共查询到20条相似文献,搜索用时 593 毫秒
1.
Solving coupled nonlinear Schrödinger equations via a direct discontinuous Galerkin method 下载免费PDF全文
In this work, we present the direct discontinuous Galerkin (DDG) method for the one-dimensional coupled nonlinear Schrödinger (CNLS) equation. We prove that the new discontinuous Galerkin method preserves the discrete mass conservations corresponding to the properties of the CNLS system. The ordinary differential equations obtained by the DDG space discretization is solved via a third-order stabilized Runge-Kutta method. Numerical experiments show that the new DDG scheme gives stable and less diffusive results and has excellent long-time numerical behaviors for the CNLS equations. 相似文献
2.
3.
Victorita Dolean Stéphane Lanteri Ronan Perrussel 《Journal of computational physics》2008,227(3):2044-2072
We present here a domain decomposition method for solving the three-dimensional time-harmonic Maxwell equations discretized by a discontinuous Galerkin method. In order to allow the treatment of irregularly shaped geometries, the discontinuous Galerkin method is formulated on unstructured tetrahedral meshes. The domain decomposition strategy takes the form of a Schwarz-type algorithm where a continuity condition on the incoming characteristic variables is imposed at the interfaces between neighboring subdomains. A multifrontal sparse direct solver is used at the subdomain level. The resulting domain decomposition strategy can be viewed as a hybrid iterative/direct solution method for the large, sparse and complex coefficients algebraic system resulting from the discretization of the time-harmonic Maxwell equations by a discontinuous Galerkin method. 相似文献
4.
In this paper we propose a new local discontinuous Galerkin method to directly solve Hamilton–Jacobi equations. The scheme is a natural extension of the monotone scheme. For the linear case with constant coefficients, the method is equivalent to the discontinuous Galerkin method for conservation laws. Thus, stability and error analysis are obtained under the framework of conservation laws. For both convex and noneconvex Hamiltonian, optimal (k + 1)th order of accuracy for smooth solutions are obtained with piecewise kth order polynomial approximations. The scheme is numerically tested on a variety of one and two dimensional problems. The method works well to capture sharp corners (discontinuous derivatives) and have the solution converges to the viscosity solution. 相似文献
5.
6.
Local discontinuous Galerkin method for solving Burgers and coupled Burgers equations 总被引:1,自引:0,他引:1 下载免费PDF全文
In the current work, we extend the local discontinuous Galerkin method to a more general application system. The Burgers and coupled Burgers equations are solved by the local discontinuous Galerkin method. Numerical experiments are given to verify the efficiency and accuracy of our method. Moreover the numerical results show that the method can approximate sharp fronts accurately with minimal oscillation. 相似文献
7.
The discontinuous Petrov-Galerkin method for one-dimensional compressible Euler equations in the Lagrangian coordinate 下载免费PDF全文
In this paper, a Petrov-Galerkin scheme named the Runge-Kutta control volume (RKCV) discontinuous finite element method is constructed to solve the one-dimensional compressible Euler equations in the Lagrangian coordinate. Its advantages include preservation of the local conservation and a high resolution. Compared with the Runge-Kutta discontinuous Galerkin (RKDG) method, the RKCV method is easier to implement. Moreover, the advantages of the RKCV and the Lagrangian methods are combined in the new method. Several numerical examples are given to illustrate the accuracy and the reliability of the algorithm. 相似文献
8.
9.
Time-Domain Numerical Solutions of Maxwell Interface Problems with Discontinuous Electromagnetic Waves 下载免费PDF全文
Ya Zhang Duc Duy Nguyen Kewei Du Jin Xu & Shan Zhao 《advances in applied mathematics and mechanics.》2016,8(3):353-385
This paper is devoted to time domain numerical solutions of two-dimensional
(2D) material interface problems governed by the transverse magnetic
(TM) and transverse electric (TE) Maxwell's equations with discontinuous electromagnetic
solutions. Due to the discontinuity in wave solutions across the interface, the
usual numerical methods will converge slowly or even fail to converge. This calls for
the development of advanced interface treatments for popular Maxwell solvers. We
will investigate such interface treatments by considering two typical Maxwell solvers
– one based on collocation formulation and the other based on Galerkin formulation. To
restore the accuracy reduction of the collocation finite-difference time-domain (FDTD)
algorithm near an interface, the physical jump conditions relating discontinuous wave
solutions on both sides of the interface must be rigorously enforced. For this purpose,
a novel matched interface and boundary (MIB) scheme is proposed in this work, in
which new jump conditions are derived so that the discontinuous and staggered features
of electric and magnetic field components can be accommodated. The resulting
MIB time-domain (MIBTD) scheme satisfies the jump conditions locally and suppresses
the staircase approximation errors completely over the Yee lattices. In the discontinuous
Galerkin time-domain (DGTD) algorithm – a popular Galerkin Maxwell solver, a
proper numerical flux can be designed to accurately capture the jumps in the electromagnetic
waves across the interface and automatically preserves the discontinuity in
the explicit time integration. The DGTD solution to Maxwell interface problems is explored
in this work, by considering a nodal based high order discontinuous Galerkin
method. In benchmark TM and TE tests with analytical solutions, both MIBTD and
DGTD schemes achieve the second order of accuracy in solving circular interfaces. In
comparison, the numerical convergence of the MIBTD method is slightly more uniform,
while the DGTD method is more flexible and robust. 相似文献
10.
A Discontinuous Galerkin Method Based on a BGK Scheme for the Navier-Stokes Equations on Arbitrary Grids 下载免费PDF全文
A discontinuous Galerkin Method based on a Bhatnagar-Gross-Krook
(BGK) formulation is presented for the solution of the compressible
Navier-Stokes equations on arbitrary grids. The idea behind this
approach is to combine the robustness of the BGK scheme with the
accuracy of the DG methods in an effort to develop a more accurate,
efficient, and robust method for numerical simulations of viscous
flows in a wide range of flow regimes. Unlike the traditional
discontinuous Galerkin methods, where a Local Discontinuous Galerkin
(LDG) formulation is usually used to discretize the viscous fluxes
in the Navier-Stokes equations, this DG method uses a BGK scheme to
compute the fluxes which not only couples the convective and
dissipative terms together, but also includes both discontinuous and
continuous representation in the flux evaluation at a cell interface
through a simple hybrid gas distribution function. The developed
method is used to compute a variety of viscous flow problems on
arbitrary grids. The numerical results obtained by this BGKDG method
are extremely promising and encouraging in terms of both accuracy
and robustness, indicating its ability and potential to become not
just a competitive but simply a superior approach than the current
available numerical methods. 相似文献
11.
针对聚合物充填过程中的裹气现象,采用一种有限元(FEM)-间断有限元(DG)耦合算法对其进行数值模拟。对于自由运动界面,采用水平集(Level Set)方法进行捕捉;用XPP(eXtended Pom-Pom)本构模型来描述黏弹性流体的流变行为。采用有限元-间断有限元耦合算法求解统一的流场方程,并采用隐式间断有限元求解XPP本构方程、Level Set及其重新初始化方程。数值结果与文献中的实验结果及模拟结果吻合较好,验证了数值方法的稳定性及准确性。分析带有非规则嵌件型腔内,注射速度及浇口尺寸对裹气现象的影响,裹气容易出现在较高注射速度及较小浇口的情形。 相似文献
12.
In this paper, central discontinuous Galerkin methods are developed for solving ideal magnetohydrodynamic (MHD) equations. The methods are based on the original central discontinuous Galerkin methods designed for hyperbolic conservation laws on overlapping meshes, and use different discretization for magnetic induction equations. The resulting schemes carry many features of standard central discontinuous Galerkin methods such as high order accuracy and being free of exact or approximate Riemann solvers. And more importantly, the numerical magnetic field is exactly divergence-free. Such property, desired in reliable simulations of MHD equations, is achieved by first approximating the normal component of the magnetic field through discretizing induction equations on the mesh skeleton, namely, the element interfaces. And then it is followed by an element-by-element divergence-free reconstruction with the matching accuracy. Numerical examples are presented to demonstrate the high order accuracy and the robustness of the schemes. 相似文献
13.
14.
Jiangxing Wang Ziqing Xie & Chuanmiao Chen 《advances in applied mathematics and mechanics.》2015,7(6):796-817
An implicit discontinuous Galerkin method is introduced to solve the time-domain
Maxwell's equations in metamaterials. The Maxwell's equations in metamaterials
are represented by integral-differential equations. Our scheme is based on discontinuous
Galerkin method in spatial domain and Crank-Nicolson method in temporal
domain. The fully discrete numerical scheme is proved to be unconditionally stable.
When polynomial of degree at most $p$ is used for spatial approximation, our scheme is
verified to converge at a rate of $\mathcal{O}(τ^2+h^{p+1/2})$. Numerical results in both 2D and 3D
are provided to validate our theoretical prediction. 相似文献
15.
考察了电、磁场分量分别基于不同近似函数空间展开的一维和二维Maxwell方程间断元求解方法。结合中心数值通量和电、磁场分量近似函数空间的不同组合,构造了各种间断元算子。通过用这些算子在规则和不规则网格上编码分析一维和二维金属腔的谐振模式,详细考察了算子的收敛和伪解支持性,并据此对基函数进行了优选。算子在时域和频域对谐振模式的计算结果彼此符合良好。优选的Maxwell方程间断元算子不仅同时具备能量守恒和免于伪解的特性,且无需引入辅助变量,为设计高品质Maxwell方程间断元算法和研发相关电磁场模拟软件提供了支撑。 相似文献
16.
We develop a smoothed aggregation-based algebraic multigrid solver for high-order discontinuous Galerkin discretizations of the Poisson problem. Algebraic multigrid is a popular and effective method for solving the sparse linear systems that arise from discretizing partial differential equations. However, high-order discontinuous Galerkin discretizations have proved challenging for algebraic multigrid. The increasing condition number of the matrix and loss of locality in the matrix stencil as p increases, in addition to the effect of weakly enforced Dirichlet boundary conditions all contribute to the challenging algebraic setting. 相似文献
17.
Posteriori Error Estimation for an Interior Penalty Discontinuous Galerkin Method for Maxwell's Equations in Cold Plasma 下载免费PDF全文
Jichun Li 《advances in applied mathematics and mechanics.》2009,1(1):107-124
In this paper, we develop a residual-based a posteriori error
estimator for the time-dependent Maxwell's equations in the cold
plasma. Here we consider a semi-discrete interior penalty
discontinuous Galerkin (DG) method for solving the governing
equations. We provide both the upper bound and lower bound analysis
for the error estimator. This is the first posteriori error analysis
carried out for the Maxwell's equations in dispersive media. 相似文献
18.
19.
20.
构造矩形网格下求解Lagrangian坐标系下气动方程组的单元中心型格式. 空间离散采用控制体积间断Petrov-Galerkin方法,时间离散采用二阶TVD Runge-Kutta方法. 利用限制器来抑制非物理震荡并保证RKCV算法的稳定性. 构造的算法可以保证物理量的局部守恒. 与Runge-Kutta间断Galerkin(RKDG)方法相比较,RKCV方法的计算公式少一项积分项使得计算较简单. 给出一些数值算例验证了算法的可靠性及效率. 相似文献