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,081 records · Page 60

Analysis of internal flows relative to the space shuttle main engine

Cooperative efforts between the Lockheed-Huntsville Computational Mechanics Group and the NASA-MSFC Computational Fluid Dynamics staff has resulted in improved capabilities for numerically simulating incompressible flows generic to the Space Shuttle Main Engine (SSME). A well established and documented CFD code was obtained, modified, and applied to laminar and turbulent flows of the type occurring in the SSME Hot Gas Manifold. The INS3D code was installed on the NASA-MSFC CRAY-XMP computer system and is currently being used by NASA engineers. Studies to perform a transient analysis of the FPB were conducted. The COBRA/TRAC code is recommended for simulating the transient flow of oxygen into the LOX manifold. Property data for modifying the code to represent LOX/GOX flow was collected. The ALFA code was developed and recommended for representing the transient combustion in the preburner. These two codes will couple through the transient boundary conditions to simulate the startup and/or shutdown of the fuel preburner. A study, NAS8-37461, is currently being conducted to implement this modeling effort.

Source record↗

Numerical Analysis of Combustion Dynamics in a Full-Scale Rotating Detonation Rocket Engine using Large Eddy Simulations

Large eddy simulations (LES) using detailed chemistry and leveraging adaptive mesh refinement (AMR) are performed to gain insights into the combustion dynamics within a full-scale methane-oxygen non-premixed rotating detonation rocket engine (RDRE) employing impinging discrete injection schemes. In particular, a comparative analysis of two operating conditions corresponding to the same global equivalence ratio but different mass flow rates is carried out to investigate the resultant impact on detonation wave characteristics and RDRE global performance. Multiple co-rotating detonation waves with spatially-distributed wave structure and preferential alignment with the inner wall of the annulus (due to asymmetry in fuel distribution) are encountered under both conditions. Both cases exhibit pre-detonation deflagrative burning in the fill region, while one of the cases shows higher susceptibility to backflow into the feed plenums due to lower plenum pressures. Furthermore, heat release analysis shows that the thrust obtained from the RDRE is closely linked to the distribution of total heat release between detonative and deflagrative combustion. On the other hand, combustion efficiency is associated with the fraction of heat release occurring in fuel-rich versus fuel-lean regions within the RDRE.

33 ADVANCED PROPULSION SYSTEMS↗

Simulating the filamentation of smoothed laser beams with three-dimensional nonlinear dynamics

In a plasma, the ponderomotive force of an inhomogeneous electromagnetic field expels plasma from regions of high intensity. When a laser propagates through a plasma, this force creates density wells that subsequently modify the index of refraction. The beam refracts and focuses into these wells and may filament. In extreme cases, the laser beam will spray due to increasing angular divergence of the beam. The threshold for ponderomotive self-focusing is well established for isolated laser hotspots or speckles. Here, we define a practical threshold for characterizing the filamentation of thousands of speckles that are found in the focal plane of high-power laser beams spatially smoothed with random phase plates as used at high energy and power laser facilities studying inertial confinement fusion. This threshold is tested against three-dimensional simulations of speckled laser light propagating through plasma. Four metrics are applied to assess filamentation: the fraction of power above five times the average intensity, an effective f-number, the mean-squared perpendicular wavenumber, and the fraction of rarefied density with deviation from the initial density exceeding | δn/ n| = 0.1. The speckled beams studied are generated by random phase plates, both with and without additional polarization smoothing, in a parameter regime of relevance to indirect drive experiments. While filamentation has been discussed extensively in the literature, we believe this to be the first published simulation study with three-dimensional nonlinear hydrodynamics that addresses the onset threshold of ponderomotive filamentation and establishes the lengths and time scales necessary to reach a statistical steady state.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Computational strategies for three-dimensional flow simulations on distributed computer systems

An increasing amount of research activity in computational fluid dynamics has been devoted to the development of efficient algorithms for parallel computing systems. The increasing performance to price ratio of engineering workstations has led to research to development procedures for implementing a parallel computing system composed of distributed workstations. This thesis proposal outlines an ongoing research program to develop efficient strategies for performing three-dimensional flow analysis on distributed computing systems. The PVM parallel programming interface was used to modify an existing three-dimensional flow solver, the TEAM code developed by Lockheed for the Air Force, to function as a parallel flow solver on clusters of workstations. Steady flow solutions were generated for three different wing and body geometries to validate the code and evaluate code performance. The proposed research will extend the parallel code development to determine the most efficient strategies for unsteady flow simulations.

