Energy conservation and super convergence analysis of the EC‐S‐FDTD schemes for Maxwell equations with periodic boundaries

This paper is concerned with new energy analysis of the two dimensional Maxwell's equations and the symmetric energy‐conserved splitting finite difference time domain (EC‐S‐FDTD) method with the periodic boundary (PB) condition. New energy identities of the Maxwell's equations in terms of H¹ and H² norms are proposed and interpreted by considering the physical meanings of the H¹ and H² semi‐norms in the identities. It is found from these new identities that the first and second curls of the electromagnetic fields are conserved in terms their magnitudes. By the energy methods, the numerical energy identities of the symmetric EC‐S‐FDTD method are derived and shown to converge to the continuous energy identities of the Maxwell's equations. This proves that the symmetric EC‐S‐FDTD scheme is unconditionally stable and energy conserved in the discrete H¹ and H² norms. Also by the energy methods, it is proved that the symmetric EC‐S‐FDTD method with PB condition is of second order (super) convergence in the discrete H¹ and H² norms. Numerical experiments are carried out and confirm the analysis on energy conservation, stability and super convergence.     

An analysis for the compressible Stokes equations by first‐order system of least‐squares finite element method

This article applies the first‐order system least‐squares (fosls) finite element method developed by Cai, Manteuffel and McCormick to the compressible Stokes equations. By introducing a new dependent velocity flux variable, we recast the compressible Stokes equations as a first‐order system. Then it is shown that the ellipticity and continuity hold for the least‐squares functionals employing the mixture of H^?1 and L², so that the fosls finite element methods yield best approximations for the velocity flux and velocity. © 2001 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 17:689–699, 2001     

Analysis of newton multilevel stabilized finite volume method for the three‐dimensional stationary Navier‐Stokes equations

This article proposes and analyzes a multilevel stabilized finite volume method(FVM) for the three‐dimensional stationary Navier–Stokes equations approximated by the lowest equal‐order finite element pairs. The method combines the new stabilized FVM with the multilevel discretization under the assumption of the uniqueness condition. The multilevel stabilized FVM consists of solving the nonlinear problem on the coarsest mesh and then performs one Newton correction step on each subsequent mesh thus only solving one large linear systems. The error analysis shows that the multilevel‐stabilized FVM provides an approximate solution with the convergence rate of the same order as the usual stabilized finite element solution solving the stationary Navier–Stokes equations on a fine mesh for an appropriate choice of mesh widths: h_j ～ h_j‐1², j = 1,…,J. Therefore, the multilevel stabilized FVM is more efficient than the standard one‐level‐stabilized FVM. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

Solving the driven flow by linearized elimination with differential quadrature

In the Navier‐Stokes equations the velocity and the pressure are coupled together by the incompressibility condition div u = 0( u = (u,v)^T) which makes the equations difficult to solve numerically. In this article, a new method named linearized elimination of unknowns with differential quadrature method is applied to the Navier‐Stokes equations. The method is of second‐order accuracy in time and of spectral accuracy in space. The numerical results show that our method is of high accuracy, of good convergence with little computational efforts. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2005     

Numerical solution of nonlinear partial quadratic integro‐differential equations of fractional order via hybrid of block‐pulse and parabolic functions

In this paper, an effective numerical approach based on a new two‐dimensional hybrid of parabolic and block‐pulse functions (2D‐PBPFs) is presented for solving nonlinear partial quadratic integro‐differential equations of fractional order. Our approach is based on 2D‐PBPFs operational matrix method together with the fractional integral operator, described in the Riemann–Liouville sense. The main characteristic behind this approach is to reduce such problems to those of solving systems of algebraic equations, which greatly simplifies the problem. By using Newton's iterative method, this system is solved, and the solution of fractional nonlinear partial quadratic integro‐differential equations is achieved. Convergence analysis and an error estimate associated with the proposed method is obtained, and it is proved that the numerical convergence order of the suggested numerical method is O(h³) . The validity and applicability of the method are demonstrated by solving three numerical examples. Numerical examples are presented in the form of tables and graphs to make comparisons with the exact solutions much easier.     

Space‐time spectral collocation method for the one‐dimensional sine‐Gordon equation

In this article, we introduce a new space‐time spectral collocation method for solving the one‐dimensional sine‐Gordon equation. We apply a spectral collocation method for discretizing spatial derivatives, and then use the spectral collocation method for the time integration of the resulting nonlinear second‐order system of ordinary differential equations (ODE). Our formulation has high‐order accurate in both space and time. Optimal a priori error bounds are derived in the L²‐norm for the semidiscrete formulation. Numerical experiments show that our formulation have exponential rates of convergence in both space and time. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 670–690, 2015     

