多尺度有限单元法求解非均质多孔介质中的三维地下水流问题

应用多尺度有限单元法模拟非均质多孔介质中的三维地下水流问题。与传统有限单元法相比，多尺度有限单元法的基函数具有能反映单元内参数变化的优点，所以这种方法能在大尺度上抓住解的小尺度特征获得较精确的解。在介绍多尺度有限单元法求解非均质多孔介质中三维地下水流问题的基本原理之后，对参数水平方向渐变垂直方向突变的非均质多孔介质中的三维地下水流和Borden实验场的三维地下水流分别用多尺度有限单元法和传统等参有限单元法进行了计算，结果表明在模拟高度非均质多孔介质中的三维地下水流问题时，多尺度有限单元法比传统有限单元法有效，既节省计算量又有较高的精度；在模拟非均质性弱的多孔介质中的三维地下水流问题时，多尺度有限单元法虽然也能在大尺度上获得较为精确的解，但效果不明显。     

Transient seepage analysis in zoned anisotropic soils based on the scaled boundary finite‐element method

The scaled boundary finite‐element method, a semi‐analytical computational scheme primarily developed for dynamic stiffness of unbounded domains, is applied to the analysis of unsteady seepage flow problems. This method is based on the finite‐element technology and gains the advantages of the boundary element method as well. Only boundary of the domain is discretized, no fundamental solution is required and singularity problems can be modeled rigorously. Anisotropic and non‐homogeneous materials satisfying similarity are modeled with no additional efforts. In this study, firstly, formulation of the method for the transient seepage flow problems is derived followed by its solution procedures. The accuracy, simplicity and applicability of the method are demonstrated via four numerical examples of transient seepage flow – three of them are available in the literature. Homogenous, non‐homogenous, isotropic and anisotropic material properties are considered to show the versatility of the technique. Excellent agreement with the finite‐element method is observed. The method out‐performs the finite‐element method in modeling singularity points. Copyright © 2014 John Wiley & Sons, Ltd.     

A practical and efficient numerical scheme for the analysis of steady state unconfined seepage flows

The scaled boundary finite‐element method (SBFEM), a novel semi‐analytical technique, is applied to the analysis of the confined and unconfined seepage flow. This method combines the advantages of the finite‐element method and the boundary element method. In this method, only the boundary of the domain is discretized; no fundamental solution is required, and singularity problems can be modeled rigorously. Anisotropic and nonhomogeneous materials satisfying similarity are modeled without additional efforts. In this paper, SBFE equations and solution procedures for the analysis of seepage flow are outlined. The accuracy of the proposed method in modeling singularity problems is demonstrated by analyzing seepage flow under a concrete dam with a cutoff at heel. As only the boundary is discretized, the variable mesh technique is advisable for modeling unconfined seepage analyses. The accuracy, effectiveness, and efficiency of the method are demonstrated by modeling several unconfined seepage flow problems. Copyright © 2011 John Wiley & Sons, Ltd.     

水气二相流方程的一种离散数值解法

提出了水气二相流方程的一种数值解法.在利用有限元方法求解水气二相流方程时,引入了离散Newton迭代方法,用于非线性有限元方程组的线性化处理,将这一步计算的收敛阶由原有研究的线性收敛提高到平方收敛,并避免了直接应用Newton迭代方法给编程带来的不便.同时在求解两相的有限元方程组时,采用两相方程组并行迭代的方法,与联立计算相比节省了大量的内存空间.     

用Gao-Yong湍流方程组数值模拟高雷诺数顶盖驱动方腔流

基于有限元方法中的虚位移原理,写出了不可压流Gao-Yong湍流方程组弱形式的微分方程表达式,利用有限元程序自动生成系统(FEPG)得到了全部有限元计算程序。对二维顶盖驱动方腔流进行了数值模拟,得到了方腔内的涡量分布和中心线上的速度变化曲线,并且成功地模拟出了高雷诺数时方腔内的双涡结构,计算结果与文献提供基准解符合良好。进一步证实了Gao-Yong湍流方程组能够准确模拟复杂湍流粘性场,并且为Gao-Yong湍流方程组的工程应用奠定了基础。     

