Moving mesh methods provide an efficient way of solving partial differential equations for which large, localised variations in the solution necessitate locally dense spatial meshes. In one-dimension, meshes are typically specified using the arclength mesh density function. This choice is well-justified for piecewise polynomial interpolants, but it is only justified for spectral methods when model solutions include localised steep gradients. In this paper, one-dimensional mesh density functions are presented which are based on a spatially localised measure of the bandwidth of the approximated model solution. In considering bandwidth, these mesh density functions are well-justified for spectral methods, but are not strictly tied to the error properties of any particular spatial interpolant, and are hence widely applicable. The bandwidth mesh density functions are illustrated in two ways. First, by applying them to Chebyshev polynomial approximation of two test functions, and second, through use in periodic spectral and finite-difference moving mesh methods applied to a number of model problems in acoustics. These problems include a heterogeneous advection equation, the viscous Burgers’ equation, and the Korteweg-de Vries equation. Simulation results demonstrate solution convergence rates that are up to an order of magnitude faster using the bandwidth mesh density functions than uniform meshes, and around three times faster than those using the arclength mesh density function.

Numerical solutions of a non-Fickian diffusion equation belonging to a hyperbolic type are presented in one space dimension. The Brownian particle modelled by this diffusion equation is subjected to a symmetric periodic potential whose spatial shape can be varied by a single parameter. We consider a numerical method which consists of applying Laplace transform in time; we then obtain an elliptic diffusion equation which is discretized using a finite difference method. We analyze some aspects of the convergence of the method. Numerical results for particle density, flux and mean-square-displacement (covering both inertial and diffusive regimes) are presented.

In this article, we present two new novel finite difference approximations of order two and four, respectively, for the three dimensional non-linear triharmonic partial differential equations on a compact stencil where the values of u, ∂2u/∂n2 and ∂4u/∂n4 are prescribed on the boundary. We introduce new ideas to handle the boundary conditions and there is no need to discretize the derivative boundary conditions. We require only 7- and 19-grid points on the compact cell for the second and fourth order approximation, respectively. The Laplacian and the biharmonic of the solution are obtained as by-product of the methods. We require only system of three equations to obtain the solution. Numerical results are provided to illustrate the usefulness of the proposed methods.

This paper establishes relations between the stability and the high-order truncated corrections for modeling of the mass conservation equation with the two-relaxation-times (TRT) collision operator. First we propose a simple method to derive the truncation errors from the exact, central-difference type, recurrence equations of the TRT scheme. They also supply its equivalent three-time-level discretization form. Two different relationships of the two relaxation rates nullify the third (advection) and fourth (pure diffusion) truncation errors, for any linear equilibrium and any velocity set. However, the two relaxation times alone cannot remove the leading-order advection-d if fusion error, because of the intrinsic fourth-order numerical diffusion. The truncation analysis is carefully verified for the evolution of concentration waves with the anisotropic diffusion tensors. The anisotropic equilibrium functions are presented in a simple but general form, suitable for the minimal velocity sets and the d2Q9, d3Q13, d3Q15 and d3Q19 velocity sets. All anisotropic schemes are complemented by their exact necessary von Neumann stability conditions and equivalent finite-difference stencils. The sufficient stability conditions are proposed for the most stable (OTRT) family, which enables modeling at any Peclet numbers with the same velocity amplitude. The heuristic stability analysis of the fourth-order truncated corrections extends the optimal stability to larger relationships of the two relaxation rates, in agreement with the exact (one-dimensional) and numerical (multi-dimensional) stability analysis. A special attention is put on the choice of the equilibrium weights. By combining accuracy and stability predictions, several strategies for selecting the relaxation and free-tunable equilibrium parameters are suggested and applied to the evolution of the Gaussian hill.

This paper discusses a new fourth-order compact off-step discretization for the solution of a system of two-dimensional nonlinear elliptic partial differential equations subject to Dirichlet boundary conditions. New methods to obtain the fourth-order accurate numerical solution of the first order normal derivatives of the solution are also derived. In all cases, we use only nine grid points to compute the solution. The proposed methods are directly applicable to singular problems and problems in polar coordinates, which is a main attraction. The convergence analysis of the derived method is discussed in detail. Several physical problems are solved to demonstrate the usefulness of the proposed methods.

We investigate the ultra weak variational formulation (UWVF) of the 2-D Helmholtz equation using a new choice of basis functions. Traditionally the UWVF basis functions are chosen to be plane waves. Here, we instead use first kind Bessel functions. We compare the performance of the two bases. Moreover, we show that it is possible to use coupled plane wave and Bessel bases in the same mesh. As test cases we shall consider propagating plane and evanescent waves in a rectangular domain and a singular 2-D Helmholtz problem in an L-shaped domain.

The emission of scission neutrons from fissioning nuclei is of high practical interest. To study this process we have used the sudden approximation and also a more realistic approach that takes into account the scission dynamics. Numerically, this implies the solution of the bi-dimensional Schrödinger equation, both stationary and time-dependent. To describe axially symmetric extremely deformed nuclear shapes, we have used the Cassini parametrization. The Hamiltonian is discretized by using finite difference approximations of the derivatives. The main computational challenges are the solution of algebraic eigenvalue problems and of linear systems with large sparse matrices. We have employed appropriate procedures (Arnoldi and bi-conjugate gradients). The numerical solutions have been used to evaluate physical quantities, like the number of emitted neutrons per scission event, the primary fragments’ excitation energy and the distribution of the emission points.