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.     

Nutrient transports in a Swedish estuary

Himmerfjärd is a Swedish estuary bordering on the Baltic. The estuary lacks astronomical tides and its circulation is driven by winds and freshwater runoff. Because of a tertiary sewage treatment plant located at its inner end, the estuary is becoming increasingly eutrophic. A field study was carried out for a 78-day period in late summer and early fall of 1977 to determine rates of nutrient transport and to construct nutrient budgets. Since physical parameters (current velocity, temperature, salinity, winds and water level changes) were measured more frequently than nutrients (phosphate, nitrate, nitrite and ammonium) it was necessary to develop a suitable method to calculate nutrient flux time series and net nutrient fluxes. Over the study period, Himmerfjärd imported phosphorus and exported nitrogen. Direction of nutrient fluxes and changes in flux direction were consistent with the structure of the baroclinic currents.     

Impact processing of chondritic planetesimals: Siderophile and volatile element fractionation in the Chico L chondrite

Abstract— A large impact event 500 Ma ago shocked and melted portions of the L‐chondrite parent body. Chico is an impact melt breccia produced by this event. Sawn surfaces of this 105 kg meteorite reveal a dike of fine‐grained, clast‐poor impact melt cutting shocked host chondrite. Coarse (1–2 cm diameter) globules of FeNi metal + sulfide are concentrated along the axis of the dike from metal‐poor regions toward the margins. Refractory lithophile element abundance patterns in the melt rock are parallel to average L chondrites, demonstrating near‐total fusion of the L‐chondrite target by the impact and negligible crystal‐liquid fractionation during emplacement and cooling of the dike. Significant geochemical effects of the impact melting event include fractionation of siderophile and chalcophile elements with increasing metal‐silicate heterogeneity, and mobilization of moderately to highly volatile elements. Siderophile and chalcophile elements ratios such as Ni/Co, Cu/Ga, and Ir/Au vary systematically with decreasing metal content of the melt. Surprisingly small (?10²) effective metal/silicate‐melt distribution coefficients for highly siderophile elements probably reflect inefficient segregation of metal despite the large degrees of melting. Moderately volatile lithophile elements such K and Rb were mobilized and heterogeneously distributed in the L‐chondrite impact breccias whereas highly volatile elements such as Cs and Pb were profoundly depleted in the region of the parent body sampled by Chico. Volatile element variations in Chico and other L chondrites are more consistent with a mechanism related to impact heating rather than condensation from a solar nebula. Impact processing can significantly alter the primary distributions of siderophile and volatile elements in chondritic planetesimals.     

Spectra of simulated lightning on Venus,Jupiter, and Titan

Laser-induced plasmas in various gas mixtures were used to simulate lightning in other planetary atmospheres. This method of simulation has the advantage of producing short-duration, high-temperature plasmas free from electrode contamination. The laser-induced plasma discharges in air are shown to accurately simulate terrestrial lightning and can be expected to simulate lightning spectra in other planetary atmospheres. Spectra from 240 to 880 nm are presented for simulated lightning in the atmospheres of Venus, Earth, Jupiter, and Titan. The spectra of lightning on the other giant planets are expected to be similar to that of Jupiter because the atmospheres of these planets are composed mainly of hydrogen and helium. The spectra of Venus and Titan show substantial amounts of radiation due to the presence of carbon atoms and ions and show CN Violet radiation. Although small amounts of CH₄ and NH₃ are present in the Jovian atmosphere, only emission from hydrogen and helium is observed. Most differences in the spectra can be understood in terms of the elemental ratios of the gas mixtures. Consequently, observations of the spectra of lightning on other planets should provide in situ estimates of the atmospheric and aerosol composition in the cloud layers in which lightning is occuring. In particular, the detection of inert gases such as helium should be possible and the relative abundance of these gases compared to major constituents might be determined.     

Resonant oscillations in sunspot umbrae

The Uchida and Sakurai (1975), Thomas and Scheuer (1982), and Scheuer and Thomas (1981) theory of umbral oscillations as resonant modes of magneto-acoustic-gravity waves is re-examined. For an isothermal atmosphere it is found that the quasi-Alfvén approximation is not a good approximation to the complete linearized wave equations. The new results presented here show that 3 min umbral oscillation periods are fairly insensitive to magnetic field strength above some critical value. For a detailed model umbra (Thomas and Scheuer, 1982) the calculations presented here show that 3 min umbral oscillations do not depend to any great extent on the level of forcing of the oscillations for those magnetic field strengths which are observed in sunspot umbras. Modes outside the 3 min range appear, as the lowest mode, as the level of forcing is placed at deeper and deeper levels in the solar atmosphere.     

Near-Infrared Imaging Spectrophotometry Of Comet C/1995 O1 (Hale-Bopp) Near Perihelion

We present 1- to 5-μm broadband and CVF images of comet Hale-Bopp taken 1997 February 10.5 UT, 50 days before perihelion. All the images exhibit a nonspherical coma with a bright “ridge” in the direction of the dust tail approximately 10″ from the coma. Synthetic aperture spectrophotometry implies that the optically important grains are of a radius ≤0.4 μm; smallest radius for any comet seen to date. The variation of the integrated surface brightness with radial distance from the coma (ρ) in all the images closely follows the “steady state” ρ⁻¹ model for comet dust ablation (Gehrz and Ney, 1992). The near-infrared colors taken along the dust tail are not constant implying the dust grain properties vary with coma distance. This revised version was published online in July 2006 with corrections to the Cover Date.