首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 15 毫秒
1.
  相似文献   

2.
The structure of shocks propagating through partially ionized hydrogen gas with characteristics typical of the atmospheres of RR Lyr, W Vir, and RV Tau type variables is analyzed in terms of a self-consistent solution of the equations of gas dynamics, atomic kinetics, and radiation transfer. The solutions were obtained for shock waves with velocities 20 km/s≤U 1≤90 km/s and unperturbed hydrogen gas with temperatures 3000 K≤T 1≤9000 K and density ρ1=10?10 g/cm3. The fraction of the energy of the gas-dynamic flux converted into radiation increases with the shock amplitude, and the ratio of the radiation flux emitted by the shock to the gas kinetic energy flux is 0.4???0.8 for the velocities U 1 considered. This ratio also increases slightly with the ambient gas temperature T 1 due to an increase in hydrogen ionization in the radiative precursor. The flux emitted by the leading edge of the shock opposite to the gas flow is several percent higher than the flux emitted in the opposite direction by the trailing edge of the shock. Radiation is mostly concentrated in the Balmer continuum, and the region of efficient Lyman radiation transfer includes gas layers located near the viscous jump (δX=±104 cm). The final gas-compression ratio in units of the limiting compression corresponding to an isothermal approximation is virtually independent of the shock amplitude, and increases with the unperturbed gas temperature from r≈0.5 at T 1=3000 K to r≈0.9 at T 1=9000 K.  相似文献   

3.
In many geostatistical applications, spatially discretized unknowns are conditioned on observations that depend on the unknowns in a form that can be linearized. Conditioning takes several matrix–matrix multiplications to compute the cross-covariance matrix of the unknowns and the observations and the auto-covariance matrix of the observations. For large numbers n of discrete values of the unknown, the storage and computational costs for evaluating these matrices, proportional to n 2, become strictly inhibiting. In this paper, we summarize and extend a collection of highly efficient spectral methods to compute these matrices, based on circulant embedding and the fast Fourier transform (FFT). These methods are applicable whenever the unknowns are a stationary random variable discretized on a regular equispaced grid, imposing an exploitable structure onto the auto-covariance matrix of the unknowns. Computational costs are reduced from O(n 2) to O(nlog2 n) and storage requirements are reduced from O(n 2) to O(n).  相似文献   

4.
We present an a priori stability and convergence analysis of a new mixed discontinuous Galerkin scheme applied to the instationary Darcy problem. The analysis accounts for a spatially and temporally varying permeability tensor in all estimates. The proposed method is stabilized using penalty terms in the primary and the flux unknowns.  相似文献   

5.
Variations in the seawater 87Sr/86Sr curve through time can be caused by fluctuations in the strontium flux or variations in the isotopic ratio from at least six different sources and sinks. Thus, 12 or more parameters control each single measurement although widely accepted assumptions allow this to be reduced to typically six unknowns. Interpreting the causes of time-variation in the seawater 87Sr/86Sr curve is therefore hampered by inherent non-uniqueness. However, this problem is under-constrained rather than unconstrained. As a result, whilst there are an infinite number of possible interpretations, these all come from a few families of very similar solutions. Using this insight, it is possible to find solutions having the smallest possible variations in source flux or source 87Sr/86Sr ratio. Thus, lower-bounds can be placed upon the source variations responsible for the observed fluctuations in the seawater 87Sr/86Sr curve. When applied to the evolution of the Early Jurassic 87Sr/86Sr seawater curve, this approach demonstrates that a short-lived Toarcian event is genuine since it is present in all models, regardless of the values chosen for the unknown source fluxes and unknown source isotope ratios. However, the variations in strontium flux or isotopic ratio necessary to explain the Toarcian event may be significantly smaller than would be predicted assuming modern values for the unknown parameters.  相似文献   

6.
A sequential solution procedure is used to simulate compositional two-phase flow in porous media. We employ a multiphysics concept that adapts the numerical complexity locally according to the underlying processes to increase efficiency. The framework is supplemented by a local refinement of the simulation grid. To calculate the fluxes on such grids, we employ a combination of the standard two-point flux approximation and a multipoint flux approximation where the grid is refined. This is then used to simulate a large-scale example related to underground CO2 storage.  相似文献   

