Bodhinanda (Nanda) Chandra

This paper presents a novel stabilized mixed material point method (MPM) designed for the unified modeling of free-surface and seepage flow. The unified formulation integrates the Navier-Stokes equation with the Darcy-Brinkman-Forchheimer equation, effectively capturing flows in both non-porous and porous domains. In contrast to the conventional Eulerian computational fluid dynamics (CFD) solver, which solves the velocity and pressure fields as unknown variables, the proposed method employs a monolithic displacement-pressure formulation adopted from the mixed-form updated-Lagrangian finite element method (FEM). To satisfy the discrete inf-sup stability condition, a stabilization strategy based on the variational multiscale method (VMS) is derived and integrated into the proposed formulation. Another distinctive feature is the implementation of blurred interfaces, which facilitate a seamless and stable transition of flows between free and porous domains, as well as across two distinct porous media. The efficacy of the proposed formulation is verified and validated through several benchmark cases in 1D, 2D, and 3D scenarios. Conducted numerical examples demonstrate enhanced accuracy and stability compared to analytical, experimental, and other numerical solutions.

This paper proposes novel and robust stabilization strategies for accurately modeling incom- pressible fluid flow problems in the material point method (MPM). To address the modeling of Newtonian fluids with incompressibility constraints, a new mixed implicit MPM formulation is proposed. Here, instead of solving the velocity and pressure fields as the unknown variables like the typical Eulerian computational fluid dynamics (CFD) solver, the proposed approach adopts a monolithic displacement–pressure formulation inspired by the mixed-form updated- Lagrangian Finite Element Method (FEM). To satisfy the inf–sup stability condition, two stabilization strategies are integrated into the formulation: the variational multiscale method (VMS) and the pressure-stabilization Petrov–Galerkin method (PSPG). By concurrently solving the displacement and pressure fields, the developed monolithic solver obviates the need for free-surface detection as well as Dirichlet and Neumann pressure imposition, in contrast to the fractional-step method. This attribute mitigates spurious pressure and velocity oscillations in simulating dynamic and transient flow problems. This study also addresses other prevalent challenges in MPM simulations, such as the pressure oscillations triggered by cell-crossing errors, particle-distribution-induced quadrature errors, and particle-grid information transfer. To resolve these issues, the quadratic B-Spline basis function, the delta-correction method, and the Taylor particle-in-cell method are incorporated into the proposed mixed MPM formulation, thereby enhancing numerical stability. The efficacy of the proposed stabilized incompressible MPM framework is validated through several benchmark cases, comparing the obtained results with other numerical methods and analytical solutions. Furthermore, the method’s capability in simulating real-world problems involving violent free-surface fluid motion is demonstrated through comparisons with experimental results of water sloshing and dam break scenarios.

An enhanced XMPM formulation is proposed to simulate the evolution of shear bands and post-failure behaviors with large deformations. In the current study, a localization search algorithm based on the theory of bifurcation is integrated into the XMPM to predict the initiation and propagation of discontinuity. In addition, a formulation of self-contact is proposed to deal with the dynamic frictional contact mechanism between the generated shear planes. In order to ensure the smoothness of the discontinuity surface during localization propagation, a hybrid implicit–explicit description of discontinuity is assumed by employing the level-set method and a point cloud approach. Several numerical examples are investigated to assess the accuracy and demonstrate the capability of the proposed XMPM approach in simulating the shear band evolution of different engineering problems in both 2D and 3D. The proposed formulation also exhibits minor sensitivity with respect to mesh refinements in predicting the shear-band path. To show the optimum performance compared to the regular MPM approach, a simulation of large-deformation heterogeneous slope failure is presented as one of the applications of the proposed method towards simulating real-scale engineering problems.

In many geomechanics applications, material boundaries are subjected to large displacements and deformation. Under these circumstances, the application of boundary conditions using particle methods, such as the material point method (MPM), becomes a challenging task since material boundaries do not coincide with the background mesh. This paper presents a formulation of penalty augmentation to impose nonhomogeneous, nonconforming Dirichlet boundary conditions in implicit MPM. The penalty augmentation is implemented utilizing boundary particles, which can move either according to or independently from the material deformation. Furthermore, releasing contact boundary condition, as well as the capability to accommodate slip boundaries, is introduced in the current work. The accuracy of the proposed method is assessed in both 2D and 3D cases, by convergence analysis reaching the analytical solution and by comparing the results of nonconforming and classical grid-conforming simulations. Keywords