SEARCH · Search NASA

Results for “numerical validation”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 217 records · Page 12

All sources NASA NTRS DOE OSTI

NASA NTRS1997

Computations of Complex Three-Dimensional Turbulent Free Jets

Three-dimensional, incompressible turbulent jets with rectangular and elliptical cross-sections are simulated with a finite-difference numerical method. The full Navier- Stokes equations are solved at low Reynolds numbers, whereas at high Reynolds numbers filtered forms of the equations are solved along with a sub-grid scale model to approximate the effects of the unresolved scales. A 2-N storage, third-order Runge-Kutta scheme is used for temporary discretization and a fourth-order compact scheme is used for spatial discretization. Although such methods are widely used in the simulation of compressible flows, the lack of an evolution equation for pressure or density presents particular difficulty in incompressible flows. The pressure-velocity coupling must be established indirectly. It is achieved, in this study, through a Poisson equation which is solved by a compact scheme of the same order of accuracy. The numerical formulation is validated and the dispersion and dissipation errors are documented by the solution of a wide range of benchmark problems. Three-dimensional computations are performed for different inlet conditions which model the naturally developing and forced jets. The experimentally observed phenomenon of axis-switching is captured in the numerical simulation, and it is confirmed through flow visualization that this is based on self-induction of the vorticity field. Statistical quantities such as mean velocity, mean pressure, two-point velocity spatial correlations and Reynolds stresses are presented. Detailed budgets of the mean momentum and Reynolds stresses are presented. Detailed budgets of the mean momentum and Reynolds stress equations are presented to aid in the turbulence modeling of complex jets. Simulations of circular jets are used to quantify the effect of the non-uniform curvature of the non-circular jets.

Wilson, Robert V.↗

NASA NTRS1996

On the Spatial and Temporal Accuracy of Overset Grid Methods for Moving Body Problems

A study of numerical attributes peculiar to an overset grid approach to unsteady aerodynamics prediction is presented. Attention is focused on the effect of spatial error associated with interpolation of intergrid boundary conditions and temporal error associated with explicit update of intergrid boundary points on overall solution accuracy. A set of numerical experiments are used to verify whether or not the use of simple interpolation for intergrid boundary conditions degrades the formal accuracy of a conventional second-order flow solver, and to quantify the error associated with explicit updating of intergrid boundary points. Test conditions correspond to the transonic regime. The validity of the numerical results presented here are established by comparison with existing numerical results of documented accuracy, and by direct comparison with experimental results.

Meakin, Robert L.↗

NASA NTRS2005

Turbulent Vortex-Flow Simulation Over a 65 deg Sharp and Blunt Leading-Edge Delta Wing at Subsonic Speeds

Turbulent thin-layer, Reynolds-Averaged Navier-Stokes solutions, based on a multi-block structured grid, are presented for a 65 deg delta wing having either a sharp leading edge (SLE) or blunt leading edge (BLE) geometry. The primary objective of the study is to assess the prediction capability of the method for simulating the leading-edge flow separation and the ensuing vortex flow characteristics. Computational results are obtained for two angles of attack of approximately 13 and 20 deg, at free-stream Mach number of 0.40 and Reynolds number of 6 million based on the wing mean aerodynamic chord. The effects of two turbulence models of Baldwin-Lomax with Degani-Schiff (BL/DS) and the Spalart-Allmaras (SA) on the numerical results are also discussed. The computations also explore the effects of two numerical flux-splitting schemes, i.e., flux difference splitting (fds) and flux vector splitting (fvs), on the solution development and convergence characteristics. The resulting trends in solution sensitivity to grid resolution for the selected leading-edge geometries, angles of attack, turbulence models and flux splitting schemes are also presented. The validity of the numerical results is evaluated against a unique set of experimental wind-tunnel data that was obtained in the National Transonic Facility at the NASA Langley Research Center.

Ghaffari, Farhad↗

NASA NTRS1977

A numerical study of the unsteady leading edge separation bubble on an oscillating airfoil

A preliminary study of the unsteady viscous flow in the region of an airfoil leading edge was performed, in which the interaction between the viscous and inviscid flow fields is neglected. The solution method uses the finite difference form of the governing equations throughout the separated flow field and incorporates a transition model based on the integral turbulence kinetic energy equation. The validity of the numerical procedure is verified by making comparisons with analytical solutions to several test problems, including unsteady flow over a plate oscillating in its own plane. The method was then applied to the problem of unsteady viscous flow over a NACA 0012 airfoil oscillating sinusoidally in pitch. The flow field characteristics were in qualitative agreement with experimental results. The bubble moved forward on the airfoil and decreased in size as incidence was increased. Viscous flow in the leading edge region was found to be quasi-steady, while bubble height varied inversely with Reynolds number.

Kreskovsky, J. P.↗

NASA NTRS1997

