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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
NASA NTRS2015

Experimental, Numerical and Analytical Characterization of Slosh Dynamics Applied to In-Space Propellant Storage, Management and Transfer

Experimental and numerical results are presented from a new cryogenic fluid slosh program at the Florida Institute of Technology (FIT). Water and cryogenic liquid nitrogen are used in various ground-based tests with an approximately 30 cm diameter spherical tank to characterize damping, slosh mode frequencies, and slosh forces. The experimental results are compared to a computational fluid dynamics (CFD) model for validation. An analytical model is constructed from prior work for comparison. Good agreement is seen between experimental, numerical, and analytical results.

slosh mode frequencies↗
NASA NTRS2016

Impact of an Exhaust Throat on Semi-Idealized Rotating Detonation Engine Performance

A computational fluid dynamic (CFD) model of a rotating detonation engine (RDE) is used to examine the impact of an exhaust throat (i.e., a constriction) on performance. The model simulates an RDE which is premixed, adiabatic, inviscid, and which contains an inlet valve that prevents backflow from the high pressure region directly behind the rotating detonation. Performance is assessed in terms of ideal net specific impulse which is computed on the assumption of lossless expansion of the working fluid to the ambient pressure through a notional diverging nozzle section downstream of the throat. Such a semi-idealized analysis, while not real-world, allows the effect of the throat to be examined in isolation from, rather than coupled to (as it actually is) various loss mechanisms. For the single Mach 1.4 flight condition considered, it is found that the addition of a throat can yield a 9.4 percent increase in specific impulse. However, it is also found that when the exit throat restriction gets too small, an unstable type of operation ensues which eventually leads to the detonation failing. This behavior is found to be somewhat mitigated by the addition of an RDE inlet restriction across which there is an aerodynamic loss. Remarkably, this loss is overcome by the benefits of the further exhaust restrictions allowed. The end result is a configuration with a 10.3 percent improvement in ideal net specific thrust.

combustion↗
NASA NTRS2016

Impact of an Exhaust Throat on Semi-Idealized Rotating Detonation Engine Performance

A computational fluid dynamic (CFD) model of a rotating detonation engine (RDE) is used to examine the impact of an exhaust throat (i.e. a constriction) on performance. The model simulates an RDE which is premixed, adiabatic, inviscid, and which contains an inlet valve that prevents backflow from the high pressure region directly behind the rotating detonation. Performance is assessed in terms of ideal net specific impulse which is computed on the assumption of lossless expansion of the working fluid to the ambient pressure through a notional diverging nozzle section downstream of the throat. Such a semi-idealized analysis, while not real-world, allows the effect of the throat to be examined in isolation from, rather than coupled to (as it actually is) various loss mechanisms. For the single Mach 1.4 flight condition considered, it is found that the addition of a throat can yield a 9.4 percent increase in specific impulse. However, it is also found that when the exit throat restriction gets too small, an unstable type of operation ensues which eventually leads to the detonation failing. This behavior is found to be somewhat mitigated by the addition of an RDE inlet restriction across which there is an aerodynamic loss. Remarkably, this loss is overcome by the benefits of the further exhaust restrictions allowed. The end result is a configuration with a 10.3 percent improvement in ideal net specific thrust.

propulsion↗
NASA NTRS2016

Zero Boil-Off Tank (ZBOT) Experiment

The Zero-Boil-Off Tank (ZBOT) experiment has been developed as a small scale ISS experiment aimed at delineating important fluid flow, heat and mass transport, and phase change phenomena that affect cryogenic storage tank pressurization and pressure control in microgravity. The experiments use a simulant transparent low boiling point fluid (PnP) in a sealed transparent Dewar to study and quantify: (a) fluid flow and thermal stratification during pressurization; (b) mixing, thermal destratification, depressurization, and jet-ullage penetration during pressure control by jet mixing. The experiment will provide valuable microgravity empirical two-phase data associated with the above-mentioned physical phenomena through highly accurate local wall and fluid temperature and pressure measurements, full-field phase-distribution and flow visualization. Moreover, the experiments are performed under tightly controlled and definable heat transfer boundary conditions to provide reliable high-fidelity data and precise input as required for validation verification of state-of-the-art two-phase CFD models developed as part of this research and by other groups in the international scientific and cryogenic fluid management communities.

