首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 343 毫秒
1.
In this paper, we systematically analyze the pressure projection stabilization method for the Darcy and coupled Darcy–Stokes flow problems in multiple dimensions. Stability results for this stabilization method are established. For the Darcy flow, optimal error estimates in the divergence norm for velocity and suboptimal error estimates in the $L^{2}$ -norm for pressure are obtained, and a superconvergence result for the pressure is derived; a local postprocessing scheme is constructed to generate optimal error estimates in the L 2-norm for pressure. For the coupled Darcy–Stokes flow, error estimates of optimal order are obtained in terms of the energy norm of velocity and pressure. Numerical results are presented to check the theory developed.  相似文献   

2.
The theory of variational inequalities enables us to formulate and solve free boundary problems in fixed domains, while most other methods assume the position of the unknown domain in solving the problem. Here the problem of seepage flow through a rectangular dam with a free boundary is formulated as a vertical inequality following the ideas of Baiocchi. In order to demonstrate the essential ideas of extending the domain of the solution of problems with free boundaries, the problem of the deflection, of a string on a rigid support is first examined. Next, variational inequalities are derived which are associated with several cases of seepage problems. An approximation theory, including a priori error estimates, is developed using finite element methods, and an associated numerical scheme is given. It is shown that for linear and quadratic finite element methods, the rates of convergence are 0(h) and 0(h1.25-δ), 0 < δ < 0.25, respectively, if the permeability is constant.  相似文献   

3.
Damage models are capable of representing crack initiation and mimicking crack propagation within a continuum framework. Thus, in principle, they do not describe crack openings. In durability analyses of concrete structures however, transfer properties are a key issue controlled by crack propagation and crack opening. We extend here a one‐dimensional approach for estimating a crack opening from a continuum‐based finite element calculation to two‐dimensional cases. The technique operates in the case of mode I cracking described in a continuum setting by a nonlocal isotropic damage model. We used the global tracking method to compute the idealized crack location as a post‐treatment procedure. The original one‐dimensional problem devised in Dufour et al. [4] is recovered as profiles of deformation orthogonal to the idealized crack direction are computed. An estimate of the crack opening and an error indicator are computed by comparing finite element deformation profiles and theoretical profiles corresponding to a displacement discontinuity. Two estimates have been considered: In the strong approach, the maxima of the profiles are assumed to be equal; in the weak approach, the integrals of each profile are set equal. Two‐dimensional numerical calculations show that the weak estimates perform better than do the strong ones. Error indicators, defined as the distance between the numerical and theoretical profiles, are less than a few percentages. In the case of a three‐point bending, test results are in good agreement with experimental data, with an error lower than 10% for widely opened crack (> 40µm). Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

4.
In this paper, we formulate a finite-element procedure for approximating the coupled fluid and mechanics in Biot’s consolidation model of poroelasticity. We approximate the flow variables by a mixed finite-element space and the displacement by a family of discontinuous Galerkin methods. Theoretical convergence error estimates are derived and, in particular, are shown to be independent of the constrained specific storage coefficient, c o . This suggests that our proposed algorithm is a potentially effective way to combat locking, or the nonphysical pressure oscillations, which sometimes arise in numerical algorithms for poroelasticity.  相似文献   

5.
We develop and analyze a mixed finite element method for the solution of an elliptic system modeling a porous medium with large cavities, called vugs. It consists of a second-order elliptic (i.e., Darcy) equation on part of the domain coupled to a Stokes equation on the rest of the domain, and a slip boundary condition (due to Beavers–Joseph–Saffman) on the interface between them. The tangential velocity is not continuous on the interface. We consider a 2-D vuggy porous medium with many small cavities throughout its extent, so the interface is not isolated. We use a certain conforming Stokes element on rectangles, slightly modified near the interface to account for the tangential discontinuity. This gives a mixed finite element method for the entire Darcy–Stokes system with a regular sparsity pattern that is easy to implement, independent of the vug geometry, as long as it aligns with the grid. We prove optimal global first-order L 2 convergence of the velocity and pressure, as well as the velocity gradient in the Stokes domain. Numerical results verify these rates of convergence and even suggest somewhat better convergence in certain situations. Finally, we present a lower dimensional space that uses Raviart–Thomas elements in the Darcy domain and uses our new modified elements near the interface in transition to the Stokes elements.  相似文献   

