Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computational fluid dynamics simulation”

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,675 records · Page 93
All sourcesNASA NTRSDOE OSTI
NASA NTRS2024

Improved Chemistry and Attenuation Models for Communication Blackout Simulation During Mars 2020 Entry

As a blunt body enters a planetary atmosphere, a plasma forms in the hypersonic shock layer and attenuates radio communication causing signal blackout for some duration of the entry sequence. In our previous work,1 computational fluid dynamics (CFD) was applied to model the entry flow around the Mars 2020 spacecraft, including ionization and electron density throughout the flow field, and predict ultra-high frequency (UHF) radio wave attenuation due to electrons. In total, 17 chemical species and their spatial profiles are modelled around the Mars 2020 spacecraft at 11 different points in time during entry. Although the simulation predicted the onset of attenuation well, the timing of the end of the predicted blackout window significantly preceded the end time observed during the 2021 landing. The present work seeks to improve the attenuation model by accounting for the fact that electrons undergo collisions with heavier species in the flow, which is an effect that was neglected in previous analyses. It is determined that including electron collisions increases the overall magnitude of attenuation predicted especially towards the end of the measured attenuation period, improving qualitative agreement between predicted and measured attenuation to both spacecraft receiving the signal from Mars 2020. To explore the remaining uncertainty in signal attenuation predictions further, a sensitivity study is performed to investigate the impact of associative ionization and electron-impact ionization rate coefficients on the electron density predicted by CFD and on the resulting attenuation predictions. These coefficients are believed to contain up to order-of-magnitude uncertainty, and therefore may significantly affect the number density of electrons throughout the flow field. Variations in associative ionization coefficients demonstrate significant impact on the magnitude of attenuation due to variation in the electron density coming from associative ionization. However, the start and end times of the predicted signal attenuation period are only slightly impacted by said variation.

Eve Papajak↗
DOE OSTI2022

Impact of Variable Gas Mixtures on Bubble Size Distribution and Mass Transfer in Gas Fermentation Reactors

Gas fermentation has emerged as a promising new technology for the generation of fuels and chemicals from mixtures of greenhouse and energy rich gas streams (CO2/CH4/H2/CO) via microbial bioreaction. Example pathways include biomethanation (CO2/H2 to CH4), biogas upgrading, CO fermentation and wet-waste conversion. Effective Gas-liquid mass-transfer is an important physical phenomenon that determines the design and scale-up of these systems. There is currently a knowledge-gap regarding bubble-size distributions when using a mixture of gases with vastly different properties, which can have a significant impact on overall mass-transfer. For example, hydrogen bubbles are more buoyant compared to other relatively heavier gases (CO2/CH4/CO), resulting in a large distribution of residence times and bubble sizes. This work therefore develops a deeper understanding of bubble dynamics and interphase mass transfer in such heterogenous gas mixtures through well-resolved computational models. We use a detailed multiphase computational-fluid-dynamics (CFD) model to study the impact of gas-mixtures on overall mass-transfer in bubble column and air-lift reactors. The CFD tool previously developed by the authors (1) for simulating aerobic fermentation reactors at scale is used in this study. The Reynolds-averaged mass, momentum, energy, and species transport equations are solved for interpenetrating gas and liquid phase in this model. We use a population balance-based bubble-size-distribution model that is validated against small-scale experiments in our solver. Results pertaining to multiple simulations of gas-fermentation reactors are presented where gas mixtures with varying compositions of CO2/CH4/CO/H2 are imposed at the sparger boundaries. The spatio-temporal variations in bubble-size distribution and mass transfer coefficient is analyzed for varying superficial velocities and gas-compositions for varying sizes of bubble-column and airlift reactors.

BIOMASS FUELS,MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Uncertainty Quantification for Multiphase Computational Fluid Dynamics Closure Relations with a Physics-Informed Bayesian Approach

