A stochastic continuum damage model for dynamic fracture analysis of quasi-brittle materials using asynchronous Spacetime Discontinuous Galerkin (aSDG) method

The microstructural design has an essential effect on the fracture response of brittle materials. We present a stochastic bulk damage formulation to model dynamic brittle fracture. This model is compared with a similar interfacial model for homogeneous and heterogeneous materials. The damage models are rate-dependent, and the corresponding damage evolution includes delay effects. The evolution equation specifies the rate at which damage tends to its quasi-static limit. The relaxation time of the model introduces an intrinsic length scale for dynamic fracture and addresses the mesh sensitivity problem of earlier damage models with much less computational efforts. The ordinary differential form of the damage equation makes this remedy quite simple and enables capturing the loading rate sensitivity of strain-stress response. A stochastic field is defined for material cohesion and fracture strength to involve microstructure effects in the proposed formulations. The statistical fields are constructed through the Karhunen-Loeve (KL) method.An advanced asynchronous Spacetime Discontinuous Galerkin (aSDG) method is used to discretize the final system of coupled equations. Local and asynchronous solution process, linear complexity of the solution versus the number of elements, local recovery of balance properties, and high spatial and temporal orders of accuracy are some of the main advantages of the aSDG method.Several numerical examples are presented to demonstrate mesh insensitivity of the method and the effect of boundary conditions on dynamic fracture patterns. It is shown that inhomogeneity greatly differentiates fracture patterns from those of a homogeneous rock, including the location of zones with maximum damage. Moreover, as the correlation length of the random field decreases, fracture patterns resemble angled-cracks observed in compressive rock fracture. The final results show that a stochastic bulk damage model produces more realistic results in comparison with a homogenizes model

Theoretical and numerical comparison of hyperelastic and hypoelastic formulations for Eulerian non-linear elastoplasticity

The aim of this paper is to compare a hyperelastic with a hypoelastic model describing the Eulerian dynamics of solids in the context of non-linear elastoplastic deformations. Specifically, we consider the well-known hypoelastic Wilkins model, which is compared against a hyperelastic model based on the work of Godunov and Romenski. First, we discuss some general conceptual differences between the two approaches. Second, a detailed study of both models is proposed, where differences are made evident at the aid of deriving a hypoelastic-type model corresponding to the hyperelastic model and a particular equation of state used in this paper. Third, using the same high order ADER Finite Volume and Discontinuous Galerkin methods on fixed and moving unstructured meshes for both models, a wide range of numerical benchmark test problems has been solved. The numerical solutions obtained for the two different models are directly compared with each other. For small elastic deformations, the two models produce very similar solutions that are close to each other. However, if large elastic or elastoplastic deformations occur, the solutions present larger differences.Comment: 14 figure

Numerical modelling of wave propagation phenomena in thermo-poroelastic media via discontinuous Galerkin methods

We present and analyze a high-order discontinuous Galerkin method for the space discretization of the wave propagation model in thermo-poroelastic media. The proposed scheme supports general polytopal grids. Stability analysis and $hp$-version error estimates in suitable energy norms are derived for the semi-discrete problem. The fully-discrete scheme is then obtained based on employing an implicit Newmark-$\beta$ time integration scheme. A wide set of numerical simulations is reported, both for the verification of the theoretical estimates and for examples of physical interest. A comparison with the results of the poroelastic model is provided too, highlighting the differences between the predictive capabilities of the two models

ADER-WENO Finite Volume Schemes with Space-Time Adaptive Mesh Refinement

We present the first high order one-step ADER-WENO finite volume scheme with Adaptive Mesh Refinement (AMR) in multiple space dimensions. High order spatial accuracy is obtained through a WENO reconstruction, while a high order one-step time discretization is achieved using a local space-time discontinuous Galerkin predictor method. Due to the one-step nature of the underlying scheme, the resulting algorithm is particularly well suited for an AMR strategy on space-time adaptive meshes, i.e.with time-accurate local time stepping. The AMR property has been implemented 'cell-by-cell', with a standard tree-type algorithm, while the scheme has been parallelized via the Message Passing Interface (MPI) paradigm. The new scheme has been tested over a wide range of examples for nonlinear systems of hyperbolic conservation laws, including the classical Euler equations of compressible gas dynamics and the equations of magnetohydrodynamics (MHD). High order in space and time have been confirmed via a numerical convergence study and a detailed analysis of the computational speed-up with respect to highly refined uniform meshes is also presented. We also show test problems where the presented high order AMR scheme behaves clearly better than traditional second order AMR methods. The proposed scheme that combines for the first time high order ADER methods with space--time adaptive grids in two and three space dimensions is likely to become a useful tool in several fields of computational physics, applied mathematics and mechanics.Comment: With updated bibliography informatio