We have numerically investigated the interaction of a side jet

positioned on the small rocket, with the supersonic cross-flow.

An open source CFD tool, OpenFOAM is used to model the

complex flow of a jet-atmosphere interaction. The flow fields

are computed by the steady 3-dimensional Navier-Stokes solver

with k- ! SST turbulence model. Our solver is validated with

the experimental pressure data available on the rocket wall and

a systematic study is done by varying parameters like jet pressure

ratio. Aerothermodynamic coefficients for various flow

conditions are reported, and pitching moments and normal

forces are found to vary linearly with the jet pressure ratios.

Possible contamination of the on-board sensor located on the

rocket wall due to impinging plumes is also examined. This

study helps in designing effective missile control by selection

of the location of jet and pressure ratios

Bhagat, A

Dongari, Nishanth

Gijare, H

Research Archive of Indian Institute of Technology Hyderabad

Proceedings of the 6th International and 43rd National Conference on Fluid Mechanics and Fluid PowerDecember 15-17, 2016, MNNITA, Allahabad, U.P., IndiaFMFP2016PAPER NO.327Computational Analysis of Side Jet Interaction With a Super-sonicCross-flowApurva BhagatResearch ScholarIndian Institute of Technology,Hyderabad, 502285Email: me13m1026@iith.ac.inHarshal GijareResearch ScholarIndian Institute of Technology,Hyderabad, 502285Email: me13m1028@iith.ac.inNishanth DongariAssistant ProfessorIndian Institute of Technology,Hyderabad, 502285Email: nishanth@iith.ac.inAbstractWe have numerically investigated the interaction of a side jetpositioned on the small rocket, with the supersonic cross-flow.An open source CFD tool, OpenFOAM is used to model thecomplex flow of a jet-atmosphere interaction. The flow fieldsare computed by the steady 3-dimensional Navier-Stokes solverwith k- ω SST turbulence model. Our solver is validated withthe experimental pressure data available on the rocket wall anda systematic study is done by varying parameters like jet pres-sure ratio. Aerothermodynamic coefficients for various flowconditions are reported, and pitching moments and normalforces are found to vary linearly with the jet pressure ratios.Possible contamination of the on-board sensor located on therocket wall due to impinging plumes is also examined. Thisstudy helps in designing effective missile control by selectionof the location of jet and pressure ratios.Keywords—CFD; OpenFOAM; rhoCentralFoam; side jet con-trol; supersonic cross-flow.I. INTRODUCTIONSide jet control plays an important role in controllingmanoeuvrability and agility of the missile. Advantage of sidejet over the aerodynamic control surface is the precise controlof forces and quick response time. Side jet control is effectiveat both, the higher stagnation pressures conditions and therarefied conditions (altitude > 70 km). Figure 1 demonstratesthe complex flow structure developed due to the interactionof jet and supersonic cross-flow. Freestream hits the jet flow,which acts as a blunt obstacle and a bow shock is developed.Interaction of the jet flow with the boundary layer aroundmissile produces a separation shock upstream of the jet andrecirculation zone in the downstream region.Stahl et al. [1] carried out extensive wall pressure mea-surements on missile models to investigate the interference ofa side jet with the cross flow. Aswin et al. [2] further carriedout numerical analysis for the same case, and a very goodagreement between experimental and numerical results wasobtained. Similar problems are analysed using DES/LES bya few researchers.Major objective of this work is to develop an open sourceCFD solver for the application of complex side jet andfreestream interaction, and parametric study to design aneffective missile control. The solver can tackle multi-species,high temperature transport modelling (T 6000 K) and workon parallel computing architecture.Fig. 1: Schematic of flow structure of jet interaction with supersoniccross-flow.II. PROBLEM SETUPFigure 2 demonstrates the schematic of the missile modelwhich consists of a cone, cylindrical fuselage and flare. Outerdiameter (D) of the cylindrical fuselage is 40 mm and exitdiameter (d) of cylindrical nozzle is 4 mm. Nozzle is locatedat distance 4.3 D and radiance sensor at 4.6 D from the origin.To study the interference of side jet flow in supersoniccross-flow atmosphere, simulations are carried out step-wise:for external flow only (without jet flow), for jet flow only(without external flow), and with both external and jet flow.Validation study is carried out for the case, where freestreamMach number (Ma∞) = 3, Reynolds number ReD = 1.9 ×106 and pressure ratio of jet pressure to ambient pressure (PR)poj/p∞ = 200. This case is referred as the baseline case of thepaper and all other cases will be compared against this data.Simulations are conducted for PR = 55, 110 and 200. Behaviorof bow shock and recirculation zones are investigated byvarying parameters like PR. Lift coefficient, drag coefficient,pitching moment and the heat load values on the radiancesensor are reported.III. COMPUTATIONAL METHODOLOGYOpenFOAM (Open Field Operation and Manipulation) isan open source CFD software, which is parallel friendlyand handles both the structured and unstructured meshes forcomplex geometries. It is based on C++ library tools and1Fig. 2: Schematic of missile model under consideration (D = 40mm). Structured mesh on the missile wall (very fine near nozzle) isdemonstrated.a collection of various applications (created using these li-braries). The rhoCentralFoam is a density-based compressibleflow solver based on central upwind schemes of Kurganov andTadmor [4].The rhoCentralFoam solver has been validated by Green-shields et al. [5] for supersonic jet experiment by Ladenburget al. [6] and various standard compressible flow cases. k - ωSST turbulence model is implemented in simulations which isthe mix of k - ω and k -  models. The k - ω SST modelis merited for its good behavior in adverse pressure gradientsand separating flow.3D structured mesh is generated which consists of 43million hexahedral elements and fig.2 demonstrates mesh onthe wall of the missile. Average Y+ value on the fuselageof missile is 22.31 and cell-size near the wall is 0.01 mm(which can capture the desired Y+ value and viscous effects)and spline law is applied for gradual growth in mesh. Meshindependence study is carried out for 3 different meshes byvarying no. of cells on the missile wall, and solution is foundto be independent of mesh refinement.Boundary conditions are given as follows: 1. Missile -isothermal wall (no-slip for velocity, no-jump for temperature).2. Nozzle - pressure inlet. 3. Outlet - pressure far-field. 4. Inlet- freestream conditions. In current simulations, a caloricallyperfect ideal gas, air is used for both the freestream flow andnozzle flow. Each test case is simulated in parallel on 32 IntelHaswell cores on the HPC facility at IIT Hyderabad.IV. RESULTS AND DISCUSSIONFigure 3 demonstrates the comparison of Cpdif data (case- PR = 200) of rhoCentralFoam solver against experimentaldata [1], where Cp is the pressure coefficient defined as,Cp = (p − p∞)/p∞kMa2∞ and Cpdif = CpwithsideJet −CpwithoutsideJet Cp values are high at the upstream of nozzlelocation (x/D = 4.3) due to bow shock wave and reduces inthe downstream wake region. The rhoCentralFoam values arein decent agreement with the experimental data except in thedownstream of jet.Fig. 3: Distribution of pressure coefficient difference on cylindricalfuselage.Effect of change in PR - To investigate the effect of changein jet pressure ratio, simulations are carried out for PR = 55,110 and 200. Freestream Ma is 3 and exit Ma of jet flow is 1for this study. Fig.4 demonstrates density distribution and flowfeatures for the baseline case in XY plane (refer fig.2).Fig. 4: Density contours and flow features (PR = 200) in X-Y plane.Fig.5 provides pressure distribution on the missile surfacefor various PRs. It can be observed that bow shock shiftsupstream of the jet as PR increases and low pressure regionin the downstream region is more pronounced at higher PR.However, change in PR does not affect the flow on the coneas well as flare region of the missile. Higher pressure zone inthe upstream of jet and low pressure zone downstream to theinjection point are the major contributors of the moment ofmissile. Desired control forces can be obtained by choosing asuitable location of jet and an optimum range of PR.Fig. 5: Pressure distribution for various jet pressure ratios on the missilewall.Fig. 6: Wall heat flux distribution for various jet pressure ratios on themissile wall.Table 1 reports values of aero-thermodynamic coefficientsfor 3 jet pressure ratio cases. Here, Cm (moment coefficient)and CL (lift coefficient) linearly increase with increase in PR,while CD (drag coefficient) inversely varies with the same.This can be explained by the fact that higher PR case acts asbigger obstruction to free-stream flow, while lower PR casesleading to increase in drag in the forward direction.Heat load values on the radiance sensor (which is locatedat x/D = 4.6) are reported. Fig.6 demonstrates wall heat fluxdistribution on fuselage missile body. It is observed that heatload on the radiance sensor reduces with increase in PRwhereas heat flux values increase in bow shock region. This isbecause recompression/barrel shock is stronger for higher PRTABLE 1: Aero-thermodynamic coefficients and Heat load.PR = 200 PR = 110 PR = 55Cm 0.1352 0.1080 0.0760CD 0.6173 0.6210 0.6252CL 0.1478 0.1214 0.0876Wall Heat Load (W/m2) -723.33 628.61 1246.9leading to intense freestream-jet flow interaction and kineticenergy of freestream flow is converted to thermal energy. Jetflow expands rapidly in the downstream region where radiancesensor is located, and a low pressure and temperature zone isformed downstream of the nozzle. Temperature of jet flow forPR = 200 case drops below wall temperature and cooling ofwall takes place, which causes the negative heat load valueat the sensor. Therefore, thermal protection system for themissile surface should be designed considering the locationof the sensor and PR value.V. CONCLUSIONSNumerical simulations are carried out to demonstrate theinterference effects of a side jet in the supersonic cross-flow atmosphere. OpenFOAM open source CFD tool is usedto model a realistic 3-dimensional missile geometry in thecontinuum regime. Solver is validated with the experimentalsurface pressure data on the missile wall. Parametric studies arecarried out to study the effect of change in jet pressure ratio onthe aerodynamic coefficients and heat load on the sensor. Heatload on the fuselage increases with PR in the upstream regionof jet, as freestream and plume interaction is intense at highPR. Heat load reduces with increase in PR, in the downstreamregion of jet, due to the formation of a low pressure wake. Bycareful selection of the location of jet and pressure ratios, itis possible to design an effective missile control. An efficientthermal protection system can be designed based on predictedheat load values.REFERENCES[1] Stahl, Bernhard, Helmut Esch, and Ali Glhan. “Experimental investi-gation of side jet interaction with a supersonic cross flow.” AerospaceScience and Technology 12.4 (2008): 269-275.[2] Aswin, G., and Debasis Chakraborty. “Numerical simulation of trans-verse side jet interaction with supersonic free stream.” AerospaceScience and Technology 14.5 (2010): 295-301.[3] Gimelshein, S. F., D. A. Levin, and A. A. Alexeenko. “Modeling ofchemically reacting flows from a side jet at high altitudes.” Journal ofspacecraft and rockets 41.4 (2004): 582-591.[4] Kurganov, Alexander, Sebastian Noelle, and Guergana Petrova.”Semidiscrete central-upwind schemes for hyperbolic conservation lawsand Hamilton–Jacobi equations.” SIAM Journal on Scientific Comput-ing 23.3 (2001): 707-740.[5] Greenshields, Christopher J., et al.“Implementation of semidiscrete,nonstaggered central schemes in a colocated, polyhedral, finite volumeframework, for highspeed viscous flows.” International journal fornumerical methods in fluids 63.1 (2010): 1-21.[6] Ladenburg, R., C. C. Van Voorhis, and J. Winckler.“Interferometricstudies of faster than sound phenomena. Part II. Analysis of supersonicair jets.” Physical Review 76.5 (1949): 662.AcknowledgementsThe research is supported by DST:SERB/F/2684/2014-15 andMHRD fellowship. We are also thankful to Dr. V. K. Sarswatfor valuable suggestions.

Computational Analysis of Side Jet Interaction With a Super-sonic Cross-flow

http://raiith.iith.ac.in/3814/1/National%20Conference%20on%20Fluid%20Mechanics%20and%20Fluid%20Power.pdf

Computational Analysis of Side Jet Interaction With a Super-sonic Cross-flow

Abstract

Similar works

Full text

Available Versions

Research Archive of Indian Institute of Technology Hyderabad