Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “CFD simulations”

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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
NASA NTRS2010

Numerical Simulations of the Boundary Layer Transition Flight Experiment

Computational Fluid Dynamics (CFD) simulations were used to study the possible effects that the Boundary Layer Transition (BLT) Flight Experiments may have on the heating environment of the Space Shuttle during its entry to Earth. To investigate this issue, hypersonic calculations using the Data-Parallel Line Relaxation (DPLR) and Langley Aerothermodynamic Upwind Relaxation (LAURA) CFD codes were computed for a 0.75 tall protuberance at flight conditions of Mach 15 and 18. These initial results showed high surface heating on the BLT trip and the areas surrounding the protuberance. Since the predicted peak heating rates would exceed the thermal limits of the materials selected to construct the BLT trip, many changes to the geometry were attempted in order to reduce the surface heat flux. The following paper describes the various geometry revisions and the resulting heating environments predicted by the CFD codes.

Tang, Chun Y.↗
NASA NTRS2003

Three Dimensional Numerical Simulation of Rocket-based Combined-cycle Engine Response During Mode Transition Events

The GTX program at NASA Glenn Research Center is designed to develop a launch vehicle concept based on rocket-based combined-cycle (RBCC) propulsion. Experimental testing, cycle analysis, and computational fluid dynamics modeling have all demonstrated the viability of the GTX concept, yet significant technical issues and challenges still remain. Our research effort develops a unique capability for dynamic CFD simulation of complete high-speed propulsion devices and focuses this technology toward analysis of the GTX response during critical mode transition events. Our principal attention is focused on Mode 1/Mode 2 operation, in which initial rocket propulsion is transitioned into thermal-throat ramjet propulsion. A critical element of the GTX concept is the use of an Independent Ramjet Stream (IRS) cycle to provide propulsion at Mach numbers less than 3. In the IRS cycle, rocket thrust is initially used for primary power, and the hot rocket plume is used as a flame-holding mechanism for hydrogen fuel injected into the secondary air stream. A critical aspect is the establishment of a thermal throat in the secondary stream through the combination of area reduction effects and combustion-induced heat release. This is a necessity to enable the power-down of the rocket and the eventual shift to ramjet mode. Our focus in this first year of the grant has been in three areas, each progressing directly toward the key initial goal of simulating thermal throat formation during the IRS cycle: CFD algorithm development; simulation of Mode 1 experiments conducted at Glenn's Rig 1 facility; and IRS cycle simulations. The remainder of this report discusses each of these efforts in detail and presents a plan of work for the next year.

Edwards, Jack R.↗
DOE OSTI2024

Development & Experimental Validation of a Generalized Resistance-Capacitance Model for Numerical Simulation of Phase-Change Material Embedded Heat Exchangers

Latent heat thermal energy storage (LHTES) using phase change material (PCM) has attracted increased attention as a viable solution for overcoming the mismatch between energy supply and demand for renewable energy-based systems. PCM-embedded heat exchangers (PCM-HX) have the potential to significantly improve thermal performance due to high storage capacity and low temperature variation during the phase change process. Most models for simulating LHTES heat transfer use Computational Fluid Dynamics (CFD) simulations, which have high computational costs resulting from considering the complex and time-dependent physics relevant to PCM-HXs. In this paper, a Generalized Resistance Capacitance-based Model (GRCM) was developed to predict the thermal performance of arbitrary PCM-HXs in a computationally efficient manner without compromising modeling accuracy. The GRCM is exercised for three case studies: (i) verification for a single-slabbed finned PCM-HX, (ii) verification and validation for a copper foam/paraffin composite PCM-HX, and (iii) validation for a straight tube annular finned PCM-HX. The copper foam PCM-HX uses an electric heater at the top of HX, while the other two configurations utilize water as heat transfer fluid. For the single-slabbed finned PCM-HX melting case, the mean deviation in average PCM temperature predicted by the GRCM compared to the CFD model was between 0.56 – 0.73 K, with maximum temperature deviation of 2.68 K. For the HTF outlet temperature, the validation results showed that GRCM prediction matches very well with experimental data, with mean temperature deviation of 0.24 K during melting case, while for solidification case was 0.34 K. These results showcase the GRCM’s capability for accurately reproducing the thermal characteristics of PCM-HXs with considerably lower computational effort.