ISS↗
NASA NTRS2018

Characterization and Analysis of Phoca Vitulina, Zalophus Californianus and Mirounga Angustirostris Vibrissae

Vibrissae of Phoca Vitulina (harbor seal) and Mirounga Angustirostris (elephant seal) possess undulations along their length. Harbor seal vibrissae have shown potential to reduce vortex induced vibrations and reduce drag compared to cylinders and ellipses. The exact geometry of the whiskers has not been well documented and the parameters that are responsible for the reduction in drag and vortex induced vibrations have not been characterized. Samples of six harbor seal vibrissae, six elephant seal vibrissae and six California sea lion (Zalophus californianus) vibrissae were collected from the Marine Mammal Center in California. The objectives of this study were to (1) Compare measurement techniques for digitizing and extracting parameters of the seal whiskers for the PeTaL (Periodic Table of Life) database. CT scanning, microscopy and 3D scanning techniques were compared. (2) Compare aerodynamic characteristics of a representative harbor seal whisker, elephant seal whisker, California sea lion whisker and ellipse at Re = 12000 and Re = 23000 based on major axis and free stream velocity. The data (in appendices) is available to compare CFD models or for further experimental validation, (3) Show close up images of whiskers and look for surface roughness effects. Variations in the seven parameters of the seal whisker were observed that may either be a feature of the vibration reduction mechanism or a result of natural variation. It is hypothesized that six parameters are sufficient to characterize seal whiskers based on analytic fitting. The drag coefficient of harbor seal whiskers examined in this study were found to be 25 percent lower than that of an ellipse with comparable major and minor axis lengths at Reynolds number of 12000. The dissipation length scale was found to be larger for seal whiskers. Potential applications of seal whisker morphology for aerospace are discussed. Roughness is not thought to play a factor in seal hydrodynamics.

Shyam, Vikram↗
NASA NTRSDate not supplied

Computational Fluid Dynamics Simulations of the Transonic Dynamics Tunnel Airstream Oscillator System

This paper presents a Computational Fluid Dynamics (CFD) model of the flow in the NASA Langley Research Center Transonic Dynamics Tunnel (TDT) with the Airstream Oscillator System (AOS) in operation. The TDT is a continuous-flow, closed circuit wind tunnel with a 16- by 16-foot slotted test section with cropped corners. The tunnel was originally built as the 19-ft Pressure Tunnel in 1938, but it was converted to the current transonic tunnel in the 1950s, with capabilities to use either air or heavy gas (R-134a) as the test medium. The TDT was also fitted with an AOS which can be used to create a gust field in the test section. To date, no computational analyses of the TDT involving the AOS have been performed. In this study, experimental data acquired of an oscillating airstream in the tunnel will be used to calibrate the computational analyses. A significant motivation of this work is to attempt to compare computational gust velocities to those recorded in the TDT literature. In addition to a validation of this model with experimental data, it may also be possible to supplement the experimental data with computational data. The experimental data is rather sparse and at only a few Mach numbers. Computational results may be able to expand the AOS data set. Another motivation for this work is that the Integrated Adaptive Wing Technology Maturation (IAWTM) semi-span model of the high-aspect-ratio CRM with 10 trailing edge control surfaces is currently being fabricated and will be delivered to the NASA Langley Transonic Dynamics Tunnel (TDT) for testing in late 2020. Among tests to be conducted, gust load alleviation (GLA) will be demonstrated at transonic conditions using the AOS.

Computational Fluid Dynamics↗
NASA NTRSDate not supplied

Development of Physics-Based Transition Models for Unstructured-Mesh CFD Codes Using Deep Learning Models