7.
Natural gas, consisting primarily of methane(CH4), has become a major source of clean energy in modern society in many parts of the globe. Recent experimental observations and discoveries of deep-sourced abiotic CH4 in cold subduction zones indicate the important ability of cold subducted slabs to generate natural gas reservoirs. However, most CH4 flux and reservoirs remain unknown and their potential is overlooked in global carbon flux estimations. Massive abiot...  相似文献   

8.
Volvach  L. N.  Volvach  A. E.  Larionov  M. G.  Wolak  P.  Kramer  B.  Menten  K.  Kraus  A.  Brand  J.  Zanichelli  A.  Poppi  S.  Rigini  S.  Ipatov  A. V.  Ivanov  D. V.  Mikhailov  A. G.  Mel’nikov  A. 《Astronomy Reports》2019,63(8):652-665

The most powerful flare ever registered in the Galactic water-maser source W49N has been detected in long-term monitoring data in the 616–523 transition with line frequency f = 22.235 GHz carried out on the 22-m Simeiz, 32-m Toruń, 100-m Effelsberg, and 32-m Medicina radio telescopes, beginning in September 2017 and continuing in 2018. Some stages of the flare were monitored daily. Detailed variations of the source spectral flux density with time have been obtained. At the flare maximum, the flux exceeded P ≈ 8 × 104 Jy, and this was record highest flux registered over the entire history of observations of this source. Important conclusions related to details of the mechanism for the H2O line emission have been drawn. An exponential increase in the flare flux density was detected during both the rise and decline of the flare. The data obtained indicate that the maser is unsaturated, and remained in this state up to the maximum observed flux densities. Additional support for the idea that the maser is unsaturated is the shape of the dependence of the line width on the flux. The characteristics of the variations of the spectral flux density are probably associated with a sharp increase in the density of the medium and the photon flux that led to an increase in the temperature from an initial level of 10–40 K to hundreds of Kelvins. Interferometric maps of the object during the increase in the spectral flux density of the flare have been obtained. A possible mechanism for the primary energy release in W49N is considered.

  相似文献   

9.
The Nu Expression for Probabilistic Data Integration   总被引：4，自引：0，他引：4  
The general problem of data integration is expressed as that of combining probability distributions conditioned to each individual datum or data event into a posterior probability for the unknown conditioned jointly to all data. Any such combination of information requires taking into account data interaction for the specific event being assessed. The nu expression provides an exact analytical representation of such a combination. This representation allows a clear and useful separation of the two components of any data integration algorithm: individual data information content and data interaction, the latter being different from data dependence. Any estimation workflow that fails to address data interaction is not only suboptimal, but may result in severe bias. The nu expression reduces the possibly very complex joint data interaction to a single multiplicative correction parameter ν 0, difficult to evaluate but whose exact analytical expression is given; availability of such an expression provides avenues for its determination or approximation. The case ν 0=1 is more comprehensive than data conditional independence; it delivers a preliminary robust approximation in presence of actual data interaction. An experiment where the exact results are known allows the results of the ν 0=1 approximation to be checked against the traditional estimators based on assumption of data independence.  相似文献   

10.
A method for history matching of an in-house petroleum reservoir compositional simulator with multipoint flux approximation is presented. This method is used for the estimation of unknown reservoir parameters, such as permeability and porosity, based on production data and inverted seismic data. The limited-memory Broyden–Fletcher–Goldfarb–Shanno method is employed for minimization of the objective function, which represents the difference between simulated and observed data. In this work, we present the key features of the algorithm for calculations of the gradients of the objective function based on adjoint variables. The test example shows that the method is applicable to cases with anisotropic permeability fields, multipoint flux approximation, and arbitrary fluid compositions.  相似文献   

