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,567 records · Page 87

Newton-Krylov-Schwarz: An implicit solver for CFD

Newton-Krylov methods and Krylov-Schwarz (domain decomposition) methods have begun to become established in computational fluid dynamics (CFD) over the past decade. The former employ a Krylov method inside of Newton's method in a Jacobian-free manner, through directional differencing. The latter employ an overlapping Schwarz domain decomposition to derive a preconditioner for the Krylov accelerator that relies primarily on local information, for data-parallel concurrency. They may be composed as Newton-Krylov-Schwarz (NKS) methods, which seem particularly well suited for solving nonlinear elliptic systems in high-latency, distributed-memory environments. We give a brief description of this family of algorithms, with an emphasis on domain decomposition iterative aspects. We then describe numerical simulations with Newton-Krylov-Schwarz methods on aerodynamics applications emphasizing comparisons with a standard defect-correction approach, subdomain preconditioner consistency, subdomain preconditioner quality, and the effect of a coarse grid.

Cai, Xiao-Chuan↗

Modeling of Supersonic Film Cooling on the J-2X Nozzle Extension

Supersonic film cooling (SSFC) of nozzles has been used in several liquid rocket engine designs, and is being applied to the nozzle extension (NE) of the J-2X upper stage engine currently under development. Turbine exhaust gas (TEG) is injected tangentially from a manifold along the NE, and provides a thermal barrier from the core nozzle flow for the NE. As the TEG stream mixes with the nozzle flow, the effectiveness of the thermal barrier is reduced. This paper documents computational fluid dynamics (CFD) analysis work performed by NASA Marshall Space Flight Center (MSFC) to model the flow of the TEG through the manifold, into the nozzle, and the subsequent mixing of the TEG stream with the core flow. The geometry and grid of the TEG manifold, structural support ribs, and the NE wall will be shown, and the CFD boundary conditions described. The Loci-CHEM CFD code used in this work will also be briefly described. A unique approach to modeling the combined TEG manifold/thrust chamber assembly (TCA) was employed, as it was not practical to model the entire 360 circumferential range in one simulation. Prior CFD validation work modeling Calspan SSFC experiments in the early 1990s, documented in a previous AIAA paper, will also be briefly discussed. The fluid dynamics of the TEG flow through the manifold, into and between the structural support ribs, and into the nozzlette that feeds the TCA will be described. Significant swirl and non-uniformities are present, which along with the wakes from the ribs, act to degrade the film cooling effectiveness compared to idealized injection of TEG gas. The effect of these flow characteristics on the adiabatic wall temperature profile on the NE will be discussed.

Ruf, Joseph H.↗

Propulsion-Airframe Integration for Conceptual Redesign of a Low-Boom Supersonic Transport

A low-boom supersonic transport was generated in a previous low-boom multidisciplinary optimization (MDO) study that used computational fluid dynamics (CFD) off-body pressure to compute the undertrack sonic boom ground signature and calibrated low-fidelity aerodynamics analyses to compute the mission performance metrics. This low-boom aircraft, referred to as the Mach 1.7 40-PAX concept, can carry 40 passengers for a low-boom overland mission with cruise Mach 1.7 for the airport pairs over the continental US and has the potential to achieve an undertrack sonic boom ground noise level below 70 perceived level of decibels at the start of overland cruise (SOC). The engine for the Mach 1.7 40-PAX concept was designed using the Numerical Propulsion System Simulation (NPSS) for the mission analysis and modeled as a flow-through nacelle for CFD-based sonic boom analysis. To understand how the engine plume affects the undertrack ground signature, the Mach 1.7 40-PAX concept is redesigned in this paper after replacing the flow-through nacelles by CFD engines for sonic boom analysis using aeropropulsive CFD simulation. A process is developed for approximation of the NPSS engine at SOC by a CFD engine for aeropropulsive CFD simulation. The generated CFD engine has the identical nozzle boundary conditions and approximately the same mass flow and thrust as those of the NPSS engine at SOC. Then, the outer mold line of the configuration with the CFD engines is redesigned to approximately restore the low-boom characteristics of the Mach 1.7 40-PAX concept at SOC. Finally, the CFD simulation data for the redesigned concept with the CFD engines is used to calibrate the low-fidelity aerodynamic analyses for mission analysis of this concept. The cyclic dependency of the involved disciplinary analyses is resolved using an iteration method for a consistent coupling of the mission analysis, NPSS engine analysis, and low-boom redesign using the aeropropulsive CFD simulation. This low-boom redesign study is used as an example to demonstrate how the propulsion-airframe integration could be implemented for conceptual design of supersonic transports that satisfy both the low-boom and mission performance requirements.