Multiphase Computational Fluid Dynamics (MCFD) based on the two-fluid model is considered a promising tool to model complex two-phase flow systems. MCFD simulation can predict local flow features without resolving interfacial information. As a result, the MCFD solver relies on closure relations to describe the interaction between the two phases. Those empirical or semi-mechanistic closure relations constitute a major source of uncertainty for MCFD predictions. In this paper, we leverage a physics-informed uncertainty quantification (UQ) approach to inversely quantify the closure relations’ model form uncertainty in a physically consistent manner. This proposed approach considers the model form uncertainty terms as stochastic fields that are additive to the closure relation outputs. Combining dimensionality reduction and Gaussian processes, the posterior distribution of the stochastic fields can be effectively quantified within the Bayesian framework with the support of experimental measurements. As this UQ approach is fully integrated into the MCFD solving process, the physical constraints of the system can be naturally preserved in the UQ results. Here, in a case study of adiabatic bubbly flow, we demonstrate that this UQ approach can quantify the model form uncertainty of the MCFD interfacial force closure relations, thus effectively improving the simulation results with relatively sparse data support.

42 ENGINEERING↗
NASA NTRS1975

Dynamic heave-pitch analysis of air cushion landing systems

A program to develop analytical tools for evaluating the dynamic performance of Air Cushion Landing Systems (ACLS) is described. The heave (vertical) motion of the ACLS was analyzed, and the analysis was extended to cover coupled heave-pitch motions. The mathematical models developed are based on a fundamental analysis of the body dynamics and fluid mechanics of the aircraft-cushion-runway interaction. The air source characteristics, flow losses in the feeding ducts, trunk and cushion, the effects of fluid compressibility, and dynamic trunk deflections, including ground contact are considered. A computer program, based on the heave-pitch analysis, was developed to simulate the dynamic behavior of an ACLS during landing impact and taxi over an irregular runway. The program outputs include ACLS motions, loadings, pressures, and flows as a function of time. To illustrate program use, three basic types of simulations were carried out. The results provide an initial indication of ACLS performance during (1) a static drop, (2) landing impact, and (3) taxi over a runway irregularity.

Captain, K. M.↗
NASA NTRS1984

Coolant passage heat transfer with rotation, a progress report

The objective of this 36-month experimental and analytical program is to develop a heat transfer and pressure drop database, computational fluid dynamic techniques, and correlations for multipass rotating coolant passages with and without flow turbulators. The experimental effort will be focused on the simulation of configurations and conditions expected in the blades of advanced aircraft high pressure turbines so that the effects of Coriolis and buoyancy forces on the coolant side flow can be rationally included in the design of turbine blades.

Kopper, F. C.↗
NASA NTRS2003

Performance Trades Study for Robust Airfoil Shape Optimization

From time to time, existing aircraft need to be redesigned for new missions with modified operating conditions such as required lift or cruise speed. This research is motivated by the needs of conceptual and preliminary design teams for smooth airfoil shapes that are similar to the baseline design but have improved drag performance over a range of flight conditions. The proposed modified profile optimization method (MPOM) modifies a large number of design variables to search for nonintuitive performance improvements, while avoiding off-design performance degradation. Given a good initial design, the MPOM generates fairly smooth airfoils that are better than the baseline without making drastic shape changes. Moreover, the MPOM allows users to gain valuable information by exploring performance trades over various design conditions. Four simulation cases of airfoil optimization in transonic viscous ow are included to demonstrate the usefulness of the MPOM as a performance trades study tool. Simulation results are obtained by solving fully turbulent Navier-Stokes equations and the corresponding discrete adjoint equations using an unstructured grid computational fluid dynamics code FUN2D.

Li, Wu↗
NASA NTRS2009

Modeling Unsteady Cavitation and Dynamic Loads in Turbopumps

