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.

We present a moving-least-square immersed boundary method for solving viscous incompressible flow involving deformable and rigid boundaries on a uniform Cartesian grid. For rigid boundaries, no-slip conditions at the rigid interfaces are enforced using the immersed-boundary direct-forcing method. We propose a reconstruction approach that utilizes moving least squares (MLS) method to reconstruct the velocity at the forcing points in the vicinity of the rigid boundaries. For deformable boundaries, MLS method is employed to construct the interpolation and distribution operators for the immersed boundary points in the vicinity of the rigid boundaries instead of using discrete delta functions. The MLS approach allows us to avoid distributing the Lagrangian forces into the solid domains as well as to avoid using the velocity of points inside the solid domains to compute the velocity of the deformable boundaries. The present numerical technique has been validated by several examples including a Poiseuille flow in a tube, deformations of elastic capsules in shear flow and dynamics of red-blood cell in microfluidic devices.

An adaptive moving mesh finite difference method is presented to solve two types of equations with dynamic capillary pressure effect in porous media. One is the non-equilibrium Richards Equation and the other is the modified Buckley-Leverett equation. The governing equations are discretized with an adaptive moving mesh finite difference method in the space direction and an implicit-explicit method in the time direction. In order to obtain high quality meshes, an adaptive monitor function with directional control is applied to redistribute the mesh grid in every time step, then a diffusive mechanism is used to smooth the monitor function. The behaviors of the central difference flux, the standard local Lax-Friedrich flux and the local Lax-Friedrich flux with reconstruction are investigated by solving a 1D modified Buckley-Leverett equation. With the moving mesh technique, good mesh quality and high numerical accuracy are obtained. A collection of one-dimensional and two-dimensional numerical experiments is presented to demonstrate the accuracy and effectiveness of the proposed method.

A second order Ghost Fluid method is proposed for the treatment of interface problems of elliptic equations with discontinuous coefficients. By appropriate use of auxiliary virtual points, physical jump conditions are enforced at the interface. The signed distance function is used for the implicit description of irregular domain. With the additional unknowns, high order approximation considering the discontinuity can be built. To avoid the ill-conditioned matrix, the interpolation stencils are selected adaptively to balance the accuracy and the numerical stability. Additional equations containing the jump restrictions are assembled with the original discretized algebraic equations to form a new sparse linear system. Several Krylov iterative solvers are tested for the newly derived linear system. The results of a series of 1-D, 2-D tests show that the proposed method possesses second order accuracy in L^∞ norm. Besides, the method can be extended to the 3-D problems straightforwardly. Numerical results reveal the present method is highly efficient and robust in dealing with the interface problems of elliptic equations.

This paper presents a GPU-accelerated implementation of the Locally Optimal Block Preconditioned Conjugate Gradient (LOBPCG) method with an inexact null-space filtering approach to find eigenvalues in electromagnetics analysis with higher-order FEM. The performance of the proposed approach is verified using the Kepler (Tesla K40c) graphics accelerator, and is compared to the performance of the implementation based on functions from the Intel MKL on the Intel Xeon (E5-2680 v3, 12 threads) central processing unit (CPU) executed in parallel mode. Compared to the CPU reference implementation based on the Intel MKL functions, the proposed GPU-based LOBPCG method with inexact null-space filtering allowed us to achieve up to 2.9-fold acceleration.

The present paper follows our previous work [Yang et al., Phys. Rev. E, 90 (2014), 063011] in which the bending modes of a symmetric flexible fiber in viscous flows were studied by using a coupling approach of smoothed particle hydrodynamics (SPH) and element bending group (EBG). It was shown that a symmetric flexible fiber can undergo four different bending modes including stable U-shape, slight swing, violent flapping and stable closure modes. For an asymmetric flexible fiber, the bending modes can be different. This paper numerically studies the fiber shape, flow field and fluid drag of an asymmetric flexible fiber immersed in a viscous fluid flow by using the SPH-EBG coupling method. An asymmetric number is defined to describe the asymmetry of a flexible fiber. The effects of the asymmetric number on the fiber shape, flow field and fluid drag are investigated.

