SEARCH · Engineering Papers
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Several turbulence models are studied and evaluated. All models are compared to each other and experimental values for flat plate boundary layer flow and channel flow. Among the models examined are: K-epsilon, k-omega, q-omega, and k-tau turbulence models. All results are presented in viewgraph format.
The constant C(sub mu) and the near-wall damping function f(sub mu) in the eddy-viscosity relation of the kappa-epsilon model are evaluated from direct numerical simulation (DNS) data for developed channel and boundary layer flow at two Reynolds numbers each. Various existing f(sub mu) model functions are compared with the DNS data, and a new function is fitted to the high-Reynolds-number channel flow data. The epsilon-budget is computed for the fully developed channel flow. The relative magnitude of the terms in the epsilon-equation is analyzed with the aid of scaling arguments, and the parameter governing this magnitude is established. Models for the sum of all source and sink terms in the epsilon-equation are tested against the DNS data, and an improved model is proposed.
A Boeing TR77 airfoil associated with the Advanced Technology Airfoil Test (ATAT) program was tested in the Langley 0.3 m Transonic Cryogenic Tunnel. Limited analysis of the data indicated that increasing Reynolds number for a fixed Mach number resulted in increased normal-force, nose-down pitching moment, and decreased drag coefficient. Increasing Mach number while keeping the Reynolds number constant yielded the expected increase in normal-force slopes, nose-down pitching moment coefficients, and decrease in angle of attack associated with maximum normal-force coefficient. Turbulent boundary layer flow was achieved over the airfoil at low Reynolds numbers for the test Mach number range using aluminum discs.
Measurements of Stanton numbers, velocity profiles, temperature profiles, and turbulence intensity profiles are reported for turbulent flat plate boundary layer flows with a step change in surface roughness. The first 0.9 m length of the test surface is roughened with 1.27 mm diameter hemispheres spaced 2 base diameters apart in a staggered array. The remaining 1.5 m length is smooth. The experiments show that the step change from a rough to a smooth surface has a dramatic effect on the convective heat transfer. In many cases, the Stanton number drops below the smooth-wall correlation immediately downstream of the change in roughness. The Stanton number measurements are compared with predictions using the discrete element method with excellent results.
The stability of a high-speed, axisymmetric boundary layer is investigated using secondary instability theory and direct numerical simulation. Parametric studies based on the temporal secondary instability theory identify subharmonic secondary instability as a likely path to transition on a cylinder at Mach 4.5. The theoretical predictions are validated by direct numerical simulation at temporally-evolving primary and secondary disturbances in an axisymmetric boundary-layer flow. At small amplitudes of the secondary disturbance, predicted growth rates agree to several significant digits with values obtained from the spectrally-accurate solution of the compressible Navier-Stokes equations. Qualitative agreement persists to large amplitudes of the secondary disturbance. Moderate transverse curvature is shown to significantly affect the growth rate of axisymmetric second mode disturbances, the likely candidates of primary instability. The influence of curvature on secondary instability is largely indirect but most probably significant, through modulation of the primary disturbance amplitude. Subharamonic secondary instability is shown to be predominantly inviscid in nature, and to account for spikes in the Reynolds stress components at or near the critical layer.
Near-wall modeling of the dissipation rate equation is investigated and its asymptotic behavior is studied in detail using a k-epsilon model. It is found that all existing modeled dissipation rate equations predict an incorrect behavior for the dissipation rate near a wall. An improvement is proposed and the resulting near-wall dissipation rate distribution is found to be similar to that given by numerical simulation data. To further validate the improved k-epsilon model, it is used to calculate flat-plate turbulent boundary-layer flows at high- as well as low-turbulence Reynolds numbers, and the results are compared with measurements, numerical simulation data, and the calculations of three different two-equation models. These comparisons show that all the models tested give essentially the same flow properties away from the wall; significant differences only occur in a region very close to the wall. In this region, the calculations of the improved k-epsilon model are in better agreement with measurements and numerical simulation data. In particular, the modeled distribution of the dissipation rate is significantly improved and a maximum is predicted at the wall instead of away from the wall. Furthermore, the improved k-epsilon model is found to be the most asymptotically consistent among the four different two-equation models examined.
Temporal direct numerical simulation of laminar breakdown via subharmonic secondary instability in high-speed axisymmetric boundary layers has been accomplished using a highly accurate, fully explicit algorithm which combines spectral collocation and high-order compact-difference techniques. Numerical test cases confirm that subharmonic secondary instability is confirmed to be a viable path to transition in high-speed boundary-layer flow. Secondary instability is shown to account for peaks in the Reynolds stresses at or near the critical layer which are not possible from the second-mode primary instability alone. Reynolds stresses spatially reconstructed from the temporal model via the Gaster transformation show a 'spreading angle' of about 12 deg, in qualitative agreement with experimental findings. The rate of broadening of the Reynolds stress peak is a strongly nonlinear phenomenon which cannot be reproduced by secondary instability theory.
The influence of convection on the development of morphological instability at the liquid-solid interface during directional solidification in a positive thermal gradient has been examined in Pb-10 wt pct Sn and succinonitrile-1.9 wt pct acetone. The onset of interfacial breakdown occurs at higher growth speeds in the presence of convection. The linear stability analysis due to Favier and Rouzaud which uses the 'deformable' mass flow boundary layer concept shows a good agreement with the experimentally observed behavior.
Report presents second part of two-part comprehensive study of coherent motions in, and structures of, turbulent boundary-layer flows. Study critically reviews current fragmented knowledge with view toward eventual unification and understanding of physical causes and statistical properties of turbulence phenomena.
A new numerical method is developed for the solution of the two-dimensional, steady Navier-Stokes equations. The method that is presented differs in significant ways from the established numerical methods for solving the Navier-Stokes equations. The major differences are described. First, the focus of the present method is on satisfying flux conservation in an integral formulation, rather than on simulating conservation laws in their differential form. Second, the present approach provides a unified treatment of the dependent variables and their unknown derivatives. All are treated as unknowns together to be solved for through simulating local and global flux conservation. Third, fluxes are balanced at cell interfaces without the use of interpolation or flux limiters. Fourth, flux conservation is achieved through the use of discrete regions known as conservation elements and solution elements. These elements are not the same as the standard control volumes used in the finite volume method. Fifth, the discrete approximation obtained on each solution element is a functional solution of both the integral and differential form of the Navier-Stokes equations. Finally, the method that is presented is a highly localized approach in which the coupling to nearby cells is only in one direction for each spatial coordinate, and involves only the immediately adjacent cells. A general third-order formulation for the steady, compressible Navier-Stokes equations is presented, and then a Newton's method scheme is developed for the solution of incompressible, low Reynolds number channel flow. It is shown that the Jacobian matrix is nearly block diagonal if the nonlinear system of discrete equations is arranged approximately and a proper pivoting strategy is used. Numerical results are presented for Reynolds numbers of 100, 1000, and 2000. Finally, it is shown that the present scheme can resolve the developing channel flow boundary layer using as few as six to ten cells per channel width, depending on the Reynolds number.
The compressible Navier-Stokes equations are solved using spectral collocation and high-order compact-difference techniques to simulate the laminar breakdown in high-speed, axisymmetric boundary-layer flow. Mach 4.5 flow along a hollow cylinder and Mach 6.8 flow along a sharp cone are considered. Data obtained replicate two previously unexplained phenomena, namely, the appearance of so-called 'rope-like waves' and 'the precursor transition effect', in which transitional flow originates near the critical layer well upstream of the transition location at the wall. The numerical data also reveal that neither of these effects can be explained, even qualitatively, by linear stability theory alone. It is shown that rope-like appearance arises from secondary instability. Certain features of the precursor transition effect also emerge from secondary instability but its nature is revealed to be fundamentally nonlinear.
The two main hypotheses concerning the origin of compositional heterogeneities in magma chambers are discussed: (1) models in which the development of compositional zonation is simultaneous with the birth and growth of the magma body and (2) models in which zonation develops within an initially homogeneous batch of magma. The paper presents an overview of the geological possibilities and evaluates them on the basis of current research. Calculations are presented for boundary-layer flow in isothermal ternary component systems, and it is demonstrated that multicomponent diffusion effects may be very significant, as was earlier suggested by Trial and Spera (1988).
In modeling the laminar-turbulent transition region, the designer depends largely on benchmark data from experiments and/or direct numerical simulations that are usually extremely expensive. An understanding of the evolution of the Reynolds stresses, turbulent kinetic energy, and quantifies in the transport equations like the dissipation and production is essential in the modeling process. The secondary instability theory and the parabolized stability equations method are used to calculate these quantities, which are then compared with corresponding quantities calculated from available direct numerical simulation data for the incompressible boundary-layer flow of laminar-turbulent transition conditions. The potential of the secondary instability theory and the parabolized stability equations approach in predicting these quantities is discussed; results indicate that inexpensive data that are useful for transition modeling in the early stages of the transition region can be provided by these tools.
The constant C sub mu and the near-wall damping function f sub mu in the eddy-viscosity relation of the k-epsilon model are evaluated from direct numerical simulation (DNS) data for developed channel and boundary layer flow at two Reynolds numbers each. Various existing f sub mu model functions are compared with the DNS data, and a new function is fitted to the high-Reynolds-number channel flow data. The epsilon-budget is computed for the fully developed channel flow. The relative magnitude of the terms in the epsilon-equation is analyzed with the aid of scaling arguments, and the parameter governing this magnitude is established. Models for the sum of all source and sink terms in the epsilon-equation are tested against the DNS data, and an improved model is proposed.
An incompressible three-dimensional laminar boundary-layer flow over a swept wing is used as a model to study both the wall-curvature and streamline-curvature effects on the stationary crossflow instability. The basic state is obtained by solving the full Navier-Stokes (N-S) equations numerically. The linear disturbance equations are cast on a fixed, body-intrinsic, curvilinear coordinate system. Those nonparallel terms which contribute mainly to the streamline-curvature effect are retained in the formulation of the disturbance equations and approximated by their local finite difference values. The resulting eigenvalue problem is solved by a Chebyshev collocation method. The present results indicate that the convex wall curvature has a stabilizing effect, whereas the streamline curvature has a destabilizing effect. A validation of these effects with an N-S solution for the linear disturbance flow is provided.
To analyze ignition transient of Space Shuttle solid rocket motor, a transient two-dimensional numerical model based on turbulent compressible Navier-Stokes equations in a generalized coordinate system was developed. One-dimensional numerical models (Peretz et al., 1973; Han, 1992; Pardue and Han, 1992) with empirical correlations data obtained from steady turbulent boundary layer flows agrees reasonably well with test rocket data by adjusting a few parameters. However, a 1D model can not provide a physical insight into the complex multidimensional thermal fields and flowfields in the chamber and the converging-diverging rocket nozzle. As an interim step, a 2D model was developed and compared with test data. A modified version of SIMPLE algorithm was used for the numerical model, and the standard k-epsilon model with a wall function was used for turbulence closure. Transient flowfields and thermal fields in the combustion chamber and the attached nozzle were obtained for a selected rocket geometry and propellant. Transient behaviors of the flow and thermal fields were analyzed, and were found to be in good agreement with physical expectations.
The subgrid scales are modeled dynamically in a large-eddy simulation of transitional boundary-layer flow along a hollow cylinder at a Mach number of 4.5. The behavior of the dynamic-model coefficients, which is determined from local information in the resolved field, is investigated through both an a priori test with direct numerical simulation data for the same case and a complete large-eddy simulation. Both contractions proposed by Germano et al. (1991) and Lilly (1992) are used for the unique determination of the coefficients of the dynamic model, and their results are compared and assessed. The behavior, as well as the energy cascade of the subgridscale field structure, is investigated at various stages of the transition process.
With the recent revitalization of high speed flow research, compressibility presents a new set of challenging problems to turbulence researchers. Questions arise as to what extent compressibility affects turbulence dynamics, structures, the Reynolds stress-mean velocity (constitutive) relation, and the accompanying processes of heat transfer and mixing. In astrophysical applications, compressible turbulence is believed to play an important role in intergalactic gas cloud dynamics and in accretion disk convection. Understanding and modeling of the compressibility effects in free shear flows, boundary layers, and boundary layer/shock interactions is discussed.