Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “CFD modeling”

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 775 records · Page 43
All sourcesNASA NTRSDOE OSTI
NASA NTRS1996

Analysis and Modeling of a Two-Phase Jet Pump of a Flow Boiling Test Facility for Aerospace Applications

Jet pumps are devices capable of pumping fluids to a higher pressure employing a nozzle/diffuser/mixing chamber combination. A primary fluid is usually allowed to pass through a converging-diverging nozzle where it can accelerate to supersonic speeds at the nozzle exit. The relatively high kinetic energy that the primary fluid possesses at the nozzle exit is accompanied by a low pressure region in order to satisfy Bernoulli's equation. The low pressure region downstream of the nozzle exit permits a secondary fluid to be entrained into and mixed with the primary fluid in a mixing chamber located downstream of the nozzle. Several combinations may exist in terms of the nature of the primary and secondary fluids in so far as whether they are single or two-phase fluids. Depending on this, the jet pump may be classified as gas/gas, gas/liquid, liquid/liquid, two-phase/liquid, or similar combinations. The mixing chamber serves to create a homogeneous single-phase or two-phase mixture which enters a diffuser where the high kinetic energy of the fluid is converted into pressure energy. If the fluid mixture entering the diffuser is in the supersonic flow regime, a normal shock wave usually develops inside the diffuser. If the fluid mixture is one that can easily change phase, a condensation shock would normally develop. Because of the overall rise in pressure in the diffuser as well as the additional rise in pressure across the shock layer, condensation becomes more likely. Associated with the pressure rise across the shock is a velocity reduction from the supersonic to the subsonic range. If the two-phase flow entering the diffuser is predominantly gaseous with liquid droplets suspended in it, it will transform into a predominantly liquid flow containing gaseous bubbles (bubbly flow) somewhere in the diffuser. While past researchers have been able to model the two-phase flow jet pump using the one-dimensional assumption with no shock waves and no phase change, there is no research known to the authors apart from that of Anand (1992) which accounted for condensation shocks. One of the objectives of this research effort is to develop a comprehensive model in which the effects of phase slip and inter-phase heat transfer as well as the wall friction and shock waves are accounted for. While this modeling effort is predominantly analytical in nature and is primarily intended to provide a parametric understanding of the jet pump performance under different operating scenarios, another parallel effort employing a commercial CFD code is also implemented. The latter effort is primarily intended to model an axisymmetric counterpart of the problem in question. The viability of using the CFD code to model a two-phase flow jet pump will be assessed by attempting to recreate some of the existing performance data of similar jet pumps. The code will eventually be used to generate the jet pump performance characteristics of several scenarios involving jet pump geometries as well as flow regimes in order to be able to determine an optimum design which would be suitable for a two-phase flow boiling test facility at NASA-Marshall. Because of the extensive nature of the analytical model developed, the following section will only provide very brief highlights of it, while leaving the details to a more complete report submitted to the NASA colleague. This report will also contain some of the simulation results obtained using the CFD code.

Sherif, S. A.↗
NASA NTRSDate not supplied

Progress on Inverse Estimation Technique of Non-Linear Pitch Damping Coefficient Curves Using Free-Flight CFD Generated Trajectories

Characterization of entry vehicle pitch damping coefficient curves is crucial to ensure appropriate re-entry and overall mission success. The pitch damping coefficient (C_(m_q )+C_(m_α ̇ )) is used to encapsulate the oscillatory growth or decay of a body during a trajectory. The inverse estimation technique utilizes an existing Free-Flight CFD (FF-CFD) dataset and wraps a reconstruction algorithm in an optimizer. The reconstruction integrates the planar equations of motion derived by Schoenenberger, Queen [1] using Python’s scipy.integrate.solve_ivp. The optimizer’s objective function is the normalized 𝐿2 residual of the angle of attack peaks between the reconstructed trajectory and the original data produced with FF-CFD. Inclusion of the peak times in this residual calculation allows for simultaneous optimization of the pitch moment coefficient, C_(m_α ). This residual equation is shown below in Eq. 1. The optimizer scipy.optimize.minimize was used with the gradient-based Powell method for the analysis presented, however the differential evolution method was investigated as means of comparison, and was found to produce marginally lower residual values with prohibitively longer run times. Further, the pitch damping curve is found by fitting a cubic interpolation function to a set of (α, (C_(m_q )+C_(m_α ̇ ))) control points, where the α points are held constant and the (C_(m_q )+C_(m_α ̇ )) values are the optimized parameters. The pitch moment curve uses a linear interpolation between the minimum and maximum α in the dataset. FF-CFD generated trajectories using the Dragonfly capsule geometry with the Genesis ballistic range model parameters were simulated and used for this analysis. These FF-CFD trajectories simulate planar motion, as restricted by the reconstructing the equations of motion, of three different cases: 1-DoF (free-to-pitch), 2-DoF (free-to-pitch and heave), and 3-DoF (free-to-pitch, heave, and decelerate). Pitch damping coefficient curves generated using this inverse estimation curve technique with FF-CFD 1-DoF Dragonfly data are found in Fig. 1. Preliminary results reconstructing ballistic range shots using these FF-CFD derived predictions of the pitch damping curve (Fig. 1) are shown in Fig. 2. It should be noted that the ballistic range shot used a Genesis model whereas the FF-CFD data used a Dragonfly geometry, however these geometries are similar.

