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 901 records · Page 50

Assessment of Near-Field Sonic Boom Simulation Tools

A recent study for the Supersonics Project, within the National Aeronautics and Space Administration, has been conducted to assess current in-house capabilities for the prediction of near-field sonic boom. Such capabilities are required to simulate the highly nonlinear flow near an aircraft, wherein a sonic-boom signature is generated. There are many available computational fluid dynamics codes that could be used to provide the near-field flow for a sonic boom calculation. However, such codes have typically been developed for applications involving aerodynamic configuration, for which an efficiently generated computational mesh is usually not optimum for a sonic boom prediction. Preliminary guidelines are suggested to characterize a state-of-the-art sonic boom prediction methodology. The available simulation tools that are best suited to incorporate into that methodology are identified; preliminary test cases are presented in support of the selection. During this phase of process definition and tool selection, parallel research was conducted in an attempt to establish criteria that link the properties of a computational mesh to the accuracy of a sonic boom prediction. Such properties include sufficient grid density near shocks and within the zone of influence, which are achieved by adaptation and mesh refinement strategies. Prediction accuracy is validated by comparison with wind tunnel data.

Casper, J. H.↗

Implicit Thermochemical Nonequilibrium Flow Simulations on Unstructured Grids using GPUs

Thermochemical nonequilibrium flow simulation capabilities have been previously implemented, verified, and validated for central processing unit (CPU) systems in NASA’s unstructured-grid computational fluid dynamics solver FUN3D. Many exascale-class high-performance computing systems will rely on graphics processing unit (GPU) architectures for high throughput and energy efficiency; thus, CPU-based scientific computing software unable to effectively utilize these systems must be updated. In this work, we present a CUDA C++ implementation of FUN3D’s thermochemical nonequilibrium flow simulation capabilities targeting NVIDIA Tesla GPUs. An overview of the porting and optimization strategy is described and performance comparisons with other recent architectures are presented. Scaling to thousands of GPUs is demonstrated, yielding computational performance equivalent to that of several million CPU cores. The implementation enables efficient, high-fidelity, scale-resolving simulations of thermochemical nonequilibrium flows for many applications including atmospheric entry, hypersonics, and combustion.

GPU↗

Toward Adaptive Mixed-Element Unstructured Grids for Simulations of Viscous Flows

Accurate prediction of surface skin friction and heat transfer for viscous flow applications using computational fluid dynamics tools necessitates thin wall-normal elements in the boundary layer. A tetrahedral adaptation process is extended to incorporate the advancing front local reconnection procedure. Thin boundary-layer grids are adapted from the solution, which governs the boundary-layer grid parameters such as first cell height off the wall. The process enables grid-converged steady-state results from a CAD file. Results are presented across the speed range from low subsonic flow over an airfoil to high-speed flows over blunt bodies and a rocket.

CFD↗

Toward Adaptive Mixed-Element Unstructured Grids for Simulations of Viscous Flows

Accurate prediction of surface skin friction and heat transfer for viscous flow applications using computational fluid dynamics tools necessitates thin wall-normal elements in the boundary layer. A tetrahedral adaptation process is extended to incorporate the advancing front local reconnection procedure. Thin boundary-layer grids are adapted from the solution, which governs the boundary-layer grid parameters such as first cell height off the wall. The process enables grid-converged steady-state results from a CAD file. Results are presented across the speed range from low subsonic flow over an airfoil to high-speed flows over blunt bodies and a rocket.

chemically reacting flow↗

Emission Spectral Measurements in the Plenum of an Arc Jet Wind Tunnel

Arc jet wind tunnel facilities are used to evaluate thermal protection system materials for re-entry vehicles. The high speed, high temperature flowfield generated by the arc jet can simulate the extreme aerodynamic heating environment experienced during re-entry so that the survivability of heat shield materials and performance of various designs options can be tested. Although the re-entry heating environment can be approximated in the arc jet facility, the flowfield only partially simulates the actual re-entry flight conditions. Reynolds numbers are not matched so that surface shear stress distributions and mass transfer rates due to ablation or other mechanisms are not modeled correctly. Unlike flight conditions the arc freestream air is in non-equilibrium because of the rapid expansion that occurs in the supersonic nozzle. To properly study the actual re-entry flow environment, computational fluid dynamics, computational chemistry and radiation models must be used. Arc jet tunnel tests serve to validate these models. To perform accurate simulations inlet and boundary-conditions are needed, which come from measurements of the flowfield. The present study is concerned with measurements in the plenum region of an arc heater. In the past, conditions in the arc heater flowfield have been predicted using simulations since conventional measurement techniques could not be used in the harsh extremely high temperature environment. The present study is part of a recent push to utilize optical techniques to help better characterize the arc jet flowfields. Emission measurements have been made in the shock layer and the constrictor section of the arc heater to determine temperatures and species number densities. LIF measurements have been made in the free stream to determine temperature and velocity.

