Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “CFD simulations”

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 253 records · Page 14

NASA ERA Integrated CFD for Wind Tunnel Testing of Hybrid Wing-Body Configuration

NASAs Environmentally Responsible Aviation (ERA) Project explores enabling technologies to reduce aviations impact on the environment. One research challenge area for the project has been to study advanced airframe and engine integration concepts to reduce community noise and fuel burn. In order to achieve this, complex wind tunnel experiments at both the NASA Langley Research Centers (LaRC) 14x22 and the Ames Research Centers 40x80 low-speed wind tunnel facilities were conducted on a Boeing Hybrid Wing Body (HWB) configuration. These wind tunnel tests entailed various entries to evaluate the propulsion airframe interference effects including aerodynamic performance and aeroacoustics. In order to assist these tests in producing high quality data with minimal hardware interference, extensive Computational Fluid Dynamic (CFD) simulations were performed for everything from sting design and placement for both the wing body and powered ejector nacelle systems to the placement of aeroacoustic arrays to minimize its impact on the vehicles aerodynamics. This paper will provide a high level summary of the CFD simulations that NASA performed in support of the model integration hardware design as well as some simulation guideline development based on post-test aerodynamic data. In addition, the paper includes details on how multiple CFD codes (OVERFLOW, STAR-CCM+, USM3D, and FUN3D) were efficiently used to provide timely insight into the wind tunnel experimental setup and execution.

Computational Fluid Dynamics↗

Aerodynamic Analysis of Multistage Turbomachinery Flows in Support of Aerodynamic Design

This paper summarizes the state of 3D CFD based models of the time average flow field within axial flow multistage turbomachines. Emphasis is placed on models which are compatible with the industrial design environment and those models which offer the potential of providing credible results at both design and off-design operating conditions. The need to develop models which are free of aerodynamic input from semi-empirical design systems is stressed. The accuracy of such models is shown to be dependent upon their ability to account for the unsteady flow environment in multistage turbomachinery. The relevant flow physics associated with some of the unsteady flow processes present in axial flow multistage machinery are presented along with procedures which can be used to account for them in 3D CFD simulations. Sample results are presented for both axial flow compressors and axial flow turbines which help to illustrate the enhanced predictive capabilities afforded by including these procedures in 3D CFD simulations. Finally, suggestions are given for future work on the development of time average flow models.

Adamczyk, John J.↗

Multi-Fidelity Heating Prediction of Adaptable, Deployable Entry Placement Technology Vehicles

The objective of this work was to investigate a multi-fidelity modeling approach to accurately and efficiently predict the laminar and turbulent convective heating on adaptable, deployable entry placement technology vehicles in Mars entry. A previously developed co-Kriging based multi-fidelity modeling approach was used to model the laminar and turbulent convective heat fluxes at several surface locations along the vehicle, including the rib sections. The laminar convective heat flux multi-fidelity model was found to have a mean convective heat rate error of approximately 3% when compared to high-fidelity CFD simulations. The turbulent convective heat flux multi-fidelity model was found to have a mean convective heat rate error of approximately 8% when compared to high-fidelity CFD simulations. Compared to a single fidelity model, the multi-fidelity model required approximately one-third the number of high-fidelity model evaluations to obtain the same accuracy level. The computational cost of evaluating the multi-fidelity model was approximately five orders of magnitude less than one high-fidelity model simulation.

Mario Santos↗

Tracking Blade Tip Vortices for Numerical Flow Simulations of Hovering Rotorcraft

Blade tip vortices generated by a helicopter rotor blade are a major source of rotor noise and airframe vibration. This occurs when a vortex passes closely by, and interacts with, a rotor blade. The accurate prediction of Blade Vortex Interaction (BVI) continues to be a challenge for Computational Fluid Dynamics (CFD). Though considerable research has been devoted to BVI noise reduction and experimental techniques for measuring the blade tip vortices in a wind tunnel, there are only a handful of post-processing tools available for extracting vortex core lines from CFD simulation data. In order to calculate the vortex core radius, most of these tools require the user to manually select a vortex core to perform the calculation. Furthermore, none of them provide the capability to track the growth of a vortex core, which is a measure of how quickly the vortex diffuses over time. This paper introduces an automated approach for tracking the core growth of a blade tip vortex from CFD simulations of rotorcraft in hover. The proposed approach offers an effective method for the quantification and visualization of blade tip vortices in helicopter rotor wakes. Keywords: vortex core, feature extraction, CFD, numerical flow visualization

