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 721 records · Page 40

Equivalent Source Method Applied to Launch Acoustic Simulations

Aeroacoustic simulations of the launch environment are described. A hybrid computational fluid dynamics (CFD)/computational aeroacoustic (CAA) approach is developed in order to accurately and efficiently predict the sound pressure level spectrum on the launch vehicle and surrounding structures. The high-fidelity CFD code LAVA (Launch Ascent and Vehicle Analysis), is used to generate pressure time history at select locations in the flow field. A 3D exterior Helmholtz solver is then used to iteratively determine a set of monopole sources which mimic the noise generating mechanisms identified by the CFD solver. The acoustic pressure field generated from the Helmholtz solver is then used to evaluate the sound pressure levels.

Housman, Jeffrey A.↗

Mitigation of Autoignition Due to Premixing in a Hypervelocity Flow Using Active Wall Cooling

Preinjection of fuel on the forebody of an airbreathing vehicle is a proposed method to gain access to hypervelocity flight Mach numbers. However, this creates the possibility of autoignition either near the wall or in the core of the flow, thereby consuming fuel prematurely as well as increasing the amount of pressure drag on the vehicle. The computational fluid dynamics code VULCAN was used to conduct three dimensional simulations of the reacting flow in the vicinity of hydrogen injectors on a flat plate at conditions relevant to a Mach 12 notional flight vehicle forebody to determine the location where autoignition occurs. Active wall cooling strategies were formulated and simulated in response to regions of autoignition. It was found that tangential film cooling using hydrogen or helium were both able to nearly or completely eliminate wall autoignition in the flow domain of interest.

Axdahl, Erik↗

Design of a Lunar Plume-Surface Interaction Measurement System

Plume-surface interactions between a rocket plume and the lunar surface will be studied in-situ during two of NASA’s upcoming Commercial Lunar Payload Services Program missions. The payload, Stereo Cameras for Lunar Plume-Surface Studies (SCALPSS), will employ a multiple-camera photogrammetry system to obtain stereo images of the lunar regolith before, during, and after significant erosion and subsequent crater formation. The evaluation of measurement system capabilities and design process of the payload is informed by computational fluid dynamics predictions, accuracy modeling based on experimental data, camera simulation, and lander design, all of which are combined in the Virtual Diagnostic Interface. The second SCALPSS mission, traveling to the moon on Firefly Aerospace’s Blue Ghost lander in early 2023, aims to build on the design and complexity of the first payload, which is launching in early 2022 on Intuitive Machines’ Nova-C lander. The Blue Ghost SCALPSS mission will include two additional cameras and a total of three different lens focal lengths which will target specific points in the vehicle’s descent to obtain quantitative and accurate 3D reconstruction of the lunar surface both prior to and after crater formation.

photogrammetry↗

Advancing xEMU Lunar Dust Mitigation Devices

Fine, electrically charged, glass like dust particles caused significant damage to the Apollo EMU9 during lunar EVAs, identified as one of the greatest challenges to future exploration. Passive Lunar Dust Mitigation Devices (LDMD) were developed, within a SBIR Phase II, to prohibit this dust from interrupting venting space suit component operation. A Computational Fluid-Dynamics and Discrete Element Method (CFD-DEM) Simulation Tool was developed at the University of Colorado, Boulder to predict venting gas flow ability to self-clean adhered dust particles from LDMD surfaces. Lunar dust properties (i.e., adhesion and cohesion strengths) required to complete Simulation Tool analysis are relatively unknown due to considerable differences between the Earth and the Moon (i.e., gravity, humidity) and due to an absence of dust particles in their native state. Analytically determining gas flow velocity, density and direction within the fluid Boundary Layer, microns from LDMD surfaces presented a second challenge. Dusty Plasma Chamber testing is being performed at Auburn University, Auburn to observe electrostatically charged dust behavior as it adheres to LDMD prototypes and specific geometry and is then blown away by metered gas flow. Test articles were developed to offer insight into the impact of different flow geometries, surface roughness and dust removal within the gas Boundary-Layer. Observed dust behavior is currently being developed to support the CFD-DEM analysis. Many Simulation Tool analytical cases have been processed to support the intention of completing sensitivity studies to assess how different dust adherence values and Boundary Layer fluid properties impact LDMD self-cleaning effectivity.

