Search NASA⌕ Search

SEARCH · Search NASA

Results for “Computational fluid dynamics code”

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.

At least 1,315 records · Page 73
All sourcesNASA NTRSDOE OSTI
NASA NTRS2000

Functionality Improvements to Overaero

The functionality of the overset, static aeroelasticity, Navier-Stokes flow solver OVERAERO was increased by adding capability to the flow solver and enhancing code performance. Improvements were made to the fluids/structure interface, an MLP version of the parallel OVERAERO code was developed, and the OVERAERO-MPI code was ported to the Cray T3E. The OVERFLOW-MPI and OVERAERO-MPI codes were tested successfully on the IPG testbed and a means of reducing communication overhead within OVERFLOW-MPI was investigated. To solve an aeroelastic problem computationally, a structures grid surface definition and a fluids grid surface definition are required. Typically, the structures grid surface has a lower fidelity than the fluids grid surface. Thus, the methods developed to transfer data between the two grid systems are vital to the accuracy and efficiency of the aeroelasticity code. The fluids/structures interface developed for the OVERAERO code was improved to more accurately treat fluids surfaces that bridge between two different structural surfaces. For example, the method allowed the forward portion of a flap track fairing to deform with the wing and the aft end of the fairing to deform with the flap. A tightly-coupled version of the code based on OVERFLOW-MLP was developed to improve code performance on the SGI Origin 2000. This required a new parallelization strategy to couple the fluids and structures codes. The OVERAERO-MPI code was ported to the Cray T3E to extend the usability of the code. The port required extensive use of dynamic memory management techniques to fit large problems within the memory limitations of the T3E. The OVERFLOW-MPI and OVERAERO-MPI codes were tested on the IPG testbed being developed within NASA. For small problems with minimal data transfer between grids, there was little to no performance penalty spreading the computation across two machines. For very large problems, methods were developed to minimize intermachine communication via the grid partitioning scheme. By minimizing the intermachine communication requirements of the problem, it may still be beneficial to run a tightly-coupled flow solver across two machines within the IPG.

Gee, Ken↗
DOE OSTI2020

Coupling of Laminar-Turbulent Transition with RANS Computational Fluid Dynamics

This project combines several new concepts to create a boundary layer transition prediction capability that is suitable for analyzing modern hypersonic flight vehicles. The first new concept is the use of ''optimization'' methods to detect the hydrodynamic instabilities that cause boundary layer transition; the use of this method removes the need for many limiting assumptions of other methods and enables quantification of the interactions between boundary layer instabilities and the flow field imperfections that generate them. The second new concept is the execution of transition analysis within a conventional hypersonics CFD code, using the same mesh and numerical schemes for the transition analysis and the laminar flow simulation. This feature enables rapid execution of transition analysis with less user oversight required and no interpolation steps needed.

42 ENGINEERING↗
NASA NTRS1995

CFD-Based Design Optimization Tool Developed for Subsonic Inlet

The traditional approach to the design of engine inlets for commercial transport aircraft is a tedious process that ends with a less-than-optimum design. With the advent of high-speed computers and the availability of more accurate and reliable computational fluid dynamics (CFD) solvers, numerical optimization processes can effectively be used to design an aerodynamic inlet lip that enhances engine performance. The designers' experience at Boeing Corporation showed that for a peak Mach number on the inlet surface beyond some upper limit, the performance of the engine degrades excessively. Thus, our objective was to optimize efficiency (minimize the peak Mach number) at maximum cruise without compromising performance at other operating conditions. Using a CFD code NPARC, the NASA Lewis Research Center, in collaboration with Boeing, developed an integrated procedure at Lewis to find the optimum shape of a subsonic inlet lip and a numerical optimization code, ADS. We used a GRAPE-based three-dimensional grid generator to help automate the optimization procedure. The inlet lip shape at the crown and the keel was described as a superellipse, and the superellipse exponents and radii ratios were considered as design variables. Three operating conditions: cruise, takeoff, and rolling takeoff, were considered in this study. Three-dimensional Euler computations were carried out to obtain the flow field. At the initial design, the peak Mach numbers for maximum cruise, takeoff, and rolling takeoff conditions were 0.88, 1.772, and 1.61, respectively. The acceptable upper limits on the takeoff and rolling takeoff Mach numbers were 1.55 and 1.45. Since the initial design provided by Boeing was found to be optimum with respect to the maximum cruise condition, the sum of the peak Mach numbers at takeoff and rolling takeoff were minimized in the current study while the maximum cruise Mach number was constrained to be close to that at the existing design. With this objective, the optimum design satisfied the upper limits at takeoff and rolling takeoff while retaining the desirable cruise performance. Further studies are being conducted to include static and cross-wind operating conditions in the design optimization procedure. This work was carried out in collaboration with Dr. E.S. Reddy of NYMA, Inc.