A fracture mapping and extended finite element scheme for coupled deformation and fluid flow in fractured porous media

This paper presents a fracture mapping (FM) approach combined with the extended finite element method (XFEM) to simulate coupled deformation and fluid flow in fractured porous media. Specifically, the method accurately represents the impact of discrete fractures on flow and deformation, although the individual fractures are not part of the finite element mesh. A key feature of FM‐XFEM is its ability to model discontinuities in the domain independently of the computational mesh. The proposed FM approach is a continuum‐based approach that is used to model the flow interaction between the porous matrix and existing fractures via a transfer function. Fracture geometry is defined using the level set method. Therefore, in contrast to the discrete fracture flow model, the fracture representation is not meshed along with the computational domain. Consequently, the method is able to determine the influence of fractures on fluid flow within a fractured domain without the complexity of meshing the fractures within the domain. The XFEM component of the scheme addresses the discontinuous displacement field within elements that are intersected by existing fractures. In XFEM, enrichment functions are added to the standard finite element approximation to adequately resolve discontinuous fields within the simulation domain. Numerical tests illustrate the ability of the method to adequately describe the displacement and fluid pressure fields within a fractured domain at significantly less computational expense than explicitly resolving the fracture within the finite element mesh. Copyright © 2013 John Wiley & Sons, Ltd.     

采用复合单元法模拟裂隙多孔介质变饱和流动

采用复合单元法建立了模拟裂隙多孔介质变饱和流动的数值模型。该模型具有以下特点：裂隙不需要离散成特定单元，而是根据几何位置插入到孔隙基质单元中形成复合单元；在复合单元中，分别建立裂隙流和孔隙基质流的计算方程，二者通过裂隙?基质界面产生联系并整合成复合单元方程；复合单元方程具有和常规有限单元方程相同的格式，因此，可以使用常规有限单元方程的求解技术。采用欠松弛迭代、集中质量矩阵以及自适应时步调节等技术，开发了裂隙多孔介质变饱和流动计算程序。通过模拟一维干土入渗和复杂裂隙含水层内的流动问题，验证了该模型的合理性和适用性。模拟结果为进一步认识非饱和裂隙含水层地下水流动特性提供了理论依据。     

二维扩散输移问题的一种新的有限体积算法

从空间离散格式上比较了有限元格式与有限体积格式的异同处,证明有限体积格式是有限元格式的一种特例,有限元格式也是一种守恒型格式。并根据有限元格式与有限体积格式的比较,提出了一种新的高精度、稳定性好的有限体积格式。通过模拟非定常的纯对流方程和二维对流扩散方程,证明了有限体积格式的优点。     

双层介质中潜水井稳定渗流的有限单元分析

均质含水层中的自由水面(潜水)问题,虽然有各种各样的分析解和近似解,但对非均质条件还缺乏合理的解答。目前虽可采用虚构流网的概念相当准确地计算侧向补给情况的潜水井流量,但由于这种方法不能计算渗流场的水头分布,仍不易把它引用到不稳定渗流问题中去。     

Analysis of ventilation shock losses by finite element method

Summary This paper presents an analysis of ventilation shock losses using the two-dimensional finite element method. It outlines the finite element formulation, in which the ventilation air flow is assumed to be a steady-state, viscous, incompressible and fully developed turbulent flow. In addition, the modified version of Van Driest model is assumed for the effective kinematic viscosity.The shock losses were estimated for the case simulating the airway configuration for (a) split, (b) junction and (c) gradual contraction. The results are compared and discussed with the shock loss data available in the ventilation literature. To estimate the shock losses, the flow solution in terms of pressure was first obtained for each airway configuration using the finite element method. The solution was then compared with that of a straight airway for estimating the difference in the head losses.     

城区渗流场中建筑物间缝隙有限元简化分析方法

