Abstract We present an adaptively updated Lagrangian Material Point Method (A‐ULMPM) to alleviate non‐physical artifacts, such as the cell‐crossing instability and numerical fracture, that plague state‐of‐the‐art Eulerian formulations of MPM, while still allowing for large deformations that arise in fluid simulations. A‐ULMPM spans MPM discretizations from total Lagrangian formulations to Eulerian formulations. We design an easy‐to‐implement physics‐based criterion that allows A‐ULMPM to update the reference configuration adaptively for measuring physical states, including stress, strain, interpolation kernels and their derivatives. For better efficiency and conservation of angular momentum, we further integrate the APIC [JSS*15] and MLS‐MPM [HFG*18] formulations in A‐ULMPM by augmenting the accuracy of velocity rasterization using both the local velocity and its first‐order derivatives. Our theoretical derivations use a nodal discretized Lagrangian, instead of the weak form discretization in MLS‐MPM [HFG*!!18], and naturally lead to a “modified” MLS‐MPM in A‐ULMPM, which can recover MLS‐MPM using a completely Eulerian formulation. A‐ULMPM does not require significant changes to traditional Eulerian formulations of MPM, and is computationally more efficient since it only updates interpolation kernels and their derivatives during large topology changes. We present end‐to‐end 3D simulations of stretching and twisting hyperelastic solids, viscous flows, splashing liquids, and multi‐material interactions with large deformations to demonstrate the efficacy of our new method.

We propose an agent-to-agent contagion-immunity formulation that can simulate detailed COVID-19 spreading within moving crowds. Specifically, we develop a diffusion-based disease contagion model for discrete systems that considers the effect of health interventions, such as social distancing, immunity, and vaccination. We integrate our contagion-immunity formulation with the governing equations of motion for crowd dynamics for investigating the distribution of disease in crowds with different numbers of people. For the same crowd simulation, our model can interactively simulate virus spread for different initial distributions of infected people. To the best of our knowledge, our work is the first to simulate the disease contagion within moving crowds in computer graphics. Our numerical results for the number of infected people in unprotected dense crowds agree with the SIS model, while our model provides richer information for disease spread and shows that vaccination is the best health intervention to prevent infection.

Abstract We propose a Lagrangian particle‐based formulation for simulating deformation, fracture, and diffusion in thin membranelike structures, such as aluminium foil, rubbery films, and seaweed flakes. We integrate our model with diffusion processes and derive a unified framework for simulating deformation‐diffusion coupled phenomena, which is applied to provide realistic heterogeneity induced by the diffusion process to fracture patterns. To the best of our knowledge, our work is the first to simulate the complex fracture patterns of single‐layered membranes in computer graphics and introduce heterogeneity induced by the diffusion process, which generates more geometrically rich fracture patterns. Our end‐to‐end 3D simulations show that our deformation‐diffusion coupling framework captures detailed fracture growth patterns in thin membranes due to both in‐plane and out‐of‐plane motions, producing realistically wrinkled slit edges, and heterogeneity introduced due to diffusion.

We assume that the viscous forces in any liquid are simultaneously local and non-local , and introduce the extended POM-POM model [McLeish and Larson 1998; Oishi et al. 2012; Verbeeten et al. 2001] to computer graphics to design a unified constitutive model for viscosity that generalizes prior models, such as Oldroyd-B, the Upper-convected Maxwell (UCM) model [Sadeghy et al. 2005], and classical Newtonian viscosity under one umbrella, recovering each of them with different parameter values. Implicit discretization of our model via backward Euler recovers the variational Stokes solver of [Larionov et al. 2017] for Newtonian viscosity. For greater accuracy, however, we introduce the second-order accurate Generalized Single Step Single Solve (GS4) scheme [Tamma et al. 2000; Zhou and Tamma 2004] to computer graphics, which recovers all prior second-order accurate time integration schemes to date. Using GS4 and our generalized constitutive model, we present a Material Point Method (MPM) for simulating various viscoelastic liquid behaviors, such as classical liquid rope coiling, buckling, folding, and shear thinning/thickening. In addition, we show how to couple our viscoelastic liquid simulator with the recently introduced non-Fourier heat diffusion solver [Xue et al. 2020] for simulating problems with phase change, such as melting chocolate and digital fabrication with 3D printing. While the discretization of heat diffusion is slightly different within GS4, we show that it can still be efficiently solved using an assembly-free Multigrid-preconditioned Conjugate Gradients solver. We present end-to-end 3D simulations to demonstrate the versatility of our framework.

We introduce the C-F diffusion model [Anderson and Tamma 2006; Xue et al. 2018] to computer graphics for diffusion-driven problems that has several attractive properties: (a) it fundamentally explains diffusion from the perspective of the non-equilibrium statistical mechanical Boltzmann Transport Equation, (b) it allows for a finite propagation speed for diffusion, in contrast to the widely employed Fick's/Fourier's law, and (c) it can capture some of the most characteristic visual aspects of diffusion-driven physics, such as hydrogel swelling, limited diffusive domain for smoke flow, snowflake and dendrite formation, that span from Fourier-type to non-Fourier-type diffusive phenomena. We propose a unified convection-diffusion formulation using this model that treats both the diffusive quantity and its associated flux as the primary unknowns, and that recovers the traditional Fourier-type diffusion as a limiting case. We design a novel semi-implicit discretization for this formulation on staggered MAC grids and a geometric Multigrid-preconditioned Conjugate Gradients solver for efficient numerical solution. To highlight the efficacy of our method, we demonstrate end-to-end examples of elastic porous media simulated with the Material Point Method (MPM), and diffusion-driven Eulerian incompressible fluids.