首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
Accurate simulation of fluid pressures in layered reservoirs with strong permeability contrasts is a challenging problem. For this purpose, the Discontinuous Galerkin (DG) method has become increasingly popular. Unfortunately, standard linear solvers are usually too inefficient for the aforementioned application. To increase the efficiency of the conjugate gradient (CG) method for linear systems resulting from symmetric interior penalty (discontinuous) Galerkin (SIPG) discretizations, we cast an existing two-level preconditioner into the deflation framework. The main idea is to use coarse corrections based on the DG solution with polynomial degree p = 0. This paper provides a numerical comparison of the performance of the original preconditioner and the resulting deflation variant in terms of scalability and overall efficiency. Furthermore, it studies the influence of the SIPG penalty parameter, weighted averages in the SIPG formulation (SWIP), the smoother, damping of the smoother, and the strategy for solving the coarse systems. We have found that the penalty parameter can best be chosen diffusion-dependent. In that case, both two-level methods yield fast and scalable convergence. Whether preconditioning or deflation is to be favored depends on the choice of the smoother and on the damping of the smoother. Altogether, both two-level methods can contribute to cheaper and more accurate fluid pressure simulations.  相似文献   

2.

Hydrochemical analysis of the high-salinity lakes in the Ishim Plain (>250–300 g/L) located at the border with the Northern Kazakhstan uranium ore province is performed. The studies have shown that the main factor of concentration and redistribution of uranium in the lake basins of the Ishim Plain are the processes of intense salt deflation causing sanding of lakes and uranium depletion in the near-surface layer of the bottom deposits. The correlation between the hydroxide forms of uranium binding in the bottom lacustrine deposits of the Ishim Plain and the coffinite composition of the Semizbai deposit makes it possible to consider this province to be promising for the discovery of hydromineral uranium deposits.

  相似文献   

3.

In this work, we measure the performance of the fixed stress split algorithm for the immiscible water-oil flow coupled with linear poromechanics. The two-phase flow equations are solved on general hexahedral elements using the multipoint flux mixed finite element method whereas the poromechanics equations are discretized using the conforming Galerkin method. We introduce a rigorous calculation of the update in poroelastic properties during the iterative solution of the coupled system equations. The effects of the coupling parameter on the performance of the fixed stress algorithm is demonstrated in two field studies: the Frio oil reservoir and the Cranfield injection site.

  相似文献   

4.
Zhu  S. R.  Wu  L. Z.  Huang  J. 《Computational Geosciences》2022,26(1):131-145

This paper studies the potential of using the successive over-relaxation iteration method with polynomial preconditioner (P(m)-SOR) to solve variably saturated flow problems described by the linearized Richards’ equation. The finite difference method is employed to numerically discretize and produce a system of linear equations. Generally, the traditional Picard method needs to re-evaluate the iterative matrix in each iteration, so it is time-consuming. And under unfavorable conditions such as infiltration into extremely dry soil, the Picard method suffers from numerical non-convergence. For linear iterative methods, the traditional Gauss-Seidel iteration method (GS) has a slow convergence rate, and it is difficult to determine the optimum value of the relaxation factor w in the successive over-relaxation iteration method (SOR). Thus, the approximate optimum value of w is obtained based on the minimum spectral radius of the iterative matrix, and the P(m)-SOR method is extended to model underground water flow in unsaturated soils. The improved method is verified using three test examples. Compared with conventional Picard iteration, GS and SOR methods, numerical results demonstrate that the P(m)-SOR has faster convergence rate, less computation cost, and good error stability. Besides, the results reveal that the convergence rate of the P(m)-SOR method is positively correlated with the parameter m. This method can serve as a reference for numerical simulation of unsaturated flow.

  相似文献   

5.

We consider the finite element (FE) approximation of the two dimensional shallow water equations (SWE) by considering discretizations in which both space and time are established using a stable FE method. Particularly, we consider the automatic variationally stable FE (AVS-FE) method, a type of discontinuous Petrov-Galerkin (DPG) method. The philosophy of the DPG method allows us to establish stable FE approximations as well as accurate a posteriori error estimators upon solution of a saddle point system of equations. The resulting error indicators allow us to employ mesh adaptive strategies and perform space-time mesh refinements, i.e., local time stepping. We establish a priori error estimates for the AVS-FE method and linearized SWE and perform numerical verifications to confirm corresponding asymptotic convergence behavior. In an effort to keep the computational cost low, we consider an alternative space-time approach in which the space-time domain is partitioned into finite sized space-time slices. Hence, we can perform adaptive mesh refinements on each individual slice to preset error tolerances as needed for a particular application. Numerical verifications comparing the two alternatives indicate the space-time slices are superior for simulations over long times, whereas the solutions are indistinguishable for short times. Multiple numerical verifications show the adaptive mesh refinement capabilities of the AVS-FE method, as well the application of the method to some commonly applied benchmarks for the SWE.

  相似文献   

