Discontinuous Galerkin and multiscale variational schemes for a coupled damped nonlinear system of Schrödinger equations

In this article, we study a streamline diffusion‐based discontinuous Galerkin approximation for the numerical solution of a coupled nonlinear system of Schrödinger equations and extend the resulting method to a multiscale variational scheme. We prove stability estimates and derive optimal convergence rates due to the maximal available regularity of the exact solution. In the weak formulation, to make the underlying bilinear form coercive, it was necessary to supply the equation system with an artificial viscosity term with a small coefficient of order proportional to a power of mesh size. We justify the theory by implementing an example of an application of the time‐dependent Schrödinger equation in the coupled ultrafast laser. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

An efficient numerical approach to solve Schrödinger equations with space fractional derivative

We design and analyze an efficient numerical approach to solve the coupled Schrödinger equations with space‐fractional derivative. The numerical scheme is based on leap‐frog in time direction and Fourier method in spatial direction. The advantage of the numerical scheme is that only a linear equation needs to be solved for each time step size, and we proved that the energy and mass of space‐fractional coupled Schrödinger equations (SFCSEs) are conserved in the case of full‐discrete scheme. Moreover, we also analyze the error estimate of the numerical scheme, and numerical solutions converge with the order in L² norm. Numerical examples are illustrated to verify the theoretical results.     

用Fourier Galerkin方法求解间断特征值问题

本文我们考虑一个带间断系数的特征值问题,使用Fourier G a lerk in方法求解.数值试验表明特征值收敛速度达到三阶,而特征函数收敛速度达到2.5阶.表明此方法对间断系数问题非常有效.     

Convergence of compact ADI method for solving linear Schrödinger equations

A compact ADI scheme of second‐order in time and fourth‐order in space is proposed for solving linear Schrödinger equations with periodic boundary conditions. By using the recently suggested discrete energy method, it is shown that the stable compact ADI method is unconditionally convergent in the maximum norm. Numerical experiments, including the comparisons with the second‐order ADI scheme and the time‐splitting Fourier pseudospectral method, are presented to support the theoretical results and show the effectiveness of our method. © 2011 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2012     

A space–time discontinuous Galerkin method for linear convection-dominated Sobolev equations

This article presents a space–time discontinuous Galerkin (DG) finite element method for linear convection-dominated Sobolev equations. The finite element method has basis functions that are continuous in space and discontinuous in time, and variable spatial meshes and time steps are allowed. In the discrete intervals of time, using properties of the Radau quadrature rule, eliminates the restriction to space–time meshes of convectional space–time Galerkin methods. The existence and uniqueness of the approximate solution are proved. An optimal priori error estimate in L^∞(H¹) is derived. Numerical experiments are presented to confirm theoretical results.     

Spatial and spectral superconvergence of discontinuous Galerkin method for hyperbolic problems

In this paper, we analyze the spatial and spectral superconvergence properties of one-dimensional hyperbolic conservation law by a discontinuous Galerkin (DG) method. The analyses combine classical mathematical arguments with MATLAB experiments. Some properties of the DG schemes are discovered using discrete Fourier analyses: superconvergence of the numerical wave numbers, Radau structure of the X spatial error.     

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.     

Two‐grid mixed finite‐element methods for nonlinear Schrödinger equations

Two‐grid mixed finite element schemes are developed for solving both steady state and unsteady state nonlinear Schrödinger equations. The schemes use discretizations based on a mixed finite‐element method. The two‐grid approach yields iterative procedures for solving the nonlinear discrete equations. The idea is to relegate all of the Newton‐like iterations to grids much coarser than the final one, with no loss in order of accuracy. Numerical tests are performed. © 2010 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 28: 63‐73, 2012     

New analysis and application of fractional order Schrödinger equation using with Atangana–Batogna numerical scheme

In this work, an analytical approximation to the solution of Schrodinger equation has been provided. The fractional derivative used in this equation is the Caputo derivative. The existence and uniqueness conditions of solutions for the proposed model are derived based on the power law. While solving the fractional order Schrodinger equation, Atangana–Batogna numerical method is presented for fractional order equation. We obtain an efficient recurrence relation for solving these kinds of equations. To illustrate the usefulness of the numerical scheme, the numerical simulations are presented. The results show that the numerical scheme is very effective and simple.     

热传导方程的间断Galerkin数值解法

本文对具有周期边界的热传导方程采用间断Galerkin(DG)方法给出数值求解方法,并利用傅里叶分析,对数值解进行L∞-误差估计,以一次分段多项式为例,得到半离散格式的误差估计.     

An application of exp‐function method to approximate general and explicit solutions for nonlinear Schrödinger equations

In this work, we implement a relatively new analytical technique, the exp‐function method, for solving nonlinear equations and absolutely a special form of generalized nonlinear Schrödinger equations, which may contain high‐nonlinear terms. This method can be used as an alternative to obtain analytical and approximate solutions of different types of fractional differential equations, which is applied in engineering mathematics. Some numerical examples are presented to illustrate the efficiency and the reliability of exp‐function method. It is predicted that exp‐function method can be found widely applicable in engineering. © 2010 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 27: 1016–1025, 2011     

On the Dirac delta as initial condition for nonlinear Schrödinger equations