entry↗
NASA NTRSDate not supplied

Rapid Aero Modeling of a Lift+Cruise UAM Configuration for Stability and Control Using Overset Grid CFD

This talk will address some aspects of an effort to develop both static and dynamic stability and control derivatives of a generic Lift+Cruise Urban Air Mobility vehicle. In particular, the use of NASA overset grid CFD packages OVERFLOW and Chimera Grid Tools allowed for rapid grid and solution generation while varying many database factors. Approach details include coarse gridding and control effectors, rotor disk modeling, convergence evaluation, and parallel execution on a midrange compute resource. A brief *description of the Rapid Aero Modeling (RAM) approach, including design of experiments and response surface modeling, will be given, and recent simulations of powered and unpowered dynamic multisine maneuvers will be presented.

Overset Grid CFD↗
NASA NTRS1987

Wind tunnel requirements for computational fluid dynamics code verification

The role of experiment in the development of Computational Fluid Dynamics (CFD) for aerodynamic flow field prediction is discussed. Requirements for code verification from two sources that pace the development of CFD are described for: (1) development of adequate flow modeling, and (2) establishment of confidence in the use of CFD to predict complex flows. The types of data needed and their accuracy differs in detail and scope and leads to definite wind tunnel requirements. Examples of testing to assess and develop turbulence models, and to verify code development, are used to establish future wind tunnel testing requirements. Versatility, appropriate scale and speed range, accessibility for nonintrusive instrumentation, computerized data systems, and dedicated use for verification were among the more important requirements identified.

Marvin, Joseph G.↗
NASA NTRS1988

Wind tunnel requirements for computational fluid dynamics code verification

The role of experiment in the development of computational fluid dynamics (CFD) for aerodynamic flow field prediction is discussed. Requirements for code verification from two sources that pace the development of CFD are described for: (1) development of adequate flow modeling, and (2) establishment of confidence in the use of CFD to predict complex flows. The types of data needed and their accuracy differs in detail and scope and leads to definite wind tunnel requirements. Examples of testing to assess and develop turbulence models, and to verify code development, are used to establish future wind tunnel testing requirements. Versatility, appropriate scale and speed range, accessibility for nonintrusive instrumentation, computerized data systems, and dedicated use for verification were among the more important requirements identified.

Marvin, Joseph G.↗
NASA NTRS2011

An Overview of Preliminary Computational and Experimental Results for the Semi-Span Super-Sonic Transport (S4T) Wind-Tunnel Model

A summary of computational and experimental aeroelastic and aeroservoelastic (ASE) results for the Semi-Span Super-Sonic Transport (S4T) wind-tunnel model is presented. A broad range of analyses and multiple ASE wind-tunnel tests of the S4T have been performed in support of the ASE element in the Supersonics Program, part of NASA's Fundamental Aeronautics Program. The computational results to be presented include linear aeroelastic and ASE analyses, nonlinear aeroelastic analyses using an aeroelastic CFD code, and rapid aeroelastic analyses using CFD-based reduced-order models (ROMs). Experimental results from two closed-loop wind-tunnel tests performed at NASA Langley's Transonic Dynamics Tunnel (TDT) will be presented as well.

Silva, Walter A.↗
NASA NTRS1993

Comparison of liquid rocket engine base region heat flux computations using three turbulence models