A Study of Fundamental Shock Noise Mechanisms

This paper investigates two mechanisms fundamental to sound generation in shocked flows: shock motion and shock deformation. Shock motion is modeled numerically by examining the interaction of a sound wave with a shock. This numerical approach is validated by comparison with results obtained by linear theory for a small-disturbance case. Analysis of the perturbation energy with Myers' energy corollary demonstrates that acoustic energy is generated by the interaction of acoustic disturbances with shocks. This analysis suggests that shock motion generates acoustic and entropy disturbance energy. Shock deformation is modeled numerically by examining the interaction of a vortex ring with a shock. These numerical simulations demonstrate the generation of both an acoustic wave and contact surfaces. The acoustic wave spreads cylindrically. The sound intensity is highly directional and the sound pressure increases with increasing shock strength. The numerically determined relationship between the sound pressure and the Mach number is found to be consistent with experimental observations of shock noise. This consistency implies that a dominant physical process in the generation of shock noise is modeled in this study.

Meadows, Kristine R.↗

NASA NTRS1994

Backward-facing step measurements at low Reynolds number, Re(sub h)=5000

An experimental study of the flow over a backward-facing step at low Reynolds number was performed for the purpose of validating a direct numerical simulation (DNS) which was performed by the Stanford/NASA Center for Turbulence Research. Previous experimental data on back step flows were conducted at Reynolds numbers and/or expansion ratios which were significantly different from that of the DNS. The geometry of the experiment and the simulation were duplicated precisely, in an effort to perform a rigorous validation of the DNS. The Reynolds number used in the DNS was Re(sub h)=5100 based on step height, h. This was the maximum possible Reynolds number that could be economically simulated. The boundary layer thickness, d, was approximately 1.0 h in the simulation and the expansion ratio was 1.2. The Reynolds number based on the momentum thickness, Re(sub theta), upstream of the step was 610. All of these parameters were matched experimentally. Experimental results are presented in the form of tables, graphs and a floppy disk (for easy access to the data). An LDV instrument was used to measure mean velocity components and three Reynolds stresses components. In addition, surface pressure and skin friction coefficients were measured. LDV measurements were acquired in a measuring domain which included the recirculating flow region.

Jovic, Srba↗

NASA NTRS2015

Toward a Nonlinear Acoustic Analogy: Turbulence as a Source of Sound and Nonlinear Propagation

An acoustic analogy is proposed that directly includes nonlinear propagation effects. We examine the Lighthill acoustic analogy and replace the Green's function of the wave equation with numerical solutions of the generalized Burgers' equation. This is justified mathematically by using similar arguments that are the basis of the solution of the Lighthill acoustic analogy. This approach is superior to alternatives because propagation is accounted for directly from the source to the far-field observer instead of from an arbitrary intermediate point. Validation of a numerical solver for the generalized Burgers' equation is performed by comparing solutions with the Blackstock bridging function and measurement data. Most importantly, the mathematical relationship between the Navier- Stokes equations, the acoustic analogy that describes the source, and canonical nonlinear propagation equations is shown. Example predictions are presented for nonlinear propagation of jet mixing noise at the sideline angle

Miller, Steven A. E.↗

NASA NTRS2015

Toward a Nonlinear Acoustic Analogy: Turbulence as a Source of Sound and Nonlinear Propagation

An acoustic analogy is proposed that directly includes nonlinear propagation effects. We examine the Lighthill acoustic analogy and replace the Green's function of the wave equation with numerical solutions of the generalized Burgers' equation. This is justified mathematically by using similar arguments that are the basis of the solution of the Lighthill acoustic analogy. This approach is superior to alternatives because propagation is accounted for directly from the source to the far-field observer instead of from an arbitrary intermediate point. Validation of a numerical solver for the generalized Burgers' equation is performed by comparing solutions with the Blackstock bridging function and measurement data. Most importantly, the mathematical relationship between the Navier-Stokes equations, the acoustic analogy that describes the source, and canonical nonlinear propagation equations is shown. Example predictions are presented for nonlinear propagation of jet mixing noise at the sideline angle.

Miller, Steven A. E.↗

NASA NTRS2008

High-precision Non-Contact Measurement of Creep of Ultra-High Temperature Materials for Aerospace

