Discontinuous Galerkin methods for dispersive and lossy Maxwell's equations and PML boundary conditions

In this paper, we will present a unified formulation of discontinuous Galerkin method (DGM) for Maxwell's equations in linear dispersive and lossy materials of Debye type and in the artificial perfectly matched layer (PML) regions. An auxiliary differential equation (ADE) method is used to handle the frequency-dependent constitutive relations with the help of auxiliary polarization currents in the computational and PML regions. The numerical flux for the dispersive lossy Maxwell's equations with the auxiliary polarization current variables is derived. Various numerical results are provided to validate the proposed formulation.     

Amplitude equations for non-linear Rayleigh waves

We investigate asymptotic equations describing small amplitude surface elastic waves in the half-plane (Rayleigh waves). For hyperelastic materials such model equations are Hamiltonian systems, and are seen to lead to the formation of singularities in the surface elastic displacement. At the time of singularity formation the Fourier spectra of the solutions exhibit power law decay, and the observed exponents suggest the existence of both differentiable and non-differentiable singular profiles.     

An energy-conserving Galerkin scheme for a class of nonlinear dispersive equations

A Galerkin scheme is presented for a class of conservative nonlinear dispersive equations, such as the Camassa–Holm equation and the regularized long wave equation. The scheme has two advantageous features: first, it is conservative in that it keeps the discrete analogue of the continuous energy conservation property in the original equations; second, it can be formulated only with cheap H¹$H^1$-elements even if the original equations include third derivative u_xxx$u_{\mathit{xxx}}$. Numerical experiments confirm the stability and effectiveness of the proposed scheme.     

Efficient stochastic Galerkin methods for random diffusion equations

We discuss in this paper efficient solvers for stochastic diffusion equations in random media. We employ generalized polynomial chaos (gPC) expansion to express the solution in a convergent series and obtain a set of deterministic equations for the expansion coefficients by Galerkin projection. Although the resulting system of diffusion equations are coupled, we show that one can construct fast numerical methods to solve them in a decoupled fashion. The methods are based on separation of the diagonal terms and off-diagonal terms in the matrix of the Galerkin system. We examine properties of this matrix and show that the proposed method is unconditionally stable for unsteady problems and convergent for steady problems with a convergent rate independent of discretization parameters. Numerical examples are provided, for both steady and unsteady random diffusions, to support the analysis.     

Smoothing properties and retarded estimates for some dispersive evolution equations

Smoothing properties, in the form of space-time integrability properties, play an important role in the study of dispersive evolution equations. A number of them follow from a combination of general arguments and specific estimates. We present a general formulation which makes the separation between the two types of ingredients as clear as possible, and we illustrate it with the examples of the Schrödinger equation, of the wave equation, and of a class of 1+1 dimensional equations related to the Benjamin-Ono equation. Of special interest for the Cauchy problem are retarded estimates expressed in terms of those properties. We derive a number of such estimates associated with the last example, and we mention briefly an application of those estimates to the Cauchy problem for the generalized Benjamin-Ono equation.Laboratoire associé au Centre National de la Recherche Scientifique     

A comparison of three perturbation methods for non-linear hyperbolic waves

Simple wave solutions of non-linear hyperbolic equations are studied by using the method of renormalization, the analytic method of characteristics, and the method of multiple scales. It is shown that the results of the method of renormalization depend on whether the potential function or the velocity is normalized. This arbitrariness does not occur when using either the analytic method of characteristics or the method of multiple scales. However, special consideration must be given in determining the potential from the velocity obtained by the analytic method of characteristics. No such consideration is needed when the method of multiple scales is applied. The first term obtained for the potential by the method of multiple scales contains a cumulative term in addition to a non-cumulative term. This first-order term is shown to yield the equal area rule for shock waves, and the slope of an equipotential line is the arithmetic mean of the slope of the characteristic in the unperturbed medium and the slope of the characteristic at the point under consideration.     

Mathematical model for non-linear standing waves in a tube

A mathematical model is developed to describe the behavior of one-dimensional, non-linear standing waves in a fluid-filled, rigid-wall tube. The tube is bounded on one end by a piston vibrating with periodic motion and terminated on the other end by an impedance boundary. A numerical procedure for calculations based on this model also is described. Calculations are carried out on a digital computer. The model and computer program are presently restricted to pre-shock conditions. Numerical results obtained for the special case of a rigid reflector and a sinusoidally moving piston are in excellent agreement with the results of Coppens and Sanders (private communication). The general approach is applicable to any one-dimensional system bounded by a velocity condition on one end and an impedance condition on the other end.     

