Search NASA⌕ Search

Engineering topics

Balakumar, P.

Publications and source records attributed to Balakumar, P..

At least 19 records

Receptivity and Forced Response to Acoustic Disturbances in High-Speed Boundary Layers

Supersonic boundary-layer receptivity to freestream acoustic disturbances is investigated by solving the Navier–Stokes equations for Mach 3.5 flow over a 7 deg half-angle cone. The freestream disturbances are generated from a wavy wall placed at the nozzle wall. The freestream acoustic disturbances radiated by the wavy wall are obtained by solving the linearized Euler equations. The results show that no noticeable instability modes are generated when the acoustic disturbances impinge the cone obliquely. The results show that the perturbations generated inside the boundary layer by the acoustic disturbances are the response of the boundary layer to the external forcing. The amplitude of the forced disturbances inside the boundary layer are about 2.5 times larger than the incoming field for zero azimuthal wave number, and they are about 1.5 times for large azimuthal wave numbers.

Balakumar, P.↗

Receptivity and Forced Response to Acoustic Disturbances in High-Speed Boundary Layers

Supersonic boundary-layer receptivity to freestream acoustic disturbances is investigated by solving the Navier-Stokes equations for Mach 3.5 flow over a sharp flat plate and a 7-deg half-angle cone. The freestream disturbances are generated from a wavy wall placed at the nozzle wall. The freestream acoustic disturbances radiated by the wavy wall are obtained by solving the linearized Euler equations. The results for the flat plate show that instability modes are generated at all the incident angles ranging from zero to highly oblique. However, the receptivity coefficient decreases by about 20 times when the incident angle increases from zero to a highly oblique angle of 68 degrees. The results for the cone show that no instability modes are generated when the acoustic disturbances impinge the cone obliquely. The results show that the perturbations generated inside the boundary layer by the acoustic disturbances are the response of the boundary layer to the external forcing. The amplitude of the forced disturbances inside the boundary layer are about 2.5 times larger than the incoming field for zero azimuthal wavenumber and they are about 1.5 times for large azimuthal wavenumbers.

Balakumar, P.↗

Transition Prediction in Hypersonic Boundary Layers Using Receptivity and Freestream Spectra

Boundary-layer transition in hypersonic flows over a straight cone can be predicted using measured freestream spectra, receptivity, and threshold values for the wall pressure fluctuations at the transition onset points. Simulations are performed for hypersonic boundary-layer flows over a 7-degree half-angle straight cone with varying bluntness at a freestream Mach number of 10. The steady and the unsteady flow fields are obtained by solving the two-dimensional Navier-Stokes equations in axisymmetric coordinates using a 5th-order accurate weighted essentially non-oscillatory (WENO) scheme for space discretization and using a third-order total-variation-diminishing (TVD) Runge-Kutta scheme for time integration. The calculated N-factors at the transition onset location increase gradually with increasing unit Reynolds numbers for flow over a sharp cone and remain almost the same for flow over a blunt cone. The receptivity coefficient increases slightly with increasing unit Reynolds numbers. They are on the order of 4 for a sharp cone and are on the order of 1 for a blunt cone. The location of transition onset predicted from the simulation including the freestream spectrum, receptivity, and the linear and the weakly nonlinear evolutions yields a solution close to the measured onset location for the sharp cone. The simulations over-predict transition onset by about twenty percent for the blunt cone.

Balakumar, P.↗

On the Effect of Rigid Swept Surface Waves on Turbulent Drag

