New energy identities and super convergence analysis of the energy conserved splitting FDTD methods for 3D Maxwell's equations

The energy‐conserved splitting finite‐difference time‐domain (EC‐S‐FDTD) method has recently been proposed to solve the Maxwell equations with second order accuracy while numerically keep the L² energy conservation laws of the equations. In this paper, the EC‐S‐FDTD scheme for the 3D Maxwell equations is proved to be energy‐conserved and unconditionally stable in the discrete H¹ norm. The EC‐S‐FDTD scheme is of second‐order accuracy both in time step and spatial steps, which suggests the super‐convergence of this scheme in the discrete H¹ norm. And the divergence of the electric field of the EC‐S‐FDTD scheme in the discrete L² norm is second‐order accurate. Numerical experiments confirm our theoretical analysis. Copyright © 2012 John Wiley & Sons, Ltd.     

General decay and blowup of solutions for coupled viscoelastic equation of Kirchhoff type with degenerate damping terms

In this work, we consider a nonlinear system of viscoelastic equations of Kirchhoff type with degenerate damping and source terms in a bounded domain. Under suitable assumptions on the initial data, the relaxation functions g_i(i = 1,2) and degenerate damping terms, we obtain global existence of solutions. Then, we prove the general decay result. Finally, we prove the finite time blow‐up result of solutions with negative initial energy. This work generalizes and improves earlier results in the literature.     

Compact difference schemes for heat equation with Neumann boundary conditions (II)

This is the further work on compact finite difference schemes for heat equation with Neumann boundary conditions subsequent to the paper, [Sun, Numer Methods Partial Differential Equations (NMPDE) 25 (2009), 1320–1341]. A different compact difference scheme for the one‐dimensional linear heat equation is developed. Truncation errors of the proposed scheme are O(τ² + h⁴) for interior mesh point approximation and O(τ² + h³) for the boundary condition approximation with the uniform partition. The new obtained scheme is similar to the one given by Liao et al. (NMPDE 22 (2006), 600–616), while the major difference lies in no extension of source terms to outside the computational domain any longer. Compared with ones obtained by Zhao et al. (NMPDE 23 (2007), 949–959) and Dai (NMPDE 27 (2011), 436–446), numerical solutions at all mesh points including two boundary points are computed in our new scheme. The significant advantage of this work is to provide a rigorous analysis of convergence order for the obtained compact difference scheme using discrete energy method. The global accuracy is O(τ² + h⁴) in discrete maximum norm, although the spatial approximation order at the Neumann boundary is one lower than that for interior mesh points. The analytical techniques are important and can be successfully used to solve the open problem presented by Sun (NMPDE 25 (2009), 1320–1341), where analyzed theoretical convergence order of the scheme by Liao et al. (NMPDE 22 (2006), 600–616) is only O(τ² + h^3.5) while the numerical accuracy is O(τ² + h⁴), and convergence order of theoretical analysis for the scheme by Zhao et al. (NMPDE 23 (2007), 949–959) is O(τ² + h^2.5), while the actual numerical accuracy is O(τ² + h³). Following the procedure used for the new obtained difference scheme in this work, convergence orders of these two schemes can be proved rigorously to be O(τ² + h⁴) and O(τ² + h³), respectively. Meanwhile, extension to the case involving the nonlinear reaction term is also discussed, and the global convergence order O(τ² + h⁴) is proved. A compact ADI difference scheme for solving two‐dimensional case is derived. Finally, several examples are given to demonstrate the numerical accuracy of new obtained compact difference schemes. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

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.     

Quasi exponential decay of a finite difference space discretization of the 1-d wave equation by pointwise interior stabilization

We consider the wave equation on an interval of length 1 with an interior damping at ξ. It is well-known that this system is well-posed in the energy space and that its natural energy is dissipative. Moreover, as it was proved in Ammari et al. (Asymptot Anal 28(3–4):215–240, 2001), the exponential decay property of its solution is equivalent to an observability estimate for the corresponding conservative system. In this case, the observability estimate holds if and only if ξ is a rational number with an irreducible fraction x = \fracpq,\xi=\frac{p}{q}, where p is odd, and therefore under this condition, this system is exponentially stable in the energy space. In this work, we are interested in the finite difference space semi-discretization of the above system. As for other problems (Zuazua, SIAM Rev 47(2):197–243, 2005; Tcheugoué Tébou and Zuazua, Adv Comput Math 26:337–365, 2007), we can expect that the exponential decay of this scheme does not hold in general due to high frequency spurious modes. We first show that this is indeed the case. Secondly we show that a filtering of high frequency modes allows to restore a quasi exponential decay of the discrete energy. This last result is based on a uniform interior observability estimate for filtered solutions of the corresponding conservative semi-discrete system.     

On decay and blow-up of the solution for a viscoelastic wave equation with boundary damping and source terms

