Search NASA⌕ Search

SEARCH · Search NASA

Results for “CFD simulation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 721 records · Page 40

Interaction of external conditions with the internal flowfield in liquid rocket engines - A computational study

Computational studies have been conducted to examine the capability of a CFD code by simulating the steady state thrust chamber internal flow. The SSME served as the sample case, and significant parameter profiles are presented and discussed. Performance predictions from TDK, the recommended JANNAF reference computer program, are compared with those from PHOENICS to establish the credibility of its results. The investigation of an overexpanded nozzle flow is particularly addressed since it plays an important role in the area ratio selection of future rocket engines. Experience gained during this uncompleted flow separation study and future steps are outlined.

Trinh, H. P.↗

Study Of Heating Of The Base Region Of A Rocket

Report describes theoretical study of heating in base region of proposed rocket called "NLS 1.5 stage reference vehicle." Study employed approach based on computational fluid dynamics (CFD). Involved numerical simulations of flow field in base region and in main exhaust plume of cluster of six engines with heat shields.

Ascoli, Edward P.↗

Marshall Space Flight Center surface modeling and grid generation applications

The Solid Rocket Motors (SRM) used by NASA to propel the Space Shuttle employ gimballing nozzles as a means for vehicular guidance during launch and ascent. Gimballing a nozzle renders the pressure field of the exhaust gases nonaxisymmetric. This has two effects: (1) it exerts a torque and side load on the nozzle; and (2) the exhaust gases flow circumferentially in the aft-dome region, thermally loading the flexible boot, case-to-nozzle joint, and casing insulation. The use of CFD models to simulate such flows is imperative in order to assess SRM design. The grids for these problems were constructed by obtaining information from drawings and tabulated coordinates. The 2D axisymmetric grids were designed and generated using the EZ-Surf and GEN2D surface and grid generation codes. These 2D grids were solved using codes such as FDNS, GASP, and MINT. These axisymmetric grids were rotated around the center-line to form 3D nongimballed grids. These were then gimballed around the pivot point and the gaps or overlaps resurfaced to obtain the final domains, which contained approximately 366,000 grid points. The 2D solutions were then rotated and manipulated as appropriate for geometry and used as initial guesses in the final solution. The analyses were used in answering questions about flight criteria.

Williams, Robert W.↗

An Adaptive Flow Solver for Air-Borne Vehicles Undergoing Time-Dependent Motions/Deformations

This report describes a concurrent Euler flow solver for flows around complex 3-D bodies. The solver is based on a cell-centered finite volume methodology on 3-D unstructured tetrahedral grids. In this algorithm, spatial discretization for the inviscid convective term is accomplished using an upwind scheme. A localized reconstruction is done for flow variables which is second order accurate. Evolution in time is accomplished using an explicit three-stage Runge-Kutta method which has second order temporal accuracy. This is adapted for concurrent execution using another proven methodology based on concurrent graph abstraction. This solver operates on heterogeneous network architectures. These architectures may include a broad variety of UNIX workstations and PCs running Windows NT, symmetric multiprocessors and distributed-memory multi-computers. The unstructured grid is generated using commercial grid generation tools. The grid is automatically partitioned using a concurrent algorithm based on heat diffusion. This results in memory requirements that are inversely proportional to the number of processors. The solver uses automatic granularity control and resource management techniques both to balance load and communication requirements, and deal with differing memory constraints. These ideas are again based on heat diffusion. Results are subsequently combined for visualization and analysis using commercial CFD tools. Flow simulation results are demonstrated for a constant section wing at subsonic, transonic, and a supersonic case. These results are compared with experimental data and numerical results of other researchers. Performance results are under way for a variety of network topologies.

Singh, Jatinder↗

Summary Report of the Orbital X-34 Wing Static Aeroelastic Study

This report documents the results of a computational study conducted on the Orbital Sciences X-34 vehicle to compute its inviscid aerodynamic characteristics taking into account the wing structural flexibility. This was a joint exercise between LaRC and SDRC of California. SDRC modeled the structural details of the wing, and provided the structural deformation for a given pressure distribution on its surfaces. This study was done for a Mach number of 1.35 and an angle of attack of 9 deg.; the freestream dynamic pressure was assumed to be 607 lb/sq ft. Only the wing and the body were simulated in the CFD computations. Two wing configurations were examined. The first had the elevons in the undeflected position and the second had the elevons deflected 20 deg. up. The results indicated that with elevon undeflected, the wing twists by about 1.5 deg. resulting in a reduction in the angle of attack at the wing tip to by 1.5 deg. The maximum vertical deflection of the wing is about 3.71 inches at the wing tip. For the wing with the undeflected elevons, the effect of this wing deformation is to reduce the normal force coefficient (C(sub N)) by 0.012 and introduce a noise up pitching moment coefficient (C(sub m)) of 0.042.

