共查询到20条相似文献，搜索用时 46 毫秒
1.
A time-splitting approach for advection–dispersion equations is considered. The dispersive and advective fluxes are split into two separate partial differential equations (PDEs), one containing the dispersive term and the other one the advective term. On triangular elements a triangle-based high resolution Finite Volume (FV) scheme for advection is combined with a Mixed Hybrid Finite Element (MHFE) technique to solve dispersion. This approach introduces an error proportional to the time step and the overall scheme is only first order accurate if special care is not taken in the definition of the numerical flux approximation for advection. By incorporating the diffusive effects into the definition of this numerical flux, near second order accuracy (up to a logh factor) can be proved theoretically and validated by numerical experiments in both one- and two-dimensional cases. 相似文献
2.
Rigoberto G. Sanabria Castro Sandra M. C. Malta Abimael F. D. Loula Luiz Landau 《Computational Geosciences》2002,5(4):301-330
A finite element formulation is proposed to approximate a nonlinear system of partial differential equations, composed by an elliptic subsystem for the pressure–velocity and a transport equation (convection–diffusion) for the concentration, which models the incompressible miscible displacement of one fluid by another in a rigid porous media. The pressure is approximated by the classical Galerkin method and the velocity is calculated by a post-processing technique. Then, the concentration is obtained by a Galerkin/least-squares space–time (GLS/ST) finite element method. A numerical analysis is developed for the concentration approximation. Then, stability, convergence and numerical results are presented confirming the a priori error estimates. 相似文献
3.
This paper is devoted to the numerical reliability and time requirements of the Mixed Finite Element (MFE) and Mixed-Hybrid Finite Element (MHFE) methods. The behavior of these methods is investigated under the influence of two factors: the mesh discretization and the medium heterogeneity. We show that, unlike the MFE, the MHFE suffers with the presence of badly shaped discretized elements. Thereat, a numerical reliability analyzing software (Aquarels) is used to detect the instability of a matrix-inversion code generated automatically by a symbolic manipulator. We also show that the spectral condition number of the algebraic systems furnished by both methods in heterogeneous media grows up linearly according to the smoothness of the hydraulic conductivity. Furthermore, it is found that the MHFE could accumulate numerical errors if large jumps in the tensor of conductivity take place. Finally, we compare running-times for both algorithms by giving various numerical experiments. 相似文献
4.
This paper is concerned with numerical methods for the modeling of flow and transport of contaminant in porous media. The
numerical methods feature the mixed finite element method over triangles as a solver to the Darcy flow equation and a conservative
finite volume scheme for the concentration equation. The convective term is approximated with a Godunov scheme over the dual
finite volume mesh, whereas the diffusion–dispersion term is discretized by piecewise linear conforming triangular finite
elements. It is shown that the scheme satisfies a discrete maximum principle. Numerical examples demonstrate the effectiveness
of the methodology for a coupled system that includes an elliptic equation and a diffusion–convection–reaction equation arising
when modeling flow and transport in heterogeneous porous media. The proposed scheme is robust, conservative, efficient, and
stable, as confirmed by numerical simulations.
相似文献
5.
浅水湖泊风生流的迎风有限元数值模型研究 总被引：7，自引：1，他引：6
从控制方程组出发,以Galerkin有限元法为基础,引入沿流线加权的权函数(沿流线加权的迎风有限元法)和选择性集中系数矩阵方法,推导、建立了一适合于浅水湖泊风生流计算的二维迎风有限元数值模型。并以太湖为例,对模型作了检验,分析了均匀、定常风场持续作用下太湖风生流场的形态特征。 相似文献
6.
On 28 February 1969, the coasts of Portugal, Spain and Morocco were affected by sea waves generated by a submarine earthquake (Ms = 7.3) with its epicenter located off Portugal. The propagation of this tsunami has been simulated by a finite element numerical model solving the Boussinesq equations. These equations have been discretized using the finite element Galerkin method and a Crank–Nicholson scheme in time. The model is validated by investigating the propagation of a solitary wave over a flat bottom. The grid sizes for the 1969 event have been determined by one-dimensional tests offshore and in shallow water regions. The two-dimensional simulation of the 1969 tsunami is carried out using the hydraulic source calculated from the geophysical model of Okada and seismic parameters of Fukao. The modeled waves are compared with the recorded ones with respect to travel times, maximum amplitudes and periods of the signal. The comparison between Boussinesq and shallow-water models shows that the effects of frequency dispersion are minor. Good agreement is found for most of the studied gauges. 相似文献
7.
David Trujillo 《Computational Geosciences》2004,8(2):173-185
A mixed finite volume method is applied in order to simulate a simplified Far Field model. We show how this method can be adapted to both the elliptic and the convection–diffusion equations describing the water flow and the variations of the concentration of the radioactive element respectively. For the water flow problem, we justify an a posteriori error estimator. Finally, we present some numerical results. 相似文献
8.
9.
In this paper, we present a numerical model for simulating two-phase (oil–water and air–water) incompressible and immiscible flow in porous media. The mathematical model which is based on a fractional flow formulation is formed of two nonlinear partial differential equations: a mean pressure equation and a water saturation equation. These two equations can be solved in a sequential manner. Two numerical methods are used to discretize the equations of the two-phase flow model: mixed hybrid finite elements are used to treat the pressure equation, h-based Richards' equation and the diffusion term in the saturation equation, the advection term in the saturation equation is treated with the discontinuous finite elements. We propose a better way to calculate the nonlinear coefficients contained in our equations on each element of the discretized domain. In heterogeneous porous media, the saturation becomes discontinuous at the interface between two porous media. We show in this paper how to use the capillary pressure–saturation relationship in order to handle the saturation jump in the mixed hybrid finite element method. The two-phase flow simulator is verified against analytical solutions for some flow problems treated by other authors. 相似文献
10.
This paper presents three-dimensional finite element simulations to evaluate diffusion and dispersion tensors in periodic
porous media in the presence of an advective velocity field. These tensors are evaluated in the framework of the double-scale
expansion technique. Two problems, a Newtonian flow and a vector-valued advection–diffusion equation, have to be sequentially
solved at the pore scale. Finite element techniques to approximate these problems are proposed and analyzed. Numerical results
in three-dimensional networks of spheres are presented to quantitatively assess the impact of the pore morphology and of the
advection velocity on the diffusion and dispersion tensors. 相似文献
11.
以重力位在场源内部满足泊松方程为依据,以重力矢量满足第三类边界条件为切入点,推导了与三度体重力矢量满足的边值问题相对应的变分问题,进而利用有限单元法实现了对变分问题的求解.立方体模型试验结果表明:文中提出的新的系数矩阵存储方式较之传统方式能够更有效地节约存储空间,且为利用预条件共轭梯度技术更加快速地求解线性方程组提供了保障;重力矢量的计算精度与边界长度及单元网格的边长息息相关,其计算效率则主要取决于所要计算的节点总数和大型稀疏线性方程组求解算法的优劣;一般情况下,当单元的边长小于场源体边长的1/10、边界长度大于场源体长度的7.5倍时,能够获得理想的结果. 相似文献
12.
基于二次插值的线源可控源有限元数值模拟 总被引：3，自引：0，他引：3
在准静态近似条件下,采用矩形网格单元和双二次函数插值就频率域二维线源边值问题进行了有限元数值模拟。在二维地电条件下,给出了边值问题和变分问题,并通过有限单元法对模型进行单元剖分、插值、积分和整体合成,最后通过求解复系数方程组得到了地表视电阻率响应。引入伪delta函数模拟线源,消弱了源带来的奇异性。通过与均匀大地以及层状介质模型的解析解对比,平均相对误差分别为0.71%和1.12%。建立了两个异常体模型,数值模拟表明异常响应比较明显,为进一步实现三维可控源电磁法有限元数值模拟提供了基础。 相似文献
13.
Numerical identification of diffusion parameters in a nonlinear convection–diffusion equation is studied. This partial differential
equation arises as the saturation equation in the fractional flow formulation of the two-phase porous media flow equations.
The forward problem is discretized with the finite difference method, and the identification problem is formulated as a constrained
minimization problem. We utilize the augmented Lagrangian method and transform the minimization problem into a coupled system
of nonlinear algebraic equations, which is solved efficiently with the nonlinear conjugate gradient method. Numerical experiments
are presented and discussed.
This work was partially supported by the Research Council of Norway (NFR), under grant 128224/431. 相似文献
14.
We give some results obtained for the Couplex test cases proposed by the ANDRA. In this paper our aim is twofold. Firstly, to compute the release of nuclides out of the repository by concentrating on the 3D near field (Couplex 2). The simulation of the transport phenomena takes into account the dissolution of the glass containers and congruent emissions of the radio-nuclides including filiation chains and some simplified chemistry. Secondly, it is to use the near field computations in order to simulate the nuclide migrations in a 2D far field (Couplex 3). Coupling in between the two simulations takes into consideration the periodicity of the disposal modules and the geometry of the repository described in Couplex 1. The mixed finite element and discontinuous Galerkin methods are used to solve the convection–diffusion equations. In order to handle the nonlinear precipitation/dissolution term, we developed a new iterative technique that combines Picard and Newton–Raphson methods. 相似文献
15.
在求解非稳定地下水溶质运移模型时,若对流项占优,则模型表现出双曲方程的特性。针对这种特性,采用非标准Galerkin有限元方法进行求解是解决这类问题的有效途径。分别采用Wavelet-Galer-kin有限元方法、迎风有限元方法和特征有限元方法对强对流溶质运移模型进行了求解,并将其结果与标准Galerkin有限元和解析解进行对比。结果表明:标准Galerkin有限元方法会产生强烈的数值振荡;Wavelet-Galerkin有限元方法的时空定位效果好;迎风有限元方法能够有效降低数值振荡现象,但迎风因子对解的影响较大,而且会带来时间延迟;特征有限元方法能够提高解的精度,故可以认为特征有限元方法是求解强对流地下水溶质运移模型的首选方法。 相似文献
16.
采用复合单元法建立了模拟裂隙多孔介质变饱和流动的数值模型。该模型具有以下特点：裂隙不需要离散成特定单元，而是根据几何位置插入到孔隙基质单元中形成复合单元；在复合单元中，分别建立裂隙流和孔隙基质流的计算方程，二者通过裂隙?基质界面产生联系并整合成复合单元方程；复合单元方程具有和常规有限单元方程相同的格式，因此，可以使用常规有限单元方程的求解技术。采用欠松弛迭代、集中质量矩阵以及自适应时步调节等技术，开发了裂隙多孔介质变饱和流动计算程序。通过模拟一维干土入渗和复杂裂隙含水层内的流动问题，验证了该模型的合理性和适用性。模拟结果为进一步认识非饱和裂隙含水层地下水流动特性提供了理论依据。 相似文献
17.
介绍了垂直钻进系统的工作原理及组成，纠斜机构是由均布在同一圆周上的3个导向块组成，导向块是自动垂钻系统中最主要的执行机构，承受较大的压力。为保证系统设计的可靠性，利用先进的有限元分析方法模拟导向块在工作状态下的受力情况，通过创建有限元模型，施加载荷并求解，分析查看结果等步骤对零件进行仿真分析，较真实地反映出零件在受力状态的应力和应变情况，预测零件的结构强度并分析其设计合理性，使对零部件关键参数的设计更合理，有效解决设计中的可靠性分析和结构优化问题。 相似文献
18.
In this paper, we are interested in the modeling of wave propagation in viscoelastic media. We present a family of models which generalize the Zeners model. We achieve its mathematical analysis: existence and uniqueness of solutions, energy decay and propagation with finite speed. For the numerical resolution, we extend a mixed finite element method proposed in [8]. This method combines mass lumping with a centered explicit scheme for time discretization. For the resulting scheme, we prove a discrete energy decay result and provide a sufficient stability condition. For the numerical simulation in open domains we adapt the perfectly matched layers techniques to viscoelastic waves [23]. Various numerical results are presented. 相似文献
19.
双排桩-锚杆支护的有限元模拟 总被引：2，自引：0，他引：2
基坑支护模型往往在计算假定的基础上建立,而依据计算假定建立的模型往往不能得到和实际相符的结果。但有限元软件可以较为真实的模拟基坑开挖的过程中应力和位移的变化情况,是土工数值分析的重要手段。本文运用大型有限元软件ABAQUS建立双排桩-锚杆结构有限元模型,全面的考虑土体的特性、桩土的共同作用及桩间土对支护结构的影响等因素,分析支护结构在土体开挖荷载作用下的内力和变形,为设计和施工提供了参考。 相似文献
20.
We present the results of a study on a posteriori error control strategies for finite volume element approximations of second order elliptic differential equations. Finite volume methods ensure local mass conservation and, combined with some up-wind strategies, give monotone solutions. We adapt the local refinement techniques known from the finite element method to the finite volume discretizations of various boundary value problems for steady-state convection–diffusion–reaction equations. In this paper we derive and study a residual type error estimator and illustrate its practical performance on a series of computational tests in 2 and 3 dimensions. Our tests show that the discussed locally conservative approximation methods with a posteriori error control can be used successfully in numerical simulation of fluid flow and transport in porous media. 相似文献