6.
This paper describes the first use of inter-particle force measurement in reworked aerosols to better understand the mechanics of dust deflation and its consequent ecological ramifications. Dust is likely to carry hydrocarbons and micro-organisms including human pathogens and cultured microbes and thereby is a threat to plants, animals and human. Present-day global aerosol emissions are substantially greater than in 1850; however, the projected influx rates are highly disputable. This uncertainty, in part, has roots in the lack of understanding of deflation mechanisms. A growing body of literature shows that whether carbon emission continues to increase, plant transpiration drops and soil water retention enhances, allowing more greenery to grow and less dust to flux. On the other hand, a small but important body of geochemistry literature shows that increasing emission and global temperature leads to extreme climates, decalcification of surface soils containing soluble carbonate polymorphs and hence a greater chance of deflation. The consistency of loosely packed reworked silt provides background data against which the resistance of dust’s bonding components (carbonates and water) can be compared. The use of macro-scale phenomenological approaches to measure dust consistency is trivial. Instead, consistency can be measured in terms of inter-particle stress state. This paper describes a semi-empirical parametrisation of the inter-particle cohesion forces in terms of the balance of contact-level forces at the instant of particle motion. We put forward the hypothesis that the loss of Ca2+-based pedogenic salts is responsible for much of the dust influx and surficial drying pays a less significant role.  相似文献   

7.

Phosphorus (P) plays an important role in the eutrophication of river and marine environments. The adsorption-desorption processes of P by estuarine sediment were studied to better understand the P behaviour in the Loughor Estuary in Carmarthen Bay, UK. Three types of models were used to estimate the P adsorption isotherm for five sediment samples from the Loughor Estuary. The Langmuir adsorption isotherm showed a good fit with the linear (data) isotherm. The P desorption from bed sediment was measured as the easy exchangeable P (MgCl2 extractable). A negative correlation was found between adsorbed P in bed sediment and median particle size. A positive correlation was also found between P adsorbed and metals such as iron (Fe), aluminium (Al) and calcium (Ca) as well as the total organic carbon (TOC).

  相似文献   

8.
ABSTRACT

Collision of Gondwana and Laurentia in the late Palaeozoic created new topography, drainages, and foreland basin systems that controlled sediment dispersal patterns on southern Laurentia. We utilize sedimentological and detrital zircon data from early Permian (Cisuralian/Leonardian) submarine-fan deposits in the Midland Basin of west Texas to reconstruct sediment dispersal pathways and palaeogeography. New sedimentological data and wire-line log correlation suggest a portion of the early Permian deposits have a southern entry point. A total of 3259 detrital zircon U-Pb and 357 εHf data from 12 samples show prominent groups of zircon grains derived from the Appalachian (500–270 Ma) and Grenville (1250–950 Ma) provinces in eastern Laurentia and the peri-Gondwana terranes (800–500 Ma) incorporated in the Alleghanian-Ouachita-Marathon orogen. Other common zircon groups of Mesoproterozoic-Archaean age are also present in the samples. The detrital zircon data suggest throughout the early Permian, Appalachia and Gondwana detritus was delivered by a longitudinal river system that flowed along the Appalachian-Ouachita-Marathon foreland into the Midland Basin. Tributary channels draining the uplifted Ouachita-Marathon hinterland brought Gondwana detritus into the longitudinal river with headwaters in the Appalachians or farther northeast. This drainage extended downstream westward and delivered sediments into the Permian Basin near the west terminus of the Laurentia-Gondwana suture. Estimated rates of deposition and proportions of zircons from more local (Grenville) versus more distal (Pan-African) sources indicate that river strength decreased throughout early Permian time. Primary sediment delivery pathway was augmented by minor input from the Ancestral Rocky Mountains and wind deflation of fluvial sediments north and east of the basin. Slope failure associated with early Permian deposition in the southeastern margin of the Midland Basin triggered gravity flows leading to submarine fan deposition.  相似文献   

9.
Pan  Yucong  Liu  Quansheng  Kong  Xiaoxuan  Liu  Jianping  Peng  Xingxin  Liu  Qi 《Acta Geotechnica》2019,14(4):1249-1268

