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

A NURBS-based Discontinuous Galerkin method for conservation laws with high-order moving meshes

International audienceThe objective of the present work is to develop a new numerical framework for simulations including moving bodies, in the specific context of high-order meshes consistent with Computer-Aided Design (CAD) representations. Thus, the proposed approach combines ideas from isogeometric analysis, able to handle exactly CAD-based geometries, and Discontinuous Galerkin (DG) methods with an Arbitrary Lagrangian-Eulerian (ALE) formulation , able to solve complex problems with moving grids. The resulting approach is a DG method based on rational Bézier elements, that can be easily constructed from Non-Uniform Rational B-Splines (NURBS), formulated in a general ALE setting. We focus here on applications in compressible aerodynamics, but the method could be applied to other models. Two verification exercises are conducted, to assess rigorously the properties of the method and the convergence rates for representations up to sixth order. Finally, two problems are analysed in depth, involving compressible Euler and Navier-Stokes equations , for an oscillating cylinder and a pitching airfoil. In particular, the convergence of flow characteristics is investigated, as well as the impact of using curved boundaries in the context of deformable domains

Heat Transfer Mechanism In Particle-Laden Turbulent Shearless Flows

Particle-laden turbulent flows are one of the complex flow regimes involved in a wide range of environmental, industrial, biomedical and aeronautical applications. Recently the interest has included also the interaction between scalars and particles, and the complex scenario which arises from the interaction of particle finite inertia, temperature transport, and momentum and heat feedback of particles on the flow leads to a multi-scale and multi-physics phenomenon which is not yet fully understood. The present work aims to investigate the fluid-particle thermal interaction in turbulent mixing under one-way and two-way coupling regimes. A recent novel numerical framework has been used to investigate the impact of suspended sub-Kolmogorov inertial particles on heat transfer within the mixing layer which develops at the interface of two regions with different temperature in an isotropic turbulent flow. Temperature has been considered a passive scalar, advected by the solenoidal velocity field, and subject to the particle thermal feedback in the two-way regime. A self-similar stage always develops where all single-point statistics of the carrier fluid and the suspended particles collapse when properly re-scaled. We quantify the effect of particle inertial, parametrized through the Stokes and thermal Stokes numbers, on the heat transfer through the Nusselt number, defined as the ratio of the heat transfer to the thermal diffusion. A scale analysis will be presented. We show how the modulation of fluid temperature gradients due to the statistical alignments of the particle velocity and the local carrier flow temperature gradient field, impacts the overall heat transfer in the two-way coupling regime

High-order gas-kinetic scheme with TENO class reconstruction for the Euler and Navier-Stokes equations

The high-order gas-kinetic scheme(HGKS) with WENO spatial reconstruction method has been extensively validated through many numerical experiments, demonstrating its superior accuracy efficiency, and robustness. Compared with WENO class schemes, TENO class schemes exhibit significantly improved robustness, low numerical dissipation and sharp discontinuity capturing. In this paper, two kinds of fifth-order HGKS with TENO class schemes are designed. One involves replacing WENO5 scheme with the TENO5 scheme in the conventional WENO5-GKS. WENO and TENO schemes only provide the non-equilibrium state values at the cell interface. The slopes of the non-equilibrium state along with the equilibrium values and slopes, are obtained by additional linear reconstruction. Another kind of TENO5-D GKS is similar to WENO5-AO GKS. Following a strong scale-separation procedure, a tailored novel ENO-like stencil selection strategy is proposed such that the high-order accuracy is restored in smooth regions by selecting the candidate reconstruction on the large stencil while the ENO property is enforced near discontinuities by adopting the candidate reconstruction from smooth small stencils. The such TENO schemes are TENO-AA and TENO-D scheme. The HGKS scheme based on WENO-AO or TENO-D reconstruction take advantage of the large stencil to provide point values and slopes of the non-equilibrium state. By dynamically merging the reconstructed non-equilibrium slopes, extra reconstruction of the equilibrium state at the beginning of each time step can be avoided. The simplified schemes have better robustness and efficiency than the conventional WENO5-GKS or TENO5-GKS. TENO-D GKS is also as easy to develop as WENO-AO GKS to high-order finite volume method for unstructured mesh.Comment: arXiv admin note: text overlap with arXiv:2304.05572; text overlap with arXiv:1905.08489 by other author