Thomas J Stapleton↗

Numerical Simulation of a High-Lift Configuration Embedded with High Momentum Fluidic Actuators

Numerical simulations have been performed for a vertical tail configuration with deflected rudder. The suction surface of the main element of this configuration, just upstream of the hinge line, is embedded with an array of 32 fluidic actuators that produce oscillating sweeping jets. Such oscillating jets have been found to be very effective for flow control applications in the past. In the current paper, a high-fidelity computational fluid dynamics (CFD) code known as the PowerFLOW R code is used to simulate the entire flow field associated with this configuration, including the flow inside the actuators. A fully compressible version of the PowerFLOW R code valid for high speed flows is used for the present simulations to accurately represent the transonic flow regimes encountered in the flow field due to the actuators operating at higher mass flow (momentum) rates required to mitigate reverse flow regions on a highly-deflected rudder surface. The computed results for the surface pressure and integrated forces compare favorably with measured data. In addition, numerical solutions predict the correct trends in forces with active flow control compared to the no control case. The effect of varying the rudder deflection angle on integrated forces and surface pressures is also presented.

Vatsa, Veer N.↗

Experience with Aero- and Fluid-Dynamic Testing for Engineering and CFD Validation

Ever since computations have been used to simulate aerodynamics the need to ensure that the computations adequately represent real life has followed. Many experiments have been performed specifically for validation and as computational methods have improved, so have the validation experiments. Validation is also a moving target because computational methods improve requiring validation for the new aspect of flow physics that the computations aim to capture. Concurrently, new measurement techniques are being developed that can help capture more detailed flow features pressure sensitive paint (PSP) and particle image velocimetry (PIV) come to mind. This paper will present various wind-tunnel tests the author has been involved with and how they were used for validation of various kinds of CFD. A particular focus is the application of advanced measurement techniques to flow fields (and geometries) that had proven to be difficult to predict computationally. Many of these difficult flow problems arose from engineering and development problems that needed to be solved for a particular vehicle or research program. In some cases the experiments required to solve the engineering problems were refined to provide valuable CFD validation data in addition to the primary engineering data. All of these experiments have provided physical insight and validation data for a wide range of aerodynamic and acoustic phenomena for vehicles ranging from tractor-trailers to crewed spacecraft.

experiment↗

Turbulence Simulations of Transonic Flows over NACA-0012 and OAT15A Airfoils

Three different simulation approaches, namely unsteady Reynolds-averaged Navier-Stokes (URANS), delayed detached-eddy simulation (DDES), and wall-modeled large-eddy simulation (WMLES) are employed to simulate transonic flow over an NACA-0012, OAT15A airfoils at different angles of attack covering pre- and post-buffet-onset regimes. The freestream Mach number is 0.75, and the Reynolds number based on the chord length is 10 and 3million, respectively. The NASA FUN3D solver is used for the simulations, which is an unstructured, compressible flow solver. The URANS simulations are performed using the Spalart-Allmaras (SA) model with the compressibility correction, the DDES predictions are based on the SA model, and the WMLES are performed using an equilibrium wall-model. The unsteady RANS simulations, only with the compressibility correction, predict the pre- and post- buffet characteristics, which compare well with the experimental results. DDES results predicted a lower buffet onset angle compared to experiment. The predicted shock locations are upstream of the locations predicted by URANS. Using a fine grid in the spanwise direction, WMLES predictions show buffeting consistent with the experiment.

Computational Fluid Dynamics↗

Turbulence Simulations of Transonic Flows over NACA-0012 and OAT15A Airfoils