为解决城区渗流场中缝隙多、区域的大尺度与物体的小尺寸相差了几个数量级,有限元分析上存在的困难,根据缝隙流动的解析解,经适当地简化,将缝隙附近的渗流场分为缝外区和缝内区两部分,将缝外区流动简化为一个半径为半缝宽的井流,缝内区流动简化为一个均匀流。引入附加阻力系数,考虑缝隙出入口附近的由于过水断面突然变化而引的局部水头损失,并得到附加阻力系数的表达式。利用该式,结合渗流场中井点水头的修正公式,得到了能有效模拟渗流场中建筑物间缝隙的修正线单元公式。利用该线单元编制了相关程序模拟城区渗流场中的缝隙,使得缝隙出入口处的网格尺寸为缝隙两边的建筑物边长1/3～1/4时就能得到较为准确的结果,避免了区域性渗流场中小尺寸物体的网格划分问题,通过算例验证了该公式的精度和边界适应性。     

城区渗流场中建筑物间缝隙有限元简化分析方法

为解决城区渗流场中缝隙多、区域的大尺度与物体的小尺寸相差了几个数量级，有限元分析上存在的困难，根据缝隙流动的解析解，经适当地简化，将缝隙附近的渗流场分为缝外区和缝内区两部分，将缝外区流动简化为一个半径为半缝宽的井流，缝内区流动简化为一个均匀流。引入附加阻力系数，考虑缝隙出入口附近的由于过水断面突然变化而引的局部水头损失，并得到附加阻力系数的表达式。利用该式，结合渗流场中井点水头的修正公式，得到了能有效模拟渗流场中建筑物间缝隙的修正线单元公式。利用该线单元编制了相关程序模拟城区渗流场中的缝隙，使得缝隙出入口处的网格尺寸为缝隙两边的建筑物边长1/3～1/4时就能得到较为准确的结果，避免了区域性渗流场中小尺寸物体的网格划分问题，通过算例验证了该公式的精度和边界适应性。     

A fast flux tube-based method for solute-transport simulation

A new method to calculate the transport of dissolved species in aquifers is presented. This approach is an extension of the stream tubes which are used for flow computation. The flux tubes defined here are conservative for solutes, but not for water mass. The flux tubes are first defined in a general domain and then calculated in a two-dimensional uniform flow field. The tubes?? computation is based on a parametric solution. The method is extended further in order to deal with heterogeneous media. A particle-tracking algorithm is used where the deviation of the flux-tube boundaries due to dispersion is accounted for. The approximate solution obtained by this approach is compared to classical numerical solutions given by a finite difference code (RT3D) and a finite element code (FEFLOW). This comparison was performed for several test cases with increasing complexity. The differences between the flux-tube approach and the other methods always remain small, even regarding mass conservation. The major advantage of the flux-tube approach is the ability to reach a solution quickly, as the method is hundreds to thousands of times faster than classical finite difference or finite element models.     

A Comparative Study of Numercial Solution of PDE inGeosciences

A number of phenomena and processes in geosciences can be summarized by second order partial differential equations. The major numerical methods for their solution include the classical finite difference method and the finite element method newly developed in the last two or three decades. Since 1977 the author has proved that for the Laplace and Poisson equations, these two methods are identical and are different only in the process of formulation. For transient problems, such as heat conduction in the earth and the groundwater and oil-gas unsteady flow in porous media, there are some differences in resulting linear algebraic euqations. In general, two methods give similar results, but when the time step is decreased to some extent, the resulting algebraic equation will be consistent with the anti-heat conduction equation rather than the original heat conduction equation. This is the reason why unrealistic potentials are produced by the finite element method. Such a problem can be overcome by using the     

A hybrid approach to flow net generation