Prabhn, Ramadas K.↗

Navier-Stokes Solutions of Rotor and Rotor-Body Flows

This paper will review the advances made recently in the Navier-Stokes CFD methods to simulate aerodynamics and aeroacoustics of helicopter rotors and rotor-body flows. Although a complete flowfield simulation of full helicopter is currently not feasible with these methods, impressive gains have been made in analyzing individual components of this complex problem in a very detailed manner. The use of the state-of-the-art numerical algorithms in solution methods, in conjunction with powerful supercomputers, like the Cray-2, have enabled noticeable progress to be made in modeling viscous-inviscid interactions, blade-vortex interactions, tip-vortex: simulation and wake effects, as well as high speed impulsive noise in hover and forward flight for isolated rotor blades. This paper will critically evaluate the presently available Euler and Navier-Stokes methods, both finite-difference and finite volume methods using structured and unstructured grids for helicopter applications for accuracy, suitability, and computational efficiency. The review will also include the recent progress made using overset grids to model rotor-body flows. All the material for this review will be drawn from the published material shown below.

Srinivasan, G. R.↗

Flow field measurements in the cell culture unit

The cell culture unit (CCU) is being designed to support cell growth for long-duration life science experiments on the International Space Station (ISS). The CCU is a perfused loop system that provides a fluid environment for controlled cell growth experiments within cell specimen chambers (CSCs), and is intended to accommodate diverse cell specimen types. Many of the functional requirements depend on the fluid flow field within the CSC (e.g., feeding and gas management). A design goal of the CCU is to match, within experimental limits, all environmental conditions, other than the effects of gravity on the cells, whether the hardware is in microgravity ( micro g), normal Earth gravity, or up to 2g on the ISS centrifuge. In order to achieve this goal, two steps are being taken. The first step is to characterize the environmental conditions of current 1g cell biology experiments being performed in laboratories using ground-based hardware. The second step is to ensure that the design of the CCU allows the fluid flow conditions found in 1g to be replicated from microgravity up to 2g. The techniques that are being used to take these steps include flow visualization, particle image velocimetry (PIV), and computational fluid dynamics (CFD). Flow visualization using the injection of dye has been used to gain a global perspective of the characteristics of the CSC flow field. To characterize laboratory cell culture conditions, PIV is being used to determine the flow field parameters of cell suspension cultures grown in Erlenmeyer flasks on orbital shakers. These measured parameters will be compared to PIV measurements in the CSCs to ensure that the flow field that cells encounter in CSCs is within the bounds determined for typical laboratory experiments. Using CFD, a detailed simulation is being developed to predict the flow field within the CSC for a wide variety of flow conditions, including microgravity environments. Results from all these measurements and analyses of the CSC flow environment are presented and discussed. The final configuration of the CSC employs magnetic stir bars with angled paddles to achieve the necessary flow requirements within the CSC.

Cell Culture↗

Demonstration, Testing and Qualification of a High Temperature, High Speed Magnetic Thrust Bearing

The gas turbine industry has a continued interest in improving engine performance and reducing net operating and maintenance costs. These goals are being realized because of advancements in aeroelasticity, materials, and computational tools such as CFD and engine simulations. These advancements aid in increasing engine thrust-to-weight ratios, specific fuel consumption, pressure ratios, and overall reliability through higher speed, higher temperature, and more efficient engine operation. Currently, rolling element bearing and squeeze film dampers are used to support rotors in gas turbine engines. Present ball bearing configurations are limited in speed (<2 million DN) and temperature (<5OO F) and require both cooling air and an elaborate lubrication system. Also, ball bearings require extensive preventative maintenance in order to assure their safe operation. Since these bearings are at their operational limits, new technologies must be found in order to take advantage of other advances. Magnetic bearings are well suited to operate at extreme temperatures and higher rotational speeds and are a promising solution to the problems that conventional rolling element bearings present. Magnetic bearing technology is being developed worldwide and is considered an enabling technology for new engine designs. Using magnetic bearings, turbine and compressor spools can be radically redesigned to be significantly larger and stiffer with better damping and higher rotational speeds. These advances, a direct result of magnetic bearing technology, will allow significant increases in engine power and efficiency. Also, magnetic bearings allow for real-time, in-situ health monitoring of the system, lower maintenance costs and down time.