Source record↗
DOE OSTI2023

Improved analytic modeling of the linear Rayleigh–Taylor instability with plasma transport

The Rayleigh–Taylor (R–T) instability with transport effects has long been studied analytically in the neutral fluid regime. How this picture may differ in plasmas has been given less consideration, however. Leveraging previous numerical simulations with plasma viscosity and diffusion of R–T at a binary plasma interface, we show how plasma transport and kinetic effects alter the linear R–T dispersion relation. Our insight informs a much improved analytical dispersion relation, which better matches the simulation data than previously established models. Additionally, we demonstrate the importance of employing asymptotically correct viscosity coefficients in fluid codes for accurately capturing plasma R–T growth.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2016

Wind-US Users Guide Version 3.0

Wind-US is a computational platform which may be used to numerically solve various sets of equations governing physical phenomena. Currently, the code supports the solution of the Euler and Navier-Stokes equations of fluid mechanics, along with supporting equation sets governing turbulent and chemically reacting flows. Wind-US is a product of the NPARC Alliance, a partnership between the NASA Glenn Research Center (GRC) and the Arnold Engineering Development Complex (AEDC) dedicated to the establishment of a national, applications-oriented flow simulation capability. The Boeing Company has also been closely associated with the Alliance since its inception, and represents the interests of the NPARC User's Association. The "Wind-US User's Guide" describes the operation and use of Wind-US, including: a basic tutorial; the physical and numerical models that are used; the boundary conditions; monitoring convergence; the files that are read and/or written; parallel execution; and a complete list of input keywords and test options. For current information about Wind-US and the NPARC Alliance, please see the Wind-US home page at http://www.grc.nasa.gov/WWW/winddocs/ and the NPARC Alliance home page at http://www.grc.nasa.gov/WWW/wind/. This manual describes the operation and use of Wind-US, a computational platform which may be used to numerically solve various sets of equations governing physical phenomena. Wind-US represents a merger of the capabilities of four CFD codes - NASTD (a structured grid flow solver developed at McDonnell Douglas, now part of Boeing), NPARC (the original NPARC Alliance structured grid flow solver), NXAIR (an AEDC structured grid code used primarily for store separation analysis), and ICAT (an unstructured grid flow solver developed at the Rockwell Science Center and Boeing).

Navier-Stokes equation↗
DOE OSTI2020

Free-Surface Liquid Lithium Flow Modeling and Stability Analysis for Fusion Applications

Liquid metal plasma facing components are considered an attractive design choice for fusion devices including pilot plants. Virtual prototyping of such devices includes modeling of free-surface flow of the electrically conductive liquid, which requires computational fluid dynamics (CFD) and magnetohydrodynamics (MHD) simulations. Numerical tools capable of simulating flows and heat transfer in the free-surface MHD flow were developed at PPPL based on the customized ANSYS CFX. Here, MHD is introduced using a magnetic vector potential approach. Free-surface flow capabilities are available in the code and were tested. Special stabilization procedures were derived and applied to improve convergence of the momentum equations with the source terms due to the Lorentz force and surface tension. Important characteristics of the fusion-relevant liquid metal flow is free surface smoothness and stability. Heat flux from the plasma impacts the liquid surface at a very acute angle, so any change of the free surface from axisymmetry can dramatically increase the local heat flux density and thus create excessive evaporation of liquid lithium into the plasma, which is detrimental to operations. Stability analysis of the liquid metal film flow was performed to determine applicable flow regimes. Thin film flow along horizontal wall is considered. Effects of gravity, magnetic field, and surface tension are included in the analysis.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Challenges in Rotor Aerodynamic Modeling for Non-Uniform Inflow Conditions

