Engineering Papers⌕ Search

Engineering topics

Prabhu, Ramadas K.

Publications and source records attributed to Prabhu, Ramadas K..

28 records · Page 2

An approximate Riemann solver for thermal and chemical nonequilibrium flows

Among the many methods available for the determination of inviscid fluxes across a surface of discontinuity, the flux-difference-splitting technique that employs Roe-averaged variables has been used extensively by the CFD community because of its simplicity and its ability to capture shocks exactly. This method, originally developed for perfect gas flows, has since been extended to equilibrium as well as nonequilibrium flows. Determination of the Roe-averaged variables for the case of a perfect gas flow is a simple task; however, for thermal and chemical nonequilibrium flows, some of the variables are not uniquely defined. Methods available in the literature to determine these variables seem to lack sound bases. The present paper describes a simple, yet accurate, method to determine all the variables for nonequilibrium flows in the Roe-average state. The basis for this method is the requirement that the Roe-averaged variables form a consistent set of thermodynamic variables. The present method satisfies the requirement that the square of the speed of sound be positive.

Prabhu, Ramadas K.↗

An implementation of a chemical and thermal nonequilibrium flow solver on unstructured meshes and application to blunt bodies

This paper presents a nonequilibrium flow solver, implementation of the algorithm on unstructured meshes, and application to hypersonic flow past blunt bodies. Air is modeled as a mixture of five chemical species, namely O2, N2, O, NO, and N, having two temperatures namely translational and vibrational. The solution algorithm is a cell centered, point implicit upwind scheme that employs Roe's flux difference splitting technique. Implementation of this algorithm on unstructured meshes is described. The computer code is applied to solve Mach 15 flow with and without a Type IV shock interference on a cylindrical body of 2.5mm radius representing a cowl lip. Adaptively generated meshes are employed, and the meshes are refined several times until the solution exhibits detailed flow features and surface pressure and heat flux distributions. Effects of a catalytic wall on surface heat flux distribution are studied. For the Mach 15 Type IV shock interference flow, present results showed a peak heat flux of 544 MW/m2 for a fully catalytic wall and 431 MW/m(exp 2) for a noncatalytic wall. Some of the results are compared with available computational data.

Prabhu, Ramadas K.↗

Application of a two-equation turbulence model for high speed compressible flows using unstructured grids

Application of a two-equation compressible turbulence model for practical hypersonic flows 18 presented. The solution algorithm is based on solving all of the governing equations simultaneously. Application of the solution procedure to several test cases for compressible flows show good agreement with theoretical predictions and/or other computational results. The solution procedure is employed to investigate the effects of turbulence in Type III and Type IV shock-shock interactions in hypersonic flows in association with adaptive unstructured grids. Computational results for these cases are compared with available experimental data.

Vemaganti, Gururaja R.↗

Computational method to predict thermodynamic, transport, and flow properties for the modified Langley 8-foot high-temperature tunnel

The Langley 8 foot high temperature tunnel (8 ft HTT) is used to test components of hypersonic vehicles for aerothermal loads definition and structural component verification. The test medium of the 8 ft HTT is obtained by burning a mixture of methane and air under high pressure; the combustion products are expanded through an axisymmetric conical contoured nozzle to simulate atmospheric flight at Mach 7. This facility was modified to raise the oxygen content of the test medium to match that of air and to include Mach 4 and Mach 5 capabilities. These modifications will facilitate the testing of hypersonic air breathing propulsion systems for a wide range of flight conditions. A computational method to predict the thermodynamic, transport, and flow properties of the equilibrium chemically reacting oxygen enriched methane-air combustion products was implemented in a computer code. This code calculates the fuel, air, and oxygen mass flow rates and test section flow properties for Mach 7, 5, and 4 nozzle configurations for given combustor and mixer conditions. Salient features of the 8 ft HTT are described, and some of the predicted tunnel operational characteristics are presented in the carpet plots to assist users in preparing test plans.

Venkateswaran, S.↗

Computational studies of a fluid spike as a leading edge protection device for shock-shock interference heating

