Engineering topics

Hassan, H. A.

Publications and source records attributed to Hassan, H. A..

At least 19 records

NASA NTRS2014

NASA Trapezoidal Wing Simulation Using Stress-w and One- and Two-Equation Turbulence Models

The Wilcox 2006 stress-omega model (also referred to as WilcoxRSM-w2006) has been implemented in the NASA Langley code CFL3D and used to study a variety of 2-D and 3-D configurations. It predicted a variety of basic cases reasonably well, including secondary flow in a supersonic rectangular duct. One- and two-equation turbulence models that employ the Boussinesq constitutive relation were unable to predict this secondary flow accurately because it is driven by normal turbulent stress differences. For the NASA trapezoidal wing at high angles of attack, the WilcoxRSM-w2006 model predicted lower maximum lift than experiment, similar to results of a two-equation model.

Rodio, J. J.↗

NASA NTRS2007

Role of Turbulent Prandtl Number on Heat Flux at Hypersonic Mach Numbers

A new turbulence model suited for calculating the turbulent Prandtl number as part of the solution is presented. The model is based on a set of two equations: one governing the variance of the enthalpy and the other governing its dissipation rate. These equations were derived from the exact energy equation and thus take into consideration compressibility and dissipation terms. The model is used to study two cases involving shock wave/boundary layer interaction at Mach 9.22 and Mach 5.0. In general, heat transfer prediction showed great improvement over traditional turbulence models where the turbulent Prandtl number is assumed constant. It is concluded that using a model that calculates the turbulent Prandtl number as part of the solution is the key to bridging the gap between theory and experiment for flows dominated by shock wave/boundary layer interactions.

Xiao, X.↗

NASA NTRS2006

Modeling Scramjet Flows with Variable Turbulent Prandtl and Schmidt Numbers

A complete turbulence model, where the turbulent Prandtl and Schmidt numbers are calculated as part of the solution and where averages involving chemical source terms are modeled, is presented. The ability of avoiding the use of assumed or evolution Probability Distribution Functions (PDF's) results in a highly efficient algorithm for reacting flows. The predictions of the model are compared with two sets of experiments involving supersonic mixing and one involving supersonic combustion. The results demonstrate the need for consideration of turbulence/chemistry interactions in supersonic combustion. In general, good agreement with experiment is indicated.

Xiao, X.↗

NASA NTRS2005

Role of Turbulent Prandtl Number on Heat Flux at Hypersonic Mach Numbers

Gaffney, R. L., Jr.↗

NASA NTRS2004

Modeling Turbulent Combustion for Variable Prandtl and Schmidt Number

This report consists of two abstracts submitted for possible presentation at the AIAA Aerospace Science Meeting to be held in January 2005. Since the submittal of these abstracts we are continuing refinement of the model coefficients derived for the case of a variable Turbulent Prandtl number. The test cases being investigated are a Mach 9.2 flow over a degree ramp and a Mach 8.2 3-D calculation of crossing shocks. We have developed an axisymmetric code for treating axisymmetric flows. In addition the variable Schmidt number formulation was incorporated in the code and we are in the process of determining the model constants.

Hassan, H. A.↗

NASA NTRS2004

Role of Turbulent Prandtl Number on Heat Flux at Hypersonic Mach Number

Present simulation of turbulent flows involving shock wave/boundary layer interaction invariably overestimates heat flux by almost a factor of two. One possible reason for such a performance is a result of the fact that the turbulence models employed make use of Morkovin's hypothesis. This hypothesis is valid for non-hypersonic Mach numbers and moderate rates of heat transfer. At hypersonic Mach numbers, high rates of heat transfer exist in regions where shock wave/boundary layer interactions are important. As a result, one should not expect traditional turbulence models to yield accurate results. The goal of this investigation is to explore the role of a variable Prandtl number formulation in predicting heat flux in flows dominated by strong shock wave/boundary layer interactions. The intended applications involve external flows in the absence of combustion such as those encountered in supersonic inlets. This can be achieved by adding equations for the temperature variance and its dissipation rate. Such equations can be derived from the exact Navier-Stokes equations. Traditionally, modeled equations are based on the low speed energy equation where the pressure gradient term and the term responsible for energy dissipation are ignored. It is clear that such assumptions are not valid for hypersonic flows. The approach used here is based on the procedure used in deriving the k-zeta model, in which the exact equations that governed k, the variance of velocity, and zeta, the variance of vorticity, were derived and modeled. For the variable turbulent Prandtl number, the exact equations that govern the temperature variance and its dissipation rate are derived and modeled term by term. The resulting set of equations are free of damping and wall functions and are coordinate-system independent. Moreover, modeled correlations are tensorially consistent and invariant under Galilean transformation. The final set of equations will be given in the paper.

Xiao, X.↗

NASA NTRS2004

A Variable Turbulent Schmidt Number Formulation for Scramjet Application

In high speed engines, thorough turbulent mixing of fuel and air is required to obtain high performance and high efficiency. Thus, the ability to predict turbulent mixing is crucial in obtaining accurate numerical simulation of an engine and its performance. Current state of the art in CFD simulation is to assume both turbulent Prandtl number and Schmidt numbers to be constants. However, since the mixing of fuel and air is inversely proportional to the Schmidt number, a value of 0.45 for the Schmidt number will produce twice as much diffusion as that with a value of 0.9. Because of this, current CFD tools and models have not been able to provide the needed guidance required for the efficient design of a scramjet engine. The goal of this investigation is to develop the framework needed to calculate turbulent Prandtl and Schmidt numbers as part of the solution. This requires four additional equations: two for the temperature variance and its dissipation rate and two for the concentration variance and its dissipation rate. In the current investigation emphasis will be placed on studying mixing without reactions. For such flows, variable Prandtl number does not play a major role in determining the flow. This, however, will have to be addressed when combustion is present. The approach to be used is similar to that used to develop the k-zeta model. In this approach, relevant equations are derived from the exact Navier-Stokes equations and each individual correlation is modeled. This ensures that relevant physics is incorporated into the model equations. This task has been accomplished. The final set of equations have no wall or damping functions. Moreover, they are tensorially consistent and Galilean invariant. The derivation of the model equations is rather lengthy and thus will not be incorporated into this abstract, but will be included in the final paper. As a preliminary to formulating the proposed model, the original k-zeta model with constant turbulent Prandtl and Schmidt numbers is used to model the supersonic coaxial jet mixing experiments involving He, O2 and air.

Xiao, X.↗

NASA NTRS2000

Theory and Experiment of Multielement Airfoils: A Comparison

A detailed comparison of computed and measured pressure distributions, velocity profiles, transition onset, and Reynolds shear stresses for multi-element airfoils is presented. It is shown that the transitional k-zeta model, which is implemented into CFL3D, does a good job of predicting pressure distributions, transition onset, and velocity profiles with the exception of velocities in the slat wake region. Considering the fact that the hot wire used was not fine enough to resolve Reynolds stresses in the boundary layer, comparisons of turbulence stresses varied from good to fair. It is suggested that the effects of unsteadiness be thoroughly evaluated before more complicated transition/turbulence models are used. Further, it is concluded that the present work presents a viable and economical method for calculating laminar/transitional/turbuient flows over complex shapes without user interface.

Czerwiec, Ryan↗

NASA NTRS1999

Study of the Transition Flow Regime using Monte Carlo Methods

This NASA Cooperative Agreement presents a study of the Transition Flow Regime Using Monte Carlo Methods. The topics included in this final report are: 1) New Direct Simulation Monte Carlo (DSMC) procedures; 2) The DS3W and DS2A Programs; 3) Papers presented; 4) Miscellaneous Applications and Program Modifications; 5) Solution of Transitional Wake Flows at Mach 10; and 6) Turbulence Modeling of Shock-Dominated Fows with a k-Enstrophy Formulation.

