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

A. Al‐Omari  E. Masad 《国际地质力学数值与分析法杂志》2004,28(13):1327-1360

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.  相似文献   

A comparison of theoretical Sv emission line strengths with extreme ultraviolet observations of a sunspot     

F. P. Keenan  J. G. Doyle 《Solar physics》1990,128(2):345-351

Electron impact excitation rates for transitions in the S v ion, calculated with theR-matrix code, are used to derive the electron temperature sensitive emission line ratiosR ₁ =I(854.8 Å)/I(786.9 Å),R ₂ =I(852.2 Å)/I(786.9 Å),R ₃ =I(849.2 Å)/I(786.9 Å), andR ₄ =I(1199.1 Å)/I(786.9 Å), which are found to be significantly different from previous estimates. A comparison of the present results with observational data for a sunspot obtained with the Harvard S-055 spectrometer on boardSkylab reveals generally good agreement between theory and experiment, except in the case ofR ₁, which is probably due to blending in the 854.8 Å feature. The possible effects of Lyman continuum absorption on the observed line ratios is briefly discussed.  相似文献   

Waveform inversion of seismic refraction data and applications to young Pacific crust   总被引：3，自引：0，他引：3  

Peter R. Shaw  John A. Orcutt 《Geophysical Journal International》1985,82(3):375-414

  相似文献   

The international law of the sea: Volume II: the late D. P. O'Connell,edited by I. A. Shearer Clarendon Press,Oxford, 1984, xxvi and 524 pp,plus 43 pp,Hardback £40     

A.V. Lowe 《Marine Policy》1985,9(1):79-80

  相似文献   

Axial convergence in a well-mixed estuary     

R.A. Nunes  J.H. Simpson 《Estuarine, Coastal and Shelf Science》1985,20(5):637-649

Transverse secondary circulations involving surface convergence, observed in a well-mixed estuary in North Wales, are made visible by the collection of surface material along an axial line which extends continuously for many kilometres through the estuary. The circulation and axial convergence, however, are seen only during the flood phase of the tide and no similar behaviour has been observed during the ebb phase.Convergent circulations in the estuary are associated with small but steady transverse density gradients in the cross-section, produced by non-uniform advection of the longitudinal gradient through the channel. A diagnostic model, using measured mean distributions of cross-sectional density, indicates surface transverse velocities (～0.1 ms^?1) similar to those observed in the estuary. The model further predicts appreciable transverse divergent currents at a fractional depth of 0.75: a prediction which has been tested in the estuary using a vertical array of accurately resolving current direction indicators.  相似文献   

The dust coma of periodic comet Churyumov-Gerasimenko (1982 VIII)     

M.S. Hanner  E. Tedesco  A.T. Tokunaga  G.J. Veeder  D.F. Lester  F.C. Witteborn  J.D. Bregman  J. Gradie  L. Lebofsky 《Icarus》1985,64(1):11-19

The dust coma of Comet P/Churyumov-Gerasimenko was monitored in the infrared (1–20 μm) from September 1982 to March 1983. Maximum dust production rate of ～2 × 10⁵ g/sec occured in December, 1 month postperihelion. The ratio of dust/gas production was higher than that in other short-period comets. No silicate feature was visible in the 8- to 13-μm spectrum on 23 October. The mean geometric albedo of the grains was ～0.04 at 1.25 μm and ～0.05 at 2.2 μm.  相似文献   

A WKBJ spectral method for computation of SV synthetic seismograms in a cylindrically symmetric medium     

Chang-Eob Baag  Charles A. Langston 《Geophysical Journal International》1985,80(2):387-417

Summary. A method of synthetic seismogram computation for teleseismic SV -waves is developed in order to treat quantitatively SV -waves in problems of body wave source inversion and source—receiver structure studies. The method employs WKBJ theory for a generalized ray in a vertically inhomogeneous half-space and the propagator matrix technique for waves in near-surface homogeneous layers. Wavenumber integration is done along the real axis of the wavenumber plane and anelasticity is included by using complex velocity in all regions of the earth model. The near-surface source structure is taken into account in the computation for the case of the shallow source by allowing a point source to be located in the homogeneous layers. Source and receiver area structures are also allowed to differ. A general moment tensor point source is considered.  相似文献   

Rapport de la commission de finance     

L. A. Adams 《Bulletin of Volcanology》1940,6(1):76-79

  相似文献   

Simulations of multimode bolometric interferometers     

George Saklatvala  Michael P. Hobson  Stafford Withington 《Monthly notices of the Royal Astronomical Society》2006,367(3):1194-1200

We describe a procedure for the numerical modelling of astronomical interferometers, with particular relevance to far-infrared and submillimetre wavelengths. The scheme is based on identifying a set of modes that carry power from the sky to the detector. The procedure is extremely general, and can be used to model scalar or vector fields, in any state of coherence and polarization, the only limitation being that the propagation of a coherent field through the system be described by an integral transform, a constraint that is in practise always met.
We present simulations of ideal, multimode two-dimensional interferometers, and show that the modal theory reproduces the correct behaviour of both Michelson and Fizeau interferometers. We calculate simulated visibility data for a multimode bolometric Michelson interferometer, with a synthesized source, and produce a dirty map, recovering the original source with the usual artefacts associated with interferometers.  相似文献   

Thick and thin models of the evolution of carbon dioxide on Mars     

Curtis V. Manning  Christopher P. McKay 《Icarus》2006,180(1):38-59

We present results from a new simulation code that accounts for the evolution of the reservoirs of carbon dioxide on Mars, from its early years to the present. We establish a baseline model parameter set that produces results compatible with the present (i.e., Patm?6.5 mbar with permanent CO₂ ice cap) for a wide range of initial inventories. We find that the initial inventory of CO₂ broadly determines the evolutionary course of the reservoirs of CO₂. The reservoirs include the atmosphere, ice cap, adsorbed CO₂ in the regolith, and carbonate rocks. We track the evolution of the free inventory: the atmosphere, ice cap and regolith. Simulations begin at 4.53 Gyr before present with a rapid loss of free inventory to space in the early Noachian. Models that assume a relatively small initial inventory (?5 bar) have pronounced minima in the free inventory of CO₂ toward the end of the Noachian. Under baseline parameters, initial inventories below ∼4.5 bar result in a catastrophic loss of the free inventory to space. The current free inventory would be then determined by the balance between outgassing, sputtering losses and chemical weathering following the end of the late bombardment. We call these “thin” models. They generically predict small current free inventories in line with expectations of a small present CO₂ ice cap. For “thick” models, with initial inventories ?5 bar, a surplus of 300-700 mbar of free CO₂ remains during the late-Noachian. The histories of free inventory in time for thick models tend to converge within the last 3.5 Gyr toward a present with an ice cap plus atmospheric inventory of about 100 mbar. For thick models, the convergence is largely due to the effects of chemical weathering, which draws down higher free inventories more rapidly than the low. Thus, thick models have ?450 mbar carbonate reservoirs, while thin models have ?200 mbar. Though both thick and thin scenarios can reproduce the current atmospheric pressure, the thick models imply a relatively large current CO₂ ice cap and thin models, little or none. While the sublimation of a massive cap at a high obliquity would create a climate swing of greenhouse warming for thick models, under the thin model, mean temperatures and pressures would be essentially unaffected by increases in obliquity.  相似文献