42 ENGINEERING↗
NASA NTRS2004

Development of the Planar Inlet Design and Analysis Process (PINDAP)

The aerodynamic development of an engine inlet requires a comprehensive program of both wind tunnel testing and Computational Fluid Dynamics (CFD) simulations. To save time and resources, much "testing" is done using CFD before any design ever enters a wind tunnel. The focus of my project this summer is on CFD analysis tool development. In particular, I am working to further develop the capabilities of the Planar Inlet Design and Analysis Process (PINDAP). "PINDAP" is a collection of computational tools that allow for efficient and accurate design and analysis of the aerodynamics about and through inlets that can make use of a planar (two-dimensional or axisymmetric) geometric and flow assumption. PINDAP utilizes the WIND CFD flow solver, which is capable of simulating the turbulent, compressible flow field. My project this summer is a continuation of work that I performed for two previous summers. Two years ago, I used basic features of the PINDAP to design a Mach 5 hypersonic scramjet engine inlet and to demonstrate the feasibility of the PINDAP. The following summer, I worked to develop its geometry and grid generation capabilities to include subsonic and supersonic inlets, complete bodies and cowls, conic leading and trailing edges, as well as airfoils. These additions allowed for much more design flexibility when using the program.

Gruber, Christopher R.↗
NASA NTRS2003

Modified Laser Flash Method for Thermal Properties Measurements and the Influence of Heat Convection

The study examined the effect of natural convection in applying the modified laser flash method to measure thermal properties of semiconductor melts. Common laser flash method uses a laser pulse to heat one side of a thin circular sample and measures the temperature response of the other side. Thermal diffusivity can be calculations based on a heat conduction analysis. For semiconductor melt, the sample is contained in a specially designed quartz cell with optical windows on both sides. When laser heats the vertical melt surface, the resulting natural convection can introduce errors in calculation based on heat conduction model alone. The effect of natural convection was studied by CFD simulations with experimental verification by temperature measurement. The CFD results indicated that natural convection would decrease the time needed for the rear side to reach its peak temperature, and also decrease the peak temperature slightly in our experimental configuration. Using the experimental data, the calculation using only heat conduction model resulted in a thermal diffusivity value is about 7.7% lower than that from the model with natural convection. Specific heat capacity was about the same, and the difference is within 1.6%, regardless of heat transfer models.

Lin, Bochuan↗
NASA NTRS2015

On the Comparison of the Long Penetration Mode (LPM) Supersonic Counterflowing Jet to the Supersonic Screech Jet

The authors provide a brief overview of the classic tonal screech noise problem created by underexpanded supersonic jets, briefly describing the fluid dynamic-acoustics feedback mechanism that has been long established as the basis for this well-known aeroacoustics problem. This is followed by a description of the Long Penetration Mode (LPM) supersonic underexpanded counterflowing jet phenomenon which has been demonstrated in several wind tunnel tests and modeled in several computational fluid dynamics (CFD) simulations. The authors provide evidence from test and CFD analysis of LPM that indicates that acoustics feedback and fluid interaction seen in LPM are analogous to the aeroacoustics interactions seen in screech jets. Finally, the authors propose applying certain methodologies to LPM which have been developed and successfully demonstrated in the study of screech jets and mechanically induced excitation in fluid oscillators for decades. The authors conclude that the large body of work done on jet screech, other aeroacoustic phenomena, and fluid oscillators can have direct application to the study and applications of LPM counterflowing supersonic cold flow jets.

Farr, Rebecca A.↗
DOE OSTI2023

Mixing in Low Reynolds Number Reacting Impinging Jets in Crossflow