Numerical↗

Unified Material-Environment Interaction Model for Binary UHTC Composites

A new model for the surface chemistry of zirconium and hafnium diborides containing silicon-carbide when exposed to high temperature conditions in air is formulated. Oxidation dynamics are simplified by considering the diffusion-limited equilibrium regime within the material at elevated temperatures. This model is first assessed against experimental mass-gain data taken from UHTC samples heated to high temperature in an oxygen environment. Then, further evaluations of the new model are conducted using hypersonic CFD simulations to analyze conditions experienced by HfB2-SiC in arc jet experiments conducted at NASA Ames Research Center. Measurements of stagnation point heat transfer and pressure are then used to calibrate the simulations. Results of coupled CFD-material response simulations are then compared to the surface temperatures measured during the Ames arc jet test. Results support the use of the proposed model while highlighting the need for improved gas-phase data on the additional constituents considered.

ultra high temperature ceramic↗

Gas-Particle Interaction Model Development in Plume Surface Interaction Erosion and Cratering

The development of a predictive simulation capability for extra-terrestrial Plume Surface Interaction (PSI) environments is undertaken by the Fluid Dynamics Branch at NASA/MSFC under the Game Changing Development (GCD) Program funded by NASA’s Space Technology Mission Directorate (STMD). Predictive simulation capability for propulsive Martian and Lunar landing PSI requires accurate modeling of the complex high-speed plume impingement and resulting gas/particle-cloud and ejecta stream formation. Physics modeling gap analysis during formulation of the PSI project scope identified two particular phenomena of first order importance in gas-particle interactions that lacked existing modeling formulations. The first arises from the lack of models capturing the unsteady drag forces imparted on particles by the rapidly expanding supersonic compressible flow. The second addresses the fluctuating forces and dispersions imparted on both the particle and the fluid resulting from the interference and wake turbulence generated by close proximity particle clouds, dubbed the Particle Turbulent Kinetic Energy(PTKE). Their first order significance has been identified in experiments, but simulation models currently do not exist for either effect. The development of models and the eventual implementation into the Eulerian Gas-Granular Flow Solver (GGFS) simulation tools was constructed as one of four major tasks of the PSI project. In this process, data on particle kinetics and gas-particle interactions are collected from carefully designed experiments of particles embedded in jets. The effects observed in the experiments are then replicated in high-fidelity particle-resolved CFD simulations to inform the formulation of improvements to particle phase drag models for implementations in the more efficient Eulerian-Lagrangian CFD simulations. The resulting models are ultimately ported to the Eulerian-Eulerian models applied for most efficient simulations in PSI production application tools. This paper describes the activities and accomplishments of the past year in the gas-particle interaction modeling task of the PSI project.

Jeff West↗

Development and Validation of a Slosh-Based Ullage Collapse Model

Ullage collapse modeling is an essential tool for the success of spacecraft powered by cryogenic propellants. The degree of ullage collapse can have a major impact on the design of cryogenic propellant tanks as well as the mission’s concept of operations. Additional commodities (e.g., helium gas for pressure control) and their quantities are dependent on the amount of heat transfer and phase change of propellants expected during flight. Currently, high-fidelity computational fluid dynamics (CFD) with heat transfer and phase change is the most accurate method for predicting ullage collapse. CFD simulations of this kind, however, can become very costly for problems requiring large domains, long simulation times, or significant liquid motion. Consequently, there is a need for a reduced order model to provide fast responses which can be used in engineering design decisions. The proposed numerical model is a nodal approach, augmented by outputs from CFD analyses that are generated rapidly and reliably. These outputs include transient values for liquid-vapor interface area, interface velocity, and liquid position using a volume of fluid (VoF) method. There are several benefits of combining a nodal approach with slosh predictions from CFD simulations: leveraging of the Loci/STREAM-VoF solver that has been validated in high-g and low-g environments, rapid modeling of VoF with isothermal propellant simplifications and relatively coarse meshes, and improved model predictions from the nodal tool. The nodal approach, informed by the slosh predictions from CFD, assumes uniform conditions in the liquid and gas domains. Modifications to account for temperature stratification in the liquid and gas may be implemented as a future improvement. Thermodynamic changes within the tank are solved at user-defined timesteps. First, impacts due to mass entering or leaving the control volume are captured. Second, heat transfer impacts are accounted for by tracking the volume of the gas which has been contacted by liquid in that timestep and rebalancing the internal energy. Finally, phase change is approximated using the area and velocity of the interface. The tank conditions are re-computed by conserving mass and energy. Calibration of the model to several test flights and experiments are required to baseline the predictive power and uncertainty. The final product will enable designers and engineers to create quick-turnaround predictions for ullage collapse and inform design feasibility.

