The Crank–Nicolson finite spectral element method and numerical simulations for 2D non-stationary Navier–Stokes equations

In this paper, we first build a semi-discretized Crank–Nicolson (CN) model about time for the two-dimensional (2D) non-stationary Navier–Stokes equations about vorticity–stream functions and discuss the existence, stability, and convergence of the time semi-discretized CN solutions. And then, we build a fully discretized finite spectral element CN (FSECN) model based on the bilinear trigonometric basic functions on quadrilateral elements for the 2D non-stationary Navier–Stokes equations about the vorticity–stream functions and discuss the existence, stability, and convergence of the FSECN solutions. Finally, we utilize two sets of numerical experiments to check out the correctness of theoretical consequences.     

Numerical solution of the heat equation with nonlinear boundary conditions in unbounded domains

The numerical solution of the heat equation in unbounded domains (for a 1D problem‐semi‐infinite line and for a 2D one semi‐infinite strip) is considered. The artificial boundaries are introduced and the exact artificial boundary conditions are derived. The original problems are transformed into problems on finite domains. The space semi‐discretization by finite element method and the full approximation by the implicit‐explicit Euler's method are presented. The solvability of the full discretization schemes is analyzed. Computational examples demonstrate the accuracy and the efficiency of the algorithms. Also, the behavior of blowing up solutions is examined numerically. © 2006 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 23: 379–399, 2007     

An explicit spectral collocation method using nonpolynomial basis functions for the time‐dependent Schrödinger equation

We propose a spectral collocation method for the numerical solution of the time‐dependent Schrödinger equation, where the newly developed nonpolynomial functions in a previous study are used as basis functions. Equipped with the new basis functions, various boundary conditions can be imposed exactly. The preferable semi‐implicit time marching schemes are employed for temporal discretization. Moreover, the new basis functions build in a free parameter λ intrinsically, which can be chosen properly so that the semi‐implicit scheme collapses to an explicit scheme. The method is further applied to linear Schrödinger equation set in unbounded domain. The transparent boundary conditions are constructed for time semidiscrete scheme of the linear Schrödinger equation. We employ spectral collocation method using the new basis functions for the spatial discretization, which allows for the exact imposition of the transparent boundary conditions. Comprehensive numerical tests both in bounded and unbounded domain are performed to demonstrate the attractive features of the proposed method.     

Asymptotic stability of the exact boundary controllability of nodal profile for 1D quasi-linear wave equations

In this paper, we consider the asymptotic stability of the exact boundary controllability of nodal profile for 1D quasi-linear wave equations. First, for 1D quasi-linear hyperbolic systems with zero eigenvalues, we establish the existence and uniqueness of semiglobal classical solution to the one-sided mixed initial-boundary value problem on a semibounded initial axis and discuss the asymptotic behavior of the corresponding solutions under different hypotheses on the initial data. Based on these results, we obtain the asymptotic stability of the exact boundary controllability of nodal profile for 1D quasi-linear wave equations on a semibounded time interval.     

BEM‐FEM coupling for the 1D Klein–Gordon equation

A transmission (bidomain) problem for the one‐dimensional Klein–Gordon equation on an unbounded interval is numerically solved by a boundary element method‐finite element method (BEM‐FEM) coupling procedure. We prove stability and convergence of the proposed method by means of energy arguments. Several numerical results are presented, confirming theoretical results. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 2042–2082, 2014     

A reduced-order extrapolating space–time continuous finite element method for the 2D Sobolev equation

This paper mainly concerns with the order reduction to the coefficient vectors of the classical space–time continuous finite element (STCFE) solutions for a two-dimensional Sobolev equation. The classical STCFE model is first constructed for the governing equation, and the theoretical results of the existence, stability, and convergence are provided for the STCFE solutions. We then employ a proper orthogonal decomposition to develop a reduced-order extrapolating STCFE (ROESTCFE) vector model with the lower dimension, and demonstrate the existence, stability, and convergence for the ROESTCFE solutions by the matrix means, resulting in the very concise and flexible theoretical analysis. Lastly, we examine the effectiveness of the developed ROESTCFE model by several numerical tests. It is shown that the ROESTCFE method is computationally very cheap in actual applications.     

A new reliable algorithm based on the sinc function for the time fractional diffusion equation

This paper deals with the numerical solution of time fractional diffusion equation. In this work, we consider the fractional derivative in the sense of Riemann-Liouville. At first, the time fractional derivative is discretized by integrating both sides of the equation with respect to the time variable and we arrive at a semi–discrete scheme. The stability and convergence of time discretized scheme are proven by using the energy method. Also we show that the convergence order of this scheme is O(τ2?α). Then we use the sinc collocation method to approximate the solution of semi–discrete scheme and show that the problem is reduced to a Sylvester matrix equation. Besides by performing some theorems, the exponential convergence rate of sinc method is illustrated. The numerical experiments are presented to show the excellent behavior and high accuracy of the proposed hybrid method in comparison with some other well known methods.     

