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.

This paper is about the construction of a unified gas-kinetic scheme (UGKS) for a coupled system of radiative transport and material heat conduction with different diffusive limits. Different from the previous approach, instead of including absorption/emission only, the current method takes both scattering and absorption/emission mechanism into account in the radiative transport process. As a result, two asymptotic limiting solutions will appear in the diffusive regime. In the strong absorption/emission case, an equilibrium diffusion limit is obtained, where the system is mainly driven by a nonlinear diffusion equation for the equilibrium radiation and material temperature. However, in the strong scattering case, a non-equilibrium limit can be obtained, where coupled nonlinear diffusion system with different radiation and material temperature is obtained. In addition to including the scattering term in the transport equation, an implicit UGKS (IUGKS) will be developed in this paper as well. In the IUGKS, the numerical flux for the radiation intensity is constructed implicitly. Therefore, the conventional CFL constraint for the time step is released. With the use of a large time step for the radiative transport, it becomes possible to couple the IUGKS with the gas dynamic equations to develop an efficient numerical method for radiative hydrodynamics. The IUGKS is a valid method for all radiative transfer regimes. A few numerical examples will be presented to validate the current implicit method for both optical thin to optical thick cases.

The computation of the radiative transfer equation is expensive mainly due to two stiff terms: the transport term and the collision operator. The stiffness in the former comes from the fact that particles (such as photons) travel at the speed of light, while that in the latter is due to the strong scattering in the optically thick region. We study the fully implicit scheme for this equation to account for the stiffness. The main challenge in the implicit treatment is the coupling between the spacial and angular coordinates that requires the large size of the to-be-inverted matrix, which is also ill-conditioned and not necessarily symmetric. Our main idea is to utilize the spectral structure of the ill-conditioned matrix to construct a pre-conditioner, which, along with an exquisite split of the spatial and angular dependence, significantly improve the condition number and allows a matrix-free treatment. We also design a fast solver to compute this pre-conditioner explicitly in advance. Our method is shown to be efficient in both diffusive and free streaming limit, and the computational cost is comparable to the state-of-the-art method. Various examples including anisotropic scattering and two-dimensional problems are provided to validate the effectiveness of our method.

In this paper, an alternative lattice Boltzmann (LB)model for incompressible flows is proposed. By modifying directly the moments of the equilibrium distribution function (EDF), the continuous expression of the EDF in tensor Hermite polynomials is derived using the moment expansion and then discretizedwith the discrete velocity vectors of the D2Q9 lattice. The present model as well as its counterpart, the incompressible LB model proposed by Guo, reproduces the incompressible Navier-Stokes (N-S) equations for both steady and unsteady flows. Besides, an alternative pressure formula, which represents the pressure as the diagonal part of the stress tensor, is adopted in the present model. Furthermore, in order to enhance the stability of the present LB model, an additional relaxation time pertaining to the non-hydrodynamic mode is added to the BGK collision operator. The present LB model is validated by two benchmark tests: the cavity flow with different Reynolds number (Re) and the flow past an impulsively started cylinder at Re=40 and 550.

We employ the lattice Boltzmann method and random walk particle tracking to simulate the time evolution of hydrodynamic dispersion in bulk, random, monodisperse, hard-sphere packings with bed porosities (interparticle void volume fractions) between the random-close and the random-loose packing limit. Using Jodrey-Tory and Monte Carlo-based algorithms and a systematic variation of the packing protocols we generate a portfolio of packings, whose microstructures differ in their degree of heterogeneity (DoH). Because the DoH quantifies the heterogeneity of the void space distribution in a packing, the asymptotic longitudinal dispersion coefficient calculated for the packings increases with the packings’ DoH. We investigate the influence of packing length (up to 150 dp, where dp is the sphere diameter) and grid resolution (up to 90 nodes per dp) on the simulated hydrodynamic dispersion coefficient, and demonstrate that the chosen packing dimensions of 10 dpx 10 dpx 70 dp and the employed grid resolution of 60 nodes per dp are sufficient to observe asymptotic behavior of the dispersion coefficient and to minimize finite size effects. Asymptotic values of the dispersion coefficients calculated for the generated packings are compared with simulated as well as experimental data from the literature and yield good to excellent agreement.

In this paper, a lattice Boltzmann BGK (LBGK) model is proposed for simulating incompressible axisymmetric flows. Unlike other existing axisymmetric lattice Boltzmann models, the present LBGK model can eliminate the compressible effects only with the small Mach number limit. Furthermore the source terms of the model are simple and contain no velocity gradients. Through the Chapman-Enskog expansion, the macroscopic equations for incompressible axisymmetric flows can be exactly recovered from the present LBGK model. Numerical simulations of the Hagen-Poiseuille flow, the pulsatile Womersley flow, the flow over a sphere, and the swirling flow in a closed cylindrical cavity are performed. The results agree well with the analytic solutions and the existing numerical or experimental data reported in some previous studies.

A cluster dynamics model based on rate theory has been developed to describe the accumulation and diffusion processes of helium in tungsten under helium implantation alone or synergistic irradiation with neutron, by involving different types of objects, adopting up-to-date parameters and complex reaction processes as well as considering the diffusion process along with depth. The calculated results under different conditions are in good agreement with experiments much well. The model describes the behavior of helium in tungsten within 2D space of defect type/size and depth on different ions incident conditions (energies and fluences) and material conditions (system temperature and existent sinks), by including the synergistic effect of helium-neutron irradiations and the influence of inherent sinks (dislocation lines and grain boundaries). The model, coded as IRadMat, would be universally applicable to the evolution of defects for ions/neutron irradiated on plasma-facing materials.

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.

We present a new Finite Volume Evolution Galerkin (FVEG) scheme for the solution of the shallow water equations (SWE) with the bottom topography as a source term. Our new scheme will be based on the FVEG methods presented in (Noelle and Kraft, J. Comp. Phys., 221 (2007)), but adds the possibility to handle dry boundaries. The most important aspect is to preserve the positivity of the water height. We present a general approach to ensure this for arbitrary finite volume schemes. The main idea is to limit the outgoing fluxes of a cell whenever they would create negative water height. Physically, this corresponds to the absence of fluxes in the presence of vacuum. Wellbalancing is then re-established by splitting gravitational and gravity driven parts of the flux. Moreover, a new entropy fix is introduced that improves the reproduction of sonic rarefaction waves.

A computational study on the enhancement of the second harmonic generation (SHG) in one-dimensional (1D) photonic crystals is presented. The mathematical model is derived from a nonlinear system of Maxwell’s equations, which partly overcomes the shortcoming of some existing models based on the undepleted pump approximation. We designed an iterative scheme coupled with the finite element method which can be applied to simulate the SHG in one dimensional nonlinear photonic band gap structures in our previous work. For the case that the nonlinearity is strong which is desirable to enhance the conversion efficiency, a continuation method is introduced to ensure the convergence of the iterative procedure. The convergence of our method is fast. Numerical experiments also indicate the conversion efficiency of SHG can be significantly enhanced when the frequencies of the fundamental and the second harmonic wave are tuned at the photonic band edges. The maximum total conversion efficiency available reaches more than 50% in all the cases studied.