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,315 records · Page 73
All sourcesNASA NTRSDOE OSTI
NASA NTRS2010

Slot Nozzle Effects for Reduced Sonic Boom on a Generic Supersonic Wing Section

NASA has conducted research programs to reduce or eliminate the operational restrictions of supersonic aircraft over populated areas. Restrictions are due to the disturbance from the sonic boom, caused by the coalescence of shock waves formed off the aircraft. Results from two-dimensional computational fluid dynamic (CFD) analyses (performed on a baseline Mach 2.0 nozzle in a simulated Mach 2.2 flow) indicate that over-expanded and under-expanded operation of the nozzle has an effect on the N-wave boom signature. Analyses demonstrate the feasibility of reducing the magnitude of the sonic boom N-wave by controlling the nozzle plume interaction with the nozzle boat tail shock structure. This work was extended to study the impact of integrating a high aspect ratio exhaust nozzle or long slot nozzle on the trailing edge of a supersonic wing. The nozzle is operated in a highly under-expanded condition, creating a large exhaust plume and a shock at the trailing edge of the wing. This shock interacts with and suppresses the expansion wave caused by the wing, a major contributor to the sonic boom signature. The goal was to reduce the near field pressures caused by the expansion using a slot nozzle located at the wing trailing edge. Results from CFD analysis on a simulated wing cross-section and a slot nozzle indicate potential reductions in sonic boom signature compared to a baseline wing with no propulsion or trailing edge exhaust. Future studies could investigate if this effect could be useful on a supersonic aircraft for main propulsion, auxiliary propulsion, or flow control.

Caster, Raymond S.↗
NASA NTRS2010

CFD Fuel Slosh Modeling of Fluid-Structure Interaction in Spacecraft Propellant Tanks with Diaphragms

Liquid sloshing within spacecraft propellant tanks causes rapid energy dissipation at resonant modes, which can result in attitude destabilization of the vehicle. Identifying resonant slosh modes currently requires experimental testing and mechanical pendulum analogs to characterize the slosh dynamics. Computational Fluid Dynamics (CFD) techniques have recently been validated as an effective tool for simulating fuel slosh within free-surface propellant tanks. Propellant tanks often incorporate an internal flexible diaphragm to separate ullage and propellant which increases modeling complexity. A coupled fluid-structure CFD model is required to capture the damping effects of a flexible diaphragm on the propellant. ANSYS multidisciplinary engineering software employs a coupled solver for analyzing two-way Fluid Structure Interaction (FSI) cases such as the diaphragm propellant tank system. Slosh models generated by ANSYS software are validated by experimental lateral slosh test results. Accurate data correlation would produce an innovative technique for modeling fuel slosh within diaphragm tanks and provide an accurate and efficient tool for identifying resonant modes and the slosh dynamic response.

Sances, Dillon J.↗
NASA NTRS2017

Validation Process for LEWICE by Use of a Navier-Stokes Solver

A research project is underway at NASA Glenn to produce computer software that can accurately predict ice growth under any meteorological conditions for any aircraft surface. This report will present results from the latest LEWICE release, version 3.5. This program differs from previous releases in its ability to model mixed phase and ice crystal conditions such as those encountered inside an engine. It also has expanded capability to use structured grids and a new capability to use results from unstructured grid flow solvers. A quantitative comparison of the results against a database of ice shapes that have been generated in the NASA Glenn Icing Research Tunnel (IRT) has also been performed. This paper will extend the comparison of ice shapes between LEWICE 3.5 and experimental data from a previous paper. Comparisons of lift and drag are made between experimentally collected data from experimentally obtained ice shapes and simulated (CFD) data on simulated (LEWICE) ice shapes. Comparisons are also made between experimentally collected and simulated performance data on select experimental ice shapes to ensure the CFD solver, FUN3D, is valid within the flight regime. The results show that the predicted results are within the accuracy limits of the experimental data for the majority of cases.

large-eddy simulations↗
NASA NTRS2017