6.

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.

  相似文献   

7.
In this paper we prove the convergence of a finite volume scheme for the discretization of an elliptic–parabolic problem, namely Richards equation β(P)t?div(K(β(P))× ?(P+z))=0, together with Dirichlet boundary conditions and an initial condition. This is done by means of a priori estimates in L2 and the use of Kolmogorov's theorem on relative compactness of subsets of L2.  相似文献   

8.
This work presents and analyzes, on unstructured grids, a discrete duality finite volume method (DDFV method for short) for 2D-flow problems in nonhomogeneous anisotropic porous media. The derivation of a symmetric discrete problem is established. The existence and uniqueness of a solution to this discrete problem are shown via the positive definiteness of its associated matrix. Properties of this matrix combined with adequate assumptions on data allow to define a discrete energy norm. Stability and error estimate results are proven with respect to this norm. L 2-error estimates follow from a discrete Poincaré inequality and an L ?∞?-error estimate is given for a P 1-DDFV solution. Numerical tests and comparison with other schemes (especially those from FVCA5 benchmark) are provided.  相似文献   

9.
This work concerns linearization methods for efficiently solving the Richards equation, a degenerate elliptic-parabolic equation which models flow in saturated/unsaturated porous media. The discretization of Richards’ equation is based on backward Euler in time and Galerkin finite elements in space. The most valuable linearization schemes for Richards’ equation, i.e. the Newton method, the Picard method, the Picard/Newton method and the L-scheme are presented and their performance is comparatively studied. The convergence, the computational time and the condition numbers for the underlying linear systems are recorded. The convergence of the L-scheme is theoretically proved and the convergence of the other methods is discussed. A new scheme is proposed, the L-scheme/Newton method which is more robust and quadratically convergent. The linearization methods are tested on illustrative numerical examples.  相似文献   

10.
In this paper, we formulate a finite element procedure for approximating the coupled fluid and mechanics in Biot’s consolidation model of poroelasticity. Here, we approximate the pressure by a mixed finite element method and the displacements by a Galerkin method. Theoretical convergence error estimates are derived in a discrete-in-time setting. Of particular interest is the case when the lowest-order Raviart–Thomas approximating space or cell-centered finite differences are used in the mixed formulation and continuous piecewise linear approximations are used for displacements. This approach appears to be the one most frequently applied to existing reservoir engineering simulators.  相似文献   

11.
In this paper, we formulate a finite element procedure for approximating the coupled fluid and mechanics in Biot’s consolidation model of poroelasticity. Here, we approximate the pressure by a mixed finite element method and the displacements by a Galerkin method. Theoretical convergence error estimates are derived in a continuous in-time setting for a strictly positive constrained specific storage coefficient. Of particular interest is the case when the lowest-order Raviart–Thomas approximating space or cell-centered finite differences are used in the mixed formulation, and continuous piecewise linear approximations are used for displacements. This approach appears to be the one most frequently applied to existing reservoir engineering simulators.  相似文献   

12.
We derive the governing equations for the dynamic response of unsaturated poroelastic solids at finite strain. We obtain simplified governing equations from the complete coupled formulation by neglecting the material time derivative of the relative velocities and the advection terms of the pore fluids relative to the solid skeleton, leading to a so‐called us ? pw ? pa formulation. We impose the weak forms of the momentum and mass balance equations at the current configuration and implement the framework numerically using a mixed finite element formulation. We verify the proposed method through comparison with analytical solutions and experiments of quasi‐static processes. We use a neo‐Hookean hyperelastic constitutive model for the solid matrix and demonstrate, through numerical examples, the impact of large deformation on the dynamic response of unsaturated poroelastic solids under a variety of loading conditions. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

