Discrete singular convolution method with perfectly matched absorbing layers for the wave scattering by periodic structures

A new computational algorithm is introduced for solving scattering problem in periodic structure. The PML technique is used to deal with the difficulty on truncating the unbounded domain while the DSC algorithm is utilized for the spatial discretization. The present study reveals that the method is efficient for solving the problem.     

New constructions of perfectly matched layers for the linearized Euler equations

Based on a PML for the advective wave equation, we propose two PML models for the linearized Euler equations. The derivation of the first model can be applied to other physical models. The second model was implemented. Numerical results are shown. To cite this article: F. Nataf, C. R. Acad. Sci. Paris, Ser. I 340 (2005).     

Development of discontinuous Galerkin methods for Maxwell’s equations in metamaterials and perfectly matched layers

The discontinuous Galerkin method has proved to be an accurate and efficient way to numerically solve many differential equations. In this paper, we extend this method to solve the time-dependent Maxwell’s equations when metamaterials and perfectly matched layers are involved. Numerical results are presented to demonstrate that our method is not only simple to implement, but also quite effective in solving Maxwell’s equations in complex media.     

Studies on some perfectly matched layers for one-dimensional time-dependent systems

We will analyze some perfectly matched layers (PMLs) for the one-dimensional time-dependent Maxwell system, acoustic equations and hyperbolic systems in unbounded domains. The exponential decays and convergence of the PML solutions are studied. Some finite difference schemes are proposed for the PML equations and their stability and convergence are established. The work of Y. Lin was fully supported by the NSERC of Canada, of K. Zhang partially supported by a Direct Grant of CUHK (2060276) and NNSF (No. 10701039 of China), whereas the work of J. Zou was fully supported by Hong Kong RGC grants (Project 404606 and Project 404407).     

Convergence analysis of some first order and second order time accurate gradient schemes for semilinear second order hyperbolic equations

This work is devoted to the convergence analysis of finite volume schemes for a model of semilinear second order hyperbolic equations. The model includes for instance the so‐called Sine‐Gordon equation which appears for instance in Solid Physics (cf. Fang and Li, Adv Math (China) 42 (2013), 441–457; Liu et al., Numer Methods Partial Differ Equ 31 (2015), 670–690). We are motivated by two works. The first one is Eymard et al. (IMA J Numer Anal 30 (2010), 1009–1043) where a recent class of nonconforming finite volume meshes is introduced. The second one is Eymard et al. (Numer Math 82 (1999), 91–116) where a convergence of a finite volume scheme for semilinear elliptic equations is provided. The mesh considered in Eymard et al. (Numer Math 82 (1999), 91–116) is admissible in the sense of Eymard et al. (Elsevier, Amsterdam, 2000, 723–1020) and a convergence of a family of approximate solutions toward an exact solution when the mesh size tends to zero is proved. This article is also a continuation of our previous two works (Bradji, Numer Methods Partial Differ Equ 29 (2013), 1278–1321; Bradji, Numer Methods Partial Differ Equ 29 (2013), 1–39) which dealt with the convergence analysis of implicit finite volume schemes for the wave equation. We use as discretization in space the generic spatial mesh introduced in Eymard et al. (IMA J Numer Anal 30 (2010), 1009–1043), whereas the discretization in time is performed using a uniform mesh. Two finite volume schemes are derived using the discrete gradient of Eymard et al. (IMA J Numer Anal 30 (2010), 1009–1043). The unknowns of these two schemes are the values at the center of the control volumes, at some internal interfaces, and at the mesh points of the time discretization. The first scheme is inspired from the previous work (Bradji, Numer Methods Partial Differ Equ 29 (2013), 1–39), whereas the second one (in which the discretization in time is performed using a Newmark method) is inspired from the work (Bradji, Numer Methods Partial Differ Equ 29 (2013), 1278–1321). Under the assumption that the mesh size of the time discretization is small, we prove the existence and uniqueness of the discrete solutions. If we assume in addition to this that the exact solution is smooth, we derive and prove three error estimates for each scheme. The first error estimate is concerning an estimate for the error between a discrete gradient of the approximate solution and the gradient of the exact solution whereas the second and the third ones are concerning the estimate for the error between the exact solution and the discrete solution in the discrete seminorm of and in the norm of . The convergence rate is proved to be for the first scheme and for the second scheme, where (resp. k) is the mesh size of the spatial (resp. time) discretization. The existence, uniqueness, and convergence results stated above do not require any relation between k and . The analysis presented in this work is also applicable in the gradient schemes framework. © 2016 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 5–33, 2017     

