We propose efficient and accurate numerical methods for computing the ground state and dynamics of the dipolar Bose-Einstein condensates utilising a newly developed dipole-dipole interaction (DDI) solver that is implemented with the non-uniform fast Fourier transform (NUFFT) algorithm. We begin with the three-dimensional (3D) Gross-Pitaevskii equation (GPE) with a DDI term and present the corresponding two-dimensional (2D) model under a strongly anisotropic confining potential. Different from existing methods, the NUFFT based DDI solver removes the singularity by adopting the spherical/polar coordinates in the Fourier space in 3D/2D, respectively, thus it can achieve spectral accuracy in space and simultaneously maintain high efficiency by making full use of FFT and NUFFT whenever it is necessary and/or needed. Then, we incorporate this solver into existing successful methods for computing the ground state and dynamics of GPE with a DDI for dipolar BEC. Extensive numerical comparisons with existing methods are carried out for computing the DDI, ground states and dynamics of the dipolar BEC. Numerical results show that our new methods outperform existing methods in terms of both accuracy and efficiency.

We propose an entropy stable high-resolution finite volume scheme to approximate systems of two-dimensional symmetrizable conservation laws on unstructured grids. In particular we consider Euler equations governing compressible flows. The scheme is constructed using a combination of entropy conservative fluxes and entropy-stable numerical dissipation operators. High resolution is achieved based on a linear reconstruction procedure satisfying a suitable sign property that helps to maintain entropy stability. The proposed scheme is demonstrated to robustly approximate complex flow features by a series of benchmark numerical experiments.

The purpose of this paper is to numerically realize the inverse scattering scheme proposed in [19] of reconstructing complex elastic objects by a single far-field measurement. The unknown elastic scatterers might consist of both rigid bodies and traction-free cavities with components of multiscale sizes presented simultaneously. We conduct extensive numerical experiments to show the effectiveness and efficiency of the imaging scheme proposed in [19]. Moreover, we develop a two-stage technique, which can significantly speed up the reconstruction to yield a fast imaging scheme.

In this paper, an approach combining the DG method in space with CG method in time (CG-DG method) is developed to solve time-dependent Maxwell's equations when meta-materials are involved. Both the unconditional L2-stability and error estimate of order are obtained when polynomials of degree at most r is used for the temporal discretization and at most k for the spatial discretization. Numerical results in 3D are given to validate the theoretical results.

This paper explores the discrete singular convolution method for Hamiltonian PDEs. The differential matrices corresponding to two delta type kernels of the discrete singular convolution are presented analytically, which have the properties of high-order accuracy, bandlimited structure and thus can be excellent candidates for the spatial discretizations for Hamiltonian PDEs. Taking the nonlinear Schrödinger equation and the coupled Schrödinger equations for example, we construct two symplectic integrators combining this kind of differential matrices and appropriate symplectic time integrations, which both have been proved to satisfy the square conservation laws. Comprehensive numerical experiments including comparisons with the central finite difference method, the Fourier pseudospectral method, the wavelet collocation method are given to show the advantages of the new type of symplectic integrators.

The purpose of this paper is to develop and test novel invariant-preserving finite difference schemes for both the Camassa-Holm (CH) equation and one of its 2-component generalizations (2CH). The considered PDEs are strongly nonlinear, admitting soliton-like peakon solutions which are characterized by a slope discontinuity at the peak in the wave shape, and therefore suitable for modeling both short wave breaking and long wave propagation phenomena. The proposed numerical schemes are shown to preserve two invariants, momentum and energy, hence numerically producing wave solutions with smaller phase error over a long time period than those generated by other conventional methods. We first apply the scheme to the CH equation and showcase the merits of considering such a scheme under a wide class of initial data. We then generalize this scheme to the 2CH equation and test this scheme under several types of initial data.