Passive turbulent drag reduction techniques are of interest as a cost effective means to improve air vehicle fuel consumption. In the past, rigid surface waves slanted at an angle from the streamwise direction were deemed ineffective to reduce skin friction drag due to the pressure drag that they generate. A recent analysis seeking similarities to the spanwise shear stress generated by spatial Stokes layers suggested that there may be a range of wavelength, amplitude, and orientation in which the wavy surface would reduce turbulent drag. The present work explores, by experiments and Direct Numerical Simulations (DNS), the effect of swept wavy surfaces on skin friction and pressure drag. Plates with shallow and deep wave patterns were rapid-prototyped and tested using a drag balance in the 7x11 inch Low-Speed Wind Tunnel at the NASA LaRC Research Center. The measured drag o set between the wavy plates and the reference at plate is found to be within the experimental repeatability limit. Oil vapor flow measurements indicate a mean spanwise flow over the deep waves. The turbulent flow in channels with at walls, swept wavy walls and spatial Stokes spanwise velocity forcing was simulated at a friction Reynolds number of two hundred. The time-averaged and dynamic turbulent flow characteristics of the three channel types are compared. The drag obtained for the channel with shallow waves is slightly larger than for the at channel, within the range of the experiments. In the case of the large waves, the simulation over predicts the drag. The shortcomings of the Stokes layer analogy model for the estimation of the spanwise shear stress and drag are discussed.

Denison, M.↗

DNS/LES Simulations of Separated Flows at High Reynolds Numbers

Direct numerical simulations (DNS) and large-eddy simulations (LES) simulations of flow through a periodic channel with a constriction are performed using the dynamic Smagorinsky model at two Reynolds numbers of 2800 and 10595. The LES equations are solved using higher order compact schemes. DNS are performed for the lower Reynolds number case using a fine grid and the data are used to validate the LES results obtained with a coarse and a medium size grid. LES simulations are also performed for the higher Reynolds number case using a coarse and a medium size grid. The results are compared with an existing reference data set. The DNS and LES results agreed well with the reference data. Reynolds stresses, sub-grid eddy viscosity, and the budgets for the turbulent kinetic energy are also presented. It is found that the turbulent fluctuations in the normal and spanwise directions have the same magnitude. The turbulent kinetic energy budget shows that the production peaks near the separation point region and the production to dissipation ratio is very high on the order of five in this region. It is also observed that the production is balanced by the advection, diffusion, and dissipation in the shear layer region. The dominant term is the turbulent diffusion that is about two times the molecular dissipation.

Balakumar, P.↗

Receptivity of Hypersonic Boundary Layers to Acoustic and Vortical Disturbances (Invited)

Boundary-layer receptivity to two-dimensional acoustic and vortical disturbances for hypersonic flows over two-dimensional and axi-symmetric geometries were numerically investigated. The role of bluntness, wall cooling, and pressure gradients on the receptivity and stability were analyzed and compared with the sharp nose cases. It was found that for flows over sharp nose geometries in adiabatic wall conditions the instability waves are generated in the leading-edge region and that the boundary layer is much more receptive to slow acoustic waves as compared to the fast waves. The computations confirmed the stabilizing effect of nose bluntness and the role of the entropy layer in the delay of boundary layer transition. The receptivity coefficients in flows over blunt bodies are orders of magnitude smaller than that for the sharp cone cases. Wall cooling stabilizes the first mode strongly and destabilizes the second mode. However, the receptivity coefficients are also much smaller compared to the adiabatic case. The adverse pressure gradients increased the unstable second mode regions.

Balakumar, P.↗

Effects of Forward- and Backward-Facing Steps on the Crossflow Receptivity and Stability in Supersonic Boundary Layers

The effects of forward- and backward-facing steps on the receptivity and stability of three-dimensional supersonic boundary layers over a swept wing with a blunt leading edge are numerically investigated for a freestream Mach number of 3 and a sweep angle of 30 degrees. The flow fields are obtained by solving the full Navier-Stokes equations. The evolution of instability waves generated by surface roughness is simulated with and without the forward- and backward-facing steps. The separation bubble lengths are about 5-10 step heights for the forward-facing step and are about 10 for the backward-facing step. The linear stability calculations show very strong instability in the separated region with a large frequency domain. The simulation results show that the presence of backward-facing steps decreases the amplitude of the stationary crossflow vortices with longer spanwise wavelengths by about fifty percent and the presence of forward-facing steps does not modify the amplitudes noticeably across the steps. The waves with the shorter wavelengths grow substantially downstream of the step in agreement with the linear stability prediction.