Computational Fluid Dynamics↗

Vortex Shedding Inside a Baffled Air Duct

Common in the operation of both segmented and un-segmented large solid rocket motors is the occurrence of vortex shedding within the motor chamber. A portion of the energy within a shed vortex is converted to acoustic energy, potentially driving the longitudinal acoustic modes of the motor in a quasi-discrete fashion. This vortex shedding-acoustic mode excitation event occurs for every Reusable Solid Rocket Motor (RSRM) operation, giving rise to subsequent axial thrust oscillations. In order to better understand this vortex shedding/acoustic mode excitation phenomena, unsteady CFD simulations were run for both a test geometry and the full scale RSRM geometry. This paper covers the results from the subscale geometry runs, which were based on work focusing on the RSRM hydrodynamics. Unsteady CFD simulation parameters, including boundary conditions and post-processing returns, are reviewed. The results were further post-processed to identify active acoustic modes and vortex shedding characteristics. Probable locations for acoustic energy generation, and subsequent acoustic mode excitation, are discussed.

Davis, Philip↗

Reduced-Order CFD Modeling to Support Waste Loading Optimization in Hanford WTP Vitrification

The U.S. DOE Hanford Site stores over 56 million gallons of radioactive liquid tank waste that must be treated and immobilized for long-term disposal The Waste Treatment and Immobilization Plant (WTP) will vitrify this waste by feeding it into Joule-heated melters, where it is incorporated into a stable borosilicate glass Computational fluid dynamics (CFD) simulations of glass melters can provide insight into the maximum achievable waste loading under varying melter operating conditions Fully resolved VOF multiphase simulations were used as the reference model to capture bubble-driven convection in the melter, including bubble formation, rise behavior, and induced glass melt circulation Effective bubble column diameter and rise velocity were extracted from the resolved simulations, compared with empirical correlations, and refit across relevant viscosity and gas flow rate conditions Explicit gas–liquid interface tracking was replaced with a single-phase momentum source term model, enabling faster steady-state CFD simulations while preserving the dominant hydrodynamic effects of bubbling New empirical correlations were developed for effective bubble column diameter and bubble rise velocity by fitting resolved simulation data across expected melter viscosity and gas flow rate ranges, providing improved inputs for the momentum source term model compared with existing literature correlations The momentum source term model reduced fluid-domain mesh size by 89% and achieved an 8.4× computational speedup relative to resolved bubbling simulations The validated momentum source term approach enables prediction of process-relevant heat transfer behavior in the integrated melter model, including heat transfer from the molten glass to the cold cap, plenum, refractory walls, and surrounding structural regions under varying melter operating conditions

12 - MGMT OF RADIOACTIVE AND NON-RADIOACTIVE WASTE↗

Isothermal Compressor Computational Fluid Dynamics Simulations (Final Report)

