Finite volume approximation for an immiscible two-phase flow in porous media with discontinuous capillary pressure     

Konstantin Brenner  Clément Cancès  Danielle Hilhorst 《Computational Geosciences》2013,17(3):573-597

We consider an immiscible incompressible two-phase flow in a porous medium composed of two different rocks so that the capillary pressure field is discontinuous at the interface between the rocks. This leads us to apply a concept of multivalued phase pressures and a notion of weak solution for the flow which have been introduced in Cancès and Pierre (SIAM J Math Anal 44(2):966–992, 2012). We discretize the problem by means of a numerical algorithm which reduces to a standard finite volume scheme in each rock and prove the convergence of the approximate solution to a weak solution of the two-phase flow problem. The numerical experiments show in particular that this scheme permits to reproduce the oil-trapping phenomenon.  相似文献   

Dust torus formed by particles ejected from a celestial body at an arbitrary point of its elliptic orbit     

Sergei A. Orlov  Konstantin V. Kholshevnikov 《Celestial Mechanics and Dynamical Astronomy》2013,116(1):35-52

Dust comnplexes make one of the components of the Solar System. The surface shape of a typical dust complex consisting of particles ejected by a celestial body is found analytically, under reasonable assumptions (the main one being the smallness of perturbations). Parametric equations of the surface are obtained. The main properties of the surface are established and studied. Singular points are found, and the topological type of the surface as a whole and in the vicinity of the singular points (one conic points and one constriction) is examined.  相似文献   

Glacial isostatic adjustment in 3-D earth models: Implications for the analysis of tide gauge records along the U.S. east coast     

Jonathan E. Davis  Konstantin Latychev  Jerry X. Mitrovica  Roblyn Kendall  Mark E. Tamisiea 《Journal of Geodynamics》2008,46(3-5):90

Tide gauge records of recent sea-level change along the U.S. east coast have received significant attention within the literature of glacial isostatic adjustment (GIA). Geographic trends in these tide gauge rates are not reduced by a GIA correction based on a commonly adopted radial viscosity profile (characterized, in particular, by a lower mantle viscosity 1−2×10²¹ Pa s), and this has led to speculation that the residual trends reflect contributions from neotectonics or oceanographic processes. While the trends can be significantly reduced by adopting an Earth model with a stiffer lower mantle, such a model appears to be incompatible with independent constraints from post-glacial decay times in Hudson Bay. We use a finite-element model of the GIA process to investigate whether 3-D viscosity variations superimposed onto the “common” radial viscosity profile may provide a route to reconciling the east coast sea-level trends. We find that the specific 3-D structure we impose has little impact on the geographic trends in the GIA-corrected rates. However, we do find that the imposed lateral variations in lower mantle viscosity introduce a nearly uniform upward shift of 0.5 mm/yr in GIA-induced sea-level rates along the U.S. east coast. Thus, inferences of regional (U.S. east coast) sea-level rise due to modern melting of ice reservoirs, based on tide gauge rates corrected using 1-D GIA models, may be significantly biased by this simplifying assumption.  相似文献   

Thermal equation of state of magnesite to 32 GPa and 2073 K     

Konstantin D. Litasov  Yingwei Fei  Eiji Ohtani  Takahiro Kuribayashi  Kenichi Funakoshi 《Physics of the Earth and Planetary Interiors》2008,168(3-4):191-203

Pressure–volume–temperature relations have been measured to 32 GPa and 2073 K for natural magnesite (Mg_0.975Fe_0.015Mn_0.006Ca_0.004CO₃) using synchrotron X-ray diffraction with a multianvil apparatus at the SPring-8 facility. A least-squares fit of the room-temperature compression data to a third-order Birch–Murnaghan equation of state (EOS) yielded K₀ = 97.1 ± 0.5 GPa and K′ = 5.44 ± 0.07, with fixed V₀ = 279.55 ± 0.02 Å³. Further analysis of the high-temperature compression data yielded the temperature derivative of the bulk modulus (∂K_T/∂T)_P = −0.013 ± 0.001 GPa/K and zero-pressure thermal expansion α = a₀ + a₁T with a₀ = 4.03 (7) × 10⁻⁵ K⁻¹ and a₁ = 0.49 (10) × 10⁻⁸ K⁻². The Anderson–Grüneisen parameter is estimated to be δ_T = 3.3. The analysis of axial compressibility and thermal expansivity indicates that the c-axis is over three times more compressible (K_Tc = 47 ± 1 GPa) than the a-axis (K_Tc = 157 ± 1 GPa), whereas the thermal expansion of the c-axis (a₀ = 6.8 (2) × 10⁻⁵ K⁻¹ and a₁ = 2.2 (4) × 10⁻⁸ K⁻²) is greater than that of the a-axis (a₀ = 2.7 (4) × 10⁻⁵ K⁻¹ and a₁ = −0.2 (2) × 10⁻⁸ K⁻²). The present thermal EOS enables us to accurately calculate the density of magnesite to the deep mantle conditions. Decarbonation of a subducting oceanic crust containing 2 wt.% magnesite would result in a 0.6% density reduction at 30 GPa and 1273 K. Using the new EOS parameters we performed thermodynamic calculations for magnesite decarbonation reactions at pressures to 20 GPa. We also estimated stability of magnesite-bearing assemblages in the lower mantle.  相似文献   

Comparison of Synchrotron Radiation X-Ray Fluorescence with Conventional Techniques for the Analysis of Sedimentary Samples   总被引：3，自引：0，他引：3  

