Numerical Simulations Of Shock-induced Gas-liquid Two-phase Oblique Interface Deformation And Turbulent Mixing

As a common and complex physical phenomenon in practical engineering applications such as supersonic combustion,the interaction between shock wave and gasliquid interface has a wide range of engineering application background and important academic research value,and has always been one of the research hotspots in the field of computational fluid dynamics.The purpose of this paper is to study the evolution law of shock-induced gas-liquid two-phase oblique interface deformation and turbulent mixing phenomenon,provide theoretical support for practical engineering applications such as optimal control of fuel mixing effect in supersonic combustion,and lay the foundation for subsequent three-dimensional high-precision numerical research.In this paper,Volume of Fluid(VOF)model and large eddy simulation(LES)were utilized to investigate the evolution of shock induced deformation and turbulent mixing at the gas-liquid two-phase inclined interface.The interaction between the classical planar shock and non-planar shock and the gas-liquid oblique interface was discussed.The effects of shock Mach number,initial dip angle of shock front,initial disturbance amplitude of shock front,initial disturbance wavelength of shock front and initial tilt angle of gas-liquid two-phase inclined interface on the interface deformation and turbulent mixing are analyzed.The main contents are as follows:(a)The gas-liquid two-phase flow model and numerical calculation method based on VOF and large eddy simulation are established.The initial flow field parameters of classical planar shock and non-planar shock are determined by using shock tube theory.(b)Numerical simulations have been carried out on the interaction between classical planar shock wave and gas-liquid two-phase oblique interface.The effects of shock Mach number,initial inclination angle of shock front and initial inclination angle of gas-liquid two-phase oblique interface on the development of interface deformation and turbulent mixing are analyzed.(c)Numerical simulations of the interaction behavior between non-planar shock wave and gas-liquid two-phase oblique interface are carried out.Relying on the MATLAB platform,the initial geometric conditions of sinusoidal shock front are constructed.The parameters of shock Mach number,initial inclination angle of gas-liquid two-phase oblique interface,amplitude and wavelength of initial disturbance of shock front have been studied.The influence mechanism of initial conditions on non-planar shock-induced oblique interface deformation and turbulent mixing is revealed.The main conclusions are as follows:(1)The width of turbulent mixing zone increases nearly linearly with time under a given condition.(2)Increasing the shock Mach number has the most significant effect on improving the interface deformation and turbulent mixing degree,and adjusting the shock Mach number is the first choice to control the interface deformation and turbulent mixing degree.(3)With the increase of the inclined angle of the gas-liquid interface,the asymmetry of the convex structure distribution on both sides of the interface becomes more obvious,and the width of the turbulent mixing zone increases at the same time.(4)In the case of classical planar shock,adjusting the initial inclination of shock front will affect the development of interface deformation and turbulent mixing by changing the intensity of normal shock.(5)In the case of non-planar shock,with the increase of the initial disturbance amplitude,the forming speed of the convex structure on the phase interface accelerates,the bending degree of the phase interface increases in the initial stage of turbulent mixing,the time of the convex structure falling off and breaking up on the phase interface is advanced in the later stage of turbulent mixing,and the width of the turbulent mixing zone increases at the same time.(6)In the case of non-planar shock,with the decrease of the initial disturbance wavelength,the width of the turbulent mixing zone increases at the same time,the convex structure on the phase interface falls off and breaks earlier,the difference of the width of the turbulent mixing zone between different working conditions increases,and the interface deformation and turbulent mixing degree increase.The innovation of this paper lies in: Based on MATLAB and fluid simulation platform,a special sinusoidal shock front geometric model is constructed;Based on the component model of multiphase fluid,a numerical method is developed to study the interaction between nonplanar shock wave and gas-liquid two-phase oblique interface.The influence of initial disturbance conditions on the evolution of interface deformation and turbulent mixing degree is revealed.