We study pseudo-arclength continuation methods for both Rydberg-dressed Bose-Einstein condensates (BEC), and binary Rydberg-dressed BEC which are governed by the Gross-Pitaevskii equations (GPEs). A divide-and-conquer technique is proposed for rescaling the range/ranges of nonlocal nonlinear term/terms, which gives enough information for choosing a proper stepsize. This guarantees that the solution curve we wish to trace can be precisely approximated. In addition, the ground state solution would successfully evolve from one peak to vortices when the affect of the rotating term is imposed. Moreover, parameter variables with different number of components are exploited in curve-tracing. The proposed methods have the advantage of tracing the ground state solution curve once to compute the contours for various values of the coefficients of the nonlocal nonlinear term/terms. Our numerical results are consistent with those published in the literatures.

We study the numerical performance of a continuous data assimilation (downscaling) algorithm, based on ideas from feedback control theory, in the context of the two-dimensional incompressible Navier-Stokes equations. Our model problem is to recover an unknown reference solution, asymptotically in time, by using continuous-in-time coarse-mesh nodal-point observational measurements of the velocity field of this reference solution (subsampling), as might be measured by an array of weather vane anemometers. Our calculations show that the required nodal observation density is remarkably less than what is suggested by the analytical study; and is in fact comparable to the number of numerically determining Fourier modes, which was reported in an earlier computational study by the authors. Thus, this method is computationally efficient and performs far better than the analytical estimates suggest.

Neighbour search (NS) is the core of any implementations of smoothed particle hydrodynamics (SPH). In this paper,we present an efficient neighbour search method based on the plane sweep (PW) algorithm with N being the number of SPH particles. The resulting method, dubbed the PWNS method, is totally independent of grids (i.e., purely meshfree) and capable of treating variable smoothing length, arbitrary particle distribution and heterogenous kernels. Several state-of-the-art data structures and algorithms, e.g., the segment tree and the Morton code, are optimized and implemented. By simply allowingmultiple lines to sweep the SPH particles simultaneously from different initial positions, a parallelization of the PWNS method with satisfactory speedup and load-balancing can be easily achieved. That is, the PWNS SPH solver has a great potential for large scale fluid dynamics simulations.

In this paper, we propose an uniformly convergent adaptive finite element method with hybrid basis (AFEM-HB) for the discretization of singularly perturbed nonlinear eigenvalue problems under constraints with applications in Bose-Einstein condensation (BEC) and quantum chemistry. We begin with the time-independent Gross-Pitaevskii equation and show how to reformulate it into a singularly perturbed nonlinear eigenvalue problem under a constraint. Matched asymptotic approximations for the problem are reviewed to confirm the asymptotic behaviors of the solutions in the boundary/interior layer regions. By using the normalized gradient flow, we propose an adaptive finite element with hybrid basis to solve the singularly perturbed nonlinear eigenvalue problem. Our basis functions and the mesh are chosen adaptively to the small parameter ε. Extensive numerical results are reported to show the uniform convergence property of our method. We also apply the AFEM-HB to compute the ground and excited states of BEC with box/harmonic/optical lattice potential in the semiclassical regime (0 <ε≪C 1). In addition, we give a detailed error analysis of our AFEM-HB to a simpler singularly perturbed two point boundary value problem, show that our method has a minimum uniform convergence order

In this paper, we study a new stabilized method based on the local pressure projection to solve the semi-linear elliptic equation. The proposed scheme combines nonconforming finite element pairs NCP1–P1 triangle element and two-level method, which has a number of attractive computational properties: parameter-free, avoiding higher-order derivatives or edge-based data structures, but have more favorable stability and less support sets. Stability analysis and error estimates have been done. Finally, numerical experiments to check estimates are presented.