13.
The multiple high‐pressure (HP), low‐temperature (LT) metamorphic units of Western and Central Anatolia offer a great opportunity to investigate the subduction‐ and continental accretion‐related evolution of the eastern limb of the long‐lived Aegean subduction system. Recent reports of the HP–LT index mineral Fe‐Mg‐carpholite in three metasedimentary units of the Gondwana‐derived Anatolide–Tauride continental block (namely the Afyon Zone, the Ören Unit and the southern Menderes Massif) suggest a more complicated scenario than the single‐continental accretion model generally put forward in previous studies. This study presents the first isotopic dates (white mica 40Ar–39Ar geochronology), and where possible are combined with P–T estimates (chlorite thermometry, phengite barometry, multi‐equilibrium thermobarometry), on carpholite‐bearing rocks from these three HP–LT metasedimentary units. It is shown that, in the Afyon Zone, carpholite‐bearing assemblages were retrogressed through greenschist‐facies conditions at c. 67–62 Ma. Early retrograde stages in the Ören Unit are dated to 63–59 Ma. In the Kurudere–Nebiler Unit (HP Mesozoic cover of the southern Menderes Massif), HP retrograde stages are dated to c. 45 Ma, and post‐collisional cooling to c. 26 Ma. These new results support that the Ören Unit represents the westernmost continuation of the Afyon Zone, whereas the Kurudere–Nebiler Unit correlates with the Cycladic Blueschist Unit of the Aegean Domain. In Western Anatolia, three successive HP–LT metamorphic belts thus formed: the northernmost Tav?anl? Zone (c. 88–82 Ma), the Ören–Afyon Zone (between 70 and 65 Ma), and the Kurudere–Nebiler Unit (c. 52–45 Ma). The southward younging trend of the HP–LT metamorphism from the upper and internal to the deeper and more external structural units, as in the Aegean Domain, points to the persistence of subduction in Western Anatolia between 93–90 and c. 35 Ma. After the accretion of the Menderes–Tauride terrane, in Eocene times, subduction stopped, leading to continental collision and associated Barrovian‐type metamorphism. Because, by contrast, the Aegean subduction did remain active due to slab roll‐back and trench migration, the eastern limb (below Southwestern Anatolia) of the Hellenic slab was dramatically curved and consequently teared. It therefore is suggested that the possibility for subduction to continue after the accretion of buoyant (e.g. continental) terranes probably depends much on palaeogeography.  相似文献   

14.
Numerical challenges occur in the simulation of groundwater flow problems because of complex boundary conditions, varying material properties, presence of sources or sinks in the flow domain, or a combination of these. In this paper, we apply adaptive isogeometric finite element analysis using locally refined (LR) B‐splines to address these types of problems. The fundamentals behind isogeometric analysis and LR B‐splines are briefly presented. Galerkin's method is applied to the standard weak formulation of the governing equation to derive the linear system of equations. A posteriori error estimates are calculated to identify which B‐splines should be locally refined. The error estimates are calculated based on recovery of the L2‐projected solution. The adaptive analysis method is first illustrated by performing simulation of benchmark problems with analytical solutions. Numerical applications to two‐dimensional groundwater flow problems are then presented. The problems studied are flow around an impervious corner, flow around a cutoff wall, and flow in a heterogeneous medium. The convergence rates obtained with adaptive analysis using local refinement were, in general, observed to be of optimal order in contrast to simulations with uniform refinement. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

15.
 The monoclinic titanite-like high-pressure form of calcium disilicate has been synthesized and quenched to ambient conditions to form the triclinic low-pressure phase containing silicon in four-, five- and sixfold coordination. The enthalpy of formation of the quench product has been measured by high-temperature oxide melt calorimetry. The value obtained from samples from a series of several synthesis experiments is ΔH f = (−26.32 ± 4.27) kJ mol−1 for the formation from the component oxides, or ΔH f  = (−2482.81 ± 4.59) kJ mol−1 for the formation from the elements. The result is identical within experimental error to available estimates, although the previously predicted energy difference between the monoclinic and triclinic phases could not be verified. Received: 16 February 2000 / Accepted: 14 July 2000  相似文献   

16.
We present a numerical scheme for reactive contaminant transport with nonequilibrium sorption in porous media. The mass conservative scheme is based on Euler implicit, mixed finite elements, and Newton method. We consider the case of a Freundlich-type sorption. In this case, the sorption isotherm is not Lipschitz but just Hölder continuous. To deal with this, we perform a regularization step. The convergence of the scheme is analyzed. An explicit order of convergence depending only on the regularization parameter, the time step, and the mesh size is derived. We give also a sufficient condition for the quadratic convergence of the Newton method. Finally, relevant numerical results are presented.  相似文献   

