Three dimensional simulation of fluid flow in X‐ray CT images of porous media

A numerical scheme is developed in order to simulate fluid flow in three dimensional (3‐D) microstructures. The governing equations for steady incompressible flow are solved using the semi‐implicit method for pressure‐linked equations (SIMPLE) finite difference scheme within a non‐staggered grid system that represents the 3‐D microstructure. This system allows solving the governing equations using only one computational cell. The numerical scheme is verified through simulating fluid flow in idealized 3‐D microstructures with known closed form solutions for permeability. The numerical factors affecting the solution in terms of convergence and accuracy are also discussed. These factors include the resolution of the analysed microstructure and the truncation criterion. Fluid flow in 2‐D X‐ray computed tomography (CT) images of real porous media microstructure is also simulated using this numerical model. These real microstructures include field cores of asphalt mixes, laboratory linear kneading compactor (LKC) specimens, and laboratory Superpave gyratory compactor (SGC) specimens. The numerical results for the permeability of the real microstructures are compared with the results from closed form solutions. Copyright © 2004 John Wiley & Sons, Ltd.     

Zircon U–Pb age, trace element and Hf isotope composition of Kongling terrane in the Yangtze Craton: refining the timing of Palaeoproterozoic high-grade metamorphism

U–Pb age, trace element and Hf isotope compositions of zircon were analysed for a metasedimentary rock and two amphibolites from the Kongling terrane in the northern part of the Yangtze Craton. The zircon shows distinct morphological and chemical characteristics. Most zircon in an amphibolite shows oscillatory zoning, high Th/U and ¹⁷⁶Lu/¹⁷⁷Hf ratios, high formation temperature, high trace element contents, clear negative Eu anomaly, as well as HREE-enriched patterns, suggesting that it is igneous. The zircon yields a weighted mean ²⁰⁷Pb/²⁰⁶Pb age of 2857 ± 8 Ma, representing the age of the magmatic protolith. The zircon in the other two samples is metamorphic. It has low Th/U ratios, low trace element concentrations, variable HREE contents (33.8 ≥ Lu_N≥2213; 14.7 ≤ Lu_N/Sm_N ≤ 354) and ¹⁷⁶Lu/¹⁷⁷Hf ratios (0.000030–0.001168). The data indicate that the zircon formed in the presence of garnet and under upper amphibolite facies conditions. The metamorphic zircon yields a weighted mean ²⁰⁷Pb/²⁰⁶Pb age of 2010 ± 13 Ma. These results combined with previously obtained Palaeoproterozoic metamorphic ages suggest a c. 2.0 Ga Palaeoproterozoic collisional event in the Yangtze Craton, which may result from the assembly of the supercontinent Columbia. The zircon in two samples yields weighted mean two-stage Hf model ( T _DM2) ages of 3217 ± 110 and 2943 ± 50 Ma, respectively, indicating that their protoliths were mainly derived from Archean crust.     

Age and early metamorphic history of the Sanbagawa belt: Lu–Hf and P–T constraints from the Western Iratsu eclogite

Two distinct age estimates for eclogite-facies metamorphism in the Sanbagawa belt have been proposed: (i) c.  120–110 Ma based on a zircon SHRIMP age for the Western Iratsu unit and (ii) c.  88–89 Ma based on a garnet–omphacite Lu–Hf isochron age from the Seba and Kotsu eclogite units. Despite the contrasting estimates of formation ages, petrological studies suggest the formation conditions of the Western Iratsu unit are indistinguishable from those of the other two units—all ∼20 kbar and 600–650 °C. Studies of the associated geological structures suggest the Seba and Western Iratsu units are parts of a larger semi-continuous eclogite unit. A combination of geochronological and petrological studies for the Western Iratsu eclogite offers a resolution to this discrepancy in age estimates. New Lu–Hf dating for the Western Iratsu eclogite yields an age of 115.9 ± 0.5 Ma that is compatible with the zircon SHRIMP age. However, petrological studies show that there was significant garnet growth in the Western Iratsu eclogite before eclogite facies metamorphism, and the early core growth is associated with a strong concentration of Lu. Pre-eclogite facies garnet (Grt1) includes epidote–amphibolite facies parageneses equilibrated at 550–650 °C and ∼10 kbar, and this is overgrown by prograde eclogite facies garnet (Grt2). The Lu–Hf age of c.  116 Ma is strongly skewed to the isotopic composition of Grt1 and is interpreted to reflect the age of the pre-eclogite phase. The considerable time gap ( c.  27 Myr) between the two Lu–Hf ages suggests they may be related to separate tectonic events or distinct phases in the evolution of the Sanbagawa subduction zone.     

