In this paper, we apply high-order finite difference (FD) schemes for multispecies and multireaction detonations (MMD). In MMD, the density and pressure are positive and the mass fraction of the ith species in the chemical reaction, say zi, is between 0 and 1, with ∑z_i=1. Due to the lack of maximum-principle, most of the previous bound-preserving technique cannot be applied directly. To preserve those bounds, we will use the positivity-preserving technique to all the z′is and enforce ∑z_i=1 by constructing conservative schemes, thanks to conservative time integrations and consistent numerical fluxes in the system. Moreover, detonation is an extreme singular mode of flame propagation in premixed gas, and the model contains a significant stiff source. It is well known that for hyperbolic equations with stiff source, the transition points in the numerical approximations near the shocks may trigger spurious shock speed, leading to wrong shock position. Intuitively, the high-order weighted essentially non-oscillatory (WENO) scheme, which can suppress oscillations near the discontinuities, would be a good choice for spatial discretization. However, with the nonlinear weights, the numerical fluxes are no longer “consistent”, leading to nonconservative numerical schemes and the bound-preserving technique does not work. Numerical experiments demonstrate that, without further numerical techniques such as subcell resolutions, the conservative FD method with linear weights can yield better numerical approximations than the nonconservative WENO scheme.

In this paper, we study the classical Allen-Cahn equations and investigate the maximum-principle-preserving (MPP) techniques. The Allen-Cahn equation has been widely used in mathematical models for problems in materials science and fluid dynamics. It enjoys the energy stability and the maximum-principle. Moreover, it is well known that the Allen-Cahn equation may yield thin interface layer, and nonuniform meshes might be useful in the numerical solutions. Therefore, we apply the local discontinuous Galerkin (LDG) method due to its flexibility on h-p adaptivity and complex geometry. However, the MPP LDG methods require slope limiters, then the energy stability may not be easy to obtain. In this paper, we only discuss the MPP technique and use numerical experiments to demonstrate the energy decay property. Moreover, due to the stiff source given in the equation, we use the conservative modified exponential Runge-Kutta methods and thus can use relatively large time step sizes. Thanks to the conservative time integration, the bounds of the unknown function will not decay. Numerical experiments will be given to demonstrate the good performance of the MPP LDG scheme.

Redirected walking (RDW) aims to reduce the collisions in the physical space for VR applications. However, most of the previous RDW methods do not consider future possibilities of collisions after imperceptibly redirecting users. In this paper, we combine the subtle RDW methods and reset strategy in our method design and propose a novel solution for RDW that can make better use of physical space and trigger fewer resets. The key idea of our method is to discretize the representation of possible user positions and orientations by a series of standard poses and rate them based on the possibilities of hitting obstacles of their reachable poses. A transfer path algorithm is proposed to measure the accessibility among standard poses and is used to support the calculation of the scores of standard poses. Using our method, the user can be redirected imperceptibly to the optimal pose with the best score among all the reachable poses from the user’s current pose during walking. Experiments demonstrate that our method outperforms state-of-the-art methods in various environment sizes and obstacle layouts.

This paper proposes a second-order pedestrian model that comprises two types of equations: continuity equation and a set of transport equations. To complete the model, we develop Eikonal equations to explicitly consider the effects of the collective decisions of individuals and crowd pressure on pedestrian dynamics. Then, the crowd movement is simulated using a set of partial differential equations under appropriate initial and boundary conditions. Based on the stability requirements derived by performing a standard linear stability analysis, suitable parameters are selected to test the model in a numerical example. The proposed second-order system is then solved using the characteristic-wise third-order weighted essentially non-oscillatory (WENO3) scheme, and the Eikonal equations are solved using the fast sweeping method. The numerical results indicate the effectiveness of the model because the derived local flow-density relationship produces a second peak in the high-density region, which is consistent with previous empirical studies. Besides, the applicability of the model to an unstable condition is verified through the simulation of complex phenomena such as stop-and-go waves. Furthermore, the estimate of crowd pressure in the simulation results can be used as a risk-level indicator for crowd management and control.

The simulation of Rayleigh waves is important in a variety of geophysical applications. The computational challenge is the fact that very fine mesh is necessary as the waves are concentrated at the free surface and decay exponentially away from the free surface. To overcome this challenge and to develop a robust high order scheme for the simulation of Rayleigh waves, we develop a mortar discontinuous Galerkin method with staggered hybridization. The use of the mortar technique allows one to use fine mesh in only a local region near the free surface, and use coarse mesh in most of the domain. This approach reduces the computational cost significantly. The staggered hybridization allows the preservation of the strong symmetry of the stress tensor without complicated construction of basis functions. In particular, the basis functions are piecewise polynomial without any continuity requirement, and the coupling of the basis functions is performed by using carefully chosen hybridized variables. The resulting scheme is explicit in time, and only local saddle point system is solved for each time step. We will present several benchmark problems to demonstrate the performance of the proposed method.