Three different simulation approaches, namely unsteady Reynolds-averaged Navier-Stokes (URANS), delayed detached-eddy simulation (DDES), and wall-modeled large-eddy simulation (WMLES) are employed to simulate transonic flow over an NACA-0012, OAT15A airfoils at different angles of attack covering pre- and post-buffet-onset regimes. The freestream Mach number is 0.75, and the Reynolds number based on the chord length is 10 and 3million, respectively. The NASA FUN3D solver is used for the simulations, which is an unstructured, compressible flow solver. The URANS simulations are performed using the Spalart-Allmaras (SA) model with the compressibility correction, the DDES predictions are based on the SA model, and the WMLES are performed using an equilibrium wall-model. The unsteady RANS simulations, only with the compressibility correction, predict the pre- and post- buffet characteristics, which compare well with the experimental results. DDES results predicted a lower buffet onset angle compared to experiment. The predicted shock locations are upstream of the locations predicted by URANS. Using a fine grid in the spanwise direction, WMLES predictions show buffeting consistent with the experiment.

Computational Fluid Dynamics↗

Simulation of Atmospheric-Entry Capsules in the Subsonic Regime

The accuracy of Computational Fluid Dynamics predictions of subsonic capsule aerodynamics is examined by comparison against recent NASA wind-tunnel data at high-Reynolds-number flight conditions. Several aspects of numerical and physical modeling are considered, including inviscid numerical scheme, mesh adaptation, rough-wall modeling, rotation and curvature corrections for eddy-viscosity models, and Detached-Eddy Simulations of the unsteady wake. All of these are considered in isolation against relevant data where possible. The results indicate that an improved predictive capability is developed by considering physics-based approaches and validating the results against flight-relevant experimental data.

Simulation↗

Wind-tunnel preflight test program for aeroassist flight experiment

