View-dependent adaptive cloth simulation

This paper describes a method for view-dependent cloth simulation using dynamically adaptive mesh refinement and coarsening. Given a prescribed camera motion, the method adjusts the criteria controlling refinement to account for visibility and apparent size in the camera's view. Objectionable dynamic artifacts are avoided by anticipative refinement and smoothed coarsening. This approach preserves the appearance of detailed cloth throughout the animation while avoiding the wasted effort of simulating details that would not be discernible to the viewer. The computational savings realized by this method increase as scene complexity grows, producing a 2× speed-up for a single character and more than 4× for a small group

Strain injection techniques in dynamic fracture modeling

A computationally affordable modeling of dynamic fracture phenomena is performed in this study by using strain injection techniques and Finite Elements with Embedded strong discontinuities (E-FEM). In the present research, classical strain localization and strong discontinuity approaches are considered by injecting discontinuous strain and displacement modes in the finite element formulation without an increase of the total number of degrees of freedom. Following the Continuum Strong Discontinuity Approach (CSDA), stress–strain constitutive laws can be employed in the context of fracture phenomena and, therefore, the methodology remains applicable to a wide number of continuum mechanics models. The position and orientation of the displacement discontinuity is obtained through the solution of a crack propagation problem, i.e. the crack path field, based on the distribution of localized strains. The combination of the above mentioned approaches is envisaged to avoid stress-locking and directional mesh bias phenomena. Dynamic simulations are performed increasing the loading rate up to the appearance of crack branching, and the variation in terms of failure modes is investigated as well as the influence of the strain injection together with the crack path field algorithm. Objectivity of the presented methodology with respect to the spatial and temporal discretization is analyzed in terms of the dissipated energy during the fracture process. The dissipation at the onset of branching is studied for different loading rate conditions and is linked to the experimental maximum velocity observed before branching takes place.Fil: Lloberas Valls, Oriol. Universidad Politecnica de Catalunya; España. Centre Internacional de Metodes Numerics en Enginyeria; EspañaFil: Huespe, Alfredo Edmundo. Centre Internacional de Metodes Numerics en Enginyeria; España. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Santa Fe. Centro de Investigaciones en Métodos Computacionales. Universidad Nacional del Litoral. Centro de Investigaciones en Métodos Computacionales; ArgentinaFil: Oliver, J.. Centre Internacional de Metodes Numerics en Enginyeria; España. Universidad Politecnica de Catalunya; EspañaFil: Dias, I.F.. Laboratório Nacional de Engenharia Civil; Portuga

Three-dimensional CFD simulations with large displacement of the geometries using a connectivity-change moving mesh approach

This paper deals with three-dimensional (3D) numerical simulations involving 3D moving geometries with large displacements on unstructured meshes. Such simulations are of great value to industry, but remain very time-consuming. A robust moving mesh algorithm coupling an elasticity-like mesh deformation solution and mesh optimizations was proposed in previous works, which removes the need for global remeshing when performing large displacements. The optimizations, and in particular generalized edge/face swapping, preserve the initial quality of the mesh throughout the simulation. We propose to integrate an Arbitrary Lagrangian Eulerian compressible flow solver into this process to demonstrate its capabilities in a full CFD computation context. This solver relies on a local enforcement of the discrete geometric conservation law to preserve the order of accuracy of the time integration. The displacement of the geometries is either imposed, or driven by fluid–structure interaction (FSI). In the latter case, the six degrees of freedom approach for rigid bodies is considered. Finally, several 3D imposed-motion and FSI examples are given to validate the proposed approach, both in academic and industrial configurations

Hierarchical path-finding for Navigation Meshes (HNA*)

Path-finding can become an important bottleneck as both the size of the virtual environments and the number of agents navigating them increase. It is important to develop techniques that can be efficiently applied to any environment independently of its abstract representation. In this paper we present a hierarchical NavMesh representation to speed up path-finding. Hierarchical path-finding (HPA*) has been successfully applied to regular grids, but there is a need to extend the benefits of this method to polygonal navigation meshes. As opposed to regular grids, navigation meshes offer representations with higher accuracy regarding the underlying geometry, while containing a smaller number of cells. Therefore, we present a bottom-up method to create a hierarchical representation based on a multilevel k-way partitioning algorithm (MLkP), annotated with sub-paths that can be accessed online by our Hierarchical NavMesh Path-finding algorithm (HNA*). The algorithm benefits from searching in graphs with a much smaller number of cells, thus performing up to 7.7 times faster than traditional A¿ over the initial NavMesh. We present results of HNA* over a variety of scenarios and discuss the benefits of the algorithm together with areas for improvement.Peer ReviewedPostprint (author's final draft

Three-dimensional dynamic rupture simulation with a high-order discontinuous Galerkin method on unstructured tetrahedral meshes

Accurate and efficient numerical methods to simulate dynamic earthquake rupture and wave propagation in complex media and complex fault geometries are needed to address fundamental questions in earthquake dynamics, to integrate seismic and geodetic data into emerging approaches for dynamic source inversion, and to generate realistic physics-based earthquake scenarios for hazard assessment. Modeling of spontaneous earthquake rupture and seismic wave propagation by a high-order discontinuous Galerkin (DG) method combined with an arbitrarily high-order derivatives (ADER) time integration method was introduced in two dimensions by de la Puente et al. (2009). The ADER-DG method enables high accuracy in space and time and discretization by unstructured meshes. Here we extend this method to three-dimensional dynamic rupture problems. The high geometrical flexibility provided by the usage of tetrahedral elements and the lack of spurious mesh reflections in the ADER-DG method allows the refinement of the mesh close to the fault to model the rupture dynamics adequately while concentrating computational resources only where needed. Moreover, ADER-DG does not generate spurious high-frequency perturbations on the fault and hence does not require artificial Kelvin-Voigt damping. We verify our three-dimensional implementation by comparing results of the SCEC TPV3 test problem with two well-established numerical methods, finite differences, and spectral boundary integral. Furthermore, a convergence study is presented to demonstrate the systematic consistency of the method. To illustrate the capabilities of the high-order accurate ADER-DG scheme on unstructured meshes, we simulate an earthquake scenario, inspired by the 1992 Landers earthquake, that includes curved faults, fault branches, and surface topography