Within an international collaboration framework, the accuracy of rotor aerodynamic models used for design load calculations of wind turbines is being assessed. Where the use of high-fidelity computation fluid dynamics (CFD) and mid-fidelity free-vortex wake (FVW) models has become commonplace within the wind energy community, these still fail to meet the requirements in terms of execution time and computational cost needed for design load calculations. The fast but engineering fidelity blade-element/momentum (BEM) method can therefore still be considered the industry workhorse for design load simulations. At the same time, upscaling of wind turbine rotors makes inflow non-uniformities (e.g. shear, veer, turbulence) more important. The objective of this work is to assess model accuracy in non-uniform inflow conditions, which violate several BEM assumptions. Thereto a comparison in turbulent inflow has been executed including a wide variety of codes, focusing on the DanAero field measurements, where a 2.3-MW turbine was equipped with, among other sensors, a pressure measurement apparatus. The results indicate that, although average load patterns are in good agreement, this does not hold for the unsteady loads that drive fatigue damage and aero-elastic stability. A simplified comparison round in vertical shear was initiated to investigate the observed differences in a more controlled manner. A consistent offset in load amplitude was observed between CFD and free-vortex codes on the one hand and BEM-type codes on the other hand. To shed more light on the observations, dedicated efforts are ongoing to pinpoint the cause for these differences, in the end leading to guidelines for an improved BEM implementation.

17 WIND ENERGY↗
NASA NTRS2023

Assessment of the Fluid Dynamics Boundary Condition in Ablating or Blowing Flows

Improved models of ablative thermal protection systems have enabled the treatment of materials and fluid behavior in a coupled manner. This paper reports a new approach to modeling the interface between fluid and material, with attention to the conservation of species mass flux and energy on the fluid side of the interface. The general equation is presented and is shown to recover the traditional uncoupled fluid/materials response interface. Including the chemical reaction terms on the CFD side of the interface makes the heat flux exchange independent of the thermodynamic reference state and, therefore, a measurable quantity. Doing so allows the material response solver to take as input the surface heat flux rather than a film coefficient. Removing the film coefficient approximation enables more direct solution of vehicle thermal response but requires consistency in the wall state. The mixing of the shock layer and pyrolysis gas is then computed with finite rate chemistry within the fluid solver. The boundary conditions described have been implemented in the DPLR v4.05.1 code. Char removal is captured using finite rate chemistry in DPLR’s gas surface interaction module. Aspects of coupling these solutions to material response are discussed.

Ablation↗
NASA NTRS2023

Assessment of the Fluid Dynamics Boundary Condition in Ablating or Blowing Flows

Improved models of ablative thermal protection systems have enabled the treatment of materials and fluid behavior in a coupled manner. This paper reports a new approach to modeling the interface between fluid and material, with attention to the conservation of species mass flux and energy on the fluid side of the interface. The general equation is presented and is shown to recover the traditional uncoupled fluid/materials response interface. Including the chemical reaction terms on the CFD side of the interface makes the heat flux exchange independent of the thermodynamic reference state and, therefore, a measurable quantity. Doing so allows the material response solver to take as input the surface heat flux rather than a film coefficient. Removing the film coefficient approximation enables more direct solution of vehicle thermal response but requires consistency in the wall state. The mixing of the shock layer and pyrolysis gas is then computed with finite rate chemistry within the fluid solver. The boundary conditions described have been implemented in the DPLR v4.05.1 code. Char removal is captured using finite rate chemistry in DPLR’s gas surface interaction module. Aspects of coupling these solutions to material response are discussed.

Ablation↗
DOE OSTI2022

Numerical Investigation of the Spark Discharge Process in a Crossflow

The present study numerically investigates the spark discharge process under crossflow conditions using a thermal equilibrium plasma solver that fully couples the electromagnetic physics and fluid dynamics in a computational framework. Numerical results are validated by the comparison with experimental data. Here, the spark discharge experiment is performed in a constant volume vessel using an inductive coil ignition system for automotive applications, and the evolution of the spark channel is measured using high-speed imaging. The crossflow in the gap between the spark-plug electrodes is generated by a rotating fan with two different fan speeds, and the flow velocity across the gap is characterized by particle image velocimetry (PIV) measurement. A computational fluid dynamics (CFD) solver is employed to simulate the crossflow and provide the flow field variables (velocity, pressure, temperature) to the plasma solver. The crossflow velocity predicted in the flow simulation agrees well with the PIV data in that the non-uniform velocity profiles at monitoring points are reproduced by the CFD code. With the crossflow initialization in the plasma solver, the simulated spark discharge process from the breakdown to spark discharge matches the experimental data, including the voltage and circuit waveforms and the high-speed images of the spark channel evolution. The stretch of spark channel captured by plasma simulations agrees with the measured data. The plasma simulation reveals that the mean temperature of the spark channel is maintained at 5000 K during the discharge phase, and the temperature varies along the spark channel so that the highest value is obtained at the spark root on the center electrode. Overall, the results presented in this paper are meant to provide valuable information about the properties of the plasma generated by the spark discharge.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