Hassan, H. A.↗

NASA NTRS1997

Transition and Turbulence Modeling for Blunt-Body Wake Flows

Aerobraking has been proposed as an efficient means of decelerating spacecraft for planetary missions. Most current aerobrake designs feature a blunt forebody shielding the payload from the intense heat generated during atmospheric entry. Although this forebody will absorb the largest portion of the heat pulse, accurate prediction of heating in the near wake is of great importance, since large local heating values can occur at points of shear-layer impingement. In order to address the various issues associated with these blunt-body wake flowfields, the Advisory Group for Aerospace Research and Development (AGARD) formed Working Group 18 in 1992. One of the objectives of this activity was to examine real-gas effects in high-speed flow fields around a 70 deg. blunted cone. To date, many researchers have conducted experiments using this geometry in various facilities, such as the Large Energy National Shock (LENS) tunnel at Cubric/Calspan and the HEG shock tunnel at DLR-Goettingen. Several computational studies have also been conducted in concert with these tests. Many of the experimental results have indicated the possible presence of a transitional shear layer through a large increase in heat transfer downstream of the reattachment point. The presence of transition could in fact lead to much higher peak heating than if the separated flow is entirely laminar or turbulent. In the shock-tunnel tests, however, it is difficult to separate such viscous-flow phenomena from real-gas effects. In order to help make this distinction, Horvath et al. recently conducted a set of experiments in the NASA Langley 20-Inch Mach 6 Tunnel, and compared the results to laminar Navier-Stokes calculations. They found heat-transfer distributions similar to those obtained in the high-enthalpy facilities, with the measured peak heating along the sting support markedly greater than that predicted by the laminar computations. These trends point to the need to find transitional and turbulent computational solutions for these flowfields.

Nance, Robert P.↗

NASA NTRS1997