DeWitt, Kenneth↗

Modeling of Slosh Dynamics in Cryogenic Propellant Tanks in Microgravity Environments

The slosh dynamics in cryogenic fuel tanks under microgravity is a pressing problem that severely affects the reliability of launching spacecraft. After reaching low Earth orbit, the propellant in a multistage rocket experiences large and cyclic changes in temperature as a result of solar heating. Tank wall heating can induce thermal stratification and propellant boiloff, particularly during slosh-inducing vehicle maneuvers. Precise understanding of the dynamic and thermodynamic effects of propellant slosh caused by these maneuvers is critical to mission performance and success. Computational fluid dynamics (CFD) analysis is used extensively within the space vehicle industry in an attempt to characterize the behavior of liquids in microgravity, yet experimental data to quantify these predictions is very limited and reduces confidence in the analytical predictions. A novel approach designed to produce high-fidelity data for correlation to CFD model predictions is being developed with the assistance of Florida Institute of Technology (FIT) and Sierra Lobo, Inc. With few exceptions, previous work in slosh dynamics was theoretical or treated the mass of fuel as a variable of inertia only; such models did not consider the viscosity, surface tension, or other important fluid effects. The challenges in this research are in the development of instrumentation able to measure the required parameters, the computational ability to quantify the fluid behaviors, and the means to assess both the measurements and predictions. The design of this experiment bridges the understanding of slosh dynamics in microgravity by a comprehensive approach that combines CFD tools, dynamic simulation tools, semianalytical models of the predominant fluid effects, and an experimental framework that includes measurement and characterization of liquid slosh in one-degree-of-freedom (DOF) and two-DOF experiments, and ultimately experiments in a NASA low-gravity aircraft.

Source record↗

A Modification and Analysis of Lagrangian Trajectory Modeling and Granular Dynamics of Lunar Dust Particles

A previously developed mathematical model is amended to more accurately incorporate the effects of lift and drag on single dust particles in order to predict their behavior in the wake of high velocity gas flow. The model utilizes output from a CFD or DSMC simulation of exhaust from a rocket nozzle hot gas jet. An extension of the Saffman equation for lift based on the research of McLaughlin (1991) and Mei (1992) is used, while an equation for the Magnus force modeled after the work of Oesterle (1994) and Tsuji et al (1985) is applied. A relationship for drag utilizing a particle shape factor (phi = 0.8) is taken from the work of Haider and Levenspiel (1989) for application to non-spherical particle dynamics. The drag equation is further adjusted to account for rarefaction and compressibility effects in rarefied and high Mach number flows according to the work of Davies (1945) and Loth (2007) respectively. Simulations using a more accurate model with the correction factor (Epsilon = 0.8 in a 20% particle concentration gas flow) given by Richardson and Zaki (1954) and Rowe (1961) show that particles have lower ejection angles than those that were previously calculated. This is more prevalent in smaller particles, which are shown through velocity and trajectory comparison to be more influenced by the flow of the surrounding gas. It is shown that particles are more affected by minor changes to drag forces than larger adjustments to lift forces, demanding a closer analysis of the shape and behavior of lunar dust particles and the composition of the surrounding gas flow.

Long, Jason M.↗

Emission Spectroscopic Measurements with an Optical Probe in the NASA Ames IHF Arc Jet Facility

