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.
Computer program predicts growth of boundary layers along any surface where air is flowing. Two integral methods, Cohen-Reshotko and Sasman-Cresici, calculate laminar boundary layers and turbulent boundary layers, respectively; Schlichting-Granville method predicts transition from laminar to turbulent flow.
Turbulence is considered as a method for converting the energy observed as mass motions during the impulsive phase into thermal energy observed during the gradual phase of solar flares. The kinetic energy of the large-scale eddies driven by the upflowing material continuously cascades to smaller scale eddies until viscosity is able to convert it into thermal energy. The general properties of steady state, homogeneous, fluid turbulence is a nonmagnetic plasma and the properties of turbulent decay are reviewed. The time-dependent behavior of the velocities and energies observed by the X-Ray Polychromator (XRP) instrument on the SMM during the November 5, 1980 flare are compared with the properties of turbulence. This study indicates that turbulence may play a role in flare energies and may account for a fraction of the total amount of thermal energy observed during the gradual phase. The rate at which the observed flare velocities decrease is consistent with the decay of turbulent energy but may be too rapid to account for the entire time delay between the impulsive and gradual phases.
The trend for decreasing the drag of aircraft is retention of laminar flow in the boundary layer over a large portion of the surface. The laminar boundary layer was studied in a low turbulence wind tunnel for low subsonic velocities. The method used and results of measurements of very low levels of turbulence are presented. Measurements were performed by a constant-resistance thermoanemometer.
Engineering codes are needed which can calculate convective heating rates accurately and expeditiously on the surfaces of high-speed vehicles. One code which has proven to meet these needs is the Langley Approximate Three-Dimensional Convective Heating (LATCH) code. It uses the axisymmetric analogue in an integral boundary-layer method to calculate laminar and turbulent heating rates along inviscid surface streamlines. It requires the solution of the inviscid flow field to provide the surface properties needed to calculate the streamlines and streamline metrics. The LATCH code has been used with inviscid codes which calculated the flow field on structured grids, Several more recent inviscid codes calculate flow field properties on unstructured grids. The present research develops a method to calculate inviscid surface streamlines, the streamline metrics, and heating rates using the properties calculated from inviscid flow fields on unstructured grids. Mr. Chris Riley, prior to his departure from NASA LaRC, developed a preliminary code in the C language, called "UNLATCH", to accomplish these goals. No publication was made on his research. The present research extends and improves on the code developed by Riley. Particular attention is devoted to the stagnation region, and the method is intended for programming in the FORTRAN 90 language.
The scalar form of the approximate factorization method was used to develop a new code for the solution of three-dimensional internal laminar and turbulent compressible flows. The Navier-Stokes equations in their Reynolds-averaged form are iterated in time until a steady solution is reached. Evidence is given to the implicit and explicit artificial damping schemes that proved to be particularly efficient in speeding up convergence and enhancing the algorithm robustness. A conservative treatment of these terms at domain boundaries is proposed in order to avoid undesired mass and/or momentum artificial fluxes. Turbulence effects are accounted for by the zero-equation Baldwin-Lomax turbulence model and the q-omega two-equation model. For the first, an investigation on the model behavior in case of multiple boundaries is performed. The flow in a developing S-duct is then solved in the laminar regime at Reynolds number (Re) 790 and in the turbulent regime at Re=40,000 using the Baldwin-Lomax model . The Stanitz elbow is then solved using an inviscid version of the same code at M(sub inlet)=0.4. Grid dependence and convergence rate are investigated showing that for this solver the implicit damping scheme may play a critical role for convergence characteristics. The same flow at Re=2.5x10(exp 6) is solved with the Baldwin-Lomax and the q-omega models. Both approaches showed satisfactory agreement with experiments, although the q-omega model is slightly more accurate.
An approximate method for calculating the noise generated by a turbulent flow within a semi-infinite duct of arbitrary cross section is developed. It is based on a previously derived high-frequency solution to Lilley's equation, which describes the sound propagation in transversely-sheared mean flow. The source term is simplified by assuming the turbulence to be axisymmetric about the mean flow direction. Numerical results are presented for the special case of a ring source in a circular duct with an axisymmetric mean flow. They show that the internally generated noise is suppressed at sufficiently large upstream angles in a hard walled duct, and that acoustic liners can significantly reduce the sound radiated in both the upstream and downstream regions, depending upon the source location and Mach number of the flow.
An approximate method for calculating the noise generated by a turbulent flow within a semi-infinite duct of arbitrary cross section is developed. It is based on a previously derived high-frequency solution to Lilley's equation, which describes the sound propagation in a transversely-sheared mean flow. The source term is simplified by assuming the turbulence to be axisymmetric about the mean flow direction. Numerical results are presented for the special case of a ring source in a circular duct with an axisymmetric mean flow. They show that the internally generated noise is suppressed at sufficiently large upstream angles in a hard walled duct, and that acoustic liners can significantly reduce the sound radiated in both the upstream and downstream regions, depending upon the source location and Mach number of the flow.
We investigate a relativistic adaptation of the lattice Boltzmann method that reproduces the equations of motion for a turbulent, two-dimensional, massless hydrodynamic system. The classical lattice Boltzmann method and its extension to relativistic fluid dynamics are described. The numeric formulation is evaluated using a zero-averaged stirring force introduced into the numerics to induce turbulence, and the flow characteristics produced are compared to properties of a classical turbulent hydrodynamic flow. The model can reasonably be expected to offer quantitative simulations of electron fluid flows in graphene or Kagome lattices.
The Johnson-King turbulence model that is a viable method for calculating two dimensional transonic separated flows was extended into three dimensions. The implementation was done for Navier Stokes flow solvers written in general curvilinear coordinates. The present approach used in turbulence modeling is based on streamwise integration of an ordinary differential equation (o.d.e.) that governs the maximum Reynolds shear stress behavior. Streamwise integration of the o.d.e. approach was found to offer great mathematical simplicity and economy for three dimensional Navier Stokes methods. Thus, the new method is quick, simple, and very cheap. The new method was first checked against the data of a well known transonic axisymmetric bump experiment, and a good agreement was obtained. Later, the new method was used to compute the flow around a low aspect ratio wing in a transonic wind tunnel. Finally it was employed to study the nonequilibrium turbulence effects on the transonic vortical flows about a 65 deg sweep round leading edge delta wing.
Predictive simulations of high-pressure sprays require accurate representation of the turbulent gaseous flow field generated by liquid jet. Typically, the accuracy that can be obtained with low-order numerical methods (e.g. finite volume, finite element) is limited by stability issues in fine grids and the order of convergence of the method. In this work, we resolve the turbulent flow field in an Eulerian manner using the high-order spectral element method, coupled with a Lagrangian parcels approach to model the atomizing liquid jet. Large eddy simulations of single-hole sprays under non-evaporative conditions were conducted and compared against experimental data from Margot et al. (2008) and Spray A data from the Engine Combustion Network. The sensitivity of liquid penetration and droplet sizes to different breakup model parameters was studied. The effect of different numerical parameters, such as polynomial order of the solution (grid resolution), on liquid penetration was also analyzed. The method achieved grid-independent results using p-refinement, achieving finer resolution (by a factor of x 1.7 - x 3.5) in the gas-phase solution than in state-of-the-art simulations using the finite-volume method. Here the results showed good agreement with experimental data, demonstrating the ability of the current method to accurately capture liquid penetration and the shape of the spray.
The roll-up of a wingtip vortex, at Reynolds number based on chord of 4.6 million was studied with an emphasis on suction side and near wake measurements. The research was conducted in a 32 in. x 48 in. low-speed wind tunnel. The half-wing model had a semi-span of 36 in. a chord of 48 in. and a rounded tip. Seven-hole pressure probe measurements of the velocity field surrounding the wingtip showed that a large axial velocity of up to 1.77 U(sub infinity) developed in the vortex core. This level of axial velocity has not been previously measured. Triple-wire probes have been used to measure all components of the Reynolds stress tensor. It was determined from correlation measurements that meandering of the vortex was small and did not appreciably contribute to the turbulence measurements. The flow was found to be turbulent in the near-field (as high as 24 percent RMS w - velocity on the edge of the core) and the turbulence decayed quickly with streamwise distance because of the nearly solid body rotation of the vortex core mean flow. A streamwise variation of the location of peak levels of turbulence, relative to the core centerline, was also found. Close to the trailing edge of the wing, the peak shear stress levels were found at the edge of the vortex core, whereas in the most downstream wake planes they occurred at a radius roughly equal to one-third of the vortex core radius. The Reynolds shear stresses were not aligned with the mean strain rate, indicating that an isotropic-eddy-viscosity based prediction method cannot accurately model the turbulence in the cortex. In cylindrical coordinates, with the origin at the vortex centerline, the radial normal stress was found to be larger than the circumferential.
A theoretical study was conducted to develop an analytical prediction method for helicopter main rotor noise due to the ingestion of atmospheric turbulence. This study incorporates an atmospheric turbulence model, a rotor mean flow contraction model and a rapid distortion turbulence model which together determine the statistics of the non-isotropic turbulence at the rotor plane. Inputs to the combined mean inflow and turbulence models are controlled by atmospheric wind characteristics and helicopter operating conditions. A generalized acoustic source model was used to predict the far field noise generated by the non-isotropic flow incident on the rotor. Absolute levels for acoustic spectra and directivity patterns were calculated for full scale helicopters, without the use of empirical or adjustable constants. Comparisons between isotropic and non-isotropic turbulence at the rotor face demonstrated pronounced differences in acoustic spectra. Turning and contraction of the flow for hover and low speed vertical ascent cases result in a 3 dB increase in the acoustic spectrum energy and a 10 dB increase in tone levels. Compared to trailing edge noise, turbulence ingestion noise is the dominant noise mechanism below approximately 30 rotor harmonics, while above 100 harmonics, trailing edge noise levels exceed turbulence ingestion noise by 25 dB.
The considered investigation is concerned with the development of an efficient computational method for obtaining a physical understanding of an internal turbulent field. The employed approach makes use of a 'two equation' type model for the turbulence to obtain the numerical solution of a two-dimensional confined turbulent flow. The mean flow governing equations are considered along with the governing equation of the mean temperature and concentrations, and the boundary conditions. The numerical procedure for solving the turbulent flow is discussed, taking into account an approximation to the nonlinear terms, and the inner and outer coupling. Attention is given to a stability convergence analysis, the stability characteristics, and computational examples.
Procedure for computing subsonic, turbulent flow in turbofan lobe mixers was extended to allow consideration of flow fields in which a swirl component of velocity may be present. Additional, an optional k-lambda turbulence model was added to the procedure. The method of specifying the initial flow field was also modified, allowing parametric specification or radial secondary flow velocities, and making it possible to consider initial flow fields which have significant inlet secondary flow vorticity. A series of example calculations was performed which demonstrate the various capabilities of the modified code. These calculations demonstrate the effects of initial secondary flows of various magnitudes, the effects of swirl, and the effects of turbulence model on the mixing process. The results of these calculations indicate that the initial secondary flows, presumed to be generated within the lobes, play a dominant role in the mixing process, and that the predicted results are relatively insensitive to the turbulence model used.
Phase-diversity techniques provide a novel observational method for overcomming the effects of turbulence and instrument-induced aberrations in ground-based astronomy. Two implementations of phase-diversity techniques that differ with regard to noise model, estimator, optimization algorithm, method of regularization, and treatment of edge effects are described. Reconstructions of solar granulation derived by applying these two implementations to common data sets are shown to yield nearly identical images. For both implementations, reconstructions from phase-diverse speckle data (involving multiple realizations of turbulence) are shown to be superior to those derived from conventional phase-diversity data (involving a single realization). Phase-diverse speckle reconstructions are shown to achieve near diffraction-limited resolution and are validated by internal and external consistency tests, including a comparison with a reconstruction using a well-accepted speckle-imaging method.
A hybrid method for the prediction of supersonic turbulent flows with combustion is developed consisting of a second order closure for the velocity field and a multi-scalar pdf method for the local thermodynamic state. It is shown that for non-premixed flames and chemical equilibrium mixture fraction, the logarithm of the (dimensionless) density, internal energy per unit mass and the divergence of the velocity have several advantages over other sets of scalars. The closure model is applied to a supersonic non-premixed flame burning hydrogen with air supplied by a supersonic coflow and the results are compared with a limited set of experimental data.
The primary objective of this research is to extend current capabilities of Large Eddy Simulations (LES) and Direct Numerical Simulations (DNS) for the computational analyses of high speed reacting flows. Our efforts in the first two years of this research have been concentrated on a priori investigations of single-point Probability Density Function (PDF) methods for providing subgrid closures in reacting turbulent flows. In the efforts initiated in the third year, our primary focus has been on performing actual LES by means of PDF methods. The approach is based on assumed PDF methods and we have performed extensive analysis of turbulent reacting flows by means of LES. This includes simulations of both three-dimensional (3D) isotropic compressible flows and two-dimensional reacting planar mixing layers. In addition to these LES analyses, some work is in progress to assess the extent of validity of our assumed PDF methods. This assessment is done by making detailed companions with recent laboratory data in predicting the rate of reactant conversion in parallel reacting shear flows. This report provides a summary of our achievements for the first six months of the third year of this program.
This document describes the reasoning and trade studies used to evaluate tools for constructing the aerodynamic lineload databases for the liftoff and transition phases of flight for the Space Launch System. Three computational fluid dynamics codes (USM3D, FUN3D, Kestrel) were investigated with various turbulence models, as well as detached eddy simulation variants for the launch vehicle in free air and in proximity to the tower. Decisions were made mostly based on results from brief developmental studies performed in response to specific, unforeseen challenges that were encountered in the analysis of a given configuration. The need to develop databases in a timely manner, as well as accurately capture the expected leeward-wake flowfield characteristics, led to the selection of the Kestrel flow solver with its delayed detached eddy simulation method, the Spalart-Allmaras turbulence model, and the adaptive mesh refinement capability in the off-body Cartesian grid region.