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 757 records · Page 42

Fluorescence Imaging of Underexpanded Jets and Comparison with CFD

An experimental study of underexpanded and highly underexpanded axisymmetric nitrogen free jets seeded with 0.5% nitric oxide (NO) and issuing from a sonic orifice was conducted at NASA Langley Research Center. Reynolds numbers based on nozzle exit conditions ranged from 770 to 35,700, and nozzle exit-to-ambient jet pressure ratios ranged from 2 to 35. These flows were non-intrusively visualized with a spatial resolution of approximately 0.14 mm x 0.14 mm x 1 mm thick and a temporal resolution of 1 s using planar laser-induced fluorescence (PLIF) of NO, with the laser tuned to the strongly-fluorescing UV absorption bands of the Q1 band head near 226.256 nm. Three laminar cases were selected for comparison with computational fluid dynamics (CFD). The cases were run using GASP (General Aerodynamic Simulation Program) Version 4. Comparisons of the fundamental wavelength of the jet flow showed good agreement between CFD and experiment for all three test cases, while comparisons of Mach disk location and Mach disk diameter showed good agreement at lower jet pressure ratios, with a tendency to slightly underpredict these parameters with increasing jet pressure ratio.

Wilkes, Jennifer A.↗

Twist Model Development and Results from the Active Aeroelastic Wing F/A-18 Aircraft

Understanding the wing twist of the active aeroelastic wing (AAW) F/A-18 aircraft is a fundamental research objective for the program and offers numerous benefits. In order to clearly understand the wing flexibility characteristics, a model was created to predict real-time wing twist. A reliable twist model allows the prediction of twist for flight simulation, provides insight into aircraft performance uncertainties, and assists with computational fluid dynamic and aeroelastic issues. The left wing of the aircraft was heavily instrumented during the first phase of the active aeroelastic wing program allowing deflection data collection. Traditional data processing steps were taken to reduce flight data, and twist predictions were made using linear regression techniques. The model predictions determined a consistent linear relationship between the measured twist and aircraft parameters, such as surface positions and aircraft state variables. Error in the original model was reduced in some cases by using a dynamic pressure-based assumption. This technique produced excellent predictions for flight between the standard test points and accounted for nonlinearities in the data. This report discusses data processing techniques and twist prediction validation, and provides illustrative and quantitative results.

Lizotte, Andrew M.↗

A Perspective on Computational Aerothermodynamics at NASA

The evolving role of computational aerothermodynamics (CA) within NASA over the past 20 years is reviewed. The presentation highlights contributions to understanding the Space Shuttle pitching moment anomaly observed in the first shuttle flight, prediction of a static instability for Mars Pathfinder, and the use of CA for damage assessment in post-Columbia mission support. In the view forward, several current challenges in computational fluid dynamics and aerothermodynamics for hypersonic vehicle applications are discussed. Example simulations are presented to illustrate capabilities and limitations. Opportunities to advance the state-of-art in algorithms, grid generation and adaptation, and code validation are identified.

Gnoffo, Peter A.↗

Development of an Effective Finite-rate Oxidation Model for NuSil-coated Charred Carbon Preform Ablators

Recently, a detailed effective finite-rate surface chemistry model was developed for the oxidation of FiberForm [1] using the molecular beam experimental data of Poovathingal et al., [2]. FiberForm is the major building block of the thermal protection system (TPS) material Phenolic Impregnated Carbon Ablator (PICA), commonly used by NASA. The surface chemistry model consists of detailed surface reaction mechanisms such as adsorption, desorption, and several types of Langmuir-Hinshelwood (LH) reactions to characterize the oxygen-carbon interactions at the surface. This model provides excellent agreement with the experimental data for oxidation product compositions and corresponding translational energy distributions. Further, an effective oxidation model was constructed that captures the equivalent interaction of oxygen inside the microstructure (including multiple surface collisions) via a modified reactivity for a smooth wall boundary condition. This enables the use of this model directly into Computational Fluid Dynamics (CFD) codes and Material Response (MR) codes to accurately simulate the gas–surface interactions within FiberForm without using the detailed micro-structure. However, this effective model is valid only for the virgin FiberForm. As the TPS material undergoes ablation, the preform carbon burns and turns into char. In addition, PICA is also coated with a protective silicon coating called NuSil for the purpose of mitigating the spread of phenolic dust, and limit contamination during clean room operations. Fig. 1 shows a X-ray microtomography image of a charred NuSil-coated PICA. The NuSil layer introduces species containing silicon into the product mixture. Further, the reactivity of the carbon within the char layer is different from the virgin FiberForm as shown in Fig. 2. Recently a new set of molecular beam experiments were performed on this NuSil-coated TPS material [3] with the same type of oxygen beam used in the previous experiments. Using the latest experimental data, the previously developed effective model will be extended to account for the charred carbon as well as the NuSil coating. Finally, this new effective model with three phases – preform carbon, char, and NuSil; will be compared and validated against the experimental product compositions.