A computational fluid dynamics (CFD) model that includes representations of effects of unsteady cavitation and associated dynamic loads has been developed to increase the accuracy of simulations of the performances of turbopumps. Although the model was originally intended to serve as a means of analyzing preliminary designs of turbopumps that supply cryogenic propellant liquids to rocket engines, the model could also be applied to turbopumping of other liquids: this can be considered to have been already demonstrated, in that the validation of the model was performed by comparing results of simulations performed by use of the model with results of sub-scale experiments in water. The need for this or a similar model arises as follows: Cavitation instabilities in a turbopump are generated as inlet pressure drops and vapor cavities grow on inducer blades, eventually becoming unsteady. The unsteady vapor cavities lead to rotation cavitation, in which the cavities detach from the blades and become part of a fluid mass that rotates relative to the inducer, thereby generating a fluctuating load. Other instabilities (e.g., surge instabilities) can couple with cavitation instabilities, thereby compounding the deleterious effects of unsteadiness on other components of the fluid-handling system of which the turbopump is a part and thereby, further, adversely affecting the mechanical integrity and safety of the system. Therefore, an ability to predict cavitation- instability-induced dynamic pressure loads on the blades, the shaft, and other pump parts would be valuable in helping to quantify safe margins of inducer operation and in contributing to understanding of design compromises. Prior CFD models do not afford this ability. Heretofore, the primary parameter used in quantifying cavitation performance of a turbopump inducer has been the critical suction specific speed at which head breakdown occurs. This parameter is a mean quantity calculated on the basis of assumed steady-state operation of the inducer; it does not account for dynamic pressure loads associated with unsteady flow caused by instabilities. Because cavitation instabilities occur well before mean breakdown in inducers, engineers have, until now, found it necessary to use conservative factors of safety when analyzing the results of numerical simulations of flows in turbopumps.

Hosangadi, Ashvin↗
NASA NTRS2015

Computational Fluid Dynamics Study on the Effects of RATO Timing on the Scale Model Acoustic Test

The Scale Model Acoustic Test (SMAT) is a 5% scale test of the Space Launch System (SLS), which is currently being designed at Marshall Space Flight Center (MSFC). The purpose of this test is to characterize and understand a variety of acoustic phenomena that occur during the early portions of lift off, one being the overpressure environment that develops shortly after booster ignition. The SLS lift off configuration consists of four RS-25 liquid thrusters on the core stage, with two solid boosters connected to each side. Past experience with scale model testing at MSFC (in ER42), has shown that there is a delay in the ignition of the Rocket Assisted Take Off (RATO) motor, which is used as the 5% scale analog of the solid boosters, after the signal to ignite is given. This delay can range from 0 to 16.5ms. While this small of a delay maybe insignificant in the case of the full scale SLS, it can significantly alter the data obtained during the SMAT due to the much smaller geometry. The speed of sound of the air and combustion gas constituents is not scaled, and therefore the SMAT pressure waves propagate at approximately the same speed as occurs during full scale. However, the SMAT geometry is much smaller allowing the pressure waves to move down the exhaust duct, through the trench, and impact the vehicle model much faster than occurs at full scale. To better understand the effect of the RATO timing simultaneity on the SMAT IOP test data, a computational fluid dynamics (CFD) analysis was performed using the Loci/CHEM CFD software program. Five different timing offsets, based on RATO ignition delay statistics, were simulated. A variety of results and comparisons will be given, assessing the overall effect of RATO timing simultaneity on the SMAT overpressure environment.

Nielsen, Tanner↗
NASA NTRS2019

Concept Evaluation: NASA Reference Model 2-Side by Side

Urban Air Mobility (UAM) is an up and coming topic in the world of aerosciences. A variety of companies, most notably Uber, has begun working towards making air taxis a regular part of day to day life. However, the implementation of this system into society will require a significant amount of work: making helicopters that can handle these flights and operate safely in urban environments, having landing zones, and working out new technologies to lower costs. NASA has also taken interest in this new idea and has formed a focus group of interns to tackle some of these problems. They have started by designing a few concept models for what these air taxis could look like, focusing on low noise, multiple passenger configurations, and the potential for electric or hybrid helicopters. One particular model is known as NR2, or the "Side by Side," which features two rotors spinning in opposite directions. Though these models are purely conceptual at this point, understanding the real-world aerodynamics and performance of this model is crucial to its future development. Using a program known as Rotorcraft Computational Fluid Dynamics (RotCFD), the performance characteristics of the rotors and aircraft, such as the lift, drag, thrust and power can be determined. These simulations can be run in a variety of flight configurations, with hover, forward flight, and climb being focused on in this study. The rotor collective angle can also be changed and used to get results about the lift and drag for the blades and the model body. These simulations can provide information for future model construction and wind tunnel testing.