Donohue, Jim↗

Parallel 3D Mortar Element Method for Adaptive Nonconforming Meshes

High order methods are frequently used in computational simulation for their high accuracy. An efficient way to avoid unnecessary computation in smooth regions of the solution is to use adaptive meshes which employ fine grids only in areas where they are needed. Nonconforming spectral elements allow the grid to be flexibly adjusted to satisfy the computational accuracy requirements. The method is suitable for computational simulations of unsteady problems with very disparate length scales or unsteady moving features, such as heat transfer, fluid dynamics or flame combustion. In this work, we select the Mark Element Method (MEM) to handle the non-conforming interfaces between elements. A new technique is introduced to efficiently implement MEM in 3-D nonconforming meshes. By introducing an "intermediate mortar", the proposed method decomposes the projection between 3-D elements and mortars into two steps. In each step, projection matrices derived in 2-D are used. The two-step method avoids explicitly forming/deriving large projection matrices for 3-D meshes, and also helps to simplify the implementation. This new technique can be used for both h- and p-type adaptation. This method is applied to an unsteady 3-D moving heat source problem. With our new MEM implementation, mesh adaptation is able to efficiently refine the grid near the heat source and coarsen the grid once the heat source passes. The savings in computational work resulting from the dynamic mesh adaptation is demonstrated by the reduction of the the number of elements used and CPU time spent. MEM and mesh adaptation, respectively, bring irregularity and dynamics to the computer memory access pattern. Hence, they provide a good way to gauge the performance of computer systems when running scientific applications whose memory access patterns are irregular and unpredictable. We select a 3-D moving heat source problem as the Unstructured Adaptive (UA) grid benchmark, a new component of the NAS Parallel Benchmarks (NPB). In this paper, we present some interesting performance results of ow OpenMP parallel implementation on different architectures such as the SGI Origin2000, SGI Altix, and Cray MTA-2.

Feng, Huiyu↗

Summary of the First AIAA Stability and Control Prediction Workshop

Results from the First AIAA Stability and Control Prediction Workshop are summarized in this paper. The workshop series was developed in support of three primary objectives: (1) to establish best practices for the prediction of stability and control derivatives using industry-standard Computational Fluid Dynamics (CFD) solvers, (2) to provide an impartial forum for evaluating the effectiveness of Reynolds-averaged-Navier-Stokes- and Detached-Eddy-Simulation-based modeling techniques, and (3) to identify areas in need of additional research and development. To address these objectives, the inaugural workshop focused on generating computational aerodynamic predictions for the ONERA version of the NASA/Boeing Common Research Model. The configuration includes the wing, body, horizontal tail, and a vertical tail designed by ONERA. While longitudinal wind tunnel test data for the model had been previously documented, unpublished lateral test data at small sideslip angles provided a unique opportunity for participants to generate blind computational predictions for wind tunnel data comparisons. Additional test cases included assessments of the Mach number effect on static lateral/directional stability and the impact of the wind tunnel sting on longitudinal stability. Participants were invited to generate solutions using two workshop-provided series of structured overset and unstructured grids, in addition to any participant custom grids of interest. Total and component-level breakdowns for the force and moment coefficients for each test case are presented, as well as sectional pressure distribution data on the wing and tail components, to assess the agreement between several different Reynolds-averaged Navier-Stokes CFD solvers.

CFD↗

Applications of CFD and visualization techniques

In this paper, three applications are presented to illustrate current techniques for flow calculation and visualization. The first two applications use a commercial computational fluid dynamics (CFD) code, FLUENT, performed on a Cray Y-MP. The results are animated with the aid of data visualization software, apE. The third application simulates a particulate deposition pattern using techniques inspired by developments in nonlinear dynamical systems. These computations were performed on personal computers.

Saunders, James H.↗