For high-temperature applications (greater than 2,000 C) such as solid rocket motors, hypersonic aircraft, nuclear electric/thermal propulsion for spacecraft, and more efficient jet engines, creep becomes one of the most important design factors to be considered. Conventional creep-testing methods, where the specimen and test apparatus are in contact with each other, are limited to temperatures approximately 1,700 C. Development of alloys for higher-temperature applications is limited by the availability of testing methods at temperatures above 2000 C. Development of alloys for applications requiring a long service life at temperatures as low as 1500 C, such as the next generation of jet turbine superalloys, is limited by the difficulty of accelerated testing at temperatures above 1700 C. For these reasons, a new, non-contact creep-measurement technique is needed for higher temperature applications. A new non-contact method for creep measurements of ultra-high-temperature metals and ceramics has been developed and validated. Using the electrostatic levitation (ESL) facility at NASA Marshall Space Flight Center, a spherical sample is rotated quickly enough to cause creep deformation due to centrifugal acceleration. Very accurate measurement of the deformed shape through digital image analysis allows the stress exponent n to be determined very precisely from a single test, rather than from numerous conventional tests. Validation tests on single-crystal niobium spheres showed excellent agreement with conventional tests at 1985 C; however the non-contact method provides much greater precision while using only about 40 milligrams of material. This method is being applied to materials including metals and ceramics for non-eroding throats in solid rockets and next-generation superalloys for turbine engines. Recent advances in the method and the current state of these new measurements will be presented.

Rogers, Jan R.↗

NASA NTRS2008

Analysis of International Space Station Materials on MISSE-3 and MISSE-4

For high-temperature applications (> 2,000 C) such as solid rocket motors, hypersonic aircraft, nuclear electric/thermal propulsion for spacecraft, and more efficient jet engines, creep becomes one of the most important design factors to be considered. Conventional creep-testing methods, where the specimen and test apparatus are in contact with each other, are limited to temperatures ~1,700 deg. C. Development of alloys for higher-temperature applications is limited by the availability of testing methods at temperatures above 2000 C. Development of alloys for applications requiring a long service life at temperatures as low as 1500 C, such as the next generation of jet turbine superalloys, is limited by the difficulty of accelerated testing at temperatures above 1700 0c. For these reasons, a new, non-contact creep-measurement technique is needed for higher temperature applications. A new non-contact method for creep measurements of ultra-high-temperature metals and ceramics has been developed and validated. Using the electrostatic levitation (ESL) facility at NASA Marshall Space Flight Center, a spherical sample is rotated quickly enough to cause creep deformation due to centrifugal acceleration. Very accurate measurement of the deformed shape through digital image analysis allows the stress exponent n to be determined very precisely from a single test, rather than from numerous conventional tests. Validation tests on single-crystal niobium spheres showed excellent agreement with conventional tests at 1985 C; however the non-contact method provides much greater precision while using only about 40 milligrams of material. This method is being applied to materials including metals and ceramics for noneroding throats in solid rockets and next-generation superalloys for turbine engines. Recent advances in the method and the current state of these new measurements will be presented.

Finckenor, Miria M.↗

NASA NTRS1992

Numerical methods for the simulation of complex multi-body flows with applications for the integrated Space Shuttle vehicle