Reference Model↗
NASA NTRSDate not supplied

Computational Analysis of a Quiet Single-Main Rotor Helicopter for Air Taxi Operations

A computational study has been done on NASA’s quiet single-main rotor helicopter concept for urban air mobility using high-fidelity computational fluid dynamics, rotorcraft comprehensive analysis tools, and computational aeroacoustics. High-order accurate schemes, dual-time stepping, and the delayed detached-eddy simulation model have been employed. A loose-coupling approach between the flow solver and the rotorcraft comprehensive code is implemented to include vehicle trim and blade motions. Acoustic simulations based on the Ffowcs-Williams and Hawkings equations have been performed to compute the rotor noise. Different blade geometries and tip speeds are analyzed, and performance and acoustics results are compared. The vehicle has been simulated in hover and cruise, with flight conditions representative of an air taxi mission. The quiet single-main rotor helicopter is one of the conceptual designs intended to focus and guide NASA’s research activities in support of aircraft development for vertical take-off and landing air taxi operations.

ARMD↗
NASA NTRSDate not supplied

Closed-Loop Simulations of Human-Scale Mars Lander Descent Trajectories on Frontier

A computational campaign was performed to run high-fidelity, free-flight simulations of a human-scale Mars lander concept vehicle decelerating under retropropulsion through the Martian atmosphere with closed-loop flight control. A novel approach is used to couple computational fluid dynamics (CFD) software with a mature flight mechanics package, where the two applications communicate in real-time across two geographically-dispersed computational facilities. The CFD is performed on the Frontier exascale system located at Oak Ridge National Laboratory, and the flight mechanics are executed on a system located at NASA Langley Research Center. In the current campaign, CFD is performed using finite-rate chemistry to account for the interactions between the LOXCH 4 engines and the CO 2 Martian atmosphere. A simulation of a closed-loop main engine throttling and RCS actuation is presented, demonstrating that the vehicle and model are able to maintain stability in a long-duration CFD-in-the-loop flight simulation. Comparisons are made to a reduced order model ignoring aero-propulsive interactions.

CFD↗
NASA NTRS1985

Large-scale computations in fluid mechanics; Proceedings of the Fifteenth Summer Seminar on Applied Mathematics, University of California, La Jolla, CA, June 27-July 8, 1983. Parts 1 & 2

Papers are presented on such topics as the use of semi-Lagrangian advective schemes in meteorological modeling; computation with high-resolution upwind schemes for hyperbolic equations; dynamics of flame propagation in a turbulent field; a modified finite element method for solving the incompressible Navier-Stokes equations; computational fusion magnetohydrodynamics; and a nonoscillatory shock capturing scheme using flux-limited dissipation. Consideration is also given to the use of spectral techniques in numerical weather prediction; numerical methods for the incorporation of mountains in atmospheric models; techniques for the numerical simulation of large-scale eddies in geophysical fluid dynamics; high-resolution TVD schemes using flux limiters; upwind-difference methods for aerodynamic problems governed by the Euler equations; and an MHD model of the earth's magnetosphere.

Engquist, B. E.↗
NASA NTRS2020

Implicit Large-Eddy Simulations of Hot and Cold Supersonic Jets in Loci-CHEM

