On error estimates for Galerkin spectral discretizations of parabolic problems with nonsmooth initial data

We analyze the Legendre and Chebyshev spectral Galerkin semidiscretizations of a one dimensional homogeneous parabolic problem with nonconstant coefficients. We present error estimates for both smooth and nonsmooth data. In the Chebyshev case a limit in the order of approximation is established. On the contrary, in the Legendre case we find an arbitrary high order of convegence.

     

Upwind difference approximations for degenerate parabolic convection-diffusion equations with a discontinuous coefficient

We analyse approximate solutions generated by an upwind differencescheme (of Engquist–Osher type) for nonlinear degenerateparabolic convection–diffusion equations where the nonlinearconvective flux function has a discontinuous coefficient (x)and the diffusion function A(u) is allowed to be strongly degenerate(the pure hyperbolic case is included in our setup). The mainproblem is obtaining a uniform bound on the total variationof the difference approximation u, which is a manifestationof resonance. To circumvent this analytical problem, we constructa singular mapping (, ·) such that the total variationof the transformed variable z = (, u) can be bounded uniformlyin . This establishes strong L¹ compactness of z and, since(, ·) is invertible, also u. Our singular mapping isnovel in that it incorporates a contribution from the diffusionfunction A(u). We then show that the limit of a converging sequenceof difference approximations is a weak solution as well as satisfyinga Krukov-type entropy inequality. We prove that the diffusionfunction A(u) is Hölder continuous, implying that the constructedweak solution u is continuous in those regions where the diffusionis nondegenerate. Finally, some numerical experiments are presentedand discussed.     

Existence of solutions for a higher-order semilinear parabolic equation with singular initial data

We establish the existence of solutions of the Cauchy problem for a higher-order semilinear parabolic equation by introducing a new majorizing kernel. We also study necessary conditions on the initial data for the existence of local-in-time solutions and identify the strongest singularity of the initial data for the solvability of the Cauchy problem.     

The heat equation with singular nonlinearity and singular initial data

We study the existence, uniqueness and regularity of positive solutions of the parabolic equation ut−Δu=a(x)uq+b(x)up in a bounded domain and with Dirichlet's condition on the boundary. We consider here a∈Lα(Ω), b∈Lβ(Ω) and 0<q?1<p. The initial data u(0)=u₀ is considered in the space Lr(Ω), r?1. In the main result (0<q<1), we assume a,b?0 a.e. in Ω and we assume that u₀?γdΩ for some γ>0. We find a unique solution in the space .     

Error estimates for a semi-implicit numerical scheme solving the Landau-Lifshitz equation with an exchange field

^** Email: Ivan.Cimrak{at}ugent.be We study the Landau–Lifshitz (LL) equation describingthe evolution of spin fields in continuum ferromagnets. We considerthe 3D case when the effective field arising in the LL equationincludes exchange interaction, the most challenging case. Thissetting corresponds to the pure isotropic case without a demagnetizingfield. We derive some regularity results for the exact solutionto the LL with Neumann-type boundary conditions. We modify thenumerical scheme studied by A. Prohl in two dimensions and weprove error estimates for this scheme in three dimensions.     

Asymptotic behavior of the plasma equation with homogeneous dirichlet boundary condition and non-negative initial data