A proposed Aeroassist Flight Experiment (AFE) will simulate return from geosynchronous orbit to low-earth orbit to provide fundamental information about the hypersonic flow conditions surrounding aeroassisted orbital transfer vehicles (AOTV's). Future AOTV design must rely heavily on computational fluid dynamic (CFD) computer codes now under development, since the entry conditions are beyond the simulation capabilities of ground-based test facilities. However, existing codes and best available ground-based test facilities must provide for design of the AFE. The AFE ground-based test program provides data for calibration of existing and developing pertinent CFD codes. Confidence in code prediction capability increases with ability to predict forces and moments, pressure and heat-transfer distributions, shock shapes, and surface streamline directions. The ground-based test program provides these data over a wide range of hypersonic test parameters. This paper describes the test program, models, facilities, and some representative results.

Wells, William L.↗

A Method for Flow Simulation About Complex Geometries Using Both Structured and Unstructured Grids

A computational fluid dynamics code which utilizes both structured and unstructured grids was developed. The objective of this study was to develop and demonstrate the ability of such a code to achieve solutions about complex geometries in two dimensions. An unstructured grid generator and flow solver were incorporated into the PARC2D structured flow solver. This new unstructured grid generator capability allows for easier generation and manipulation of complex grids. Several examples of the grid generation capabilities are provided. The coupling of different grid topologies and the manipulation of individual grids is shown. Also, grids for realistic geometries, a NACA 0012 airfoil and a wing/nacelle installation, were created. The flow over a NACA 0012 airfoil was used as a test case for the flow solver. Eight separate cases were run. They were both the inviscid and viscous solutions for two freestream Mach numbers and airfoil angle of attacks of 0 to 3.86 degrees. The Mach numbers chosen were for a subsonic case, Mach 0.6, and a case where supersonic regions and a shock wave exists, Mach 0.8. These test case conditions were selected to match experimentally obtained data for code comparison. The results show that the code accurately predicts the flow field for all cases.

Debonis, James R.↗

Wall-modeled LES of Turbulent Flow Over a Two-dimensional Gaussian Bump

We perform wall-modeled large eddy simulations (WMLES) of turbulent flow over a nominally two-dimensional Gaussian-shaped bump geometry to assess its performance in the accelerating and separation regions of the flow. The flow conditions are based on high-fidelity numerical simulation of Uzun & Malik (AIAA Journal 2022). The oncoming flow Mach number is 0.2, with the bump length-based Reynolds number (ReL) of 2 million. In our previous study, while WMLES with the constant coefficient Vreman subgrid scale model and equilibrium wall model performed satisfactorily at lower ReL = 1 million, it failed to predict flow separation at the higher Re, contrary to the experimental observations. We investigate the sensitivity of WMLES to different subgrid scale models, wall models, and grid resolution and topology on flow separation by comparing with available data.

Computational Fluid Dynamics↗

High-Density Plasma Reactors: Simulations for Design

The development of improved and more efficient plasma reactors is a costly process for the semiconductor industry. Until five years ago, the Industry made most of its advancements through a trial and error approach. More recently, the role of computational modeling in the design process has increased. Both conventional computational fluid dynamics (CFD) techniques like Navier-Stokes solvers as well as particle simulation methods are used to model plasma reactor flowfields. However, since high-density plasma reactors generally operate at low gas pressures on the order of 1 to 10 mTorr, a particle simulation may be necessary because of the failure of CFD techniques to model rarefaction effects. The direct simulation Monte Carlo method is the most widely accepted and employed particle simulation tool and has previously been used to investigate plasma reactor flowfields. A plasma DSMC code is currently under development at NASA Ames Research Center with its foundation as the object-oriented parallel Cornell DSMC code, MONACO. The present investigation is a follow up of a neutral flow investigation of the effects of process parameters as well as reactor design on etch rate and etch rate uniformity. The previous work concentrated on silicon etch of a chlorine flow in a configuration typical of electron cyclotron resonance (ECR) or helical resonator type reactors. The effects of the plasma on the dissociation chemistry were modeled by making assumptions about the electron temperature and number density. The electrons or ions themselves were not simulated.The present work extends these results by simulating the charged species.The electromagnetic fields are calculated such that power deposition is modeled self-consistently. Electron impact reactions are modeled along with mechanisms for charge exchange. An bipolar diffusion assumption is made whereby electrons remain tied to the ions. However, the velocities of tile electrons are allowed to be modified during collisions and are not confined to a Maxwellian distribution. The interaction between the neutral flow and plasma is examined, and results for etch rate uniformity from the previous research and the present plasma simulations are compared.

Hash, David B.↗

A static data flow simulation study at Ames Research Center

Demands in computational power, particularly in the area of computational fluid dynamics (CFD), led NASA Ames Research Center to study advanced computer architectures. One architecture being studied is the static data flow architecture based on research done by Jack B. Dennis at MIT. To improve understanding of this architecture, a static data flow simulator, written in Pascal, has been implemented for use on a Cray X-MP/48. A matrix multiply and a two-dimensional fast Fourier transform (FFT), two algorithms used in CFD work at Ames, have been run on the simulator. Execution times can vary by a factor of more than 2 depending on the partitioning method used to assign instructions to processing elements. Service time for matching tokens has proved to be a major bottleneck. Loop control and array address calculation overhead can double the execution time. The best sustained MFLOPS rates were less than 50% of the maximum capability of the machine.

Barszcz, Eric↗

Adjustable Autonomy Testbed

The Adjustable Autonomy Testbed (AAT) is a simulation-based testbed located in the Intelligent Systems Laboratory in the Automation, Robotics and Simulation Division at NASA Johnson Space Center. The purpose of the testbed is to support evaluation and validation of prototypes of adjustable autonomous agent software for control and fault management for complex systems. The AA T project has developed prototype adjustable autonomous agent software and human interfaces for cooperative fault management. This software builds on current autonomous agent technology by altering the architecture, components and interfaces for effective teamwork between autonomous systems and human experts. Autonomous agents include a planner, flexible executive, low level control and deductive model-based fault isolation. Adjustable autonomy is intended to increase the flexibility and effectiveness of fault management with an autonomous system. The test domain for this work is control of advanced life support systems for habitats for planetary exploration. The CONFIG hybrid discrete event simulation environment provides flexible and dynamically reconfigurable models of the behavior of components and fluids in the life support systems. Both discrete event and continuous (discrete time) simulation are supported, and flows and pressures are computed globally. This provides fast dynamic simulations of interacting hardware systems in closed loops that can be reconfigured during operations scenarios, producing complex cascading effects of operations and failures. Current object-oriented model libraries support modeling of fluid systems, and models have been developed of physico-chemical and biological subsystems for processing advanced life support gases. In FY01, water recovery system models will be developed.

Malin, Jane T.↗

A Step Towards CO2-Neutral Aviation

An approximation method for evaluation of the caloric equations used in combustion chemistry simulations is described. The method is applied to generate the equations of specific heat, static enthalpy, and Gibb's free energy for fuel mixtures of interest to gas turbine engine manufacturers. Liquid-phase fuel properties are also derived. The fuels include JP-8, synthetic fuel, and two fuel blends consisting of a mixture of JP-8 and synthetic fuel. The complete set of fuel property equations for both phases are implemented into a computational fluid dynamics (CFD) flow solver database, and multi-phase, reacting flow simulations of a well-tested liquid-fueled combustor are performed. The simulations are a first step in understanding combustion system performance and operational issues when using alternate fuels, at practical engine operating conditions.

Brankovic, Andreja↗

CFD Simulations in Support of Shuttle Orbiter Contingency Abort Aerodynamic Database Enhancement

Modern Computational Fluid Dynamics (CFD) techniques were used to compute aerodynamic forces and moments of the Space Shuttle Orbiter in specific portions of contingency abort trajectory space. The trajectory space covers a Mach number range of 3.5-15, an angle-of-attack range of 20deg-60deg, an altitude range of 100-190 kft, and several different settings of the control surfaces (elevons, body flap, and speed brake). Presented here are details of the methodology and comparisons of computed aerodynamic coefficients against the values in the current Orbiter Operational Aerodynamic Data Book (OADB). While approximately 40 cases have been computed, only a sampling of the results is provided here. The computed results, in general, are in good agreement with the OADB data (i.e., within the uncertainty bands) for almost all the cases. However, in a limited number of high angle-of-attack cases (at Mach 15), there are significant differences between the computed results, especially the vehicle pitching moment, and the OADB data. A preliminary analysis of the data from the CFD simulations at Mach 15 shows that these differences can be attributed to real-gas/Mach number effects. The aerodynamic coefficients and detailed surface pressure distributions of the present simulations are being used by the Shuttle Program in the evaluation of the capabilities of the Orbiter in contingency abort scenarios.

Papadopoulos, Periklis E.↗

Evaluation of Wall-Modeled Les for Flow Over A Multi-Element Airfoil

We perform wall-modeled large-eddy simulations (WMLES) of turbulent flow over a 30P30N multielement airfoil at various angles of attack to assess its ability to accurately predict lift using an unstructured node-based finite-volume methodology. The flow conditions are based on the experiments reported in Klausmeyer & Lin [NASA /TM-112858]. The oncoming flow Mach number is 0.2, and the stowed chord-based Reynolds number (𝑅𝑒𝑐) is 9 million. The angle of attack (𝛼) is varied between 8◦ and 23◦ with the expected stall angle of around 21◦. The simulated geometry is periodic in the span, thus any three-dimensional effects present in the experiment are ignored. We first analyze the predictions with baseline isotropic and anisotropic grids containing over 16 grid points per boundary-layer thickness. The near-wall streamwise/nominal spacing ratio is 0.4 at the leading/trailing edges of the slat, main element and flap, and 1 elsewhere, and the first wall-normal/nominal spacing ratio is 0.375 over most of the airfoil for the anisotropic grid, while these are unity for the isotropic grid. It is found that the predicted lift and wall pressures, near-wall velocity profiles, and the location of onset of resolved turbulence are sensitive to the near-wall grid anisotropy. Further sensitivities were explored at lower angles of attack by independently varying the streamwise and wall-normal grid anisotropies. The effects of both off-wall and near-wall grid refinement were also quantified. For the grid resolutions used here, our results suggest that the first wall-normal spacing, which coincides with the WMLES exchange location, has a dominant effect on the predictions for this flow configuration.

Computational Fluid Dynamics↗