The flow in the base region of launch vehicles is characterized by flow separation, flow reversals, and reattachment. Computation of the convective heat flux in the base region and on the nozzle external surface of Space Shuttle Main Engine and Space Transportation Main Engine (STME) is an important part of defining base region thermal environments. Several turbulence models were incorporated in a CFD code and validated for flow and heat transfer computations in the separated and reattaching regions associated with subsonic and supersonic flows over backward facing steps. Heat flux computations in the base region of a single STME engine and a single S1C engine were performed using three different wall functions as well as a renormalization-group based k-epsilon model. With the very limited data available, the computed values are seen to be of the right order of magnitude. Based on the validation comparisons, it is concluded that all the turbulence models studied have predicted the reattachment location and the velocity profiles at various axial stations downstream of the step very well.

Kumar, Ganesh N.↗
NASA NTRS2011

Gas Generator Feedline Orifice Sizing Methodology: Effects of Unsteadiness and Non-Axisymmetric Flow

Engine LH2 and LO2 gas generator feed assemblies were modeled with computational fluid dynamics (CFD) methods at 100% rated power level, using on-center square- and round-edge orifices. The purpose of the orifices is to regulate the flow of fuel and oxidizer to the gas generator, enabling optimal power supply to the turbine and pump assemblies. The unsteady Reynolds-Averaged Navier-Stokes equations were solved on unstructured grids at second-order spatial and temporal accuracy. The LO2 model was validated against published experimental data and semi-empirical relationships for thin-plate orifices over a range of Reynolds numbers. Predictions for the LO2 square- and round-edge orifices precisely match experiment and semi-empirical formulas, despite complex feedline geometry whereby a portion of the flow from the engine main feedlines travels at a right-angle through a smaller-diameter pipe containing the orifice. Predictions for LH2 square- and round-edge orifice designs match experiment and semi-empirical formulas to varying degrees depending on the semi-empirical formula being evaluated. LO2 mass flow rate through the square-edge orifice is predicted to be 25 percent less than the flow rate budgeted in the original engine balance, which was subsequently modified. LH2 mass flow rate through the square-edge orifice is predicted to be 5 percent greater than the flow rate budgeted in the engine balance. Since CFD predictions for LO2 and LH2 square-edge orifice pressure loss coefficients, K, both agree with published data, the equation for K has been used to define a procedure for orifice sizing.

Rothermel, Jeffry↗
NASA NTRS2008

Rolls-Royce Low Noise Highly Variable Cycle Nozzle for Next Generation Supersonic Aircraft

An overview of the work performed by Rolls-Royce under contract NNL08AA29C is presented. The work includes computational fluid dynamic (CFD) analysis for, and design of, a highly variable cycle exhaust model for the Supersonic project (NRA NN06ZEA001N). The CFD analysis shows that the latest design improvements to the clam shell doors have increased flow through the ejector over that achieved with previous designs.

Sokhey, Jack S.↗
NASA NTRS2012

Quasi One-Dimensional Unsteady Modeling of External Compression Supersonic Inlets

The AeroServoElasticity task under the NASA Supersonics Project is developing dynamic models of the propulsion system and the vehicle in order to conduct research for integrated vehicle dynamic performance. As part of this effort, a nonlinear quasi 1-dimensional model of an axisymmetric external compression supersonic inlet is being developed. The model utilizes compressible flow computational fluid dynamics to model the internal inlet segment as well as the external inlet portion between the cowl lip and normal shock, and compressible flow relations with flow propagation delay to model the oblique shocks upstream of the normal shock. The external compression portion between the cowl-lip and the normal shock is also modeled with leaking fluxes crossing the sonic boundary, with a moving CFD domain at the normal shock boundary. This model has been verified in steady state against tunnel inlet test data and it s a first attempt towards developing a more comprehensive model for inlet dynamics.

Kopasakis, George↗
NASA NTRS2018

Helium Pressurization Modeling of the EDU Submerged Diffuser Using Thermal Desktop

NASA is currently working on developing its predictive modeling capabilities to support future missions that require long term storage of cryogenic propellants. Through the Evolvable Cryogenics project, Development and Validation of Analysis Tools (DVAT) task, both CFD and multi-node fluid modeling tools are being validated against to increase the following cryogenic fluid management technologies: Self-Pressurization, Pressure Control, Helium Pressurization, de-stratification, jet mixing, Tank and transfer line chill-down, tank filling and transfer. NASA is using previous and future micro-gravity experiments to validate both types of models. Recent activities include: Validation of the Tank Pressure Control Experiment (TPCE), Hydrogen and Helium pressurization, and Completion of the Zero Boil Off Tank (ZBOT) experiment on ISS. Future activities include validation of the ZBOT experiment with both CFD and multi-node analysis tools, validation of the Cryogenic Demonstration System of the Robotic Refueling Mission (RRM3).