Carnot Compression is a startup company developing an innovative technology for air and gas compression. This technology is inherently oil-free and isothermal due to the use of water to simultaneously compress and cool the gas throughout the process. Isothermal gas compression eliminates the need to cool the compressed gas so it may lead to significant energy savings. The development of Carnot’s isothermal compressor is limited by the lack of insight into the flow and detailed behavior of the fluids (water and air/gas) inside the air end. Previous and current attempts at Computational Fluid Dynamics (CFD) simulations by Carnot have been unable to provide the level of accuracy required to use simulations for technology development. In this Cooperative Research and Development Agreement (CRADA) project, Oak Ridge National Laboratory (ORNL) used the CFD package Star-CCM+ to successfully develop a CFD simulation, providing the much-needed insight required to speed up the development. The work is intended to enable Carnot Compression to unlock its technology potential to a level sufficient for commercialization of the intended first product. The work will also prepare Carnot to scale the technology to much larger and more energy intensive applications, positioning it to broaden the product applications. Successful development of the technology has the potential to result in 20% or more efficiency gains for compression processes, or about $3 billion in annual energy savings potential for the U.S. alone.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Airfoil Computational Fluid Dynamics - 2k shapes, 25 AoA's, 3 Re numbers

This dataset contains aerodynamic quantities - including flow field values (momentum, energy, and vorticity) and summary values (coefficients of lift, drag, and momentum) - for 1,830 airfoil shapes computed using the HAM2D CFD (computational fluid dynamics) model. The airfoil shapes were designed using the separable shape tensor parameterization that encodes two-dimensional shapes as elements of the Grassmann manifold. This data-driven approach learns two independent spaces of parameter from a collection of sample airfoils. The first captures large-scale, linear perturbations, and the second defines small-scale, higher-order perturbations. For this dataset, we used the G2Aero database of over 19,000 airfoil shapes to learn a parameter space that captured a wide array of shape characteristics. We sampled airfoil designs over both parameter spaces to explore the full range of possible shape variations. The aerodynamic quantities for the generated airfoil were obtained using the HAM2D code, which is a finite-volume Reynolds-averaged Navier-Stokes (RANS) flow solver. We employ a fifth-order WENO scheme for spatial reconstruction with Roe's flux difference scheme for inviscid flux and second-order central differencing for viscous flux. A preconditioned GMRES method is applied for implicit integration. The Spalart-Allmaras 1-eq turbulence model is used for the turbulence closure, and the Medida-Baeder 2-eq transition model is applied to account for the effects of laminar turbulent transition. The airfoil grid is generated with a total of 400 points on the airfoil surface, the initial wall-normal spacing of y+ = 1, and an outer boundary located at 300 chord lengths away from the wall. The CFD simulations are performed at a freestream Mach number of 0.1, for or three different Reynolds' numbers (3M, 6M, and 9M), and for 25 angles of attack from -4 deg. to 20 deg. with 1 degree increments. Across all these various parameters, this dataset includes the results from over 250,000 CFD simulations. The simulations were performed using the Bridges-2 system at the Pittsburgh Supercomputing Center in February 2023 as part of the INTEGRATE project funded by the Advanced Research Projects Agency - Energy, in the U.S. Department of Energy. The data was collected, reformatted, and preprocessed for this OEDI submission in July 2023 under the Foundational AI for Wind Energy project funded by the U.S. Department of Energy Wind Energy Technologies Office. This dataset is intended to serve as a benchmark against which new artificial intelligence (AI) or machine learning (ML) tools may be tested. Baseline AI/ML methods for analyzing this dataset have been implemented, and a link to their repository containing those models has been provided. The .h5 data file structure can be found in the GitHub Repository resource under explore_airfoil_2k_data.ipynb.

2k↗

Dynamic separation of gases using microsieves

Separation of light weight molecules, such as nitrogen, argon, and oxygen, from heavier compounds can have significant impacts on energy capture, environmental monitoring, or isotopic applications. Large-scale gas separation techniques, like gas centrifugation and membrane mitigation, can be problematic as they impart tremendous energy and induce high mechanical stress onto the instrumentation. Microsieves, also known as micronozzles or microfunnels, are developed to create physical barriers to separate specific isotopes and gases. Separation is achieved using a converging and diverging micronozzle to impose supersonic gas flow around a curved wall, and it has been used for the separation of heavy actinide isotopes in low weight gas as well as separation of low weight gas compositions of nitrogen and argon back in 1900s. However, systematic reviews of this unique technology are lacking. The application of the Laval style nozzle, which has a converging/diverging entrance fundamental to the micronozzle, is included in this review due to its importance in industrial applications in uranium (U) isotope refinement. Using advanced computational fluid dynamic (CFD) simulations, the extent of gas separation can be modelled. Herein, we first examine the literature and survey recent advances on fabrication techniques for creating curved micronozzles, methods and separation principles used to design devices. Furthermore, we then follow with highlights of CFD simulations applied to evaluate the separation effects using microsieves. Finally, identification of the gap and recommendation for future development and applications are suggested for using intrinsic molecular features and fluidic dynamics in formulating separation strategies.

