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 will develop a new staggered hybridization technique for discontinuous Galerkin methods to discretize linear elastodynamic equations. The idea of hybridization is used extensively in many discontinuous Galerkin methods, but the idea of staggered hybridization is new. Our new approach offers several advantages, namely energy conservation, high-order optimal convergence, preservation of symmetry for the stress tensor, block diagonal mass matrices as well as low dispersion error. The key idea is to use two staggered hybrid variables to enforce the continuity of the velocity and the continuity of the normal component of the stress tensor on a staggered mesh. We prove the stability and the convergence of the proposed scheme in both the semi-discrete and the fully-discrete settings. Numerical results confirm the optimal rate of convergence and show that the method has a superconvergent property for dispersion. Furthermore, an application of this method to Rayleigh wave propagation is presented.

This paper proposes a two-level continuum transportation system approach to modeling a dynamic taxi traffic assignment (DTTA) problem in a dense network with real-time traffic information provision and three types of vehicles, including private cars, occupied taxis, and vacant taxis. The proposed approach treats the dense network as a continuum in the first level, in which private car and occupied taxi drivers are free to choose their paths in a two-dimensional continuous space. The proposed approach also divides the modeling region into many identical squares to form a cell-based network in the second level, in which the cells are classified into two categories: target cell with an acceptable expected rate of return (EROR) to vacant taxi drivers and non-target cell with an unacceptable EROR. The EROR associated with a cell is the ratio of the cumulative expected profit of a taxi driver who successfully picks up a customer during the customer search that starts from that cell to the sum of expected search time for this customer and expected occupied travel time to serve this customer. Based on the cell-based network, we develop a cell-based intervening opportunity model to capture the fact that vacant taxi drivers can meet a customer on the way to their destination zones and estimate the EROR. Each vacant taxi driver has a mixed strategy to determine his/her customer-search direction according to the EROR: Each vacant taxi driver in a target cell selects its neighbor cells with maximum EROR, and each vacant taxi driver in a non-target cell selects the travel time-based shortest path to his/her target cell. Meanwhile, each private car driver chooses the path that minimizes his/her own generalized travel cost, and each occupied taxi driver chooses the path that minimizes his/her customer's generalized in-vehicle travel cost. In our model, traffic density in the system is governed by the conservation law (CL), and the flow directions of different vehicles are determined by the path-choice strategies of their drivers, which are captured by Hamilton–Jacobi (HJ) equations. Both the proposed CL and HJ equations can be solved by the Lax–Friedrichs scheme, which forms the backbone of the developed solution algorithm. Finally, numerical examples and a case study are used to demonstrate the properties of the model, the performance of the solution algorithm, and the value of using our methodology for estimating network performance.

We present a stochastic Justh–Krishnaprasad flocking model describing interactions among individuals in a planar domain with their positions and heading angles. The deterministic counterpart of the proposed model describes the formation of nematic alignment in an ensemble of planar particles moving with a unit speed. When the noise is turned off, we show that the nematic alignment state, in which all heading angles are either same or the opposite, is nonlinearly stable using a Lyapunov functional approach. We employed a diameter-like functional via the rearrangement of heading angles in the 2π-interval. In contrast, under the additive noise, a continuous angle configuration will be deviated asymptotically from the nematic state. Nevertheless, in any finite-time interval, we will see that some part of angle configuration will stay close to the nematic state with a positive probability, where we call this phenomenon as stochastic persistency. We provide a quantitative estimate on the probability for stochastic persistency and compare several numerical examples with analytical results.