low-gravity fluid modeling↗
NASA NTRSDate not supplied

Multiobjective Multidisciplinary Optimization of Low-Boom Supersonic Transports Using Multifidelity Models

A multidisciplinary optimization (MDO) method has been developed to design a computational fluid dynamics (CFD) based low-boom configuration that can be obtained from a Pareto solution of a low-fidelity multiobjective MDO problem with mission constraints. This paper refines the developed MDO method using multifidelity models for CFD-based multiobjective MDO. The refined MDO method can generate a low-boom configuration that satisfies the mission requirements, has the lowest takeoff gross weight and the longest range for the low-boom mission, trims the low-boom cruise flight with fuel redistributions, and has a reversed equivalent area distribution closely matching a low-boom target with ground noise level below 70 PLdB. The validity of the refined MDO method is demonstrated by a design study of a low-boom supersonic transport that carries 40 passengers, flies a low-boom mission with cruise Mach of 1.7 and range of 3500 nm, and cruises overwater at Mach 1.8 with range of 3882 nm. Moreover, the refined MDO method eliminates the difference between the assumed cruise weight for CFD-based low-boom inverse design optimization and the estimated cruise weight of the optimal inverse design solution with respect to the mission requirements.

multidisciplinary optimization↗
NASA NTRSDate not supplied

CFD Validation Study of a Hypersonic Cone-Slice-Flap Variable Geometry Configuration

Model validation is the process of determining the degree of accuracy between physical reality and the model. The result of model validation can either be used to improve the model through calibration or quantify the model-form uncertainty. This work focuses on providing the model-form uncertainty through an area metric for a hypersonic cone-slice-flap variable geometry configuration given uncertainty in both the simulation and experimental data. The research here compares two different turbulence models for the simulations. For a variable geometry, performing uncertainty quantification to capture the model-form uncertainty on every configuration is computationally challenging. This work lays out a procedure that can give an accurate representation of the model-form uncertainty using a small number of high-fidelity runs and many low-fidelity runs on multiple configurations. Running this comparison provides a quantifiable measurement for the accuracy of each turbulence model for this type of design. The high-fidelity CFD solver used was VULCAN-CFD and the low-fidelity results came from Cart3D. The experimental data came from the 20-Inch Mach 6 Tunnel located at NASA Langley Research Center. The present work showed that the using both the Spalart and Allamaras and Menter Shear-Stress Transport turbulence models overpredicted the drag and lift coefficient, while underpredicting the pitching moment coefficient. The model-form uncertainty estimate resulted in up to a 13.6% change in the total uncertainty for the drag coefficient, up to a 57.4% change in total uncertainty for the lift coefficient, and up to a 100% change in total uncertainty for the pitching moment coefficient.

Laura M. White↗
NASA NTRSDate not supplied

CFD Validation Study of a Hypersonic Cone-Slice-Flap Variable Geometry Configuration

Model validation is the process of determining the degree of accuracy between physical reality and the model. The result of model validation can either be used to improve the model through calibration or quantify the model-form uncertainty. This work focuses on providing the model-form uncertainty through an area metric for a hypersonic cone-slice-flap variable geometry configuration given uncertainty in both the simulation and experimental data. For a variable geometry, performing uncertainty quantification to capture the model-form uncertainty on every configuration is computationally challenging. This work lays out a procedure that can give an accurate representation of the model-form uncertainty using a small number of high-fidelity runs and many low-fidelity runs on multiple configurations. Running this comparison provides a quantifiable measurement for the accuracy of each turbulence model for this type of design. The high-fidelity CFD solver used was VULCAN-CFD and the low-fidelity results came from Cart3D. The experimental data came from the 20-Inch Mach 6 Tunnel located at NASA Langley Research Center. The present work showed that the turbulence simulation overpredicted the drag and lift coefficient, while underpredicting the pitching moment coefficient. The model-form uncertainty estimate resulted up to a 13.6% change in the total uncertainty for the drag coefficient, up to a 57.4% change in total uncertainty for the lift coefficient, and up to a 100% change in total uncertainty for the pitching moment coefficient.

Laura White↗
NASA NTRS2005

Stirling Analysis Comparison of Commercial Versus High-Order Methods