Weed, Richard Allen↗

Training Knowledge Bots for Physics-Based Simulations Using Artificial Neural Networks

Millions of complex physics-based simulations are required for design of an aerospace vehicle. These simulations are usually performed by highly trained and skilled analysts, who execute, monitor, and steer each simulation. Analysts rely heavily on their broad experience that may have taken 20-30 years to accumulate. In addition, the simulation software is complex in nature, requiring significant computational resources. Simulations of system of systems become even more complex and are beyond human capacity to effectively learn their behavior. IBM has developed machines that can learn and compete successfully with a chess grandmaster and most successful jeopardy contestants. These machines are capable of learning some complex problems much faster than humans can learn. In this paper, we propose using artificial neural network to train knowledge bots to identify the idiosyncrasies of simulation software and recognize patterns that can lead to successful simulations. We examine the use of knowledge bots for applications of computational fluid dynamics (CFD), trajectory analysis, commercial finite-element analysis software, and slosh propellant dynamics. We will show that machine learning algorithms can be used to learn the idiosyncrasies of computational simulations and identify regions of instability without including any additional information about their mathematical form or applied discretization approaches.

Samareh, Jamshid A.↗

Preparing an Incompressible-Flow Fluid Dynamics Code for Exascale-Class Wind Energy Simulations

The U.S. Department of Energy has identified exascale-class wind farm simulation as critical to wind energy scientific discovery. A primary objective of the ExaWind project is to build high-performance, predictive computational fluid dynamics (CFD) tools that satisfy these modeling needs. GPU accelerators will serve as the computational thoroughbreds of next-generation, exascale-class supercomputers. Here, we report on our efforts in preparing the ExaWind unstructured mesh solver, Nalu-Wind, for exascale-class machines. For computing at this scale, a simple port of the incompressible-flow algorithms to GPUs is insufficient. To achieve high performance, one needs novel algorithms that are application aware, memory efficient, and optimized for the latest-generation GPU devices. The result of our efforts are unstructured-mesh simulations of wind turbines that can effectively leverage thousands of GPUs. In particular, we demonstrate a first-of-its-kind, incompressible-flow simulation using Algebraic Multigrid solvers that strong scales to more than 4000 GPUs on the Summit supercomputer.

algebraic multigrid↗

Computational Fluid Dynamics Analysis for the Orbiter LH2 Feedline Flowliner

In phase II, additional inducer rotations are simulated in order to understand the root cause of the flowliner crack problem. CFD results confirmed that there is a strong unsteady interaction between the backflow regions caused by the LPFTP inducer and secondary flow regions in the bellows cavity through the flowliner slots. It is observed that the swirl on the duct side of the downstream flowliner is stronger than on the duct side of the upstream flowliner. Due to this swirl, there are more significant unsteady flow interactions through the downstream slots than those observed in the upstream slots. Averaged values of the local velocity at the slots were provided to the NESC-ITA flow physics acoustics team to guide them in designing the acoustics experiment. A parametric study was performed to compare the flow field in the flowliner area when one upstream slot and one corresponding downstream slot were enlarged. No significant differences were observed between the flow field obtained from the enlarged slot configuration when compared with the original configuration. More cases must be analyzed with various enlarged slot configurations to generalize this observation. The flow through the A1 test stand and the flow through the orbiter fuel feedline manifold were simulated without the LPFTP. It was observed that incoming flow to the flowliner and inducer was more uniform in the A1 test stand then in the orbiter manifold. Additionally, each engine LPFTP in the orbiter receives significantly different velocity distributions. Because of the differences observed in the computed results, it is not possible for the A1 test stand to represent the three different engine feedlines simultaneously.

Kiris, Cetin C.↗

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 investigated include JP-8, synthetic fuel, and two blends 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 multiphase, 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 modeling of natural circulation in LiCl-KCl molten salt closed loop