Approximation of the incompressible non‐Newtonian fluid equations by the artificial compressibility method

This paper studies the approximation of the non‐Newtonian fluid equations by the artificial compressibility method. We first introduce a family of perturbed compressible non‐Newtonian fluid equations (depending on a positive parameter ε) that approximates the incompressible equations as ε → 0⁺. Then, we prove the unique existence and convergence of solutions for the compressible equations to the solutions of the incompressible equations. Copyright © 2012 John Wiley & Sons, Ltd.     

Orthogonal spline collocation methods for the stream function‐vorticity formulation of the Navier–Stokes equations

An orthogonal spline collocation (OSC) spatial discretization is proposed for the solution of the fully coupled stream function‐vorticity formulation of the Navier–Stokes equations in two dimensions. For the time‐stepping, a three‐level leapfrog scheme is employed. This method is algebraically linear, and, at each time step, gives rise to a system of linear equations of the form arising in the OSC approximation of the biharmonic Dirichlet problem and can be solved by a fast direct method. Error estimates in the H^l–norm in space, l = 1,2, are derived for the semi‐discrete method and the fully‐discrete leapfrog scheme which is also shown to be second order accurate in time. Numerical results are presented which confirm the theoretical analysis and exhibit superconvergence phenomena, which provide superconvergent approximations to the components of the velocity. © John Wiley & Sons, Inc. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2008     

Nonoverlapping domain decomposition characteristic finite differences for three‐dimensional convection‐diffusion equations

One domain decomposition method modified with characteristic differences is presented for non‐periodic three‐dimensional equations by multiply‐type quadratic interpolation and variant time‐step technique. This method consists of reduced‐scale, two‐dimensional computation on subdomain interface boundaries and fully implicit subdomain computation in parallel. A computational algorithm is outlined and an error estimate in discrete l²‐ norm is established by introducing new inner products and norms. Finally, numerical examples are given to illustrate the theoretical results, efficiency and parallelism of this method. © 2010 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 28: 17‐37, 2012     

Bivariate spline method for numerical solution of time evolution Navier‐Stokes equations over polygons in stream function formulation

We use a bivariate spline method to solve the time evolution Navier‐Stokes equations numerically. The bivariate splines we use in this article are in the spline space of smoothness r and degree 3r over triangulated quadrangulations. The stream function formulation for the Navier‐Stokes equations is employed. Galerkin's method is applied to discretize the space variables of the nonlinear fourth‐order equation, Crank‐Nicholson's method is applied to discretize the time variable, and Newton's iterative method is then used to solve the resulting nonlinear system. We show the existence and uniqueness of the weak solution in L₂(0, T; H²(Ω)) ∩ L_∞(0, T; H¹(Ω)) of the 2D nonlinear fourth‐order problem and give an estimate of how fast the numerical solution converges to the weak solution. The C¹ cubic splines are implemented in MATLAB for solving the Navier‐Stokes equations numerically. Our numerical experiments show that the method is effective and efficient. © 2003 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 19: 776–827, 2003.     

On the regularity of the axially symmetric solutions to the three‐dimensional MHD equations

We present the improved three‐dimensional axially symmetric incompressible magnetohydrodynamics (MHD) equations with nonzero swirl. We consider three kinds of smooth axially symmetric particular solutions to the MHD equations: (1) u^θ=0,B^r=B^z=0, (2) B^r=B^z=0, and (3) B^θ=0. In particular, we derive new regularity criteria for these three kinds of the three‐dimensional axially symmetric smooth solutions to the MHD equations. Our results also reveal some interesting dynamic behavior of the interaction by the angular vorticity field ω^θ and the angular current density field j^θ. Copyright © 2016 John Wiley & Sons, Ltd.     

Large time behavior for the incompressible magnetohydrodynamic equations in half‐spaces

The weighted L^r‐asymptotic behavior of the strong solution and its first‐order spacial derivatives to the incompressible magnetohydrodynamic (MHD) equations is established in a half‐space. Further, the L^∞‐decay rates of the second‐order spatial derivatives of the strong solution are derived by using the Stokes solution formula and employing a decomposition for the nonlinear terms in MHD equations. Copyright © 2014 John Wiley & Sons, Ltd.     

A new fourth‐order numerical algorithm for a class of three‐dimensional nonlinear evolution equations