30 Microfluidics↗

Understanding the Flow Physics of Shock Boundary-Layer Interactions Using CFD and Numerical Analyses

Computational fluid dynamic (CFD) analyses of the University of Michigan (UM) Shock/Boundary-Layer Interaction (SBLI) experiments were performed as an extension of the CFD SBLI Workshop held at the 48th AIAA Aerospace Sciences Meeting in 2010. In particular, the UM Mach 2.75 Glass Tunnel with a semi-spanning 7.75deg wedge was analyzed in attempts to explore key physics pertinent to SBLI's, including thermodynamic and viscous boundary conditions as well as turbulence modeling. Most of the analyses were 3D CFD simulations using the OVERFLOW flow solver, with additional quasi-1D simulations performed with an in house MATLAB code interfacing with the NIST REFPROP code to explore perfect verses non-ideal air. A fundamental exploration pertaining to the effects of particle image velocimetry (PIV) on post-processing data is also shown. Results from the CFD simulations showed an improvement in agreement with experimental data with key contributions including adding a laminar zone upstream of the wedge and the necessity of mimicking PIV particle lag for comparisons. Results from the quasi-1D simulation showed that there was little difference between perfect and non-ideal air for the configuration presented.

Friedlander, David J.↗

Validation of High-Resolution CFD Method for Slosh Damping Extraction of Baffled Tanks

Determination of slosh damping is a very challenging task as there is no analytical solution. The damping physics involve the vorticity dissipation which requires the full solution of the nonlinear Navier-Stokes equations. As a result, previous investigations and knowledge were mainly carried out by extensive experimental studies. A Volume-Of-Fluid (VOF) based CFD program developed at NASA MSFC was applied to extract slosh damping in a baffled tank from the first principle. First, experimental data using water with subscale smooth wall tank were used as the baseline validation. CFD simulation was demonstrated to be capable of accurately predicting natural frequency and very low damping value from the smooth wall tank at different fill levels. The damping due to a ring baffle at different liquid fill levels from barrel section and into the upper dome was then investigated to understand the slosh damping physics due to the presence of a ring baffle. Based on this study, the Root-Mean-Square error of our CFD simulation in estimating slosh damping was less than 4.8%, and the maximum error was less than 8.5%. Scalability of subscale baffled tank test using water was investigated using the validated CFD tool, and it was found that unlike the smooth wall case, slosh damping with baffle is almost independent of the working fluid and it is reasonable to apply water test data to the full scale LOX tank when the damping from baffle is dominant. On the other hand, for the smooth wall, the damping value must be scaled according to the Reynolds number. Comparison of experimental data, CFD, with the classical and modified Miles equations for upper dome was made, and the limitations of these semi-empirical equations were identified.

Yang, H. Q.↗

Transonic Correction to Theodorsen's Theory for Oscillating Airfoil in Pitch and Plunge Toward Flutter

This paper presents a transonic correction method for an oscillating airfoil in pitch and plunge. The proposed method applies correction functions to the Theodorsen’s theory to capture the transonic nonlinear aero- dynamics. These correction functions apply necessary corrections to the amplitudes and the phase angles of the unsteady lift and pitching moment coefficients to account for transonic aerodynamics. The proposed method also postulates a correction for the motion of the aerodynamic center which could be induced by moving shocks. A series of unsteady RANS CFD simulations of the airfoil at the mean aerodynamic chord of the Transonic Truss-Braced Wing aircraft are conducted using FUN3D to provide data to construct these transonic correction functions. The computed responses of the unsteady lift and pitching moment coefficients using these transonic correction functions match the CFD simulation results very well even when the pitching moment coefficient is highly nonlinear. A flutter analysis of an airfoil in pitch and plunge illustrates the potential use of the proposed transonic correction method.