In this paper we consider the decay and blow-up properties of a viscoelastic wave equation with boundary damping and source terms. We first extend the decay result (for the case of linear damping) obtained by Lu et al. (On a viscoelastic equation with nonlinear boundary damping and source terms: Global existence and decay of the solution, Nonlinear Analysis: Real World Applications 12 (1) (2011), 295-303) to the nonlinear damping case under weaker assumption on the relaxation function g(t). Then, we give an exponential decay result without the relation between g^′(t) and g(t) for the linear damping case, provided that ‖g‖_L¹(0,∞) is small enough. Finally, we establish two blow-up results: one is for certain solutions with nonpositive initial energy as well as positive initial energy for both the linear and nonlinear damping cases, the other is for certain solutions with arbitrarily positive initial energy for the linear damping case.     

Quadri-tilings of the Plane

We introduce quadri-tilings and show that they are in bijection with dimer models on a family of graphs R ^* arising from rhombus tilings. Using two height functions, we interpret a sub-family of all quadri-tilings, called triangular quadri-tilings, as an interface model in dimension 2+2. Assigning “critical" weights to edges of R ^*, we prove an explicit expression, only depending on the local geometry of the graph R ^*, for the minimal free energy per fundamental domain Gibbs measure; this solves a conjecture of Kenyon (Invent Math 150:409–439, 2002). We also show that when edges of R ^* are asymptotically far apart, the probability of their occurrence only depends on this set of edges. Finally, we give an expression for a Gibbs measure on the set of all triangular quadri-tilings whose marginals are the above Gibbs measures, and conjecture it to be that of minimal free energy per fundamental domain.     

A new stability number of the Bresse‐Cattaneo system

In this paper, we consider the Bresse‐Cattaneo system with a frictional damping term and prove some optimal decay results for the L²‐norm of the solution and its higher order derivatives. In fact, we show that there is a completely new stability number δ that controls the decay rate of the solution. To prove our results, we use the energy method in the Fourier space to build some very delicate Lyapunov functionals that give the desired results. We also prove the optimality of the results by using the eigenvalues expansion method. In addition, we show that for the absence of the frictional damping term, the solution of our problem does not decay at all. This result improves some early results     

A Richardson scheme of the decomposition method for solving singularly perturbed parabolic reaction-diffusion equation

For the one-dimensional singularly perturbed parabolic reaction-diffusion equation with a perturbation parameter ɛ, where ɛ ∈ (0, 1], the grid approximation of the Dirichlet problem on a rectangular domain in the (x, t)-plane is examined. For small ɛ, a parabolic boundary layer emerges in a neighborhood of the lateral part of the boundary of this domain. A new approach to the construction of ɛ-uniformly converging difference schemes of higher accuracy is developed for initial boundary value problems. The asymptotic construction technique is used to design the base decomposition scheme within which the regular and singular components of the grid solution are solutions to grid subproblems defined on uniform grids. The base scheme converges ɛ-uniformly in the maximum norm at the rate of O(N ⁻²ln² N + N ₀⁻¹), where N + 1 and N ₀ + 1 are the numbers of nodes in the space and time meshes, respectively. An application of the Richardson extrapolation technique to the base scheme yields a higher order scheme called the Richardson decomposition scheme. This higher order scheme convergesɛ-uniformly at the rate of O(N ⁻⁴ln⁴ N + N ₀⁻²). For fixed values of the parameter, the convergence rate is O(N ⁻⁴ + N ₀⁻²).     

Optimal H2‐error estimates of conservative compact difference scheme for the Zakharov equation in two‐space dimension

In this paper, a conservative compact difference scheme is proposed for the two‐dimensional nonlinear Zakharov equation with periodic boundary condition and initial condition. The proposed scheme not only conserve the mass and energy in the discrete level but also are efficient in practical experiments because the Fast Fourier transform (FFT) can be used to speed up the numerical computation. By using the standard energy method and induction argument, we can establish rigorously the unconditional and optimal H²‐error estimates. Some numerical examples are provided to support our theoretical results and show the accuracy and efficiency of the new scheme.     

Optimal damping of selected eigenfrequencies using dimension reduction

We consider linear vibrational systems described by a system of second‐order differential equations of the form , where M and K are positive definite matrices, representing mass and stiffness, respectively. The damping matrix D is assumed to be positive semidefinite. We are interested in finding an optimal damping matrix that will damp a certain (critical) part of the eigenfrequencies. For this, we use an optimization criterion based on the minimization of the average total energy of the system. This is equivalent to the minimization of the trace of the solution of the corresponding Lyapunov equation AX + XA^T = ?GG^T, where A is the matrix obtained from linearizing the second‐order differential equation, and G depends on the critical part of the eigenfrequencies to be damped. The main result is the efficient approximation and the corresponding error bound for the trace of the solution of the Lyapunov equation obtained through dimension reduction, which includes the influence of the right‐hand side GG^T and allows us to control the accuracy of the trace approximation. This trace approximation yields a very accelerated optimization algorithm for determining the optimal damping. Copyright © 2011 John Wiley & Sons, Ltd.     

Sharp regularity estimates for solutions of the wave equation and their traces with prescribed neumann data