Hybridizable discontinuous Galerkin methods for partial differential equations in continuum mechanics

We present hybridizable discontinuous Galerkin methods for solving steady and time-dependent partial differential equations (PDEs) in continuum mechanics. The essential ingredients are a local Galerkin projection of the underlying PDEs at the element level onto spaces of polynomials of degree k to parametrize the numerical solution in terms of the numerical trace; a judicious choice of the numerical flux to provide stability and consistency; and a global jump condition that enforces the continuity of the numerical flux to arrive at a global weak formulation in terms of the numerical trace. The HDG methods are fully implicit, high-order accurate and endowed with several unique features which distinguish themselves from other discontinuous Galerkin methods. First, they reduce the globally coupled unknowns to the approximate trace of the solution on element boundaries, thereby leading to a significant reduction in the degrees of freedom. Second, they provide, for smooth viscous-dominated problems, approximations of all the variables which converge with the optimal order of k + 1 in the L²-norm. Third, they possess some superconvergence properties that allow us to define inexpensive element-by-element postprocessing procedures to compute a new approximate solution which may converge with higher order than the original solution. And fourth, they allow for a novel and systematic way for imposing boundary conditions for the total stress, viscous stress, vorticity and pressure which are not naturally associated with the weak formulation of the methods. In addition, they possess other interesting properties for specific problems. Their approximate solution can be postprocessed to yield an exactly divergence-free and H(div)-conforming velocity field for incompressible flows. They do not exhibit volumetric locking for nearly incompressible solids. We provide extensive numerical results to illustrate their distinct characteristics and compare their performance with that of continuous Galerkin methods.     

Hybridizable discontinuous Galerkin methods for the time-harmonic Maxwell’s equations

We present two hybridizable discontinuous Galerkin (HDG) methods for the numerical solution of the time-harmonic Maxwell’s equations. The first HDG method explicitly enforces the divergence-free condition and thus necessitates the introduction of a Lagrange multiplier. It produces a linear system for the degrees of freedom of the approximate traces of both the tangential component of the vector field and the Lagrange multiplier. The second HDG method does not explicitly enforce the divergence-free condition and thus results in a linear system for the degrees of freedom of the approximate trace of the tangential component of the vector field only. For both HDG methods, the approximate vector field converges with the optimal order of k + 1 in the L²-norm, when polynomials of degree k are used to represent all the approximate variables. We propose elementwise postprocessing to obtain a new H^curl-conforming approximate vector field which converges with order k + 1 in the H^curl-norm. We present extensive numerical examples to demonstrate and compare the performance of the HDG methods.     

Discontinuous Galerkin finite element methods for hyperbolic nonconservative partial differential equations

We present space- and space–time discontinuous Galerkin finite element (DGFEM) formulations for systems containing nonconservative products, such as occur in dispersed multiphase flow equations. The main criterium we pose on the weak formulation is that if the system of nonconservative partial differential equations can be transformed into conservative form, then the formulation must reduce to that for conservative systems. Standard DGFEM formulations cannot be applied to nonconservative systems of partial differential equations. We therefore introduce the theory of weak solutions for nonconservative products into the DGFEM formulation leading to the new question how to define the path connecting left and right states across a discontinuity. The effect of different paths on the numerical solution is investigated and found to be small. We also introduce a new numerical flux that is able to deal with nonconservative products. Our scheme is applied to two different systems of partial differential equations. First, we consider the shallow water equations, where topography leads to nonconservative products, in which the known, possibly discontinuous, topography is formally taken as an unknown in the system. Second, we consider a simplification of a depth-averaged two-phase flow model which contains more intrinsic nonconservative products.     

Exact travelling wave solutions to some fifth order dispersive nonlinear wave equations

Exact travelling wave solutions to some nonlinear equations of fifth order derivatives are derived by using some accurate ansatz methods.     

A derivation of the Lorenz equations for some unstable dispersive physical systems

The amplitude equations which are valid in the neighbourhood of the bifurcation point of a class of dispersively unstable physical systems when small dissipation is included are shown to be transformable to the Lorenz equations. There is a strong connection with systems which yield equations solvable by the inverse scattering transform when damping is excluded and spatial variation included. Two examples are given in very brief detail: (1) a two-layer model for baroclinic instability, (2) the laser equations giving rise to the SIT equations.     