3D Simulations and MLT. I. Renzini’s Critique

Renzini wrote an influential critique of "overshooting" in mixing-length theory (MLT), as used in stellar evolution codes, and concluded that three-dimensional fluid dynamical simulations were needed. Such simulations are now well tested. Implicit large eddy simulations connect large-scale stellar flow to a turbulent cascade at the grid scale, and allow the simulation of turbulent boundary layers, with essentially no assumptions regarding flow except the number of computational cells. Buoyant driving balances turbulent dissipation for weak stratification, as in MLT, but with the dissipation length replacing the mixing length. The turbulent kinetic energy in our computational domain shows steady pulses after 30 turnovers, with no discernible diminution; these are caused by the necessary lag in turbulent dissipation behind acceleration. Interactions between coherent turbulent structures give multi-modal behavior, which drives intermittency and fluctuations. These cause mixing, which may justify use of the instability criterion of Schwarzschild rather than the Ledoux. Chaotic shear flow of turning material at convective boundaries causes instabilities that generate waves and sculpt the composition gradients and boundary layer structures. The flow is not anelastic; wave generation is necessary at boundaries. A self-consistent approach to boundary layers can remove the need for ad hoc procedures of "convective overshooting" and "semi-convection." In Paper II, we quantify the adequacy of our numerical resolution in a novel way, determine the length scale of dissipation—the "mixing length"—without astronomical calibration, quantify agreement with the four-fifths law of Kolmogorov for weak stratification, and deal with strong stratification.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1982

Experiences with explicit finite-difference schemes for complex fluid dynamics problems on STAR-100 and CYBER-203 computers

Several two- and three-dimensional external and internal flow problems solved on the STAR-100 and CYBER-203 vector processing computers are described. The flow field was described by the full Navier-Stokes equations which were then solved by explicit finite-difference algorithms. Problem results and computer system requirements are presented. Program organization and data base structure for three-dimensional computer codes which will eliminate or improve on page faulting, are discussed. Storage requirements for three-dimensional codes are reduced by calculating transformation metric data in each step. As a result, in-core grid points were increased in number by 50% to 150,000, with a 10% execution time increase. An assessment of current and future machine requirements shows that even on the CYBER-205 computer only a few problems can be solved realistically. Estimates reveal that the present situation is more storage limited than compute rate limited, but advancements in both storage and speed are essential to realistically calculate three-dimensional flow.

Kumar, A.↗
NASA NTRS1983

Improved method for transonic airfoil design-by-optimization

An improved method for use of optimization techniques in transonic airfoil design is demonstrated. FLO6QNM incorporates a modified quasi-Newton optimization package, and is shown to be more reliable and efficient than the method developed previously at NASA-Ames, which used the COPES/CONMIN optimization program. The design codes are compared on a series of test cases with known solutions, and the effects of problem scaling, proximity of initial point to solution, and objective function precision are studied. In contrast to the older method, well-converged solutions are shown to be attainable in the context of engineering design using computational fluid dynamics tools, a new result. The improvements are due to better performance by the optimization routine and to the use of problem-adaptive finite difference step sizes for gradient evaluation.

Kennelly, R. A., Jr.↗
NASA NTRS1983

Improved method for transonic airfoil design-by-optimization

An improved method for use of optimization techniques in transonic airfoil design is demonstrated. FLO6QNM incorporates a modified quasi-Newton optimization package, and is shown to be more reliable and efficient than the method developed previously at NASA-Ames, which used the COPES/CONMIN optimization problem. The design codes are compared on a series of test cases with known solutions, and the effects of problem scaling, proximity of initial point to solution, and objective function precision are studied. In contrast to the older method, well-converged solutions are shown to be attainable in the context of engineering design using computational fluid dynamics tools, a new result. The improvements are due to better performance by the optimization routine and to the use of problem-adaptive finite difference step sizes for gradient evaluation.

Kennelly, R. A., Jr.↗
NASA NTRS1998

Reduction and Analysis of Phosphor Thermography Data With the IHEAT Software Package