We study the asymptotic behaviour of the plasma equation at its extinction time (The first time T' after which the solution vanishes). We establish the Berryman-Holland result [3] for weak solutions of the equation eliminating the strong regularity assumptions of [3]     

Fractional Crank–Nicolson–Galerkin finite element scheme for the time‐fractional nonlinear diffusion equation

This article presents a finite element scheme with Newton's method for solving the time‐fractional nonlinear diffusion equation. For time discretization, we use the fractional Crank–Nicolson scheme based on backward Euler convolution quadrature. We discuss the existence‐uniqueness results for the fully discrete problem. A new discrete fractional Gronwall type inequality for the backward Euler convolution quadrature is established. A priori error estimate for the fully discrete problem in L²(Ω) norm is derived. Numerical results based on finite element scheme are provided to validate theoretical estimates on time‐fractional nonlinear Fisher equation and Huxley equation.     

Time‐splitting methods with charge conservation for the nonlinear Dirac equation

In this work, four numerical time‐splitting methods are proposed for the (1 + 1)‐dimensional nonlinear Dirac equation. All of these methods (or schemes) are proved to satisfy the charge conservation in the discrete level. To enhance the computation efficiency, the block Thomas algorithm is adopted. Numerical experiments are given to test the accuracy order for these schemes, to simulate numerically the binary collision including two standing waves and two moving solitons, meanwhile, the dynamic properties for the nonlinear Dirac equation are discussed. © 2017 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 1582–1602, 2017     

Sharp ill-posedness result for the periodic Benjamin-Ono equation

We prove the discontinuity for the weak L²(T)-topology of the flow-map associated with the periodic Benjamin-Ono equation. This ensures that this equation is ill-posed in Hs(T) as soon as s<0 and thus completes exactly the well-posedness result obtained in Molinet (2008) [12].     

Maximum norm a posteriori error estimates for a 1D singularly perturbed semilinear reaction-diffusion problem

^{A singularly perturbed semilinear two-point boundary-value problemis discretized on arbitrary non-uniform meshes. We present second-ordermaximum norm a posteriori error estimates that hold true uniformlyin the small parameter. Their application to monitor-functionequidistribution and a posteriori mesh refinement are discussed.Numerical results are presented that support our theoreticalestimates.}     

On cauchy problem of the Benney-Lin equation with low regularity initial data

In this work we prove that the initial value problem of the Benney-Lin equation ut ＋ uxxx ＋ β（uxx ＋ u xxxx） ＋ ηuxxxxx ＋ uux = 0 （x ∈ R, t ≥0 0）, where β 〉 0 and η∈R, is locally well-posed in Sobolev spaces HS（R） for s ≥ -7/5. The method we use to prove this result is the bilinear estimate method initiated by Bourgain.     

Sharp threshold of global existence for non-linear Gross-Pitaevskii equation in RN

^** Email: chenguanggan{at}hotmail.com This paper is concerned with the non-linear Gross–Pitaevskiiequation which describes the attractive Bose–Einsteincondensate under a magnetic trap. By an intricate variationalargument we derive out a sharp threshold of blowing up and globalexistence by applying the potential well argument and the concavitymethod. Furthermore, we answer the question: How small are theinitial data, the global solutions of the Cauchy problem ofthe equation exist for [graphic: see PDF]     

Approximation of the solution and its derivative for the singularly perturbed Black-Scholes equation with nonsmooth initial data

A problem for the black-Scholes equation that arises in financial mathematics is reduced, by a transformation of variables, to the Cauchy problem for a singularly perturbed parabolic equation with the variables x, t and a perturbation parameter ɛ, ɛ ∈ (0, 1]. This problem has several singularities such as the unbounded domain, the piecewise smooth initial function (its first-order derivative in x has a discontinuity of the first kind at the point x = 0), an interior (moving in time) layer generated by the piecewise smooth initial function for small values of the parameter ɛ, etc. In this paper, a grid approximation of the solution and its first-order derivative is studied in a finite domain including the interior layer. On a uniform mesh, using the method of additive splitting of a singularity of the interior layer type, a special difference scheme is constructed that allows us to ɛ-uniformly approximate both the solution to the boundary value problem and its first-order derivative in x with convergence orders close to 1 and 0.5, respectively. The efficiency of the constructed scheme is illustrated by numerical experiments. The text was submitted by the authors in English.     

High accuracy error estimates of a Galerkin finite element method for nonlinear time fractional diffusion equation

In this work, an effective and fast finite element numerical method with high-order accuracy is discussed for solving a nonlinear time fractional diffusion equation. A two-level linearized finite element scheme is constructed and a temporal–spatial error splitting argument is established to split the error into two parts, that is, the temporal error and the spatial error. Based on the regularity of the time discrete system, the temporal error estimate is derived. Using the property of the Ritz projection operator, the spatial error is deduced. Unconditional superclose result in H¹-norm is obtained, with no additional regularity assumption about the exact solution of the problem considered. Then the global superconvergence error estimate is obtained through the interpolated postprocessing technique. In order to reduce storage and computation time, a fast finite element method evaluation scheme for solving the nonlinear time fractional diffusion equation is developed. To confirm the theoretical error analysis, some numerical results are provided.     

Attractors for non-autonomous parabolic problems with singular initial data

In this paper, we study the asymptotic behavior of solutions of non-autonomous parabolic problems with singular initial data. We first establish the well-posedness of the equation when the initial data belongs to Lr(Ω) (1<r<∞) and W1,r(Ω) (1<r<N), respectively. When the initial data belongs to Lr(Ω), we establish the existence of uniform attractors in Lr(Ω) for the family of processes with external forces being translation bounded but not translation compact in . When we consider the existence of uniform attractors in , the solution of equation lacks the higher regularity, so we introduce a new type of solution and prove the existence result. For the long time behavior of solutions of the equation in W1,r(Ω), we only obtain the uniform attracting property in the weak topology.     

Local error estimates of the finite element method for an elliptic problem with a Dirac source term

The solutions of elliptic problems with a Dirac measure right‐hand side are not in dimension and therefore the convergence of the finite element solutions is suboptimal in the ‐norm. In this article, we address the numerical analysis of the finite element method for the Laplace equation with Dirac source term: we consider, in dimension 3, the Dirac measure along a curve and, in dimension 2, the punctual Dirac measure. The study of this problem is motivated by the use of the Dirac measure as a reduced model in physical problems, for which high accuracy of the finite element method at the singularity is not required. We show a quasioptimal convergence in the ‐norm, for on subdomains which exclude the singularity; in the particular case of Lagrange finite elements, an optimal convergence in ‐norm is shown on a family of quasiuniform meshes. Our results are obtained using local Nitsche and Schatz‐type error estimates, a weak version of Aubin‐Nitsche duality lemma and a discrete inf‐sup condition. These theoretical results are confirmed by numerical illustrations.     

WKB method for the Dirac equation with a scalar-vector coupling

We outline a recursive method for obtaining WKB expansions of solutions of the Dirac equation in an external centrally symmetric field with a scalar-vector Lorentz structure of the interaction potentials. We obtain semiclassical formulas for radial functions in the classically allowed and forbidden regions and find conditions for matching them in passing through the turning points. We generalize the Bohr-Sommerfeld quantization rule to the relativistic case where a spin-1/2 particle interacts simultaneously with a scalar and an electrostatic external field. We obtain a general expression in the semiclassical approximation for the width of quasistationary levels, which was earlier known only for barrier-type electrostatic potentials (the Gamow formula). We show that the obtained quantization rule exactly produces the energy spectrum for Coulomb- and oscillatory-type potentials. We use an example of the funnel potential to demonstrate that the proposed version of the WKB method not only extends the possibilities for studying the spectrum of energies and wave functions analytically but also ensures an appropriate accuracy of calculations even for states with n_r 1.Translated from Teoreticheskaya i Matematicheskaya Fizika, Vol. 143, No. 1, pp. 83–111, April, 2005.     

Simultaneous inversion of two initial values for a time‐fractional diffusion‐wave equation

This study is devoted to recovering two initial values for a time‐fractional diffusion‐wave equation from boundary Cauchy data. We provide the uniqueness result for recovering two initial values simultaneously by the method of Laplace transformation and analytic continuation. And then we use a nonstationary iterative Tikhonov regularization method to solve the inverse problem and propose a finite dimensional approximation algorithm to find good approximations to the initial values. Numerical examples in one‐ and two‐dimensional cases are provided to show the effectiveness of the proposed method.