K Swaminathan Gopalan↗

Fluid-Structure Interaction Simulations of the ASPIRE SR03 Supersonic Parachute Flight Test

Research into parachute performance continues to be a source of significant investment from the National Aeronautics and Space Administration to mitigate risks and to enable a variety of exploration missions, including landing on Mars as well as returning to Earth. The cost of flight tests to certify any changes to the current state-of-the-art parachute designs limits the development of next generation parachute systems. Fluid-structure interaction simulations could help accelerate this process once validated. The Launch, Ascent, and Vehicle Aerodynamics team is developing the capability to perform such fluid-structure interaction simulations by coupling a higher-order Cartesian immersed boundary computational fluid dynamics solver with adaptive mesh refinement to a finite element structural dynamics solver in space and time. We continue the effort to validate this tool with the Advanced Supersonic Parachute Inflation Research Experiments SR03 flight test featuring a strengthened parachute akin to the Mars 2020 mission that landed the Perseverance rover on Mars, and a higher freestream dynamic pressure prior to inflation. The effect of the flow conditions’ angle of attack and of the initial parachute shape are quantified. The impact of relaxing modeling assumptions with regards to radial stiffeners on the parachute canopy is also investigated. Results demonstrate improvements in agreement with the pull force recorded during the SR03 flight test as the initial conditions of the flow and parachute are brought closer to those experienced in flight, and further improved when the radial stiffener modeling assumptions are relaxed.

ESM↗

Fluid-Structure Interaction Simulations of the ASPIRE SR03 Supersonic Parachute Flight Test

Research into parachute performance continues to be a source of significant investment from the National Aeronautics and Space Administration to mitigate risks and to enable a variety of exploration missions, including landing on Mars as well as returning to Earth. The cost of flight tests to certify any changes to the current state-of-the-art parachute designs limits the development of next generation parachute systems. Fluid-structure interaction simulations could help accelerate this process once validated. The Launch, Ascent, and Vehicle Aerodynamics team is developing the capability to perform such fluid-structure interaction simulations by coupling a higher-order Cartesian immersed boundary computational fluid dynamics solver with adaptive mesh refinement to a finite element structural dynamics solver in space and time. We continue the effort to validate this tool with the Advanced Supersonic Parachute Inflation Research Experiments SR03 flight test featuring a strengthened parachute akin to the Mars 2020 mission that landed the Perseverance rover on Mars, and a higher freestream dynamic pressure prior to inflation. The effect of the flow conditions’ angle of attack and of the initial parachute shape are quantified. The impact of relaxing modeling assumptions with regards to radial stiffeners on the parachute canopy is also investigated. Results demonstrate improvements in agreement with the pull force recorded during the SR03 flight test as the initial conditions of the flow and parachute are brought closer to those experienced in flight, and further improved when the radial stiffener modeling assumptions are relaxed.

ESM↗

Simplified jet-A kinetic mechanism for combustor application

Successful modeling of combustion and emissions in gas turbine engine combustors requires an adequate description of the reaction mechanism. For hydrocarbon oxidation, detailed mechanisms are only available for the simplest types of hydrocarbons such as methane, ethane, acetylene, and propane. These detailed mechanisms contain a large number of chemical species participating simultaneously in many elementary kinetic steps. Current computational fluid dynamic (CFD) models must include fuel vaporization, fuel-air mixing, chemical reactions, and complicated boundary geometries. To simulate these conditions a very sophisticated computer model is required, which requires large computer memory capacity and long run times. Therefore, gas turbine combustion modeling has frequently been simplified by using global reaction mechanisms, which can predict only the quantities of interest: heat release rates, flame temperature, and emissions. Jet fuels are wide-boiling-range hydrocarbons with ranges extending through those of gasoline and kerosene. These fuels are chemically complex, often containing more than 300 components. Jet fuel typically can be characterized as containing 70 vol pct paraffin compounds and 25 vol pct aromatic compounds. A five-step Jet-A fuel mechanism which involves pyrolysis and subsequent oxidation of paraffin and aromatic compounds is presented here. This mechanism is verified by comparing with Jet-A fuel ignition delay time experimental data, and species concentrations obtained from flametube experiments. This five-step mechanism appears to be better than the current one- and two-step mechanisms.