First benchmark of the Unstructured Grid Adaptation Working Group

Unstructured grid adaptation is a technology that holds the potential to improve the automation and accuracy of computational fluid dynamics and other computational disciplines. Difficulty producing the highly anisotropic elements necessary for simulation on complex curved geometries that satisfies a resolution request has limited this technology's widespread adoption. The Unstructured Grid Adaptation Working Group is an open gathering of researchers working on adapting simplicial meshes to conform to a metric field. Current members span a wide range of institutions including academia, industry, and national laboratories. The purpose of this group is to create a common basis for understanding and improving mesh adaptation. We present our first major contribution: a common set of benchmark cases, including input meshes and analytic metric specifications, that are publicly available to be used for evaluating any mesh adaptation code. We also present the results of several existing codes on these benchmark cases, to illustrate their utility in identifying key challenges common to all codes and important differences between available codes. Future directions are defined to expand this benchmark to mature the technology necessary to impact practical simulation workflows.

Ibanez, Daniel↗
NASA NTRS2019

ASPIRE Aerodynamic Models and Flight Performance

The Advanced Supersonic Parachute Inflation Research Experiments (ASPIRE) project waslaunched to develop the capability for testing supersonic parachutes at Mars-relevant conditions.Three initial parachute tests, targeted as a risk-reduction activity for NASA's upcoming Mars2020 mission, successfully tested two candidate parachute designs and provided valuable data onparachute inflation, forces, and aerodynamic behavior. Design of the flight tests depended onflight mechanics simulations which in turn required aerodynamic models for the payload, andthe parachute. Computational Fluid Dynamics (CFD) was used to generate these models preflightand they are compared against the flight data after the tests. For the payload, thereconstructed aerodynamic behavior is close to the pre-flight predictions, but the uncertaintiesin the reconstructed data are high due to the low dynamic pressures and accelerations duringthe flight period of comparison. For the parachute, the predicted time to inflation agrees wellwith the pre-flight model; the peak aerodynamic force and the steady state drag on the parachuteare within the bounds of the pre-flight models, even as the models over-predict the parachutedrag at supersonic Mach numbers. Notably, the flight data do not show the transonic dragdecrease predicted by the pre-flight model. The ASPIRE flight tests provide previouslyunavailable valuable data on the performance of a large full-scale parachute behind a slenderleading body at Mars-relevant Mach number, dynamic pressure and parachute loads. TheASPIRE results are used to propose a new model for the parachute drag behind slender bodiesto aid future experiments.

Muppidi, Suman↗
NASA NTRS2019

ASPIRE Aerodynamic Models and Flight Performance

The Advanced Supersonic Parachute Inflation Research Experiments (ASPIRE) project waslaunched to develop the capability for testing supersonic parachutes at Mars-relevant conditions.Three initial parachute tests, targeted as a risk-reduction activity for NASA's upcomingMars2020 mission, successfully tested two candidate parachute designs and provided valuabledata on parachute inflation, forces, and aerodynamic behavior. Design of the flight tests dependedon flight mechanics simulations which in turn required aerodynamic models for the payload, andthe parachute. Computational Fluid Dynamics (CFD) was used to generate these models preflightand are compared against the flight data after the tests. For the payload, the reconstructedaerodynamic behavior is close to the pre-flight predictions, but the uncertainties in thereconstructed data are high due to the low dynamic pressures and accelerations during the flightperiod of comparison. For the parachute, the predicted time to inflation agrees well with the preflightmodel; the peak aerodynamic force and the steady state drag on the parachute are withinthe bounds of the pre-flight models, even as the models over-predict the parachute drag atsupersonic Mach numbers. Notably, the flight data does not show the transonic drag decreasepredicted by the pre-flight model. The ASPIRE flight tests provide previously unavailablevaluable data on the performance of a large full-scale parachute behind a slender leading bodyat Mars-relevant Mach number, dynamic pressure and parachute loads. This data is used topropose a new model for the parachute drag behind slender bodies to aid future experiments.

