A commercial CFD code based on a pressure-based algorithm (SIMPLEC) has been used to compute high Mach number flow over a cavity, which has application in supersonic combustors. The difficulties encountered in obtaining converged solution for compressible flow by the pressure-based method are discussed. The predicted solution is compared with the available experimental results. Though some quantitative differences exist between them, the qualitative agreement is generally

good

Khan, AA

Shembharkar, TR

English

National Aerospace Laboratories Institutional Repository

31' Natioml Conference on Fluid MecI&n&s &Fluid Power 16-18 Dec 2004, Jadavpur UniversiLjr, E- India Computation of High S@ Cold Cavity Flow with a Commercial CFD Code Ashfaque A. Khan and T. R. Shembharkar Propulsion Division, N & d  Aerospace Laboratories, Bangalore. ABSTRACT A commercial CFD code based on a pressure-based algorithm (SIMPLEC) has been used to compute'high Mach number flow over a cavity, which has application in supersonic combustors; difficulties encountered in obtaining converged solution for compressible flow by the pressure-based method are discussed. The predicted solution is compared with the available experimental results. Though some quantitative &fferences, exisr between them, the qualitative agreement is generally good. Keywords: CFD, cavity, supersonic flow I., INTRODUCTION Supersonic air-breathing engines are vital elements of future high-speed transportation vehicles. These engines invariably include supersonic combustors. In a scramjet, flame holding with minimum total pressure loss is a key to obtain good performance. In recent years, cavities have been proposed as a new concept for flame holding and stabilization in supersonic combustors. Therefore, it is crucial to understand the basic flow features of cavities in high-speed flows. There have been some excellent review papers [Ben-Yakar et all on the aspects of high-speed flow behaviour over the cavity. In general, the boundary layer ahead of the cavity separates at the leading edge and forms a free shear layer across the cavity. Inside the cavity, flow recirculation takes place. The shear layer reattaches at some other point downstream. The reattachment point depends upon the geometry of the cavity and the external flow conditions. Depending upon the reattachment point, the cavities are classified as "open" or "closed". In open cavity, the reattachment takes place at the back face of the cavity and in closed cavity it takes place on the lower wall. The open cavities have aspect ratio less than about 7-10 while the closed cavity have higher aspect ratios. The leading edge expansion wave, flow separation, flow recirculation, reattachment of free @hear layer and miling edge shock contribute to the pressure: loss,in the caGsy which neenis to be kept low for acceptable performance. '?'he present investigation is dated to flow analysis thrcruj~h a slant cavity (the back : 'face is inclined), which. can be used in a scramjet to get &k flatne holding and recirculdon properties. T b  analysis has bmn c&ed out wing a comrpercial CljT) &de called CFD-ACE.+. [cFDaC], which is a general puqbc  CFD. code b& in a pf~ssw-based algarithm, Besides understanding the -cavity flow, them is a second there appears to be a sart of m-bas& or density-bas4 In purely i~compressiblk flow - - aomain, density being constant and not a function of pressure, there ie no separate ,equation for p d r e .  The Poisson-type plressure ar p m u s e  cart-n equarion is -@.eiived from continuity and rnomnmm equations, Such mtft& are broadly called < ptssure-based methods. In compmsibb -Elow domain, the continuity equation ved mparately for densiry. ity and p w u r e  e.g. ideal sda are bwdly c&ld density-bawd nyrhods. Both these methods b v e  Wit plus and minus pain&. The pressure-bd methods are very robust in kstzampressible and low Mach number flaws whiIe the dmsity-based methods are well e flows. The preksllre-based methods seem to get into &wbb in flow r&me while the density-based methods face severe nce i ~ i  incompressible flow-regime. It has k e n  a very active area the working flow regimes of.thae two approaches by &vising In the p m t  Cm) scenario; when conmercial CFD codes we &masingly being used for solving p ~ ~ l ~ i n d u s t r i a l  problem, it may be noted that ; &St of these codes ape based on one or the other method mentionai above. The code method. Our goal was W algoti-thm for high speed ion has been performed for a ion at t h ~  inlet Mach number of 2.46 and the predicted ith the experimentat wu1t;s [Saumy et all. n and ttre geometry of the slant cavity. A 2-D nerated. The grid distribution was prescribed s and flow gradients. The 2-I3 steady state RMG k-E turbulence model. .The present have b n  d m  with a pressure-bd method (SWLEC algorithm) on me approach. 4 1  the enter. The mnvwive tams e diffusion terms by central diffarenee is employed, whereid the equation sets for each variable are solved sequentially and repeatedly by conjugate squared with pre-conditioning linear solver until a converged solution The relevant velocity, pressure and density .at the air inlet plane were taken Reference [Samimy et all. The inlet profiles were taken as uni walls of the cavity were specified no-slip wall condition. The however, given symmetry condition, which effectively means impervious wall. The exit boundary conditions were satisfied by corresponding variables from inside. 3. RESULTS AND DISCUSSION During the computations, it was realized that the high-speed flow computation code was fraught with difficulties. Unlike the incompressible and low Ma compressible flow regimes for which the code has shown extremely robus [Muralidhara et all, the high Mach number flow regime makes it vul could not be obtained. However, the equivalent inlet conditions in term and density at inlet stabilizes the convergence. Further, the no-slip wall thk top boundary also made the solution prone tp divergence. The differencing for convective terms also posed problems in obtaining converg In short, it required extensive experimentation with different options to meaningful solution for the problem. This experience supports about the pressure-based algorithms that they are susceptible to non-conv divergence of solution in high Mach number flow regime. The final s obtained with prescribed velocity and density conditions at the inlet compared with the experimental results [Samimy et all. The predicted solution, in spite of slight difference in top wall boundary in qualitative agreement with the experimental observations. The I expansion wave, low Mach number flow region inside the cavity, above the cavity, reattachment of shear layer on the slant wall and an downstream (recompression and shock) are clearly visible in Figure-2. the recirculation zone can be seen in the Figure-3. The position of recirculation compares favourably with the experimental ob A comparison of predicted and measured location of M=l line (sonic line (dividing line, some measure of recirculating vortex) in the flow do in Figure-4. The predicted sonic line is slightly below the experimental v cavity but matches well along the slant wall. Similarly the dividing line the experimental line but the difference is considerably reduced near Figure-5 shows the non-dimensional velocity profiles at different sueamwise locations in the flow field. These profiles also show good qualitative agreement with the experimental results shown here. Figure-6 compares the wall static pressure in the cavity. The present predictions are the static pressure at the bottom of the cavity. They seem to be slightly lower than the experimental values inside the cavity and slightly higher along the slant wall. It may be noted here that the exact location of the measurement points for static pressure is not quite clear in the reference. 4. CONCLUSIONS The commercial CFD code, namely CFD-ACE+, based on a pressure-based algorithm has been used to compute flow over a cavity with free stream Mach number of 2.46. It is concluded that the pressure- based algorithms do face problems of convergence in high Mach number flows. It is not very straightforward to use this algorithm for compressible flows. It needs quite some adjustment 'with various solution controlling options and boundary conditions to get meaningful solutions. The predicted solution is compared with the available experimental results. Though some quantitative differences exist between them, the qualitative agreement is generally good. The expansion wave at the leading edge, flow separation, flow recirculation inside the cavity, reattachment of free shear layer at the slant wall and the trailing edge shock have been qualitatively well predicted. 5. ACKNOWLEDGEMENT The authors thank the CSIR Centre for Mathematical modelling and Computer Simulation (C-MMACS), Bangalore for providing the computing faciiities. REFERENCES Adela Ben-Yakar and Ronald K. Hanson, 'Cavity flame-holders for ignition and flame stabilization in scramjets: An overview', Journal of Propulsion and Power. 17(4), 869- 877.2001. CFD-ACE+ code from CFD Research Corporation, USA (www.cfdrc.com) - available to the authors courtesy Centre for Mathematical Modelling and Computer Simulation (C-MMACS), Bangalore. Samimy, M., Patrie, H. L. and Addy, A. L., "A study of compressible turbulent reattaching free shear layers", AIAA Journal, 24(2), 261 (1986) H.S.Muralidhara, B.B.Chinmaya Kumar and T.R.Shembharkar, "CFD analysis of gas turbine combustor: A CAD to Mesh to CFD approach", 61h National Conference on Air Breathing Engines and Aerospace Propulsion (NCABE), Paper B1.l, 10-1 1 January 2003. Figure 1. Computational Domain for the Slant Cavity. Much - - - 2.891 Mach number umtnur (m 1 ) Figure 2. Mach Number Contours (Prbent Predictions). u ~ v r * r w o ~ l )  ' F Present Predictions. / Samimy et al. Figure 3. Recirculating  lo& in the Cavity (Present Predictions and s&nimy et al). Figure 4. Comparison between Present Predictions and Experimental Results for Sonic and U=O line. UIOLine 30. . + Samimy et al - prediction r a:, .;:10 d& > 0 1 0 60 100 1BO ' m *(mm) . velocity' profile at different sections (rc measured from end of step; UN1 scale for xc=IOmm profile) ~. Present Predictions. Samirny et al. Figure 5. Velocity Profiles at different Sections (hesent Predictions and Samirny et al). "11 static Pr-um + Samimy et a1 - prediction 3.5 3 2.5 r 2 % 1.5 1 0.5 0 0 0.05 0.1 0.1 5 0.2 0.25 x (m) Figure 6. Comparison between Present Predictions and Experimental Results for Wall Static Pressure. 

Computation of High Speed Cold Cavity Flow with a

Commercial CFD Code

http://nal-ir.nal.res.in/11597/1/computation.pdf

Computation of High Speed Cold Cavity Flow with a Commercial CFD Code

Abstract

Similar works

Full text

Available Versions

National Aerospace Laboratories Institutional Repository