共查询到20条相似文献,搜索用时 15 毫秒
1.
High-Order Runge-Kutta Discontinuous Galerkin Methods with a New Type of Multi-Resolution WENO Limiters on Tetrahedral Meshes
下载免费PDF全文
![点击此处可从《Communications In Computational Physics》网站下载免费的PDF全文](/ch/ext_images/free.gif)
In this paper, the second-order and third-order Runge-Kutta discontinuous
Galerkin (RKDG) methods with multi-resolution weighted essentially non-oscillatory
(WENO) limiters are proposed on tetrahedral meshes. The multi-resolution WENO
limiter is an extension of a finite volume multi-resolution WENO scheme developed
in [81], which serves as a limiter for RKDG methods on tetrahedral meshes. This new
WENO limiter uses information of the DG solution essentially only within the troubled cell itself which is identified by a new modified version of the original KXRCF
indicator [42], to build a sequence of hierarchical $L^2$ projection polynomials from zeroth degree to the second or third degree of the DG solution. The second-order and
third-order RKDG methods with the associated multi-resolution WENO limiters are
developed as examples for general high-order RKDG methods, which could maintain
the original order of accuracy in smooth regions and keep essentially non-oscillatory
property near strong discontinuities by gradually degrading from the optimal order
to the first order. The linear weights inside the procedure of the new multi-resolution
WENO limiters can be set as any positive numbers on the condition that they sum
to one. A series of polynomials of different degrees within the troubled cell itself
are applied in a WENO fashion to modify the DG solutions in the troubled cell on
tetrahedral meshes. These new WENO limiters are very simple to construct, and can
be easily implemented to arbitrary high-order accuracy on tetrahedral meshes. Such
spatial reconstruction methodology improves the robustness in the simulation on the
same compact spatial stencil of the original DG methods on tetrahedral meshes. Extensive one-dimensional (run as three-dimensional problems on tetrahedral meshes)
and three-dimensional tests are performed to demonstrate the good performance of
the RKDG methods with new multi-resolution WENO limiters. 相似文献
2.
A Runge Kutta Discontinuous Galerkin Method for Lagrangian Compressible Euler Equations in Two-Dimensions
下载免费PDF全文
![点击此处可从《Communications In Computational Physics》网站下载免费的PDF全文](/ch/ext_images/free.gif)
Zhenzhen Li Xijun Yu Jiang Zhu & Zupeng Jia 《Communications In Computational Physics》2014,15(4):1184-1206
This paper presents a new Lagrangian type scheme for solving the Euler
equations of compressible gas dynamics. In this new scheme the system of equations
is discretized by Runge-Kutta Discontinuous Galerkin (RKDG) method, and the mesh
moves with the fluid flow. The scheme is conservative for the mass, momentum and
total energy and maintains second-order accuracy. The scheme avoids solving the geometrical
part and has free parameters. Results of some numerical tests are presented
to demonstrate the accuracy and the non-oscillatory property of the scheme. 相似文献
3.
Runge-Kutta Discontinuous Galerkin Method Using WENO-Type Limiters: Three-Dimensional Unstructured Meshes
下载免费PDF全文
![点击此处可从《Communications In Computational Physics》网站下载免费的PDF全文](/ch/ext_images/free.gif)
This paper further considers weighted essentially non-oscillatory (WENO) and Hermite weighted essentially non-oscillatory (HWENO) finite volume methods as limiters for Runge-Kutta discontinuous Galerkin (RKDG) methods to solve problems involving nonlinear hyperbolic conservation laws. The application discussed here is the solution of 3-D problems on unstructured meshes. Our numerical tests again demonstrate this is a robust and high order limiting procedure, which simultaneously achieves high order accuracy and sharp non-oscillatory shock transitions. 相似文献
4.
On the Order of Accuracy and Numerical Performance of Two Classes of Finite Volume WENO Schemes
下载免费PDF全文
![点击此处可从《Communications In Computational Physics》网站下载免费的PDF全文](/ch/ext_images/free.gif)
Rui Zhang Mengping Zhang & Chi-Wang Shu 《Communications In Computational Physics》2011,9(3):807-827
In this paper we consider two commonly used classes of finite volume
weighted essentially non-oscillatory (WENO) schemes in two dimensional Cartesian
meshes. We compare them in terms of accuracy, performance for smooth and shocked
solutions, and efficiency in CPU timing. For linear systems both schemes are high
order accurate, however for nonlinear systems, analysis and numerical simulation results
verify that one of them (Class A) is only second order accurate, while the other
(Class B) is high order accurate. The WENO scheme in Class A is easier to implement
and costs less than that in Class B. Numerical experiments indicate that the resolution
for shocked problems is often comparable for schemes in both classes for the same
building blocks and meshes, despite of the difference in their formal order of accuracy.
The results in this paper may give some guidance in the application of high order finite
volume schemes for simulating shocked flows. 相似文献
5.
Adaptive Order WENO Reconstructions for the Semi-Lagrangian Finite Difference Scheme for Advection Problem
下载免费PDF全文
![点击此处可从《Communications In Computational Physics》网站下载免费的PDF全文](/ch/ext_images/free.gif)
Jiajie Chen Xiaofeng Cai Jianxian Qiu & Jing-Mei Qiu 《Communications In Computational Physics》2021,30(1):67-96
We present a new conservative semi-Lagrangian finite difference weighted
essentially non-oscillatory scheme with adaptive order. This is an extension of the
conservative semi-Lagrangian (SL) finite difference WENO scheme in [Qiu and Shu,
JCP, 230 (4) (2011), pp. 863-889], in which linear weights in SL WENO framework
were shown not to exist for variable coefficient problems. Hence, the order of accuracy is not optimal from reconstruction stencils. In this paper, we incorporate a recent
WENO adaptive order (AO) technique [Balsara et al., JCP, 326 (2016), pp. 780-804]
to the SL WENO framework. The new scheme can achieve an optimal high order of
accuracy, while maintaining the properties of mass conservation and non-oscillatory
capture of solutions from the original SL WENO. The positivity-preserving limiter is
further applied to ensure the positivity of solutions. Finally, the scheme is applied to
high dimensional problems by a fourth-order dimensional splitting. We demonstrate
the effectiveness of the new scheme by extensive numerical tests on linear advection
equations, the Vlasov-Poisson system, the guiding center Vlasov model as well as the
incompressible Euler equations. 相似文献
6.
Jun Zhu Xinghui Zhong Chi-Wang Shu & Jianxian Qiu 《Communications In Computational Physics》2016,19(4):944-969
In this paper, we propose a new type of weighted essentially non-oscillatory
(WENO) limiter, which belongs to the class of Hermite WENO (HWENO) limiters, for
the Runge-Kutta discontinuous Galerkin (RKDG) methods solving hyperbolic conservation
laws. This new HWENO limiter is a modification of the simple WENO limiter
proposed recently by Zhong and Shu [29]. Both limiters use information of the DG
solutions only from the target cell and its immediate neighboring cells, thus maintaining
the original compactness of the DG scheme. The goal of both limiters is to obtain
high order accuracy and non-oscillatory properties simultaneously. The main novelty
of the new HWENO limiter in this paper is to reconstruct the polynomial on the target
cell in a least square fashion [8] while the simple WENO limiter [29] is to use the entire
polynomial of the original DG solutions in the neighboring cells with an addition of
a constant for conservation. The modification in this paper improves the robustness
in the computation of problems with strong shocks or contact discontinuities, without
changing the compact stencil of the DG scheme. Numerical results for both one and
two dimensional equations including Euler equations of compressible gas dynamics
are provided to illustrate the viability of this modified limiter. 相似文献
7.
Marie Billaud Friess Jé rô me Breil Pierre-Henri Maire & Mikhail Shashkov 《Communications In Computational Physics》2014,15(2):330-364
In this paper we present recent developments concerning a Cell-Centered
Arbitrary Lagrangian Eulerian (CCALE) strategy using the Moment Of Fluid (MOF)
interface reconstruction for the numerical simulation of multi-material compressible
fluid flows on unstructured grids in cylindrical geometries. Especially, our attention
is focused here on the following points. First, we propose a new formulation of the
scheme used during the Lagrangian phase in the particular case of axisymmetric geometries. Then, the MOF method is considered for multi-interface reconstruction in
cylindrical geometry. Subsequently, a method devoted to the rezoning of polar meshes
is detailed. Finally, a generalization of the hybrid remapping to cylindrical geometries
is presented. These explorations are validated by mean of several test cases using unstructured grid that clearly illustrate the robustness and accuracy of the new method. 相似文献
8.
High-Order Gas-Kinetic Scheme in Curvilinear Coordinates for the Euler and Navier-Stokes Solutions
下载免费PDF全文
![点击此处可从《Communications In Computational Physics》网站下载免费的PDF全文](/ch/ext_images/free.gif)
Liang Pan & Kun Xu 《Communications In Computational Physics》2020,28(4):1321-1351
The high-order gas-kinetic scheme (HGKS) has achieved success in simulating compressible flows with Cartesian meshes. To study the flow problems in general
geometries, such as the flow over a wing-body, the development of HGKS in general
curvilinear coordinates becomes necessary. In this paper, a two-stage fourth-order gas-kinetic scheme is developed for the Euler and Navier-Stokes solutions in the curvilinear coordinates from one-dimensional to three-dimensional computations. Based on
the coordinate transformation, the kinetic equation is transformed first to the computational space, and the flux function in the gas-kinetic scheme is obtained there and
is transformed back to the physical domain for the update of flow variables inside
each control volume. To achieve the expected order of accuracy, the dimension-by-dimension reconstruction based on the WENO scheme is adopted in the computational domain, where the reconstructed variables are the cell averaged Jacobian and the
Jacobian-weighted conservative variables. In the two-stage fourth-order gas-kinetic
scheme, the point values as well as the spatial derivatives of conservative variables at
Gaussian quadrature points have to be used in the evaluation of the time dependent
flux function. The point-wise conservative variables are obtained by ratio of the above
reconstructed data, and the spatial derivatives are reconstructed through orthogonalization in physical space and chain rule. A variety of numerical examples from the
accuracy tests to the solutions with strong discontinuities are presented to validate the
accuracy and robustness of the current scheme for both inviscid and viscous flows.
The precise satisfaction of the geometrical conservation law in non-orthogonal mesh is
also demonstrated through the numerical example. 相似文献
9.
In this paper, a high-order moment-based multi-resolution Hermite
weighted essentially non-oscillatory (HWENO) scheme is designed for hyperbolic conservation laws. The main idea of this scheme is derived from our previous work [J.
Comput. Phys., 446 (2021) 110653], in which the integral averages of the function and
its first order derivative are used to reconstruct both the function and its first order
derivative values at the boundaries. However, in this paper, only the function values at
the Gauss-Lobatto points in the one or two dimensional case need to be reconstructed
by using the information of the zeroth and first order moments. In addition, an extra
modification procedure is used to modify those first order moments in the troubled-cells, which leads to an improvement of stability and an enhancement of resolution
near discontinuities. To obtain the same order of accuracy, the size of the stencil required by this moment-based multi-resolution HWENO scheme is still the same as the
general HWENO scheme and is more compact than the general WENO scheme. Moreover, the linear weights are not unique and are independent of the node position, and
the CFL number can still be 0.6 whether for the one or two dimensional case, which has
to be 0.2 in the two dimensional case for other HWENO schemes. Extensive numerical
examples are given to demonstrate the stability and resolution of such moment-based
multi-resolution HWENO scheme. 相似文献
10.
Arbitrary-Lagrangian-Eulerian One-Step WENO Finite Volume Schemes on Unstructured Triangular Meshes
下载免费PDF全文
![点击此处可从《Communications In Computational Physics》网站下载免费的PDF全文](/ch/ext_images/free.gif)
In this article we present a new class of high order accurate ArbitraryEulerian-Lagrangian (ALE) one-step WENO finite volume schemes for solving nonlinear hyperbolic systems of conservation laws on moving two dimensional unstructured triangular meshes. A WENO reconstruction algorithm is used to achieve high
order accuracy in space and a high order one-step time discretization is achieved by
using the local space-time Galerkin predictor proposed in [25]. For that purpose, a
new element-local weak formulation of the governing PDE is adopted on moving
space-time elements. The space-time basis and test functions are obtained considering Lagrange interpolation polynomials passing through a predefined set of nodes.
Moreover, a polynomial mapping defined by the same local space-time basis functions
as the weak solution of the PDE is used to map the moving physical space-time element onto a space-time reference element. To maintain algorithmic simplicity, the
final ALE one-step finite volume scheme uses moving triangular meshes with straight
edges. This is possible in the ALE framework, which allows a local mesh velocity that
is different from the local fluid velocity. We present numerical convergence rates for
the schemes presented in this paper up to sixth order of accuracy in space and time and
show some classical numerical test problems for the two-dimensional Euler equations
of compressible gas dynamics. 相似文献
11.
A New Approach of High Order Well-Balanced Finite Volume WENO Schemes and Discontinuous Galerkin Methods for a Class of Hyperbolic Systems with Source Terms
下载免费PDF全文
![点击此处可从《Communications In Computational Physics》网站下载免费的PDF全文](/ch/ext_images/free.gif)
Hyperbolic balance laws have steady state solutions in which the flux gradients are
nonzero but are exactly balanced by the source terms. In our earlier work [31–33], we designed
high order well-balanced schemes to a class of hyperbolic systems with separable source terms.
In this paper, we present a different approach to the same purpose: designing high order
well-balanced finite volume weighted essentially non-oscillatory (WENO) schemes and RungeKutta
discontinuous Galerkin (RKDG) finite element methods. We make the observation that
the traditional RKDG methods are capable of maintaining certain steady states exactly, if a
small modification on either the initial condition or the flux is provided. The computational
cost to obtain such a well balanced RKDG method is basically the same as the traditional
RKDG method. The same idea can be applied to the finite volume WENO schemes. We
will first describe the algorithms and prove the well balanced property for the shallow water
equations, and then show that the result can be generalized to a class of other balance laws.
We perform extensive one and two dimensional simulations to verify the properties of these
schemes such as the exact preservation of the balance laws for certain steady state solutions,
the non-oscillatory property for general solutions with discontinuities, and the genuine high
order accuracy in smooth regions. 相似文献
12.
A New Family of High Order Unstructured MOOD and ADER Finite Volume Schemes for Multidimensional Systems of Hyperbolic Conservation Laws
下载免费PDF全文
![点击此处可从《Communications In Computational Physics》网站下载免费的PDF全文](/ch/ext_images/free.gif)
Rapha& euml l Loub& egrave re Michael Dumbser & Steven Diot 《Communications In Computational Physics》2014,16(3):718-763
In this paper, we investigate the coupling of the Multi-dimensional Optimal
Order Detection (MOOD) method and the Arbitrary high order DERivatives (ADER)
approach in order to design a new high order accurate, robust and computationally
efficient Finite Volume (FV) scheme dedicated to solving nonlinear systems of hyperbolic conservation laws on unstructured triangular and tetrahedral meshes in two and
three space dimensions, respectively. The Multi-dimensional Optimal Order Detection
(MOOD) method for 2D and 3D geometries has been introduced in a recent series of
papers for mixed unstructured meshes. It is an arbitrary high-order accurate Finite
Volume scheme in space, using polynomial reconstructions with a posteriori detection
and polynomial degree decrementing processes to deal with shock waves and other
discontinuities. In the following work, the time discretization is performed with an
elegant and efficient one-step ADER procedure. Doing so, we retain the good properties of the MOOD scheme, that is to say, the optimal high-order of accuracy is reached
on smooth solutions, while spurious oscillations near singularities are prevented. The
ADER technique not only reduces the cost of the overall scheme as shown
on a set of numerical tests in 2D and 3D, but also increases the stability of the overall scheme. A systematic comparison between classical unstructured ADER-WENO
schemes and the new ADER-MOOD approach has been carried out for high-order
schemes in space and time in terms of cost, robustness, accuracy and efficiency. The
main finding of this paper is that the combination of ADER with MOOD generally
outperforms the one of ADER and WENO either because at given accuracy MOOD isless expensive (memory and/or CPU time), or because it is more accurate for a given
grid resolution. A large suite of classical numerical test problems has been solved
on unstructured meshes for three challenging multi-dimensional systems of conservation laws: the Euler equations of compressible gas dynamics, the classical equations
of ideal magneto-Hydrodynamics (MHD) and finally the relativistic MHD equations
(RMHD), which constitutes a particularly challenging nonlinear system of hyperbolic
partial differential equation. All tests are run on genuinely unstructured grids composed of simplex elements. 相似文献
13.
H. W. Zheng N. Qin F. C. G. A. Nicolleau & C. Shu 《Communications In Computational Physics》2011,9(1):68-88
An anisotropic solution adaptive method based on unstructured quadrilateral
meshes for inviscid compressible flows is proposed. The data structure, the directional
refinement and coarsening, including the method for initializing the refined
new cells, for the anisotropic adaptive method are described. It provides efficient high
resolution of flow features, which are aligned with the original quadrilateral mesh
structures. Five different cases are provided to show that it could be used to resolve
the anisotropic flow features and be applied to model the complex geometry as well as
to keep a relative high order of accuracy on an efficient anisotropic mesh. 相似文献
14.
Angelo L. Scandaliato & Meng-Sing Liou 《Communications In Computational Physics》2012,12(4):1096-1120
In this paper we demonstrate the accuracy and robustness of combining the
advection upwind splitting method (AUSM), specifically AUSM+-UP [9], with high-order upwind-biased interpolation procedures, the weighted essentially non-oscillatory
(WENO-JS) scheme [8] and its variations [2, 7], and the monotonicity preserving (MP)
scheme [16], for solving the Euler equations. MP is found to be more effective than the
three WENO variations studied. AUSM+-UP is also shown to be free of the so-called "carbuncle" phenomenon with the high-order interpolation. The characteristic variables are preferred for interpolation after comparing the results using primitive and
conservative variables, even though they require additional matrix-vector operations.
Results using the Roe flux with an entropy fix and the Lax-Friedrichs approximate Riemann solvers are also included for comparison. In addition, four reflective boundary
condition implementations are compared for their effects on residual convergence and
solution accuracy. Finally, a measure for quantifying the efficiency of obtaining high
order solutions is proposed; the measure reveals that a maximum return is reached
after which no improvement in accuracy is possible for a given grid size. 相似文献
15.
A High-Order Central ENO Finite-Volume Scheme for Three-Dimensional Low-Speed Viscous Flows on Unstructured Mesh
下载免费PDF全文
![点击此处可从《Communications In Computational Physics》网站下载免费的PDF全文](/ch/ext_images/free.gif)
Marc R. J. Charest Clinton P. T. Groth & Pierre Q. Gauthier 《Communications In Computational Physics》2015,17(3):615-656
High-order discretization techniques offer the potential to significantly reduce
the computational costs necessary to obtain accurate predictions when compared
to lower-order methods. However, efficient and universally-applicable high-order
discretizations remain somewhat illusive, especially for more arbitrary unstructured
meshes and for incompressible/low-speed flows. A novel, high-order, central essentially
non-oscillatory (CENO), cell-centered, finite-volume scheme is proposed for
the solution of the conservation equations of viscous, incompressible flows on three-dimensional
unstructured meshes. Similar to finite element methods, coordinate transformations
are used to maintain the scheme's order of accuracy even when dealing
with arbitrarily-shaped cells having non-planar faces. The proposed scheme is applied
to the pseudo-compressibility formulation of the steady and unsteady Navier-Stokes
equations and the resulting discretized equations are solved with a parallel implicit
Newton-Krylov algorithm. For unsteady flows, a dual-time stepping approach
is adopted and the resulting temporal derivatives are discretized using the family of
high-order backward difference formulas (BDF). The proposed finite-volume scheme
for fully unstructured mesh is demonstrated to provide both fast and accurate solutions
for steady and unsteady viscous flows. 相似文献
16.
High Order Finite Difference WENO Methods with Unequal-Sized Sub-Stencils for the Degasperis-Procesi Type Equations
下载免费PDF全文
![点击此处可从《Communications In Computational Physics》网站下载免费的PDF全文](/ch/ext_images/free.gif)
Jianfang Lin Yan Xu Huiwen Xue & Xinghui Zhong 《Communications In Computational Physics》2022,31(3):913-946
In this paper, we develop two finite difference weighted essentially
non-oscillatory (WENO) schemes with unequal-sized sub-stencils for solving the
Degasperis-Procesi (DP) and $\mu$-Degasperis-Procesi ($\mu$DP) equations, which contain
nonlinear high order derivatives, and possibly peakon solutions or shock waves. By
introducing auxiliary variable(s), we rewrite the DP equation as a hyperbolic-elliptic
system, and the $\mu$DP equation as a first order system. Then we choose a linear finite
difference scheme with suitable order of accuracy for the auxiliary variable(s), and
two finite difference WENO schemes with unequal-sized sub-stencils for the primal
variable. One WENO scheme uses one large stencil and several smaller stencils, and
the other WENO scheme is based on the multi-resolution framework which uses a series of unequal-sized hierarchical central stencils. Comparing with the classical WENO
scheme which uses several small stencils of the same size to make up a big stencil, both
WENO schemes with unequal-sized sub-stencils are simple in the choice of the stencil
and enjoy the freedom of arbitrary positive linear weights. Another advantage is that
the final reconstructed polynomial on the target cell is a polynomial of the same degree as the polynomial over the big stencil, while the classical finite difference WENO
reconstruction can only be obtained for specific points inside the target interval. Numerical tests are provided to demonstrate the high order accuracy and non-oscillatory
properties of the proposed schemes. 相似文献
17.
Improvement on Spherical Symmetry in Two-Dimensional Cylindrical Coordinates for a Class of Control Volume Lagrangian Schemes
下载免费PDF全文
![点击此处可从《Communications In Computational Physics》网站下载免费的PDF全文](/ch/ext_images/free.gif)
In [14], Maire developed a class of cell-centered Lagrangian schemes for
solving Euler equations of compressible gas dynamics in cylindrical coordinates. These
schemes use a node-based discretization of the numerical fluxes. The control volume version has several distinguished properties, including the conservation of mass,
momentum and total energy and compatibility with the geometric conservation law
(GCL). However, it also has a limitation in that it cannot preserve spherical symmetry
for one-dimensional spherical flow. An alternative is also given to use the first order
area-weighted approach which can ensure spherical symmetry, at the price of sacrificing conservation of momentum. In this paper, we apply the methodology proposed in
our recent work [8] to the first order control volume scheme of Maire in [14] to obtain
the spherical symmetry property. The modified scheme can preserve one-dimensional
spherical symmetry in a two-dimensional cylindrical geometry when computed on an
equal-angle-zoned initial grid, and meanwhile it maintains its original good properties
such as conservation and GCL. Several two-dimensional numerical examples in cylindrical coordinates are presented to demonstrate the good performance of the scheme
in terms of symmetry, non-oscillation and robustness properties. 相似文献
18.
A Hermite WENO Method with Modified Ghost Fluid Method for Compressible Two-Medium Flow Problems
下载免费PDF全文
![点击此处可从《Communications In Computational Physics》网站下载免费的PDF全文](/ch/ext_images/free.gif)
Zhuang Zhao Yong-Tao Zhang Yibing Chen & Jianxian Qiu 《Communications In Computational Physics》2021,30(3):851-873
In this paper, we develop a novel approach by combining a new robust finite difference Hermite weighted essentially non-oscillatory (HWENO) method [51]
with the modified ghost fluid method (MGFM) [25] to simulate the compressible two-medium flow problems. The main idea is that we first use the technique of the MGFM
to transform a two-medium flow problem to two single-medium cases by defining the
ghost fluids status based on the predicted interface status. Then the efficient and robust
HWENO finite difference method is applied for solving the single-medium flow cases.
By using immediate neighbor information to deal with both the solution and its derivatives, the fifth order finite difference HWENO scheme adopted in this paper is more
compact and has higher resolution than the classical fifth order finite difference WENO
scheme of Jiang and Shu [14]. Furthermore, by combining the HWENO scheme with
the MGFM to simulate the two-medium flow problems, less ghost point information
is needed than that in using the classical WENO scheme in order to obtain the same
numerical accuracy. Various one-dimensional and two-dimensional two-medium flow
problems are solved to illustrate the good performances of the proposed method. 相似文献
19.
High-Order Positivity-Preserving Well-Balanced Discontinuous Galerkin Methods for Euler Equations with Gravitation on Unstructured Meshes
下载免费PDF全文
![点击此处可从《Communications In Computational Physics》网站下载免费的PDF全文](/ch/ext_images/free.gif)
Weijie Zhang Yulong Xing Yinhua Xia & Yan Xu 《Communications In Computational Physics》2022,31(3):771-815
In this paper, we propose a high-order accurate discontinuous Galerkin
(DG) method for the compressible Euler equations under gravitational fields on unstructured meshes. The scheme preserves a general hydrostatic equilibrium state and
provably guarantees the positivity of density and pressure at the same time. Comparing with the work on the well-balanced scheme for Euler equations with gravitation
on rectangular meshes, the extension to triangular meshes is conceptually plausible
but highly nontrivial. We first introduce a special way to recover the equilibrium state
and then design a group of novel variables at the interface of two adjacent cells, which
plays an important role in the well-balanced and positivity-preserving properties. One
main challenge is that the well-balanced schemes may not have the weak positivity
property. In order to achieve the well-balanced and positivity-preserving properties
simultaneously while maintaining high-order accuracy, we carefully design DG spatial discretization with well-balanced numerical fluxes and suitable source term approximation. For the ideal gas, we prove that the resulting well-balanced scheme, coupled with strong stability preserving time discretizations, satisfies a weak positivity
property. A simple existing limiter can be applied to enforce the positivity-preserving
property, without losing high-order accuracy and conservation. Extensive one- and
two-dimensional numerical examples demonstrate the desired properties of the proposed scheme, as well as its high resolution and robustness. 相似文献
20.
A. Ferrari M. Dumbser E. F. Toro & A. Armanini 《Communications In Computational Physics》2008,4(2):378-404
The purpose of this paper is to solve some of the trouble spots of the classical SPH method by proposing an alternative approach. First, we focus on the problem of the stability for two different SPH schemes, one is based on the approach of
Vila [25] and the other is proposed in this article which mimics the classical 1D Lax Wendroff scheme. In both approaches the classical SPH artificial viscosity term is removed preserving nevertheless the linear stability of the methods, demonstrated via
the von Neumann stability analysis. Moreover, the issue of the consistency for the
equations of gas dynamics is analyzed. An alternative approach is proposed that consists of using Godunov-type SPH schemes in Lagrangian coordinates. This not only
provides an improvement in accuracy of the numerical solutions, but also assures that
the consistency condition on the gradient of the kernel function is satisfied using an
equidistant distribution of particles in Lagrangian mass coordinates. Three different
Riemann solvers are implemented for the first-order Godunov type SPH schemes in
Lagrangian coordinates, namely the Godunov flux based on the exact Riemann solver,
the Rusanov flux and a new modified Roe flux, following the work of Munz [17]. Some
well-known numerical 1D shock tube test cases [22] are solved, comparing the numerical solutions of the Godunov-type SPH schemes in Lagrangian coordinates with the
first-order Godunov finite volume method in Eulerian coordinates and the standard
SPH scheme with Monaghan's viscosity term. 相似文献