1 / 29

BAROCLINIC TURBULENCE OF DISSOCIATING GAS AS A REASON FOR INSTABILITY OF BOW SHOCK WAVE

BAROCLINIC TURBULENCE OF DISSOCIATING GAS AS A REASON FOR INSTABILITY OF BOW SHOCK WAVE. A.S.Baryshnikov Physico-Technical Institute of RAS, Russia,. The history. The gases with great molecular weight

drake
Télécharger la présentation

BAROCLINIC TURBULENCE OF DISSOCIATING GAS AS A REASON FOR INSTABILITY OF BOW SHOCK WAVE

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. BAROCLINIC TURBULENCE OF DISSOCIATING GAS AS A REASON FOR INSTABILITY OF BOW SHOCK WAVE A.S.Baryshnikov Physico-Technical Institute of RAS, Russia,

  2. The history • The gases with great molecular weight is often used for high Much number modeling in free flight of body in ballistic ranges. One of them is gaseous CF2Cl2 , which has very high energy of dissociation what is impotent for modeling flight in ideal not reacting gas.

  3. Scheme of ballistic range receiver Powder gun One of visualization windows with Shlieren set up

  4. Problem • When velocity of body riches the amount at which Chemical dissociation (not exothermic reaction) takes the place • Just a disturbances on bow shock wave are arising. • The flow behind body becomes disturbed in such manner as for turbulent flow. • It’s important that such ”turbulence” is outside the wake flow region

  5. Fig.1.Instability of bow shock wave in front of segmental body. Flow Mach number: М0=3,9; Gas: CF2Cl2. Pressure in flow: Р0=3,8·104Pа. Wake turbulence Disturbance on bow shock wave Single vortex

  6. Fig.2.Instability of bow shock wave in front of segmental body. Flow Mach number: М0=3,9; Gas: CF2Cl2. Pressure in flow: Р0=4,3·104Pа.

  7. Applications • Reducing a drag of air vehicles (destabilization anddestruction of bow shock wave), At dropping of Temperature and so of heat flux to surface!!! • Controlling a flight of them, • Without any energy deposition (Compare with case of plasma control),

  8. Usual study region

  9. Fig.3.Estimation of one dimension spectrum of density fluctuations in investigated region of turbulent flow behind segmental body higher than axis. Body diameter d=27 mm, CF2Cl2; М0=6,8 , Р0=3,8·104Pа, length of realizations Θ =18 mm. I(ξ) Db -20 -40 100 200 ξ bin ξ is normalized wave number of fluctuations

  10. Fig.4.Estimation of one dimension spectrum of density fluctuations in investigated region of turbulent flow behind sphere body higher and lower than axis. Body diameter d=27 mm, CF2Cl2; М0=6,1 , Р0=8·104Pа, length of realizations Θ =41 mm. I(ξ) Db -10 -20 -30 -40 ξ bin 50 100 150 200

  11. Estimation of energy spectra • For isotropy turbulence we have

  12. Fig.5. Estimation of three dimension energy spectrum of turbulence in investigated flow region behind segmental body higher than axis. Body diameter d=27 mm, CF2Cl2; М0=6,8, Р0=3,8·104Pа, length of realizations Θ =18 mm. 1/rc

  13. Fig.6. Estimation of three dimension energy spectrum of turbulence in investigated flow region behind sphere body higher and lower than axis. Body diameter d=27 мм, CF2Cl2; М0=6,1, Р0=8·104Па, Θ =41 mm. 1/rc

  14. Some conclusions from spectral estimations • Maximum amount of three dimensional spectrum in the investigated region of waves numbers (Fig.5-6) corresponds to vortexes with character size 2rcθξmax , rc ≈ 0,3 mm (θ - length of realization, ξmax -wave number corresponds to maximum of spectra intercity) • This size is local dimension of investigatedturbulence. It can be used for obtaining estimations of energy of turbulent pulses near100 kDj/kg [experiment]. • As it will be shown this energy is well compared with character energy of physical-chemical processes determining the turbulence pulses 90,8 kDj/kg • which is the energy of dissociation of molecule of CF2Cl2.

  15. How to explaine the turbulence arising?Let’s write equation for vorticety Ω • The baroclinic term is conditioned by change the entropy - grad S = B

  16. Estimation of energy amount conditioned the flow instability For this aim the relation connected entropy changing, -ΔS, with a vortices in gas particle in homogeneous stationary flow. The amount of energy participated in thermal exchange, - q, for chemical processes is simply the same as T·S. Stationary equation of fluid dynamics will be:

  17. Deposition of chemical processes in change of entropy ΔS for great dissociation energy ΔH, little concentration lnζkand great gradζ k of products becomes very large: ∆Hk, Mk, ζk – creation energy, molecular weight and concentration of spesies; RΓ- gas constant This situation is just realized for 1. large Damkeler number ratio of reacting time to the time of flow ( so called almost frozen flow) and 2. for such polyatomic gas as CF2Cl2 with high energy of dissociation

  18. Use the model of rolling swirl (consider that vortex of characteristic size r is moving along the surface as if it is rolling along the surface without slide) ∆S = almost ∆q = 6 102 Here is Vm ≈ 250 - the middle flow velocity near the lateral surface m/s. [ found from estimation of full mass consumption of gas stream (≈ 3,4 kg/s) between the shock wave and lateral surface].

  19. The second step in turbulence investigation, which could give the answer for a question about a nature of phenomena, could be made, consider the dissipation function: Maximum of it corresponds to dimension of vortexes, on which main part of pulses energy in gas is dissipated • energy diffusion due to the dissipation take the place if their dimension about of order of amount or larger than Kolmogorov's micro dimension of turbulence: • εx- mass density of power of dissipation, • ν - the viscosity of gas

  20. Table 1.η - Kolmogorov's microdimension. ηd - Dimension of vortexes which is most intensely dissipated ξd - amount of wave number at which dissipation function have an absolute maximum, ε - mass density of power of dissipation.

  21. Some conclusions • estimation of power of dissipation of investigated turbulent flow (second row) • is close and almost higher than • estimation made for the case of release of all equilibrium energy of oscillatory degrees of freedom third row (accuracy of ε in frames of order of value). • Then its amount is greater than estimation for usual viscose stationary turbulence in three orders of value.

  22. So • appearance of turbulent flow in experiments is not conditioned by viscose dissipation of flow energy on the lateral surface • and is connected with excitation of more energy capable degrees of freedom than oscillatory one. • Such degrees of freedom can be, for instance, the chemical degrees or degrees of electron excitation.

  23. Fig.8. Scheme of flow around blunt body at the “transformation to reacting” regime. Dashed region is high entropy layer of gas flow. Reactions are only at stagnation point and products flow near the surface of body grad S is not coincides with vector of V, so rotV is not zero

  24. The gradient is conditioned by difference of entropy between its value on surface and on shock wave – ΔS, and by the depart of shock wave from surface - δ. Assume once more the following model of rolling swirl: let's consider that vortex of characteristic size is moving along the surface as if it is rolling along the surface without slide alike cylinder of radius – r(t). Then: =2V/r

  25. 1. Interaction of vortex with shock wave takes the place if r≥0,5δ; i.e. at 2. On the other hand for perturbation of shock wave it's necessary that V≤Vsound γ* , μ – effective ratio of specific thermo capacities of gas and its molecular weight in compression zone

  26. Two conditions take the place if: • this criteria is valid • for regime at which first (at shock wave Mach number arising) disturbances on shock wave were observed at flow of CF2Cl2 around the blunt body. • ( M∞=4, P∞=0,04 MPa ) • INDEED: Temperature and pressure in critical point will be: T=750 K, P=0,1 MPa

  27. There is experimental data of equilibrium chemical composition of dissociation of CF2Cl2 about at this temperature and pressure: • Trenwith A.B., Watson R.H. Thermal decomposition of chlorofluoromethanes.// J.Chem.Soc. (London). 1957. N 5. P.2368-2372 • Using it we can find T∆S=90,8 kDj/kg. • Then (ТS)=150,7 m/s. • Since sound speed in CF2Cl2 at this regimes 239 m/s (at γ* ≈1,1) criterion (14) take the place. • At rising of M∞ and P∞ they reach the values at which criterion (14) is not available (for instance for P∞=0,04 MPa it is M∞=6) • and center of disturbances creation is moving from region of lateral flow to region of bottom flow.

  28. Comparison with experiment • In experiments we observe that perturbations of shock wave front (and of the flow behind it) indeed • at P∞=0,04 MPa at M∞=4 is of find scale 0.3 mm • But at M∞=6 they are almost neglected.

  29. Summary • We can say from experimental studyand made estimations • that effect of so called instability of bow shock wave and flow behind it • Is conditioned by the baroclinic type of flow of reacting gas in compression zone between bow shock wave and bow surface • I am very thanks for collaboration first of all to A.Bedin, then to V.Derevyanko,V.Maslennikov, I.Basargin and M.Chistyakova. Thank you for attention

More Related