Metamorphic phase relations in orthopyroxene-bearing granitoids: implication for high-pressure metamorphism and prograde melting in the continental crust

In this work, the factors controlling the formation and preservation of high-pressure mineral assemblages in the metamorphosed orthopyroxene-bearing metagranitoids of the Sandmata Complex, Aravalli-Delhi Mobile Belt (ADMB), northwestern India have been modelled. The rocks range in composition from farsundite through quartz mangerite to opdalite, and with varying K₂O, Ca/(Ca + Na)_rock and FeO_tot + MgO contents. A two stage metamorphic evolution has been recorded in these rocks.
An early hydration event stabilized biotite with or without epidote at the expense of magmatic orthopyroxene and plagioclase. Subsequent high-pressure granulite facies metamorphism (∼15 kbar, ∼800 °C) of these hydrated rocks produced two rock types with contrasting mineralogy and textures. In the non-migmatitic metagranitoids, spectacular garnet ± K-feldspar ± quartz corona was formed around reacting biotite, plagioclase, quartz and/or pyroxene. In contrast, biotite ± epidote melting produced migmatites, containing porphyroblastic garnet incongruent solids and leucosomes.
Applying NCKFMASHTO T–M (H₂O) and P–T pseudosection modelling techniques, it is demonstrated that the differential response of these magmatic rocks to high-pressure metamorphism is primarily controlled by the scale of initial hydration. Rocks, which were pervasively hydrated, produced garnetiferous migmatites, while for limited hydration, the same metamorphism formed sub-solidus garnet-bearing coronae. Based on the sequence of mineral assemblage evolution and the mineral compositional zoning features in the two metagranitoids, a clockwise metamorphic P–T path is constrained for the high-pressure metamorphic event. The finding has major implications in formulating geodynamic model of crustal amalgamation in the ADMB.     

Merging of low-mass systems and the origin of the fundamental plane

We present the preliminary results of a study of how small stellar systems merge to form larger ones. As we display the families of galaxies in the μ_e - R_e plane (effective surface brightness versus effective radius) we realize that different morphological types occupy different loci, evidencing the different physical mechanisms operating in each family. As proposed by Capaccioli et al. (1992) this diagram is the logical equivalent of the HR diagram for stars. Here we take some initial steps in understanding of how we can establish the evolutionary tracks, solely due to dynamical processes, in the μ_e - R_e plane, ultimately making a dwarf elliptical to turn into a normal elliptical galaxy. This revised version was published online in August 2006 with corrections to the Cover Date.     

Ultraviolet reflectance properties of asteroids

An analysis of the UV spectra of 28 asteroids obtained with the Internal Ultraviolet Explorer (IUE) satellite is presented. The spectra lie within the range 2100–3200 Å. Our results are examined in terms of both asteroid classification and of current ideas concerning the surface mineralogy of asteroids. For all the asteroids examined, UV reflectivity declines approximately linearly toward shorter wavelengths. In general, the same taxonomic groups are seen in the UV as in the visible and IR, although there is some evidence for asteroids with anomalous UV properties and for UV subclasses within the S class. No mineral absorption features are reported of strength similar to the strongest features in the visible and IR regions, but a number of shallow absorptions do occur and may provide valuable information on the surface composition of many asteroids.     

Electron microscopy study of the iron meteorite Santa Catharina