Muppidi, Suman↗
NASA NTRSDate not supplied

Aerodynamic Implications of Aerocapture Systems for Uranus Orbiters

Exploration of the Ice Giants, Uranus in particular, has been indicated as the highest priority new flagship-class mission by the 2022 Planetary Science Decadal Survey. However, due to it’s location in the solar system, extensive travel time and resources are required for a mission to reach Uranus. Aerocapture offers shorter interplanetary transit times and reduced propellant requirements for such a mission when compared to a traditional propulsive-capture orbiter. Considerations for the aerodynamics of an aerocapture system enabling flagship-class science mission to Uranus are presented. The applicability of a low lift-to-drag entry vehicle aeroshell with flight heritage to a new planet and atmosphere is evaluated. A preliminary aerodatabase is formulated from low-fidelity methods and heritage flight mission data to provide initial estimates of trajectory space to be considered for inserting a spacecraft into orbit at Uranus using aerocapture. Computational fluid dynamics calculations in the free molecular, transitional, and continuum regimes using Direct Simulation Monte Carlo and Navier-Stokes solutions are used to assess the validity of this preliminary aerodatabase. The considered entry vehicle is found to varying levels of agreement with work done by previous flight missions when evaluated for flight in Uranus’ atmosphere. Implications for the aerodynamics of this vehicle as applied in an aerocapture system in flight at Uranus as found by the current work are discussed.

Eli R. Shellabarger↗
NASA NTRSDate not supplied

Aerodynamic Implications of Aerocapture Systems for Uranus Orbiters

Exploration of the Ice Giants, Uranus in particular, has been indicated as the highest priority new flagship-class mission by the 2022 Planetary Science Decadal Survey. However, due to it's location in the solar system, extensive travel time and resources are required for a mission to reach Uranus. Aerocapture offers shorter interplanetary transit times and reduced propellant requirements for such a mission when compared to a traditional propulsive-capture orbiter. Considerations for the aerodynamics of an aerocapture system enabling flagship-class science mission to Uranus are presented. The applicability of a low lift-to-drag entry vehicle aeroshell with flight heritage to a new planet and atmosphere is evaluated. A preliminary aerodatabase is formulated from low-fidelity methods and heritage flight mission data to provide initial estimates of trajectory space to be considered for inserting a spacecraft into orbit at Uranus using aerocapture. Computational fluid dynamics calculations in the free molecular, transitional, and continuum regimes using Direct Simulation Monte Carlo and Navier-Stokes solutions are used to assess the validity of this preliminary aerodatabase. The considered entry vehicle is found to varying levels of agreement with work done by previous flight missions when evaluated for flight in Uranus' atmosphere. Implications for the aerodynamics of this vehicle as applied in an aerocapture system in flight at Uranus as found by the current work are discussed.

Uranus↗
NASA NTRS2024

Particle Impact Simulation and Ignition Prediction

An experimentally calibrated tool is needed to predict if a system is susceptible to failure by particle impact ignition (PI) based on use conditions, materials, and flow geometry. This tool will accelerate new components, evaluating existing hardware, and help disposition anomalies. - Conduct particle impact testing with in-situ diagnostics and complementary simulations on subset of key engineering materials (IN718, M400, 316L, 6061, Ti64, Zr) to develop a proof-of-concept predictive tool for assessing the risk of PI for idealized geometries (spherical particles) in realistic environments. - Assess particle/target interactions (coefficient of restitution, ignition, kindling) using instrumented particle impact rigs while systematically varying key parameters (materials, particle size, environment, target configuration). - Determine key field variables (temperature, strain, stress) in particle impacts using Multiphysics finite element and hydrocode simulations validated through comparison with experimental measurements and observations. - Synthesize experiments and simulations into constitutive models for PI that can be integrated with existing computational fluid dynamics (CFD) and Debris Transport Analysis (DTA) tools in future efforts.

Jonathan Tylka↗
DOE OSTI2022