Previous efforts to model uranyl fluoride formation in an impinging jet gas reactor underpredicted spatial mixing and overpredicted chemical conversion into particulates. The previous fluid dynamics model was based on the solution of the Reynolds Averaged Navier Stokes equations. After simulating fluid dynamics, aerosol dynamics were superimposed onto CFD-simulated gas reactant species concentrations. The current work explores the influence of complex unsteady flow features on the overall flow physics and chemistry for a low Reynolds number, opposed flow, impinging jet gas reactor where there is a low Reynolds number cross flow. The objective of this study was to assess the impact of model formulation on scalar mixing and transport. Here, transient flow simulations were performed using Scale Resolving Simulations. Large-Eddy Simulations with the dynamic Smagorinsky turbulence model were performed along with simulations which directly resolved the flow. Average and root-mean-square (RMS) velocities and species concentrations were computed along with modeled and resolved turbulence kinetic energy (TKE), modeled turbulence dissipation, and modeled turbulent viscosity. Lagrangian flow tracers were also used to quantify species concentrations along path lines emanating from the jet tips. Transient simulation data were compared to results from RANS simulations using the k-ω shear stress transport (SST) model and Reynolds Stress Model (RSM). Transient simulations showed spatial mixing patterns which were more consistent with experimental data and helped elucidate the process of particle formation observed in experiments.

42 ENGINEERING↗
NASA NTRS1992

A status of the activities of the NASA. Marshall Space Flight Center Combustion Devices Technology Team

The Consortium for Computational Fluid Dynamics (CFD) Applications in Propulsion Technology was established to focus on computational fluid dynamics applications in propulsion. Specific areas of effort include developing the CFD technology required to address rocket propulsion issues, validating the technology, and applying the validated technology to design problems. The Combustion Devices Technology Team was formed to implement the above objectives in the broad area of combustion driven flows. In an effort to bring CFD to bear in the design environment, the team has focused its efforts on the Space Transportation Main Engine nozzle. The main emphasis has been on the film cooling scheme used to cool the nozzle wall. Benchmark problems have been chosen to validate CFD film cooling capabilities. CFD simulations of the subscale nozzle have been made. Also, CFD predictions of the base flow resulting from this type of nozzle have been made. The status of these calculations is presented along with future plans. Information is given in viewgraph form.

Tucker, Kevin↗
DOE OSTI2022

Computational Methods to Characterize Panel Loading Conditions for Accelerated Testing

Panel cracking and degradation due to wind loading are known to have a detrimental effect on power output. Recreating these damaging conditions in a controlled experimental setting requires an understanding of the loads being generated at the panel surface as a function of both wind speed, wind direction, and panel orientation. To better understand these relationships, a computational fluid dynamics (CFD) simulation package was constructed using an open-source Python library for solving partial differential equations with the finite element method. This CFD package allows the simulated wind speed and panel orientation of a multi-panel array to be easily changed and provides the capability to measure the traction forces at discrete points along the sun-facing and ground-facing surface of an interior panel residing in the wake of one or more upstream panels. A parameter exploration was performed in which the panel angle was varied in even increments from -70 deg to 70 deg and the wind speed was varied from 2 to 30 m/s. During post processing, the measured traction along the panel surface was averaged spatially and interpreted as a time-varying signal, where further processing of this signal yielded the root mean square amplitude and a characteristic frequency associated with the loading. This study found that higher wind speeds are generally associated with increased amplitude of loading and that the panel orientation angle can significantly exacerbate or mitigate this loading. These outcomes are presented along with current work on higher-fidelity verification simulations and recommendations for performing accelerated experimental testing.

loading↗
NASA NTRS2014

Computational Simulations of Convergent Nozzles for the AIAA 1st Propulsion Aerodynamics Workshop

Computational Fluid Dynamics (CFD) simulations were completed for a series of convergent nozzles in participation of the American Institute of Aeronautics and Astronautics (AIAA) 1st Propulsion Aerodynamics Workshop. The simulations were performed using the Wind-US flow solver. Discharge and thrust coefficients were computed for four axisymmetric nozzles with nozzle pressure ratios (NPR) ranging from 1.4 to 7.0. The computed discharge coefficients showed excellent agreement with available experimental data; the computed thrust coefficients captured trends observed in the experimental data, but over-predicted the thrust coefficient by 0.25 to 1.0 percent. Sonic lines were computed for cases with NPR >= 2.0 and agreed well with experimental data for NPR >= 2.5. Simulations were also performed for a 25 deg. conic nozzle bifurcated by a flat plate at NPR = 4.0. The jet plume shock structure was compared with and without the splitter plate to the experimental data. The Wind-US simulations predicted the shock structure well, though lack of grid resolution in the plume reduced the sharpness of the shock waves. Unsteady Reynolds-Averaged Navier-Stokes (URANS) simulations and Detached Eddy Simulations (DES) were performed at NPR = 1.6 for the 25 deg conic nozzle with splitter plate. The simulations predicted vortex shedding from the trailing edge of the splitter plate. However, the vortices of URANS and DES solutions appeared to dissipate earlier than observed experimentally. It is believed that a lack of grid resolution in the region of the vortex shedding may have caused the vortices to break down too soon