Radiative Heating of Large Meteoroids During Atmospheric Entry

A high-fidelity approach for simulating the aerothermodynamic environments of meteor entries was developed, which allows the commonly assumed heat transfer coefficient of 0.1 to be assessed. This model uses chemically reacting computational fluid dynamics (CFD), coupled with radiation transport and surface ablation. Coupled radiation accounts for the impact of radiation on the flowfield energy equations, while coupled ablation explicitly models the injection of ablation products within the flowfield and radiation simulations. For a meteoroid with a velocity of 20 km/s, coupled radiation is shown to reduce the stagnation point radiative heating by over 60%. The impact of coupled ablation (with coupled radiation) is shown to provide at least a 70% reduction in the radiative heating relative to cases with only coupled radiation. This large reduction is partially the result of the low ionization energies of meteoric ablation products relative to air species. The low ionization energies of ablation products, such as Mg and Ca, provide strong photoionization and atomic line absorption in regions of the spectrum that air species do not. MgO and CaO are also shown to provide significant absorption. Turbulence is shown to impact the distribution of ablation products through the shock-layer, which results in up to a 100% increase in the radiative heating downstream of the stagnation point. To create a database of heat transfer coefficients, the developed model was applied to a range of cases. This database considered velocities ranging from 14 to 20 km/s, altitudes ranging from 20 to 50 km, and nose radii ranging from 1 to 100 m. The heat transfer coefficients from these simulations are below 0.045 for the range of cases, for both laminar and turbulent, which is significantly lower than the canonical value of 0:1. When the new heat transfer model is applied to a Tunguska-like 15 Mt entry, the effect of the new model is to lower the height of burst by up to 2 km, depending on assumed entry angle. This, in turn, results in a significantly larger ground damage footprint than when the canonical heating assumption is used.

Christopher O Johnston↗

Prediction of vortex breakdown on a delta wing

Recent studies of leading-edge vortex flows with computational fluid dynamics codes using Euler or Navier-Stokes formulations have shown fair agreement with experimental data. These studies have concentrated on simulating the flowfields associated with a sharp-edged flat plate 70 deg delta wing at angles of attack where vortex breakdown or burst is observed over the wing. There are, however, a number of discrepancies between the experimental data and the computed flowfields. The location of vortex breakdown in the computational solutions is seen to differ from the experimental data and to vary with changes in the computational grid and freestream Mach number. There also remain issues as to the validity of steady-state computations for cases which contain regions of unsteady flow, such as in the post-breakdown regions. As a partial response to these questions, a number of laminar Navier-Stokes solutions were examined for the 70 deg delta wing. The computed solutions are compared with an experimental database obtained at low subsonic speeds. The convergence of forces, moments and vortex breakdown locations are also analyzed to determine if the computed flowfields actually reach steady-state conditions.

Agrawal, S.↗

Scale Resolving Simulations of Contra Rotating Open Rotor Noise Prediction

The development of breakthrough open-rotor propulsion systems that meet performance, fuel efficiency, and noise requirements for commercial viability relies heavily on the availability of reliable and computationally efficient aerodynamic performance and noise prediction tools.This investigation showcases the proficiency of the structured curvilinear overset grid Computational Fluid Dynamics solver implemented in the Launch, Ascent, and Vehicle Aerodynamics framework for predicting the aerodynamic and aeroacoustic characteristics of a Contra-RotatingOpen Rotor propulsion system. Additionally, the study includes comprehensive sensitivity analyses employing a Hybrid Reynolds-Averaged Navier-Stokes/Large Eddy Simulation - based scale-resolving turbulence closure treatment. The computed results are validated through comparison with existing experimental measurements.

ARMD↗

Performance of Coupled Physics Solvers for Multidisciplinary Hypersonic Flow Simulations on Several Classes of Computer Architectures

The application of hypersonic flow simulation tools to realistic flight scenarios will require the coupling of multiple physical effects to the baseline fluid dynamics. Such multiphysics effects can include the aerooelastic response of the airframe or engine components, dynamic transport of atmospheric particles, the deformation of solid-fluid interfaces that can ablate, pyrolyze, or erode, as well as a host of other processes, all of which are governed by unique sets of physical equations and models. Coupling multiple (and potentially disparate) physics solvers to a robust compressible flow solver poses additional challenges related to the stability, performance and scalability of the combined solver. The choices made during the software design process can therefore lead to a variation in simulation efficiency across different computer architectures. In this paper, we will consider two representative multiphysics hypersonic flow scenarios: the interaction of solid particulates with the flow field created by a hypersonic lifting body and the aerooelastic deformation of a model airframe under high-Mach-number flow conditions. For these simulations we explore the behavior of several hypersonic simulation tools, including Kestrel, FUN3D, US3D, and JENRE multiphysics framework, on several high performance computing systems containing various CPU and GPU architectures.

