A parallel well-balanced finite volume method for shallow water equations with topography on the cubed-sphere

A finite volume scheme for the global shallow water model on the cubed-sphere mesh is proposed and studied in this paper. The new cell-centered scheme is based on Osher’s Riemann solver together with a high-order spatial reconstruction. On each patch interface of the cubed-sphere only one layer of ghost cells is needed in the scheme and the numerical flux is calculated symmetrically across the interface to ensure the numerical conservation of total mass. The discretization of the topographic term in the equation is properly modified in a well-balanced manner to suppress spurious oscillations when the bottom topography is non-smooth. Numerical results for several test cases including a steady-state nonlinear geostrophic flow and a zonal flow over an isolated mountain are provided to show the flexibility of the scheme. Some parallel implementation details as well as some performance results on a parallel supercomputer with more than one thousand processor cores are also provided.     

Semi-Lagrangian semi-implicit time-splitting scheme for a regional model of the atmosphere

Semi-Lagrangian semi-implicit (SLSI) method is currently one of the most efficient approaches for numerical solution of the atmosphere dynamics equations. In this research we apply splitting techniques in the context of a two-time-level SLSI scheme in order to simplify the treatment of the slow physical modes and optimize the solution of the elliptic equations related to implicit part of the scheme. The performed numerical experiments show the accuracy and computational efficiency of the scheme.     

Global conservative solutions for a model equation for shallow water waves of moderate amplitude

In this paper, we study the continuation of solutions to an equation for surface water waves of moderate amplitude in the shallow water regime beyond wave breaking (in [11], Constantin and Lannes proved that this equation accommodates wave breaking phenomena). Our approach is based on a method proposed by Bressan and Constantin [2]. By introducing a new set of independent and dependent variables, which resolve all singularities due to possible wave breaking, the evolution problem is rewritten as a semilinear system. Local existence of the semilinear system is obtained as fixed points of a contractive transformation. Moreover, this formulation allows one to continue the solution after collision time, giving a global conservative solution where the energy is conserved for almost all times. Finally, returning to the original variables, we obtain a semigroup of global conservative solutions, which depend continuously on the initial data.     

Godunov-type approximation for a general resonant balance law with large data

We consider the Cauchy problem for the 2×2 nonstrictly hyperbolic system

     

Global weak solution of planetary geostrophic equations with inviscid geostrophic balance

A reformulation of the planetary geostrophic equations (PGEs) with the inviscid balance equation is proposed and the existence of global weak solutions is established, provided that the mechanical force satisfies an integral constraint. There is only one prognostic equation for the temperature field, and the velocity field is statically determined by the planetary geostrophic balance combined with the incompressibility condition. Furthermore, the velocity profile can be accurately represented as a function of the temperature gradient. In particular, the vertical velocity depends only on the first-order derivative of the temperature. As a result, the bound for the L∞ (0, t ₁ ; L ²) ∩ L ² (0, t ₁ ; H ¹) norm of the temperature field is sufficient to show the existence of the weak solution.     

High order compact Alternating Direction Implicit method for the generalized sine-Gordon equation

High order compact Alternating Direction Implicit scheme is given for solving the generalized sine-Gordon equation in a two-dimensional rectangular domain. We apply the compact finite difference operators to obtain a fourth order discretization for the second order space derivatives, and we give a linearized three time level algorithm for solving the original nonlinear equation. Error estimate is given by the energy method. Numerical results are provided to verify the accuracy and efficiency of this algorithm.     

High-order finite element methods for time-fractional partial differential equations

The aim of this paper is to develop high-order methods for solving time-fractional partial differential equations. The proposed high-order method is based on high-order finite element method for space and finite difference method for time. Optimal convergence rate O((Δt)^2−α+N^−r) is proved for the (r−1)th-order finite element method (r≥2).     

A modified accelerated monotone iterative method for finite difference reaction-diffusion-convection equations

This paper is concerned with monotone algorithms for the finite difference solutions of a class of nonlinear reaction-diffusion-convection equations with nonlinear boundary conditions. A modified accelerated monotone iterative method is presented to solve the finite difference systems for both the time-dependent problem and its corresponding steady-state problem. This method leads to a simple and yet efficient linear iterative algorithm. It yields two sequences of iterations that converge monotonically from above and below, respectively, to a unique solution of the system. The monotone property of the iterations gives concurrently improving upper and lower bounds for the solution. It is shown that the rate of convergence for the sum of the two sequences is quadratic. Under an additional requirement, quadratic convergence is attained for one of these two sequences. In contrast with the existing accelerated monotone iterative methods, our new method avoids computing local maxima in the construction of these sequences. An application using a model problem gives numerical results that illustrate the effectiveness of the proposed method.     

STABILITY ANALYSIS OF NONLINEAR GALERKIN ＆ GALERKIN METHODS FOR NONLINEAR EVOLUTION EQUATIONS

Abstract. Ogr object in this artlcle is to describe tbe Galerkln scheme and nonlin-eax Galerkin scheme for the approximation of nonlinear evolution equations, and tostudy the stability of these schemes. Spatial discretizatlon can be pedormed by eitherGalerkln spectral method or nonlinear Galerldn spectral method; time discretizatlort isdone hy Euler sin.heine wklch is explicit or implicit in the nonlinear terms. According tothe stability analysis of the above schemes, the stability of nonllneex Galerkln methodis better than that of Galexkln method.     

Existence of the universal attractor for the 3-D viscous primitive equations of large-scale moist atmosphere

