共查询到20条相似文献,搜索用时 46 毫秒
1.
A Weighted Runge-Kutta Discontinuous Galerkin Method for 3D Acoustic and Elastic Wave-Field Modeling
下载免费PDF全文
![点击此处可从《Communications In Computational Physics》网站下载免费的PDF全文](/ch/ext_images/free.gif)
Numerically solving 3D seismic wave equations is a key requirement for
forward modeling and inversion. Here, we propose a weighted Runge-Kutta discontinuous Galerkin (WRKDG) method for 3D acoustic and elastic wave-field modeling. For this method, the second-order seismic wave equations in 3D heterogeneous anisotropic media are transformed into a first-order hyperbolic system, and
then we use a discontinuous Galerkin (DG) solver based on numerical-flux formulations for spatial discretization. The time discretization is based on an implicit diagonal Runge-Kutta (RK) method and an explicit iterative technique, which avoids
solving a large-scale system of linear equations. In the iterative process, we introduce
a weighting factor. We investigate the numerical stability criteria of the 3D method in
detail for linear and quadratic spatial basis functions. We also present a 3D analysis of
numerical dispersion for the full discrete approximation of acoustic equation, which
demonstrates that the WRKDG method can efficiently suppress numerical dispersion
on coarse grids. Numerical results for several different 3D models including homogeneous and heterogeneous media with isotropic and anisotropic cases show that the 3D
WRKDG method can effectively suppress numerical dispersion and provide accurate
wave-field information on coarse mesh. 相似文献
2.
A Preconditioned Implicit Free-Surface Capture Scheme for Large Density Ratio on Tetrahedral Grids
下载免费PDF全文
![点击此处可从《Communications In Computational Physics》网站下载免费的PDF全文](/ch/ext_images/free.gif)
Xin Lv Qingping Zou D. E. Reeve & Yong Zhao 《Communications In Computational Physics》2012,11(1):215-248
We present a three dimensional preconditioned implicit free-surface capture scheme on tetrahedral grids. The current scheme improves our recently reported method [10] in several aspects. Specifically, we modified the original eigensystem by applying a preconditioning matrix so that the new eigensystem is virtually independent of density ratio, which is typically large for practical two-phase problems. Further, we replaced the explicit multi-stage Runge-Kutta method by a fully implicit Euler integration scheme for the Navier-Stokes (NS) solver and the Volume of Fluids (VOF) equation is now solved with a second order Crank-Nicolson implicit scheme to reduce the numerical diffusion effect. The preconditioned restarted Generalized Minimal RESidual method (GMRES) is then employed to solve the resulting linear system. The validation studies show that with these modifications, the method has improved stability and accuracy when dealing with large density ratio two-phase problems. 相似文献
3.
An Implicit Algorithm for High-Order DG/FV Schemes for Compressible Flows on 2D Arbitrary Grids
下载免费PDF全文
![点击此处可从《Communications In Computational Physics》网站下载免费的PDF全文](/ch/ext_images/free.gif)
Laiping Zhang Ming Li Wei Liu & Xin He 《Communications In Computational Physics》2015,17(1):287-316
A Newton/LU-SGS (lower-upper symmetric Gauss-Seidel) iteration implicit
method was developed to solve two-dimensional Euler and Navier-Stokes equations
by the DG/FV hybrid schemes on arbitrary grids. The Newton iteration was employed
to solve the nonlinear system, while the linear system was solved with LU-SGS iteration.
The effect of several parameters in the implicit scheme, such as the CFL number,
the Newton sub-iteration steps, and the update frequency of Jacobian matrix, was investigated
to evaluate the performance of convergence history. Several typical test
cases were simulated, and compared with the traditional explicit Runge-Kutta (RK)
scheme. Firstly the Couette flow was tested to validate the order of accuracy of the
present DG/FV hybrid schemes. Then a subsonic inviscid flow over a bump in a channel
was simulated and the effect of parameters was alsoinvestigated. Finally, the implicit
algorithm was applied to simulate a subsonic inviscid flow over a circular cylinder
and the viscous flow in a square cavity. The numerical results demonstrated that
the present implicit scheme can accelerate the convergence history efficiently. Choosing
proper parameters would improve the efficiency of the implicit scheme. Moreover,
in the same framework, the DG/FV hybrid schemes are more efficient than the same
order DG schemes. 相似文献
4.
Fourth Order Exponential Time Differencing Method with Local Discontinuous Galerkin Approximation for Coupled Nonlinear Schrödinger Equations
下载免费PDF全文
![点击此处可从《Communications In Computational Physics》网站下载免费的PDF全文](/ch/ext_images/free.gif)
This paper studies a local discontinuous Galerkin method combined with
fourth order exponential time differencing Runge-Kutta time discretization and a fourth
order conservative method for solving the nonlinear Schrödinger equations. Based on
different choices of numerical fluxes, we propose both energy-conserving and energy-dissipative
local discontinuous Galerkin methods, and have proven the error estimates
for the semi-discrete methods applied to linear Schrödinger equation. The numerical
methods are proven to be highly efficient and stable for long-range soliton computations.
Extensive numerical examples are provided to illustrate the accuracy, efficiency
and reliability of the proposed methods. 相似文献
5.
On the Stability and CPU Time of the Implicit Runge-Kutta Schemes for Steady State Simulations
下载免费PDF全文
![点击此处可从《Communications In Computational Physics》网站下载免费的PDF全文](/ch/ext_images/free.gif)
Implicit time integration schemes are popular because their relaxed stability
constraints can result in better computational efficiency. For time-accurate unsteady
simulations, it has been well recognized that the inherent dispersion and dissipation
errors of implicit Runge-Kutta schemes will reduce the computational accuracy for
large time steps. Yet for steady state simulations using the time-dependent governing
equations, these errors are often overlooked because the intermediate solutions are of
less interest. Based on the model equation dy/dt = (µ+iλ)y of scalar convection diffusion
systems, this study examines the stability limits, dispersion and dissipation errors
of four diagonally implicit Runge-Kutta-type schemes on the complex (µ+iλ)∆t
plane. Through numerical experiments, it is shown that, as the time steps increase,
the A-stable implicit schemes may not always have reduced CPU time and the computations
may not always remain stable, due to the inherent dispersion and dissipation
errors of the implicit Runge-Kutta schemes. The dissipation errors may decelerate the
convergence rate, and the dispersion errors may cause large oscillations of the numerical
solutions. These errors, especially those of high wavenumber components, grow
at large time steps. They lead to difficulty in the convergence of the numerical computations,
and result in increasing CPU time or even unstable computations as the time
step increases. It is concluded that an optimal implicit time integration scheme for
steady state simulations should have high dissipation and low dispersion. 相似文献
6.
Development of an Explicit Symplectic Scheme that Optimizes the Dispersion-Relation Equation of the Maxwell's Equations
下载免费PDF全文
![点击此处可从《Communications In Computational Physics》网站下载免费的PDF全文](/ch/ext_images/free.gif)
Tony W. H. Sheu L. Y. Liang & J. H. Li 《Communications In Computational Physics》2013,13(4):1107-1133
In this paper an explicit finite-difference time-domain scheme for solving
the Maxwell's equations in non-staggered grids is presented. The proposed scheme
for solving the Faraday's and Ampère's equations in a theoretical manner is aimed to
preserve discrete zero-divergence for the electric and magnetic fields. The inherent local conservation laws in Maxwell's equations are also preserved discretely all the time
using the explicit second-order accurate symplectic partitioned Runge-Kutta scheme.
The remaining spatial derivative terms in the semi-discretized Faraday's and Ampère's
equations are then discretized to provide an accurate mathematical dispersion relation
equation that governs the numerical angular frequency and the wavenumbers in two
space dimensions. To achieve the goal of getting the best dispersive characteristics, we
propose a fourth-order accurate space centered scheme which minimizes the difference
between the exact and numerical dispersion relation equations. Through the computational exercises, the proposed dual-preserving solver is computationally demonstrated
to be efficient for use to predict the long-term accurate Maxwell's solutions. 相似文献
7.
An Implicit LU-SGS Scheme for the Spectral Volume Method on Unstructured Tetrahedral Grids
下载免费PDF全文
![点击此处可从《Communications In Computational Physics》网站下载免费的PDF全文](/ch/ext_images/free.gif)
Takanori Haga Keisuke Sawada & Z. J. Wang 《Communications In Computational Physics》2009,6(5):978-996
An efficient implicit lower-upper symmetric Gauss-Seidel (LU-SGS) solution
approach has been applied to a high order spectral volume (SV) method for unstructured
tetrahedral grids. The LU-SGS solver is preconditioned by the block element
matrix, and the system of equations is then solved with a LU decomposition.
The compact feature of SV reconstruction facilitates the efficient solution algorithm
even for high order discretizations. The developed implicit solver has shown more
than an order of magnitude of speed-up relative to the Runge-Kutta explicit scheme
for typical inviscid and viscous problems. A convergence to a high order solution for
high Reynolds number transonic flow over a 3D wing with a one equation turbulence
model is also indicated. 相似文献
8.
Kinetic Energy Preserving and Entropy Stable Finite Volume Schemes for Compressible Euler and Navier-Stokes Equations
下载免费PDF全文
![点击此处可从《Communications In Computational Physics》网站下载免费的PDF全文](/ch/ext_images/free.gif)
Praveen Chandrashekar 《Communications In Computational Physics》2013,14(5):1252-1286
Centered numerical fluxes can be constructed for compressible Euler equations
which preserve kinetic energy in the semi-discrete finite volume scheme. The essential
feature is that the momentum flux should be of the form where are any consistent approximations to the
pressure and the mass flux. This scheme thus leaves most terms in the numerical
flux unspecified and various authors have used simple averaging. Here we enforce
approximate or exact entropy consistency which leads to a unique choice of all the
terms in the numerical fluxes. As a consequence, a novel entropy conservative flux that
also preserves kinetic energy for the semi-discrete finite volume scheme has been proposed.
These fluxes are centered and some dissipation has to be added if shocks are
present or if the mesh is coarse. We construct scalar artificial dissipation terms which
are kinetic energy stable and satisfy approximate/exact entropy condition. Secondly,
we use entropy-variable based matrix dissipation flux which leads to kinetic energy
and entropy stable schemes. These schemes are shown to be free of entropy violating
solutions unlike the original Roe scheme. For hypersonic flows a blended scheme is
proposed which gives carbuncle free solutions for blunt body flows. Numerical results
for Euler and Navier-Stokes equations are presented to demonstrate the performance
of the different schemes. 相似文献
9.
Parameter Identification in Uncertain Scalar Conservation Laws Discretized with the Discontinuous Stochastic Galerkin Scheme
下载免费PDF全文
![点击此处可从《Communications In Computational Physics》网站下载免费的PDF全文](/ch/ext_images/free.gif)
Louisa Schlachter & Claudia Totzeck 《Communications In Computational Physics》2020,28(4):1585-1608
We study an identification problem which estimates the parameters of the
underlying random distribution for uncertain scalar conservation laws. The hyperbolic equations are discretized with the so-called discontinuous stochastic Galerkin
method, i.e., using a spatial discontinuous Galerkin scheme and a Multielement stochastic Galerkin ansatz in the random space. We assume an uncertain flux or uncertain
initial conditions and that a data set of an observed solution is given. The uncertainty is assumed to be uniformly distributed on an unknown interval and we focus on
identifying the correct endpoints of this interval. The first-order optimality conditions
from the discontinuous stochastic Galerkin discretization are computed on the time-continuous level. Then, we solve the resulting semi-discrete forward and backward
schemes with the Runge-Kutta method. To illustrate the feasibility of the approach,
we apply the method to a stochastic advection and a stochastic equation of Burgers' type. The results show that the method is able to identify the distribution parameters
of the random variable in the uncertain differential equation even if discontinuities are
present. 相似文献
10.
Simulation of Incompressible Free Surface Flow Using the Volume Preserving Level Set Method
下载免费PDF全文
![点击此处可从《Communications In Computational Physics》网站下载免费的PDF全文](/ch/ext_images/free.gif)
This study aims to develop a numerical scheme in collocated Cartesian grids
to solve the level set equation together with the incompressible two-phase flow equations.
A seventh-order accurate upwinding combined compact difference (UCCD7)
scheme has been developed for the approximation of the first-order spatial derivative
terms shown in the level set equation. Developed scheme has a higher accuracy with a
three-point grid stencil to minimize phase error. To preserve the mass of each phase all
the time, the temporal derivative term in the level set equation is approximated by the
sixth-order accurate symplectic Runge-Kutta (SRK6) scheme. All the simulated results
for the dam-break, Rayleigh-Taylor instability, bubble rising, two-bubble merging, and
milkcrown problems in two and three dimensions agree well with the available numerical
or experimental results. 相似文献
11.
A Numerical Methodology for Enforcing Maximum Principles and the Non-Negative Constraint for Transient Diffusion Equations
下载免费PDF全文
![点击此处可从《Communications In Computational Physics》网站下载免费的PDF全文](/ch/ext_images/free.gif)
K. B. Nakshatrala H. Nagarajan & M. Shabouei 《Communications In Computational Physics》2016,19(1):53-93
Transient diffusion equations arise in many branches of engineering and
applied sciences (e.g., heat transfer and mass transfer), and are parabolic partial differential
equations. It is well-known that these equations satisfy important mathematical
properties like maximum principles and the non-negative constraint, which have implications
in mathematical modeling. However, existing numerical formulations for
these types of equations do not, in general, satisfy maximum principles and the non-negative
constraint. In this paper, we present a methodology for enforcing maximum
principles and the non-negative constraint for transient anisotropic diffusion equation.
The proposed methodology is based on the method of horizontal lines in which
the time is discretized first. This results in solving steady anisotropic diffusion equation
with decay equation at every discrete time-level. We also present other plausible
temporal discretizations, and illustrate their shortcomings in meeting maximum principles
and the non-negative constraint. The proposed methodology can handle general
computational grids with no additional restrictions on the time-step. We illustrate the
performance and accuracy of the proposed methodology using representative numerical
examples. We also perform a numerical convergence analysis of the proposed
methodology. For comparison, we also present the results from the standard single-field
semi-discrete formulation and the results from a popular software package, which
all will violate maximum principles and the non-negative constraint. 相似文献
12.
In this paper, we will develop a fast iterative solver for the system of linear
equations arising from the local discontinuous Galerkin (LDG) spatial discretization
and additive Runge-Kutta (ARK) time marching method for the KdV type equations.
Being implicit in time, the severe time step ($∆t$=$\mathcal{O}(∆x^k)$, with the $k$-th order of the
partial differential equations (PDEs)) restriction for explicit methods will be removed.
The equations at the implicit time level are linear and we demonstrate an efficient,
practical multigrid (MG) method for solving the equations. In particular, we numerically
show the optimal or sub-optimal complexity of the MG solver and a two-level
local mode analysis is used to analyze the convergence behavior of the MG method.
Numerical results for one-dimensional, two-dimensional and three-dimensional cases
are given to illustrate the efficiency and capability of the LDG method coupled with
the multigrid method for solving the KdV type equations. 相似文献
13.
An Energy-Preserving Wavelet Collocation Method for General Multi-Symplectic Formulations of Hamiltonian PDEs
下载免费PDF全文
![点击此处可从《Communications In Computational Physics》网站下载免费的PDF全文](/ch/ext_images/free.gif)
In this paper, we develop a novel energy-preserving wavelet collocation
method for solving general multi-symplectic formulations of Hamiltonian PDEs. Based
on the autocorrelation functions of Daubechies compactly supported scaling functions,
the wavelet collocation method is conducted for spatial discretization. The obtained
semi-discrete system is shown to be a finite-dimensional Hamiltonian system, which
has an energy conservation law. Then, the average vector field method is used for
time integration, which leads to an energy-preserving method for multi-symplectic
Hamiltonian PDEs. The proposed method is illustrated by the nonlinear Schrödinger
equation and the Camassa-Holm equation. Since differentiation matrix obtained by
the wavelet collocation method is a cyclic matrix, we can apply Fast Fourier transform
to solve equations in numerical calculation. Numerical experiments show the high
accuracy, effectiveness and conservation properties of the proposed method. 相似文献
14.
Two-Grid Method for Miscible Displacement Problem by Mixed Finite Element Methods and Mixed Finite Element Method of Characteristics
下载免费PDF全文
![点击此处可从《Communications In Computational Physics》网站下载免费的PDF全文](/ch/ext_images/free.gif)
The miscible displacement of one incompressible fluid by another in a porous
medium is governed by a system of two equations. One is elliptic form equation for
the pressure and the other is parabolic form equation for the concentration of one of
the fluids. Since only the velocity and not the pressure appears explicitly in the concentration
equation, we use a mixed finite element method for the approximation of
the pressure equation and mixed finite element method with characteristics for the
concentration equation. To linearize the mixed-method equations, we use a two-grid
algorithm based on the Newton iteration method for this full discrete scheme problems.
First, we solve the original nonlinear equations on the coarse grid, then, we
solve the linearized problem on the fine grid used Newton iteration once. It is shown
that the coarse grid can be much coarser than the fine grid and achieve asymptotically
optimal approximation as long as the mesh sizes satisfy $h=H^2$ in this paper. Finally,
numerical experiment indicates that two-grid algorithm is very effective. 相似文献
15.
A Simple Solver for the Two-Fluid Plasma Model Based on PseudoSpectral Time-Domain Algorithm
下载免费PDF全文
![点击此处可从《Communications In Computational Physics》网站下载免费的PDF全文](/ch/ext_images/free.gif)
Benoit Morel Remo Giust Kazem Ardaneh & Francois Courvoisier 《Communications In Computational Physics》2021,29(3):955-978
We present a solver of 3D two-fluid plasma model for the simulation of
short-pulse laser interactions with plasma. This solver resolves the equations of the
two-fluid plasma model with ideal gas closure. We also include the Bhatnagar-Gross-Krook collision model. Our solver is based on PseudoSpectral Time-Domain (PSTD)
method to solve Maxwell's curl equations. We use a Strang splitting to integrate Euler
equations with source term: while Euler equations are solved with a composite scheme
mixing Lax-Friedrichs and Lax-Wendroff schemes, the source term is integrated with
a fourth-order Runge-Kutta scheme. This two-fluid plasma model solver is simple to
implement because it only relies on finite difference schemes and Fast Fourier Transforms. It does not require spatially staggered grids. The solver was tested against
several well-known problems of plasma physics. Numerical simulations gave results
in excellent agreement with analytical solutions or with previous results from the literature. 相似文献
16.
Jie Zhou Long Chen Yunqing Huang & Wansheng Wang 《Communications In Computational Physics》2015,17(1):127-145
A two-grid method for solving the Cahn-Hilliard equation is proposed in
this paper. This two-grid method consists of two steps. First, solve the Cahn-Hilliard
equation with an implicit mixed finite element method on a coarse grid. Second, solve
two Poisson equations using multigrid methods on a fine grid. This two-grid method
can also be combined with local mesh refinement to further improve the efficiency. Numerical
results including two and three dimensional cases with linear or quadratic elements
show that this two-grid method can speed up the existing mixed finite method
while keeping the same convergence rate. 相似文献
17.
We propose an a-posteriori error/smoothness indicator for standard semi-discrete
finite volume schemes for systems of conservation laws, based on the numerical
production of entropy. This idea extends previous work by the first author limited
to central finite volume schemes on staggered grids. We prove that the indicator converges
to zero with the same rate of the error of the underlying numerical scheme on
smooth flows under grid refinement. We construct and test an adaptive scheme for
systems of equations in which the mesh is driven by the entropy indicator. The adaptive
scheme uses a single nonuniform grid with a variable timestep. We show how
to implement a second order scheme on such a space-time non uniform grid, preserving
accuracy and conservation properties. We also give an example of a p-adaptive
strategy. 相似文献
18.
A fully discrete discontinuous Galerkin method is introduced for solving
time-dependent Maxwell's equations. Distinguished from the Runge-Kutta discontinuous Galerkin method (RKDG) and the finite element time domain method (FETD), in
our scheme, discontinuous Galerkin methods are used to discretize not only the spatial
domain but also the temporal domain. The proposed numerical scheme is proved to be
unconditionally stable, and a convergent rate $\mathcal{O}((∆t)^{r+1}+h^{k+1/2})$ is established under the $L^2$ -norm when polynomials of degree at most $r$ and $k$ are used for temporal and
spatial approximation, respectively. Numerical results in both 2-D and 3-D are provided to validate the theoretical prediction. An ultra-convergence of order $(∆t)^{2r+1}$ in
time step is observed numerically for the numerical fluxes w.r.t. temporal variable at
the grid points. 相似文献
19.
This paper presents a new finite-volume discretization of a generalised Lattice Boltzmann equation (LBE) on unstructured grids. This equation is the continuum LBE, with the addition of a second order time derivative term (memory), and is derived from a second-order differential form of the semi-discrete Boltzmann equation in its implicit form. The new scheme, named unstructured lattice Boltzmann equation with memory (ULBEM), can be advanced in time with a larger time-step than the previous unstructured LB formulations, and a theoretical demonstration of the improved stability is provided. Taylor vortex simulations show that the viscosity is the same as with standard ULBE and demonstrates that the new scheme improves both stability and accuracy. Model validation is also demonstrated by simulating backward-facing step flow at low and moderate Reynolds numbers, as well as by comparing the reattachment length of the recirculating eddy behind the step against experimental and numerical data available in literature. 相似文献
20.
Remapping-Free Adaptive GRP Method for Multi-Fluid Flows I: One Dimensional Euler Equations
下载免费PDF全文
![点击此处可从《Communications In Computational Physics》网站下载免费的PDF全文](/ch/ext_images/free.gif)
In this paper, a remapping-free adaptive GRP method for one dimensional
(1-D) compressible flows is developed. Based on the framework of finite volume
method, the 1-D Euler equations are discretized on moving volumes and the resulting
numerical fluxes are computed directly by the GRP method. Thus the remapping
process in the earlier adaptive GRP algorithm [17,18] is omitted. By adopting a flexible
moving mesh strategy, this method could be applied for multi-fluid problems. The interface
of two fluids will be kept at the node of computational grids and the GRP solver
is extended at the material interfaces of multi-fluid flows accordingly. Some typical numerical
tests show competitive performances of the new method, especially for contact
discontinuities of one fluid cases and the material interface tracking of multi-fluid
cases. 相似文献