Lee, Chi-Ming↗

Conjugate Heat Transfer Analyses on the Manifold for Ramjet Fuel Injectors

Three-dimensional conjugate heat transfer analyses on the manifold located upstream of the ramjet fuel injector are performed using CFdesign, a finite-element computational fluid dynamics (CFD) software. The flow field of the hot fuel (JP-7) flowing through the manifold is simulated and the wall temperature of the manifold is computed. The three-dimensional numerical results of the fuel temperature are compared with those obtained using a one-dimensional analysis based on empirical equations, and they showed a good agreement. The numerical results revealed that it takes around 30 to 40 sec to reach the equilibrium where the fuel temperature has dropped about 3 F from the inlet to the exit of the manifold.

Wang, Xiao-Yen J.↗

Flight test and numerical simulation of transonic flow around YAV-8B Harrier II wing

A computational fluid dynamics (CFD) method is used to study the aerodynamics of the YAV-8B Harrier II wing in the transonic region. A numerical procedure is developed to compute the flow field around the complicated wing-pylon-fairing geometry. The surface definition of the wing and pylons were obtained from direct measurement using theodolite triangulation. A thin-layer Navier-Stokes code with the Chimera technique is used to compute flow solutions. The computed pressure distributions at several span stations are compared with flight test data and show good agreement. Computed results are correlated with flight test data that show the flow is severely separated in the vicinity of the wing-pylon junction. Analysis shows that shock waves are induced by pylon swaybrace fairings, that the flow separation is much stronger at the outboard pylon and that the separation is caused mainly by the crossflow passing the geometry of wing-pylon junction.

Gea, Lie-Mine↗

Development of the Glenn Heat-Transfer (Glenn-HT) Computer Code to Enable Time-Filtered Navier-Stokes (TFNS) Simulations and Application to Film Cooling on a Flat Plate Through Long Cooling Tubes

Computational fluid dynamics (CFD) analysis using Reynolds-averaged Navier-Stokes (RANS) formulation for turbomachinery-related flows has enabled improved engine component designs. RANS methodology has limitations that are related to its inability to accurately describe the spectrum of flow phenomena encountered in engines. Examples of flows that are difficult to compute accurately with RANS include phenomena such as laminar/turbulent transition, turbulent mixing due to mixing of streams, and separated flows. Large eddy simulation (LES) can improve accuracy but at a considerably higher cost. In recent years, hybrid schemes that take advantage of both unsteady RANS and LES have been proposed. This study investigated an alternative scheme, the time-filtered Navier-Stokes (TFNS) method applied to compressible flows. The method developed by Shih and Liu was implemented in the Glenn-Heat-Transfer (Glenn-HT) code and applied to film-cooling flows. In this report the method and its implementation is briefly described. The film effectiveness results obtained for film cooling from a row of 30deg holes with a pitch of 3.0 diameters emitting air at a nominal density ratio of unity and two blowing ratios of 0.5 and 1.0 are shown. Flow features under those conditions are also described.

Ameri, Ali↗

Development of the Glenn-HT Computer Code to Enable Time-Filtered Navier-Stokes (TFNS) Simulations and Application to Film Cooling on a Flat Plate Through Long Cooling Tubes

Computational fluid dynamics (CFD) analysis using Reynolds-averaged Navier-Stokes (RANS) formulation for turbomachinery-related flows has enabled improved engine component designs. RANS methodology has limitations which are related to its inability to accurately describe the spectrum of flow phenomena encountered in engines. Examples of flows that are difficult to compute accurately with RANS include phenomena such as laminarturbulent transition, turbulent mixing due to mixing of streams, and separated flows. Large eddy simulation (LES) can improve accuracy but at a considerably higher cost. In recent years, hybrid schemes which take advantage of both unsteady RANS and LES have been proposed. This study investigated an alternative scheme, the time-filtered Navier-Stokes (TFNS) method applied to compressible flows. The method developed by Shih and Liu was implemented in the Glenn-HT code and applied to film cooling flows. In this report the method and its implementation is briefly described. The film effectiveness results obtained for film cooling from a row of 30 holes with a pitch of 3.0 diameters emitting air at a nominal density ratio of unity and four blowing ratios of 0.5, 1.0, 1.5 and 2.0 are shown. Flow features under those conditions are also described.