architecture↗

Performance Evaluation of the Gravity Probe B Design

This report documents the simulation of the Lockheed Martin designed Gravity Probe B (GPB) spacecraft developed tool by bd Systems Inc using the TREETOPS simulation. This study quantifies the effects of flexibility and liquid helium slosh on GPB spacecraft control performance. The TREETOPS simulation tool permits the simulation of flexible structures given that a flexible body model of the structure is available. For purposes of this study, a flexible model of the GPB spacecraft was obtained from Lockheed Martin. To model the liquid helium slosh effects, computational fluid dynamics (CFD) results' were obtained, and used to develop a dynamic model of the slosh effects. The flexible body and slosh effects were incorporated separately into the TREETOPS simulation, which places the vehicle in a 650 km circular polar orbit and subjects the spacecraft to realistic environmental disturbances and sensor error quantities. In all of the analysis conducted in this study the spacecraft is pointed at an inertially fixed guide star (GS) and is rotating at a constant rate about this line of sight.

Francis, Ronnie↗

Unsteady Analysis of Inlet-Compressor Acoustic Interactions Using Coupled 3-D and 1-D CFD Codes

It is well known that the dynamic response of a mixed compression supersonic inlet is very sensitive to the boundary condition imposed at the subsonic exit (engine face) of the inlet. In previous work, a 3-D computational fluid dynamics (CFD) inlet code (NPARC) was coupled at the engine face to a 3-D turbomachinery code (ADPAC) simulating an isolated rotor and the coupled simulation used to study the unsteady response of the inlet. The main problem with this approach is that the high fidelity turbomachinery simulation becomes prohibitively expensive as more stages are included in the simulation. In this paper, an alternative approach is explored, wherein the inlet code is coupled to a lesser fidelity 1-D transient compressor code (DYNTECC) which simulates the whole compressor. The specific application chosen for this evaluation is the collapsing bump experiment performed at the University of Cincinnati, wherein reflections of a large-amplitude acoustic pulse from a compressor were measured. The metrics for comparison are the pulse strength (time integral of the pulse amplitude) and wave form (shape). When the compressor is modeled by stage characteristics the computed strength is about ten percent greater than that for the experiment, but the wave shapes are in poor agreement. An alternate approach that uses a fixed rise in duct total pressure and temperature (so-called 'lossy' duct) to simulate a compressor gives good pulse shapes but the strength is about 30 percent low.

Suresh, A.↗

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↗

CYBER 200 Applications Seminar

Applications suited for the CYBER 200 digital computer are discussed. Various areas of application including meteorology, algorithms, fluid dynamics, monte carlo methods, petroleum, electronic circuit simulation, biochemistry, lattice gauge theory, economics and ray tracing are discussed.

Gary, J. P.↗

Transonic Symposium: Theory, Application, and Experiment, volume 1, part 2

In order to assess the state of the art in transonic flow disciplines and to glimpse at future directions, NASA-Langley held a Transonic Symposium. Emphasis was placed on steady, three dimensional external, transonic flow and its simulation, both numerically and experimentally. The symposium included technical sessions on wind tunnel and flight experiments; computational fluid dynamic applications; inviscid methods and grid generation; viscous methods and boundary layer stability; and wind tunnel techniques and wall interference. This, being volume 1, is unclassified.

Foughner, Jerome T., Jr.↗

Response of gravity level fluctuations on the Gravity Probe-B spacecraft propellant system

Time-dependent computations have been performed to investigate the dynamical behaviors of fluid under a microgravity environment. A computer algorithm is introduced which can be used to simulate the fluid behavior in that environment, in particular the excitation of sloshing waves due to different gravity environments and rotation speeds. A suggestion on the proper handling and managing of cryogenic fluid propellant to be used in the Gravity Probe-B spacecraft propulsion is made.

Hung, R. J.↗