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

Multi-scale control variate methods for uncertainty quantification in kinetic equations

Kinetic equations play a major rule in modeling large systems of interacting particles. Uncertainties may be due to various reasons, like lack of knowledge on the microscopic interaction details or incomplete informations at the boundaries. These uncertainties, however, contribute to the curse of dimensionality and the development of efficient numerical methods is a challenge. In this paper we consider the construction of novel multi-scale methods for such problems which, thanks to a control variate approach, are capable to reduce the variance of standard Monte Carlo techniques

A realizable filtered intrusive polynomial moment method

Intrusive uncertainty quantification methods for hyperbolic problems exhibit spurious oscillations at shocks, which leads to a significant reduction of the overall approximation quality. Furthermore, a challenging task is to preserve hyperbolicity of the gPC moment system. An intrusive method which guarantees hyperbolicity is the intrusive polynomial moment (IPM) method, which performs the gPC expansion on the entropy variables. The method, while still being subject to oscillations, requires solving a convex optimization problem in every spatial cell and every time step. The aim of this work is to mitigate oscillations in the IPM solution by applying filters. Filters reduce oscillations by damping high order gPC coefficients. Naive filtering, however, may lead to unrealizable moments, which means that the IPM optimization problem does not have a solution and the method breaks down. In this paper, we propose and analyze two separate strategies to guarantee the existence of a solution to the IPM problem. First, we propose a filter which maintains realizability by being constructed from an underlying Fokker-Planck equation. Second, we regularize the IPM optimization problem to be able to cope with non-realizable gPC coefficients. Consequently, standard filters can be applied to the regularized IPM method. We demonstrate numerical results for the two strategies by investigating the Euler equations with uncertain shock structures in one- and two-dimensional spatial settings. We are able to show a significant reduction of spurious oscillations by the proposed filters