Balakumar, P.↗

Receptivity of Hypersonic Boundary Layers to Distributed Roughness and Acoustic Disturbances

Boundary-layer receptivity and stability of Mach 6 flows over smooth and rough seven-degree half-angle sharp-tipped cones are numerically investigated. The receptivity of the boundary layer to slow acoustic disturbances, fast acoustic disturbances, and vortical disturbances is considered. The effects of three-dimensional isolated roughness on the receptivity and stability are also simulated. The results for the smooth cone show that the instability waves are generated in the leading edge region and that the boundary layer is much more receptive to slow acoustic waves than to the fast acoustic waves. Vortical disturbances also generate unstable second modes, however the receptivity coefficients are smaller than that of the slow acoustic wave. Distributed roughness elements located near the nose region decreased the receptivity of the second mode generated by the slow acoustic wave by a small amount. Roughness elements distributed across the continuous spectrum increased the receptivity of the second mode generated by the slow and fast acoustic waves and the vorticity wave. The largest increase occurred for the vorticity wave. Roughness elements distributed across the synchronization point did not change the receptivity of the second modes generated by the acoustic waves. The receptivity of the second mode generated by the vorticity wave increased in this case, but the increase is lower than that occurred with the roughness elements located across the continuous spectrum. The simulations with an isolated roughness element showed that the second mode waves generated by the acoustic disturbances are not influenced by the small roughness element. Due to the interaction, a three-dimensional wave is generated. However, the amplitude is orders of magnitude smaller than the two-dimensional wave.

Balakumar, P.↗

Crossflow Instability on a Wedge-Cone at Mach 3.5

As a follow-on activity to the HyBoLT flight experiment, a six degree half angle wedge-cone model at zero angle of attack has been employed to experimentally and computationally study the boundary layer crossflow instability at Mach 3.5 under low disturbance freestream conditions. Computed meanflow and linear stability analysis results are presented along with corresponding experimental Pitot probe data. Using a model-mounted probe survey apparatus, data acquired to date show a well defined stationary crossflow vortex pattern on the flat wedge surface. This effort paves the way for additional detailed, calibrated flow field measurements of the crossflow instability, both stationary and traveling modes, and transition-to-turbulence under quiet flow conditions as a means of validating existing stability theory and providing a foundation for dynamic flight instrumentation development.

Beeler, George B.↗

Stability of Hypersonic Boundary Layers on a Cone at an Angle of Attack

The stability and receptivity of a three-dimensional hypersonic boundary layer over a 7deg half-angle straight cone at an angle of attack of 6deg is numerically investigated at a freestream Mach number of 6.0 and a Reynolds number of 10.4x10(exp 6)/m. The generation and evolution of stationary crossflow vortices are also investigated by performing simulations with three-dimensional roughness elements located on the surface of the cone. The flow fields with and without the roughness elements are obtained by solving the full Navier- Stokes equations in cylindrical coordinates using a fifth-order accurate weighted essentially non-oscillatory (WENO) scheme for spatial discretization and a third-order total-variation-diminishing (TVD) Runge-Kutta scheme for temporal integration. Stability computations produced azimuthal wavenumbers in the range of m approx. 20-50 for the most amplified traveling disturbances and in the range of m approx.30-70 for the stationary disturbances. The frequency of the unstable second-mode ranges from 400 kHz to 900 kHz along the windward ray. The N-Factor computations predicted transition would occur more forward on the sides of the cone as compared to the transition fronts near the windward and the leeward rays. The simulations also show the crossflow vortices originating from the nose region propagate towards the leeward ray. No perturbations were observed toward the windward half of the cone.

Balakumar, P.↗

Boundary Layer Receptivity due to Roughness and Freestream Sound for Supersonic Flows over Axisymmetric Cones