Predicting transition locations over a vehicle surface is of fundamental importance for many engineering applications. With the transition information, the Reynolds-averaged Navier-Stokes (RANS) computations can turn on the turbulence model at the right locations so that drag, lift and other aerodynamic quantities can be accurately predicted. In contrast to the popularity of RANS-based transition modeling in which transition onset is governed by the turbulence equations, physics-based transition models that account for instability waves within the boundary layer, thus more compliant to flow physics, only gained more attention in recent years. This paper describes the development of a new physics-based transition model based on either the linear stability theory (LST) or parabolized stability equations (PSE). The model is designed to communicate with a structured or unstructured-mesh RANS solver back and forth in order to more accurately compute transition fronts over a three-dimensional body. In the developed model, the Python suite of interface codes in conjunction with the LASTRAC software can be executed autonomously to produce transition onset locations for a given laminar or RANS-computed transitional state. In addition, as a proof of concept, the tool set consists of a deep learning neural network model that has been designed and trained to predict instability wave evolutions inside the boundary layer for various instability wave mechanisms across a selected speed range. A machine-learned intelligent profile interpolation model has also been devised to enable reliable instability-wave spectra predictions with just a few points in the mean flow profiles.

Transition Modeling↗
NASA NTRSDate not supplied

Interfacial thermodynamics of cryogenic fluids: The effect on non-condensable gas on fluid storage

Propellant tanks that contain cryogenic fluids (CFs) in low gravity conditions are usually pressurized with non-condensable (NC) gases for fast extraction. Unfortunately, the presence of NC gases causes CFs to exhibit higher boil-offs compared to pure CF systems. For optimal utilization of the cryogenic fuels, these higher boil-off rates must be minimized. Our goal is to quantify the effects NC gases have on the evaporation and condensation dynamics of CFs to develop strategies that mitigate the higher boil-offs. The hypothesis is that the NC gas accumulates around the CF’s liquid vapor interface in the Knudsen layer, creating a kinetic barrier for both evaporation and condensation. Testing this hypothesis with experimental methods is difficult due to the transient nature of the Knudsen layer, which is only a few nanometers thick. Furthermore, experimental capabilities are limited in low gravity conditions and extremely expensive. Accordingly, we approach this problem from a theoretical perspective. In this study, we employ molecular dynamics (MD) simulations to probe the interfacial mechanisms that affect CF’s evaporation and condensation at varying concentrations of NC gases. Specifically, we use nitrogen (N2) and oxygen (O2) as our CFs and neon (Ne) as our NC gas. Using MD simulations, we show that Ne accumulates at N2 liquid vapor interface across a wide range of Ne concentrations, thereby impeding mass transport of N2. Our simulations allow for direct computation of the molar flux as well as the mass accommodation coefficient (MAC) which can then be used as input parameters to continuum fluid dynamics (CFD) models for optimal storage tank design.

Michael Robert DeLyser↗
NASA NTRS2022

Flow Visualization to Complement Modern Wing Tunnel Testing

In recent times, wind tunnel testing at the NASA Ames Unitary Plan Wind Tunnel (UPWT) seeks to validate computational fluid dynamics (CFD) tools. To that end, we find that now more than ever, flow visualization is a powerful capability that provides context to the experimental results bringing real-time insight to the researcher. It enables us to answer questions like: “Why does this drag polar look different?” or “What is the source of this tone in the microphone data?”. Flow visualization can verify that assumptions made in the CFD model are correct, addressing questions like: “Are our trip dots the correct size and location?” or “Does the position of this shock match our prediction?”. Techniques such as infrared flow visualization and schlieren / shadowgraph are powerful capabilities that when taken together can provide results as to the nature of the fluid structure interaction at a test condition in a format that pairs well with computational results. When properly integrated these measurements can be taken with no impact on the productivity of the test program, supplementing traditional measurements and often discovering something unforeseen, helping the researcher to answer those slippery aerodynamic questions and ultimately build confidence in the accuracy of a model.

wind tunnel↗
NASA NTRSDate not supplied

Measurements of Temperature and Enthalpy in NASA Ames 60-MW Arcjet using Atomic Oxygen and Atomic Nitrogen Absorption