This project forms part of the long term computational effort to simulate the time dependent flow over the integrated Space Shuttle vehicle (orbiter, solid rocket boosters (SRB's), external tank (ET), and attach hardware) during its ascent mode for various nominal and abort flight conditions. Due to the limitations of experimental data such as wind tunnel wall effects and the difficulty of safely obtaining valid flight data, numerical simulations are undertaken to supplement the existing data base. This data can then be used to predict the aerodynamic behavior over a wide range of flight conditions. Existing computational results show relatively good overall comparison with experiments but further refinement is required to reduce numerical errors and to obtain finer agreements over a larger parameter space. One of the important goals of this project is to obtain better comparisons between numerical simulations and experiments. In the simulations performed so far, the geometry has been simplified in various ways to reduce the complexity so that useful results can be obtained in a reasonable time frame due to limitations in computer resources. In this project, the finer details of the major components of the Space Shuttle are modeled better by including more complexity in the geometry definition. Smaller components not included in early Space Shuttle simulations will now be modeled and gridded.

Chan, William M.↗

NASA NTRS2003

Pulsed Inductive Thruster (PIT): Modeling and Validation Using the MACH2 Code

Numerical modeling of the Pulsed Inductive Thruster exercising the magnetohydrodynamics code, MACH2 aims to provide bilateral validation of the thruster's measured performance and the code's capability of capturing the pertinent physical processes. Computed impulse values for helium and argon propellants demonstrate excellent correlation to the experimental data for a range of energy levels and propellant-mass values. The effects of the vacuum tank wall and massinjection scheme were investigated to show trivial changes in the overall performance. An idealized model for these energy levels and propellants deduces that the energy expended to the internal energy modes and plasma dissipation processes is independent of the propellant type, mass, and energy level.

Schneider, Steven↗

NASA NTRS2008

An Experimental and Numerical Study of a Supersonic Burner for CFD Model Development

A laboratory scale supersonic burner has been developed for validation of computational fluid dynamics models. Detailed numerical simulations were performed for the flow inside the combustor, and coupled with finite element thermal analysis to obtain more accurate outflow conditions. A database of nozzle exit profiles for a wide range of conditions of interest was generated to be used as boundary conditions for simulation of the external jet, or for validation of non-intrusive measurement techniques. A set of experiments was performed to validate the numerical results. In particular, temperature measurements obtained by using an infrared camera show that the computed heat transfer was larger than the measured value. Relaminarization in the convergent part of the nozzle was found to be responsible for this discrepancy, and further numerical simulations sustained this conclusion.

Magnotti, G.↗

NASA NTRS1989

Thermal analysis of engine inlet anti-icing systems

A hot air anti-icing system of a gas turbine engine inlet is analyzed numerically. A three-dimensional potential flow code, which accounts for compressibility effects, is used to determine the flowfield in and around the inlet. A particle trajectory code is developed using a local linearization technique. The trajectory code is used to calculate local water impingement rates. Energy balances are performed on both the surface runback water and the metallic skin to determine their temperature distributions. A variety of test cases are considered in order to validate the various numerical components of the process as well as to demonstrate the procedure.

Keith, Theo G., Jr.↗

NASA NTRS2011

Toward Supersonic Retropropulsion CFD Validation

This paper begins the process of verifying and validating computational fluid dynamics (CFD) codes for supersonic retropropulsive flows. Four CFD codes (DPLR, FUN3D, OVERFLOW, and US3D) are used to perform various numerical and physical modeling studies toward the goal of comparing predictions with a wind tunnel experiment specifically designed to support CFD validation. Numerical studies run the gamut in rigor from code-to-code comparisons to observed order-of-accuracy tests. Results indicate that this complex flowfield, involving time-dependent shocks and vortex shedding, design order of accuracy is not clearly evident. Also explored is the extent of physical modeling necessary to predict the salient flowfield features found in high-speed Schlieren images and surface pressure measurements taken during the validation experiment. Physical modeling studies include geometric items such as wind tunnel wall and sting mount interference, as well as turbulence modeling that ranges from a RANS (Reynolds-Averaged Navier-Stokes) 2-equation model to DES (Detached Eddy Simulation) models. These studies indicate that tunnel wall interference is minimal for the cases investigated; model mounting hardware effects are confined to the aft end of the model; and sparse grid resolution and turbulence modeling can damp or entirely dissipate the unsteadiness of this self-excited flow.

Kleb, Bil↗

NASA NTRS2008

Measurement and Prediction of the Thermomechanical Response of Shape Memory Alloy Hybrid Composite Beams

An experimental and numerical investigation into the static and dynamic responses of shape memory alloy hybrid composite (SMAHC) beams is performed to provide quantitative validation of a recently commercialized numerical analysis/design tool for SMAHC structures. The SMAHC beam specimens consist of a composite matrix with embedded pre-strained SMA actuators, which act against the mechanical boundaries of the structure when thermally activated to adaptively stiffen the structure. Numerical results are produced from the numerical model as implemented into the commercial finite element code ABAQUS. A rigorous experimental investigation is undertaken to acquire high fidelity measurements including infrared thermography and projection moire interferometry for full-field temperature and displacement measurements, respectively. High fidelity numerical results are also obtained from the numerical model and include measured parameters, such as geometric imperfection and thermal load. Excellent agreement is achieved between the predicted and measured results of the static and dynamic thermomechanical response, thereby providing quantitative validation of the numerical tool.

Davis, Brian↗

NASA NTRS1979

CO2 and CO electron vibrational cooling rates

Electron cooling rates for electron vibrational excitation of CO2 and CO and rotational excitation of CO are calculated. Results of the calculation are expressed as closed form analytic relations that are both convenient for numerical calculations and valid over a wide range of electron and neutral temperatures. It is found that the cooling rates and their temperature dependences differ significantly from other calculations. For CO2, these differences are related mainly to the assumptions about the approach to equilibrium and not to the actual cross-section choices. For electron temperatures below 1000 K and gas temperatures near 300 K (representative of conditions in the Venus atmosphere, where collisional cooling of the electrons dominates), the CO2 cooling rate is as much as an order of magnitude larger than previously considered in ionospheric models.

Porter, H. S.↗

NASA NTRS1983

Calculation of boundary layers near the stagnation point of an oscillating airfoil

The results of an investigation of boundary layers close to the stagnation point of an oscillating airfoil are reported. Two procedures for generating initial conditions, the characteristics box scheme and a quasi-static approach, were investigated, and the quasi-static approach was shown to be appropriate provided the initial region was far from any flow separation. With initial conditions generated in this way, the unsteady boundary layer equations were solved for the flow in the leading edge region of a NACA 0012 airfoil oscillating from 0 to 5 deg. Results were obtained for both laminar and turbulent flow, and, in the latter case, the effect of transition was assessed by specifying its occurrence at different locations. The results demonstrate the validity of the numerical scheme and suggest that the procedures should be applied to calculation of the entire flow around oscillating airfoils.

Cebeci, T.↗