In this study, determination of some machine parameters and performance prediction for tunnel boring machine (TBM) are conducted based on laboratory rock cutting test. Firstly, laboratory full-scale linear cutting test is carried out using 432-mm CCS (constant cross section) disc cutter in Chongqing Sandstone. Then, the input parameters for TBM cutterhead design are extracted; some TBM specifications are determined and then compared to the manufactured values. Finally, laboratory full-scale linear cutting test results are compared with the field TBM excavation performance data collected in Chongqing Yangtze River Tunnel. Results show that laboratory full-scale linear cutting test results, combined with some engineering considerations, can be used for the preliminary and rough design of TBM machine capacity. Meanwhile, combined with some modification factors, it can also well predict the field TBM excavation performance.

  相似文献   

10.
《Sedimentology》2018,65(6):1859-1874
Ripples are prevalent in aeolian landscapes. Many researchers have focused on the shape and formation of sand ripples, but few have studied the differences in the particle size of sand on crests and in troughs along bed, especially the variations caused by changes in friction velocity and the wind‐blowing duration. A particle size of 158 μm (d ) was used to create aeolian ripples in a wind tunnel under four friction velocities (u *) with different wind duration times (t ). Samples were collected from the surfaces of ripple crests and troughs, respectively, at seven sites, and particle sizes were measured using a Malvern Mastersizer 2000. The main results were: (i) The particle size distributions of sand in troughs are unimodal with slight variations of particle size parameters, including mean particle size, standard deviation, skewness and kurtosis, etc., under different conditions, while these particle size parameters of sand on crests change with friction velocity and deflation time. Moreover, some of the particle distribution curves for the sand on crests do not follow typical unimodal curves. (ii) With increasing friction velocity or deflation duration, the sand on the crests shows a coarsening process relative to those on the bed surface. The particle size of sand on crests at a 1 m bed increases linearly with friction velocity (=  344·27 + 34·54 u *) at a given wind‐blowing duration. The particle sizes of sand on crests at 1 m, 2 m and 4 m beds increase with a power‐law relationship (= a + t b, where a and b are fitting parameters) with deflation time at a given friction velocity. (iii) The probability cumulative curves of sand showed a three‐section pattern in troughs and on most of the crests but a four‐section pattern at crest locations due to increased influence by friction velocity and deflation time. The proportions of the sediment moved by suspension, saltation and creep in the three‐section pattern were within the ranges of 0·2% to 2·0%, 97·0% to 98·9%, and 0·8% to 3·0%, respectively. For the four‐section pattern, suspension accounted for 0·3% and 3·0%, and the proportion of creep increased with friction velocity and deflation time, while saltation decreased accordingly. Although these results require additional validation, they help to advance current understanding of the grain‐size characteristics of aeolian ripples.  相似文献   

11.
Wu  Hui-ming  Ma  Ning  Ma  Quan-kun  Lin  Xiao-fei  Song  Ci 《Acta Geotechnica》2021,16(9):2997-3004

We present an aerosol injection technique (AIT) to accelerate the consolidation of soft soils for ground improvement. We employ high-pressure aerosol injections at different depths to enhance the drainage in soft soils for faster consolidation. The technique is briefly described. A well-instrumented field test is carried out to demonstrate its performance. Compared to the traditional methods, our approach gives rise to faster dissipation of excess pore pressure and larger ground settlement. This method is particularly attractive for the improvement in soft ground in medium depths.

  相似文献   

12.
Yadav  Ganesh  Singh  R. B.  Anand  Subhash  Pandey  B. W.  Mohanty  Ashutosh  Dash  Sushree Sangita 《GeoJournal》2021,87(4):469-483

Ambient air pollution, particularly in the urban environment of developing countries, has turned out to be a major health risk factor. We explore the compounded impact of age sensitivity, exposure, poverty, co-morbidity, etc., along with composite air pollution in determining morbidity and health burden of people in Lucknow, India. This cross-sectional study is confined to analyse respiratory health status across different socio-economic and geographic locations using n = 140 in-depth questionnaire method. We used mean daily ambient air pollution data of PM10, PM2.5, SO2, and NO2 for the 2008–2018 period. We used the ecological model framework to assess the risk at different hierarchical levels and compounded severity on a spatial scale. We also used Logistic regression model with log odds and odds ratio to analyze the association of risks outcomes with composite air pollution scores calculated using the principal component analysis method. There is a strong association of location-specific respiratory disease prevalence with an overall 32 percent prevalence. The prevalence of ecological model 1 (individual domain) is 4.3 percent, while ecological model 2 (community domain) has the highest prevalence at 32.4 percent. The logistic regression model shows that respiratory disease load is positively associated with age sensitivity (P < .001) and composite pollution level (P < .001). For another model with suffocation as the outcome variable, composite pollution level (P < .001) and exposure (P < .001) are positively associated. Optimum interventions are required at Ecological models 1, 2, and 3 levels for better respiratory health outcomes.

  相似文献   