11.
In this work, lowest-order Raviart–Thomas and Brezzi–Douglas–Marini mixed methods are considered for groundwater flow simulations. Typically, mixed methods lead to a saddle-point problem, which is expensive to solve. Two approaches are numerically compared here to allow an explicit velocity elimination: (1) the well-known hybrid formulation leading to a symmetric positive definite system where the only unknowns are the Lagrange multipliers and (2) a more recent approach, inspired from the multipoint flux approximation method, reducing low-order mixed methods to cell-centered finite difference schemes. Selected groundwater flow scenarios are used for the comparison between hybrid and multipoint approaches. The simulations are performed in the bidimensional case with a general triangular discretization because of its practical interest for hydrogeologists.  相似文献   

12.
The types and structures of adsorption complexes formed by oxalate at boehmite (γ-AlOOH)/water and corundum (α-Al2O3)/water interfaces were determined using in situ attenuated total reflectance fourier transform infrared (ATR-FTIR) spectroscopy and quantum chemical simulation methods. At pH 5.1, at least four different oxalate species were found at or near the boehmite/water interface for oxalate surface coverages (Γox) ranging from 0.25 to 16.44 μmol/m2. At relatively low coverages (Γox < 2.47), strongly adsorbed inner-sphere oxalate species (IR peaks at 1286, 1418, 1700, and 1720 cm−1) replace weakly adsorbed carbonate species, and a small proportion of oxalate anions are adsorbed in an outer-sphere mode (IR peaks at 1314 and 1591 cm−1). IR peaks indicative of inner-sphere adsorbed oxalate are also observed for oxalate at the corundum/water interface at Γox = 1.4 μmol/m2. With increasing oxalate concentration (Γox > 2.47 μmol/m2), the boehmite surface binding sites for inner-sphere adsorbed oxalate become saturated, and excess oxalate ions are present dominantly as aqueous species (IR peaks at 1309 and 1571 cm−1). In addition to these adsorption processes, oxalate-promoted dissolution of boehmite following inner-sphere oxalate adsorption becomes increasingly pronounced with increasing Γox and results in an aqueous Al(III)-oxalate species, as indicated by shifted IR peaks (1286 → 1297 cm−1 and 1418 → 1408 cm−1). At pH 2.5, no outer-sphere adsorbed oxalate or aqueous oxalate species were observed. The similarity of adsorbed oxalate spectral features at pH 2.5 and 5.1 implies that the adsorption mechanism of aqueous HOx− species involves loss of protons from this species during the ligand-exchange reaction. As a consequence, adsorbed inner-sphere oxalate and aqueous Al(III)-oxalate complexes formed at pH 2.5 have coordination geometries very similar to those formed at pH 5.1.The coordination geometry of inner-sphere adsorbed oxalate species was also predicted using quantum chemical geometry optimization and IR vibrational frequency calculations. Geometry-optimized Al8O12 and Al14O22 clusters with the reactive surface Al site coordinated by three oxygens were used as model substrates for corundum and boehmite surfaces. Among the models considered, calculated IR frequencies based on a bidentate side-on structure with a 5-membered ring agree best with the observed frequencies for boehmite/oxalate/water samples at Γox = 0.25 to 16.44 μmol/m2 and pH 2.5 and 5.1, and for a corundum/oxalate/water sample at Γox = 1.4 μmol/m2 and pH 5.1. Based on these results, we suggest that oxalate bonding on boehmite and corundum surfaces results in 5-coordinated rather than 4- or 6-coordinated Al surface sites.  相似文献   

13.
We present a general compositional formulation using multi-point flux mixed finite element (MFMFE) method on general hexahedral grids. The mixed finite element framework allows for local mass conservation, accurate flux approximation, and a more general treatment of boundary conditions. The multi-point flux inherent in MFMFE scheme allows the usage of a full permeability tensor. The proposed formulation is an extension of single and two-phase flow formulations presented by Wheeler and Yotov, SIAM J. Numer. Anal. 44(5), 2082–2106 (35) with similar convergence properties. Furthermore, the formulation allows for black oil, single-phase and multi-phase incompressible, slightly and fully compressible flow models utilizing the same design for different fluid systems. An accurate treatment of diffusive/dispersive fluxes owing to additional velocity degrees of freedom is also presented. The applications areas of interest include gas flooding, CO 2 sequestration, contaminant removal, and groundwater remediation.  相似文献   

