共查询到20条相似文献,搜索用时 46 毫秒
1.
We present a new second-order accurate monotone finite volume (FV) method for the steady-state advection–diffusion equation. The method uses a nonlinear approximation for both diffusive and advective fluxes and guarantees solution non-negativity. The interpolation-free approximation of the diffusive flux uses the nonlinear two-point stencil proposed in Lipnikov [23]. Approximation of the advective flux is based on the second-order upwind method with a specially designed minimal nonlinear correction. The second-order convergence rate and monotonicity are verified with numerical experiments. 相似文献
2.
3.
Qian-Yong Chen Jing Wan Yahan Yang Rick T. Mifflin 《Journal of computational physics》2008,227(3):1701-1721
It is well known that the two-point flux approximation, a numerical scheme used in most commercial reservoir simulators, has O(1) error when grids are not K-orthogonal. In the last decade, the multi-point flux approximations have been developed as a remedy. However, non-physical oscillations can appear when the anisotropy is really strong. We found out the oscillations are closely related to the poor approximation of pressure gradient in the flux computation.In this paper, we propose the control volume enriched multi-point flux approximation (EMPFA) for general diffusion problems on polygonal and polyhedral meshes. Non-physical oscillations are not observed for realistic and strongly anisotropic heterogeneous material properties described by a full tensor. Exact linear solutions are recovered for grids with non-planar interfaces, and a first and second order convergence are achieved for the flux and scalar unknowns, respectively. 相似文献
4.
5.
In this paper, we describe a second-order accurate cell-centered finite volume method for solving anisotropic diffusion on two-dimensional unstructured grids. The resulting numerical scheme, named CCLAD (Cell-Centered LAgrangian Diffusion), is characterized by a local stencil and cell-centered unknowns. It is devoted to the resolution of diffusion equation on distorted grids in the context of Lagrangian hydrodynamics wherein a strong coupling occurs between gas dynamics and diffusion. The space discretization relies on the introduction of two half-edge normal fluxes and two half-edge temperatures per cell interface using the partition of each cell into sub-cells. For each cell, the two half-edge normal fluxes attached to a node are expressed in terms of the half-edge temperatures impinging at this node and the cell-centered temperature. This local flux approximation can be derived through the use of either a sub-cell variational formulation or a finite difference approximation, leading to the two variants CCLADS and CCLADNS. The elimination of the half-edge temperatures is performed locally at each node by solving a small linear system which is obtained by enforcing the continuity condition of the normal heat flux across sub-cell interface impinging at the node. The accuracy and the robustness of the present scheme is assessed by means of various numerical test cases. 相似文献
6.
F. Hermeline 《Journal of computational physics》2000,160(2):844
A new finite volume method is presented for discretizing general linear or nonlinear elliptic second-order partial-differential equations with mixed boundary conditions. The advantage of this method is that arbitrary distorted meshes can be used without the numerical results being altered. The resulting algorithm has more unknowns than standard methods like finite difference or finite element methods. However, the matrices that need to be inverted are positive definite, so the most powerful linear solvers can be applied. The method has been tested on a few elliptic and parabolic equations, either linear, as in the case of the standard heat diffusion equation, or nonlinear, as in the case of the radiation diffusion equation and the resistive diffusion equation with Hall term. 相似文献
7.
In this paper, we propose a new positivity-preserving finite volume scheme with fixed stencils for the nonequilibrium radiation diffusion equations on distorted meshes. This scheme is used to simulate the equations on meshes with both the cell-centered and cell-vertex unknowns. The cell-centered unknowns are the primary unknowns, and the element vertex unknowns are taken as the auxiliary unknowns, which can be calculated by interpolation algorithm. With the nonlinear two-point flux approximation, the interpolation algorithm is not required to be positivity-preserving. Besides, the scheme has a fixed stencil and is locally conservative. The Anderson acceleration is used for the Picard method to solve the nonlinear systems efficiently. Several numerical results are also given to illustrate the efficiency and strong positivity-preserving quality of the scheme. 相似文献
8.
研究二维散乱点集上数值求解非线性扩散方程的有限方向差分方法。利用五个邻点信息构造具有最小模板的离散格式,并且离散系数具有显式表达式。另外,利用五点公式获得了间断问题物质界面的离散格式,该格式对界面流的计算具有近似二阶精度。不同计算区域及不同类型的离散点集上的计算结果验证了方法的有效性。 相似文献
9.
We consider a control volume discretization with a multi-point flux approximation to model Discrete Fracture–Matrix systems for anisotropic and fractured porous media in two and three spatial dimensions. Inspired by a recently introduced approach based on a two-point flux approximation, we explicitly account for the fractures by representing them as hybrid cells between the matrix cells. As well as simplifying the grid generation, our hybrid approach excludes small cells in the intersection of the fractures and hence avoids severe time-step restrictions associated with small cells. Excluding the small cells also reduces the condition number of the discretization matrix. For examples involving realistic anisotropy ratios in the permeability, numerical results show significant improvement compared to existing methods based on two-point flux approximations. We also investigate the hybrid method by studying the convergence rates for different apertures and fracture/matrix permeability ratios. Finally, the effect of removing the cells in the intersections of the fractures are studied. Together, these examples demonstrate the efficiency, flexibility and robustness of our new approach. 相似文献
10.
11.
《Journal of computational physics》2008,227(1):492-512
We consider a non-linear finite volume (FV) scheme for stationary diffusion equation. We prove that the scheme is monotone, i.e. it preserves positivity of analytical solutions on arbitrary triangular meshes for strongly anisotropic and heterogeneous full tensor coefficients. The scheme is extended to regular star-shaped polygonal meshes and isotropic heterogeneous coefficients. 相似文献
12.
《Journal of computational physics》2008,227(2):1046-1074
This paper presents a class of kernel-free boundary integral (KFBI) methods for general elliptic boundary value problems (BVPs). The boundary integral equations reformulated from the BVPs are solved iteratively with the GMRES method. During the iteration, the boundary and volume integrals involving Green’s functions are approximated by structured grid-based numerical solutions, which avoids the need to know the analytical expressions of Green’s functions. The KFBI method assumes that the larger regular domain, which embeds the original complex domain, can be easily partitioned into a hierarchy of structured grids so that fast elliptic solvers such as the fast Fourier transform (FFT) based Poisson/Helmholtz solvers or those based on geometric multigrid iterations are applicable. The structured grid-based solutions are obtained with standard finite difference method (FDM) or finite element method (FEM), where the right hand side of the resulting linear system is appropriately modified at irregular grid nodes to recover the formal accuracy of the underlying numerical scheme. Numerical results demonstrating the efficiency and accuracy of the KFBI methods are presented. It is observed that the number of GMRES iterations used by the method for solving isotropic and moderately anisotropic BVPs is independent of the sizes of the grids that are employed to approximate the boundary and volume integrals. With the standard second-order FEMs and FDMs, the KFBI method shows a second-order convergence rate in accuracy for all of the tested Dirichlet/Neumann BVPs when the anisotropy of the diffusion tensor is not too strong. 相似文献
13.
A non-linear finite volume method with monotone matrix for the diffusion equation is presented. It does not extrapolate the primary variable to Neumann boundaries, as this was previously done in similar methods. This change results in faster convergence. Computation time is significantly shortened further using the reduced rank extrapolation method (RRE), and imposing an upper limit on the number of linear iterations per non-linear step. Second-order accuracy and performance improvement are demonstrated by numerical examples. 相似文献
14.
Annamaria Mazzia Gianmarco Manzini Mario Putti 《Journal of computational physics》2011,230(23):8410-8426
We study the performance of Godunov mixed methods, which combine a mixed-hybrid finite element solver and a Godunov-like shock-capturing solver, for the numerical treatment of the advection–dispersion equation with strong anisotropic tensor coefficients. It turns out that a mesh locking phenomenon may cause ill-conditioning and reduce the accuracy of the numerical approximation especially on coarse meshes. This problem may be partially alleviated by substituting the mixed-hybrid finite element solver used in the discretization of the dispersive (diffusive) term with a linear Galerkin finite element solver, which does not display such a strong ill conditioning. To illustrate the different mechanisms that come into play, we investigate the spectral properties of such numerical discretizations when applied to a strongly anisotropic diffusive term on a small regular mesh. A thorough comparison of the stiffness matrix eigenvalues reveals that the accuracy loss of the Godunov mixed method is a structural feature of the mixed-hybrid method. In fact, the varied response of the two methods is due to the different way the smallest and largest eigenvalues of the dispersion (diffusion) tensor influence the diagonal and off-diagonal terms of the final stiffness matrix. One and two dimensional test cases support our findings. 相似文献
15.
16.
This paper presents a class of kernel-free boundary integral (KFBI) methods for general elliptic boundary value problems (BVPs). The boundary integral equations reformulated from the BVPs are solved iteratively with the GMRES method. During the iteration, the boundary and volume integrals involving Green’s functions are approximated by structured grid-based numerical solutions, which avoids the need to know the analytical expressions of Green’s functions. The KFBI method assumes that the larger regular domain, which embeds the original complex domain, can be easily partitioned into a hierarchy of structured grids so that fast elliptic solvers such as the fast Fourier transform (FFT) based Poisson/Helmholtz solvers or those based on geometric multigrid iterations are applicable. The structured grid-based solutions are obtained with standard finite difference method (FDM) or finite element method (FEM), where the right hand side of the resulting linear system is appropriately modified at irregular grid nodes to recover the formal accuracy of the underlying numerical scheme. Numerical results demonstrating the efficiency and accuracy of the KFBI methods are presented. It is observed that the number of GMRES iterations used by the method for solving isotropic and moderately anisotropic BVPs is independent of the sizes of the grids that are employed to approximate the boundary and volume integrals. With the standard second-order FEMs and FDMs, the KFBI method shows a second-order convergence rate in accuracy for all of the tested Dirichlet/Neumann BVPs when the anisotropy of the diffusion tensor is not too strong. 相似文献
17.
This paper focuses on flux-continuous pressure equation approximation for strongly anisotropic media. Previous work on families of flux-continuous schemes for solving the general geometry–permeability tensor pressure equation has focused on single-parameter families. These schemes have been shown to remove the O(1) errors introduced by standard two-point flux reservoir simulation schemes when applied to full-tensor flow approximation. Improved convergence of the schemes has also been established for specific quadrature points. However these schemes have conditional M-matrices depending on the strength of the off-diagonal tensor coefficients. When applied to cases involving full-tensors arising from strongly anisotropic media, the point-wise continuous schemes can fail to satisfy the maximum principle and induce severe spurious oscillations in the numerical pressure solution.New double-family flux-continuous locally conservative schemes are presented for the general geometry–permeability tensor pressure equation. The new double-family formulation is shown to expand on the current single-parameter range of existing conditional M-matrix schemes. The conditional M-matrix bounds on a double-family formulation are identified for both quadrilateral and triangle cell grids. A quasi-positive QM-matrix analysis is presented that classifies the behaviour of the new schemes with respect to double-family quadrature in regions beyond the M-matrix bounds. The extension to double-family quadrature is shown to be beneficial, resulting in novel optimal anisotropic quadrature schemes. The new methods are applied to strongly anisotropic full-tensor field problems and yield results with sharp resolution, with only minor or practically zero spurious oscillations. 相似文献
18.
Jean-Luc Guermond Richard Pasquetti Bojan Popov 《Journal of computational physics》2011,230(11):4248-4267
A new class of high-order numerical methods for approximating nonlinear conservation laws is described (entropy viscosity method). The novelty is that a nonlinear viscosity based on the local size of an entropy production is added to the numerical discretization at hand. This new approach does not use any flux or slope limiters, applies to equations or systems supplemented with one or more entropy inequalities and does not depend on the mesh type and polynomial approximation. Various benchmark problems are solved with finite elements, spectral elements and Fourier series to illustrate the capability of the proposed method. 相似文献
19.
Heterogeneous anisotropic diffusion problems arise in the various areas of science and engineering including plasma physics, petroleum engineering, and image processing. Standard numerical methods can produce spurious oscillations when they are used to solve those problems. A common approach to avoid this difficulty is to design a proper numerical scheme and/or a proper mesh so that the numerical solution validates the discrete counterpart (DMP) of the maximum principle satisfied by the continuous solution. A well known mesh condition for the DMP satisfaction by the linear finite element solution of isotropic diffusion problems is the non-obtuse angle condition that requires the dihedral angles of mesh elements to be non-obtuse. In this paper, a generalization of the condition, the so-called anisotropic non-obtuse angle condition, is developed for the finite element solution of heterogeneous anisotropic diffusion problems. The new condition is essentially the same as the existing one except that the dihedral angles are now measured in a metric depending on the diffusion matrix of the underlying problem. Several variants of the new condition are obtained. Based on one of them, two metric tensors for use in anisotropic mesh generation are developed to account for DMP satisfaction and the combination of DMP satisfaction and mesh adaptivity. Numerical examples are given to demonstrate the features of the linear finite element method for anisotropic meshes generated with the metric tensors. 相似文献
20.
The initial value problem solution of the nonlinear shallow water-wave equations is developed under initial waveforms with and without velocity. We present a solution method based on a hodograph-type transformation to reduce the nonlinear shallow water-wave equations into a second-order linear partial differential equation and we solve its initial value problem. The proposed solution method overcomes earlier limitation of small waveheights when the initial velocity is nonzero, and the definition of the initial conditions in the physical and transform spaces is consistent. Our solution not only allows for evaluation of differences in predictions when specifying an exact initial velocity based on nonlinear theory and its linear approximation, which has been controversial in geophysical practice, but also helps clarify the differences in runup observed during the 2004 and 2005 Sumatran tsunamigenic earthquakes. 相似文献