DES↗
NASA NTRS2019

Isolated Rotor Forward Flight Testing at Martian Atmospheric Density Data Report

With the recent interest in Martian exploration using Unmanned Aerial Vehicles (UAVs), an experimental study was conducted to investigate rotor performance at Martian atmospheric conditions. Both simulation and testing of rotors is vital for evaluating rotor performance and behavior, especially for operations at Martian atmospheric densities and pressures. One critical test that has not been performed to date is helicopter forward flight testing at Martian atmospheric densities. To achieve this, a test was conducted in a tunnel facility that could be evacuated to the atmospheric pressure and density of Mars. A 40-inch-diameter rotor, roughly approximating the scale of the proposed Mars Helicopter (MH) design by the NASA Jet Propulsion Laboratory (JPL), was tested in forward flight at Mars’ atmospheric pressure at the NASA Ames Planetary Aeolian Laboratory (PAL). In this forward flight testing, the drive system of the Martian Surface Wind Tunnel (MARSWIT) was never turned on; all wind speeds read were generated either through the rotor spinning or from facility effects. The goal of this experiment was to collect rotor thrust, rotational speed, power, torque, and acoustics measurements. Subsequently, these results can be used for correlation with simulated cases using a mid-fidelity Computational Fluid Dynamics (CFD) simulation. As expected, rotor thrust and power are drastically reduced at air densities 100 times lower than at sea level on Earth. In addition, Reynolds number effects seem to play a vital role at reduced pressure that cannot be neglected in the simulation.

One Atmosphere↗
NASA NTRS2019

Experimental Forward Flight Rotor Performance Testing from Terrestrial to Martian Atmospheric Densities

With the recent interest in Martian exploration using Unmanned Aerial Vehicles (UAV), an experimental study was conducted to investigate rotor performance at Martian atmospheric conditions. Both simulation and testing of rotors is vital for the evaluation of performance and behavior of rotor, especially for operations at Martian atmospheric densities and pressures. Testing and measuring rotor forward flight performance at Martian atmospheric conditions is a relatively unexplored area. Therefore, an experimental study was performed in a wind tunnel to investigate helicopter forward flight performance and to demonstrate successful rotor operation at Martian atmospheric densities. This work was a continuation of the first ever wind tunnel test of a simulated rotorcraft in forward flight at Martian atmospheric densities. A test was conducted in a facility, which could be evacuated to the atmospheric pressure and density of Mars. A 40-in diameter rotor, roughly approximating the scale of the proposed Mars Helicopter design by the NASA Jet Propulsion Laboratory (JPL), was tested in forward flight at Mars atmospheric pressure at the NASA Ames Planetary Aeolian Laboratory (PAL). In this forward flight testing, the drive system of the Martian Surface Wind Tunnel (MARSWIT) was never turned on. The goal of this experiment was to collect rotor thrust, rotational speed, power, torque, and wind speed measurements. Subsequently, these results can be used for correlation with simulated cases using a mid-fidelity Computational Fluid Dynamics (CFD) simulation. Rotor thrust and power seem to decrease approximately proportional to the decrease in density. However, the Reynolds number has an effect on rotor performance that might also be contributing to the change in thrust and power. This effect plays a vital role in rotor performance at reduced pressure that cannot be neglected in the simulation. Despite the challenges involved in testing at a large difference of atmospheric densities between Earth and Mars, repeatable data is obtained in all the measurements at Martian atmospheric conditions.

Martian↗
NASA NTRS2020

Nitric Oxide Laser-Induced Fluorescence Rotational Thermometry in a Hypersonic Non-Equilibrium Flow

