A Simulation of the Liquid Shock and Cavitation Based on a Multi-Equation Model

2016 ◽  
Vol 13 (04) ◽  
pp. 1641010
Author(s):  
Yang-Yao Niu

In this paper, an unsteady preconditioning formulation for multi-phase flows with arbitrary equation of state based on the approximated Riemann solver is developed for multi-phase flows at all speed. This paper considers a homogeneous two-phase multi-equation mixture model with the assumption of kinematics and thermodynamics equilibriums. The thermodynamics behaviors of liquid phase, vapor phase and their phase transitional process are described by a temperature-dependent hybrid equation of state. Benchmark test cases, including one-dimensional (1D) condensation shock in the cavitated nozzle and two-dimensional (2D) cavitated blunt body problem, demonstrate accurate capturing of interfaces, shock waves and cavitation zones.

Author(s):  
Lucio Raimondo ◽  
Lorenzo Iannucci ◽  
Paul Robinson ◽  
Paul T. Curtis ◽  
Garry M. Wells

This paper presents part of an ongoing programme of work on high velocity impact modelling on composite targets. The modelling approach aims to link existing low velocity constitutive failure models, including delamination modelling, with relevant orthotropic Equations Of State models. A methodology for predicting the Hugoniot states (shock velocity vs. particle velocity) of multi-phase materials at high compression is presented. The Gruneisen parameter of the mixture is also derived. The proposed approach is a step toward a full thermodynamic virtual characterisation of untested multi-phase materials, when tabulated shock data for the constituents is available [1]. Other approaches have been proposed [2], [3]; however, they require complex Finite Element coding and iterative procedures and are limited to two-phase materials. The approach is critically discussed in relation to shock data derived from existing flyer plate impact test data. An orthotropic Equation of State [4] has also been implemented into the LS-DYNA3D code. A flyer plate test is simulated using the implemented model, and with material parameters derived using the theory of mixture approach. The current orthotropic Equation of State formulation is discussed, within the limitation of classical Lagrangian FE techniques. Additionally, conclusions are drawn on the logical next step to model high velocity angled impacts onto orthotropic targets.


1971 ◽  
Vol 38 (2) ◽  
pp. 363-370 ◽  
Author(s):  
L. W. Woodruff ◽  
W. H. Giedt ◽  
J. L. Hesse

The present study was undertaken to investigate the applicability of one-dimensional computer programs for predicting the response of materials exposed to rapid surface heating produced by a pulsed electron beam. A specific objective was to determine modifications necessary to account for surface melting and spall. Measured values of mass loss, impulse, and internal stresses in lead and gold were satisfactorily predicted using a one-dimensional finite-difference program which simultaneously solved the conservation equations and a hydrodynamic equation of state. Required program modifications included: (a) specifying spall to occur to the depth below the target surface at which the energy deposited was sufficient to initiate melting, and (b) revising the equation of state for the material in the two-phase region.


2004 ◽  
Vol 15 (1) ◽  
pp. 79-102 ◽  
Author(s):  
OLIVIER BACONNEAU ◽  
JAN BOUWE VAN DEN BERG ◽  
CLAUDE-MICHEL BRAUNER ◽  
JOSEPHUU HULSHOF

We study travelling wave solutions of a one-dimensional two-phase Free Boundary Problem, which models premixed flames propagating in a gaseous mixture with dust. The model combines diffusion of mass and temperature with reaction at the flame front, the reaction rate being temperature dependent. The radiative effects due to the presence of dust account for the divergence of the radiative flux entering the equation for temperature. This flux is modelled by the Eddington equation. In an appropriate limit the divergence of the flux takes the form of a nonlinear heat loss term. The resulting reduced model is able to capture a hysteresis effect that appears if the amount of fuel in front of the flame, or equivalently, the adiabatic temperature is taken as a control parameter.


2015 ◽  
Vol 25 (9) ◽  
pp. 795-817 ◽  
Author(s):  
Mika P. Jarvinen ◽  
A. E. P. Kankkunen ◽  
R. Virtanen ◽  
P. H. Miikkulainen ◽  
V. P. Heikkila

Sign in / Sign up

Export Citation Format

Share Document