Low-boom supersonic transport↗

Validation of a Two-Phase CFD Model for Predicting Tank Self-Pressurization in the Ground-Based K-Site Experiment

A two-phase CFD model for self-pressurization of a cryogenic storage tank partially filled with liquid hydrogen is presented using the Volume-Of-Fluid approach for modeling two-phase flow, as well as interfacial heat, mass and momentum transfer between the liquid and the vapor regions. The CFD model is validated against self-pressurization experiment performed using the K-site flightweight hydrogen storage tank at NASA Glenn Research Center1. Laminar and turbulent simulations are performed together with conjugated heat transfer analysis. Effects of turbulence, the value of accommodation coefficient used for predicting phase change rates, as well as tank wall geometry are presented and discussed. Predicted tank pressures, fluid and wall temperatures are compared with the experimental data at 49% fill level and two different heat loads for validating this CFD model.

Computational Fluid Dynamics↗

Validation of a Two-Phase CFD Model for Predicting Tank Self-Pressurization in the Ground-Based K-Site Experiment

A two-phase CFD model for self-pressurization of a cryogenic storage tank partially filled with liquid hydrogen is presented using the Volume-Of-Fluid approach for modeling two-phase flow, as well as interfacial heat, mass and momentum transfer between the liquid and vapor regions. The CFD model is validated against self-pressurization experiment performed using the K-site flightweight hydrogen storage tank at NASA Glenn Research Center 1 . Laminar and turbulent simulations are performed together with conjugate heat transfer analysis. Effects of turbulence, the value of accommodation coefficient used for predicting phase change rates, as well as tank wall geometry are presented and discussed. Predicted tank pressures, fluid and wall temperatures are compared with the experimental data at 49% fill level and two different heat loads for validating this CFD model.

Computational Fluid Dynamics↗

Verification and Validation Studies for the LAVA CFD Solver

The verification and validation of the Launch Ascent and Vehicle Aerodynamics (LAVA) computational fluid dynamics (CFD) solver is presented. A modern strategy for verification and validation is described incorporating verification tests, validation benchmarks, continuous integration and version control methods for automated testing in a collaborative development environment. The purpose of the approach is to integrate the verification and validation process into the development of the solver and improve productivity. This paper uses the Method of Manufactured Solutions (MMS) for the verification of 2D Euler equations, 3D Navier-Stokes equations as well as turbulence models. A method for systematic refinement of unstructured grids is also presented. Verification using inviscid vortex propagation and flow over a flat plate is highlighted. Simulation results using laminar and turbulent flow past a NACA 0012 airfoil and ONERA M6 wing are validated against experimental and numerical data.

Validation↗

Large Vehicle Lunar Landing Surface Interaction and In-Situ Resource Based Risk Mitigation: Landing & Launch Pads