Some temporal second order difference schemes for fractional wave equations

In this article, motivated by Alikhanov's new work (Alikhanov, J Comput Phys 280 (2015), 424–438), some difference schemes are proposed for both one‐dimensional and two‐dimensional time‐fractional wave equations. The obtained schemes can achieve second‐order numerical accuracy both in time and in space. The unconditional convergence and stability of these schemes in the discrete H¹‐norm are proved by the discrete energy method. The spatial compact difference schemes with the results on the convergence and stability are also presented. In addition, the three‐dimensional problem is briefly mentioned. Numerical examples illustrate the efficiency of the proposed schemes. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 970–1001, 2016     

Impulsive stability for systems of second order retarded differential equations

This paper is concerned with systems of impulsive second order delay differential equations. We prove that unstable systems can be stabilized by imposition of impulsive controls. The main tools used are Lyapunov functionals, stability theory and control by impulses.     

On exponential stability of second order delay differential equations

We propose a new method for studying stability of second order delay differential equations. Results we obtained are of the form: the exponential stability of ordinary differential equation implies the exponential stability of the corresponding delay differential equation if the delays are small enough. We estimate this smallness through the coefficients of this delay equation. Examples demonstrate that our tests of the exponential stability are essentially better than the known ones. This method works not only for autonomous equations but also for equations with variable coefficients and delays.     

A second order scheme for the navier-stokes equations: stability and convergence

In this paper, a fully discretized projection method is introduced. It contains a parameter operator. Depending on this operator, we can obtain a first-order scheme, which is appropriate for theoretical analysis, and a second-order scheme, which is more suitable for actual computations. In this method, the boundary conditions of the intermediate velocity field and pressure are not needed. We give the proof of the stability and convergence for the first-order case. For the higher order cases, the proof were different, and we will present it elsewhere.

In a forthcoming article [7], we apply this scheme to the driven-cavity problem and compare it with other schemes     

Boundedness and stability criteria for linear ordinary differential equations of the second order

We establish some correlations for solutions of ordinary differential equations and the imaginary part of the complex solution of the corresponding Riccati equation. On the basis of these correlations and the I. M. Sobol’ theorem we prove some new stability and boundedness criteria for linear equations of the second order.     

Asymptotic behavior and stability of second order neutral delay differential equations

We study the asymptotic behavior of a class of second order neutral delay differential equations by both a spectral projection method and an ordinary differential equation method approach. We discuss the relation of these two methods and illustrate some features using examples. Furthermore, a fixed point method is introduced as a third approach to study the asymptotic behavior. We conclude the paper with an application to a mechanical model of turning processes.     

Oscillation criteria for first and second order forced difference equations with mixed nonlinearities

Some new criteria for the oscillation of certain difference equations with mixed nonlinearities are established. The main tool in the proofs is an inequality due to Hardy, Littlewood, and Pólya.     

On accuracy and unconditional stability of linear multistep methods for second order differential equations

Linear multistep methods for solution of the equationy=f(t, y) are studied by means of the test equationy=–² y, with real. It is shown that the order of accuracy cannot exceed 2 for an unconditionally stable method.This work was supported by the NASA-Ames Research Center, Moffett Field, California, under Interchange No. NCA2-OR745-712, while the author was a visitor at the Computer Science Department, Stanford University, Stanford, California.     

Hyers-Ulam stability of linear differential equations of first order, III

Let X be a complex Banach space and let I=(a,b) be an open interval. In this paper, we will prove the generalized Hyers-Ulam stability of the differential equation ty^′(t)+αy(t)+βtrx₀=0 for the class of continuously differentiable functions , where α, β and r are complex constants and x₀ is an element of X. By applying this result, we also prove the Hyers-Ulam stability of the Euler differential equation of second order.