In this article we will study the initial value problem for some Schrödinger equations with Dirac-like initial data and therefore with infinite L² mass, obtaining positive results for subcritical nonlinearities. In the critical case and in one dimension we prove that after some renormalization the corresponding solution has finite energy. This allows us to conclude a stability result in the defocusing setting. These problems are related to the existence of a singular dynamics for Schrödinger maps through the so-called Hasimoto transformation.     

An efficient spectral‐Galerkin method based on a dimension reduction scheme for eigenvalue problems of Schrödinger equations

In this paper, we propose an efficient spectral‐Galerkin method based on a dimension reduction scheme for eigenvalue problems of Schrödinger equations. Firstly, we carry out a truncation from a three‐dimensional unbounded domain to a bounded spherical domain. By using spherical coordinate transformation and spherical harmonic expansion, we transform the original problem into a series of one‐dimensional eigenvalue problem that can be solved effectively. Secondly, we introduce a weighted Sobolev space to treat the singularity in the effective potential. Using the property of orthogonal polynomials in weighted Sobolev space, the error estimate for the approximate eigenvalues and corresponding eigenfunctions are proved. Error estimates show that our numerical method can achieve spectral accuracy for approximate eigenvalues and eigenfunctions. Finally, we give some numerical examples to demonstrate the efficiency of our algorithms and the correctness of the theoretical results.     

Local energy‐ and momentum‐preserving schemes for Klein‐Gordon‐Schrödinger equations and convergence analysis

In this article, we obtain local energy and momentum conservation laws for the Klein‐Gordon‐Schrödinger equations, which are independent of the boundary condition and more essential than the global conservation laws. Based on the rule that the numerical methods should preserve the intrinsic properties as much as possible, we propose local energy‐ and momentum‐preserving schemes for the equations. The merit of the proposed schemes is that the local energy/momentum conservation law is conserved exactly in any time‐space region. With suitable boundary conditions, the schemes will be charge‐ and energy‐/momentum‐preserving. Nonlinear analysis shows LEP schemes are unconditionally stable and the numerical solutions converge to the exact solutions with order . The theoretical properties are verified by numerical experiments. © 2017 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 1329–1351, 2017     

Numerical solutions of nonlinear fractional Schrödinger equations using nonstandard discretizations

In this article, numerical study for both nonlinear space‐fractional Schrödinger equation and the coupled nonlinear space‐fractional Schrödinger system is presented. We offer here the weighted average nonstandard finite difference method (WANSFDM) as a novel numerical technique to study such kinds of partial differential equations. The space fractional derivative is described in the sense of the quantum Riesz‐Feller definition. Stability analysis of the proposed method is studied. To show that this method is reliable and computationally efficient different numerical examples are provided. We expect that the proposed schemes can be applicable to different systems of fractional partial differential equations. © 2016 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 1399–1419, 2017     

Efficient energy‐preserving scheme of the three‐coupled nonlinear Schrödinger equation

An energy‐preserving scheme is proposed for the three‐coupled nonlinear Schrödinger (T‐CNLS) equation. The T‐CNLS equation is rewritten into the classical Hamiltonian form. Then the spatial variable is discretized by using high‐order compact method to convert it into a finite‐dimensional Hamiltonian system. Next, a second‐order averaged vector field (AVF) method is employed in time which results in an energy‐preserving scheme. Some theoretical results such as convergence are investigated. In addition, it provides some numerical examples to illustrate the robustness and reliability of the theoretical results. It also explores the role of the parameters in the model and initial condition on the wave propagation.     

Finite element analysis of parabolic integro-differential equations of Kirchhoff type

The aim of this paper is to study parabolic integro-differential equations of Kirchhoff type. We prove the existence and uniqueness of the solution for this problem via Galerkin method. Semidiscrete formulation for this problem is presented using conforming finite element method. As a consequence of the Ritz–Volterra projection, we derive error estimates for both semidiscrete solution and its time derivative. To find the numerical solution of this class of equations, we develop two different types of numerical schemes, which are based on backward Euler–Galerkin method and Crank–Nicolson–Galerkin method. A priori bounds and convergence estimates in spatial as well as temporal direction of the proposed schemes are established. Finally, we conclude this work by implementing some numerical experiments to confirm our theoretical results.     

An energy‐preserving Crank–Nicolson Galerkin method for Hamiltonian partial differential equations

A semidiscretization based method for solving Hamiltonian partial differential equations is proposed in this article. Our key idea consists of two approaches. First, the underlying equation is discretized in space via a selected finite element method and the Hamiltonian PDE can thus be casted to Hamiltonian ODEs based on the weak formulation of the system. Second, the resulting ordinary differential system is solved by an energy‐preserving integrator. The relay leads to a fully discretized and energy‐preserved scheme. This strategy is fully realized for solving a nonlinear Schrödinger equation through a combination of the Galerkin discretization in space and a Crank–Nicolson scheme in time. The order of convergence of our new method is if the discrete L²‐norm is employed. An error estimate is acquired and analyzed without grid ratio restrictions. Numerical examples are given to further illustrate the conservation and convergence of the energy‐preserving scheme constructed.© 2016 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 1485–1504, 2016     

Weak Galerkin finite element methods for Parabolic equations

A newly developed weak Galerkin method is proposed to solve parabolic equations. This method allows the usage of totally discontinuous functions in approximation space and preserves the energy conservation law. Both continuous and discontinuous time weak Galerkin finite element schemes are developed and analyzed. Optimal‐order error estimates in both H¹ and L² norms are established. Numerical tests are performed and reported. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013