The effectiveness of a fluid spike as a device to protect leading edges of hypersonic atmospheric flight vehicles from high aerothermal loads produced by complex shock-shock interference is studied. The two-dimensional Navier-Stokes equations are solved using an unstructured cell-centered, fully implicit, flux-difference split algorithm. Adaptively generated unstructured meshes are employed. A type IV shock-shock interference for Mach 8 flow on a cylindrical leading edge with and without a small contraflow supersonic jet (fluid spike) placed at two different locations on the body is solved. A typical flow past a blunt body with a type IV shock-shock interference produces very high pressures and heat fluxes on the leading edge. Present results indicate that a fluid spike displaces the bow shock further in front of the body and modifies the shock-shock interference pattern. This leads to reduced peak pressures and heat fluxes on the body.

Prabhu, Ramadas K.↗

Shock interference studies on a circular cylinder at Mach 16

This paper presents the results of a computational study on the shock interference problems on a cylindrical body typical of an engine inlet cowl leading edge at a nominal Mach number of 16. The two-dimensional Navier-Stokes equations are solved assuming the flow to be in chemical and thermal equilibrium and using a finite element method. The algorithm employs a cell-centered fully implicit upwind scheme. Adaptively generated unstructured meshes of triangles and quadrilaterals are employed. Under certain conditions the finite element code resulted in oscillatory solutions for shock interference at Mach 16. Some of the causes of the unsteady behavior are identified, and to the extent possible, such situations are avoided in the present application. Two shock interference conditions involving a Type IV (supersonic jet) and a Type III (attaching shear layer) are solved. The results are compared with available experimental data and reasonable agreement is observed. A semi-empirical method is also used to estimate the maximum surface heat flux and static pressure.

Prabhu, Ramadas K.↗

A Navier-Stokes solver for high speed equilibrium flows and application to blunt bodies

This paper presents a finite element method for the solution of Navier-Stokes equations with the assumption of thermodynamic and chemical equilibrium. The method employs an upwind finite element technique with an implicit time-marching scheme for the solution, and uses an adaptively generated unstructured triangular mesh with several layers of quadrilateral elements near solid walls. The complexity associated with the assumption that the flow is in equilibrium is treated consistently, and the inviscid flux Jacobian matrices are derived. Several problems involving inviscid and viscous hypersonic flow past blunt are solved. Results are compared with other numerical results and experimental data, and generally good agreement is observed.

Prabhu, Ramadas K.↗

Euler and Navier-Stokes solutions for hypersonic flows

An upwind finite-element technique that uses cell-centered quantities and implicit and/or explicit time marching has been developed for computing hypersonic laminar viscous flows using adaptive unstructured grids in two and three dimensions. A perfect gas model as well as an equilibrium air model is implemented for solving high-speed flows. A first-order basic scheme and a higher-order flux-corrected transport (FCT) scheme have been implemented. This technique has been used to predict 'Type III and IV' shock interactions on a cylinder in two dimensions and a swept cylinder in three dimensions, with a view to determine the pressure and heating rate augmentation caused by an impinging shock on the leading edge of a cowl lip of an engine inlet. The predictions of wall pressure and heating rates compare very well with experimental data. The flow features are very distinctly captured with a sequence of adaptively-generated grids. Three-dimensional corner flow, typically encountered in engine inlets due to compression of the flow by ramps in the walls, is also modeled. This procedure is the first step in developing an integrated fluid, thermal, structural analysis capability for hypersonic flight vehicles like the National Aero-Space Plane.

Thareja, Rajiv R.↗

Finite-rate water condensation in combustion-heated wind tunnels

A quasi-one-dimensional method for computing finite rate nucleation and droplet growth of water in a supersonic expansion of combustion products is presented. Sample computations are included for the Langley 8 foot High Temperature Tunnel, but the method can also be applied to other combustion heated wind tunnels. The sample results indicate that the free stream static pressure can be in the range of 25 to 60 percent greater than that computed for isentropic nozzle flow without water condensation. The method provides a tool for examining the effects of water condensation on static state properties and velocity of the supersonic stream in combustion heated wind tunnels.

Erickson, Wayne D.↗

A rapid method for the computation of equilibrium chemical composition of air to 15000 K

A rapid computational method has been developed to determine the chemical composition of equilibrium air to 15000 K. Eleven chemically reacting species, i.e., O2, N2, O, NO, N, NO+, e-, N+, O+, Ar, and Ar+ are included. The method involves combining algebraically seven nonlinear equilibrium equations and four linear elemental mass balance and charge neutrality equations. Computational speeds for determining the equilibrium chemical composition are significantly faster than the often used free energy minimization procedure. Data are also included from which the thermodynamic properties of air can be computed. A listing of the computer program together with a set of sample results are included.

Prabhu, Ramadas K.↗