Turbulence is a critical factor in most aero- and propulsion-related flows. For many years, the industry-standard approach to incorporating the effect of turbulence in computational fluid dynamics (CFD)calculations has been the Reynolds Averaged Navier Stokes (RANS) modeling approach. RANS implicitly time-averages turbulent motion, and models its effect on the mixing of species, momentum, and energy in a flow field. The approach is economical and often allows a CFD calculation to proceed to steady-state. However, it makes several simplifying assumptions, and the results will only be as good as the RANS turbulence model. Meanwhile, the last two decades have witnessed rapid and explosive growth in computational power, the development of CFD codes which scale well over thousands of processors, and the refinement of unstructured grid-generation tools which facilitate rapid surface and volume gridding of complex geometries. CFD engineering calculations of 108– 109 finite-volume cells have become routine. Thus, large eddy simulation (LES), which attempts to directly resolve the unsteady motion of the largest scales of turbulence, is increasingly of interest for many fluids engineering problems. Additionally, hybrid RANS-LES approaches such as Spalart’s Detached Eddy Simulation (DES), which attempt to marry the best strengths of both methods, are in active use and development. Because it seeks to accurately resolve the unsteady motion of the largest scales of turbulence, LES places certain requirements on the numerical method used to discretize and advance the Navier Stokes equations. In particular, LES requires high time accuracy and minimal intrinsic numerical dispersion and dissipation over a wide range of length scales.

C I Morris↗
DOE OSTI2021

Monte Carlo Simulations in Support of NorthStar Irradiation Facility Commissioning Procedure

NorthStar Medical Radioisotopes LLC is constructing an irradiation facility at which electron accelerators will induce photo-transmutation (i.e., neutron knockout) of molybdenum-100 (Mo-100) to produce medically useful Mo-99. The facility will house two high-intensity electron RhodotronTM accelerators developed by IBA Industrial. The nominal power of each Rhodotron is 125kW at 40 MeV beam energy. The production target is made of metallic molybdenum enriched with Mo-100 isotope. It will be irradiated by electron beams and will itself play a role in electron conversion to intense X-rays. A full description of the NorthStar target as well as the concept of radiation shielding of the facility can be found in our previous report. The NorthStar irradiation facility is under commissioning now. The goal of the present work is to support the commissioning process by: (1) Studying residual radiation and activation to minimize personnel exposure and inform radiation waste management; (2) Using computational fluid dynamics (CFD) of the irradiated target to ensure the facility design is safe and reliable; and (3) Helping to determine radiation protection requirements during operation. Monte Carlo simulations are used for the analyses. Most of the results are presented as 3D arrays representing spatial distribution of values across simulated geometries: radiation energy depositions, ambient doses, residual isotopes accumulated, and residual doses.

07 ISOTOPE AND RADIATION SOURCES↗
NASA NTRS2019

Parallel Anisotropic Unstructured Grid Adaptation

Computational Fluid Dynamics (CFD) has become critical to the design and analysis of aerospace vehicles. Parallel grid adaptation that resolves multiple scales with anisotropy is identified as one of the challenges in the CFD Vision 2030 Study to increase the capacity and capability of CFD simulation. The Study also cautions that computer architectures are undergoing a radical change and dramatic increases in algorithm concurrency will be required to exploit full performance. This paper reviews four different methods to parallel anisotropic grid generation. They cover both ends of the spectrum: (i) using existing state-of-the-art software optimized for a single core and modifying it for parallel platforms and (ii) designing and implementing scalable software with incomplete, but rapidly maturating functionality. A brief overview for each grid adaptation system is presented in the context of a telescopic approach for multilevel concurrency. These methods employ different approaches to enable parallel execution, which provides a unique opportunity to illustrate the relative behavior of each approach. Qualitative and quantitative metric evaluations are used to draw lessons for future developments in this critical area for parallel CFD simulation.

Tsolakis, Christos↗
NASA NTRS1998

A Geometry Based Infra-Structure for Computational Analysis and Design