Characterizing flow within a molten salt closed-loop system is crucial for assessing system requirements, evaluating performance, and identifying potential flaws. Direct flow measurement using instrumentation is challenging due to extreme environmental conditions and the limitations associated with measuring molten salt flow under natural convection. Here, this study aims to provide comprehensive insights into the thermal-hydraulic behavior of a closed loop, with a particular focus on temperature distribution and velocity prediction. The Computational Fluid Dynamics (CFD) model demonstrated the capability to effectively simulate and predict both temperature distributions and flow velocities within the molten salt loop. The CFD model's predictive capability was validated by its ability to replicate temperature measurements under varying boundary conditions. The analysis revealed that the CFD model tends to underpredict temperatures in the cold leg and overpredict them in the hot leg, highlighting the need for continuous model refinement and acknowledging the limitations of using a steady-state approach. Furthermore, the potential of using external temperature measurements to estimate internal molten salt temperatures and predict flow velocity was explored, revealing that this approach could introduce up to a 5.5% error in flow velocity calculations. Line probes mapping temperature distributions across the tube's cross-section and molten salt provided valuable insights into temperature gradients, emphasizing the need for a thermal conductivity equation for molten salt with lower uncertainty to achieve more accurate temperature predictions of the system.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

Performance analysis and tank test validation of a hybrid ocean wave-current energy converter with a single power takeoff

This paper introduces a hybrid ocean wave-current energy converter (HWCEC) that harvests energy from ocean waves and current simultaneously with a single power takeoff (PTO). Specifically, the wave energy is extracted through relative heaving motion between a floating buoy and a submerged second body, while the current energy is extracted using a marine current turbine (MCT). Energy from both sources are integrated by a hybrid PTO whose concept is based on a mechanical motion rectifier (MMR). The hybrid PTO with two one-way clutches converts bidirectional, up-and-down motion from the waves into unidirectional rotation of the generator. Meanwhile a third one-way clutch couples the MCT with the same generator. The wave and current can simultaneously or separately drive the same generator through different engagement and disengagement statuses of the one-way clutches. Time-domain simulation is conducted with hydrodynamic coefficients obtained from computational fluid dynamics analysis and boundary element method. Tank tests were conducted for a HWCEC under co-existing wave and current inputs. For comparison, separate baseline tests of a turbine and a two-body point absorber, each acting in isolation, are conducted. Experimental results validate the dynamic modeling and show that a HWCEC can increase the output power with a range between 29-87% over either current turbine and wave energy converter acting individually, and it can reduce by up to 70% the peak-to-average power ratio compared with the wave energy converter on the tested conditions.

16 TIDAL AND WAVE POWER↗

Assessment of Lost Circulation Material Particle-Size Distribution on Fracture Sealing: A Numerical Study

Summary Lost circulation materials (LCMs) are essential to combat fluid loss while drilling and may put the whole operation at risk if a proper LCM design is not used. The focus of this research is understanding the function of LCMs in sealing fractures to reduce fluid loss. One important consideration in the success of fracture sealing is the particle-size distribution (PSD) of LCMs. Various studies have suggested different guidelines for obtaining the best size distribution of LCMs for effective fracture sealing based on limited laboratory experiments or field observations. Hence, there is a need for sophisticated numerical methods to improve the LCM design by providing some predictive capabilities. In this study, computational fluid dynamics (CFD) and discrete element methods (DEM) numerical simulations are coupled to investigate the influence of PSD of granular LCMs on fracture sealing. Dimensionless variables were introduced to compare cases with different PSDs. We validated the CFD-DEM model in reproducing specific laboratory observations of fracture-sealing experiments within the model boundary parameters. Our simulations suggested that a bimodally distributed blend would be the most effective design in comparison to other PSDs tested here.

Engineering↗

Calculation of confined swirling flows

Swirling jets are an important constituent flow of many types of combustion equipment. Flame size, shape, stability and combustion intensity are all favorably influenced by swirling some portion of the incoming air with the fuel. Fluid dynamics computer codes are being developed for the mathematical simulation of flows in practical combustors. It is important that these codes have the ability to accurately and reliably calculate swirling flows. Turbulent flow calculations with swirl can be inaccurate due to at least three major reasons: limitations of the turbulence modeling used, incorrect or inadequate specifications of inlet boundary conditions, and, error introduced through the numerics. These sources of error are described, and examples provided of each. The state-of-the-art in such calculations is reviewed. The controlling character of current numerics is demonstrated, and it is concluded that as a first priority, better numerics must be arrived at, and that an improved understanding of the discharge flow from swirl generators is essential as a second priority. When these priorities are satisfied an improved turbulence model is desirable. It is also concluded that an improved understanding of swirl-induced recirculation would be an asset.