infrared↗

Development of the Glenn-Heat-Transfer (Glenn-HT) Computer Code to Enable Time-Filtered Navier Stokes (TFNS) Simulations and Application to Film Cooling on a Flat Plate Through Long Cooling Tubes

Computational fluid dynamics (CFD) analysis using Reynolds-averaged Navier-Stokes (RANS) formulation for turbomachinery-related flows has enabled improved engine component designs. RANS methodology has limitations that are related to its inability to accurately describe the spectrum of flow phenomena encountered in engines. Examples of flows that are difficult to compute accurately with RANS include phenomena such as laminar/turbulent transition, turbulent mixing due to mixing of streams, and separated flows. Large eddy simulation (LES) can improve accuracy but at a considerably higher cost. In recent years, hybrid schemes that take advantage of both unsteady RANS and LES have been proposed. This study investigated an alternative scheme, the time-filtered Navier-Stokes (TFNS) method applied to compressible flows. The method developed by Shih and Liu was implemented in the Glenn-Heat-Transfer (Glenn-HT) code and applied to film-cooling flows. In this report the method and its implementation is briefly described. The film effectiveness results obtained for film cooling from a row of 30deg holes with a pitch of 3.0 diameters emitting air at a nominal density ratio of unity and two blowing ratios of 0.5 and 1.0 are shown. Flow features under those conditions are also described.

heat transfer↗

Numerical simulation of compressor endwall and casing treatment flow phenomena

A numerical study is presented of the flow in the endwall region of a compressor blade row, in conditions of operation with both smooth and grooved endwalls. The computations are first compared to velocity field measurements in a cantilevered stator/rotating hub configuration to confirm that the salient features are captured. Computations are then interrogated to examine the tip leakage flow structure since this is a dominant feature of the endwall region. In particular, the high blockage that can exist near the endwalls at the rear of a compressor blade passage appears to be directly linked to low total pressure fluid associated with the leakage flow. The fluid dynamic action of the grooved endwall, representative of the casing treatments that have been most successful in suppressing stall, is then simulated computationally and two principal effects are identified. One is suction of the low total pressure, high blockage fluid at the rear of the passage. The second is energizing of the tip leakage flow, most notably in the core of the leakage vortex, thereby suppressing the blockage at its source.

Crook, A. J.↗

Controls/CFD Interdisciplinary Research Software Generates Low-Order Linear Models for Control Design From Steady-State CFD Results

The NASA Lewis Research Center is developing analytical methods and software tools to create a bridge between the controls and computational fluid dynamics (CFD) disciplines. Traditionally, control design engineers have used coarse nonlinear simulations to generate information for the design of new propulsion system controls. However, such traditional methods are not adequate for modeling the propulsion systems of complex, high-speed vehicles like the High Speed Civil Transport. To properly model the relevant flow physics of high-speed propulsion systems, one must use simulations based on CFD methods. Such CFD simulations have become useful tools for engineers that are designing propulsion system components. The analysis techniques and software being developed as part of this effort are an attempt to evolve CFD into a useful tool for control design as well. One major aspect of this research is the generation of linear models from steady-state CFD results. CFD simulations, often used during the design of high-speed inlets, yield high resolution operating point data. Under a NASA grant, the University of Akron has developed analytical techniques and software tools that use these data to generate linear models for control design. The resulting linear models have the same number of states as the original CFD simulation, so they are still very large and computationally cumbersome. Model reduction techniques have been successfully applied to reduce these large linear models by several orders of magnitude without significantly changing the dynamic response. The result is an accurate, easy to use, low-order linear model that takes less time to generate than those generated by traditional means. The development of methods for generating low-order linear models from steady-state CFD is most complete at the one-dimensional level, where software is available to generate models with different kinds of input and output variables. One-dimensional methods have been extended somewhat so that linear models can also be generated from two- and three-dimensional steady-state results. Standard techniques are adequate for reducing the order of one-dimensional CFD-based linear models. However, reduction of linear models based on two- and three-dimensional CFD results is complicated by very sparse, ill-conditioned matrices. Some novel approaches are being investigated to solve this problem.