The computational steps traditionally taken for most engineering analysis suites (computational fluid dynamics (CFD), structural analysis, heat transfer and etc.) are: (1) Surface Generation -- usually by employing a Computer Assisted Design (CAD) system; (2) Grid Generation -- preparing the volume for the simulation; (3) Flow Solver -- producing the results at the specified operational point; (4) Post-processing Visualization -- interactively attempting to understand the results. For structural analysis, integrated systems can be obtained from a number of commercial vendors. These vendors couple directly to a number of CAD systems and are executed from within the CAD Graphical User Interface (GUI). It should be noted that the structural analysis problem is more tractable than CFD; there are fewer mesh topologies used and the grids are not as fine (this problem space does not have the length scaling issues of fluids). For CFD, these steps have worked well in the past for simple steady-state simulations at the expense of much user interaction. The data was transmitted between phases via files. In most cases, the output from a CAD system could go to Initial Graphics Exchange Specification (IGES) or Standard Exchange Program (STEP) files. The output from Grid Generators and Solvers do not really have standards though there are a couple of file formats that can be used for a subset of the gridding (i.e. PLOT3D data formats). The user would have to patch up the data or translate from one format to another to move to the next step. Sometimes this could take days. Specifically the problems with this procedure are:(1) File based -- Information flows from one step to the next via data files with formats specified for that procedure. File standards, when they exist, are wholly inadequate. For example, geometry from CAD systems (transmitted via IGES files) is defined as disjoint surfaces and curves (as well as masses of other information of no interest for the Grid Generator). This is particularly onerous for modern CAD systems based on solid modeling. The part was a proper solid and in the translation to IGES has lost this important characteristic. STEP is another standard for CAD data that exists and supports the concept of a solid. The problem with STEP is that a solid modeling geometry kernel is required to query and manipulate the data within this type of file. (2) 'Good' Geometry. A bottleneck in getting results from a solver is the construction of proper geometry to be fed to the grid generator. With 'good' geometry a grid can be constructed in tens of minutes (even with a complex configuration) using unstructured techniques. Adroit multi-block methods are not far behind. This means that a million node steady-state solution can be computed on the order of hours (using current high performance computers) starting from this 'good' geometry. Unfortunately, the geometry usually transmitted from the CAD system is not 'good' in the grid generator sense. The grid generator needs smooth closed solid geometry. It can take a week (or more) of interaction with the CAD output (sometimes by hand) before the process can begin. One way Communication. (3) One-way Communication -- All information travels on from one phase to the next. This makes procedures like node adaptation difficult when attempting to add or move nodes that sit on bounding surfaces (when the actual surface data has been lost after the grid generation phase). Until this process can be automated, more complex problems such as multi-disciplinary analysis or using the above procedure for design becomes prohibitive. There is also no way to easily deal with this system in a modular manner. One can only replace the grid generator, for example, if the software reads and writes the same files. Instead of the serial approach to analysis as described above, CAPRI takes a geometry centric approach. This makes the actual geometry (not a discretized version) accessible to all phases of the analysis. The connection to the geometry is made through an Application Programming Interface (API) and NOT a file system. This API isolates the top-level applications (grid generators, solvers and visualization components) from the geometry engine. Also this allows the replacement of one geometry kernel with another, without effecting these top-level applications. For example, if UniGraphics is used as the CAD package then Parasolid (UG's own geometry engine) can be used for all geometric queries so that no solid geometry information is lost in a translation. This is much better than STEP because when the data is queried, the same software is executed as used in the CAD system. Therefore, one analyzes the exact part that is in the CAD system. CAPRI uses the same idea as the commercial structural analysis codes but does not specify control. Software components of the CAD system are used, but the analysis suite, not the CAD operator, specifies the control of the software session. This also means that the license issues (may be) minimized and individuals need not have to know how to operate a CAD system in order to run the suite.

Haimes, Robert↗
NASA NTRS1985

Coolant passage heat transfer with rotation

The objective is to develop a heat transfer and pressure drop data base, computational fluid dynamic techniques, and correlations for multi-pass rotating coolant passages with and without flow turbulators. The experimental effort is focused on the simulation of configurations and conditions expected in the blades of advanced aircraft high pressure turbines. With the use of this data base, the effects of Coriolis and buoyancy forces on the coolant side flow can be included in the design of turbine blades.

Hajek, T. J.↗