The receptivity of supersonic boundary layers due to the interaction of a two-dimensional acoustic wave with a three-dimensional roughness is numerically investigated over a 5-degree sharp tipped straight cone at a free stream Mach number of 3.5 and a high Reynolds number of 10(exp 6)/inch. Both the steady and unsteady solutions are obtained by solving the full Navier-Stokes equations using the fifth-order accurate weighted essentially non-oscillatory (WENO) scheme for spatial discretization and using third-order total-variation-diminishing (TVD) Runge-Kutta scheme for temporal integration. The simulations are performed with different roughness heights and roughness locations. The steady computations with the roughness reveal that the perturbations induced by the roughness remain very close to the wall and do not penetrate to the outer part of the boundary layer. The effects of the roughness persist for about 10 boundary layer thicknesses downstream of each element. The unsteady simulation shows the generation of unstable three-dimensional first mode disturbances within a short distance from the roughness. However, the amplitudes of the instability waves are very small near the neutral point, on the order of 10(exp -4) compared to the free stream disturbances. The receptivity coefficients due to the roughness elements located very close to the tip are one order of magnitude smaller than those due to the roughness elements located closer to the neutral point.

Balakumar, P.↗

Effects of Wall Cooling on Hypersonic Boundary Layer Receptivity Over a Cone

Effects of wall cooling on the receptivity process induced by the interaction of slow acoustic disturbances in the free-stream are numerically investigated for a boundary layer flow over a 5-degrees straight cone. The free-stream Mach number is 6.0 and the Reynolds number is 7.8x10(exp 6)/ft. Both the steady and unsteady solutions are obtained by solving the full Navier-Stokes equations using 5th-order accurate weighted essentially non-oscillatory (WENO) scheme for space discretization and using 3rd-order total variation diminishing (T VD) Runge-K utta scheme for time integration. Computations are performed for a cone with nose radius of 0.001 inch for adiabatic wall temperature (T(sub aw)), 0.75*T(sub aw), 0.5*T(sub aw), 0.40*T(sub aw), 0.30*T(sub aw), and 0.20*T(sub aw). Once the mean flow field is computed, disturbances are introduced at the upstream end of the computational domain. Generation of instability waves from leading edge region and receptivity of boundary layer to slow acoustic waves are investigated. Computations showed that wall cooling has strong stabilization effect on the first mode disturbances as was observed in the experiments. T ransition location moved to upstream when wall cooling was applied It is also found that the boundary layer is much more receptive to fast acoustic wave (by almost a factor of 50). When simulations performed using the same forcing frequency growth of the second mode disturbances are delayed with wall cooling and they attained values two times higher than that of adiabatic case. In 0.20*T(sub aw) case the transition Reynolds number is doubled compared to adiabatic conditions. The receptivity coefficient for adiabatic wall case (804 R) is 1.5225 and for highly cooled cones (241, and 161 R); they are in the order of 10(exp -3).

Kara, K.↗

Receptivity of Supersonic Boundary Layers Due To Acoustic Disturbances Over Blunt Cones

Receptivity and stability of supersonic boundary layers over a 5-degree straight cone with a blunt tip are numerically investigated at a free stream Mach number of 3.5 and at a high Reynolds number of 106/inch. Both the steady and unsteady solutions are obtained by solving the full Navier-Stokes equations using the 5th-order accurate weighted essentially non-oscillatory (WENO) scheme for space discretization and using third-order total-variation-diminishing (TVD) Runge-Kutta scheme for time integration. The linear stability results showed that bluntness has less stabilizing effects on the stability of boundary layers over cones than on flat plates and wedges. The unsteady simulations of the interaction of plane threedimensional acoustic waves with the cone showed that the modulation of wavelength and the generation of instability waves first occurred near the leading edge in the plane where the constant acoustic phase lines are perpendicular to the cone axis. Further downstream, this instability region spreads in the azimuthal direction from this plane.

Balakumar, P.↗

Effects of Nose Bluntness on Stability of Hypersonic Boundary Layers over Blunt Cone