Spatially-resolved nitric oxide (NO) planar laser-induced fluorescence (PLIF) rotational thermometry was performed on a Mach 5 non-equilibrium flow around a blunt-body specimen in the Hypersonic Materials Environmental Test System (HyMETS) arc-heat wind tunnel at the NASA Langley Research Center. Transitions within the (0,0) band in the A(exp⁡ 2)∑(sup⁡ +) - X(exp⁡ 2)Π system of NO were excited with a 10 Hz pulsed ultra-violet laser. A multi-line fitting algorithm was applied to the PLIF spectra and transitions were filtered based on an iterative method that maximizes thermometry accuracy. NO PLIF thermometry results are compared to non-equilibrium computational fluid dynamics (CFD) simulations reported in the literature. Differences between the experiment and CFD are quantified for the free-stream, post-shock, shear-layer, and expansion-fan regions of the flow. Potential sources of error, applicable to both the CFD and experiment, are discussed.

Connor C McDougall↗
DOE OSTI2023

An Integrated Approach to Predicting Ash Deposition and Heat Transfer in Coal-Fired Boilers

The overall goal of this project is to develop via measurements and simulations an advanced online technology to predict, monitor and manage fireside ash deposition in a coal-fired boiler allowing for more efficient operations under a range of load conditions and fuel property variability. With this in place fuel sorting and blending can be done upstream and operations can be optimized to compensate for load and fuel properties. In support of this objective, three experimental campaigns were undertaken during the course of the project to measure ash deposition rates within the boiler at different fuel flow rates and its ash composition variability. Simulations of the experimental conditions representing actual geometry, operational scenarios in terms of air flow rates, coal flow rates as well as coal compositions, heating values, and particle size distributions were also carried out. Deposition rates were predicted using a unique particle kinetic energy and viscosity based ash deposition methodology whose validity was ascertained by comparing against deposition rate measurements for widely varying operating conditions and ash compositions in a lab-scale furnace. With a unique end-to-end combustion modeling methodology established and different simulation scenarios carried out, the results from our computational fluid dynamic (CFD) simulations in conjunction with the plant data summarized in this report were used to refine Microbeam Technology Incorporated’s MTI CSPI-CT Tool to predict and monitor fire-side ash deposition under a range of load conditions and fuel property variability in real time.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2022

Design of a gas-solid-solid separator to remove ash from circulating fluidized bed reactors

Cyclones are one of the most common types of gas-solid separators used in circulating fluidized bed boilers. However, cyclones typically do not allow ash to leave the system through the cyclone exit, causing ash to build up in the fluidized bed and necessitating additional systems to remove ash that builds up in the bed. In this study, an alternative “disengager” gas-solid separator is proposed as a way of inherently separating small and large solids, resulting in a gas-solid-solid separation system where ash is allowed to leave the system along with gas while the desired fluidized bed material is retained. Unlike cyclones, which rely on centrifugal force to separate solids and gas, the disengager separates based on entrainment velocity of the particles. Upwards-flowing gas and particles strike a deflection plate and enter the disengaging chamber where particles with low terminal velocity such as ash fines flow with the gas, while larger particles such as sand fall to the bottom of the separator and are returned to the fluidized bed. In this study, several different proposed disengager configurations are simulated and compared to a typical cyclone using computational fluid dynamic (CFD) simulations. Furthermore, it was found that separation efficiency in the disengager is strongly influenced by the size of the deflection plate, rather than by the size of the unit itself. The predicted separation efficiency showed that compared to a cyclone, the disengager design allows significantly more ash to exit the system but retains a similar amount of desirable material. Additionally, the disengager was predicted to not suffer significantly more erosion that a cyclone.

42 ENGINEERING↗
NASA NTRS2017

Progress of High Efficiency Centrifugal Compressor Simulations Using TURBO

