9,403 research outputs found
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
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
Discretely exact derivatives for hyperbolic PDE-constrained optimization problems discretized by the discontinuous Galerkin method
This paper discusses the computation of derivatives for optimization problems
governed by linear hyperbolic systems of partial differential equations (PDEs)
that are discretized by the discontinuous Galerkin (dG) method. An efficient
and accurate computation of these derivatives is important, for instance, in
inverse problems and optimal control problems. This computation is usually
based on an adjoint PDE system, and the question addressed in this paper is how
the discretization of this adjoint system should relate to the dG
discretization of the hyperbolic state equation. Adjoint-based derivatives can
either be computed before or after discretization; these two options are often
referred to as the optimize-then-discretize and discretize-then-optimize
approaches. We discuss the relation between these two options for dG
discretizations in space and Runge-Kutta time integration. Discretely exact
discretizations for several hyperbolic optimization problems are derived,
including the advection equation, Maxwell's equations and the coupled
elastic-acoustic wave equation. We find that the discrete adjoint equation
inherits a natural dG discretization from the discretization of the state
equation and that the expressions for the discretely exact gradient often have
to take into account contributions from element faces. For the coupled
elastic-acoustic wave equation, the correctness and accuracy of our derivative
expressions are illustrated by comparisons with finite difference gradients.
The results show that a straightforward discretization of the continuous
gradient differs from the discretely exact gradient, and thus is not consistent
with the discretized objective. This inconsistency may cause difficulties in
the convergence of gradient based algorithms for solving optimization problems
- …