We report on measurements of temperature and enthalpy in the 60-MW Interaction Heating Facility (IHF) Arcjet at the NASA Ames Research Center. These quantities were measured at sub-second temporal resolution using tunable diode laser absorption spectroscopy, employing four Distributed Feedback (DFB) lasers to simultaneously probe the arcjet flow. Optical access was achieved through a specially-designed optical disk mounted in series with the modules of the arc heater in the add-air plenum of the column, which is positioned between the downstream electrodes and the nozzle. Gas at this measurement location serves as the reservoir condition for the ensuing supersonic expansion. Conditions in the reservoir are extreme, with temperatures between 5000-8000 K and pressures between 1-9 atm. Absorption from the 5S state of atomic oxygen at 777.2 nm and the 4P state of atomic nitrogen at 868.0 nm, two electronically-excited quantum states, formed the basis of the flow characterization measurements. Fiber-coupled laser light was pitched through the IHF via 1/8" ports in the optical disks, permitting optical access with minimal intrusion to the flow. Data was collected at 3 nominal conditions, a Low condition (1600 A, 2000 V, 0.13 kg/s), a Medium condition (3500 A, 5700 V, 0.61 kg/s), and a High condition (5900 A, 6800 V, 0.84 kg/s). Measurements suggest agreement within 1-2% between excited oxygen and excited nitrogen thermometry, indicating approximate thermal and chemical equilibration of the arcjet reservoir state, a result in keeping with long-standing assumptions of CFD modelling. Spatial gradients are found in the flow due to enthalpy loss to the facility walls, but the measured spatial non-uniformities are consistent with axi-symmetry. Measurements of specific enthalpy in the reservoir are also made utilizing excited oxygen and excited nitrogen absorption under the assumption of chemical equilibration.

Peter M. Finch↗
NASA NTRSDate not supplied

Energy of Fluid (EOF) Simulation of Large-Scale Densification and Solidification of Hydrogen

The current research develops a finite volume-based Computation Fluid Dynamic (CFD) model utilizing an Energy of Fluid (EOF) approach to generate a simulation of the densification of liquid hydrogen (LH2) in the Integrated Refrigeration and Storage (IRAS) experimental tank for the Ground Operations Demonstration Unit for Liquid Hydrogen Project (GODU-LH2) at the Kennedy Space Center (KSC). The computational code will integrate a commercial software pressure-based model with User Defined Functions (UDF) generated and added to the model. The modified model will solve the energy equation in terms of internal energy and provide temperature and pressure calculations for a given tank geometry. This method will decrease computational time when compared to the built-in commercially available Enthalpy solvers. Currently the research focused on utilizing User Defined Scalars (UDS) compared to the enthalpy model and a lumped node analysis to show the written UDS is as accurate as the built-in solvers, and the models will be verified using known solutions for solidification and validated with the experimental results from the Cryogenics Test Laboratory (CTL) at KSC for the IRAS densification process.

C. Mahony↗
NASA NTRS2023

Simulations of Zero Boiloff, Densification, and Solidification in a Large-Scale Liquid Hydrogen Tank

The current research utilizes an Energy of Fluid (EOF) approach to develop a finite volume-based Computation Fluid Dynamic (CFD) model to create a benchmark simulation of the densification of liquid hydrogen (LH2) in an experimental Integrated Refrigeration and Storage (IRAS) tank at the Kennedy Space Center (KSC). The computational code will incorporate the FLUENT pressure-based model with User Defined Functions (UDF) to implement the EOF model. The enhanced model will solve the energy equation in terms of internal energy and provide temperature and pressure calculations for a given tank geometry. Specifically, a convection term will be added to an already existing code for a more accurate simulation of LH2 densification. However, currently the research focuses on utilizing User Defined Scalars (UDS) with FLUENT to show the built-in solvers are as accurate as the written UDS. This will reduce the computational time down significantly by utilizing an already built in solver for the momentum equation. A custom UDF has been written to solve transient and steady state conduction, and results obtained using the UDF agree with results obtained using the built-in FLUENT solver. Currently, a convection UDF is being added and will be verified using known solutions for solidification and validated with the experimental results from the Cryogenics Test Laboratory (CTL) at KSC for the IRAS densification process.

Colin Philip Mahony↗
NASA NTRSDate not supplied

Cryogenic Flow Boiling Parabolic Flight Experiment and Universal Correlations for Cryogens