Spectral methods based on Hermite functions for linear hyperbolic equations

We consider the approximation by spectral and pseudo‐spectral methods of the solution of the Cauchy problem for a scalar linear hyperbolic equation in one space dimension posed in the whole real line. To deal with the fact that the domain of the equation is unbounded, we use Hermite functions as orthogonal basis. Under certain conditions on the coefficients of the equation, we prove the spectral convergence rate of the approximate solutions for regular initial data in a weighted space related to the Hermite functions. Numerical evidence of this convergence is also presented. © 2011 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2012     

A parareal approach of semi‐linear parabolic equations based on general waveform relaxation

We present a parareal approach of semi‐linear parabolic equations based on general waveform relaxation (WR) at the partial differential equation (PDE) level. An algorithm for initial‐boundary value problem and two algorithms for time‐periodic boundary value problem are constructed. The convergence analysis of three algorithms are provided. The results show that the algorithm for initial‐boundary value problem is superlinearly convergent while both algorithms for the time‐periodic boundary value problem linearly converge to the exact solutions at most. Numerical experiments show that the parareal algorithms based on general WR at the PDE level, compared with the parareal algorithm based on the classical WR at the ordinary differential equations (ODEs) level (the PDEs is discretized into ODEs), require much fewer number of iterations to converge.     

Interaction of faults under slip‐dependent friction. Non‐linear eigenvalue analysis

We analyse the evolution of a system of finite faults by considering the non‐linear eigenvalue problems associated to static and dynamic solutions on unbounded domains. We restrict our investigation to the first eigenvalue (Rayleigh quotient). We point out its physical significance through a stability analysis and we give an efficient numerical algorithm able to compute it together with the corresponding eigenfunction. We consider the anti‐plane shearing on a system of finite faults under a slip‐dependent friction in a linear elastic domain, not necessarily bounded. The static problem is formulated in terms of local minima of the energy functional. We introduce the non‐linear (static) eigenvalue problem and we prove the existence of a first eigenvalue/eigenfunction characterizing the isolated local minima. For the dynamic problem, we discuss the existence of solutions with an exponential growth, to deduce a (dynamic) non‐linear eigenvalue problem. We prove the existence of a first dynamic eigenvalue and we analyse its behaviour with respect to the friction parameter. We deduce a mixed finite element discretization of the non‐linear spectral problem and we give a numerical algorithm to approach the first eigenvalue/eigenfunction. Finally we give some numerical results which include convergence tests, on a single fault and a two‐faults system, and a comparison between the non‐linear spectral results and the time evolution results. Copyright © 2004 John Wiley & Sons, Ltd.     

带有线性外力场的双曲平均曲率流Cauchy问题经典解的生命跨度

本文提出并研究带有线性外力场的双曲平均曲率流,通过凸曲线的支撑函数,导出一个双曲型Monge-Ampère 方程并将其转化成Riemann 不变量满足的拟线性双曲方程组。利用拟线性双曲方程组Cauchy 问题的局部解理论,讨论带有线性外力场的双曲平均曲率流Cauchy 问题经典解的生命跨度（即局部解存在的最大时间区间）。     

Solvability,asymptotic analysis and regularity results for a mixed type interaction problem of acoustic waves and piezoelectric structures

In the paper, we investigate the mixed type transmission problem arising in the model of fluid–solid acoustic interaction when a piezoceramic elastic body (Ω⁺) is embedded in an unbounded fluid domain (Ω^?). The corresponding physical process is described by the boundary‐transmission problem for second‐order partial differential equations. In particular, in the bounded domain Ω⁺, we have a 4×4 dimensional matrix strongly elliptic second‐order partial differential equation, while in the unbounded complement domain Ω^?, we have a scalar Helmholtz equation describing acoustic wave propagation. The physical kinematic and dynamic relations mathematically are described by appropriate boundary and transmission conditions. With the help of the potential method and theory of pseudodifferential equations based on the Wiener–Hopf factorization method, the uniqueness and existence theorems are proved in Sobolev–Slobodetskii spaces. We derive asymptotic expansion of solutions, and on the basis of asymptotic analysis, we establish optimal Hölder smoothness results for solutions. Copyright © 2017 John Wiley & Sons, Ltd.     

Stability analysis of block boundary value methods for neutral pantograph equation

This paper deals with the convergence and stability properties of block boundary value methods (BBVMs) for the neutral pantograph equation. Due to its unbounded time lags and limited computer memory, a change in the independent variable is used to transform a pantograph equation into a non-autonomous differential equation with a constant delay but variable coefficients. It is shown under the classical Lipschitz condition that a BBVM is convergent of order p if the underlying boundary value method is consistent with order p. Furthermore, it is proved under a certain condition that BBVMs can preserve the asymptotic stability of exact solutions for the neutral pantograph equation. Meanwhile, some numerical experiments are given to confirm the main conclusions.     