In this article, a new compact alternating direction implicit finite difference scheme is derived for solving a class of 3‐D nonlinear evolution equations. By the discrete energy method, it is shown that the new difference scheme has good stability and can attain second‐order accuracy in time and fourth‐order accuracy in space with respect to the discrete H¹ ‐norm. A Richardson extrapolation algorithm is applied to achieve fourth‐order accuracy in temporal dimension. Numerical experiments illustrate the accuracy and efficiency of the extrapolation algorithm. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

Lq‐estimates for the stationary Oseen equations on the exterior of a rotating obstacle

We study the Dirichlet problem for the stationary Oseen equations around a rotating body in an exterior domain. Our main results are the existence and uniqueness of weak and very weak solutions satisfying appropriate L^q‐estimates. The uniqueness of very weak solutions is shown by the method of cut‐off functions with an anisotropic decay. Then our existence result for very weak solutions is deduced by a duality argument from the existence and estimates of strong solutions. From this and interior regularity of very weak solutions, we finally establish the complete D^1,r‐result for weak solutions of the Oseen equations around a rotating body in an exterior domain, where 4/3<r <4. Here, D^1,r is the homogeneous Sobolev space.     

Stabilized mixed finite element methods for the Navier‐Stokes equations with damping

In this paper, the stabilized mixed finite element methods are presented for the Navier‐Stokes equations with damping. The existence and uniqueness of the weak solutions are proven by use of the Brouwer fixed‐point theorem. Then, optimal error estimates for the H¹‐norm and L²‐norm of the velocity and the L²‐norm of the pressure are derived. Moreover, on the basis of the optimal L²‐norm error estimate of the velocity, a stabilized two‐step method is proposed, which is more efficient than the usual stabilized methods. Finally, two numerical examples are implemented to confirm the theoretical analysis.     

Combination of standard Galerkin and subspace methods for the time‐dependent Navier‐Stokes equations with nonsmooth initial data

We consider a combination of the standard Galerkin method and the subspace decomposition methods for the numerical solution of the two‐dimensional time‐dependent incompressible Navier‐Stokes equations with nonsmooth initial data. Because of the poor smoothness of the solution near t = 0, we use the standard Galerkin method for time interval [0, 1] and the subspace decomposition method time interval [1, ∞). The subspace decomposition method is based on the solution into the sum of a low frequency component integrated using a small time step Δt and a high frequency integrated using a larger time step pΔt with p > 1. From the H¹‐stability and L²‐error analysis, we show that the subspace decomposition method can yield a significant gain in computing time. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2009     

Chebyshev‐Legendre spectral method for solving the two‐dimensional vorticity equations with homogeneous Dirichlet conditions

The Chebyshev‐Legendre spectral method for the two‐dimensional vorticity equations is considered. The Legendre Galerkin Chebyshev collocation method is used with the Chebyshev‐Gauss collocation points. The numerical analysis results under the L²‐norm for the Chebyshev‐Legendre method of one‐dimensional case are generalized into that of the two‐dimensional case. The stability and optimal order convergence of the method are proved. Numerical results are given. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2009     

Superconvergence of a stabilized finite element approximation for the Stokes equations using a local coarse mesh L2 projection

This article first recalls the results of a stabilized finite element method based on a local Gauss integration method for the stationary Stokes equations approximated by low equal‐order elements that do not satisfy the inf‐sup condition. Then, we derive general superconvergence results for this stabilized method by using a local coarse mesh L² projection. These supervergence results have three prominent features. First, they are based on a multiscale method defined for any quasi‐uniform mesh. Second, they are derived on the basis of a large sparse, symmetric positive‐definite system of linear equations for the solution of the stationary Stokes problem. Third, the finite elements used fail to satisfy the inf‐sup condition. This article combines the merits of the new stabilized method with that of the L² projection method. This projection method is of practical importance in scientific computation. Finally, a series of numerical experiments are presented to check the theoretical results obtained. © 2010 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 28: 115‐126, 2012     

Global Well‐Posedness of the Three‐Dimensional Primitive Equations with Only Horizontal Viscosity and Diffusion

In this paper, we consider the initial boundary value problem of the three‐dimensional primitive equations for planetary oceanic and atmospheric dynamics with only horizontal eddy viscosity in the horizontal momentum equations and only horizontal diffusion in the temperature equation. Global well‐posedness of the strong solution is established for any H² initial data. An N‐dimensional logarithmic Sobolev embedding inequality, which bounds the L^∞‐norm in terms of the L^q‐norms up to a logarithm of the L^p‐norm for p > N of the first‐order derivatives, and a system version of the classic Grönwall inequality are exploited to establish the required a~priori H² estimates for global regularity.© 2016 Wiley Periodicals, Inc.