Transonic Flutter↗

Characterization of Slosh Damping for Ortho-Grid and Iso-Grid Internal Tank Structures

Grid stiffened tank structures such as Ortho-Grid and Iso-Grid are widely used in cryogenic tanks for providing stiffening to the tank while reducing mass, compared to tank walls of constant cross-section. If the structure is internal to the tank, it will positively affect the fluid dynamic behavior of the liquid propellant, in regard to fluid slosh damping. As NASA and commercial companies endeavor to explore the solar system, vehicles will by necessity become more mass efficient, and design margin will be reduced where possible. Therefore, if the damping characteristics of the Ortho-Grid and Iso-Grid structure is understood, their positive damping effect can be taken into account in the systems design process. Historically, damping by internal structures has been characterized by rules of thumb and for Ortho-Grid, empirical design tools intended for slosh baffles of much larger cross-section have been used. There is little or no information available to characterize the slosh behavior of Iso-Grid internal structure. Therefore, to take advantage of these structures for their positive damping effects, there is much need for obtaining additional data and tools to characterize them. Recently, the NASA Marshall Space Flight Center conducted both sub-scale testing and computational fluid dynamics (CFD) simulations of slosh damping for Ortho-Grid and Iso-Grid tanks for cylindrical tanks containing water. Enhanced grid meshing techniques were applied to the geometrically detailed and complex Ortho-Grid and Iso-Grid structures. The Loci-STREAM CFD program with the Volume of Fluid Method module for tracking and locating the water-air fluid interface was used to conduct the simulations. The CFD simulations were validated with the test data and new empirical models for predicting damping and frequency of Ortho-Grid and Iso-Grid structures were generated.

Douglas G. Westra↗

A summary of computational experience at GE Aircraft Engines for complex turbulent flows in gas turbines

This viewgraph presentation summarizes some CFD experience at GE Aircraft Engines for flows in the primary gaspath of a gas turbine engine and in turbine blade cooling passages. It is concluded that application of the standard k-epsilon turbulence model with wall functions is not adequate for accurate CFD simulation of aerodynamic performance and heat transfer in the primary gas path of a gas turbine engine. New models are required in the near-wall region which include more physics than wall functions. The two-layer modeling approach appears attractive because of its computational complexity. In addition, improved CFD simulation of film cooling and turbine blade internal cooling passages will require anisotropic turbulence models. New turbulence models must be practical in order to have a significant impact on the engine design process. A coordinated turbulence modeling effort between NASA centers would be beneficial to the gas turbine industry.

Zerkle, Ronald D.↗

Numerical Simulations of an Unsteady Rocket Launch from the AH-64D Apache Longbow Helicopter

Rocket and missile firings from attack helicopters can cause main engine compressor stall. Studies of this phenomenon suggest that the main engine ingests either the plume from the rockets or the rocket blast waves. This creates surges at the inlet face, causing a loss of power in the main engine. The objective of this project is to set-up a computational fluid dynamics (CFD) simulation of the AH-64D Apache Longbow helicopter during a rocket launch, in order to qualitatively study the fluid dynamics of the problem. This project presents a progression of three unsteady Navier-Stokes solutions. The first unsteady solution involves only a rocket launch from its launch canister. The second solution is a launch from a canister mounted on the Apache's wing-pylon assembly. The last solution includes the Apache main engine and fuselage. The computations use a series of structured, overset grid systems, which allow for a rocket moving in a prescribed path. The method implements a Roe upwind scheme with LU-SGS (lower-upper factored symmetric Gauss-Seidel). A rotor pressure disk model approximates the helicopter rotor, while the rocket engine exit properties are applied as a prescribed boundary condition. Although the project is only at the half-way point, the first and second CFD simulations suggest the possibility of pressure wave interference. Sudden surges in pressure occur from two sources: at rocket start-up, and as the rocket leaves the canister. Wave patterns set-up by these sources appear to propagate to the location of the engine inlet. However the simplified geometry simulation with the main engine needs to be performed before coming to a conclusion.

Okamoto, Kevin↗