A key capability required for the exploration of planetary bodies is the ability to land on the surface. Previous work performed by NASA and other institutions has primarily focused on landing small spacecraft on planetary surfaces and the associated small-to-medium thrusters required for the soft landing. In the case of human exploration—particularly the establishment of long duration exploration and habitation outposts—the ability to land large landers, such as the SpaceX Starship, is necessary. These larger landing systems require the use of more powerful engines, with higher engine exhaust temperatures and higher landing loads. Understanding the excavation of material by the engines, as well as the potential for the landing legs to sink into the subsurface, is key in ensuring reliable and safe landings. A further improvement in landing reliability can be achieved by constructing landing / launch pads, especially with in-situ resources. Some material excavation by the plume is inevitable, leaving at least a portion of the surface scoured and uneven under the lander and ejecting regolith particles and rocks at very high velocities. One possible solution would be to robotically build landing / launch pads (ideally autonomously) at the destination using in-situ materials. In this case, the first one or few landers will need to land on unimproved surfaces at higher risk; however, they would bring the required equipment to build the landing pads with mostly local resources, thus increasing the reliability of safe landing for subsequent larger landers. A number of methods to build in-situ landing and launch pads have already been developed. These methods include no, or some, addition of required binder additives to the local regolith material, different processing approaches and result in varying landing pad strengths. A sub-scale rocket engine plume, was used to simulate some of the conditions of a landing on the Moon to assess the effectiveness of various materials for an in-situ built landing pad, The GO2/GCH4rocket engine fired on a 1m2area coupons of representative pad materials. The results will allow continued development towards materials that satisfy the landing pad properties required for the effective risk reduction and increased reliability for landing people and equipment on the lunar surface. This work contained two parts: (1) computer modeling of a large rocket engine plume interacting with regolith on the Moon, using the Granular Gas Flow Solver (GGFS) provided by CFD Research Corporation as well as other computational fluid dynamics codes (CFD) such as Loci/CHEM. (2) Developing landing/launch pad materials that could be used for in-situ construction on the lunar surface in the future, to mitigate the calculated effects of a large vehicle rocket engine landing and launching on the Moon.

Lunar↗

Improved Models for Simulation of Cryogenic Spray Cooling for Propellant Tank Chilldown

This paper presents advances in spray boiling models for cryogenic fluids that exhibit strong variation in thermodynamic properties; it accounts for the interaction of the spray droplets with walls and the complex heat transfer and boiling processes that play a vital role in chilling solid surfaces. The model has been developed by corroborating and validating it with test data in liquid nitrogen obtained at the University of Connecticut. The test data includes detailed visualization and plate temperature measurements for liquid nitrogen sprays impinging on a warm test plate and subsequent transient chilldown of the plate. Simulations of transient plate chilldown in liquid nitrogen were carried out and temperature histories at various locations on the plate were compared with test measurements to demonstrate the validity of the spray boiling models that have been developed for cryogens.

Boiling↗

Dynamical behavior of surface tension on rotating fluids in low and microgravity environments

Consideration is given to the time-dependent evolutions of the free surface profile (bubble shapes) of a cylindrical container, partially filled with a Newtonian fluid of constant density, rotating about its axis of symmetry in low and microgravity environments. The dynamics of the bubble shapes are calculated for four cases: linear time-dependent functions of spin-up and spin-down in low and microgravity, linear time-dependent functions of increasing and decreasing gravity at high and low rotating cylinder speeds, time-dependent step functions of spin-up and spin-down in low gravity, and sinusoidal function oscillation of the gravity environment in high and low rotating cylinder speeds. It is shown that the computer algorithms developed by Hung et al. (1988) may be used to simulate the profile of time-dependent bubble shapes under variations of centrifugal, capillary, and gravity forces.

Hung, R. J.↗

Computer Modeling of Sand Transport on Mars Using a Compart-Mentalized Fluids Algorithm (CFA)