Detailed aeroheating information is critical to the successful design of a thermal protection system (TPS) for an aerospace vehicle. This report describes NASA Langley Research Center's (LaRC) two-color relative-intensity phosphor thermography method and the IHEAT software package which is used for the efficient data reduction and analysis of the phosphor image data. Development of theory is provided for a new weighted two-color relative-intensity fluorescence theory for quantitatively determining surface temperatures on hypersonic wind tunnel models; an improved application of the one-dimensional conduction theory for use in determining global heating mappings; and extrapolation of wind tunnel data to flight surface temperatures. The phosphor methodology at LaRC is presented including descriptions of phosphor model fabrication, test facilities and phosphor video acquisition systems. A discussion of the calibration procedures, data reduction and data analysis is given. Estimates of the total uncertainties (with a 95% confidence level) associated with the phosphor technique are shown to be approximately 8 to 10 percent in the Langley's 31-Inch Mach 10 Tunnel and 7 to 10 percent in the 20-Inch Mach 6 Tunnel. A comparison with thin-film measurements using two-inch radius hemispheres shows the phosphor data to be within 7 percent of thin-film measurements and to agree even better with predictions via a LATCH computational fluid dynamics solution (CFD). Good agreement between phosphor data and LAURA CFD computations on the forebody of a vertical takeoff/vertical lander configuration at four angles of attack is also shown. In addition, a comparison is given between Mach 6 phosphor data and laminar and turbulent solutions generated using the LAURA, GASP and LATCH CFD codes. Finally, the extrapolation method developed in this report is applied to the X-34 configuration with good agreement between the phosphor extrapolation and LAURA flight surface temperature predictions. The phosphor process outlined in the paper is believed to provide the aerothermodynamic community with a valuable capability for rapidly obtaining (4 to 5 weeks) detailed heating information needed in TPS design.

Merski, N. Ronald↗
NASA NTRS1991

Computation of supersonic jet mixing noise for an axisymmetric CD nozzle using k-epsilon turbulence model

The turbulent mixing noise of a supersonic jet is calculated for a round convergent-divergent nozzle at the design pressure ratio. Aerodynamic computations are performed using the PARC code with a k-epsilon turbulence model. Lighthill's acoustic analogy combined with Ribner's assumption is adopted. The acoustics solution is based upon the methodology followed by GE in the MGB code. The source correlation function is expressed as a linear combination of second-order tensors. Assuming separable second-order correlations and incorporating Batchelor's isotropic turbulence model, the source term was calculated from the kinetic energy of turbulence. A Gaussian distribution for the time-delay of correlation was introduced. The computational fluid dynamics (CFD) solution was used to obtain the source strength as well as the characteristic time-delay of correlation. The effect of sound/flow interaction was incorporated using the high frequency asymptotic solution to Lilley's equation for axisymmetric geometries. Acoustic results include sound pressure level directivity and spectra at different polar angles. The aerodynamic and acoustic results demonstrate favorable agreement with experimental data.

Khavaran, Abbas↗
NASA NTRS1992

Computation of supersonic jet mixing noise for an axisymmetric CD nozzle using k-epsilon turbulence model

The turbulent mixing noise of a supersonic jet is calculated for a round convergent-divergent nozzle at the design pressure ratio. Aerodynamic computations are performed using the PARC code with a k-epsilon turbulence model. Lighthill's acoustic analogy combined with Ribner's assumption is adopted. The acoustics solution is based upon the methodology followed by GE in the MGB code. The source correlation function is expressed as a linear combination of second-order tensors. Assuming separable second-order correlations and incorporating Batchelor's isotropic turbulence model, the source term was calculated from the kinetic energy of turbulence. A Gaussian distribution for the time-delay of correlation was introduced. The computational fluid dynamics (CFD) solution was used to obtain the source strength as well as the characteristic time-delay of correlation. The effect of sound/flow interaction was incorporated using the high frequency asymptotic solution to Lilley's equation for axisymmetric geometries. Acoustic results include sound pressure level directivity and spectra at different polar angles. The aerodynamic and acoustic results demonstrate favorable agreement with experimental data.

Khavaran, A.↗
NASA NTRS1995

Computational aeroacoustics on massively parallel computers

A computational fluid dynamics (CFD) technique is employed to solve aeroacoustics problems on massively parallel computers. The algorithm is based on a 4th order accurate central finite difference and 4th order accurate 4-stage Runge-Kutta time integration method that solves the 3-D full Navier-Stokes/Euler equations on a curvilinear coordinate system. The code has been developed to predict noise radiation from ducted fans. A specific application of this code is made here to the oscillating circular piston problem (Workshop Category 4, Problem 2) by solving the full Euler equations with nonreflecting boundary conditions in the axisymmetric mode. The features of the code that apply to this particular problem are described in this paper.

Oezyoeruek, Yusuf↗