Three-dimensional, time-accurate, and phase-lagged computational fluid dynamics (CFD) simulations of the High Efficiency Centrifugal Compressor (HECC) stage were generated using the TURBO solver. Changes to the TURBO Parallel Version 4 source code were made in order to properly model the no-slip boundary condition along the spinning hub region for centrifugal impellers. A startup procedure was developed to generate a converged flow field in TURBO. This procedure initialized computations on a coarsened mesh generated by the Turbomachinery Gridding System (TGS) and relied on a method of systematically increasing wheel speed and backpressure. Baseline design-speed TURBO results generally overpredicted total pressure ratio, adiabatic efficiency, and the choking flow rate of the HECC stage as compared with the design-intent CFD results of Code Leo. Including diffuser fillet geometry in the TURBO computation resulted in a 0.6 percent reduction in the choking flow rate and led to a better match with design-intent CFD. Diffuser fillets reduced annulus cross-sectional area but also reduced corner separation, and thus blockage, in the diffuser passage. It was found that the TURBO computations are somewhat insensitive to inlet total pressure changing from the TURBO default inlet pressure of 14.7 pounds per square inch (101.35 kilopascals) down to 11.0 pounds per square inch (75.83 kilopascals), the inlet pressure of the component test. Off-design tip clearance was modeled in TURBO in two computations: one in which the blade tip geometry was trimmed by 12 mils (0.3048 millimeters), and another in which the hub flow path was moved to reflect a 12-mil axial shift in the impeller hub, creating a step at the hub. The one-dimensional results of these two computations indicate non-negligible differences between the two modeling approaches.

turbomachinery↗
NASA NTRS2011

Aerodynamic Testing of the Orion Launch Abort Tower Separation with Jettison Motor Jet Interactions

The aerodynamic database for the Orion Launch Abort System (LAS) was developed largely from wind tunnel tests involving powered jet simulations of the rocket exhaust plumes, supported by computational fluid dynamics (CFD) simulations. The LAS contains three solid rocket motors used in various phases of an abort to provide propulsion, steering, and Launch Abort Tower (LAT) jettison from the Crew Module (CM). This paper describes a pair of wind tunnel experiments performed at transonic and supersonic speeds to determine the aerodynamic effects due to proximity and jet interactions during LAT jettison from the CM at the end of an abort. The tests were run using two different scale models at angles of attack from 150deg to 200deg , sideslip angles from -10deg to +10deg , and a range of powered thrust levels from the jettison motors to match various jet simulation parameters with flight values. Separation movements between the CM and LAT included axial and vertical translations as well as relative pitch angle between the two bodies. The paper details aspects of the model design, nozzle scaling methodology, instrumentation, testing procedures, and data reduction. Sample data are shown to highlight trends seen in the results.

Rhode, Matthew N.↗
NASA NTRS2006

Estimating Resolution Lengths of Hybrid Turbulence Models

A two-stage procedure has been devised for estimating the spatial resolution achievable in the simulation of a given flow on a given computational grid by a computational fluid dynamics (CFD) code that incorporates a hybrid model of turbulence. The hybrid models to which this procedure is especially relevant are those of the Reynolds-averaged Navier-Stokes (RANS) and the partial-averaged Navier-Stokes (PANS) approaches. This procedure represents the first step toward adding variable-resolution turbulence-modeling capabilities to CFD codes as part of a continuing effort to increase the accuracy and robustness of CFD simulations of unsteady flows. Some background information is prerequisite to a meaningful summary of the procedure. Among experts in CFD, it is well known that combination of the Reynolds-averaged Navier-Stokes (RANS) approach and eddy-viscosity turbulence models offers limited capability for simulating unsteady and complex flows. The RANS approach includes an assumption that most of the energy in a given flow is modeled through turbulence-transport equations and is resolved in a computational grid used to simulate the flow. RANS also overpredicts eddy viscosity, thereby yielding excessive damping of unsteady motion. The eddy viscosity attains an unphysically large value because of unresolved scales, and suppresses most temporal and spatial fluctuations in the resolved flow field. One approach used to overcome this deficiency is to provide a mechanism for the RANS equations to resolve motion only on the largest scales and to use a hybrid model to represent effects at smaller scales. The RANS approach involves the use of a standard two-equation turbulence model in which the effect of turbulence is summarized by a viscosity that is a function of (1) the time-averaged kinetic- energy density (k) associated with the local fluctuating (turbulent) component of flow and (2) the time-averaged rate of dissipation of the turbulent-kinetic- energy density ( ). In PANS, which was developed to overcome the grid dependency associated with other hybrid turbulence models (including that of RANS), the standard two-equation turbulence model is replaced by another two-equation model in which one solves for the previously unresolved k and , which are now allowed to vary in space and time.

Abdol-Hamid, Khaled S.↗