An Efficient Time-Domain Model to Simulate Parametric Resonances in a Floating Body Free to Move in Six Degrees of Freedom: Preprint

We present a computationally efficient time-domain model capable of simulating parametric resonances in a floating body in waves. The model assumes all wave forces to be linear, but the inertia and restoring forces acting on the body are expanded to second order in body motions. The simulation speed on a standard computer is approximately 40 times faster than real time. The model is applied to a soft-moored floating axisymmetric body which absorbs energy through heave, but is otherwise free to move in six degrees of freedom. Under certain conditions, we show that the body responds parametrically with large amplitudes not only in surge and pitch, but also in sway, roll, and yaw, provided it is given some small initial displacement in one of these out-of-plane modes. The predictions are confirmed by simulations using state-of-the-art nonlinear Froude-Krylov and computational fluid dynamics models.

parametric resonance↗
DOE OSTI2024

Membrane-based carbon capture process optimization using CFD modeling

Carbon capture is a promising option to mitigate CO2 emissions from existing coal-fired power plants, cement and steel industries, and petrochemical complexes. Among the available technologies, membrane-based carbon capture presents the lowest energy consumption, operating costs, and carbon footprint. In addition, membrane processes have important operational flexibility and response times. On the other hand, the major challenges to widespread application of this technology are related to reducing capital costs and improving membrane stability and durability. To upscale the technology into stacked flat sheet configurations, high fidelity computational fluid dynamics (CFD) that describes the separation process accurately are required. High fidelity simulations have been shown to be effective in studying the complex transport phenomena in membrane systems. In addition, obtaining high CO2 recovery percentages and product purity requires a multi-stage membrane process, where the optimal network configuration of the membrane modules must be studied in a systematic way. In order to address the design problem at process scale, we formulate a superstructure for the membrane-based carbon capture, including up to three separation stages. In the formulation of the optimization problem, we include reduced models, based on rigorous CFD simulations of the membrane modules. Numerical results indicate that the optimal design includes three membrane stages, and the capture cost is 45.4 $/t-CO2.

Pedrozo, Hector A.↗
DOE OSTI2022

The role of breakup and coalescence in fine-scale bubble-induced turbulence. I. Dynamics

We study the effect of bubble breakup and coalescence on fine-scale dynamics of bubbly turbulent flows using direct numerical simulations. We perform two different simulations of dilute bubbly flows of void fraction 0.5%: one with bubbles breaking up and coalescing and the other without these physical processes. The volume of the fluid method is used for simulating bubbles undergoing breakup and coalescence while the bubbles are treated as rigid spheres in the immersed boundary method simulation. The energy spectrum in both types of simulation, consistent with previous studies, exhibits a -3 slope. We follow a single infinitesimal fluid element as it evolves to understand velocity gradient dynamics using conditional mean trajectories. We note finite-time divergence when the fluid element evolves under the action of inertial and pressure dynamics. The inertial, pressure, and viscous velocity gradient dynamics, when considered individually, produce the same results for bubble-induced turbulence (BIT) as with the classical homogeneous isotropic turbulence (HIT). Yet when the overall velocity gradient dynamics is considered, BIT results in non-cyclic trajectories moving toward stable node and unstable saddle while classical HIT shows cyclic behavior in their trajectories that move toward the origin. Interestingly, both the volume of fluid and immersed boundary simulations produce similar results. Therefore, there are two main takeaways from this research. First, new velocity gradient models are needed for BIT as their velocity gradients behavior is entirely different from the HIT. Second, we can neglect the bubble topology, breakup, and coalescence while studying or modeling the fine-scale dynamics of BIT.

42 ENGINEERING↗
DOE OSTI2020

Hydrodynamic computations of high-power laser drives generating metal ejecta jets from surface grooves

