共查询到20条相似文献,搜索用时 31 毫秒
1.
High Order Finite Difference Hermite WENO Fixed-Point Fast Sweeping Method for Static Hamilton-Jacobi Equations
下载免费PDF全文
![点击此处可从《Communications In Computational Physics》网站下载免费的PDF全文](/ch/ext_images/free.gif)
In this paper, we combine the nonlinear HWENO reconstruction in [43] and
the fixed-point iteration with Gauss-Seidel fast sweeping strategy, to solve the static
Hamilton-Jacobi equations in a novel HWENO framework recently developed in [22].
The proposed HWENO frameworks enjoys several advantages. First, compared with
the traditional HWENO framework, the proposed methods do not need to introduce
additional auxiliary equations to update the derivatives of the unknown function $\phi$.
They are now computed from the current value of $\phi$ and the previous spatial derivatives of $\phi$. This approach saves the computational storage and CPU time, which greatly
improves the computational efficiency of the traditional HWENO scheme. In addition,
compared with the traditional WENO method, reconstruction stencil of the HWENO
methods becomes more compact, their boundary treatment is simpler, and the numerical errors are smaller on the same mesh. Second, the fixed-point fast sweeping method
is used to update the numerical approximation. It is an explicit method and does
not involve the inverse operation of nonlinear Hamiltonian, therefore any Hamilton-Jacobi equations with complex Hamiltonian can be solved easily. It also resolves some
known issues, including that the iterative number is very sensitive to the parameter $ε$ used in the nonlinear weights, as observed in previous studies. Finally, to further
reduce the computational cost, a hybrid strategy is also presented. Extensive numerical experiments are performed on two-dimensional problems, which demonstrate the
good performance of the proposed fixed-point fast sweeping HWENO methods. 相似文献
2.
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. 相似文献
3.
Yifei Wan & Yinhua Xia 《Communications In Computational Physics》2023,33(5):1270-1331
For steady Euler equations in complex boundary domains, high-order shockcapturing schemes usually suffer not only from the difficulty of steady-state convergence but also from the problem of dealing with physical boundaries on Cartesian
grids to achieve uniform high-order accuracy. In this paper, we utilize a fifth-order
finite difference hybrid WENO scheme to simulate steady Euler equations, and the
same fifth-order WENO extrapolation methods are developed to handle the curved
boundary. The values of the ghost points outside the physical boundary can be obtained by applying WENO extrapolation near the boundary, involving normal derivatives acquired by the simplified inverse Lax-Wendroff procedure. Both equivalent expressions involving curvature and numerical differentiation are utilized to transform
the tangential derivatives along the curved solid wall boundary. This hybrid WENO
scheme is robust for steady-state convergence and maintains high-order accuracy in
the smooth region even with the solid wall boundary condition. Besides, the essentially non-oscillation property is achieved. The numerical spectral analysis also shows
that this hybrid WENO scheme has low dispersion and dissipation errors. Numerical
examples are presented to validate the high-order accuracy and robust performance of
the hybrid scheme for steady Euler equations in curved domains with Cartesian grids. 相似文献
4.
In this paper, we introduce a new type of troubled-cell indicator to improve
hybrid weighted essentially non-oscillatory (WENO) schemes for solving the hyperbolic conservation laws. The hybrid WENO schemes selectively adopt the high-order
linear upwind scheme or the WENO scheme to avoid the local characteristic decompositions and calculations of the nonlinear weights in smooth regions. Therefore,
they can reduce computational cost while maintaining non-oscillatory properties in
non-smooth regions. Reliable troubled-cell indicators are essential for efficient hybrid
WENO methods. Most of troubled-cell indicators require proper parameters to detect
discontinuities precisely, but it is very difficult to determine the parameters automatically. We develop a new troubled-cell indicator derived from the mean value theorem
that does not require any variable parameters. Additionally, we investigate the characteristics of indicator variable; one of the conserved properties or the entropy is considered as indicator variable. Detailed numerical tests for 1D and 2D Euler equations are
conducted to demonstrate the performance of the proposed indicator. The results with
the proposed troubled-cell indicator are in good agreement with pure WENO schemes.
Also the new indicator has advantages in the computational cost compared with the
other indicators. 相似文献
5.
Michael Dumbser Ariunaa Uuriintsetseg & Olindo Zanotti 《Communications In Computational Physics》2013,14(2):301-327
In this article we present a new family of high order accurate Arbitrary
Lagrangian-Eulerian one-step WENO finite volume schemes for the solution of stiff
hyperbolic balance laws. High order accuracy in space is obtained with a standard
WENO reconstruction algorithm and high order in time is obtained using the local
space-time discontinuous Galerkin method recently proposed in [20]. In the Lagrangian
framework considered here, the local space-time DG predictor is based on a weak
formulation of the governing PDE on a moving space-time element. For the space-time basis and test functions we use Lagrange interpolation polynomials defined by
tensor-product Gauss-Legendre quadrature points. The moving space-time elements
are mapped to a reference element using an isoparametric approach, i.e. the space-time mapping is defined by the same basis functions as the weak solution of the PDE.
We show some computational examples in one space-dimension for non-stiff and for
stiff balance laws, in particular for the Euler equations of compressible gas dynamics,
for the resistive relativistic MHD equations, and for the relativistic radiation hydrodynamics equations. Numerical convergence results are presented for the stiff case up to
sixth order of accuracy in space and time and for the non-stiff case up to eighth order
of accuracy in space and time. 相似文献
6.
Jun Zhu & Jianxian Qiu 《Communications In Computational Physics》2020,27(3):897-920
In this paper, a new type of third-order and fourth-order weighted essentially non-oscillatory (WENO) schemes is designed for simulating the Hamilton-Jacobi equations on triangular meshes. We design such schemes with the use of the nodal information defined on five unequal-sized spatial stencils, the application of monotone Hamiltonians as a building block, the artificial set of positive linear weights to make up high-order approximations in smooth regions simultaneously avoiding spurious oscillations nearby discontinuities of the derivatives of the solutions. The spatial reconstructions are convex combinations of the derivatives of a modified cubic/quartic polynomial defined on a big spatial stencil and four quadratic polynomials defined on small spatial stencils, and a third-order TVD Runge-Kutta method is used for the time discretization. The main advantages of these WENO schemes are their efficiency, simplicity, and can be easily implemented to higher dimensional unstructured meshes. Extensive numerical tests are performed to illustrate the good performance of such new WENO schemes. 相似文献
7.
Adaptive Bayesian Inference for Discontinuous Inverse Problems,Application to Hyperbolic Conservation Laws
下载免费PDF全文
![点击此处可从《Communications In Computational Physics》网站下载免费的PDF全文](/ch/ext_images/free.gif)
Alexandre Birolleau Gaë l Poë tte & Didier Lucor 《Communications In Computational Physics》2014,16(1):1-34
Various works from the literature aimed at accelerating Bayesian inference
in inverse problems. Stochastic spectral methods have been recently proposed as surrogate approximations of the forward uncertainty propagation model over the support
of the prior distribution. These representations are efficient because they allow affordable simulation of a large number of samples from the posterior distribution. Unfortunately, they do not perform well when the forward model exhibits strong nonlinear
behavior with respect to its input.In this work, we first relate the fast (exponential) $L^2$-convergence of the forward
approximation to the fast (exponential) convergence (in terms of Kullback-Leibler divergence) of the approximate posterior. In particular, we prove that in case the prior
distribution is uniform, the posterior is at least twice as fast as the convergence rate of
the forward model in those norms. The Bayesian inference strategy is developed in the
framework of a stochastic spectral projection method. The predicted convergence rates
are then demonstrated for simple nonlinear inverse problems of varying smoothness.We then propose an efficient numerical approach for the Bayesian solution of inverse problems presenting strongly nonlinear or discontinuous system responses. This
comes with the improvement of the forward model that is adaptively approximated by
an iterative generalized Polynomial Chaos-based representation. The numerical approximations and predicted convergence rates of the former approach are compared
to the new iterative numerical method for nonlinear time-dependent test cases of varying dimension and complexity, which are relevant regarding our hydrodynamics motivations and therefore regarding hyperbolic conservation laws and the apparition of
discontinuities in finite time. 相似文献
8.
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. 相似文献
9.
Piecewise Polynomial Mapping Method and Corresponding WENO Scheme with Improved Resolution
下载免费PDF全文
![点击此处可从《Communications In Computational Physics》网站下载免费的PDF全文](/ch/ext_images/free.gif)
The method of mapping function was first proposed by Henrick et al. [J.
Comput. Phys. 207:542-547 (2005)] to adjust nonlinear weights in [0,1] for the fifth-order
WENO scheme, and through which the requirement of convergence order is
satisfied and the performance of the scheme is improved. Different from Henrick's
method, a concept of piecewise polynomial function is proposed in this study and
corresponding WENO schemes are obtained. The advantage of the new method is
that the function can have a gentle profile at the location of the linear weight (or the
mapped nonlinear weight can be close to its linear counterpart), and therefore is favorable
for the resolution enhancement. Besides, the function also has the flexibility
of quick convergence to identity mapping near two endpoints of [0,1], which is favorable
for improved numerical stability. The fourth-, fifth- and sixth-order polynomial
functions are constructed correspondingly with different emphasis on aforementioned
flatness and convergence. Among them, the fifth-order version has the flattest profile.
To check the performance of the methods, the 1-D Shu-Osher problem, the 2-D
Riemann problem and the double Mach reflection are tested with the comparison of
WENO-M, WENO-Z and WENO-NS. The proposed new methods show the best resolution
for describing shear-layer instability of the Riemann problem, and they also
indicate high resolution in computations of double Mach reflection, where only these
proposed schemes successfully resolved the vortex-pairing phenomenon. Other investigations
have shown that the single polynomial mapping function has no advantage
over the proposed piecewise one, and it is of no evident benefit to use the proposed
method for the symmetric fifth-order WENO. Overall, the fifth-order piecewise polynomial
and corresponding WENO scheme are suggested for resolution improvement. 相似文献
10.
A Comparative Study of Rosenbrock-Type and Implicit Runge-Kutta Time Integration for Discontinuous Galerkin Method for Unsteady 3D Compressible Navier-Stokes equations
下载免费PDF全文
![点击此处可从《Communications In Computational Physics》网站下载免费的PDF全文](/ch/ext_images/free.gif)
Xiaodong Liu Yidong Xia Hong Luo & Lijun Xuan 《Communications In Computational Physics》2016,20(4):1016-1044
A comparative study of two classes of third-order implicit time integration
schemes is presented for a third-order hierarchical WENO reconstructed discontinuous
Galerkin (rDG) method to solve the 3D unsteady compressible Navier-Stokes
equations: — 1) the explicit first stage, single diagonally implicit Runge-Kutta (ESDIRK3)
scheme, and 2) the Rosenbrock-Wanner (ROW) schemes based on the differential
algebraic equations (DAEs) of Index-2. Compared with the ESDIRK3 scheme,
a remarkable feature of the ROW schemes is that, they only require one approximate
Jacobian matrix calculation every time step, thus considerably reducing the overall
computational cost. A variety of test cases, ranging from inviscid flows to DNS of
turbulent flows, are presented to assess the performance of these schemes. Numerical
experiments demonstrate that the third-order ROW scheme for the DAEs of index-2
can not only achieve the designed formal order of temporal convergence accuracy in
a benchmark test, but also require significantly less computing time than its ESDIRK3
counterpart to converge to the same level of discretization errors in all of the flow
simulations in this study, indicating that the ROW methods provide an attractive alternative
for the higher-order time-accurate integration of the unsteady compressible
Navier-Stokes equations. 相似文献
11.
High Order Schemes on Three-Dimensional General Polyhedral Meshes — Application to the Level Set Method
下载免费PDF全文
![点击此处可从《Communications In Computational Physics》网站下载免费的PDF全文](/ch/ext_images/free.gif)
In this article, we detail the methodology developed to construct arbitrarily high order schemes — linear and WENO — on 3D mixed-element unstructured
meshes made up of general convex polyhedral elements. The approach is tailored
specifically for the solution of scalar level set equations for application to incompressible two-phase flow problems. The construction of WENO schemes on 3D unstructured meshes is notoriously difficult, as it involves a much higher level of complexity
than 2D approaches. This due to the multiplicity of geometrical considerations introduced by the extra dimension, especially on mixed-element meshes. Therefore, we
have specifically developed a number of algorithms to handle mixed-element meshes
composed of convex polyhedra with convex polygonal faces. The contribution of this
work concerns several areas of interest: the formulation of an improved methodology
in 3D, the minimisation of computational runtime in the implementation through the
maximum use of pre-processing operations, the generation of novel methods to handle complex 3D mixed-element meshes and finally the application of the method to
the transport of a scalar level set. 相似文献
12.
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. 相似文献
13.
Multi-Valued Solution and Level Set Methods in Computational High Frequency Wave Propagation
下载免费PDF全文
![点击此处可从《Communications In Computational Physics》网站下载免费的PDF全文](/ch/ext_images/free.gif)
We review the level set methods for computing multi-valued solutions to a class of nonlinear first order partial differential equations, including Hamilton-Jacobi equations, quasi-linear hyperbolic equations, and conservative transport equations with multi-valued transport speeds. The multivalued solutions are embedded as the zeros of a set of scalar functions that solve the initial value problems of a time dependent partial differential equation in an augmented space. We discuss the essential ideas behind the techniques, the coupling of these techniques to the projection of the interaction of zero level sets and a collection of applications including the computation of the semiclassical limit for Schrödinger equations and the high frequency geometrical optics limits of linear wave equations. 相似文献
14.
Ruo Li & Wei Zhong 《Communications In Computational Physics》2021,30(5):1545-1588
We propose a new family of mapped WENO schemes by using several adaptive control functions and a smoothing approximation of the signum function. The
proposed schemes admit an extensive permitted range of the parameters in the mapping functions. Consequently, they have the capacity to achieve optimal convergence
rates, even near critical points. Particularly, the new schemes with fine-tuned parameters illustrates a significant advantage when solving problems with discontinuities. It
produces numerical solutions with high resolution without generating spurious oscillations, especially for long output times. 相似文献
15.
This paper presents a new and better suited formulation to implement the
limiting projection to high-order schemes that make use of high-order local reconstructions
for hyperbolic conservation laws. The scheme, so-called MCV-WENO4 (multi-moment
Constrained finite Volume with WENO limiter of 4th order) method, is an
extension of the MCV method of Ii & Xiao (2009) by adding the 1st order derivative
(gradient or slope) at the cell center as an additional constraint for the cell-wise local
reconstruction. The gradient is computed from a limiting projection using the WENO
(weighted essentially non-oscillatory) reconstruction that is built from the nodal values
at 5 solution points within 3 neighboring cells. Different from other existing methods
where only the cell-average value is used in the WENO reconstruction, the present
method takes account of the solution structure within each mesh cell, and thus minimizes
the stencil for reconstruction. The resulting scheme has 4th-order accuracy and
is of significant advantage in algorithmic simplicity and computational efficiency. Numerical
results of one and two dimensional benchmark tests for scalar and Euler conservation
laws are shown to verify the accuracy and oscillation-less property of the
scheme. 相似文献
16.
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. 相似文献
17.
Convergence Detection in Direct Simulation Monte Carlo Calculations for Steady State Flows
下载免费PDF全文
![点击此处可从《Communications In Computational Physics》网站下载免费的PDF全文](/ch/ext_images/free.gif)
A new criterion is presented to detect global convergence to steady state,
and to identify local transient characteristics, during rarefied gas flow simulations performed
using the direct simulation Monte Carlo (DSMC) method. Unlike deterministic
computational fluid dynamics (CFD) schemes, DSMC is generally subject to large statistical
scatter in instantaneous flow property evaluations, which prevents the use of
residual tracking procedures as are often employed in CFD simulations. However,
reliable prediction of the time to reach steady state is necessary for initialization of
DSMC sampling operations. Techniques currently used in DSMC to identify steady
state convergence are usually insensitive to weak transient behavior in small regions
of relatively low density or recirculating flow. The proposed convergence criterion is
developed with the goal of properly identifying such weak transient behavior, while
adding negligible computational expense and allowing simple implementation in any
existing DSMC code. Benefits of the proposed technique over existing convergence
detection methods are demonstrated for representative nozzle/plume expansion flow,
hypersonic blunt body flow and driven cavity flow problems. 相似文献
18.
Conservative Semi-Lagrangian Finite Difference WENO Formulations with Applications to the Vlasov Equation
下载免费PDF全文
![点击此处可从《Communications In Computational Physics》网站下载免费的PDF全文](/ch/ext_images/free.gif)
In this paper, we propose a new conservative semi-Lagrangian (SL) finite
difference (FD) WENO scheme for linear advection equations, which can serve as a
base scheme for the Vlasov equation by Strang splitting [4]. The reconstruction procedure
in the proposed SL FD scheme is the same as the one used in the SL finite volume
(FV) WENO scheme [3]. However, instead of inputting cell averages and approximate
the integral form of the equation in a FV scheme, we input point values and approximate
the differential form of equation in a FD spirit, yet retaining very high order
(fifth order in our experiment) spatial accuracy. The advantage of using point values,
rather than cell averages, is to avoid the second order spatial error, due to the shearing
in velocity (v) and electrical field (E) over a cell when performing the Strang splitting
to the Vlasov equation. As a result, the proposed scheme has very high spatial accuracy,
compared with second order spatial accuracy for Strang split SL FV scheme for
solving the Vlasov-Poisson (VP) system. We perform numerical experiments on linear
advection, rigid body rotation problem; and on the Landau damping and two-stream
instabilities by solving the VP system. For comparison, we also apply (1) the conservative
SL FD WENO scheme, proposed in [22] for incompressible advection problem, (2)
the conservative SL FD WENO scheme proposed in [21] and (3) the non-conservative
version of the SL FD WENO scheme in [3] to the same test problems. The performances
of different schemes are compared by the error table, solution resolution of sharp interface,
and by tracking the conservation of physical norms, energies and entropies,
which should be physically preserved. 相似文献
19.
Ruo Li & Wei Zhong 《Communications In Computational Physics》2022,31(5):1362-1401
A new type of finite volume WENO schemes for hyperbolic problems was
devised in [33] by introducing the order-preserving (OP) criterion. In this continuing
work, we extend the OP criterion to the WENO-Z-type schemes. We firstly rewrite the
formulas of the Z-type weights in a uniform form from a mapping perspective inspired
by extensive numerical observations. Accordingly, we build the concept of the locally
order-preserving (LOP) mapping which is an extension of the order-preserving (OP)
mapping and the resultant improved WENO-Z-type schemes are denoted as LOP-GMWENO-X. There are four major advantages of the LOP-GMWENO-X schemes superior to the existing WENO-Z-type schemes. Firstly, the new schemes can amend
the serious drawback of the existing WENO-Z-type schemes that most of them suffer
from either producing severe spurious oscillations or failing to obtain high resolutions
in long calculations of hyperbolic problems with discontinuities. Secondly, they can
maintain considerably high resolutions on solving problems with high-order critical
points at long output times. Thirdly, they can obtain evidently higher resolution in the
region with high-frequency but smooth waves. Finally, they can significantly decrease
the post-shock oscillations for simulations of some 2D problems with strong shock
waves. Extensive benchmark examples are conducted to illustrate these advantages. 相似文献
20.
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. 相似文献