Cryogenic fluids are employed in many industries and play crucial roles in space applications, including nuclear thermal propulsion systems, LOX/LCH4 ascent and descent stage feedlines, LOX/LH2 transfer lines for in-space fuel depots, and LHe transfer lines that cool space experiments. Due to the ultra-low boiling point of cryogens, phase change is inevitable in these systems, which leads to complex boiling behavior in reduced gravity. For cost-effective design and sizing of the aforementioned systems, accurate knowledge of two-phase parameters such as pressure drop, heat transfer coefficient (HTC), and critical heat flux (CHF) is desired. Since 2018, researchers from Purdue University and NASA Glenn Research Center (GRC) have been collaborating to investigate the effects of reduced gravity on cryogenic pool boiling, flow boiling ins tubes, and spray cooling. Using a variety of both terrestrial and parabolic flight experimental rigs, extensive LN2 data have been amassed and combined with cryogenic data from World literature to develop both ‘universal’ empirical correlations and CFD models, which are expected to serve as foundation for design of future space systems involving cryogenic fluids. This presentation will provide an overview of both the experimental work and predictive tools that have resulted from the joined Purdue-Glenn studies.

Flow Boiling↗
NASA NTRS2008

Development of a Common Research Model for Applied CFD Validation Studies

The development of a wing/body/nacelle/pylon/horizontal-tail configuration for a common research model is presented, with focus on the aerodynamic design of the wing. Here, a contemporary transonic supercritical wing design is developed with aerodynamic characteristics that are well behaved and of high performance for configurations with and without the nacelle/pylon group. The horizontal tail is robustly designed for dive Mach number conditions and is suitably sized for typical stability and control requirements. The fuselage is representative of a wide/body commercial transport aircraft; it includes a wing-body fairing, as well as a scrubbing seal for the horizontal tail. The nacelle is a single-cowl, high by-pass-ratio, flow-through design with an exit area sized to achieve a natural unforced mass-flow-ratio typical of commercial aircraft engines at cruise. The simplicity of this un-bifurcated nacelle geometry will facilitate grid generation efforts of subsequent CFD validation exercises. Detailed aerodynamic performance data has been generated for this model; however, this information is presented in such a manner as to not bias CFD predictions planned for the fourth AIAA CFD Drag Prediction Workshop, which incorporates this common research model into its blind test cases. The CFD results presented include wing pressure distributions with and without the nacelle/pylon, ML/D trend lines, and drag-divergence curves; the design point for the wing/body configuration is within 1% of its max-ML/D. Plans to test the common research model in the National Transonic Facility and the Ames 11-ft wind tunnels are also discussed.

Vassberg, John C.↗
NASA NTRSDate not supplied

Reduced-Order Aerodynamic Modeling Based on CFD Frequency Responses from Multisine Inputs

A system identification analysis was performed to determine a reduced-order model (ROM) of a computational fluid dynamics (CFD) solver in support of linear aeroservoelastic model development and feedback control design. The approach was applied to the FUN3D code for the half-span wind tunnel test article used in the NASA-Boeing collaboration called the Integrated Adaptive Wing Technology Maturation (IAWTM) project. In a transonic flow condition, multiple inputs (11 structural mode displacements and 3 control surface deflections) were simultaneously excited with orthogonal phase-optimized multisines while multiple outputs (the corresponding 14 generalized aerodynamic forces) were recorded. From these recorded times series, the matrix of frequency responses was computed and subsequently fit using rational function approximations (RFAs). It was found that the entire (14 x 14) matrix of frequency responses could be determined from a single CFD run and that results generally followed trends predicted using other methods. Differences were attributed to the modeling fidelity and nonlinearities from structural mode and control surface interactions at higher reduced frequencies. More accurate fits of the RFAs to the frequency response data were obtained by making two CFD runs, one with only structural mode excitations and one with only control surface excitations, which reduced the degree of nonlinearity in the modeling data.

Aeroservoelasticity↗
NASA NTRS2021

Pretest Computational Assessment of Boundary Layer Transition in the NASA Juncture Flow Model with an NACA 0015-Based Wing