Abstract— Characterization of the microstructural features of the metal of the Santa Catharina meteorite was performed using a variety of electron optical techniques. Sample USNM#6293 is chemically homogeneous on the micron scale and has a Ni content of 28.2 wt.%. Its microstructure is similar to that of the Twin City ataxite and contains clear taenite II, i.e., fcc taenite with domains of tetrataenite, < 10 nm in size. Sample USNM#3043 is a more typical Santa Catharina specimen with dark and light regions as observed with the light optical microscope. The dark regions are inhomogeneous and contain 45–50 wt.% Ni and 7–12 wt.% O. The light regions are homogeneous and contain 35 wt.% Ni and no detectable oxygen. The microstructure is that of cloudy zone, i.e., islands of tetrataenite, ～20 nm in size, in a honeycomb matrix. The honeycomb phase contains Ni rich oxide in the dark regions and contains metal, fcc taenite, in the light regions. The original metal structure of USNM#3043 is cloudy zone which formed during cooling into the low temperature miscibility gap of the Fe-Ni phase diagram. The dark regions were developed from the metal by selective corrosion of the honeycomb structure, transforming it into Ni containing oxides, possibly non-stoichiometric Fe₂NiO₄ while retaining the tetrataenite islands. Using the results of this study, many of the existing discrepancies concerning the microstructure of Santa Catharina can be explained.     

Effect of particle acceleration process on the flare characteristics of blazars

Following the kinetic equation approach, we study the flare processes in blazars in the optical-to-X-ray region, considering energy dependent acceleration time-scale of electrons and synchrotron and adiabatic cooling as their dominant energy loss processes.     

Eastern Aphrodite Terra on Venus: Characteristics,structure, and mode of origin

Eastern Aphrodite Terra and Western Aphrodite form an altimetrically prominent 14,000 km long part of the equatorial highlands on Venus. Several parallel linear discontinuities striking northwest across the general east-west regional strike of the highlands are mapped in the altimetric and radar image data of Eastern Aphrodite and identified on the basis of abrupt termination of rift-like central chasma, offset and segmentation of the center of the highlands, and radar image discontinuities in the lowlands to the north. These characteristics are similar to those of linear discontinuities previously mapped in Western Aphrodite in terms of length, orientation, and influence on the central highlands and adjacent lowlands.Altimetric profiles in directions parallel to the discontinuities are regionally symmetric, more ridge-like in Eastern Aphrodite compared to the plateau-dominated form of topography in Western Aphrodite, and are characterized by alternating paired ridge-and-trough forms near their crests and on their flanks. By mapping the center of symmetry in multiple profiles, the prominent segmentation of the highland is shown to be imparted by an offset of the regional symmetry along the mapped discontinuities. These characteristics are morphologically similar to several of the large-scale characteristics of divergent plate boundaries of Earth, including mid-ocean rise crests and rifts, offset at fracture zones and transform faults, and symmetric thermal boundary layer topography.The altitude of the surface in profiles parallel to the discontinuities decreases as the square root of distance from the symmetry axes and with a slope similar to that predicted for thermal boundary layer topography associated with rates of divergence on Venus of ~ 1 ± 0.5 cm/yr. In order to test the hypothesis that the linear discontinuities are analogous to fracture zones, the predicted altitude of the surface at great distance from the centers of symmetry of the central highland and in directions across the discontinuities was calculated on the basis of a thermal boundary layer topography model with offset of altimetric symmetry at each discontinuity. Similarity of observed Arecibo high-resolution altimetric profiles across the discontinuities with that calculated for thermal boundary layer topography offset by transform faults reveals that in terms of the sense and magnitude of regional steps in altimetry across discontinuities and the altitude of the surface, Eastern Aphrodite is similar to the known characteristics of crustal spreading at divergent boundaries. The plateau-like form of Western Aphrodite and the ridge-like form of Eastern Aphrodite are analogous respectively to the difference between areas of anomalous (Iceland) and normal crustal production along rise crests on Earth. Estimates of volumetric differences in crustal production in the environment of Venus and as it would be influenced by differences in mantle temperature beneath Western and Eastern Aphrodite imply that Eastern Aphrodite represents normal crustal production. On this basis, Western Aphrodite may be characterized by a mantle temperature that is warmer than the mantle beneath Eastern Aphrodite Terra, perhaps in association with deep convective mantle upwelling.'Geology and Tectonics of Venus', special issue edited by Alexander T. Basilevsky (USSR Acad. of Sci., Moscow), James W. Head (Brown University, Providence), Gordon H. Pettengill (MIT. Cambridge, Massachusetts) and R. S. Saunders (J.P.L., Pasadena).