It has been postulated that aeolian transport on Mars may be significantly different from that on Earth. From laboratory experiments simulating martian grain transport [2], it has been observed that (saltating) grains striking the bed can cause hundreds of secondary reptation trajectories when impact occurs at speeds postulated for Mars. Some of the ballistically induced trajectories "die ouf' and effectively join the ranks on the creep population that is merely nudged along by impact. Many of the induced reptation trajectories, however, are sufficiently high for the grains to become part of the saltation load (it is irrelevant to the boundary layer how a grain attained its initial lift force). When these grains, in turn, strike the surface, they too are capable of inducing more reptating grains. This cascading effect has been discussed in connection with terrestrial aeolian transport in an attempt to dispel the notion that sand motion is divisible only into creep and saltation loads. On Earth, only a few grains are splashed by impact. On Mars, it may be hundreds. We developed a computer model to address this phenomenon because there are some important ramifications: First, this ratio may mean that martian aeolian transport is dominated by reptation flux rather than saltation. On Earth, the flux would be a roughly balanced mixture between reptation/creep and saltation. On Venus, there would be no transport other than by saltation. In other words, an understanding of planetary aeolian processes may not be necessarily understood by extrapolating from the "Earth case", with only gravity and atmospheric density/viscosity being considered as variables. Second, the reptation flux on Mars may be self sustaining, so that little input is required by the wind once transport has been initiated. The number of grains saturating the boundary layer near the bed may mean that average grain speed on Mars might conceivably be less than that on Earth. This would say much for models of sand comminution on Mars. A multiple-grain transport model using just the equations of grain motion describing lift and drag is impossible to develop owing to stochastic effects --the very effects we wish to model. Also, unless we were to employ supercomputing techniques and extremely complex computer codes that could deal with millions of grains simultaneously, it would also be difficult to model grain transport if we attempted to consider every grain in motion. No existing computer models were found that satisfactorily used the equations of motion to arrive at transport flux numbers for the different populations of saltation and reptation. Modeling all the grains in a transport system was an intractable problem within our resources, and thus we developed what we believe to be a new modeling approach to simulating grain transport. The CFA deals with grain populations, but considers them to belong to various compartmentalized fluid units in the boundary layer. In this way, the model circumvents the multigrain problem by dealing primarily with the consequences of grain transport --momentum transfer between air and grains, which is the physical essence of a dynamic grain-fluid mixture. We thus chose to model the aeolian transport process as a superposition of fluids. These fluids include the air as well as particle populations of various properties. The prime property distinguishing these fluids is upward and downward grain motion. In a normal saltation trajectory, a grain's downwind velocity increases with time, so a rising grain will have a smaller downwind velocity than a failing grain. Because of this disparity in rising and falling grain proper-ties, it seemed appropriate to track these as two separate grain populations within the same physical space. The air itself can be considered a separate fluid superimposed within and interacting with the various grain-cloud "fluids". Additional informaiton is contained in the original.

Marshall, J.↗

Relaxation and Preconditioning for High Order Discontinuous Galerkin Methods with Applications to Aeroacoustics and High Speed Flows

This project is about the investigation of the development of the discontinuous Galerkin finite element methods, for general geometry and triangulations, for solving convection dominated problems, with applications to aeroacoustics. Other related issues in high order WENO finite difference and finite volume methods have also been investigated. methods are two classes of high order, high resolution methods suitable for convection dominated simulations with possible discontinuous or sharp gradient solutions. In [18], we first review these two classes of methods, pointing out their similarities and differences in algorithm formulation, theoretical properties, implementation issues, applicability, and relative advantages. We then present some quantitative comparisons of the third order finite volume WENO methods and discontinuous Galerkin methods for a series of test problems to assess their relative merits in accuracy and CPU timing. In [3], we review the development of the Runge-Kutta discontinuous Galerkin (RKDG) methods for non-linear convection-dominated problems. These robust and accurate methods have made their way into the main stream of computational fluid dynamics and are quickly finding use in a wide variety of applications. They combine a special class of Runge-Kutta time discretizations, that allows the method to be non-linearly stable regardless of its accuracy, with a finite element space discretization by discontinuous approximations, that incorporates the ideas of numerical fluxes and slope limiters coined during the remarkable development of the high-resolution finite difference and finite volume schemes. The resulting RKDG methods are stable, high-order accurate, and highly parallelizable schemes that can easily handle complicated geometries and boundary conditions. We review the theoretical and algorithmic aspects of these methods and show several applications including nonlinear conservation laws, the compressible and incompressible Navier-Stokes equations, and Hamilton-Jacobi-like equations.