An optical probe was designed to measure radiation (from inside the arc heater) incident on a test sample immersed in the arc-heated stream. Currently, only crude estimates are available for this incident radiation. Unlike efforts of the past, where the probe line of sight was inclined to the nozzle centerline, the present development focuses on having the probe line of sight coincide with the nozzle centerline. A fiber-coupled spectrometer was used to measure the spectral distribution of incident radiation in the wavelength range of 225 to 900 nm. The radiation heat flux in this wavelength range was determined by integration of measured emission spectral intensity calibrated to incident irradiance from an integrating sphere. Two arc-heater conditions, corresponding to stream bulk enthalpy levels of 12 and 22 MJ/kg, were investigated in the 13-inch diameter nozzle of the Interaction Heating Facility at NASA Ames Research Center. With the probe placed at a distance of 10 inches from the nozzle exit plane, total radiative heat fluxes were measured to be 3.3 and 8.4 W/sq cm for the 12 and 22 MJ/kg conditions, respectively. About 17% of these radiative fluxes were due to bound-bound radiation from atoms and molecules, while the remaining 83% could be attributed to continua (bound-free and/or free-free). A comparison with spectral simulation based on CFD solutions for the arc-heater flow field and with spectroscopic measurements in the plenum region indicates that more than 95% of the measured radiation is generated in the arc region. The total radiative heat flux from the line radiation could increase by a factor of two through contributions from wavelengths outside the measured range, i.e., from the vacuum ultraviolet (wavelengths less than 225 nm) and the infrared (wavelengths greater than 900 nm). An extrapolation of the continuum radiation to these two wavelength regions was not attempted. In the tested configuration, the measured radiative heat flux accounts for only about 1.4% of the nominal heat flux on a flat face model and therefore is considered negligible. In the 6-inch diameter nozzle, on account of shorter path lengths, the radiation heat flux could be significant. Therefore, future tests in the 6-inch nozzle will have radiometers in addition to the optical probe.

Winter, Michael↗

Performance Impact of Deflagration to Detonation Transition Enhancing Obstacles

A sub-model is developed to account for the drag and heat transfer enhancement resulting from deflagration-to-detonation (DDT) inducing obstacles commonly used in pulse detonation engines (PDE). The sub-model is incorporated as a source term in a time-accurate, quasi-onedimensional, CFD-based PDE simulation. The simulation and sub-model are then validated through comparison with a particular experiment in which limited DDT obstacle parameters were varied. The simulation is then used to examine the relative contributions from drag and heat transfer to the reduced thrust which is observed. It is found that heat transfer is far more significant than aerodynamic drag in this particular experiment.

Paxson, Daniel E.↗

Development of a Robust and Efficient Parallel Solver for Unsteady Turbomachinery Flows

The traditional design and analysis practice for advanced propulsion systems relies heavily on expensive full-scale prototype development and testing. Over the past decade, use of high-fidelity analysis and design tools such as CFD early in the product development cycle has been identified as one way to alleviate testing costs and to develop these devices better, faster and cheaper. In the design of advanced propulsion systems, CFD plays a major role in defining the required performance over the entire flight regime, as well as in testing the sensitivity of the design to the different modes of operation. Increased emphasis is being placed on developing and applying CFD models to simulate the flow field environments and performance of advanced propulsion systems. This necessitates the development of next generation computational tools which can be used effectively and reliably in a design environment. The turbomachinery simulation capability presented here is being developed in a computational tool called Loci-STREAM [1]. It integrates proven numerical methods for generalized grids and state-of-the-art physical models in a novel rule-based programming framework called Loci [2] which allows: (a) seamless integration of multidisciplinary physics in a unified manner, and (b) automatic handling of massively parallel computing. The objective is to be able to routinely simulate problems involving complex geometries requiring large unstructured grids and complex multidisciplinary physics. An immediate application of interest is simulation of unsteady flows in rocket turbopumps, particularly in cryogenic liquid rocket engines. The key components of the overall methodology presented in this paper are the following: (a) high fidelity unsteady simulation capability based on Detached Eddy Simulation (DES) in conjunction with second-order temporal discretization, (b) compliance with Geometric Conservation Law (GCL) in order to maintain conservative property on moving meshes for second-order time-stepping scheme, (c) a novel cloud-of-points interpolation method (based on a fast parallel kd-tree search algorithm) for interfaces between turbomachinery components in relative motion which is demonstrated to be highly scalable, and (d) demonstrated accuracy and parallel scalability on large grids (approx 250 million cells) in full turbomachinery geometries.

West, Jeff↗

Modeling Droplet Heat and Mass Transfer during Spray Bar Pressure Control of the Multipurpose Hydrogen Test Bed (MHTB) Tank in Normal Gravity

