Drucker-prager elastoplasticity for sand animation

SummaryTo improve the computational efficiency of the numerical manifold method for discontinuous deformation simulations, a spatial‐domain coupled explicit‐implicit time integration algorithm is proposed. A subdomain partition algorithm based on a super manifold element is developed for the numerical manifold method to simulate dynamic motions of blocky rock mass. In different subdomains, explicit or implicit time integration method is employed respectively based on its contact and motion status. These subdomains interact through assembling the corresponding explicit or implicit time integration‐based matrices of different rock blocks. The computational efficiency of the discontinuity system under dynamic loading is improved by partially diagonalizing the global matrices. Two verification examples of a sliding block along an inclined plane under a horizontal acceleration excitation and a multiblock system acted on by dynamic forces are studied to examine the accuracy of the proposed numerical method, respectively. A highly fractured rock mass situated on an inclined slope subjected to seismic excitations is then studied to show the computational efficiency of the developed algorithm. The simulated results are in good agreement with those from the versions using purely implicit or explicit time integration algorithm for the numerical manifold method. The computational efficiency is shown to be higher using the proposed algorithm, which demonstrates its potential for application in dynamic analysis of highly fractured rock masses.

A coupled implicit-explicit time integration method for compressible unsteady flows

This study investigates the full penetration simulation of piles from the ground surface, focusing on frictional contact modeling without mesh distortion. To overcome issues related to mesh distortion and improve solution convergence, the Arbitrary Lagrangian–Eulerian (ALE) adaptive mesh technique was implemented within both explicit and implicit time integration schemes. The numerical model was validated against field experiments conducted at Bothkennar, Scotland, using the Imperial College instrumented displacement pile (ICP) in soft clay, where the soil behavior was effectively represented using the modified Cam-Clay model and the Mohr–Coulomb model. The primary objectives of this study are to evaluate the ALE method performance in handling mesh distortion; analyze the effects of soil–pile interface friction, pile dimensions, and various dilation angles on pile resistance; and compare the effectiveness of explicit and implicit time integration schemes in terms of stability, computational efficiency, and solution accuracy. The ALE method effectively modeled pile penetration in Bothkennar clay, validating the numerical model against field experiments. Comparative analysis revealed the explicit time integration method as more robust and computationally efficient, particularly for complex soil–pile interactions with higher friction coefficients.

Implicit–explicit time integration of a high-order particle-in-cell method with hyperbolic divergence cleaning

Two- and three-dimensional fluid-structure interaction computer programs for the simulation of nonlinear dynamics were developed and applied to a number of problems. The programs were created by coupling Arbitrary Lagrangian-Eulerian finite volume fluid dynamics programs with strictly Lagrangian finite element structural dynamics programs. The resulting coupled programs can use either fully explicit or implicit time integration. The implicit time integration is accomplished by iterations of the fluid dynamics pressure solver and the structural dynamics system solver. The coupled programs have been used to solve problems involving incompressible fluids, membrane and shell elements, compressible multiphase flows, explosions in both air and water, and large displacements. In this paper, we present the approach used for the coupling and describe test problems that verify the two-dimensional programs against an experiment and an analytical linear problem. The experiment involves an explosion underwater near an instrumented thin steel plate. The analytical linear problem is the vibration of an infinite cylinder surrounded by an incompressible fluid to a given radius.

This work presents the development of an explicit/implicit particle finite element method (PFEM) for the 2D modeling of metal cutting processes. The purpose is to study the efficiency of implicit and explicit time integration schemes in terms of precision, accuracy and computing time. The formulation for implicit and explicit time marching schemes is developed, and a detailed study on the explicit solution steps is presented. The PFEM remeshing procedures for insertion and removal of particles have been improved to model the multiple scales of time and/or space of the solution. The detection and treatment of the rigid tool contact are presented for both, implicit and explicit schemes. The performance of explicit/implicit integration is studied with a set of different two-dimensional orthogonal cutting tests of AISI 4340 steel at cutting speeds ranging from 1 m/s up to 30 m/s. It was shown that if the correct selection of the time integration scheme is made, the computing time can decrease up to 40 times. It allows us to affirm that the computing time of the PFEM simulations can be excessive due to the used time marching scheme independently of the meshing process. As a practical result, a set of recommendations to select the time integration schemes for a given cutting speed are given. This is intended to minimize one of the negative constraints pointed out by the industry when using metal cutting simulators.

