TWO-DIMENSIONAL NUMERICAL SIMULATION OF THE SHOCK WAVE-PARTICLE LAYER INTERACTION USING CARTESIAN GRID METHOD

Author(s):  
D. A. SIDORENKO ◽  
◽  
P. S. UTKIN ◽  

The process of dispersion of particles behind a shock wave (SW) propagating along the dust layer is one of the stages of the complex process of dust-layered detonation [1]. The paper presents the results of direct numerical simulation of the problem of interaction of a sliding SW with a layer of moving and interacting particles on an impermeable surface.

Author(s):  
YA. E. POROSHYNA ◽  
◽  
P. S. UTKIN ◽  

The problem of shock wave - dense particle layer interaction is a fundamental basis for the study of a more complex process of dust explosion or dust-layered detonation. The work presents results of numerical simulation of the experiment on interaction of an SW with particles layer deposited on the impermeable wall.


Materials ◽  
2021 ◽  
Vol 14 (13) ◽  
pp. 3553
Author(s):  
Dengwang Wang ◽  
Yong Gao ◽  
Sheng Wang ◽  
Jie Wang ◽  
Haipeng Li

Carbon/Phenolic (C/P), a typical anisotropic material, is an important component of aerospace and often used to protect the thermodynamic effects of strong X-ray radiation. In this paper, we establish the anisotropic elastic-plastic constitutive model, which is embedded in the in-house code “RAMA” to simulate a two-dimensional thermal shock wave induced by X-ray. Then, we compare the numerical simulation results with the thermal shock wave stress generated by the same strong current electron beam via experiment to verify the correctness of the numerical simulation. Subsequently, we discuss and analyze the rules of thermal shock wave propagation in C/P material by further numerical simulation. The results reveal that the thermal shock wave represents different shapes and mechanisms by the radiation of 1 keV and 3 keV X-rays. The vaporization recoil phenomenon appears as a compression wave under 1 keV X-ray irradiation, and X-ray penetration is caused by thermal deformation under 3 keV X-ray irradiation. The thermal shock wave propagation exhibits two-dimensional characteristics, the energy deposition of 1 keV and 3 keV both decays exponentially, the energy deposition of 1 keV-peak soft X-ray is high, and the deposition depth is shallow, while the energy deposition of 3 keV-peak hard X-ray is low, and the deposition depth is deep. RAMA can successfully realize two-dimensional orthotropic elastoplastic constitutive relation, the corresponding program was designed and checked, and the calculation results for inspection are consistent with the theory. This study has great significance in the evaluation of anisotropic material protection under the radiation of intense X-rays.


2019 ◽  
Vol 221 ◽  
pp. 01021
Author(s):  
Aleksandr Kraus ◽  
Evgeny Kraus ◽  
Ivan Shabalin

A two-dimensional and three-dimensional non-stationary problem of the interaction of a homogeneous impactor and a heterogeneous structure made of steel and ceramics and placed in a Kevlar pocket is considered. The model of the human body is a plate of gelatine with cylindrical inserts-imitators of human bones. The results of numerical simulation using different approaches for describing heterogeneous media are compared. On the basis of direct numerical simulation, it is shown that the gradient armor plate (steel + B4C) has the best weight and size parameters.


2019 ◽  
Vol 877 ◽  
pp. 167-195 ◽  
Author(s):  
Feng-Yuan Zuo ◽  
Antonio Memmolo ◽  
Guo-ping Huang ◽  
Sergio Pirozzoli

Direct numerical simulation of the Navier–Stokes equations is carried out to investigate the interaction of a conical shock wave with a turbulent boundary layer developing over a flat plate at free-stream Mach number $M_{\infty }=2.05$ and Reynolds number $Re_{\unicode[STIX]{x1D703}}\approx 630$, based on the upstream boundary layer momentum thickness. The shock is generated by a circular cone with half opening angle $\unicode[STIX]{x1D703}_{c}=25^{\circ }$. As found in experiments, the wall pressure exhibits a distinctive N-wave signature, with a sharp peak right past the precursor shock generated at the cone apex, followed by an extended zone with favourable pressure gradient, and terminated by the trailing shock associated with recompression in the wake of the cone. The boundary layer behaviour is strongly affected by the imposed pressure gradient. Streaks are suppressed in adverse pressure gradient (APG) zones, but re-form rapidly in downstream favourable pressure gradient (FPG) zones. Three-dimensional mean flow separation is only observed in the first APG region associated with the formation of a horseshoe vortex, whereas the second APG region features an incipient detachment state, with scattered spots of instantaneous reversed flow. As found in canonical geometrically two-dimensional wedge-generated shock–boundary layer interactions, different amplification of the turbulent stress components is observed through the interacting shock system, with approach to an isotropic state in APG regions, and to a two-component anisotropic state in FPG. The general adequacy of the Boussinesq hypothesis is found to predict the spatial organization of the turbulent shear stresses, although different eddy viscosities should be used for each component, as in tensor eddy-viscosity models, or in full Reynolds stress closures.


AIP Advances ◽  
2020 ◽  
Vol 10 (2) ◽  
pp. 025034
Author(s):  
Balaji Jayaraman ◽  
Saadbin Khan

2013 ◽  
Vol 5 (3) ◽  
pp. 435-445
Author(s):  
M. S. I. Mallik ◽  
M. A. Uddin ◽  
M. A. Rahman

Direct numerical simulation (DNS) in two-dimensional homogeneous isotropic turbulence is performed by using the Spectral method at a Reynolds number Re = 1000 on a uniformly distributed grid points. The Reynolds number is low enough that the computational grid is capable of resolving all the possible turbulent scales. The statistical properties in the computed flow field show a good agreement with the qualitative behavior of decaying turbulence. The behavior of the flow structures in the computed flow field also follow the classical idea of the fluid flow in turbulence. Keywords: Direct numerical simulation, Isotropic turbulence, Spectral method. © 2013 JSR Publications. ISSN: 2070-0237 (Print); 2070-0245 (Online). All rights reserved. doi:http://dx.doi.org/10.3329/jsr.v5i3.12665 J. Sci. Res. 5 (3), 435-445 (2013)  


1997 ◽  
Vol 334 ◽  
pp. 353-379 ◽  
Author(s):  
KRISHNAN MAHESH ◽  
SANJIVA K. LELE ◽  
PARVIZ MOIN

Direct numerical simulation and inviscid linear analysis are used to study the interaction of a normal shock wave with an isotropic turbulent field of vorticity and entropy fluctuations. The role of the upstream entropy fluctuations is emphasized. The upstream correlation between the vorticity and entropy fluctuations is shown to strongly influence the evolution of the turbulence across the shock. Negative upstream correlation between u′ and T′ is seen to enhance the amplification of the turbulence kinetic energy, vorticity and thermodynamic fluctuations across the shock wave. Positive upstream correlation has a suppressing effect. An explanation based on the relative effects of bulk compression and baroclinic torque is proposed, and a scaling law is derived for the evolution of vorticity fluctuations across the shock. The validity of Morkovin's hypothesis across a shock wave is examined. Linear analysis is used to suggest that shock-front oscillation would invalidate the relation between urms and Trms, as expressed by the hypothesis.


Sign in / Sign up

Export Citation Format

Share Document