In this work, we concern with the numerical approach for delay differential equations with random coefficients. We first show that the exact solution of the problem considered admits good regularity in the random space, provided that the given data satisfy some reasonable assumptions. A stochastic collocation method is proposed to approximate the solution in the random space, and we use the Legendre spectral collocation method to solve the resulting deterministic delay differential equations. Convergence property of the proposed method is analyzed. It is shown that the numerical method yields the familiar exponential order of convergence in both the random space and the time space. Numerical examples are given to illustrate the theoretical results.

In this paper, we study a high-order compact difference scheme for the fourth-order fractional subdiffusion system. We consider the situation in which the unknown function and its first-order derivative are given at the boundary. The scheme is shown to have high order convergence. Numerical examples are given to verify the theoretical results.

A lattice Boltzmann flux solver (LBFS) is presented in this work for simulation of incompressible viscous and inviscid flows. The new solver is based on Chapman-Enskog expansion analysis, which is the bridge to link Navier-Stokes (N-S) equations and lattice Boltzmann equation (LBE). The macroscopic differential equations are discretized by the finite volume method, where the flux at the cell interface is evaluated by local reconstruction of lattice Boltzmann solution from macroscopic flow variables at cell centers. The new solver removes the drawbacks of conventional lattice Boltzmann method such as limitation to uniform mesh, tie-up of mesh spacing and time interval, limitation to viscous flows. LBFS is validated by its application to simulate the viscous decaying vortex flow, the driven cavity flow, the viscous flow past a circular cylinder, and the inviscid flow past a circular cylinder. The obtained numerical results compare very well with available data in the literature, which show that LBFS has the second order of accuracy in space, and can be well applied to viscous and inviscid flow problems with non-uniform mesh and curved boundary.

Interaction of vortex rings with solid is an important research topic of hydrodynamic. In this study, a multiple-relaxation time (MRT) lattice Boltzmann method (LBM) is used to investigate the flow of a vortex ring impacting spheroidal particles. The MRT-LBM is validated through the cases of vortex ring impacting a flat wall. The vortex evolution due to particle size, the aspect ratio of a prolate particle, as well as Reynolds $(Re)$ number are discussed in detail. When the vortex ring impacting a stationary sphere, the primary and secondary vortex rings wrap around each other, which is different from the situation of the vortex ring impacting a plate. For the vortex ring impacting with a prolate spheroid, the secondary vortex ring stretches mainly along the long axis of the ellipsoid particle. However, it is found that after the vortex wrapping stage, the primary vortex recovers along the short axis of the particle faster than that in the long axis, i.e., the primary vortex ring stretches mainly along the short axis of the particle. That has never been addressed in the literature.

Fluid-structure-interaction problems are ubiquitous, complicated, and not yet well understood. In this paper we investigate the interaction of a leading rigid circular cylinder and a trailing compliant filament and analyze the dynamic responses of the filament in the wake of the cylinder. It is revealed that there exist two flapping states of the filament depending on the cylinder-filament separation distance and the relevant critical distance distinguishing the two states is associated with the Reynolds number and the filament length. It is also found that the drag coefficient of the cylinder is reduced but that of the filament may be increased or decreased depending on its length. Compared with a single filament in a uniform flow, the filament of the same mechanical properties flapping in the wake of the cylinder has a lower frequency and a greater amplitude.

In this paper, we propose two new explicit multi-symplectic splitting methods for the nonlinear Dirac (NLD) equation. Based on its multi-symplectic formulation, the NLD equation is split into one linear multi-symplectic system and one nonlinear infinite Hamiltonian system. Then multi-symplectic Fourier pseudospectral method and multi-symplectic Preissmann scheme are employed to discretize  the linear subproblem, respectively. And the nonlinear subsystem is solved by a symplectic scheme. Finally, a composition method is applied to obtain the final schemes for the NLD equation. We find that the two proposed schemes preserve the total symplecticity and can be solved explicitly. Numerical experiments are presented to show the effectiveness of the proposed methods.

We consider an inverse problem of determining unknown coefficients for a one-dimensional analogue of radiative transport equation. We show that some combination of the unknown coefficients can be uniquely determined by giving pulse-like inputs at the boundary and observing the corresponding outputs. Our result can be applied for determination of absorption and scattering properties of an optically turbid medium if the radiative transport equation is appropriate for describing the propagation of light in the medium.

It is very important to predict the bypass transition in the simulation of flows through turbomachinery. This paper presents a four-equation eddy-viscosity turbulence transition model for prediction of bypass transition. It is based on the SST turbulence model and the laminar kinetic energy concept. A transport equation for the non-turbulent viscosity is proposed to predict the development of the laminar kinetic energy in the pre-transitional boundary layer flow which has been observed in experiments. The turbulence breakdown process is then captured with an intermittency transport equation in the transitional region. The performance of this new transition model is validated through the experimental cases of T3AM, T3A and T3B. Results in this paper show that the new transition model can reach good agreement in predicting bypass transition, and is compatible with modern CFD software by using local variables.