共查询到20条相似文献,搜索用时 280 毫秒
1.
Zhijun Shen Wei Yan & Guangwei Yuan 《Communications In Computational Physics》2014,15(5):1320-1342
The carbuncle phenomenon has been regarded as a spurious solution produced by most of contact-preserving methods. The hybrid method of combining high
resolution flux with more dissipative solver is an attractive attempt to cure this kind
of non-physical phenomenon. In this paper, a matrix-based stability analysis for 2-D
Euler equations is performed to explore the cause of instability of numerical schemes.
By combining the Roe with HLL flux in different directions and different flux components, we give an interesting explanation to the linear numerical instability. Based on
such analysis, some hybrid schemes are compared to illustrate different mechanisms in
controlling shock instability. Numerical experiments are presented to verify our analysis results. The conclusion is that the scheme of restricting directly instability source
is more stable than other hybrid schemes. 相似文献
2.
A Class of Hybrid DG/FV Methods for Conservation Laws III: Two-Dimensional Euler Equations
下载免费PDF全文
![点击此处可从《Communications In Computational Physics》网站下载免费的PDF全文](/ch/ext_images/free.gif)
Laiping Zhang Wei Liu Lixin He & Xiaogang Deng 《Communications In Computational Physics》2012,12(1):284-314
A concept of "static reconstruction" and "dynamic reconstruction" was introduced for higher-order (third-order or more) numerical methods in our previous
work. Based on this concept, a class of hybrid DG/FV methods had been developed
for one-dimensional conservation law using a "hybrid reconstruction" approach, and
extended to two-dimensional scalar equations on triangular and Cartesian/triangular
hybrid grids. In the hybrid DG/FV schemes, the lower-order derivatives of the piecewise polynomial are computed locally in a cell by the traditional DG method (called
as "dynamic reconstruction"), while the higher-order derivatives are reconstructed by
the "static reconstruction" of the FV method, using the known lower-order derivatives
in the cell itself and in its adjacent neighboring cells. In this paper, the hybrid DG/FV
schemes are extended to two-dimensional Euler equations on triangular and Cartesian/triangular hybrid grids. Some typical test cases are presented to demonstrate
the performance of the hybrid DG/FV methods, including the standard vortex evolution problem with exact solution, isentropic vortex/weak shock wave interaction,
subsonic flows past a circular cylinder and a three-element airfoil (30P30N), transonic
flow past a NACA0012 airfoil. The accuracy study shows that the hybrid DG/FV
method achieves the desired third-order accuracy, and the applications demonstrate
that they can capture the flow structure accurately, and can reduce the CPU time and
memory requirement greatly than the traditional DG method with the same order of
accuracy. 相似文献
3.
A conservative modification to the ghost fluid method (GFM) is developed
for compressible multiphase flows. The motivation is to eliminate or reduce the conservation
error of the GFM without affecting its performance. We track the conservative
variables near the material interface and use this information to modify the numerical
solution for an interfacing cell when the interface has passed the cell. The modification
procedure can be used on the GFM with any base schemes. In this paper we use the
fifth order finite difference WENO scheme for the spatial discretization and the third
order TVD Runge-Kutta method for the time discretization. The level set method is
used to capture the interface. Numerical experiments show that the method is at least
mass and momentum conservative and is in general comparable in numerical resolution
with the original GFM. 相似文献
4.
A High Order Sharp-Interface Method with Local Time Stepping for Compressible Multiphase Flows
下载免费PDF全文
![点击此处可从《Communications In Computational Physics》网站下载免费的PDF全文](/ch/ext_images/free.gif)
Angela Ferrari Claus-Dieter Munz & Bernhard Weigand 《Communications In Computational Physics》2011,9(1):205-230
In this paper, a new sharp-interface approach to simulate compressible
multiphase flows is proposed. The new scheme consists of a high order WENO finite volume scheme for solving the Euler equations coupled with a high order path-conservative
discontinuous Galerkin finite element scheme to evolve an indicator function
that tracks the material interface. At the interface our method applies ghost cells
to compute the numerical flux, as the ghost fluid method. However, unlike the original
ghost fluid scheme of Fedkiw et al. [15], the state of the ghost fluid is derived
from an approximate-state Riemann solver, similar to the approach proposed in [25],
but based on a much simpler formulation. Our formulation leads only to one single
scalar nonlinear algebraic equation that has to be solved at the interface, instead of
the system used in [25]. Away from the interface, we use the new general Osher-type
flux recently proposed by Dumbser and Toro [13], which is a simple but complete Riemann
solver, applicable to general hyperbolic conservation laws. The time integration
is performed using a fully-discrete one-step scheme, based on the approaches recently
proposed in [5, 7]. This allows us to evolve the system also with time-accurate local
time stepping. Due to the sub-cell resolution and the subsequent more restrictive
time-step constraint of the DG scheme, a local evolution for the indicator function is
applied, which is matched with the finite volume scheme for the solution of the Euler
equations that runs with a larger time step. The use of a locally optimal time step
avoids the introduction of excessive numerical diffusion in the finite volume scheme.
Two different fluids have been used, namely an ideal gas and a weakly compressible
fluid modeled by the Tait equation. Several tests have been computed to assess the
accuracy and the performance of the new high order scheme. A verification of our
algorithm has been carefully carried out using exact solutions as well as a comparison
with other numerical reference solutions. The material interface is resolved sharply
and accurately without spurious oscillations in the pressure field. 相似文献
5.
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. 相似文献
6.
B. Costa W. S. Don D. Gottlieb & R. Sendersky 《Communications In Computational Physics》2006,1(3):548-574
The multi-domain hybrid Spectral-WENO method (Hybrid) is introduced for the numerical solution of two-dimensional nonlinear hyperbolic systems in a Cartesian physical domain which is partitioned into a grid of rectangular subdomains. The main idea of the Hybrid scheme is to conjugate the spectral and WENO methods for solving problems with shock or high gradients such that the scheme adapts its solver spatially and temporally depending on the smoothness of the solution in a given subdomain. Built as a multi-domain method, an adaptive algorithm is used to keep the solutions parts exhibiting high gradients and discontinuities always inside WENO subdomains while the smooth parts of the solution are kept inside spectral ones, avoiding oscillations related to the well-known Gibbs phenomenon and increasing the numerical efficiency of the overall scheme. A higher order version of the multi-resolution analysis proposed by Harten is used to determine the smoothness of the solution in each subdomain. We also discuss interface conditions for the two-dimensional problem and the switching procedure between WENO and spectral subdomains. The Hybrid method is applied to the two-dimensional Shock-Vortex Interaction and the Richtmyer-Meshkov Instability (RMI) problems. 相似文献
7.
Comparison of Fifth-Order WENO Scheme and Finite Volume WENO-Gas-Kinetic Scheme for Inviscid and Viscous Flow Simulation
下载免费PDF全文
![点击此处可从《Communications In Computational Physics》网站下载免费的PDF全文](/ch/ext_images/free.gif)
The development of high-order schemes has been mostly concentrated on
the limiters and high-order reconstruction techniques. In this paper, the effect of the
flux functions on the performance of high-order schemes will be studied. Based on the
same WENO reconstruction, two schemes with different flux functions, i.e., the fifth-order WENO method and the WENO-Gas-Kinetic scheme (WENO-GKS), will be compared. The fifth-order finite difference WENO-SW scheme is a characteristic variable
reconstruction based method which uses the Steger-Warming flux splitting for inviscid terms, the sixth-order central difference for viscous terms, and three stages Runge-Kutta time stepping for the time integration. On the other hand, the finite volume
WENO-GKS is a conservative variable reconstruction based method with the same
WENO reconstruction. But it evaluates a time dependent gas distribution function
along a cell interface, and updates the flow variables inside each control volume by
integrating the flux function along the boundary of the control volume in both space
and time. In order to validate the robustness and accuracy of the schemes, both methods are tested under a wide range of flow conditions: vortex propagation, Mach 3
step problem, and the cavity flow at Reynolds number 3200. Our study shows that
both WENO-SW and WENO-GKS yield quantitatively similar results and agree with
each other very well provided a sufficient grid resolution is used. With the reduction of mesh points, the WENO-GKS behaves to have less numerical dissipation and
present more accurate solutions than those from the WENO-SW in all test cases. For
the Navier-Stokes equations, since the WENO-GKS couples inviscid and viscous terms
in a single flux evaluation, and the WENO-SW uses an operator splitting technique, it
appears that the WENO-SW is more sensitive to the WENO reconstruction and boundary treatment. In terms of efficiency, the finite volume WENO-GKS is about 4 times
slower than the finite difference WENO-SW in two dimensional simulations. The current study clearly shows that besides high-order reconstruction, an accurate gas evolution model or flux function in a high-order scheme is also important in the capturing of physical solutions. In a physical flow, the transport, stress deformation, heat conduction, and viscous heating are all coupled in a single gas evolution process. Therefore,
it is preferred to develop such a scheme with multi-dimensionality, and unified treatment of inviscid and dissipative terms. A high-order scheme does prefer a high-order
gas evolution model. Even with the rapid advances of high-order reconstruction techniques, the first-order dynamics of the Riemann solution becomes the bottleneck for
the further development of high-order schemes. In order to avoid the weakness of the
low order flux function, the development of high-order schemes relies heavily on the
weak solution of the original governing equations for the update of additional degree
of freedom, such as the non-conservative gradients of flow variables, which cannot be
physically valid in discontinuous regions. 相似文献
8.
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. 相似文献
9.
10.
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. 相似文献
11.
Linear Scaling Discontinuous Galerkin Density Matrix Minimization Method with Local Orbital Enriched Finite Element Basis: 1-D Lattice Model System
下载免费PDF全文
![点击此处可从《Communications In Computational Physics》网站下载免费的PDF全文](/ch/ext_images/free.gif)
Tiao Lu Wei Cai Jianguo Xin & Yinglong Guo 《Communications In Computational Physics》2013,14(2):276-300
In the first of a series of papers, we will study a discontinuous Galerkin (DG)
framework for many electron quantum systems. The salient feature of this framework
is the flexibility of using hybrid physics-based local orbitals and accuracy-guaranteed
piecewise polynomial basis in representing the Hamiltonian of the many body system.
Such a flexibility is made possible by using the discontinuous Galerkin method
to approximate the Hamiltonian matrix elements with proper constructions of numerical
DG fluxes at the finite element interfaces. In this paper, we will apply the DG
method to the density matrix minimization formulation, a popular approach in the
density functional theory of many body Schrödinger equations. The density matrix
minimization is to find the minima of the total energy, expressed as a functional of the
density matrix ρ(r,r′), approximated by the proposed enriched basis, together with
two constraints of idempotency and electric neutrality. The idempotency will be handled
with the McWeeny's purification while the neutrality is enforced by imposing the
number of electrons with a penalty method. A conjugate gradient method (a Polak-Ribiere
variant) is used to solve the minimization problem. Finally, the linear-scaling
algorithm and the advantage of using the local orbital enriched finite element basis in
the DG approximations are verified by studying examples of one dimensional lattice
model systems. 相似文献
12.
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. 相似文献
13.
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. 相似文献
14.
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. 相似文献
15.
A Novel Strong-Coupling Pseudo-Spectral Method for Numerical Studies of Two-Layer Turbulent Channel Flows
下载免费PDF全文
![点击此处可从《Communications In Computational Physics》网站下载免费的PDF全文](/ch/ext_images/free.gif)
Tong Wu Zixuan Yang Shizhao Wang & Guowei He 《Communications In Computational Physics》2020,28(3):1133-1146
A new method is proposed to simulate a coupled air-water two-layer turbulent channel flow. A numerically effective dynamic viscosity is implemented to calculate the viscous momentum flux at the interface, leading to a strong-coupling scheme
for the evolution of air and water motions. The direct numerical simulation results
are compared with those in the literature obtained from a weak-coupling scheme. It
is discovered that while the turbulence statistics of the air phase based on the strong- and weak-coupling schemes are close to each other, the results on the water side are
influenced by the coupling approach, especially near the water surface. 相似文献
16.
Implicit Quadrature-Free Direct Reconstruction Method for Efficient Scale-Resolving Simulations
下载免费PDF全文
![点击此处可从《Communications In Computational Physics》网站下载免费的PDF全文](/ch/ext_images/free.gif)
The present study develops implicit physical domain-based discontinuous Galerkin (DG) methods for efficient scale-resolving simulations on mixed-curved
meshes. Implicit methods are essential to handle stiff systems in many scale-resolving
simulations of interests in computational science and engineering. The physical
domain-based DG method can achieve high-order accuracy using the optimal bases set
and preserve the required accuracy on non-affine meshes. When using the quadrature-based DG method, these advantages are overshadowed by severe computational costs
on mixed-curved meshes, making implicit scale-resolving simulations unaffordable.
To address this issue, the quadrature-free direct reconstruction method (DRM) is extended to the implicit DG method. In this approach, the generalized reconstruction
approximates non-linear flux functions directly in the physical domain, making the
computing-intensive numerical integrations precomputable at a preprocessing step.
The DRM operator is applied to the residual computation in the matrix-free method.
The DRM operator can be further extended to the system matrix computation for the
matrix-explicit Krylov subspace method and preconditioning. Finally, the A-stable
Rosenbrock-type Runge–Kutta methods are adopted to achieve high-order accuracy
in time. Extensive verification and validation from the manufactured solution to implicit large eddy simulations are conducted. The computed results confirm that the
proposed method significantly improves computational efficiency compared to the
quadrature-based method while accurately resolving detailed unsteady flow features
that are hardly captured by scale-modeled simulations. 相似文献
17.
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. 相似文献
18.
A hybrid lattice-Boltzmann finite-difference method is presented to simulate incompressible, resistive magnetohydrodynamic (MHD) flows. The lattice Boltzmann equation (LBE) with the Lorentz force term is solved to update the flow field
while the magnetic induction equation is solved using the finite difference method to
calculate the magnetic field. This approach is methodologically intuitive because the
governing equations for MHD are solved in their respective original forms. In addition, the extension to 3-D is straightforward. For validation purposes, this approach
was applied to simulate the Hartmann flow, the Orszag-Tang vortex system (2-D and
3-D) and the magnetic reconnection driven by doubly periodic coalescence instability. The obtained results agree well with analytical solutions and simulation results
available in the literature. 相似文献
19.
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. 相似文献
20.
Guoliang Zhang & Tao Xiong 《Communications In Computational Physics》2022,32(1):126-155
We propose a high order finite difference linear scheme combined with a
high order bound preserving maximum-principle-preserving (MPP) flux limiter to
solve the incompressible flow system. For such problem with highly oscillatory structure but not strong shocks, our approach seems to be less dissipative and much less
costly than a WENO type scheme, and has high resolution due to a Hermite reconstruction. Spurious numerical oscillations can be controlled by the weak MPP flux
limiter. Numerical tests are performed for the Vlasov-Poisson system, the 2D guiding-center model and the incompressible Euler system. The comparison between the linear
and WENO type schemes, with and without the MPP flux limiter, will demonstrate the
good performance of our proposed approach. 相似文献