In this paper, the (G'/G)-expansion method is suggested to establish new exact solutions for fractional differential-difference equations in the sense of modified Riemann-Liouville derivative. The fractional complex transform is proposed to convert a fractional partial differential difference equation into its differential difference equation of integer order. With the aid of symbolic computation, we choose nonlinear lattice equations to illustrate the validity and advantages of the algorithm. It is shown that the proposed algorithm is effective and can be used for many other nonlinear lattice equations in mathematical physics and applied mathematics.

In this work, the Bishop and Love models for longitudinal vibrations are adopted to study the dynamics of isotropic rods with conical and exponential cross-sections. Exact solutions of both models are derived, using appropriate transformations. The analytical solutions of these two models are obtained in terms of generalised hypergeometric functions and Legendre spherical functions respectively. The exact solution of Love model for a rod with exponential cross-section is expressed as a sum of Gauss hypergeometric functions. The models are solved numerically by using the method of lines to reduce the original PDE to a system of ODEs. The accuracy of the numerical approximations is studied in the case of special solutions.

We investigate the decay for |x|→∞ of weak Sobolev-type solutions of semilinear nonlocal equations Pu = F(u). We consider the case when P = p(D) is an elliptic Fourier multiplier with polyhomogeneous symbol p(ξ), and we derive algebraic decay estimates in terms of weighted Sobolev norms. Our basic example is the celebrated Benjamin–Ono equation

for internal solitary waves of deep stratified fluids. Their profile presents algebraic decay, in strong contrast with the exponential decay for KdV shallow water waves.

We prove the linear and nonlinear instability of periodic traveling wave solutions for a generalized version of the symmetric regularized long wave (SRLW) equation. Using analytic and asymptotic perturbation theory, we establish sufficient conditions for the existence of exponentially growing solutions to the linearized problem and so the linear instability of periodic profiles is obtained. An application of this approach is made to obtain the linear/nonlinear instability of cnoidal wave solutions for the modified SRLW (mSRLW) equation. We also prove the stability of dnoidal wave solutions associated to the equation just mentioned.

We consider an elliptic self-adjoint first-order differential operator $L$ acting on pairs (2-columns) of complex-valued half-densities over a connected compact three-dimensional manifold without boundary. The principal symbol of the operator $L$ is assumed to be trace-free and the subprincipal symbol is assumed to be zero. Given a positive scalar weight function, we study the weighted eigenvalue problem for the operator $L$ . The corresponding counting function (number of eigenvalues between zero and a positive $\unicode[STIX]{x1D706}$ ) is known to admit, under appropriate assumptions on periodic trajectories, a two-term asymptotic expansion as $\unicode[STIX]{x1D706}\rightarrow +\infty$ and we have recently derived an explicit formula for the second asymptotic coefficient. The purpose of this paper is to establish the geometric meaning of the second asymptotic coefficient. To this end, we identify the geometric objects encoded within our eigenvalue problem—metric, non-vanishing spinor field and topological charge—and express our asymptotic coefficients in terms of these geometric objects. We prove that the second asymptotic coefficient of the counting function has the geometric meaning of the massless Dirac action.

Numerical oscillation of the total energy can be observed when the Kohn- Sham equation is solved by real-space methods to simulate the translational move of an electronic system. Effectively remove or reduce the unphysical oscillation is crucial not only for the optimization of the geometry of the electronic structure, but also for the study of molecular dynamics. In this paper, we study such unphysical oscillation based on the numerical framework in [G. Bao, G. H. Hu, and D. Liu, An h-adaptive finite element solver for the calculations of the electronic structures, Journal of Computational Physics, Volume 231, Issue 14, Pages 4967–4979, 2012], and deliver some numerical methods to constrain such unphysical effect for both pseudopotential and all-electron calculations, including a stabilized cubature strategy for Hamiltonian operator, and an a posteriori error estimator of the finite element methods for Kohn-Sham equation. The numerical results demonstrate the effectiveness of our method on restraining unphysical oscillation of the total energies.

We prove the quasiinvariance of Gaussian measures (supported by functions of increasing Sobolev regularity) under the flow of one-dimensional Hamiltonian partial differential equations such as the regularized long wave, also known as the Benjamin–Bona–Mahony (BBM) equation.

We prove global well-posedness of the time-dependent degenerate thermistor problem by establishing a uniform-in-time bounded mean ocsillation (BMO) estimate of inhomogeneous parabolic equations. Applying this estimate to the temperature equation, we derive a BMO bound of the temperature uniform with respect to time, which implies that the electric conductivity is an $A_{2}$ weight. The Hölder continuity of the electric potential is then proved by applying the De Giorgi–Nash–Moser estimate for degenerate elliptic equations with an $A_{2}$ coefficient. The uniqueness of the solution is proved based on the established regularity of the weak solution. Our results also imply the existence of a global classical solution when the initial and boundary data are smooth.