Understanding dynamic fragmentation in shock-loaded metals and predicting properties of the resulting ejecta are of considerable importance for both basic and applied science. The nature of material ejection has been shown to change drastically when the free surface melts on compression or release. In this work, we present hydrodynamic simulations of laser-driven microjetting from micron-scale grooves on a tin surface. We study microjet formation across a range of shock strengths from drives that leave the target solid after release to drives that induce shock melting in the target. The shock-state particle velocity (Up) varies from 0.3 to 3 km/s and the shock breakout pressure is 3–120 GPa. The microjet tip velocity is 1–8 km/s and the free-surface velocity varies from 0.1 to 5 km/s. Two tin equations of state are examined: a “soft” model (LEOS 501) where the target melts for U p > 1 km/s and a more detailed multiphase model (SESAME 2161) that melts for U p > 1.4 km/s. We use these two models to examine the influence of phase change and the choice of the material model on microjet formation and evolution. We observe in our computational results that jet formation can be classified into three regimes: a low-energy regime where material strength affects jet formation, a moderate-energy regime dominated by the changing phase of tin material, and a high-energy regime where results are insensitive to the material model and jet formation is described by an idealized steady-jet theory. Using an ensemble of 2D simulations, we show that these trends hold across a wide range of drive energies and groove angles.

36 MATERIALS SCIENCE↗
NASA NTRS2013

CFD Lagrangian Modeling of Water Droplet Transport for ISS Hygiene Activity Application

The goal of this study was to assess the impacts of free water propagation in the Waste and Hygiene Compartment (WHC) installed in Node 3. Free water can be generated inside the WHC in small quantities due to crew hygiene activity. To mitigate potential impact of free water in Node 3 cabin the WHC doorway is enclosed by a waterproof bump-out, Kabin, with openings at the top and bottom. At the overhead side of the rack, there is a screen that prevents large drops of water from exiting. However, as the avionics fan in the WHC causes airflow toward the deck side of the rack, small quantities of free water may exit at the bottom of the Kabin. A Computational Fluid Dynamics (CFD) analysis of Node 3 cabin airflow enable identifying the paths of water transport. To simulate the droplet transport the Lagrangian discrete phase approach was used. Various initial droplet distributions were considered in the study. The droplet diameter was varied in the range of 5-20 mm. The results of the computations showed that most of the drops fall to the rack surface not far from the WHC curtain.

Son, Chang H.↗
NASA NTRSDate not supplied

Modeling Autogenous Pressurization and Draining of a Cryogenic Storage Tank in Normal Gravity

A two-phase CFD model for autogenous pressurization and draining of a cryogenic storage tank is presented using both the Sharp Interface and Volume-Of-Fluid (VOF) approaches for capturing the front and the associated interfacial heat, mass and momentum transfer between the liquid and vapor regions. Both models are validated against data provided by the Cryogenic Propellant Storage and Transfer (CPST) Engineering Development Unit (EDU) experiment. The results of the autogenous pressurization are presented first, focusing on the phase change and turbulence effects on the tank pressure and temperature predictions. Both the Sharp Interface (SI-CFD) and VOF (VOF-CFD) multiphase models predict tank pressure during pressurization within 3% of the measured values. The sensitivity of key physical and numerical parameters of the problem are tested using the Sharp Interface model. Effects of the accommodation coefficient (AC), and the computational grid structure are studied. The second part of this paper is devoted to validating the SI-VOF model, with an enhanced capability of moving the liquid-vapor interface, against the draining segment of the EDU experiment. The VOF-CFD model was also used to simulate tank draining and its results were compared with the results of the Sharp Interface model with the moving interface. Both models predict tank pressure during draining within 3.5% of the measured values. Pressure decrease rate is underpredicted by both models during the first 100 seconds of draining but matches the experimental rate for the rest of the simulation.

Computational Fluid Dynamics↗
NASA NTRSDate not supplied

Modeling Autogenous Pressurization and Draining of a Cryogenic Storage Tank in Normal Gravity

