共查询到20条相似文献，搜索用时 31 毫秒
1.
针对溃坝水流数值模拟面临的复杂地形和不规则边界等问题,基于结构网格建立了适应复杂地形和不规则边界的溃坝水流数值模拟有限体积模型（HydroM2D）。模型基于具有守恒特性的二维浅水方程,利用HLLC格式的近似Riemann解计算网格界面通量,利用MUSCL-Hancock法不断向前积分,使模型在时空上具有二阶精度;对源项进行离散处理确保模型的稳定性;模型引入有效干湿边界和不规则地形边界处理方法,准确模拟了干湿单元的动态交替和复杂边界上的水流特性。最后分别利用水槽试验、物理模型和实际算例对模型进行验证。结果表明,该模型对不同情景下的溃坝洪水模拟结果和实测资料以及现有模型模拟结果具有较高的一致性,模拟精度较高,稳定性较好,具有推广应用价值。 相似文献
2.
We present a new nonlinear monotone finite volume method for diffusion equation and its application to two-phase flow model. We consider full anisotropic discontinuous diffusion or permeability tensors on conformal polyhedral meshes. The approximation of the diffusive flux uses the nonlinear two-point stencil which provides the conventional seven-point stencil for the discrete diffusion operator on cubic meshes. We show that the quality of the discrete flux in a reservoir simulator has great effect on the front behavior and the water breakthrough time. We compare two two-point flux approximations (TPFA), the proposed nonlinear TPFA and the conventional linear TPFA, and multipoint flux approximation (MPFA). The new nonlinear scheme has a number of important advantages over the traditional linear discretizations. Compared to the linear TPFA, the nonlinear TPFA demonstrates low sensitivity to grid distortions and provides appropriate approximation in case of full anisotropic permeability tensor. For nonorthogonal grids or full anisotropic permeability tensors, the conventional linear TPFA provides no approximation, while the nonlinear flux is still first-order accurate. The computational work for the new method is higher than the one for the conventional TPFA, yet it is rather competitive. Compared to MPFA, the new scheme provides sparser algebraic systems and thus is less computational expensive. Moreover, it is monotone which means that the discrete solution preserves the nonnegativity of the differential solution. 相似文献
3.
Eirik Keilegavlen Jeremy E. Kozdon Bradley T. Mallison 《Computational Geosciences》2012,16(4):1021-1042
The governing equations for multiphase flow in porous media have a mixed character, with both nearly elliptic and nearly hyperbolic variables. The flux for each phase can be decomposed into two parts: (1) a geometry- and rock-dependent term that resembles a single-phase flux; and (2) a mobility term representing fluid properties and rock–fluid interactions. The first term is commonly discretized by two- or multipoint flux approximations (TPFA and MPFA, respectively). The mobility is usually treated with single-point upstream weighting (SPU), also known as dimensional or donor cell upstream weighting. It is well known that when simulating processes with adverse mobility ratios, SPU suffers from grid orientation effects. An important example of this, which will be considered in this work, is the displacement of a heavy oil by water. For these adverse mobility ratio flows, the governing equations are unstable at the modeling scale, rendering a challenging numerical problem. These challenges must be addressed in order to avoid systematic biasing of simulation results. In this work, we present a framework for multidimensional upstream weighting for multiphase flow with buoyancy on general two-dimensional grids. The methodology is based on a dual grid, and the resulting transport methods are provably monotone. The multidimensional transport methods are coupled with MPFA methods to solve the pressure equation. Both explicit and fully implicit approaches are considered for time integration of the transport equations. The results show considerable reduction of grid orientation effects compared to SPU, and the explicit multidimensional approach allows larger time steps. For the implicit method, the total number of non-linear iterations is also reduced when multidimensional upstream weighting is used. 相似文献
4.
The trend toward unstructured grids in subsurface flow modeling has prompted interest in the issue of streamline or pathline tracing on unstructured grids. Streamline tracing on unstructured grids is problematic because a continuous velocity field is required for the calculation, while numerical solutions to the groundwater flow equations provide velocity in discretized form only. A method for calculating flow streamlines or pathlines from a finite-volume flow solution is presented. The method uses an unconstrained least squares method on interior cells and a constrained least squares method on boundary cells to approximate cell-centered velocities, which can then be continuously interpolated to any point in the domain of interest. Two-dimensional tests demonstrate that the method correctly reproduces uniform and corner-to-corner flow on fully unstructured grids. In three dimensions using regular hexahedral grids, the method agrees well with established semianalytical methods. Tests also demonstrate that the method produces physically realistic results on fully unstructured three-dimensional grids. 相似文献
5.
Vasiliy Kramarenko Kirill Nikitin Yuri Vassilevski 《Computational Geosciences》2017,21(5-6):1023-1033
We present the latest enhancement of the nonlinear monotone finite volume method for the near-well regions. The original nonlinear method is applicable for diffusion, advection-diffusion, and multiphase flow model equations with full anisotropic discontinuous permeability tensors on conformal polyhedral meshes. The approximation of the diffusive flux uses the nonlinear two-point stencil which reduces to the conventional two-point flux approximation (TPFA) on cubic meshes but has much better accuracy for the general case of non-orthogonal grids and anisotropic media. The latest modification of the nonlinear method takes into account the nonlinear (e.g., logarithmic) singularity of the pressure in the near-well region and introduces a correction to improve accuracy of the pressure and the flux calculation. In this paper, we consider a linear version of the nonlinear method waiving its monotonicity for sake of better accuracy. The new method is generalized for anisotropic media, polyhedral grids and nontrivial cases such as slanted, partially perforated wells or wells shifted from the cell center. Numerical experiments show noticeable reduction of numerical errors compared to the original monotone nonlinear FV scheme with the conventional Peaceman well model or with the given analytical well rate. 相似文献
6.
7.
This paper presents a finite-volume method for hexahedral multiblock grids to calculate multiphase flow in geologically complex reservoirs. Accommodating complex geologic and geometric features in a reservoir model (e.g., faults) entails non-orthogonal and/or unstructured grids in place of conventional (globally structured) Cartesian grids. To obtain flexibility in gridding as well as efficient flow computation, we use hexahedral multiblock grids. These grids are locally structured, but globally unstructured. One major advantage of these grids over fully unstructured tetrahedral grids is that most numerical methods developed for structured grids can be directly used for dealing with the local problems. We present several challenging examples, generated via a commercially available tool, that demonstrate the capabilities of hexahedral multiblock gridding. Grid quality is discussed in terms of uniformity and orthogonality. The presence of non-orthogonal grid and full permeability tensors requires the use of multi-point discretization methods. A flux-continuous finite-difference (FCFD) scheme, previously developed for stratigraphic hexahedral grid with full-tensor permeability, is employed for numerical flow computation. We extend the FCFD scheme to handle exceptional configurations (i.e. three- or five-cell connections as opposed to the regular four), which result from employing multiblock gridding of certain complex objects. In order to perform flow simulation efficiently, we employ a two-level preconditioner for solving the linear equations that results from the wide stencil of the FCFD scheme. The individual block, composed of cells that form a structured grid, serves as the local level; the higher level operates on the global block configuration (i.e. unstructured component). The implementation uses an efficient data structure where each block is wrapped with a layer of neighboring cells. We also examine splitting techniques [14] for the linear systems associated with the wide stencils of our FCFD operator. We present three numerical examples that demonstrate the method: (1) a pinchout, (2) a faulted reservoir model with internal surfaces and (3) a real reservoir model with multiple faults and internal surfaces. 相似文献
8.
This paper extends the multipoint flux-approximation (MPFA) control-volume method to quadrilateral grids for which the adjacent cells do not necessarily share corners. Examples are grids with faults and locally refined grids. This paper gives a derivation of the method for such grids. The difference between two-point flux-approximation (TPFA) results and MPFA results for faults and local grid refinements is demonstrated for synthetic problems. Further, the results are compared with results from uniform fine-grid simulations. The effect of repeated fault patterns as well as anisotropy is investigated. Large errors may be found for the TPFA method for flow through a series of faults in an anisotropic medium. Finally, a comparison is done for a reservoir field application. 相似文献
9.
In this article we present a series of tests to study how well suited the TPFA coefficient matrix is as a preconditioner for the MPFA discrete system of equations in an iterative solver, using a flux splitting method. These tests have been conducted for single-phase flow for a wide range of anisotropy, heterogeneity, and grid skewness (mainly parallelogram grids). We use the K-orthogonal part of the MPFA transmissibilities for a parallelogram grid to govern the TPFA transmissibilities. The convergence of the flux splitting method is for each test case measured by the spectral radius of the iteration matrix. 相似文献
10.
Numerically modeling groundwater flow on finely discretized two- and three-dimensional domains requires solution algorithms appropriate for distributed memory multiprocessor architectures. Multilevel and domain decomposition algorithms are appropriate for preconditioning or solving linear systems in parallel and have, therefore, been applied to linear models for saturated groundwater flow. These algorithms have also been incorporated into more complex nonlinear multiphase flow models in the context of a linearization procedure such as Newton's method. In this work, we study a class of parallel preconditioners based on two-level Schwarz domain decomposition applied in a nonlinear two-phase flow numerical model. The restriction and interpolation operators are based on an aggregation approach that has a straightforward implementation for a variety of applications arising in subsurface modeling: structured and unstructured discretizations, finite elements and finite differences, and multicomponent model equations. We present model formulations, results from numerical experiments, and a comparison of a standard one-level Schwarz method to three two-level aggregation-based methods. 相似文献
11.
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. 相似文献
12.
A Boundary Element based Discontinuous Deformation Analysis (BE‐DDA) method is developed by implementing the improved dual reciprocity boundary element method into the open close iterations based DDA. This newly developed BE‐DDA is capable of simulating both the deformation and movement of blocks in a blocky system. Based on geometry updating, it adopts an incremental dynamic formulation taking into consideration initial stresses and dealing with external concentrated and contact forces conveniently. The boundaries of each block in the discrete blocky system are discretized with boundary elements while the domain of each block is divided into internal cells only for the integration of the domain integral of the initial stress term. The contact forces among blocks are treated as concentrated forces and the open–close iterations are applied to ensure the computational accuracy of block interactions. In the current method, an implicit time integration scheme is adopted for numerical stability. Three examples are used to show the effectiveness of the algorithm in simulating block movement, sliding, deformation and interaction of blocks. At last, block toppling and tunnel stability examples are conducted to demonstrate that the BE‐DDA is applicable for simulation of blocky systems. Copyright © 2016 John Wiley & Sons, Ltd. 相似文献
13.
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.
相似文献
14.
Multiscale methods can in many cases be viewed as special types of domain decomposition preconditioners. The localisation
approximations introduced within the multiscale framework are dependent upon both the heterogeneity of the reservoir and the
structure of the computational grid. While previous works on multiscale control volume methods have focused on heterogeneous
elliptic problems on regular Cartesian grids, we have tested the multiscale control volume formulations on two-dimensional
elliptic problems involving heterogeneous media and irregular grid structures. Our study shows that the tangential flow approximation
commonly used within multiscale methods is not suited for problems involving rough grids. We present a more robust mass conservative
domain decomposition preconditioner for simulating flow in heterogeneous porous media on general grids. 相似文献
15.
将追踪自由表面的流体体积(VOF)法应用于曲线坐标系下水流控制方程的求解中,计入流线弯曲对水流紊动特性的影响,建立了垂向二维强紊动水流的曲率修正的紊流模型,并对溢流坝反弧段的紊流特性进行了数值模拟。数值计算时,采用有限体积法离散水流的控制方程;物理变量,如:压力P、紊动参量κ、ε、γt等,采用交错方式排列(交错网格布置),用SIMPLEC算法求解离散方程。计算结果表明,得到的溢流坝反弧段的自由表面位置、速度场、压力场、剪应力分布和紊动能分布与实验结果吻合良好。 相似文献
16.
岩土体的渗透破坏、地下工程的防渗设计等无不与渗流计算有关。针对渗流自由面问题,提出一种重心拉格朗日插值的配点型无网格方法。由于渗流自由面问题的求解区域是不规则区域,该方法通过将不规则求解区域嵌入一个正则矩形区域,在正则区域上采用重心拉格朗日插值近似未知函数,利用配点法离散渗流问题的控制方程,将重心拉格朗日插值的微分矩阵离散成代数方程表达的矩阵形式。将自由面上的边界条件通过重心拉格朗日插值离散,通过置换方程法和附加方程法施加边界条件,利用正则区域上的重心插值配点法,通过迭代确定最终自由面的位置。数值算例表明所提出的无网格方法对于求解渗流自由面问题的正确性和高精度。 相似文献
17.
H. Hajinezhad Ali R. Soheili Mohammad R. Rasaei F. Toutounian 《Computational Geosciences》2018,22(6):1433-1444
In this paper, the numerical methods for solving the problem of steam injection in the heavy oil reservoirs are presented. We consider a 3-dimensional model of 3-phase flow, oil, water, and steam, with the effect of 3-phase relative permeability. Interphase mass transfer of water and steam is considered; oil is assumed nonvolatile. We apply the simultaneous solution approach to solve the corresponding nonlinear discretized partial differential equation in the fully implicit form. The convergence of finite difference scheme is proved by the Rosinger theorem. The heuristic Jacobian-Free-Newton-Krylov (HJFNK) method is proposed for solving the system of algebraic equations. The result of this proposed numerical method is well compared with some experimental results. Our numerical results show that the first iteration of the full approximation scheme (FAS) provides a good initial guess for the Newton method. Therefore, we propose a new hybrid-FAS-HJFNK method while there is no steam in the reservoir. The numerical results show that the hybrid-FAS-HJFNK method converges faster than the HJFNK method. 相似文献
18.
Katerina Georgiou John Harte Ali Mesbah William J. Riley 《Computational Geosciences》2018,22(3):851-865
We present a numerical method for solving a class of systems of partial differential equations (PDEs) that arises in modeling environmental processes undergoing advection and biogeochemical reactions. The salient feature of these PDEs is that all partial derivatives appear in linear expressions. As a result, the system can be viewed as a set of ordinary differential equations (ODEs), albeit each one along a different characteristic. The method then consists of alternating between equations and integrating each one step-wise along its own characteristic, thus creating a customized grid on which solutions are computed. Since the solutions of such PDEs are generally smoother along their characteristics, the method offers the potential of using larger time steps while maintaining accuracy and reducing numerical dispersion. The advantages in efficiency and accuracy of the proposed method are demonstrated in two illustrative examples that simulate depth-resolved reactive transport and soil carbon cycling. 相似文献
19.
An iterative method is presented for solving a fully coupled and implicit formulation of fluid flow in a porous medium. The mathematical model describes a set of fully coupled three-phase flow of compressible and immiscible fluids in a saturated oil reservoir. The finite element method is applied to obtain the simultaneous solution (SS) for the resulting highly non-linear partial differential equations where fluid pressures are the primary unknowns. The final discretized equations are solved iteratively by using a fully implicit numerical scheme. Several examples, illustrating the use of the present model, are described. The increased stability achieved with this scheme has permitted the use of larger time steps with smaller material balance errors. 相似文献
20.
This paper presents a numerical model for simulating free surface flow in porous media with spatially varying porosity. The governing equations are based on the mixture theory. The resistance forces between solid and fluid is assumed to be nonlinear. A multiphase SPH approach is presented to solve the governing equations. In the multiphase SPH, water is modeled as a weakly compressible fluid, and solid phase is discretized by fixed solid particles carrying information of porosity. The model is validated by several numerical examples including seepage through specimen, fast flow through rockfill dam and wave interaction with porous structure. Good agreements between numerical results and experimental data are obtained in terms of flow rate and evolution of free surface. Parameter study shows that (1) the nonlinear resistance law provides more accurate results; (2) particle size and porosity have significant influence on the porous flow. 相似文献