In this paper the regularity properties of second-order hyperbolic equations defined over a rectangular domain Θ with boundary Γ under the action of a Neumann boundary forcing term inL ² (0,T;H ^1/4 (Γ)) are investigated. With this given boundary input, we prove by a cosine operator/functional analytical approach that not only is the solution of the wave equation and its derivatives continuous in time, with their pointwise values in a basic energy space (in the interior of Ω), but also that a trace regularity thereof can be assigned for the solution’s time derivative in an appropriate (negative) Sobolev space. This new-found information on the solution and its traces is crucial in handling a mathematical model derived for a particular fluid/structure interaction system.     

Uniform exponential long time decay for the space semi-discretization of a locally damped wave equation via an artificial numerical viscosity

Summary. We consider the finite-difference space semi-discretization of a locally damped wave equation, the damping being supported in a suitable subset of the domain under consideration, so that the energy of solutions of the damped wave equation decays exponentially to zero as time goes to infinity. The decay rate of the semi-discrete systems turns out to depend on the mesh size h of the discretization and tends to zero as h goes to zero. We prove that adding a suitable vanishing numerical viscosity term leads to a uniform (with respect to the mesh size) exponential decay of the energy of solutions. This numerical viscosity term damps out the high frequency numerical spurious oscillations while the convergence of the scheme towards the original damped wave equation is kept. We discuss this problem in 1D and 2D in the interval and the square respectively. Our method of proof relies on discrete multiplier techniques. Mathematics Subject Classification (1991):65M06     

Numerical simulation for the initial-boundary value problem of the Klein-Gordon-Zakharov equations

In this paper,a new conservative finite difference scheme with a parameter θ is proposed for the initial-boundary problem of the Klein-Gordon-Zakharov (KGZ) equations.Convergence of the numerical solutions are proved with order O(h 2 + τ 2) in the energy norm.Numerical results show that the scheme is accurate and efficient.     

A new difference scheme with high accuracy and absolute stability for solving convection–diffusion equations

In this paper, we use a semi-discrete and a padé approximation method to propose a new difference scheme for solving convection–diffusion problems. The truncation error of the difference scheme is O(h⁴+τ⁵). It is shown through analysis that the scheme is unconditionally stable. Numerical experiments are conducted to test its high accuracy and to compare it with Crank–Nicolson method.     

Local energy decay for the wave equation with a nonlinear time-dependent damping

This article addresses a wave equation on a exterior domain in ? ^d (d odd) with nonlinear time-dependent dissipation. Under a microlocal geometric condition we prove that the decay rates of the local energy functional are obtained by solving a nonlinear non-autonomous differential equation     

A new approximation for the generalized fractional derivative and its application to generalized fractional diffusion equation

In this paper, we derive a fourth order approximation for the generalized fractional derivative that is characterized by a scale function z(t) and a weight function w(t) . Combining the new approximation with compact finite difference method, we develop a numerical scheme for a generalized fractional diffusion problem. The stability and convergence of the numerical scheme are proved by the energy method, and it is shown that the temporal and spatial convergence orders are both 4. Several numerical experiments are provided to illustrate the efficiency of our scheme.     

A numerical method for computing radially symmetric solutions of a dissipative nonlinear modified Klein‐Gordon equation

In this article we develop a finite‐difference scheme to approximate radially symmetric solutions of a dissipative nonlinear modified Klein‐Gordon equation subject to smooth initial conditions ? and ψ in an open sphere D around the origin, with constant internal and external damping coefficients—β and γ, respectively—, and nonlinear term of the form G′(w) = w^p, with p > 1 an odd number. The functions ? and ψ are radially symmetric in D, and ?, ψ, r?, and rψ are assumed to be small at infinity. We prove that our scheme is consistent order ??(Δt²) + ??(Δr²) for G′ identically equal to zero and provide a necessary condition for it to be stable order n. Part of our study will be devoted to compare the physical effects of β and γ. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2005     

and Norms Error Estimates in Finite Element Method for Linear Parabolic Interface Problems

We derive new a priori error estimates for linear parabolic equations with discontinuous coefficients. Due to low global regularity of the solutions the error analysis of the standard finite element method for parabolic problems is difficult to adopt for parabolic interface problems. A finite element procedure is, therefore, proposed and analyzed in this paper. We are able to show that the standard energy technique of finite element method for non-interface parabolic problems can be extended to parabolic interface problems if we allow interface triangles to be curved triangles. Optimal pointwise-in-time error estimates in the L ²(Ω) and H ¹(Ω) norms are shown to hold for the semidiscrete scheme. A fully discrete scheme based on backward Euler method is analyzed and pointwise-in-time error estimates are derived. The interfaces are assumed to be arbitrary shape but smooth for our purpose.     

A posteriori error estimates for finite element approximations to the wave equation with discontinuous coefficients

We derive residual‐based a posteriori error estimates of finite element method for linear wave equation with discontinuous coefficients in a two‐dimensional convex polygonal domain. A posteriori error estimates for both the space‐discrete case and for implicit fully discrete scheme are discussed in L^∞(L²) norm. The main ingredients used in deriving a posteriori estimates are new Clément type interpolation estimates in conjunction with appropriate adaption of the elliptic reconstruction technique of continuous and discrete solutions. We use only an energy argument to establish a posteriori error estimates with optimal order convergence in the L^∞(L²) norm.