A two-phase CFD model for autogenous pressurization and draining of a cryogenic storage tank is presented using both the Sharp Interface and Volume-Of-Fluid (VOF) approaches for capturing the front and the associated interfacial heat, mass and momentum transfer between the liquid and vapor regions. Both models are validated against data provided by the Cryogenic Propellant Storage and Transfer (CPST) Engineering Development Unit (EDU) experiment1. The results of the autogenous pressurization are presented first, focusing on the phase change and turbulence effects on the tank pressure and temperature predictions. Both the Sharp Interface (SI-CFD) and VOF (VOF-CFD) multiphase models predict tank pressure during pressurization within 3% of the measured values. The sensitivity of key physical and numerical parameters of the problem are tested using the Sharp Interface model. Effects of the accommodation coefficient (AC), and the computational grid structure are studied. The second part of this paper is devoted to validating the SI-VOF model, with an enhanced capability of moving the liquid-vapor interface, against the draining segment of the EDU experiment. The VOF-CFD model was also used to simulate tank draining and its results were compared with the results of the Sharp Interface model with the moving interface. Both models predict tank pressure during draining within 3.5% of the measured values. Pressure decrease rate is underpredicted by both models during the first 100 seconds of draining but matches the experimental rate for the rest of the simulation.

Computational Fluid Dynamics↗
DOE OSTI2026

Two-Way Coupled CFD/DEM Simulation of Electromagnetic Heating in a Fluidized Bed Reactor

A set of tools has been developed to couple a multiphase reacting computational fluid dynamics solver (MFIX) with two computational electromagnetic solvers (COMSOL and Elmer). In contrast to several recent investigations, which calculate a single static EM field prior to the CFD calculation, in the present work the electromagnetic field is periodically updated during the simulation. This enables the EM field and associated EM heating predictions to reflect the effects of evolving material properties (temperature, composition, size etc.) and the evolving distribution of material in the system. The present work uses effective medium theory (EMT) for the electromagnetic field calculations since the goal is to apply the tool for industrial scale reacting systems. Predictions using MFIX/COMSOL and MFIX/Elmer are compared to in bed thermocouple measurements from a fluidized bed experiment and recent simulation predictions and recent experiments performed at NETL. The NETL experiments also include packed bed operating conditions. Coupled MFIX/COMSOL simulations at these conditions are compared with similar simulations using only COMSOL to measure errors due to the solver coupling approach.

CFD↗
DOE OSTI2023

Demonstration of RANS models with wall functions in the spectral element code Nek5000

The spectral element based computational fluid dynamics (CFD) code Nek5000 has been traditionally used for high-fidelity applications, such as direct numerical simulation (DNS) and large eddy simulation (LES). These techniques require very fine numerical resolution to accurately capture turbulent fluctuations which can be prohibitively expensive for users without access to leadership class computing facilities. For broader application and adoption, significant effort has been invested to develop Reynolds-averaged Navier–Stokes (RANS) capabilities in Nek5000. Here, this work presents details of the implementation and demonstration of the standard wall functions for the κ–τ model in Nek5000. Results using the wall-modeled approach are compared to a wall-resolved approach for cases with negligible pressure gradient, viz., channel flow, pipe flow and flow in a reactor subchannel. Results show reasonably good agreement between the two approaches for friction factor and Nusselt number. Some expected differences are identified near the wall. These cases demonstrate the potential for significant computational savings by using much coarser meshes for the wall-modeled approach, with only minor differences between the predicted result. Additionally, several Reynolds numbers up to 1,000,000 are demonstrated for pipe flow and predicted friction factors and Nusselt numbers compared well to available correlations, with the worst below 10%. As the Reynolds number is increased, better agreement is observed between the correlations and the wall-modeled approach. In addition, flow in a molten salt fast reactor (MSFR) core is considered which features an adverse pressure gradient and flow separation. It showcases the inability of standard wall functions to accurately predict flows with adverse pressure gradients. The results, however, match reasonably well in trend in regions of the flow where the boundary layer is attached. Ongoing research is dedicated to include a pressure gradient correction to wall functions to improve the accuracy of flows with separation or reattachment and adverse or favorable pressure gradients.

42 ENGINEERING↗