By following a simple set of rules, a flow net can be manually constructed to obtain a graphical solution to the Laplace equation for simple two‐dimensional (2‐D) flow problems. With the advent of numerical solutions such as the finite difference and finite element methods, it is more common to generate a flow net automatically using the nodal head and flow values output by the computer program. Two methods have been published for automatically generating flow nets from finite element solutions: the stream‐function method and the particle‐tracking method. The stream‐function method works well for many cases, but it does not work for problems with holes in the mesh or internal sources or sinks. The particle‐tracking method works for all cases, but previously published algorithms that utilize this method do not result in the properly sized flow channels. A new approach is presented in this paper that is a hybrid of the stream‐function and particle‐tracking approaches. This method works for all cases and generates the properly sized flow channels. Copyright © 2001 John Wiley & Sons, Ltd.     

交叉裂隙水流N-S方程有限元分析

运用有限元理论从N-S方程(Navier-Stokesequation)出发,研究了交叉裂隙水流的基本特征.通过有限元计算结果与传统网络水力学算法及实验成果的对比分析,论述了在裂隙渗流分析中,传统网络计算方法的近似程度完全能为工程要求所接受.网络裂隙渗流分析计算时,不必再单独考虑裂隙水的偏流问题.     

无压渗流问题分析的多节点有限元方法

采用基于平均值插值的多节点有限元方法分析有自由面渗流问题。在自由面附近采用多节点单元逼近自由面,利用平均值插值建立多节点单元的形函数,在远离自由面的区域采用四边形/三角形单元剖分。给定一个初始自由面位置,通过对渗流控制方程的多节点有限元求解,根据自由面上节点水头值判断自由面节点的调整方向和大小,最终迭代求出自由面的位置。土坝渗流问题的数值分析表明了所提方法的有效性和足够的计算精度。     

Numerical modelling of large deformations associated with driving of open‐ended piles

This paper presents a ‘Eulerian‐like’ finite element technique to simulate the large accumulated displacements of piles subjected to multiple hammer blows. For each hammer blow, results are obtained using a standard small strain finite element model and, at the end of each hammer blow, material flow is taken into account with reference to a fixed finite element mesh. Residual stresses calculated at the Gauss integration points of the deformed finite element mesh are mapped on to the fixed finite element mesh, and these stresses are used as initial stresses for the next hammer blow. At the end of each hammer blow, stiffness and mass matrices are recalculated for the volume of material remaining inside the fixed finite element mesh. Results obtained with and without allowing material to flow through the fixed mesh are compared for several hammer blows. Build up of residual stresses, soil flow and yielded points around the pile are presented for plugged, partially‐plugged and unplugged piles. Using the new finite element technique, the driving of a pile from the soil surface is studied. The ability to analyse this and other large deformation problems is the main advantage of the new finite element technique. Copyright © 2000 John Wiley & Sons, Ltd.     

非稳定温度场的三维p型有限元方法研究

非稳定温度场的边界往往是温度梯度较大的区域,需要划分较密的网格。但在三维温度应力的计算中太密的网格往往使得计算成本很高。将p型阶谱有限元方法引入温度场的计算,则可以开拓一条解决这类问题的有效途径。研究了p型有限元中的温度边界条件和初值,编制了有关程序,对同一个算例,采用不同网格密度的常规有限元法和p型有限元法进行了温度场计算分析和对比。研究结果表明,用p型有限元进行非稳定温度场计算,在网格稀疏的条件下能够达到较高的精度。     

基于塑性极限分析上限法理论的土质边坡可靠度分析

将极限分析的上限定理、有限元离散思想、随机规划理论和蒙特卡洛方法这四者结合起来,提出了一种土质边坡可靠度分析的上限数值方法。首先采用三节点有限单元离散土质边坡,然后将土体的抗剪参数设为随机变量,根据上限定理构建同时满足三角形单元的塑性流动约束条件、单元公共边的塑性流动约束条件和单元速度边界条件的机动许可速度场,并根据内功功率等于外功功率条件建立目标函数,构建土质边坡可靠度分析的上限法随机规划模型。采用蒙特卡洛方法求解上限法随机规划模型,同时提出了一种基于上限法速度场的边坡失效风险系数估算方法,该方法特别适用于具有多种失效模式的边坡风险分析。对2个经典算例进行了深入分析,验证了方法的正确性。