Secular-free solution up to third order of a model nonlinear equation for dispersive waves

The perturbation method of Lindstedt is applied to study the non linear effect of a nonlinear equation $$\nabla ^2 {\rm E} - \frac{1}{{c^2 }}\frac{{\partial ^2 {\rm E}}}{{\partial t^2 }} - \frac{{\omega _0^2 }}{{c^2 }}{\rm E} + \frac{{2v}}{{c^2 }}\frac{{\partial {\rm E}}}{{\partial t}} + E^2 \left[ {\frac{{\partial {\rm E}}}{{\partial t}} \times A} \right] = 0,$$ where (A. E)=0 andA,c, ω ₀ andν are constants in space and time. Amplitude dependent frequency shifts and the solution up to third order are derived.     

一类非线性色散方程中的新型奇异孤立波

研究一类非线性色散广义DGH方程的新型奇异孤立波及其Painlevé可积性.利用Painlevé分析发现当对流项强度m=2时广义DGH方程是可积的,这是一个新的可积方程.通过构造新的变量代换以及auto-Backlund变换获得该方程丰富的奇异孤立波解,如紧孤立波（compacton）、尖峰孤立波（peakon）、新型带尖点的双孤立波和带爆破点的双孤立波等. 关键词： 非线性色散方程 可积性 奇异孤立波     

Note on coefficient matrices from stochastic Galerkin methods for random diffusion equations

In a recent work by Xiu and Shen [D. Xiu, J. Shen, Efficient stochastic Galerkin methods for random diffusion equations, J. Comput. Phys. 228 (2009) 266–281], the Galerkin methods are used to solve stochastic diffusion equations in random media, where some properties for the coefficient matrix of the resulting system are provided. They also posed an open question on the properties of the coefficient matrix. In this work, we will provide some results related to the open question.     

Another possible model equation for long waves in nonlinear dispersive systems

In the same spirit in which Benjamin, Bona, and Mahoney modified the Korteweg-de Vries equation (U_x+U_t+UU_x+U_xxx=0) to obtain the so-called BBM equation, U_x+U_t+UU_x?U_xxt=0, we propose a different modification: U_x+U_t+UU_x+U_xtt=0. The advantages in this equation are 1) the system is conservative since it can be derived from the Lagrangian density $\text{L=}\text{1}\text{2}\text{θ}{}_{\mathrm{x}}\text{θ}{}_{\mathrm{t}}\text{+}\text{1}\text{2}\text{θ}{}^2{}_{\mathrm{x}}\text{+}\text{1}\text{6}\text{θ}{}^3{}_{\mathrm{x}}\text{?}\text{1}\text{2}\text{θ}{}^2{}_{\mathrm{xt}}\text{,}\text{where}\text{θ}{}_{\mathrm{x}}\text{≡ U;2)}$ for large wave-numbers |k|, the infinitesimal-wave phase speed falls off like $\text{1}\text{|k|}$, in accord with physical intuition; 3) since the equation is of second order in t, both U and U_t can be independently specified for t = 0. Several conservation laws satisfied by solutions to this equation are given.     

Arbitrary order exactly divergence-free central discontinuous Galerkin methods for ideal MHD equations

Ideal magnetohydrodynamic (MHD) equations consist of a set of nonlinear hyperbolic conservation laws, with a divergence-free constraint on the magnetic field. Neglecting this constraint in the design of computational methods may lead to numerical instability or nonphysical features in solutions. In our recent work [F. Li, L. Xu, S. Yakovlev, Central discontinuous Galerkin methods for ideal MHD equations with the exactly divergence-free magnetic field, Journal of Computational Physics 230 (2011) 4828–4847], second and third order exactly divergence-free central discontinuous Galerkin methods were proposed for ideal MHD equations. In this paper, we further develop such methods with higher order accuracy. The novelty here is that the well-established H(div)-conforming finite element spaces are used in the constrained transport type framework, and the magnetic induction equations are extensively explored in order to extract sufficient information to uniquely reconstruct an exactly divergence-free magnetic field. The overall algorithm is local, and it can be of arbitrary order of accuracy. Numerical examples are presented to demonstrate the performance of the proposed methods especially when they are fourth order accurate.     

Soliton solutions in some non-linear Schrödinger-like equations

Soliton solutions are obtained for a class of non-linear Schrödinger-like equations. The parameters of the soliton solutions are written out explicitly.