14.
We consider the slightly compressible two-phase flow problem in a porous medium with capillary pressure. The problem is solved using the implicit pressure, explicit saturation (IMPES) method, and the convergence is accelerated with iterative coupling of the equations. We use discontinuous Galerkin to discretize both the pressure and saturation equations. We apply two improvements, which are projecting the flux to the mass conservative H(div)-space and penalizing the jump in capillary pressure in the saturation equation. We also discuss the need and use of slope limiters and the choice of primary variables in discretization. The methods are verified with two- and three-dimensional numerical examples. The results show that the modifications stabilize the method and improve the solution.  相似文献   

15.
Disorder models of oxygen positions in P63/ mmc, C2221 and P212121 tridymites were given in applying geometrical and lattice dynamical calculations. Sixmembered rings of rigid SiO4 units are all collapsed in these forms; with silicon atoms fixed, SiO4 units can take six different orientations in forming tridymite frameworks in both the P63/mmc and C2221 forms, and three orientations in the P212121 form. Atomic distances and angles obtained from the distance least-squares method are about equal for the three forms: 〈Si-O〉 (mean Si-O) = 1.611 Å, 〈O-O〉 = 2.629 Å, and 〈Si-O-Si〉 = 147°. Domain formation models are given for the three forms. The tridymite framework structures may possibly undergo lattice vibrations with low frequencies in two kinds of pair-wise rotational modes of SiO4 units joined by the apical oxygen atoms, at the Γ-point: one is around 〈100〉 (or 〈210〉 for the hexagonal case), and the other is around 〈010〉. As temperature approaches the hexagonal-orthorhombic transition from below, the rotational mode around 〈100〉 remarkably softens at the Γ-point. The behavior of the atoms at the hexagonal-orthorhombic transition is explained in terms of a coupled softening of the two rotational modes of neighboring local domains in different orientations.  相似文献   

16.
Proper characterizations of background soil CO2 respiration rates are critical for interpreting CO2 leakage monitoring results at geologic sequestration sites. In this paper, a method is developed for determining temperature-dependent critical values of soil CO2 flux for preliminary leak detection inference. The method is illustrated using surface CO2 flux measurements obtained from the AmeriFlux network fit with alternative models for the soil CO2 flux versus soil temperature relationship. The models are fit first to determine pooled parameter estimates across the sites, then using a Bayesian hierarchical method to obtain both global and site-specific parameter estimates. Model comparisons are made using the deviance information criterion (DIC), which considers both goodness of fit and model complexity. The hierarchical models consistently outperform the corresponding pooled models, demonstrating the need for site-specific data and estimates when determining relationships for background soil respiration. A hierarchical model that relates the square root of the CO2 flux to a quadratic function of soil temperature is found to provide the best fit for the AmeriFlux sites among the models tested. This model also yields effective prediction intervals, consistent with the upper envelope of the flux data across the modeled sites and temperature ranges. Calculation of upper prediction intervals using the proposed method can provide a basis for setting critical values in CO2 leak detection monitoring at sequestration sites.  相似文献   

17.
为了探讨接触变质带内变碳酸盐岩变质过程CO2释放的数量和排放CO2的物理、化学及地质条件,根据递进变质反应和时间积分流体通量模型,定量分析和定量计算了双山地区变碳酸盐岩在接触变质作用中释放CO2的通量。计算结果得到CO2的通量值为0.729×104~2.446×104 mol/cm2,CO2的来源以接触变质反应释放为主;CO2的生成释放与变质程度呈正相关关系。自白云石带至方解石带变质流体中XCO2不断升高,但钙铝榴石带由于岩浆水影响,流体通量最高而XCO2急剧下降。   相似文献   