13.

Minimization of a stochastic cost function is commonly used for approximate sampling in high-dimensional Bayesian inverse problems with Gaussian prior distributions and multimodal posterior distributions. The density of the samples generated by minimization is not the desired target density, unless the observation operator is linear, but the distribution of samples is useful as a proposal density for importance sampling or for Markov chain Monte Carlo methods. In this paper, we focus on applications to sampling from multimodal posterior distributions in high dimensions. We first show that sampling from multimodal distributions is improved by computing all critical points instead of only minimizers of the objective function. For applications to high-dimensional geoscience inverse problems, we demonstrate an efficient approximate weighting that uses a low-rank Gauss-Newton approximation of the determinant of the Jacobian. The method is applied to two toy problems with known posterior distributions and a Darcy flow problem with multiple modes in the posterior.

  相似文献   

14.

This paper introduces a simplified method to investigate the influence of thermal loads on the shaft friction and tip resistance of energy piles. The method is based on the influence factors (λ and η) which are back-calculated drawing on a large number of field and model tests. Values for λ and η during heating and cooling are suggested. Moreover, a new equation is proposed to calculate total shaft friction. The equations concerning the relationship between η and temperature difference are recommended to investigate the impacts of the thermal load on the pile tip resistance. The slope of the linear equation of an end-bearing pile is 2.14 times that of a floating pile indicating that the pile tip resistance of an end-bearing pile is much more affected by the same thermal load.

  相似文献   

15.
Hasan  Md. Rakib  Yu  Xinbao  Abu-Farsakh  Murad 《Acta Geotechnica》2021,16(2):491-506

Load tests of drilled shafts are often performed using Osterberg cell (O-cell) testing, a popular load test method for drilled shafts, which measures both side and tip resistance. However, it is common that only one of the resistance components can be fully mobilized. Therefore, extrapolation of the partially mobilized resistance is often required to determine the total resistance or the equivalent top-down curve. The extrapolation tends to introduce errors to the constructed total resistance values, which subsequently affect the calibrated resistance factors required for the LRFD design of drilled shafts. In this study, eight O-cell tests of drilled shafts with total measured resistances close to the failure criteria defined by FHWA, 5% of the shaft diameter (B), were collected among 64 drilled shaft load tests from Louisiana and Mississippi. For each of the eight cases, extrapolation was performed on both tip and side movement curves for the construction of the equivalent top-down load-settlement (ELT) curves. Data points from the measured side or tip movement curve were removed systematically to create a total of 80 cases with partially mobilized movement curves, and extrapolation exercises were performed on each fabricated case to obtain its equivalent top-down curve. The error of bias for each fabricated case was determined for statistical analyses. Multiple linear regression analysis was performed on the bias errors to model the bias errors caused by extrapolation. Calibrated resistance factors were determined and compared between the original database and fabricated database needing extrapolation. A correction method is proposed, based on a linear regression relationship, to estimate and minimize the extrapolation error of bias for less mobilized databases.

  相似文献   

16.
Gong  Hongmei  Guo  Pengyuan  Chen  Shuo  Duan  Meng  Sun  Pu  Wang  Xiaohong  Niu  Yaoling 《中国地球化学学报》2020,39(3):355-364

Element doping has been proved to be a useful method to correct for the mass bias fractionation when analyzing iron isotope compositions. We present a systematic re-assessment on how the doped nickel may affect the iron isotope analysis in this study by carrying out several experiments. We find three important factors that can affect the analytical results, including the Ni:Fe ratio in the analyte solutions, the match of the Ni:Fe ratio between the unknown sample and standard solutions, and the match of the Fe concentration between the sample and standard solutions. Thus, caution is required when adding Ni to the analyte Fe solutions before analysis. Using our method, the δ56Fe and δ57Fe values of the USGS standards W-2a, BHVO-2, BCR-2, AGV-2 and GSP-2 are consistent with the recommended literature values, and the long-term (one year) external reproducibility is better than 0.03 and 0.05‰ (2SD) for δ56Fe and δ57Fe, respectively. Therefore, the analytical method established in our laboratory is a method of choice for high quantity Fe isotope data in geological materials.

  相似文献   