The fractional derivatives include nonlocal information and thus their calculation requires huge storage and computational cost for long time simulations. We present an efficient and high-order accurate numerical formula to speed up the evaluation of the Caputo fractional derivative based on the L2-1_σ formula proposed in [A. Alikhanov, J. Comput. Phys., 280 (2015), pp. 424-438], and employing the sum-of-exponentials approximation to the kernel function appeared in the Caputo fractional derivative. Both theoretically and numerically, we prove that while applied to solving time fractional diffusion equations, our scheme not only has unconditional stability and high accuracy but also reduces the storage and computational cost.

In this study, a numerical framework of the high order well-balanced weighted compact nonlinear (WCN) schemes is proposed for the shallow water equations based on the work in [S. Zhang, S. Jiang, C.-W Shu, J. Comput. Phys. 227 (2008) 7294-7321]. We employ a special splitting technique for the source term proposed in [Y. Xing, C.-W Shu, J. Comput. Phys. 208 (2005) 206-227] to maintain the exact C-property, which can be proved theoretically. In the meantime, the genuine high order accuracy of the numerical scheme can be observed successfully, and small perturbation of the stationary state can be resolved and evolved well. In order to capture the strong discontinuities and large gradients, the fifth-order upwind weighted nonlinear interpolations together with the fourth/sixth order cell-centered compact scheme are used to construct different WCN schemes. In addition, the local characteristic projections are considered to further restrain the potential numerical oscillations. A variety of representative one- and two-dimensional examples are tested to demonstrate the good performance of the proposed schemes.

In this paper we are concerned with numerical methods for nonlinear time-dependent problem coupled by electron, ion and photon temperatures in two dimensions, which is called the 2D-3T heat conduction equations. We propose discontinuous Galerkin (DG) methods for the discretization of the equations. For solving the resulting discrete system, we employ two domain decomposition (DD) preconditioners, one of which is associated with the non-overlapping DDM and the other is based on DDM with small overlap. The preconditioners are constructed by dropping the couplings between particles and each preconditioner consists of three preconditioners with smaller matrix size. To gauge the efficiency of the preconditioners, we test two examples and make different settings of parameters. Numerical results show that the proposed preconditioners are very effective to the 2D-3T problem.

Electrophoresis of a cylinder suspended in a cylindrical tube is analytically studied in the limit of thin electric double layer approximation. The electric and fluid flow fields within the annulus, and the cylinder velocities are analytically obtained in bipolar coordinates. The results are analyzed with various values of dimensionless parameters: eccentricity, cylinder-to-tube radius ratio and tube-to-cylinder zeta potential ratio (i.e., tube-to-cylinder velocity scale ratio). The analysis shows that microvortices are generated within the annulus. By changing the parameters, different flow patterns can be created, which shows potential for mixing enhancement in micro/nanofluidics. Moreover, the cylinder not only translates but also rotates when the cylinder and tube are eccentric. The cylinder rotation might be utilized as a micromotor or an electric field detector. The cylinder trajectories show that the cylinder may approach the tube wall or rest within the tube depending on the zeta potential ratio.

With the use of temporal derivative of flux function, a two-stage temporal discretization has been recently proposed in the design of fourth-order schemes based on the generalized Riemann problem (GRP) [21] and gas-kinetic scheme (GKS) [28]. In this paper, the fourth-order gas-kinetic scheme will be extended to solve the compressible multicomponent flow equations, where the two-stage temporal discretization and fifth-order WENO reconstruction will be used in the construction of the scheme. Based on the simplified two-species BGK model [41], the coupled Euler equations for individual species will be solved. Each component has its individual gas distribution function and the equilibrium states for each component are coupled by the physical requirements of total momentum and energy conservation in particle collisions. The second-order flux function is used to achieve the fourth-order temporal accuracy, and the robustness is as good as the second-order schemes. At the same time, the source terms, such as the gravitational force and the chemical reaction, will be explicitly included in the two-stage temporal discretization through their temporal derivatives. Many numerical tests from the shock-bubble interaction to ZND detonative waves are presented to validate the current approach.

In this work, we introduce an IMEX discontinuous Galerkin solver for the weakly compressible isentropic Euler equations. The splitting needed for the IMEX temporal integration is based on the recently introduced reference solution splitting [32, 52], which makes use of the incompressible solution. We show that the overall method is asymptotic preserving. Numerical results show the performance of the algorithm which is stable under a convective CFL condition and does not show any order degradation.