Sturgess, G. J.↗

Shuttle entry aerothermodynamic flight research - The Orbiter Experiments (OEX) Program

Results of the OEX program are summarized with emphasis on the information on entry aerothermodynamic phenomena derived from Space Shuttle operations. The discussion focuses on OEX experiment complement and operational history, freestream environment and vehicle attitude data, aerodynamic force and moment data, aerodynamic surface data, and vehicle configuration data. Attention is also given to orbiter aerodynamic performance, stability and control, high-altitude atmospheric density variability, direct simulation Monte Carlo validation, orbital drag variation, and computational fluid dynamic technique validation.

Throckmorton, David A.↗

Transonic Investigation of Two-Dimensional Nozzles Designed for Supersonic Cruise

An experimental and computational investigation has been conducted to determine the off-design uninstalled drag characteristics a of a two-dimensional convergent-divergent nozzle designed for a supersonic cruise civil transport. The main objective of this investigation was to determine the effects of varying nozzle external flap curvature and sidewall boattail angle and curvature on nozzle drag The experimental investigation was conducted in the Langley 16-Foot Transonic Tunnel at Mach numbers from 0.80 to 1.20 at nozzle pressure ratios up to nine. Three-dimensional simulations of nozzle performance were obtained with the computational fluid dynamics code PAB using turbulence closure and nonlinear Reynolds stress modeling. The results of this investigation indicate that excellent correlation between experimental and predicted results was obtained for the nozzle with a moderate amount of boattail curvature. The nozzle with an external flap having a sharp shoulder (no curvature) had the lowest nozzle pressure drag. At a Mach number of 1.2, sidewall pressure drag doubled as sidewall boattail angle was increased from 4 to 8 deg. Reducing the height of the sidewall caused large decreases in both the sidewall and flap pressure drags.

Capone, Francis J.↗

Analysis of X-33 Linear Aerospike Plume Induced Base-Heating Environment During Power-Pack Out

The objectives of this viewgraph presentation are to predict the following: (1) dual-engine base-heating at 57% PL at sea level, and (2) dual-engine base-heating during PPO at three ascent abort trajectories. A systematically anchored computational fluid dynamics and heat transfer three-dimensional transfer simulation is being used to study the effect of reduced power levels on base-heating environments during sea level testing and during PPO. Preliminary results show the following: (1) convective heating is higher for the 57% PL than for 100% PL on most of the pillows and flex seals during sea level testing; and (2) convective heating on pillows and flex seals on the 'off' engine side is higher than that on the 'on' engine side.

Wang, Ten-See↗

Method for CFD Simulation of Propellant Slosh in a Spherical Tank

Propellant sloshing can impart unwanted disturbances to spacecraft, especially if the spacecraft controller is driving the system at the slosh frequency. This paper describes the work performed by the authors in simulating propellant slosh in a spherical tank using computational fluid dynamics (CFD). ANSYS-CFX is the CFD package used to perform the analysis. A 42 in spherical tank is studied with various fill fractions. Results are provided for the forces on the walls and the frequency of the slosh. Snapshots of slosh animation give a qualitative understanding of the propellant slosh. The results show that maximum slosh forces occur at a tank fill fraction of 0.4 and 0.6 due to the amount of mass participating in the slosh and the room available for sloshing to occur. The slosh frequency increases as the tank fill fraction increases.

Benson, David J.↗

Simulating Effects of High Angle of Attack on Turbofan Engine Performance

A method of investigating the effects of high angle of attack (AOA) flight on turbofan engine performance is presented. The methodology involves combining a suite of diverse simulation tools. Three-dimensional, steady-state computational fluid dynamics (CFD) software is used to model the change in performance of a commercial aircraft-type inlet and fan geometry due to various levels of AOA. Parallel compressor theory is then applied to assimilate the CFD data with a zero-dimensional, nonlinear, dynamic turbofan engine model. The combined model shows that high AOA operation degrades fan performance and, thus, negatively impacts compressor stability margins and engine thrust. In addition, the engine response to high AOA conditions is shown to be highly dependent upon the type of control system employed.

Liu, Yuan↗