17.
Small-scale faults with associated drag folds in brittle-ductile rocks can retain detailed information on the kinematics and amount of deformation the host rock experienced. Measured fault orientation (α), drag angle (β) and the ratio of the thickness of deflected layers at the fault (L) and further away (T) can be compared with α, β and L/T values that are calculated with a simple analytical model. Using graphs or a numerical best-fit routine, one can then determine the kinematic vorticity number and initial fault orientation that best fits the data. The proposed method was successfully tested on both analogue experiments and numerical simulations with BASIL. Using this method, a kinematic vorticity number of one (dextral simple shear) and a minimum finite strain of 2.5–3.8 was obtained for a population of antithetic faults with associated drag folds in a case study area at Mas Rabassers de Dalt on Cap de Creus in the Variscan of the easternmost Pyrenees, Spain.  相似文献   

18.
A four‐node plane parametric element AQGβ6‐I is constructed on the basis of the quadrilateral area coordinate, the generalized conforming principle and the projection technique with a penalty factor β within an interval of 0–1. When β = 0, the element has excellent bending performance. When β = 1, the element can pass patch test strictly; its performance is as good as many famous elements. When β value is between 0 and 1, such as β = 0.5, the element can arrive at a compromise between (relatively) low sensitivity to mesh distortion and perfect convergence. The work provides an illuminating method to alleviate a difficult problem in finite element modelling using the four‐node quadrilateral element, which can pass the strict patch test, but has poor performance in bending dominated problem; on the contrary, it has excellent performance in bending dominated problem but cannot pass the strong patch test. The AQGβ6‐I with the convergence formulation (β = 1) is then applied to coupled solid‐deformation/fluid‐flow simulation for porous geomaterials. The computational examples are carried out to demonstrate that the AQGβ6‐I (β = 1) element is not only stable, reliable and efficient but also of high accuracy. The present study provides a good applicable element for finite element simulations of solid‐deformation/fluid‐flow for porous geomaterials. Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

19.
An analytical solution to 1D coupled water infiltration and deformation in layered soils is derived using a Laplace transformation. Coupling between seepage and deformation, and initial conditions defined by arbitrary continuous pore‐water pressure distributions are considered. The analytical solutions describe the transient pore‐water pressure distributions during 1D, vertical infiltration toward the water table through two‐layer unsaturated soils. The nonlinear coupled formulations are first linearized and transformed into a form that is solvable using a Laplace transformation. The solutions provide a reliable means of comparing the accuracy of various numerical methods. Parameters considered in the coupled analysis include the saturated permeability (ks), desaturation coefficient (α), and saturated volumetric water content (θs) of each soil layer, and antecedent and subsequent rainfall infiltration rates. The analytical solution demonstrates that the coupling of seepage and deformation plays an important role in water infiltration in layered unsaturated soils. A smaller value of α or a smaller absolute value of the elastic modulus of the soil with respect to a change in soil suction (H) for layered unsaturated soils means more marked coupling effect. A smaller absolute value of H of the upper layer soil also tends to cause more marked coupling effect. A large difference between the saturated coefficients of permeability for the top and bottom soil layers leads to reduced rainfall infiltration into the deep soil layer. The initial conditions also play a significant role in the pore‐water pressure redistribution and coupling effect. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

20.
Strain rates from snowball garnet   总被引：3，自引：0，他引：3  
Spiral inclusion trails in garnet porphyroblasts are likely to have formed due to simultaneous growth and rotation of the crystals, during syn‐metamorphic deformation. Thus, they contain information on the strain rate of the rock. Strain rates may be interpreted from such inclusion trails if two functions are known: (1) The relationship between rotation rate and shear strain rate; (2) the growth rate of the crystal. We have investigated details of both functions using a garnetiferous mica schist from the eastern European Alps as an example. The rotation rate of garnet porphyroblasts was determined using finite element modelling of the geometrical arrangement of the crystals in the rock. The growth rate of the porphyroblasts was determined by using the major and trace element distributions in garnet crystals, thermodynamic pseudosections and information on the grain size distribution. For the largest porphyroblast size fraction (size L=12 mm) we constrain a growth interval between 540 and 590 °C during the prograde evolution of the rock. Assuming a reasonable heating rate and using the angular geometry of the spiral inclusion trails we are able to suggest that the mean strain rate during crystal growth was of the order of =6.6 × 10?14 s?1. These estimates are consistent with independent estimates for the strain rates during the evolution of this part of the Alpine orogen.  相似文献   

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

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