Shu, Chi-Wang↗

Compressible 2D Flow Field Interaction of Two Contra-Rotating Blades

The blades of coaxial, contra-rotating rotor systems cross each other in close proximity and at high relative speeds. This crossing event is a potential source of noise and severe blade loads. Effects of compressibility can aggravate the interaction and significantly alter the pressure field signature and phase relationships. A 2-D analysis of this phenomenon is performed by simulating two airfoils passing each other at specified speeds and vertical separation distances. Several test cases spanning a relevant range of Reynolds numbers, angles of attack, and relative Mach number are considered. The Mach number is varied to simulate the radial variation of velocity from the root to tip of a rotor blade to capture the pressure signature, lift, and drag of the airfoils. The velocity and pressure distributions on the airfoils, and in the space between the airfoils are computed before, at, and after airfoil crossing. The variations of lift and drag coefficients through the interaction are captured. The upper airfoil experiences an increase in lift followed by a very sharp drop in lift during the interaction. When relative Mach numbers are transonic, the region of interaction is greatly extended, with shock interactions occurring. The results show the complex nature of the aerodynamic and fluid dynamic impulses generated by blade-blade interactions, with implications to aeroelastic loads and aeroacoustic sources.

Barbely, Natasha L.↗

Improved Models for Simulation of Cryogenic Spray Cooling for Propellant Tank Chilldown

This paper presents advances in spray boiling models for cryogenic fluids that exhibit strong variation in thermodynamic properties; it accounts for the interaction of the spray droplets with walls and the complex heat transfer and boiling processes that play a vital role in chilling solid surfaces. The model has been developed by corroborating and validating it with test data in liquid nitrogen obtained at the University of Connecticut. The test data includes detailed visualization and plate temperature measurements for liquid nitrogen sprays impinging on a warm test plate and subsequent transient chilldown of the plate. Simulations of transient plate chilldown in liquid nitrogen were carried out and temperature histories at various locations on the plate were compared with test measurements to demonstrate the validity of the spray boiling models that have been developed for cryogens.

Boiling↗

Lattice Boltzmann Method for Spacecraft Propellant Slosh Simulation

A scalable computational approach to the simulation of propellant tank sloshing dynamics in microgravity is presented. In this work, we use the lattice Boltzmann equation (LBE) to approximate the behavior of two-phase, single-component isothermal flows at very low Bond numbers. Through the use of a non-ideal gas equation of state and a modified multiple relaxation time (MRT) collision operator, the proposed method can simulate thermodynamically consistent phase transitions at temperatures and density ratios consistent with typical spacecraft cryogenic propellants, for example, liquid oxygen. Determination of the tank forces and moments is based upon a novel approach that relies on the global momentum conservation of the closed fluid domain, and a parametric wall wetting model allows tuning of the free surface contact angle. Development of the interface is implicit and no interface tracking approach is required. A numerical example illustrates the method's application to prediction of bulk fluid behavior during a spacecraft ullage settling maneuver.

Orr, Jeb S.↗

Microgravity Diode Laser Spectroscopy Measurements in a Reacting Vortex Ring