Initial boundary value problem for the 3D quasilinear hyperbolic equations with nonlinear damping

We investigate global existence and asymptotic behavior of the 3D quasilinear hyperbolic equations with nonlinear damping on a bounded domain with slip boundary condition, which describes the propagation of heat waves for rigid solids at very low temperature, below about 20 K. The global existence and uniqueness of classical solutions are obtained when the initial data are near its equilibrium. Time asymptotically, the internal energy is conjectured to satisfy the porous medium equation and the heat flux obeys the classical Darcy’s-type law. Based on energy estimates, we show that the classical solution converges to steady state exponentially fast in time. Moreover, we also verify that the same is true for the corresponding initial boundary value problem of porous medium equation and thus justifies the validity of Darcy’s-type law in large time.     

Convergence analysis for parallel‐in‐time solution of hyperbolic systems

Parallel‐in‐time algorithms have been successfully employed for reducing time‐to‐solution of a variety of partial differential equations, especially for diffusive (parabolic‐type) equations. A major failing of parallel‐in‐time approaches to date, however, is that most methods show instabilities or poor convergence for hyperbolic problems. This paper focuses on the analysis of the convergence behavior of multigrid methods for the parallel‐in‐time solution of hyperbolic problems. Three analysis tools are considered that differ, in particular, in the treatment of the time dimension: (a) space–time local Fourier analysis, using a Fourier ansatz in space and time; (b) semi‐algebraic mode analysis, coupling standard local Fourier analysis approaches in space with algebraic computation in time; and (c) a two‐level reduction analysis, considering error propagation only on the coarse time grid. In this paper, we show how insights from reduction analysis can be used to improve feasibility of the semi‐algebraic mode analysis, resulting in a tool that offers the best features of both analysis techniques. Following validating numerical results, we investigate what insights the combined analysis framework can offer for two model hyperbolic problems, the linear advection equation in one space dimension and linear elasticity in two space dimensions.     

On the combination of some semi‐discretization methods and boundary integral equations for the numerical solution of initial boundary value problems

We consider initial boundary value problems for the homogeneous differential equation of hyperbolic or parabolic type in the unbounded two‐ or three‐dimensional spatial domain with the homogeneous initial conditions and with Dirichlet or Neumann boundary condition. The numerical solution is realized in two steps. At first using the Laguerre transformation or Rothe's method with respect to the time variable the non‐stationary problem is reduced to the sequence of boundary value problems for the non‐homogeneous Helmholtz equation. Further we construct the special integral representation for solutions and obtain the sequence of boundary integral equations (without volume integrals). For the full‐discretization of integral equations we propose some projection methods.     

Exact boundary controllability for a kind of second‐order quasilinear hyperbolic systems and its applications

In this paper, by means of a constructive method based on the existence and uniqueness of the semi‐global C² solution, we establish the local exact boundary controllability for a kind of second‐order quasilinear hyperbolic systems. As an application, we obtain the one‐sided local exact boundary controllability for the first‐order quasilinear hyperbolic systems of diagonal form with boundary conditions in which the diagonal variables corresponding to the positive eigenvalues and those corresponding to the negative eigenvalues are decoupled. Copyright © 2009 John Wiley & Sons, Ltd.     

Hyperbolic–parabolic singular perturbation for quasilinear equations of Kirchhoff type with weak dissipation

We consider a hyperbolic–parabolic singular perturbation problem for a quasilinear hyperbolic equation of Kirchhoff type with dissipation weak in time. The purpose of this paper is to give time‐decay convergence estimates of the difference between the solutions of the hyperbolic equation above and those of the corresponding parabolic equation, together with the unique existence of the global solutions of the hyperbolic equation above. Copyright © 2009 John Wiley & Sons, Ltd.     

对流扩散方程的混合时间间断时空有限元方法

构造并分析二阶对流扩散方程的混合时间间断时空有限元格式．利用混合有限元方法将二阶方程降阶,利用空间连续而时间允许间断的时空有限元方法离散低阶方程．证明数值解的稳定性、存在唯一性和收敛性．最后通过数值结果验证该算法的有效性和可行性．     

Convergence of a misanthrope process to the entropy solution of 1D problems

We prove the convergence, in some strong sense, of a Markov process called “a misanthrope process” to the entropy weak solution of a one-dimensional scalar nonlinear hyperbolic equation. Such a process may be used for the simulation of traffic flows. The convergence proof relies on the uniqueness of entropy Young measure solutions to the nonlinear hyperbolic equation, which holds for both the bounded and the unbounded cases. In the unbounded case, we also prove an error estimate. Finally, numerical results show how this convergence result may be understood in practical cases.