Solutions of the interfacial shear and normal stresses in plate flexural-strengthened beams based on different complementary strain energy assumptions

In this paper high order time integration schemes are investigated for the integration of uid-structure interaction (FSI) model problem with deforming domains. The integration algorithm consists of a partitioned scheme in which a combination of implicit and explicit (IMEX) high order Runge-Kutta time integration schemes is used: the uid and structure systems are decoupled at their interface and both systems are integrated using the implicit system under a given boundary condition at the interface. The inuence of this boundary condition (the coupling term) is integrated using the explicit Runge-Kutta scheme. The uid dynamics are written in the arbitrary Lagrange Eulerian formulation to cope with the deforming uid mesh. An expression is proposed for the mesh face velocities of the deforming mesh such that the discrete geometric conservation law is satised. For the model problem numerical investigations are made into the eciency of these high order time integration schemes compared to lower order methods. The IMEX scheme retains the order of the implicit and explicit schemes and obtains the same accuracy as a second order staggered scheme for less computational work.

In this collaborative research project between Pennsylvania State University, Colorado State University and Florida State University, we mainly focused on developing multi-resolution algorithms which are suitable to regional ocean modeling. We developed hybrid implicit and explicit adaptive multirate time integration method to solve systems of time-dependent equations that present two signi cantly di erent scales. We studied the e ects of spatial simplicial meshes on the stability and the conditioning of fully discrete approximations. We also studies adaptive nite element method (AFEM) based upon the Centroidal Voronoi Tessellation (CVT) and superconvergent gradient recovery. Some of these techniques are now being used by geoscientists(such as those at LANL).

In this paper a robust and effective 4–node shell element for the structural analysis of thin structures is described. A Hu–Washizu functional with independent displacements, stress resultants and shell strains is the variational basis of the formulation. With application of three basic assumptions within the FE formulation a stabilized one point integrated element comes out automatically. Besides static computations it can be used for dynamic problems with explicit or implicit time integration without further adaptions. The element possesses the correct rank and fulfills the membrane and bending patch test. The stabilization terms can be integrated analytically and do not require an input of problem dependent control parameters. In contrast to several approaches the results are not path dependent for hyperelastic material behavior. This unphysical behavior follows when the material matrix of the last load increment is used in the stabilization matrix to preserve quadratic convergence in the Newton-Raphson iteration. The element possesses the well-known robustness of Hu–Washizu elements with the possibility of very large load steps in nonlinear applications. Another essential advantage of present development is the application to FE2 computations where the computing times are dominated by the solution of the micro problems at the shell integration points. Depending on the used interpolation for the shell strains the new element is approximately by a factor 4 or 9 faster than the corresponding version with full integration.

Stabilization of the Lax-Wendroff method and a generalized one-step Runge-Kutta method for hyperbolic initial-value problems

Dynamic burst pressure analysis of cylindrical shells based on average shear stress yield criterion