In this paper, we consider the initial-boundary value problem for the three-dimensional viscous primitive equations of large-scale moist atmosphere which are used to describe the turbulent behavior of long-term weather prediction and climate changes. By obtaining the existence and uniqueness of global strong solutions for the problem and studying the long-time behavior of strong solutions, we prove the existence of the universal attractor for the dynamical system generated by the primitive equations of large-scale moist atmosphere.     

Solving parabolic and hyperbolic equations by the generalized finite difference method

Classical finite difference schemes are in wide use today for approximately solving partial differential equations of mathematical physics. An evolution of the method of finite differences has been the development of generalized finite difference (GFD) method, that can be applied to irregular grids of points.     

Construction of solutions to the two-dimensional stationary Euler equations by the pseudo-advection method

The finite difference approximation of a nonstationary pseudo-advected vorticity equation is proved to yield generalized solutions to the two-dimensional stationary Euler equations with nonvanishing vorticity. This result is obtained by the simultaneous limiting of lattice scale and time.Received: 15 May 2002     

A numerical scheme for the Maxwell equations in the quasi-static regime with a non-local source

In this paper we study a non-linear evolution equation, based on quasi-static electromagnetic fields, with a non-local field-dependent source. This model occurs in transformer driven active magnetic shielding. We present a numerical scheme for both time and space discretization and prove convergence of this scheme. We also derive the corresponding error estimates.     

Error and extrapolation of a compact LOD method for parabolic differential equations

This paper is concerned with a compact locally one-dimensional (LOD) finite difference method for solving two-dimensional nonhomogeneous parabolic differential equations. An explicit error estimate for the finite difference solution is given in the discrete infinity norm. It is shown that the method has the accuracy of the second-order in time and the fourth-order in space with respect to the discrete infinity norm. A Richardson extrapolation algorithm is developed to make the final computed solution fourth-order accurate in both time and space when the time step equals the spatial mesh size. Numerical results demonstrate the accuracy and the high efficiency of the extrapolation algorithm.     

Global conservative solutions of the generalized hyperelastic-rod wave equation

We prove existence of global and conservative solutions of the Cauchy problem for the nonlinear partial differential equation where f is strictly convex or concave and g is locally uniformly Lipschitz. This includes the Camassa-Holm equation (f(u)=u²/2 and g(u)=κu+u²) as well as the hyperelastic-rod wave equation (f(u)=γu²/2 and g(u)=(3−γ)u²/2) as special cases. It is shown that the problem is well-posed for initial data in H¹(R) if one includes a Radon measure that corresponds to the energy of the system with the initial data. The solution is energy preserving. Stability is proved both with respect to initial data and the functions f and g. The proof uses an equivalent reformulation of the equation in terms of Lagrangian coordinates.     

Fractional evolution Dirac-like equations: Some properties and a discrete Von Neumann-type analysis

A system of fractional evolution equations results from employing the tool of the Fractional Calculus and following the method used by Dirac to obtain his well-known equation from Klein–Gordon’s one. It represents a possible interpolation between Dirac and diffusion and wave equations in one space dimension.     

Least-squares spectral element method for non-linear hyperbolic differential equations

The least-squares spectral element method has been applied to the one-dimensional inviscid Burgers equation which allows for discontinuous solutions. In order to achieve high order accuracy both in space and in time a space–time formulation has been applied. The Burgers equation has been discretized in three different ways: a non-conservative formulation, a conservative system with two variables and two equations: one first order linear PDE and one linearized algebraic equation, and finally a variant on this conservative formulation applied to a direct minimization with a QR-decomposition at elemental level. For all three formulations an h/p-convergence study has been performed and the results are discussed in this paper.     

A self-adaptive moving mesh method for the Camassa-Holm equation

A self-adaptive moving mesh method is proposed for the numerical simulations of the Camassa-Holm equation. It is an integrable scheme in the sense that it possesses the exact N-soliton solution. It is named a self-adaptive moving mesh method, because the non-uniform mesh is driven and adapted automatically by the solution. Once the non-uniform mesh is evolved, the solution is determined by solving a tridiagonal linear system. Due to these two superior features of the method, several test problems give very satisfactory results even if by using a small number of grid points.     

Convergence analysis of the FEM approximation of the first order projection method for incompressible flows with and without the inf-sup condition

In this paper we obtain convergence results for the fully discrete projection method for the numerical approximation of the incompressible Navier–Stokes equations using a finite element approximation for the space discretization. We consider two situations. In the first one, the analysis relies on the satisfaction of the inf-sup condition for the velocity-pressure finite element spaces. After that, we study a fully discrete fractional step method using a Poisson equation for the pressure. In this case the velocity-pressure interpolations do not need to accomplish the inf-sup condition and in fact we consider the case in which equal velocity-pressure interpolation is used. Optimal convergence results in time and space have been obtained in both cases.     

The solution of the two-dimensional sine-Gordon equation using the method of lines

The method of lines is used to transform the initial/boundary-value problem associated with the two-dimensional sine-Gordon equation in two space variables into a second-order initial-value problem. The finite-difference methods are developed by replacing the matrix-exponential term in a recurrence relation with rational approximants. The resulting finite-difference methods are analyzed for local truncation error, stability and convergence. To avoid solving the nonlinear system a predictor–corrector scheme using the explicit method as predictor and the implicit as corrector is applied. Numerical solutions for cases involving the most known from the bibliography line and ring solitons are given.