A CFD model for simulating pressure control in cryogenic storage tanks through the injection of a subcooled liquid into the ullage is presented and applied to the 1g MHTB spray bar cooling experiments. An Eulerian-Lagrangian approach is utilized to track the spray droplets and capture the interaction between the discrete droplets and the continuous ullage phase. The spray model is coupled with the VOF model by performing particle tracking in the ullage, removing particles from the ullage when they reach the interface, and then adding their contributions to the liquid. A new model for calculating the droplet-ullage heat and mass transfer is developed. In this model, a droplet is allowed to warm up to the saturation temperature corresponding to the ullage vapor pressure, after which it evaporates while remaining at the saturation temperature. The droplet model is validated against the results of the MHTB spray-bar cooling experiments with 50% and 90% tank fill ratios. The predictions of the present T-sat based model are compared with those of a previously developed kinetic-based droplet mass transfer model. The predictions of the two models regarding the evolving tank pressure and temperature distributions, as well as the droplets' trajectories and temperatures, are examined and compared in detail. Finally, the ullage pressure and local vapor and liquid temperature evolutions are validated against the corresponding data provided by the MHTB spray bar mixing experiment.

mass transfer↗

Computational Modeling Approaches to Multiscale Design of Icephobic Surfaces

To aid in the design of surfaces that prevent icing, a model and computational simulation of impact ice formation at the single droplet scale was implemented. The nucleation of a single supercooled droplet impacting on a substrate, in rime ice conditions, was simulated. Open source computational fluid dynamics (CFD) software was used for the simulation. To aid in the design of surfaces that prevent icing, a model of impact ice formation at the single droplet scale was proposed•No existing model simulates simultaneous impact and freezing of a single super-cooled water droplet•For the 10-week project, a low-fidelity feasibility study was the goal.

adhesion↗

RANS Predictions of Inlet/Isolator Unstart at Mach 4

This paper describes several approaches to back pressuring an inlet/isolator model at Mach 4.0 and the differences in unstart predictions using common RANS simulation methods. The CFD codes OVERFLOW and CREATE-AV Kestrel were used for this work. Verification and validation was conducted on a relevant shock-wave-boundary-layer-interaction case with both solvers. Using aftmounted back pressure flaps did not produce a stable pre-combustion shock train. An unstart comparison between CFD codes OVERFLOW and CREATE-AV Kestrel as well as common 2-equation turbulence models was conducted with numerical back pressure. Unstart predictions using a novel louver back pressure mechanism was done with OVERFLOW. There is notable differences in shock train travel and predicted unstart limit due to turbulence modeling, temporal resolution, and use of the 2D simplification.

CFD↗

Multi-d CFD Modeling of a Free-piston Stirling Convertor at NASA Glenn

A high efficiency Stirling Radioisotope Generator (SRG) is being developed for possible use in long duration space science missions. NASA s advanced technology goals for next generation Stirling convertors include increasing the Carnot efficiency and percent of Carnot efficiency. To help achieve these goals, a multidimensional Computational Fluid Dynamics (CFD) code is being developed to numerically model unsteady fluid flow and heat transfer phenomena of the oscillating working gas inside Stirling convertors. Simulations of the Stirling convertors for the SRG will help characterize the thermodynamic losses resulting from fluid flow and heat transfer between the working gas and solid walls. The current CFD simulation represents approximated 2-dimensional convertor geometry. The simulation solves the Navier Stokes equations for an ideal helium gas oscillating at low speeds. The current simulation results are discussed.

Wilson, Scott D.↗

Development, Verification and Use of Gust Modeling in the NASA Computational Fluid Dynamics Code FUN3D

This paper presents the implementation of gust modeling capability in the CFD code FUN3D. The gust capability is verified by computing the response of an airfoil to a sharp edged gust. This result is compared with the theoretical result. The present simulations will be compared with other CFD gust simulations. This paper also serves as a users manual for FUN3D gust analyses using a variety of gust profiles. Finally, the development of an Auto-Regressive Moving-Average (ARMA) reduced order gust model using a gust with a Gaussian profile in the FUN3D code is presented. ARMA simulated results of a sequence of one-minus-cosine gusts is shown to compare well with the same gust profile computed with FUN3D. Proper Orthogonal Decomposition (POD) is combined with the ARMA modeling technique to predict the time varying pressure coefficient increment distribution due to a novel gust profile. The aeroelastic response of a pitch/plunge airfoil to a gust environment is computed with a reduced order model, and compared with a direct simulation of the system in the FUN3D code. The two results are found to agree very well.

Bartels, Robert E.↗