The present thesis describes the development of a computational method for the numerical simulation of compressible flows in aerodynamics. The main goal is to study high-order discretization schemes that can potentially deliver a better computational efficiency. The spatial part of the conservation equations is discretized with the finite difference method. In particular, Summation By Parts (SBP) finite difference operators and treatment of boundary conditions via Simultaneous Approximation Terms (SAT) have been investigated. An SBP finite difference operator is essentially a (high-order) centered finite difference scheme with a specific closure at the boundary. The SAT are penalty-like terms that enforce the boundary conditions weakly and are used to augment the SBP schemes. The combination of SBP operators with SAT boundary terms constitutes schemes that, for the linearized equations governing a smooth flow, are provably stable (thus convergent) schemes. The discretization in time is carried out with high-order schemes as well. Both implicit and explicit time integration schemes have been investigated. The system of non-linear equations that arises when using the implicit time integration schemes is solved via (a damped) Newton's method. On the target grid, the initial guess required for Newton's method is obtained by performing grid sequencing combined a series of explicit relaxation iterations on each grid level. The flow equations and the equation of the turbulence model (if present) are solved in a fully-coupled manner. Numerical simulations for inviscid flow, laminar flow, as well as simulations based on DNS, LES and RANS have been performed. Results show that the described computational method is particularly well-suited for smooth external flows about relatively simple geometries, solved with the DNS and LES approach. For these types of flows the gain in efficiency obtained by using high-order schemes is clear (more than a factor of 5). For RANS, however, the comparison with a state-of-the-art second order finite volume code suggests that solving the turbulence and the flow equations fully coupled with implicit schemes may not be the most efficient iteration strategy. On the other hand, the numerical simulations of non-smooth/internal flows highlighted the need of further development of the SBP-SAT framework.

In this paper, we study 3-D multiparameter full waveform inversion (FWI) in viscoelastic media based on the generalized Maxwell/Zener body including arbitrary number of attenuation mechanisms. We present a frequency-domain energy analysis to establish the stability condition of a full anisotropic viscoelastic system, according to zero-valued boundary condition and the elastic–viscoelastic correspondence principle: the real-valued stiffness matrix becomes a complex-valued one in Fourier domain when seismic attenuation is taken into account. We develop a least-squares optimization approach to linearly relate the quality factor with the anelastic coefficients by estimating a set of constants which are independent of the spatial coordinates, which supplies an explicit incorporation of the parameter <it>Q</it> in the general viscoelastic wave equation. By introducing the Lagrangian multipliers into the matrix expression of the wave equation with implicit time integration, we build a systematic formulation of multiparameter FWI for full anisotropic viscoelastic wave equation, while the equivalent form of the state and adjoint equation with explicit time integration is available to be resolved efficiently. In particular, this formulation lays the foundation for the inversion of the parameter <it>Q</it> in the time domain with full anisotropic viscoelastic properties. In the 3-D isotropic viscoelastic settings, the anelastic coefficients and the quality factors using bulk and shear moduli parametrization can be related to the counterparts using <it>P</it> and <it>S</it> velocity. Gradients with respect to any other parameter of interest can be found by chain rule. Pioneering numerical validations as well as the real applications of this most generic framework will be carried out to disclose the potential of viscoelastic FWI when adequate high-performance computing resources and the field data are available.

Understanding the flow of fluid, either liquid or gas, through and around solid bodies has challenged man since the dawn of scientific inquiry. Many of the great minds of science and math have progressively built up a hierarchy of fluid models. This report is concerned with the computational modeling of turbulent flow around aerodynamic bodies such as planes and wind turbines. In this case, viscous effects near the solid bodies create very thin boundary layers that yield highly anisotropic (gradients normal to the surface may be 10^6 times larger than gradients along the surface) solutions to the governing non-linear partial differential equations (PDE); the Navier-Stokes equations. Furthermore, turbulent flows develop extremely broad ranges of length and time scales.This disparity motivates the use of discretization methods capable of employing adaptivity and implicit time integration. The combination of these features (non-linear, anisotropy, adaptivity, andimplicit) dramatically raise the complexity of the discretization, posing large challenges to efficient scalable parallel implementation. However, through careful design, the more complex algorithms can provide great reductions in computational cost relative to simpler methods (e.g., Cartesian grids with explicit time integration) that are easier to mate efficiently to hardware. In this report, we not only describe our approach but we also address the fact that while complex algorithms may never be as efficient flop-for-flop as simple methods, in the important measure of science-per-core-hour, they can still win big by making complex features like adaptivity and implicit methods as efficient and scalable as possible.