Receptivity and stability of hypersonic boundary layers are numerically investigated for boundary layer flows over a 5-degree straight cone at a free-stream Mach number of 6.0. To compute the shock and the interaction of shock with the instability waves, we solve the Navier-Stokes equations in axisymmetric coordinates. The governing equations are solved using the 5th-order accurate weighted essentially non-oscillatory (WENO) scheme for space discretization and using third-order total-variation-diminishing (TVD) Runge-Kutta scheme for time integration. After the mean flow field is computed, disturbances are introduced at the upstream end of the computational domain. Generation of instability waves from leading edge region and receptivity of boundary layer to slow acoustic waves are investigated. Computations are performed for a cone with nose radii of 0.001, 0.05 and 0.10 inches that give Reynolds numbers based on the nose radii ranging from 650 to 130,000. The linear stability results showed that the bluntness has a strong stabilizing effect on the stability of axisymmetric boundary layers. The transition Reynolds number for a cone with the nose Reynolds number of 65,000 is increased by a factor of 1.82 compared to that for a sharp cone. The receptivity coefficient for a sharp cone is about 4.23 and it is very small, approx.10(exp -3), for large bluntness.

Kara, K.↗

Computations of Flow Over the Hump Model Using Higher-Order Method With Turbulence Modeling

The flow over the two-dimensional hump model is computed by solving the RANS equations with kappa-omega (SST) model. The governing equations, the flow equations and the turbulent equations, are solved using the 5th order accurate weighted essentially non-oscillatory (WENO) scheme for space discretization and using explicit third order total-variation-diminishing (TVD) Runge-Kutta scheme for time integration. The WENO and the TVD methods and the formulas are explained in [1] and the application of ENO method to N-S equations is given in [2]. The solution method implemented in this computation is described in detail in [3].

Balakumar, P.↗

Receptivity of Hypersonic Boundary Layers Due to Acoustic Disturbances over Blunt Cone

The transition process induced by the interaction of acoustic disturbances in the free-stream with boundary layers over a 5-degree straight cone and a wedge with blunt tips is numerically investigated at a free-stream Mach number of 6.0. To compute the shock and the interaction of shock with the instability waves the Navier-Stokes equations are solved in axisymmetric coordinates. The governing equations are solved using the 5th -order accurate weighted essentially non-oscillatory (WENO) scheme for space discretization and using third-order total-variation-diminishing (TVD) Runge-Kutta scheme for time integration. After the mean flow field is computed, acoustic disturbances are introduced at the outer boundary of the computational domain and unsteady simulations are performed. Generation and evolution of instability waves and the receptivity of boundary layer to slow and fast acoustic waves are investigated. The mean flow data are compared with the experimental results. The results show that the instability waves are generated near the leading edge and the non-parallel effects are stronger near the nose region for the flow over the cone than that over a wedge. It is also found that the boundary layer is much more receptive to slow acoustic wave (by almost a factor of 67) as compared to the fast wave.

Kara, K.↗

Control of Supersonic Boundary Layers Using Steady Suction

Control of supersonic boundary layers using steady suction through a series of very small two-dimensional strips is numerically investigated at a free stream Mach number of 1.8. The mean flow induced by rows of suction holes is also computed. Both the steady and unsteady solutions are obtained by solving the full Navier-Stokes equations using the 5th-order accurate weighted essentially non-oscillatory (WENO) scheme for space discretization and using third-order total-variationdiminishing (TVD) Runge-Kutta scheme for time integration. Computations for the two-dimensional cases are performed at suction coefficients 0.001 and 0.002 to investigate the stabilizing effects of suction. The simulation showed that a series of shock waves are generated at the slots. The stability results showed that the total amplification is reduced up to the end of the computational domain. However, the growth rates become larger at downstream distances away from the suction region. The computations for the suction holes showed the generation of Mach waves from each hole and the formation of longitudinal vortices.

Balakumar, P.↗