The technique of Diode Laser Spectroscopy (DLS) with wavelength modulation is utilized to measure the concentration of methane in reacting vortex rings under microgravity conditions. From the measured concentration of methane, other major species such as water, carbon dioxide, nitrogen, and oxygen can be easily computed under the assumption of equilibrium chemistry with the method of Interactive Temperature with Assumed Chemistry (ITAC). The conserved scalar approach in modelling the coupling between fluid dynamics and combustion is utilized to represent the unknown variables in terms of the mixture fraction and scalar dissipation rate in conjunction with ITAC. Post-processing of the DLS measurements and the method of ITAC used in computing the species concentration are discussed. From the flame luminosity results, the increase in ring circulation appears to increase the fuel consumption rate inside the reacting vortex ring and the flame height for cases with similar fuel volumes. Preliminary results and application of ITAC show some potential capabilities of ITAC in DLS. The measured concentration of methane, and computed concentrations of water and carbon dioxide agree well with available results from numerical simulations.

Chen, Shin-Juh↗

High Performance Parallel Analysis of Coupled Problems for Aircraft Propulsion

In order to predict the dynamic response of a flexible structure in a fluid flow, the equations of motion of the structure and the fluid must be solved simultaneously. In this paper, we present several partitioned procedures for time-integrating this focus coupled problem and discuss their merits in terms of accuracy, stability, heterogeneous computing, I/O transfers, subcycling, and parallel processing. All theoretical results are derived for a one-dimensional piston model problem with a compressible flow, because the complete three-dimensional aeroelastic problem is difficult to analyze mathematically. However, the insight gained from the analysis of the coupled piston problem and the conclusions drawn from its numerical investigation are confirmed with the numerical simulation of the two-dimensional transient aeroelastic response of a flexible panel in a transonic nonlinear Euler flow regime.

Felippa, C. A.↗

Experimentation and Modeling of Jet A Thermal Stability in a Heated Tube

High performance aircraft typically use hydrocarbon fuel to regeneratively cool the airframe and engine components. As the coolant temperatures increase, the fuel may react with dissolved oxygen forming deposits that limit the regenerative cooling system performance. This study investigates the deposition of Jet A using a thermal stability experiment and computational fluid dynamics (CFD) modeling. The experimental portion of this study is performed with a high Reynolds number thermal stability (HiRets) tester in which fuel passes though an electrically heated tube and the fuel outlet temperature is held constant. If the thermal stability temperature of the fuel is exceeded, deposits form and adhere to the inside of the tube creating an insulating layer between the tube and the fuel. The HiRets tester measures the tube outer wall temperatures near the fuel outlet to report the effect of deposition occurring inside the tube. Final deposits are also estimated with a carbon burn off analysis. The CFD model was developed and used to simulate the fluid dynamics, heat transfer, chemistry, and transport of the deposit precursors. The model is calibrated to the experiment temperature results and carbon burn-off deposition results. The model results show that the dominant factor in deposition is the heated wall temperature and that most of the deposits are formed in the laminar sublayer. The models predicted a 7.0E-6 kilograms per square meter-sec deposition rate, which compared well to the carbon burn-off analysis deposition rate of 1.0E-6 kilograms per square meter-sec.

Khodabandeh, Julia W.↗

An Overview of Unsteady Three-Dimensional Navier-Stokes Simulations of Multistage Turbomachinery Flows

Fluid dynamics of turbomachines are complicated by inherently three dimensional structures such as endwall boundary layers, hub corner separation bubbles and tip-leakage flows. In addition, the relative motion between rotors and stators causes unsteady aerodynamic interactions to occur between blade rows. It is necessary to understand the aerodynamics associated with these interactions in order to design turbomachines that are both light and compact as well as reliable and efficient. An unsteady, three-dimensional, thin-layer, Navier-Stokes zonal algorithm is used to investigate the unsteady aerodynamics of multi-stage turbines and compressors. Relative motion between rotors and stators is made possible by the use of systems of patched and overlaid grids. Time-averaged pressures and pressure envelopes have been computed for several two- and three-dimensional single- and multi-stage configurations. Flow visualizations and computed results are in good agreement with experimental data.

Gundy-Burlet, Karen L.↗