Role of Boundary Conditions in Monte Carlo Simulation of MEMS Devices

A study is made of the issues surrounding prediction of microchannel flows using the direct simulation Monte Carlo method. This investigation includes the introduction and use of new inflow and outflow boundary conditions suitable for subsonic flows. A series of test simulations for a moderate-size microchannel indicates that a high degree of grid under-resolution in the streamwise direction may be tolerated without loss of accuracy. In addition, the results demonstrate the importance of physically correct boundary conditions, as well as possibilities for reducing the time associated with the transient phase of a simulation. These results imply that simulations of longer ducts may be more feasible than previously envisioned.

Nance, Robert P.↗

NASA NTRS1997

Numerical Solutions of the Complete Navier-Stokes Equations

This report details the development of a new two-equation turbulence closure model based on the exact turbulent kinetic energy k and the variance of vorticity, zeta. The model, which is applicable to three dimensional flowfields, employs one set of model constants and does not use damping or wall functions, or geometric factors.

Robinson, David F.↗

NASA NTRS1997

A New Approach for Determining Onset of Transition

The final report consists of three papers which outline and demonstrate the new method for determining transition onset. The procedure developed under this grant requires specification of the instability mechanism, i.e., Tollmien-Schlichting or crossflow, that leads to transition. The attached papers are entitled: 'An Alternative to the e(sup n) Method for Determining Onset of Transition', 'Transition Model for Swept Wing Flows', and 'A Transition Closure Model for Predicting Transition Onset'.

Hassan, H. A.↗

NASA NTRS1997

Transition and Turbulence Modeling for Blunt-Body Wake Flows

This study attempts t o improve the modeling and computational prediction of high- speed transitional wake flows. The recently developed kappa - zeta (Enstrophy) turbulence model is coupled with a newly developed transition prediction method and implemented in an implicit flow solver well-suited to hypersonic flows. In this model, transition onset is determined as part of the solution. Results obtained using the new model for a 70- deg blunted cone/sting geometry demonstrate better agreement with experimental heat- transfer measurements when compared to laminar calculations as well as solutions using the kappa - omega model. Results are also presented for the situation where transition onset is preselected. It is shown that, in this case, results are quite sensitive to location of the transition point.

Nance, Robert P.↗

NASA NTRS1996

Numerical solutions of the complete Navier-Strokes equations

This report describes the development of an enstrophy model capable of predicting turbulence separation and its application to two airfoils at various angles of attack and Mach numbers. In addition, a two equation kappa-xi model with a tensor eddy viscosity was developed. Plans call for this model to be used in calculating three dimensional turbulent flows.

Hassan, H. A.↗

NASA NTRS1995

Numerical solutions of the complete Navier-Stokes equations

A compressible turbulence model based on the kappa-zeta (enstrophy) model of Robinson et al. is developed. The model incorporates a new compressibility correction model that is quite different from previous models. Results are compared to the mixing layer data compiled by Settles and Dodson. It is shown that the predictions of the theory are in good agreement with experiment. Moreover, the present model is superior to the kappa-epsilon model with traditional compressibility corrections.

Hassan, H. A.↗

NASA NTRS1995

Effects of Chemistry on Blunt-Body Wake Structure

Results of a numerical study are presented for hypersonic low-density flow about a 70-deg blunt cone using direct simulation Monte Carlo (DSMC) and Navier-Stokes calculations. Particular emphasis is given to the effects of chemistry on the near-wake structure and on the surface quantities and the comparison of the DSMC results with the Navier-Stokes calculations. The flow conditions simulated are those experienced by a space vehicle at an altitude of 85 km and a velocity of 7 km/s during Earth entry. A steady vortex forms in the near wake for these freestream conditions for both chemically reactive and nonreactive air gas models. The size (axial length) of the vortex for the reactive air calculations is 25% larger than that of the nonreactive air calculations. The forebody surface quantities are less sensitive to the chemistry than the base surface quantities. The presence of the afterbody has no effect on the forebody flow structure or the surface quantities. The comparisons of DSMC and Navier-Stokes calculations show good agreement for the wake structure and the forebody surface quantities.

Dogra, Virendra K.↗

NASA NTRS1994

Parallel DSMC Solution of Three-Dimensional Flow Over a Finite Flat Plate

This paper describes a parallel implementation of the direct simulation Monte Carlo (DSMC) method. Runtime library support is used for scheduling and execution of communication between nodes, and domain decomposition is performed dynamically to maintain a good load balance. Performance tests are conducted using the code to evaluate various remapping and remapping-interval policies, and it is shown that a one-dimensional chain-partitioning method works best for the problems considered. The parallel code is then used to simulate the Mach 20 nitrogen flow over a finite-thickness flat plate. It is shown that the parallel algorithm produces results which compare well with experimental data. Moreover, it yields significantly faster execution times than the scalar code, as well as very good load-balance characteristics.

Nance, Robert P.↗