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.

In this paper, we present a discontinuous Galerkin method with staggered hybridization to discretize a class of nonlinear Stokes equations in two dimensions. The utilization of staggered hybridization is new and this approach combines the features of traditional hybridization method and staggered discontinuous Galerkin method. The main idea of our method is to use hybrid variables to impose the staggered continuity conditions instead of enforcing them in the approximation space. Therefore, our method enjoys some distinctive advantages, including mass conservation, optimal convergence and preservation of symmetry of the stress tensor. We will also show that, one can obtain superconvergent and strongly divergence-free velocity by applying a local postprocessing technique on the approximate solution. We will analyze the stability and derive a priori error estimates of the proposed scheme. The resulting nonlinear system is solved by using the Newton’s method, and some numerical results will be demonstrated to confirm the theoretical rates of convergence and superconvergence.