The self-similar growth-fragmentation equation describes the evolution of a medium in which particles grow and divide as time proceeds, with the growth and splitting of each particle depending only upon its size. The critical case of the equation, in which the growth and division rates balance one another, was considered in Doumic and Escobedo (2015) for the homogeneous case where the rates do not depend on the particle size. Here, we study the general self-similar case, using a probabilistic approach based on Lévy processes and positive self-similar Markov processes which also permits us to analyse quite general splitting rates. Whereas existence and uniqueness of the solution are rather easy to establish in the homogeneous case, the equation in the nonhomogeneous case has some surprising features. In particular, using the fact that certain self-similar Markov processes can enter (0,∞) continuously from either 0 or ∞, we exhibit unexpected spontaneous generation of mass in the solutions.

We consider the zero-resistivity limit for Hasegawa–Wakatani equations in a cylindrical domain when the initial data are Stepanov almost-periodic in the axial direction. First, we prove the existence of a solution to Hasegawa–Wakatani equations with zero resistivity; second, we obtain uniform a priori estimates with respect to resistivity. Such estimates can be obtained in the same way as for our previous results; therefore, the most important contribution of this paper is the proof of the existence of a local-in-time solution to Hasegawa–Wakatani equations with zero resistivity. We apply the theory of Bohr–Fourier series of Stepanov almost-periodic functions to such a proof.

We construct and analyze conservative local discontinuous Galerkin (LDG) methods for the Generalized Korteweg-de-Vries equation. LDG methods are designed by writing the equation as a system and performing separate approximations to the spatial derivatives. The main focus is on the development of conservative methods which can preserve discrete versions of the first two invariants of the continuous solution, and a posteriori error estimates for a fully discrete approximation that is based on the idea of dispersive reconstruction. Numerical experiments are provided to verify the theoretical estimates.

Integrable couplings of the Boiti-Pempinelli-Tu hierarchy are constructed by a class of non-semisimple block matrix loop algebras. Further, through using the variational identity theory, the Hamiltonian structures of those integrable couplings are obtained. The method can be applied to obtain other integrable hierarchies.

It is known that large time-stepping method are useful for simulating phase field models. In this work, an adaptive time-stepping strategy is proposed based on numerical energy stability and equi-distribution principle. The main idea is to use the energy variation as an indicator to update the time step, so that the resulting algorithm is free of user-defined parameters, which is different from several existing approaches. Some numerical experiments are presented to illustrate the effectiveness of the algorithms.

In the first part, we study the convergence of discrete solutions to splitting schemes for two-phase flow with different mass densities suggested in [Guillen-Gonzalez, Tierra, J.Comput.Math. (6)2014]. They have been formulated for the diffuse interface model in [Abels, Garcke, Grün, M3AS, 2012, DOI:10.1142/S0218202511500138] which is consistent with thermodynamics. Our technique covers various discretization methods for phase-field energies, ranging from convex-concave splitting to difference quotient approaches for the double-well potential. In the second part of the paper, numerical experiments are presented in two space dimensions to identify discretizations of Cahn-Hilliard energies which are ϕ-stable and which do not reduce the acceleration of falling droplets. Finally, 3d simulations in axial symmetric geometries are shown to underline even more the full practicality of the approach.

Consider the time-harmonic acoustic scattering from an extended elastic body surrounded by a finite number of point-like obstacles in a fluid. We assume point source waves are emitted from arrayed transducers and the signals of scattered near-field data are recorded by receivers not far away from the scatterers (compared to the incident wavelength). The forward scattering can be modeled as an interaction problem between acoustic and elastic waves together with a multiple scattering problem between the extend solid and point scatterers. We prove a necessary and sufficient condition that can be used simultaneously to recover the shape of the extended elastic solid and to locate the positions of point scatterers. The essential ingredient in our analysis is the outgoing-to-incoming (OtI) operator applied to the resulting near-field response matrix (or operator). In the first part, we justify the MUSIC algorithm for locating point scatterers from near-field measurements. In the second part, we apply the factorization method, the continuous analogue of MUSIC, to the two-scale scattering problem for determining both extended and point scatterers. Numerical examples in 2D are demonstrated to show the validity and accuracy of our inversion algorithms.

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.