Melcher, Kevin J.↗

Axisymmetric Plume Simulations with NASA's DSMC Analysis Code

A comparison of axisymmetric Direct Simulation Monte Carlo (DSMC) Analysis Code (DAC) results to analytic and Computational Fluid Dynamics (CFD) solutions in the near continuum regime and to 3D DAC solutions in the rarefied regime for expansion plumes into a vacuum is performed to investigate the validity of the newest DAC axisymmetric implementation. This new implementation, based on the standard DSMC axisymmetric approach where the representative molecules are allowed to move in all three dimensions but are rotated back to the plane of symmetry by the end of the move step, has been fully integrated into the 3D-based DAC code and therefore retains all of DAC s features, such as being able to compute flow over complex geometries and to model chemistry. Axisymmetric DAC results for a spherically symmetric isentropic expansion are in very good agreement with a source flow analytic solution in the continuum regime and show departure from equilibrium downstream of the estimated breakdown location. Axisymmetric density contours also compare favorably against CFD results for the R1E thruster while temperature contours depart from equilibrium very rapidly away from the estimated breakdown surface. Finally, axisymmetric and 3D DAC results are in very good agreement over the entire plume region and, as expected, this new axisymmetric implementation shows a significant reduction in computer resources required to achieve accurate simulations for this problem over the 3D simulations.

Stewart, B. D.↗

Parametric Modeling Investigation of a Radially-Staged Low-Emission Aviation Combustor

Aviation gas-turbine combustion demands high efficiency, wide operability and minimal trace gas emissions. Performance critical design parameters include injector geometry, combustor layout, fuel-air mixing and engine cycle conditions. The present investigation explores these factors and their impact on a radially staged low-emission aviation combustor sized for a next-generation 24,000-lbf-thrust engine. By coupling multi-fidelity computational tools, a design exploration was performed using a parameterized annular combustor sector at projected 100% takeoff power conditions. Design objectives included nitrogen oxide emission indices and overall combustor pressure loss. From the design space, an optimal configuration was selected and simulated at 7.1, 30 and 85% part-power operation, corresponding to landing-takeoff cycle idle, approach and climb segments. All results were obtained by solution of the steady-state Reynolds-averaged Navier-Stokes equations. Species concentrations were solved directly using a reduced 19-step reaction mechanism for Jet-A. Turbulence closure was obtained using a nonlinear K-epsilon model. This research demonstrates revolutionary combustor design exploration enabled by multi-fidelity physics-based simulation.

design analysis↗

Unstructured Grid Adaptation: Status, Potential Impacts, and Recommended Investments Toward CFD Vision 2030

Unstructured grid adaptation is a powerful tool to control discretization error for Computational Fluid Dynamics (CFD). It has enabled key increases in the accuracy, automation, and capacity of some fluid simulation applications. Slotnick et al. provides a number of case studies in the CFD Vision 2030 Study: A Path to Revolutionary Computational Aerosciences to illustrate the current state of CFD capability and capacity. The authors forecast the potential impact of emerging High Performance Computing (HPC) environments forecast in the year 2030 and identify that mesh generation and adaptivity continue to be significant bottlenecks in the CFD work flow. These bottlenecks may persist because very little government investment has been targeted in these areas. To motivate investment, the impacts of improved grid adaptation technologies are identified. The CFD Vision 2030 Study roadmap and anticipated capabilities in complementary disciplines are quoted to provide context for the progress made in grid adaptation in the past fifteen years, current status, and a forecast for the next fifteen years with recommended investments. These investments are specific to mesh adaptation and impact other aspects of the CFD process. Finally, a strategy is identified to diffuse grid adaptation technology into production CFD work flows.

Park, Michael A.↗

Vibration dissociation coupling in nonequilibrium flows

The final report on research between North Carolina State University and the NASA Ames Research Center is presented. The research was aimed at using the Schwartz, Slawsky, Herzfeld (SSH) theory to simulate the vibrational relaxation of nitrogen molecules undergoing dissociation or recombination over a wide range of conditions. The results of these simulations were then treated as exact, and they were used to develop a model for the coupled vibration-dissociation process. This new model is simple enough to be used in computational fluid dynamics codes, but still captures the physics of the complex process. The model is used to simulate the flow over typical geometries to test it and to determine how much impact it has on the flow field. The key elements of this research are summarized.

Candler, Graham V.↗