Recently, three-dimensional Stirling engine simulations have been accomplished utilizing commercial Computational Fluid Dynamics software. The validations reported can be somewhat inconclusive due to the lack of precise time accurate experimental results from engines, export control/proprietary concerns, and the lack of variation in the methods utilized. The last issue may be addressed by solving the same flow problem with alternate methods. In this work, a comprehensive examination of the methods utilized in the commercial codes is compared with more recently developed high-order methods. Specifically, Lele's compact scheme and Dyson's Ultra Hi-Fi method will be compared with the SIMPLE and PISO methods currently employed in CFD-ACE, FLUENT, CFX, and STAR-CD (all commercial codes which can in theory solve a three-dimensional Stirling model with sliding interfaces and their moving grids limit the effective time accuracy). We will initially look at one-dimensional flows since the current standard practice is to design and optimize Stirling engines with empirically corrected friction and heat transfer coefficients in an overall one-dimensional model. This comparison provides an idea of the range in which commercial CFD software for modeling Stirling engines may be expected to provide accurate results. In addition, this work provides a framework for improving current one-dimensional analysis codes.

Dyson, Rodger W.↗
NASA NTRS2007

Stirling Analysis Comparison of Commercial vs. High-Order Methods

Recently, three-dimensional Stirling engine simulations have been accomplished utilizing commercial Computational Fluid Dynamics software. The validations reported can be somewhat inconclusive due to the lack of precise time accurate experimental results from engines, export control/ proprietary concerns, and the lack of variation in the methods utilized. The last issue may be addressed by solving the same flow problem with alternate methods. In this work, a comprehensive examination of the methods utilized in the commercial codes is compared with more recently developed high-order methods. Specifically, Lele's Compact scheme and Dyson s Ultra Hi-Fi method will be compared with the SIMPLE and PISO methods currently employed in CFD-ACE, FLUENT, CFX, and STAR-CD (all commercial codes which can in theory solve a three-dimensional Stirling model although sliding interfaces and their moving grids limit the effective time accuracy). We will initially look at one-dimensional flows since the current standard practice is to design and optimize Stirling engines with empirically corrected friction and heat transfer coefficients in an overall one-dimensional model. This comparison provides an idea of the range in which commercial CFD software for modeling Stirling engines may be expected to provide accurate results. In addition, this work provides a framework for improving current one-dimensional analysis codes.

Dyson, Rodger W.↗
NASA NTRS2012

Numerical Simulation of Noise from Supersonic Jets Passing Through a Rigid Duct

The generation, propagation and radiation of sound from a perfectly expanded Mach 2.5 cold supersonic jet flowing through an enclosed rigid-walled duct with an upstream J-deflector have been numerically simulated with the aid of OVERFLOW Navier-Stokes CFD code. A one-equation turbulence model is considered. While the near-field sound sources are computed by the CFD code, the far-field sound is evaluated by Kirchhoff surface integral formulation. Predictions of the farfield directivity of the OASPL (Overall Sound Pressure Level) agree satisfactorily with the experimental data previously reported by the author. Calculations also suggest that there is significant entrainment of air into the duct, with the mass flow rate of entrained air being about three times the jet exit mass flow rate.

Kandula, Max↗
NASA NTRSDate not supplied

CFD Validation Study of a Hypersonic Cone-Slice-Flap Configuration

Model validation is the process of determining the degree of accuracy between the real world and the model. The result of model validation can be used to either improve the model through calibration or quantify the model-form uncertainty. This work focuses on estimating the model-form uncertainty through an area metric for a hypersonic cone-slice-flap geometry configuration given the uncertainty in both the simulation and experimental data. A procedure that can give an accurate representation of the model-form uncertainty using a small number of high-fidelity runs is outlined. The work also assesses the impact of using different high-fidelity solvers and different turbulence models. The goal of performing this comparison is to provide a quantifiable measurement of the accuracy of each solver and turbulence model for this type of design. A low-fidelity analysis is also performed to get a model-form uncertainty for this type of analysis. The two high-fidelity CFD solvers used here are VULCAN-CFD and FUN3D, and the low-fidelity results comes from Cart3D.The experimental data comes from the 20-inch Mach 6 Tunnel located at NASA Langley Research Center. The work here shows that the models used tend to under-predict the drag and moment aerodynamic coefficients for this type of design leading to a model-form uncertainty estimate which can account for up to an 86% of the total uncertainty in the model predictions.

Laura M White↗