Michael A. Phedorin  Vladislav A. Bobrov  Evgeny P. Chebykin  Evgeny L. Goldberg  Michael S. Melgunov  Svetlana V. Filippova  Konstantin V. Zolotarev 《Geostandards and Geoanalytical Research》2000,24(2):205-216

A comparison is made of elemental analyses of sediment samples (lake bottom sediments and sediments from the petroliferous Bazhenovka Formation) by five competitive methods (synchrotron radiation X-ray fluorescence (SR-XRF), INAA, ICP-MS, AAS and traditional XRF). The object of this study was to identify the most suitable analytical techniques for applications in sedimentary geochemistry. Advantages and shortcomings of the five techniques were considered with respect to applications related to studies of palaeo-climate changes recorded in sedimentary cores from Lake Baikal and the geochemical stratigraphy of the Bazhenovka Formation, which is rich in organic material. We have concluded that SR-XRF was the best technique for producing extensive analytical data series, from the point of view of its speed, ease of application, cost, non-destructive nature and sensitivity, which allowed for the simultaneous determination of more than twenty elements of geochemical interest in small samples of sediments. An additional benefit of the method was demonstrated when used for the independent certification of geological reference materials. International reference samples BIL-1 and SDO-1 were certified by means of this method without use of reference materials. Analyses were based on the physics of the theory of energy transport.  相似文献   

地球物理学及大地构造学中关于地壳均衡的新模型(英文)     

Tyap.  KF 《地学前缘》2000,7(1):69-78

讨论了一些有关地壳均衡概念的模型。业已证明 ,这些模型既不能用以正确评估地球的平衡状态 ,也不能用以解释构造成因。介绍一种新的地壳均衡模型 ,它与地球的某一种旋转方式相适应。地球体与球的偏差 ,被用来作为地球平衡的判据。地球被认为是在偏差趋近于零的那些点上得到平衡的。地球体与球体之间标志的差异 ,是由地球外壳中地质作用的方向所确定的。所提出的是大地地壳均衡模型 ,是地壳中构造形成新旋转假说的基础。  相似文献   

Reply to “Discussion of ‘Study of Air-Ventilated Cavity under Model Hull on Water Surface’ by K.I. Matveev,T.J. Burnett,and A.E. Ockfen [Ocean Engineering 36 (12–13) 891–1038 (September 2009)]”     

Konstantin I. Matveev 《Ocean Engineering》2010,37(14-15):1392

  相似文献   

Static performance of power-augmented ram vehicle model on water   总被引：1，自引：0，他引：1  

Konstantin I. Matveev  Ryan K. Soderlund 《Ocean Engineering》2008,35(10):1060-1065

An experimental parametric study of a novel air-assisted platform-type model called power-augmented ram vehicle is described. The zero-speed regimes of the model over the water are investigated. The recovered thrust, pressure underneath the platform, and the model attitude are recorded for variable system geometry, loading conditions, and propulsor thrust. The stern flap under the model platform provides an effective mechanism for controlling the thrust recovery and the air-jet-induced lift. Unstable behavior of the model is found at sufficiently high levels of the propulsor thrust.  相似文献   

Preindustrial, historical, and fertilization simulations using a global ocean carbon model with new parameterizations of iron limitation, calcification, and N₂ fixation     

Konstantin Zahariev  James R. Christian  Kenneth L. Denman 《Progress in Oceanography》2008,77(1):56-82

The Canadian Model of Ocean Carbon (CMOC) has been developed as part of a global coupled climate carbon model. In a stand-alone integration to preindustrial equilibrium, the model ecosystem and global ocean carbon cycle are in general agreement with estimates based on observations. CMOC reproduces global mean estimates and spatial distributions of various indicators of the strength of the biological pump; the spatial distribution of the air-sea exchange of CO₂ is consistent with present-day estimates. Agreement with the observed distribution of alkalinity is good, consistent with recent estimates of the mean rain ratio that are lower than historic estimates, and with calcification occurring primarily in the lower latitudes. With anthropogenic emissions and climate forcing from a 1850-2000 climate model simulation, anthropogenic CO₂ accumulates at a similar rate and with a similar spatial distribution as estimated from observations. A hypothetical scenario for complete elimination of iron limitation generates maximal rates of uptake of atmospheric CO₂ of less than 1 PgC y⁻¹, or about 11% of 2004 industrial emissions. Even a ‘perfect’ future of sustained fertilization would have a minor impact on atmospheric CO₂ growth. In the long term, the onset of fertilization causes the ocean to take up an additional 77 PgC after several thousand years, compared with about 84 PgC thought to have occurred during the transition into the last glacial maximum due to iron fertilization associated with increased dust deposition.  相似文献   

Distance Between Two Arbitrary Unperturbed Orbits     

Roman V. Baluyev  Konstantin V. Kholshevnikov 《Celestial Mechanics and Dynamical Astronomy》2005,91(3-4):287-300

In the paper by Kholshevnikov and Vassilie, 1999, (see also references therein) the problem of finding critical points of the distance function between two confocal Keplerian elliptic orbits (hence finding the distance between them in the sense of set theory) is reduced to the determination of all real roots of a trigonometric polynomial of degree eight. In non-degenerate cases a polynomial of lower degree with such properties does not exist. Here we extend the results to all possible cases of ordered pairs of orbits in the Two–Body–Problem. There are nine main cases corresponding to three main types of orbits: ellipse, hyperbola, and parabola. Note that the ellipse–hyperbola and hyperbola–ellipse cases are not equivalent as we exclude the variable marking the position on the second curve. For our purposes rectilinear trajectories can be treated as particular (not limiting) cases of elliptic or hyperbolic orbits.  相似文献