17.
The present study examines the geoarchaeological history of an oasis in Kharga Depression in central Egypt. El‐Deir is renowned for its Ptolemaic temple and Roman fortress on the road from former Hibis (Kharga) to the Nile Valley. During the survey, spring mounds and irrigation soils belonging to an ancient agricultural zone were discovered, and further documented by ceramics found on the site. Our methodology combines the geomorphological interpretation of landforms (especially yardangs) with ceramics and 14C‐dated charcoal to distinguish and date former agricultural areas in El‐Deir. The results show that the oasis experienced several phases of soil accretion and destruction through time. Playa sediments were deposited in the humid early Holocene and severely eroded by deflation before the onset of irrigated agriculture between Pharaonic and Persian times. Very fast vertical soil accretion occurred in the Ptolemaic period, but irrigation soils were later destroyed during the Roman period by a combination of wind deflation and flash floods (second to fourth century A.D.), suggesting a period of climate instability. The case of El‐Deir invites reevaluation of constructive agencies for the development of irrigated land and destructive agencies as limiting factors for the sustainability of agricultural practices in late antiquity.  相似文献   

18.
Artyukh  V. S. 《Astronomy Reports》2019,63(3):167-173

An expression for the intensity of synchrotron emission from a radio source (in the optically thin regime) in terms of the energy densities in the magnetic field and particles is obtained, based on a definition of a relativistic electron related to its rest energy. A relationship is obtained between the energy densities in particles Ee and the magnetic field EH for a physical system containing a magnetic field and relativistic electrons in a minimum-energy state. A method for estimating the Doppler factors of the relativistic electrons has been developed. This method does not requires that all radio sources have the same radiation energies (brightness temperatures): it is sufficient that the energies of the magnetic fields and relativistic particles in the source be approximately equal. The method yields Doppler-factor estimates with reasonably good accuracy, even when there are modest deviations from energy equipartition in the radio source,making it applicable to many radio sources. The method is used to estimate the Doppler factor of the radio jet in CTA 21.

  相似文献   

19.
di Prisco  Claudio  Flessati  Luca  Porta  Davide 《Acta Geotechnica》2020,15(4):1013-1030

The fronts of tunnels excavated under particularly difficult ground conditions by employing conventional tunnelling methods are commonly supported: the stabilization is usually achieved either by improving the mechanical properties of the soil (injections, jet grouting, soil freezing, etc.) or by introducing linear inclusions. This last technique, consisting in the introduction of pipes (usually made of fibreglass reinforced polymers) in the front, is particularly popular since it is very simple to adapt the reinforcement geometry, length and number to the different conditions encountered during the excavation. The design of this reinforcement technique is nowadays based on very simplified approaches: on either empirical formula or the limit equilibrium method. In a previous paper, the authors numerically studied the mechanical response of unreinforced fronts in cohesive soils and defined a non-dimensional front characteristic curve. In this paper, the authors intend to take into consideration the role of reinforcements by following the same approach. A procedure allowing the definition of the reinforced non-dimensional front characteristic curve, once the reinforcement pattern is assigned, is introduced. The practical use of this curve is described.

  相似文献   

20.
Kadeethum  T.  Lee  S.  Nick  H. M. 《Mathematical Geosciences》2020,52(8):977-1015

We study and compare five different combinations of finite element spaces for approximating the coupled flow and solid deformation system, so-called Biot’s equations. The permeability and porosity fields are heterogeneous and depend on solid displacement and fluid pressure. We provide detailed comparisons among the continuous Galerkin, discontinuous Galerkin, enriched Galerkin, and two types of mixed finite element methods. Several advantages and disadvantages for each of the above techniques are investigated by comparing local mass conservation properties, the accuracy of the flux approximation, number of degrees of freedom (DOF), and wall and CPU times. Three-field formulation methods with fluid velocity as an additional primary variable generally require a larger number of DOF, longer wall and CPU times, and a greater number of iterations in the linear solver in order to converge. The two-field formulation, a combination of continuous and enriched Galerkin function space, requires the fewest DOF among the methods that conserve local mass. Moreover, our results illustrate that three out of the five methods conserve local mass and produce similar flux approximations when conductivity alteration is included. These comparisons of the key performance indicators of different combinations of finite element methods can be utilized to choose the preferred method based on the required accuracy and the available computational resources.

  相似文献   

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

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