共查询到20条相似文献，搜索用时 78 毫秒
1.
2.
An unconditionally stable, fully explicit and highly precise multiple timescale finite element modeling scheme is described for a fully coupled hydro-mechanical (FCHM) analysis of saturated poroelastic media. The finite element method (FEM) is used for the discretization of the FCHM differential equation in the space domain. Direct integration is performed based on the precise time step integration method (PTSIM) for the time derivatives. Two configurations for the proposed scheme are constructed (abbreviated as PTSIM-f1 and -f2, respectively). The stability and convergence of the PTSIM-f1 and -f2 are proved using a matrix-based spectral analysis in the time domain. It is demonstrated that the explicit scheme proposed in this paper is unconditionally stable and independent of the time-step size. The algorithmic error estimation results indicate that the numerical modeling performed using PTSIM-f1 and -f2 in the time domain match the computer precision. Theoretically, the algorithmic error is caused by only the mesh discretization. Therefore, the proposed modeling scheme is a semi-analytical scheme. The applicability and accuracy of the proposed scheme are examined using sample calculations. By comparing with the analytical solutions, it is indicated that the modeling results have significant advantages over the standard FEM in terms of precision and computational efficiency for large timescales. 相似文献
3.
In this paper a finite volume (FV) numerical method is implemented to solve a Biot consolidation model with discontinuous coefficients. Our studies show that the FV scheme leads to a locally mass conservative approach which removes pressure oscillations especially along the interface between materials with different properties and yields higher accuracy for the flow and mechanics parameters. Then this numerical discretization is utilized to investigate different sequential strategies with various degrees of coupling including: iteratively, explicitly and loosely coupled methods. A comprehensive study is performed on the stability, accuracy and rate of convergence of all of these sequential methods. In the iterative and explicit solutions four splits of drained, undrained, fixed-stress and fixed-strain are studied. In loosely coupled methods three techniques of the local error method, the pore pressure method, and constant step size are considered and results are compared with other types of coupling methods. It is shown that the fixed-stress method is the best operator split in comparison with other sequential methods because of its unconditional stability, accuracy and the rate of convergence. Among loosely coupled schemes, the pore pressure and local error methods which are, respectively, based on variation of pressure and displacement, show consistency with the physics of the problem. In these methods with low number of total mechanical iterations, errors within acceptance range can be achieved. As in the pore pressure method mechanics time step increases more uniformly, this method would be less costly in comparison with the local error method. These results are likely to be useful in decision making regarding choice of solution schemes. Moreover, the stability of the FV method in multilayered media is verified using a numerical example. 相似文献
4.
Coupled formulation and algorithms for the simulation of non‐planar three‐dimensional hydraulic fractures using the generalized finite element method 
下载免费PDF全文
下载免费PDF全文 This paper presents a coupled hydro‐mechanical formulation for the simulation of non‐planar three‐dimensional hydraulic fractures. Deformation in the rock is modeled using linear elasticity, and the lubrication theory is adopted for the fluid flow in the fracture. The governing equations of the fluid flow and elasticity and the subsequent discretization are fully coupled. A Generalized/eXtended Finite Element Method (G/XFEM) is adopted for the discretization of the coupled system of equations. A Newton–Raphson method is used to solve the resulting system of nonlinear equations. A discretization strategy for the fluid flow problem on non‐planar three‐dimensional surfaces and a computationally efficient strategy for handling time integration combined with mesh adaptivity are also presented. Several three‐dimensional numerical verification examples are solved. The examples illustrate the generality and accuracy of the proposed coupled formulation and discretization strategies. Copyright © 2015 John Wiley & Sons, Ltd. 相似文献
5.
A new formulation of the element‐free Galerkin (EFG) method is developed for solving coupled hydro‐mechanical problems. The numerical approach is based on solving the two governing partial differential equations of equilibrium and continuity of pore water simultaneously. Spatial variables in the weak form, i.e. displacement increment and pore water pressure increment, are discretized using the same EFG shape functions. An incremental constrained Galerkin weak form is used to create the discrete system equations and a fully implicit scheme is used for discretization in the time domain. Implementation of essential boundary conditions is based on a penalty method. Numerical stability of the developed formulation is examined in order to achieve appropriate accuracy of the EFG solution for coupled hydro‐mechanical problems. Examples are studied and compared with closed‐form or finite element method solutions to demonstrate the validity of the developed model and its capabilities. The results indicate that the EFG method is capable of handling coupled problems in saturated porous media and can predict well both the soil deformation and variation of pore water pressure over time. Some guidelines are proposed to guarantee the accuracy of the EFG solution for coupled hydro‐mechanical problems. Copyright © 2008 John Wiley & Sons, Ltd. 相似文献
6.
The Crank–Nicolson scheme has second‐order accuracy, but often leads to oscillations affecting numerical stability. On the other hand, the implicit scheme is free from oscillation, but it has only first‐order accuracy. In this work, a three‐point discretization scheme with variable time step is presented for the time marching of parabolic partial differential equations. The method proposed has second‐order accuracy, is unconditionally stable and dampens spurious oscillations of the numerical results. The application and effectiveness of the new method are demonstrated through several numerical examples. It is shown that, unlike the Crank–Nicolson method, the approach proposed produces no oscillatory response irrespective of the time step adopted. Copyright © 2010 John Wiley & Sons, Ltd. 相似文献
7.
An alternative approach for quasi‐static large deformation analysis of saturated porous media using meshfree method 下载免费PDF全文
下载免费PDF全文 An alternative coupled large deformation formulation combined with a meshfree approach is proposed for flow–deformation analysis of saturated porous media. The formulation proposed is based on the Updated Lagrangian (UL) approach, except that the spatial derivatives are defined with respect to the configuration of the medium at the last time step rather than the configuration at the last iteration. In this way, the Cauchy stresses are calculated directly, rendering the second Piola–Kirchhoff stress tensor not necessary for the numerical solution of the equilibrium equations. Moreover, in contrast with the UL approach, the nodal shape function derivatives are calculated once in each time step and stored for use in subsequent iterations, which reduces the computational cost of the algorithm. Stress objectivity is satisfied using the Jaumann stress rate, and the spatial discretisation of the governing equations is achieved using the standard Galerkin method. The equations of equilibrium are satisfied directly, and the nonlinear parts of the system matrix are derived independent of the stresses of the medium resulting in a stable numerical algorithm. Temporal discretisation is effected based on a three‐point approximation technique that avoids spurious ripple effects and has second‐order accuracy. The radial point interpolation method is used to construct the shape functions. The application of the formulation and the significance of large deformation effects on the numerical results are demonstrated through several numerical examples. Copyright © 2015 John Wiley & Sons, Ltd. 相似文献
8.
9.
Philip Cardiff 《国际地质力学数值与分析法杂志》2015,39(13):1410-1430
Accurate prediction of the interactions between the nonlinear soil skeleton and the pore fluid under loading plays a vital role in many geotechnical applications. It is therefore important to develop a numerical method that can effectively capture this nonlinear soil‐pore fluid coupling effect. This paper presents the implementation of a new finite volume method code of poro‐elasto‐plasticity soil model. The model is formulated on the basis of Biot's consolidation theory and combined with a perfect plasticity Mohr‐Coulomb constitutive relation. The governing equation system is discretized in a segregated manner, namely, those conventional linear and uncoupled terms are treated implicitly, while those nonlinear and coupled terms are treated explicitly by using any available values from previous time or iteration step. The implicit–explicit discretization leads to a linearized and decoupled algebraic system, which is solved using the fixed‐point iteration method. Upon the convergence of the iterative method, fully nonlinear coupled solutions are obtained. Also explored in this paper is the special way of treating traction boundary in finite volume method compared with FEM. Finally, three numerical test cases are simulated to verify the implementation procedure. It is shown in the simulation results that the implemented solver is capable of and efficient at predicting reasonable soil responses with pore pressure coupling under different loading situations. Copyright © 2015 John Wiley & Sons, Ltd. 相似文献
10.
为揭示尾迹区添加横隔板对圆柱绕流流场特性影响,把多步格式引入到特征线算子分裂有限元法中,建立了基于多步格式的特征线算子分裂有限元法:在每个时间步内将Navier-Stokes方程分裂成对流项和扩散项,对流项时间离散采用多步格式,在每一子时间步内沿特征线展开并显式求解。方腔流数值模拟结果表明该算法既可降低对整体时间步长的要求又可提高计算精度。对比有无横隔板圆柱绕流流场和圆柱表面压力变化表明,横隔板可以有效地抑制绕流尾迹区涡旋脱落,提高圆柱背流面压力,减少圆柱上下表面的压力差。 相似文献
11.
Piyang Liu Gary Douglas Couples Jun Yao Zhaoqin Huang Wenhui Song Jingsheng Ma 《Computational Geosciences》2018,22(5):1187-1201
The two-scale continuum model is widely used in simulating the reactive dissolution process and predicting the optimum injection rate for carbonate reservoir acidizing treatment. The numerical methods of this model are currently based on structured grids, which are not applicable for complicated geometries. In this study, a general numerical scheme for simulating a reactive flow problem on both structured and unstructured grids is presented based on the finite volume method (FVM). The convection and diffusion terms involved in the reactive flow model are discretized by using the upwind scheme and two-point flux approximation (TPFA), respectively. The location of the centroid node inside each control volume is moved by using an optimization algorithm to make the connections with the surrounding elements as orthogonal as possible, which systematically improves the accuracy of the TPFA scheme. Additionally, in order to avoid the computational complexity resulting from the discretization of the non-linear term, the mass balance equation is only discretized in the spatial domain to get a set of ordinary differential equations (ODEs). These ODEs are coupled with the reaction equations and then solved using the numerical algorithm on ODEs. The accuracy and efficiency of the proposed method are studied by comparing the results obtained from the proposed numerical method with previous experimental and numerical results. This comparison indicates that, compared with the previous methods, the proposed method predicts the wormhole structure more accurately. Finally, the presented method is used to check the effect of the domain geometry, and it is found that the geometry of the flow domain has no effect on the optimum injection velocity, but the radial domain requires a larger breakthrough volume than the linear domain when other parameters are fixed. 相似文献
12.
在岩土工程分析中求解精度控制常常是必需的,在数值流形法中可以通过控制数学覆盖网格的稀疏和覆盖位移的阶数来达到精度的要求。提出了基于等几何分析的数值流形方法,定义了相应的数学覆盖的构造形式,推导了基于二次B样条的9节点数值流形方法分析格式;针对基于Lagrange插值函数的4节点数值流形方法提出了基于T样条思想的数学覆盖网格的局部加密方法。算例计算结果表明,相对于4节点的数值流形方法,基于非均匀有理B样条的9节点数值流形方法具有更高的精度;基于T样条思想的加密网格在保持计算精度的前提下降低了自由度的数量,表明T样条加密是一种自然的局部加密算法。 相似文献
13.
Richards方程在非饱和渗流模拟及其他相关领域应用广泛。在数值求解过程中,可以采用有限差分方法进行数值离散并迭代求解,为了获得较可靠的数值解,常规的均匀网格空间步长往往是较小的。在一些不利数值条件下,如入渗于干燥土壤,迭代计算费时甚至精度也不能得到很好改善。因此,文章提出Chebyshev空间网格改进方法,结合有限差分方法对Richards方程进行数值离散以获得线性方程组,并通过经典的Picard迭代方法进行迭代求解线性方程组以得到Richards方程的数值解。通过均质土和分层土2个不利情况下的非饱和渗流算例,又结合模型解析解和软件Hydrus-1D,对比研究了改进网格方法与均匀网格方法获得数值解的精度。结果表明,提出的Chebyshev网格方法相较于传统的均匀网格,可以在较少的节点数下获得较高的数值精度,又具有较小的计算开销,有较好的应用前景。 相似文献
14.
Unsaturated flow problems in porous media often described by Richards’ equation are of great importance in many engineering applications. In this contribution, we propose a new numerical flow approach based on isogeometric analysis (IGA) for modeling the unsaturated flow problems. The non-uniform rational B-spline (NURBS) basis is utilized for spatial discretization whereas the stable implicit backward Euler method for time discretization. The nonlinear Richards’ equation is iteratively solved with the aid of the Newton–Raphson scheme. Owing to some desirable features of an efficient numerical flow approach, major advantages of the present formulation involve: (a) numerical oscillation at the wetting front can be avoided or facilitated, simply by using either an h-refinement or a lumped mass matrix technique; (b) higher-order exactness can be obtained due to the nature of the IGA features; (c) the approach is straightforward to implement and it does not need any transformation, e.g., Kirchhoff transformation or filter algorithm; and (d) in contrast to the Picard iteration scheme, which forms linear convergences, the proposed approach can however yield quadratic convergences by using the Newton–Raphson method for solving resultant nonlinear equations. Numerical model validation is analyzed by solving a three-dimensional unsaturated flow problem in soil, and its derived results are verified against analytical solutions. Numerical applications are then studied by considering three extensive examples with simple and complex configurations to further show the accuracy and applicability of the present IGA. 相似文献
15.
In this paper, the numerical model of the transverse vibrations of a thin poroelastic plate saturated by a fluid was proposed. Two coupled dynamic equations of equilibrium related to the plate deflection and the equivalent moment were established for an isotropic porous medium with uniform porosity. The fundamental solutions for a porous plate were derived both in the Laplace transform domain and in the time domain. A meshless method was developed and demonstrated in the Laplace transform domain for solving two coupled dynamic equations. Numerical examples demonstrated the accuracy of the method of the fundamental solutions and comparisons were made with analytical solutions. The proposed meshless method was shown to be simple to implement and gave satisfactory results for a poroelastic plate dynamic analysis. Copyright © 2009 John Wiley & Sons, Ltd. 相似文献
16.
This paper presents a numerical scheme for fluid‐particle coupled discrete element method (DEM), which is based on poro‐elasticity. The motion of the particles is resolved by means of DEM. While within the proposition of Darcian regime, the fluid is assumed as a continuum phase on a Eulerian mesh, and the continuity equation on the fluid mesh for a compressible fluid is solved using the FEM. Analytical solutions of traditional soil mechanics examples, such as the isotropic compression and one‐dimensional upward seepage flow, were used to validate the proposed algorithm quantitatively. The numerical results showed very good agreement with the analytical solutions, which show the correctness of this algorithm. Sensitivity studies on the effect of some influential factors of the coupling scheme such as pore fluid bulk modulus, volumetric strain calculation, and fluid mesh size were performed to display the accuracy, efficiency, and robustness of the numerical algorithm. It is revealed that the pore fluid bulk modulus is a critical parameter that can affect the accuracy of the results. Because of the iterative coupling scheme of these algorithms, high value of fluid bulk modulus can result in instability and consequently reduction in the maximum possible time‐step. Furthermore, the increase of the fluid mesh size reduces the accuracy of the calculated pore pressure. This study enhances our current understanding of the capacity of fluid‐particle coupled DEM to simulate the mechanical behavior of saturated granular materials. Copyright © 2014 John Wiley & Sons, Ltd. 相似文献
17.
Delfim Soares Jr. 《国际地质力学数值与分析法杂志》2014,38(4):391-405
This work proposes an iterative procedure to analyze dynamic linear/nonlinear fully saturated porous media considering time‐domain finite element discretization. In this iterative approach, each phase of the coupled problem is treated separately, uncoupling the governing equations of the model. Thus, simpler, smaller, and better conditioned systems of equations are obtained, rendering more attractive techniques. A relaxation parameter is introduced in order to improve the efficiency and robustness of the iterative solution, and an expression to compute optimal values for the relaxation parameter is discussed. At the end of the paper, numerical examples are presented, illustrating the effectiveness and potentialities of the proposed methodology. Copyright © 2013 John Wiley & Sons, Ltd. 相似文献
18.
In this paper, we present a semi-implicit method for the incompressible three-phase flow equations in two dimensions. In particular, a high-order discontinuous Galerkin spatial discretization is coupled with a backward Euler discretization in time. We consider a pressure-saturation formulation, decouple the pressure and saturation equations, and solve them sequentially while still keeping each equation implicit in its respective unknown. We present several numerical examples on both homogeneous and heterogeneous media, with varying permeability and porosity. Our results demonstrate the robustness of the scheme. In particular, no slope limiters are required and a relatively large time step may be taken. 相似文献
19.
下限分析有限单元法将下限定理这一数学变分问题转化为一个数学规划问题,克服了人为构造可静应力场的困难,在实际工程中具有广阔的应用前景。通过有限元离散得到的非线性下限规划模型中包含大量的优化变量与约束条件,常规优化算法难以求解。为此,在分析非线性下限规划模型自身特点的基础上,引入可行弧技术和Wolfe非精确搜索技术改进其优化求解效率。算例分析表明,基于可行弧技术和Wolfe非精确搜索技术,下限分析有限单元法优化求解程序的收敛速度和步长搜索效率得到明显的提升,并且其数值稳定性良好、计算精度较高,可以较好地适应实际工程问题的计算。 相似文献