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.
Explore the source record for details and available documents.
The present paper gives the results of measurements of the resistance of a plate placed parallel to the flow in the range of Reynolds numbers from 10 to 2300; in this range the resistance deviates from the formula of Blasius. The lower limit of validity of the Blasius formula is determined and also the increase in resistance at the edges parallel to the flow in the case of a plate of finite width.
Previous work at NASA LeRC has shown that flow distortions in aircraft engine inlet ducts can be significantly reduced by mounting vortex generators, or small wing sections, on the inside surface of the engine inlet. The placement of the vortex generators is an important factor in obtaining the optimal effect over a wide operating envelope. In this regard, the only alternative to a long and expensive test program which would search out this optimal configuration is a good prediction procedure which could narrow the field of search. Such a procedure has been developed in collaboration with NASA LeRC, and results obtained by NASA personnel indicate that it shows considerable promise for predicting the viscous turbulent flow in engine inlet ducts in the presence of vortex generators. The prediction tool is a computer code which numerically solves the reduced Navier-Stokes equations and so is commonly referred to as RNS3D. Obvious deficiencies in RNS3D have been addressed in previous work. Primarily, it is known that the predictions of the mean velocity field of a turbulent boundary layer flow approaching separation are not in good agreement with data. It was suggested that the use of an algebraic mixing-length turbulence model in RNS3D is at least partly to blame for this. Additionally, the current turbulence model includes an assumption of isotropy which will ultimately fail to capture turbulence-driven secondary flow known to exist in noncircular ducts.
Chemical and thermal nonequilibrium phenomena are studied in the stagnation region of a hypervelocity blunt body. This investigation is motivated by the need to predict the heat-transfer rate to the leading edge of aeromaneuvering orbital transfer vehicles. Flight speeds of approximately 10 km/s at altitudes of approximately 80 km are considered for body radii of 1-50 cm. The analysis is based on continuum theory and is applicable to the viscous and incipient merged layer regimes of rarefied flow. A two-species, two-temperature gas model is assumed. Comparisons are made with previous theories, experimental data, and results based on the thermodynamic equilibrium assumption. The equation accounting for vibrational nonequilibrium is presented and its effects on flow properties are discussed. Parameters requiring further investigation are identified. Preliminary results indicate that the inclusion of vibrational relaxation has little effect on the heat-transfer rate for a fully catalytic surface. However, vibrational nonequilibrium may increase the heat-transfer rate to a noncatalytic surface, depending on the degree of nonequilibrium.
A sailplane is currently being developed at NASA's Dryden Flight Research Center to support a high altitude flight experiment. The purpose of the experiment is to measure the performance characteristics of an airfoil at altitudes between 100,000 and 70,000 feet at Mach numbers between 0.65 and 0.5. The airfoil lift and drag are measured from pilot and static pressures. The location of the separation bubble and vortex shedding are measured from a hot film strip. The details of the flight experiment are presented. A comparison of several estimates of the airfoil performance is also presented. The airfoil, APEX-16, was designed by Drela (MIT) with his MSES code. A two dimensional Navier-Stokes analysis has been performed by Tatineni and Zhong (UCLA) and another at the Dryden Flight Research Center. The role these analysis served to define the experiment is discussed.
The mechanisms of sound generation in a Mach 0.9, Reynolds number 3600 turbulent jet are investigated by direct numerical simulation. Details of the numerical method are briefly outlined and results are validated against an experiment at the same flow conditions. Lighthill's theory is used to define a nominal acoustic source in the jet, and a numerical solution of Lighthill's equation is compared to the simulation to verify the computational procedures. The acoustic source is Fourier transformed in the axial coordinate and time and then filtered in order to identify and separate components capable of radiating to the far field. This procedure indicates that the peak radiating component of the source is coincident with neither the peak of the full unfiltered source nor that of the turbulent kinetic energy. The phase velocities of significant components range from approximately 5% to 50% of the ambient sound speed which calls into question the commonly made assumption that the noise sources convect at a single velocity. Space-time correlations demonstrate that the sources are not acoustically compact in the streamwise direction and that the portion of the source that radiates at angles greater than 45 deg. is stationary. Filtering non-radiating wavenumber components of the source at single frequencies reveals that a simple modulated wave forms for the source, as might be predicted by linear stability analysis. At small angles from the jet axis the noise from these modes is highly directional, better described by an exponential than a standard Doppler factor.
Airfoils at high Reynolds numbers, in general, have small separation bubbles that are usually confined to the leading edge. Since the Reynolds number is large, the turbulence model for the transition region between the laminar and turbulent flow is not important. Furthermore, the onset of transition occurs either at separation or prior to separation and can be predicted satisfactorily by empirical correlations when the incident angle is small and can be assumed to correspond to laminar separation when the correlations do not apply, i.e., at high incidence angles.
NASA, FAA, ONERA, and other partner organizations have embarked on a significant, collaborative research effort to address the technical challenges associated with icing on large scale, three-dimensional swept wings. These are extremely complex phenomena important to the design, certification and safe operation of small and large transport aircraft. There is increasing demand to balance trade-offs in aircraft efficiency, cost and noise that tend to compete directly with allowable performance degradations over an increasing range of icing conditions. Computational fluid dynamics codes have reached a level of maturity that they are being proposed by manufacturers for use in certification of aircraft for flight in icing. However, sufficient high-quality data to evaluate their performance on iced swept wings are not currently available in the public domain and significant knowledge gaps remain.
Aerodynamic assessment of icing effects on swept wings is an important component of a larger effort to improve three-dimensional icing simulation capabilities. An understanding of ice-shape geometric fidelity and Reynolds and Mach number effects on the iced-wing aerodynamics is needed to guide the development and validation of ice-accretion simulation tools. To this end, wind-tunnel testing and computational flow simulations were carried out for an 8.9%-scale semispan wing based upon the Common Research Model airplane configuration. The wind-tunnel testing was conducted at the Wichita State University 7 ft x 10 ft Beech wind tunnel from Reynolds numbers of 0.8×10(exp 6) to 2.4×10(exp 6) and corresponding Mach numbers of 0.09 to 0.27. This paper presents the results of initial studies investigating the model mounting configuration, clean-wing aerodynamics and effects of artificial ice roughness. Four different model mounting configurations were considered and a circular splitter plate combined with a streamlined shroud was selected as the baseline geometry for the remainder of the experiments and computational simulations. A detailed study of the clean-wing aerodynamics and stall characteristics was made. In all cases, the flow over the outboard sections of the wing separated as the wing stalled with the inboard sections near the root maintaining attached flow. Computational flow simulations were carried out with the ONERA elsA software that solves the compressible, three-dimensional RANS equations. The computations were carried out in either fully turbulent mode or with natural transition. Better agreement between the experimental and computational results was obtained when considering computations with free transition compared to turbulent solutions. These results indicate that experimental evolution of the clean wing performance coefficients were due to the effect of three-dimensional transition location and that this must be taken into account for future data analysis. This research also confirmed that artificial ice roughness created with rapid-prototype manufacturing methods can generate aerodynamic performance effects comparable to grit roughness of equivalent size when proper care is exercised in design and installation. The conclusions of this combined experimental and computational study contributed directly to the successful implementation of follow-on test campaigns with numerous artificial ice-shape configurations for this 8.9% scale model.
Aerodynamic assessment of icing effects on swept wings is an important component of a larger effort to improve three-dimensional icing simulation capabilities. An understanding of ice-shape geometric fidelity and Reynolds and Mach number effects on the iced-wing aerodynamics is needed to guide the development and validation of ice-accretion simulation tools. To this end, wind-tunnel testing and computational flow simulations were carried out for an 8.9 percent-scale semispan wing based upon the Common Research Model airplane configuration. The wind-tunnel testing was conducted at the Wichita State University 7 by 10 ft Beech wind tunnel from Reynolds numbers of 0.8×10(exp 6) to 2.4×10(exp 6) and corresponding Mach numbers of 0.09 to 0.27. This paper presents the results of initial studies investigating the model mounting configuration, clean-wing aerodynamics and effects of artificial ice roughness. Four different model mounting configurations were considered and a circular splitter plate combined with a streamlined shroud was selected as the baseline geometry for the remainder of the experiments and computational simulations. A detailed study of the clean-wing aerodynamics and stall characteristics was made. In all cases, the flow over the outboard sections of the wing separated as the wing stalled with the inboard sections near the root maintaining attached flow. Computational flow simulations were carried out with the ONERA elsA software that solves the compressible, threedimensional RANS equations. The computations were carried out in either fully turbulent mode or with natural transition. Better agreement between the experimental and computational results was obtained when considering computations with free transition compared to turbulent solutions. These results indicate that experimental evolution of the clean wing performance coefficients were due to the effect of three-dimensional transition location and that this must be taken into account for future data analysis. This research also confirmed that artificial ice roughness created with rapid-prototype manufacturing methods can generate aerodynamic performance effects comparable to grit roughness of equivalent size when proper care is exercised in design and installation. The conclusions of this combined experimental and computational study contributed directly to the successful implementation of follow-on test campaigns with numerous artificial ice-shape configurations for this 8.9 percent scale model.
In this paper we report the outcome of selected workshops organized as part of the NATO Applied Vehicle Technology (AVT)-313 activity Incompressible Laminar-to-Turbulent Flow Transition Study that focused on assessing the numerical and modeling accuracy of the γ−Reθ and γ transition models coupled to the k−ω Shear-Stress Transport (SST) two-equation eddy-viscosity model. Three different test cases involving nominally 2D flow configurations were selected: flow over a flat plate with two different levels of turbulence intensity at the inlet; flow around the Eppler 387 foil at a Reynolds number of 3×10^5 and angles of attack of 1 deg. and 7 deg. flow around the NACA 0015 foil at a Reynolds number of 1.8×10^5 and angles of attack of 5 deg. and10 deg. The flat plate flow conditions correspond to natural and by-pass transition, whereas the other two test cases include laminar separation bubbles that lead to separation-induced transition. For each test case, the selected quantities of interest include both integral and local flow quantities. Geometrically similar grids with a wide range of grid refinement ratios were generated for each of the test cases to allow the estimation of numerical uncertainties for all quantities of interest selected for this study. Several RANS flow solvers were used, employing common grids with the same boundary conditions and mathematical models. Therefore, it is possible to analyze the consistency of the results, i.e., to check if the intervals defined by the different numerical solutions with their respective uncertainties overlap with each other. Modeling errors can also be addressed for the selected flow quantities that have experimental data available. However, the experimental information available in these cases is not sufficient to guarantee that experiments and simulations are performed with the same settings. Nonetheless, the available experimental data is sufficient to guarantee that modeling errors are significantly reduced with the use of the transition models when compared to simulations performed using only the k−ω SST model.
Explore the source record for details and available documents.
Implicit large-eddy simulations of the self-similar regions of two compressible mixing layers at high Reynolds number and convective Mach numbers of 0.381 and 0.690 were carried out. Experimental data was used as the input into a synthetic eddy method turbulent inflow boundary condition to initiate simulations. Computational grids with nearly isotropic cells were used and the grid counts ranged from 24 to 270 million points. Agreement with experimental velocity and Reynolds stress data was good given the simplifications of the simulation. The decrease in mixing layer growth rate with increasing compressibility compares well with the literature. Budgets for the Reynolds stress transport equation were extracted from the simulations. It was demonstrated that the grids used did not resolve the dissipation term at the experimental Reynolds number. However, low Reynolds number simulations indicate that the numerical dissipation replaces the physical dissipation and the dissipation term in the transport equation can be represented by the summation term that quantifies the imbalance in the budget. The budgets for the high Reynolds number mixing layers agree well with previous low Reynolds number studies; indicating that the primary mechanism for the reduced growth rate and increased anistropy of the Reynolds stress tensor due to compressibility is a reduction in the pressure-strain and production terms. A new scaling based on the magnitude of the Reynolds stress tensor is proposed. This scaling provides better relative comparisons between the data and shows that the decrease in pressure-strain in the shear component is the primary driver in reducing the transverse normal stress.
An efficient and accurate method for solving the viscous shock layer equations for hypersonic flows over long slender bodies is presented. The two first order equations, continuity and normal momentum, are solved simultaneously as a coupled set. The flow conditions included are from high Reynolds numbers at low altitudes to low Reynolds numbers at high altitudes. For high Reynolds number flows, both chemical nonequilibrium and perfect gas cases are analyzed with surface catalytic effects and different turbulence models, respectively. At low Reynolds number flow conditions, corrected slip models are implemented with perfect gas case. Detailed comparisons are included with other predictions and experimental data.
The major hypersonic aerodynamic results obtained experimentally from the 0.010-scale (51-0) 140A/B orbiter vehicle are presented. The test was conducted over an angle-of-attack range from 20 to 35 degrees for nominal Mach numbers of 16 and 20. The Mach 16 data were obtained at Reynolds numbers, based on model length, of nominally 0.54 x 106 and 1.18 x 106. The Mach number 20 data were obtained at nominal Reynolds numbers of 0.08 x 106 and 0.27 x 106. Various elevon deflections and body flap deflections were tested for determination of control effectiveness. The test objectives were: (1) verify December 1973 aerodynamic Design Data Book estimates for Mach number greater than 10; (2) determine configuration forward and aft trim boundaries and control characteristics; (3) determine high Mach number - low Reynolds number effects; (4) determine if positive control surface flow separation is present and magnitude of separation effects.
Preconditioning methods have become the method of choice for the solution of flowfields involving the simultaneous presence of low Mach and transonic regions. It is well known that these methods are important for insuring accurate numerical discretization as well as convergence efficiency over various operating conditions such as low Mach number, low Reynolds number and high Strouhal numbers. For unsteady problems, the preconditioning is introduced within a dual-time framework wherein the physical time-derivatives are used to march the unsteady equations and the preconditioned time-derivatives are used for purposes of numerical discretization and iterative solution. In this paper, we describe the implementation of the preconditioned dual-time methodology in the OVERFLOW code. To demonstrate the performance of the method, we employ both simple and practical unsteady flowfields, including vortex propagation in a low Mach number flow, flowfield of an impulsively started plate (Stokes' first problem) arid a cylindrical jet in a low Mach number crossflow with ground effect. All the results demonstrate that the preconditioning algorithm is responsible for improvements to both numerical accuracy and convergence efficiency and, thereby, enables low Mach number unsteady computations to be performed at a fraction of the cost of traditional time-marching methods.
The immersed boundary method (IBM) with direct forcing is very popular in the simulation of rigid particulate flows. In the IBM, an interaction force is introduced at the interface between fluid and particle in order to approximate the no-slip boundary condition. The interaction force is calculated through dividing the velocity difference (or error) between fluid and particle at the interface by the time step. Here, a dynamic equation for the velocity difference is derived. Additionally, analyses on the dynamic equation provide a few new findings: (i) The interaction force is the solution of a least-squares error problem, with the direct implication that the Lagrangian marker distribution has no effect on the large scale flow structure once the distribution of Lagrangian markers become saturated along the interface (i.e., each marker remains properly correlated with all its neighbors); (ii) The Lagrangian volume-weight is a relaxation factor to control how fast the velocity error decays to the ideal value of zero; (iii) The optimal choice of the Lagrangian volume-weight is the largest value permissible by a stability condition. A comprehensive convergence analysis with regard to the spatial and temporal resolution is presented for the velocity error and also for the shear-stress and surface pressure. In three simple canonical problems, it is analytically and numerically shown that the IBM results converge to the theoretical solutions obtained with precise imposition of no-slip and no-penetration boundary conditions. It is observed that it is not necessary to match the Lagrangian marker volume-weight to that of the local Eulerian cell volume and in fact this matching leads to lower than optimal computational efficiency. However, it is found that extremely high Eulerian grid resolution and small time step have to be used to obtain high precision simulation results. Especially, the time step should be inversely proportional to the particle Reynolds number for low Reynolds number flows. For high frequency oscillation problems, the grid size needs to be reduced by a factor of the square root of the frequency, and the time step to be reduced by a factor of the frequency. The theoretical findings here can be used to alleviate the technical difficulties in simulating non-spherical particles by not requiring the Lagrangian marker distribution to match the Eulerian grids and also in the implementation of IBM on non-uniform Eulerian grids. The present work also provides simple practical guidance on the choice of temporal and spatial resolution so as to control the simulation error a priori.
The field of computational aeroacoustics (CAA), in which sound generation by unsteady flows is calculated from first principles, has recently emerged as a viable discipline. This paper describes one CAA technique which has been validated for low Mach number, low Reynolds number, two-dimensional flows. Results of its application to geometries producing both dipole and quadrupole type sound sources are presented along with theoretical and experimental comparisons. In addition, the future promise of this and other CAA techniques in removing the present limitations is evaluated.