18.
A systematic analysis ofmethods for computing the trajectories of solid-phase particles applied in modern astrophysics codes designed for modeling gas–dust circumstellar disks has been carried out for the first time. Themotion of grains whose velocities are determinedmainly by the gas drag, that is, for which the stopping time or relaxation time for the velocity of the dust to the velocity of the gas tstop is less than or comparable to the rotation period, are considered. The methods are analyzed from the point of view of their suitability for computing the motions of small bodies, including dust grains less than 1 μm in size, which are strongly coupled to the gas. Two test problems are with analytical solutions. Fast first order accurate methods that make it possible to avoid additional restrictions on the time step size τ due to gas drag in computations of the motion of grains of any size are presented. For the conditions of a circumstellar disk, the error in the velocity computations obtained when using some stable methods becomes unacceptably large when the time step size is τ > tstop. For the radial migration of bodies that exhibit drifts along nearly Keplerian orbits, an asymptotic approximation, sometimes called the short friction time approximation or drift flux model, gives a relative error for the radial-velocity computations equals to St2, where St is the Stokes number, the ratio of the stopping time of the body to some fraction of the rotation period (dynamical time scale) in the disk.  相似文献   

19.
Lakes worldwide are commonly oversaturated with CO2, however the source of this CO2 oversaturation is not well understood. To examine the magnitude of the C flux to the atmosphere and determine if an excess of respiration (R) over gross primary production (GPP) is sufficient to account for this C flux, metabolic parameters and stable isotopes of dissolved O2 and C were measured in 23 Québec lakes. All of the lakes sampled were oversaturated with CO2 over the sampling period, on average 221 ± 25%. However, little evidence was found to conclude that this CO2 oversaturation was the result of an excess of pelagic R over GPP. In lakes Croche and à l’Ours, where CO2 flux, R and GPP were measured weekly, the annual difference between pelagic GPP and R, or net primary production (NPP), was not sufficient to account for the size of the CO2 flux to the atmosphere. In Lac Croche average annual NPP was 14.4 mg C m−2 d−1 while the average annual flux of CO2 to the atmosphere was 34 mg C m−2 d−1. In Lac à l’Ours average annual NPP was −9.1 mg C m−2 d−1 while the average annual flux of CO2 to the atmosphere was 55 mg C m−2 d−1. In all of the lakes sampled, O2 saturation averaged 104.0 ± 1.7% during the ice-free season and the isotopic composition of dissolved O2 (δ18ODO) was 22.9 ± 0.3‰, lower than atmospheric values and indicative of net autotrophy. Carbon evasion was not a function of R, nor did the isotopic signature of dissolved CO2 in the lakes present evidence of excess R over GPP. External inputs of C must therefore subsidize the lake to explain the continued CO2 oversaturation. The isotopic composition of dissolved inorganic C (δ13CDIC) indicates that the CO2 oversaturation cannot be attributed to in situ aerobic respiration. δ13CDIC reveals a source of excess C enriched in 13C, which may be accounted for by anaerobic sediment respiration or groundwater inputs followed by kinetic isotope fractionation during degassing under open system conditions.  相似文献   

20.
The carbon cycle of global inland waters is quantitatively comparable to other components in the global carbon budget. Among inland waters, a significant part is man-made lakes formed by damming rivers. Man-made lakes are undergoing a rapid increase in number and size. Human impacts and frequent algae blooms lead to it necessary to make a better constraint on their carbon cycles. Here, we make a primary estimation on the air–water CO2 transfer flux through an algae bloom year for a subtropical man-made lake—Hongfeng Lake, Southwest China. To do this a new type of glass bottles was designed for content and isotopic analysis of DIC and other environmental parameters. At the early stage of algae bloom, CO2 was transferred from the atmosphere to the lake with a net flux of 1.770 g·C·m?2. Later, the partial pressure (pCO2) of the aqueous CO2 increased rapidly and the lake outgassed to the atmosphere with a net flux of 95.727 g·C·m?2. In the remaining days, the lake again took up CO2 from the atmosphere with a net flux of 14.804 g·C·m?2. As a whole, Lake Hongfeng released 4527 t C to the atmosphere, accounting for one-third of the atmosphere/soil CO2 sequestered by chemical weathering in the whole drainage. With an empirical mode decomposition method, we found air temperature plays a major role in controlling water temperature, aqueous pCO2 and hence CO2 flux. This work indicates a necessity to make detailed and comprehensive carbon budgets in man-made lakes.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号