The first two phases of the NASA Juncture Flow experiment were carried out on a DLR-F6 swept-wing model and were designed to provide “CFD validation-quality” data toward the assessment and improvement of existing CFD turbulence models in predicting onset and extent of three-dimensional separated flow near the wing-juncture trailing-edge region. The next phase of experiments will involve an NACA 0015-based swept wing, as prior risk reduction experiments had indicated that this wing shape resulted in reduced separation near the juncture region than the DLR-F6 wing, thus providing a better option to evaluate the ability of CFD models to predict incipient turbulent separation. The NACA 0015 measurements will also include IR thermography to infer the variation of transition front with respect to an increasing angle of attack. The primary objective of this work is to computationally make a preliminary assessment of the transition front on both surfaces of the NACA 0015 wing at a crank-chord-based Reynolds number of 2.4 x 106 for four different angles of attack, (0°, 2.5°, 5°, and 7.5°) and to determine the dominant mechanisms responsible for transition. This assessment includes both RANS-based transition models from NASA’s OVERFLOW 2.3b flow solver and linear parabolized stability equations (PSE) stability analysis based on the Langley Stability and Transition Analysis code, LASTRAC. Linear PSE results indicate that the upper surface of the wing is dominated by Tollmien- Schlichting (TS) instabilities, and that the laminar flow region shrinks from about 50% chord to a very small region just downstream of the attachment line as the angle of attack is increased from 0° to 7.5°. Consequently, the transition fronts predicted by the Spalart- Allmaras-based amplification factor transport (AFT-2017b) equation model (which accounts for the TS instabilities alone) and the Menter’s shear-stress transport equation (SST2003)- based Langtry-Menter transition model with ability to account for both TS and crossflow effects (LM2015) compare well with those predicted using linear PSE. On the lower surface of the wing, stationary crossflow (CF) instabilities begin to appear on the inboard portion of the wing in addition to the TS-instabilities for the larger angles of attack (5° and 7.5°), further reducing the laminar flow extent within the inboard region. The LM2015 model that accounts for CF effects is able to replicate this trend but appears to predict a slightly earlier transition. The outcome of this effort will inform the experiment and, when the actual experimental data become available, provide further opportunity to assess and improve the various transition models.

CFD modeling↗
NASA NTRS2021

Pretest Computational Assessment of Boundary Layer Transition in the NASA Juncture Flow Model with an NACA 0015-Based Wing

The first two phases of the NASA Juncture Flow experiment were carried out on a DLR-F6 swept-wing model and were designed to provide “CFD validation-quality” data toward the assessment and improvement of existing CFD turbulence models in predicting onset and extent of three-dimensional separated flow near the wing-juncture trailing-edge region. The next phase of experiments will involve an NACA 0015-based swept wing, as prior risk reduction experiments had indicated that this wing shape resulted in reduced separation near the juncture region than the DLR-F6 wing, thus providing a better option to evaluate the ability of CFD models to predict incipient turbulent separation. The NACA 0015 measurements will also include IR thermography to infer the variation of transition front with respect to an increasing angle of attack. The primary objective of this work is to computationally make a preliminary assessment of the transition front on both surfaces of the NACA 0015 wing at a crank-chord-based Reynolds number of 2.4 x 106 for four different angles of attack, (0°, 2.5°, 5°, and 7.5°) and to determine the dominant mechanisms responsible for transition. This assessment includes both RANS-based transition models from NASA’s OVERFLOW 2.3b flow solver and linear parabolized stability equations (PSE) stability analysis based on the Langley Stability and Transition Analysis code, LASTRAC. Linear PSE results indicate that the upper surface of the wing is dominated by Tollmien- Schlichting (TS) instabilities, and that the laminar flow region shrinks from about 50% chord to a very small region just downstream of the attachment line as the angle of attack is increased from 0° to 7.5°. Consequently, the transition fronts predicted by the Spalart- Allmaras-based amplification factor transport (AFT-2017b) equation model (which accounts for the TS instabilities alone) and the Menter’s shear-stress transport equation (SST2003)- based Langtry-Menter transition model with ability to account for both TS and crossflow effects (LM2015) compare well with those predicted using linear PSE. On the lower surface of the wing, stationary crossflow (CF) instabilities begin to appear on the inboard portion of the wing in addition to the TS-instabilities for the larger angles of attack (5° and 7.5°), further reducing the laminar flow extent within the inboard region. The LM2015 model that accounts for CF effects is able to replicate this trend but appears to predict a slightly earlier transition. The outcome of this effort will inform the experiment and, when the actual experimental data become available, provide further opportunity to assess and improve the various transition models.

CFD modeling↗