A Fourier spectral-discontinuous Galerkin method for time-dependent 3-D Schrödinger–Poisson equations with discontinuous potentials

In this paper, we propose a high order Fourier spectral-discontinuous Galerkin method for time-dependent Schrödinger–Poisson equations in 3-D spaces. The Fourier spectral Galerkin method is used for the two periodic transverse directions and a high order discontinuous Galerkin method for the longitudinal propagation direction. Such a combination results in a diagonal form for the differential operators along the transverse directions and a flexible method to handle the discontinuous potentials present in quantum heterojunction and supperlattice structures. As the derivative matrices are required for various time integration schemes such as the exponential time differencing and Crank Nicholson methods, explicit derivative matrices of the discontinuous Galerkin method of various orders are derived. Numerical results, using the proposed method with various time integration schemes, are provided to validate the method.     

Flows formed after the passage of a discontinuous wave over a bottom drop

The solvability of the problem of the flow formed after a discontinuous wave has passed over a bottom drop is studied within the framework of the first approximation of shallow water theory. Solutions in which the total energy of the flow is either conserved or lost at the drop are considered. Stable self-similar solutions of five qualitatively different types are derived and their domains of existence are determined in the dimensionless parameter plane.     

High‐order 𝒞1 finite‐element interpolating schemes—Part I: Semi‐Lagrangian linear advection

This paper is devoted to the development of accurate high‐order interpolating schemes for semi‐Lagrangian advection. The characteristic‐Galerkin formulation is obtained by using a semi‐Lagrangian temporal discretization of the total derivative. The semi‐Lagrangian method requires high‐order interpolators for accuracy. A class of ??¹ finite‐element interpolating schemes is developed and two semi‐Lagrangian methods are considered by tracking the feet of the characteristic lines either from the interpolation or from the integration nodes. Numerical stability and analytical results quantifying the amount of artificial viscosity induced by the two methods are presented in the case of the one‐dimensional linear advection equation, based on the modified equation approach. Results of test problems to simulate the linear advection of a cosine hill illustrate the performance of the proposed approach. Copyright © 2007 John Wiley & Sons, Ltd.     

A fundamental solution to the Cauchy problem for a fourth-order pseudoparabolic equation

The Cauchy problem for a fourth-order pseudoparabolic equation describing liquid filtration problems in fissured media, moisture transfer in soil, etc., is studied. Under certain summability and boundedness conditions imposed on the coefficients, the operator of this problem and its adjoint operator are proved to be homeomorphism between certain pairs of Banach spaces. Introduced under the same conditions, the concept of a θ-fundamental solution is introduced, which naturally generalizes the concept of the Riemann function to the equations with discontinuous coefficients; the new concept makes it possible to find an integral form of the solution to a nonhomogeneous problem.     

Flow of viscoelastic fluids between plates rotating about distinct axes

We discuss the flow of BKZ fluids in an orthogonal rheometer. Some analytical results are proved, and numerical solutions are obtained for the Currie model. These solutions show a boundary layer behavior at high Reynolds numbers and the possibility of discontinuous solutions or nonexistence at high Weissenberg numbers.     

Remarks on the links between low‐order DG methods and some finite‐difference schemes for the Stokes problem

In this paper we demonstrate that some well‐known finite‐difference schemes can be interpreted within the framework of the local discontinuous Galerkin (LDG) methods using the low‐order piecewise solenoidal discrete spaces introduced in (SIAM J. Numer. Anal. 1990; 27 (6): 1466–1485). In particular, it appears that it is possible to derive the well‐known MAC scheme using a first‐order Nédélec approximation on rectangular cells. It has been recently interpreted within the framework of the Raviart–Thomas approximation by Kanschat (Int. J. Numer. Meth. Fluids 2007; published online). The two approximations are algebraically equivalent to the MAC scheme, however, they have to be applied on grids that are staggered on a distance h/2 in each direction. This paper also demonstrates that both discretizations allow for the construction of a divergence‐free basis, which yields a linear system with a ‘biharmonic’ conditioning. Both this paper and Kanschat (Int. J. Numer. Meth. Fluids 2007; published online) demonstrate that the LDG framework can be used to generalize some popular finite‐difference schemes to grids that are not parallel to the coordinate axes or that are unstructured. Copyright © 2008 John Wiley & Sons, Ltd.     

The n-order iterative schemes for a nonlinear Kirchhoff-Carrier wave equation associated with the mixed inhomogeneous conditions

In this paper, a high-order iterative scheme is established in order to get a convergent sequence at a rate of order N(N?1) to a local unique weak solution of a nonlinear Kirchhoff-Carrier wave equation associated with mixed nonhomogeneous conditions - the boundary conditions are Dirichlet in one part and Robin in other part of boundary. On the other hand, some numerical results were presented.     

A high‐order element‐based Galerkin method for the barotropic vorticity equation

A high‐order element‐based Galerkin method is developed to solve the non‐divergent barotropic vorticity equation (BVE). The solution process involves solving a conservative transport equation for the vorticity fields and a Poisson equation for the stream function fields. The discontinuous Galerkin method is employed for solving the transport equation and a spectral element method (continuous Galerkin) is used for the Poisson equation. A third‐order strong stability preserving explicit Runge–Kutta scheme is used for time integration. A series of tests have been performed to validate the model, which include the evolution of an idealized tropical cyclone and interaction of dual vortices in close proximity. The numerical convergence study is performed by solving the BVE on the sphere where the analytic solution is known. The test results are consistent with physical observations, and the model exhibits exponential convergence. Copyright © 2008 John Wiley & Sons, Ltd.     

High-Order Leap-Frog Based Discontinuous Galerkin Method for the Time-Domain Maxwell Equations on Non-Conforming Simplicial Meshes

A high-order leap-flog based non-dissipative discontinuous Galerkin timedomain method for solving Maxwell's equations is introduced and analyzed. The proposed method combines a centered approximation for the evaluation of fluxes at the interface between neighboring elements, with a Nth-order leap-frog time scheme. Moreover, the interpolation degree is defined at the element level and the mesh is refined locally in a non-conforming way resulting in arbitrary level hanging nodes. The method is proved to be stable under some CFL-like condition on the time step. The convergence of the semi-discrete approximation to Maxwell's equations is established rigorously and bounds on the global divergence error are provided. Numerical experiments with highorder elements show the potential of the method.     

多维非定常中子迁移方程Galerkin有限元法近似解的收敛性和广义解的存在性

使用Galerkin有限元法研究了多维非定常中子迁移方程，证明了Galerkin有限元法近似解的收敛性和广义解的存在性．