We consider the Schrodinger equation with a logarithmic nonlinearity and a repulsive harmonic potential. Depending on the parameters of the equation, the solution may or may not be dispersive. When dispersion occurs, it does with an exponential rate in time. To control this, we change the unknown function through a generalized lens transform. This approach neutralizes the possible boundary effects, and could be used in the case of the Schrodinger equation without potential. We then employ standard splitting methods on the new equation via a nonuniform grid, after the logarithmic nonlinearity has been regularized. We also discuss the case of a power nonlinearity and give some results concerning the error estimates of the first-order Lie-Trotter splitting method for both cases of nonlinearities. Finally extensive numerical experiments are reported to investigate the dynamics of the equations.

In this paper, an arbitrary Lagrangian–Eulerian local discontinuous Galerkin (ALE-LDG) method for Hamilton–Jacobi equations will be developed, analyzed and numerically tested. This method is based on the time-dependent approximation space defined on the moving mesh. A priori error estimates will be stated with respect to the L∞ 􏰒(0, T ; L2)norm. In particular, the optimal (k + 1) convergence in one dimension and the suboptimal (k + 1 ) convergence in two dimensions will be proven for the semi-discrete method, when a local Lax–Friedrichs flux and piecewise polynomials of degree k on the reference cell are used. Furthermore, the validity of the geometric conservation law will be proven for the fully-discrete method. Also, the link between the piecewise constant ALE-LDG method and the monotone scheme, which converges to the unique viscosity solution, will be shown. The capability of the method will be demonstrated by a variety of one and two dimensional numerical examples with convex and noneconvex Hamiltonian.

In this paper, we introduce a general approach for proving optimal $L^2$ error estimates for the semi-discrete local discontinuous Galerkin (LDG) methods solving linear high order wave equations. The optimal order of error estimates hold not only for the solution itself but also for the auxiliary variables in the LDG method approximating the various order derivatives of the solution. Several examples including the one-dimensional third order wave equation, one-dimensional fifth order wave equation, and multi-dimensional Schr\"{o}dinger equation are explored to demonstrate this approach. The main idea is to derive energy stability for the various auxiliary variables in the LDG discretization, via using the scheme and its time derivatives with different test functions. Special projections are utilized to eliminate the jump terms at the cell boundaries in the error estimate in order to achieve the optimal order of accuracy.

No abstract uploaded!

We present a new formulaiton of the incompressible Navier-Stokes equation in terms of an auxiliary field that differs from the velocity by a gauge transformation. The gauge freedom allows us to assign simple and specific boundary conditions for both the auxiliary field and the gauge field, thus eliminating the issue of pressure boundary condition in the usual primitive variable formulation. The resulting dynamic and kinematic equations can then be solved by standard methods for heat and Poisson equations. A normal mode analysis suggests that in contrast to the classical projection method, the gauge method does not suffer from the problem of numerical boundary layers. Thus the subtleties on the spatial discretization for the projection methods are removed. Consequently, the projection step in the gauge method can be accomplished by standard Poissin solves. We demonstrate the efficiency and accuracy of the gauge method by several numerical examples, including the flow past cylinder.