Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Computational Fluid Dynamics Conduction”

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 505 records · Page 28

Nonlinear Dynamic Modeling of a Supersonic Commercial Transport Turbo-Machinery Propulsion System for Aero-Propulso-Servo-Elasticity Research

This paper covers the development of an integrated nonlinear dynamic model for a variable cycle turbofan engine, supersonic inlet, and convergent-divergent nozzle that can be integrated with an aeroelastic vehicle model to create an overall Aero-Propulso-Servo-Elastic (APSE) modeling tool. The primary focus of this study is to provide a means to capture relevant thrust dynamics of a full supersonic propulsion system by using relatively simple quasi-one dimensional computational fluid dynamics (CFD) methods that will allow for accurate control algorithm development and capture the key aspects of the thrust to feed into an APSE model. Previously, propulsion system component models have been developed and are used for this study of the fully integrated propulsion system. An overview of the methodology is presented for the modeling of each propulsion component, with a focus on its associated coupling for the overall model. To conduct APSE studies the described dynamic propulsion system model is integrated into a high fidelity CFD model of the full vehicle capable of conducting aero-elastic studies. Dynamic thrust analysis for the quasi-one dimensional dynamic propulsion system model is presented along with an initial three dimensional flow field model of the engine integrated into a supersonic commercial transport.

Control Theory↗

Unsteady Navier-Stokes computations on a wing-body configuration in ramp motions

Unsteady Navier-Stokes computations are conducted for transonic flows over a wing-body configuration undergoing prescribed ramp motions. The ramp motion from 0 to 15 deg includes angles of attack where vortex breakdown is observed experimentally. The vortex breakdown is found to be delayed until after the ramp motion ends. The dynamic effects on the loads are also demonstrated. The moment coefficient is found to be sensitive to the effect of the virtual mass of the fluid. To verify the numerical results, a grid refinement study has been performed for both steady and unsteady flow conditions, by using up to approximately one million grid points. The results show that the coarse grid can give integrated quantities reasonably well, whereas the finer grids give a more detailed flow structure. Comparisons are also made with available steady-state experimental data.

Obayashi, Shigeru↗

Advancing the State-of-the-Practice for Liquid Rocket Engine Injector Design

Current shortcomings in both the overall injector design process and its underlying combustion stability assessment methodology are rooted in the use of empirically based or low fidelity representations of complex physical phenomena and geometry details that have first order effects on performance, thermal environments and combustion stability. The result is a design and analysis capability that is often inadequate to reliably arrive at a suitable injector design in an efficient manner. Specifically, combustion instability has been particularly difficult to predict and mitigate. Large hydrocarbon-fueled booster engines have been especially problematic in this regard. Where combustion instability has been a problem, costly and time-consuming redesign efforts have often been an unfortunate consequence. This paper presents an overview of a recently completed effort at NASA Marshall Space Flight Center to advance the state-of-the-practice for liquid rocket engine injector design. Multiple perturbations of a gas-centered swirl coaxial (GCSC) element that burned gaseous oxygen and RP-1 were designed, assessed for combustion stability, and tested. Three designs, one stable, one marginally unstable and one unstable, were used to demonstrate both an enhanced overall injector design process and an improved combustion stability assessment process. High-fidelity results from state-of-the-art computational fluid dynamics CFD simulations were used to substantially augment and improve the injector design methodology. The CFD results were used to inform and guide the overall injector design process. They were also used to upgrade selected empirical or low-dimensional quantities in the ROCket Combustor Interactive Design (ROCCID) stability assessment tool. Hot fire single element injector testing was used to verify both the overall injector designs and the stability assessments. Testing was conducted at the Air Force Research Laboratory and at Purdue University. Companion papers provide details of the overall injector design process, full- and sub-scale testing, ROCCID-based stability assessments and the CFD simulations.

Tucker, P. K.↗

Fluid-Structure Interaction Simulations of the ASPIRE SR01 Supersonic Parachute

High-fidelity computational fluid dynamics (CFD) simulations have so far only played exploratory and supporting roles in the study and qualification of parachutes for planetary entry. The last few years have seen the maturation of coupled computational methods that are capable of modeling the complex fluid-structure interactions between a parachute canopy and the supersonic flow in the wake of an entry vehicle in flight conditions. One of the primary goals of these methods is to predict the peak opening load experienced by the parachute during inflation. The Launch, Ascent, and Vehicle Aerodynamics (LAVA)team is developing efficient, high-fidelity numerical methods to perform such challenging fluid-structure interaction simulations. A loose coupling approach is used to advance the solutions of a Cartesian ghost cell immersed boundary method CFD solver and a finite element computational structural dynamics (CSD) solver in space and time. The coupled solver is employed to simulate the ASPIRE SR01 flight test, where a build-to-print version of the Mars Science Laboratory parachute was inflated in supersonic conditions in the upper terrestrial atmosphere. The simulations conducted in the current paper predict a peak opening load that is within 10% of that from the flight test. Grid convergence with respect to the volume and structural domains is demonstrated, and less than 1% variation in the peak opening load is predicted between all grid resolutions.

ARMD↗

Multi-dimensional modeling of a thermal energy storage canister

The Solar Dynamic Power Module being developed for Space Station Freedom uses a eutectic mixture of LiF-CaF2 phase change material (PCM) contained in toroidal canisters for thermal energy storage. Presented are the results from heat transfer analyses of a PCM containment canister. One and two dimensional finite difference computer models are developed to analyze heat transfer in the canister walls, PCM, void, and heat engine working fluid coolant. The modes of heat transfer considered include conduction in canister walls and solid PCM, conduction and pseudo-free convection in liquid PCM, conduction and radiation across PCM vapor filled void regions, and forced convection in the heat engine working fluid. Void shape, location, growth or shrinkage (due to density difference between the solid and liquid PCM phases) are prescribed based on engineering judgment. The PCM phase change process is analyzed using the enthalpy method. The discussion of the results focuses on how canister thermal performance is affected by free convection in the liquid PCM and void heat transfer. Characterizing these effects is important for interpreting the relationship between ground-based canister performance (in 1-g) and expected on-orbit performance (in micro-g). Void regions accentuate canister hot spots and temperature gradients due to their large thermal resistance. Free convection reduces the extent of PCM superheating and lowers canister temperatures during a portion of the PCM thermal charge period. Surprisingly small differences in canister thermal performance result from operation on the ground and operation on-orbit. This lack of a strong gravity dependency is attributed to the large contribution of container walls in overall canister energy redistribution by conduction.

Kerslake, Thomas W.↗

Thermal Analysis of an In-Space Heat Shield Exposed to a Rocket Plume During Stage Separation

Thermal analyses were conducted to evaluate the thermal response of various individual materials and multi-layer configurations for an in-space heat shield exposed to a rocket plume during stage separation. Frequently used, readily available thermal protection materials consisting of cork, carbon cloth phenolic (CCP), and silica cloth phenolic (SCP) were selected for an initial screening analysis. SCP was chosen as the primary material for further evaluation using an updated thermal environment calculated using computational fluid dynamics (CFD) results. A multi-layer configuration consisting of SCP providing erosion resistance and a Nomex honeycomb providing thermal protection was then evaluated as a solution for more mass-efficient performance. Pyrolysis depth, bondline temperature, back-side temperature, and a thermal factor of safety were used to define the necessary material thickness profile. Adhesives were included in the analyses with considerations for maximum temperature and application method. One dimensional thermal analyses at several stations along the heat shield surface were performed for the initial screening and follow-on analyses. The methodology used standard processes and computer programs applicable to solid rocket motor internal insulation and nozzle thermal analysis such as Chemical Equilibrium and Applications (CEA), Aerotherm Chemical Equilibrium (ACE), Momentum/Energy Integral Technique (MEIT), and Insulation Thermal Response and Ablation Code (ITRAC). The baseline SCP/Nomex honeycomb multi-layer configuration thermal performance predictions satisfied the objectives with an acceptable mass estimate for the design maturity. The EA9673 film adhesive and SCP resin were the best performing bonding agents considered for the SCP/Nomex honeycomb and Nomex honeycomb/structural substrate interfaces. The baseline mass estimate was lowest when the resin was the bonding agent due to the higher maximum service temperature. These thermal analyses evaluated typical internal and external nozzle materials using a calculated thermal environment and heat loads that were lower than a typical rocket motor internal environment but greater than standard external aeroheating.

Andrew T Hiatt↗

Experimental results for a hypersonic nozzle/afterbody flow field

This study was conducted to experimentally characterize the flow field created by the interaction of a single-expansion ramp-nozzle (SERN) flow with a hypersonic external stream. Data were obtained from a generic nozzle/afterbody model in the 3.5 Foot Hypersonic Wind Tunnel at the NASA Ames Research Center, in a cooperative experimental program involving Ames and McDonnell Douglas Aerospace. The model design and test planning were performed in close cooperation with members of the Ames computational fluid dynamics (CFD) team for the National Aerospace Plane (NASP) program. This paper presents experimental results consisting of oil-flow and shadow graph flow-visualization photographs, afterbody surface-pressure distributions, rake boundary-layer measurements, Preston-tube skin-friction measurements, and flow field surveys with five-hole and thermocouple probes. The probe data consist of impact pressure, flow direction, and total temperature profiles in the interaction flow field.

Spaid, Frank W.↗

Hypersonic Nozzle/Afterbody Experiment: Flow Visualization and Boundary Layer Experiments

This study was conducted to experimentally characterize the flow field created by the interaction of a single-expansion-ramp-nozzle (SERN) flow with a hypersonic external stream Data were obtained from a generic nozzle/afterbody model in the 3.5-Foot Hypersonic Wind Tunnel of the NASA Ames Research Center in a cooperative experimental program involving Ames and the McDonnell Douglas Aerospace. The model design and test planning were performed in close cooperation with members of the Ames computational fluid dynamics (CFD) team for the National Aero-Space Plane (NASP) program. This paper presents experimental results consisting of oil-flow and shadowgraph flow-visualization photographs, afterbody surface-pressure distributions, boundary-layer rake measurements, and Preston-tube skin-friction measurements.

Keener, Earl R.↗

Impact of Flow Configurations on Response Time and Data Quality in Real-Time, In-Line Fourier Transform Infrared (FTIR) Monitoring of Viscous Flows

The real-time, in-line monitoring of continuous flow concentrations is widely conducted via infrared (IR) spectrometry by using a flow cell connected to a reactive flow stream. For protective purposes, the IR sensor tip is typically offset from the flow. This offset can cause the formation of a stagnant boundary layer above the sensor, especially when dealing with high-viscosity fluids. As a result, the IR signal response time is often controlled by the slow diffusional exchange of fluid in the boundary layer, as confirmed via 2D computational fluid dynamics (CFD) simulations. We evaluated several flow configuration modifications in a typical IR flow cell in order to identify the changes to the flow dynamics that enable improved response times with minimal changes to the cell configuration: the use of (i) vertical flow, where the standard horizontal flow over the sensor is redirected to contact vertically with the sensor, (ii) a static mixer to create radial flow momentum above the IR sensor, and (iii) horizontal or vertical nozzles to direct the flow toward the IR sensor. The vertical flow configuration did not show any significant improvement over the standard horizontal flow configuration. However, the static mixer, horizontal nozzle, and vertical nozzle configurations all resulted in markedly improved response times and signal quality, albeit at the expense of a higher pressure drop across the flow cell. These results point toward straightforward, user-accessible modifications of in-line IR flow cells that result in significant improvements in signal stability and acquisition times.

42 ENGINEERING↗

Ion properties of the high-latitude middle-atmosphere

Investigations of high-latitude, middle-atmosphere ion properties have been conducted using rocket probe techniques. Electrodynamics studies have addressed: mesospheric dynamics and small-scale structure (the MAC/EPSILON Campaign during October-November 1977); ionization effects associated with highly energetic greater than or equal to 1 MeV) precipitating electrons (the MAC/REP Program during May 1990); and the characteristics of noctilucent cloud (NLC) and polar mesosphere summer echo (PMSE) regions (the NLC-91 Campaign during July-August 1991). Recent developments in Gerdien condenser probe theory, using computational fluid dynamics to model ion collection, have resulted in more accurate evaluation of the ion mobilities and number densities. Measurements of very small ion mobilities are thought to indicate an aerosol ion presence in the upper mesosphere.

Li, C.↗

Multiscale graph neural network autoencoders for interpretable scientific machine learning

The goal of this work is to address two limitations in autoencoder-based models: latent space interpretability and compatibility with unstructured meshes. This is accomplished here with the development of a novel graph neural network (GNN) autoencoding architecture with demonstrations on complex fluid flow applications. To address the first goal of interpretability, the GNN autoencoder achieves reduction in the number nodes in the encoding stage through an adaptive graph reduction procedure. Further, this reduction procedure essentially amounts to flowfieldconditioned node sampling and sensor identification, and produces interpretable latent graph representations tailored to the flowfield reconstruction task in the form of so-called masked fields. These masked fields allow the user to (a) visualize where in physical space a given latent graph is active, and (b) interpret the time-evolution of the latent graph connectivity in accordance with the time-evolution of unsteady flow features (e.g. recirculation zones, shear layers) in the domain. To address the goal of unstructured mesh compatibility, the autoencoding architecture utilizes a series of multi-scale message passing (MMP) layers, each of which models information exchange among node neighborhoods at various lengthscales. The MMP layer, which augments standard single-scale message passing with learnable coarsening operations, allows the decoder to more efficiently reconstruct the flowfield from the identified regions in the masked fields. Analysis of latent graphs produced by the autoencoder for various model settings are conducted using unstructured snapshot data sourced from large-eddy simulations in a backward-facing step (BFS) flow configuration with an OpenFOAM-based flow solver at high Reynolds numbers.

97 MATHEMATICS AND COMPUTING↗

Investigation of thermal hydraulic behavior of the High Temperature Test Facility's lower plenum via large eddy simulation

A high-fidelity computational fluid dynamics (CFD) analysis was performed using the Large Eddy Simulation (LES) model for the lower plenum of the High–Temperature Test Facility (HTTF), a ¼ scale test facility of the modular high temperature gas-cooled reactor (MHTGR) managed by Oregon State University. In most next–generation nuclear reactors, thermal stress due to thermal striping is one of the risks to be curiously considered. This is also true for HTGRs, especially since the exhaust helium gas temperature is high. In order to evaluate these risks and performance, organizations in the United States led by the OECD NEA are conducting a thermal hydraulic code benchmark for HTGR, and the test facility used for this benchmark is HTTF. HTTF can perform experiments in both normal and accident situations and provide high-quality experimental data. However, it is difficult to provide sufficient data for benchmarking through experiments, and there is a problem with the reliability of CFD analysis results based on Reynolds–averaged Navier–Stokes to analyze thermal hydraulic behavior without verification. To solve this problem, high-fidelity 3-D CFD analysis was performed using the LES model for HTTF. It was also verified that the LES model can properly simulate this jet mixing phenomenon via a unit cell test that provides experimental information. As a result of CFD analysis, the lower the dependency of the sub-grid scale model, the closer to the actual analysis result. In the case of unit cell test CFD analysis and HTTF CFD analysis, the volume-averaged sub-grid scale model dependency was calculated to be 13.0% and 9.16%, respectively. As a result of HTTF analysis, quantitative data of the fluid inside the HTTF lower plenum was provided in this paper. As a result of qualitative analysis, the temperature was highest at the center of the lower plenum, while the temperature fluctuation was highest near the edge of the lower plenum wall. The power spectral density of temperature was analyzed via fast Fourier transform (FFT) for specific points on the center and side of the lower plenum. FFT results did not reveal specific frequency-dominant temperature fluctuations in the center part. It was confirmed that the temperature power spectral density (PSD) at the top increased from the center to the wake. The vortex was visualized using the well-known scalar Q-criterion, and as a result, the closer to the outlet duct, the greater the influence of the mainstream, so that the inflow jet vortex was dissipated and mixed at the top of the lower plenum. Additionally, FFT analysis was performed on the support structure near the corner of the lower plenum with large temperature fluctuations, and as a result, it was confirmed that the temperature fluctuation of the flow did not have a significant effect near the corner wall. In addition, the vortices generated from the lower plenum to the outlet duct were identified in this paper. It is considered that the quantitative and qualitative results presented in this paper will serve as reference data for the benchmark.

97 MATHEMATICS AND COMPUTING↗

CFD analysis of pump consortium impeller

Current design of high performance turbopumps for rocket engines requires effective and robust analytical tools to provide design impact in a productive manner. The main goal of this study is to develop a robust and effective computational fluid dynamics (CFD) pump model for general turbopump design and analysis applications. A Navier-Stokes flow solver, FDNS, embedded with the extended k-epsilon turbulence model and with appropriate moving interface boundary conditions, is developed to analyze turbulent flows in the turbomachinery devices. The FDNS code was benchmarked with its numerical predictions of the pump consortium inducer, and provides satisfactory results. In the present study, a CFD analysis of the pump consortium impeller will be conducted with the application of the FDNS code. The pump consortium impeller, with partial blades, is the new design concept of the advanced rocket engine.

Cheng, Gary C.↗

Computational Fluid Dynamics Analysis Method Developed for Rocket-Based Combined Cycle Engine Inlet

Renewed interest in hypersonic propulsion systems has led to research programs investigating combined cycle engines that are designed to operate efficiently across the flight regime. The Rocket-Based Combined Cycle Engine is a propulsion system under development at the NASA Lewis Research Center. This engine integrates a high specific impulse, low thrust-to-weight, airbreathing engine with a low-impulse, high thrust-to-weight rocket. From takeoff to Mach 2.5, the engine operates as an air-augmented rocket. At Mach 2.5, the engine becomes a dual-mode ramjet; and beyond Mach 8, the rocket is turned back on. One Rocket-Based Combined Cycle Engine variation known as the "Strut-Jet" concept is being investigated jointly by NASA Lewis, the U.S. Air Force, Gencorp Aerojet, General Applied Science Labs (GASL), and Lockheed Martin Corporation. Work thus far has included wind tunnel experiments and computational fluid dynamics (CFD) investigations with the NPARC code. The CFD method was initiated by modeling the geometry of the Strut-Jet with the GRIDGEN structured grid generator. Grids representing a subscale inlet model and the full-scale demonstrator geometry were constructed. These grids modeled one-half of the symmetric inlet flow path, including the precompression plate, diverter, center duct, side duct, and combustor. After the grid generation, full Navier-Stokes flow simulations were conducted with the NPARC Navier-Stokes code. The Chien low-Reynolds-number k-e turbulence model was employed to simulate the high-speed turbulent flow. Finally, the CFD solutions were postprocessed with a Fortran code. This code provided wall static pressure distributions, pitot pressure distributions, mass flow rates, and internal drag. These results were compared with experimental data from a subscale inlet test for code validation; then they were used to help evaluate the demonstrator engine net thrust.

Source record↗

MPEX High Heat Flux Plasma Dump Design, Manufacture, and Articles Test

Here, the Material Plasma Exposure eXperiment (MPEX) device is a steady-state linear plasma device currently in the final design phase at the Oak Ridge National Laboratory. This device will reach ion fluences up to 10 31 m -2 and will be used to expose neutron-irradiated materials to divertor-relevant plasmas and to study the effects of plasma-material interactions. These studies will elucidate the complex effects of plasmas with divertor candidate materials capable of withstanding high heat flux and high fluences for next-generation fusion devices. Bidirectional plasma will be generated using a high-power (200 kW) helicon source. Plasma will be confined by superconducting magnets. The last plasma-facing component on the upstream side of the MPEX device is the dump, which has been designed to intercept plasma and energetic particles. The dump will have a total heat load of 9.2 kW. A copper alloy (Glidcop AL-15) was selected for use in the water-cooled flange design because of its high thermal conductivity, its retaining strength at elevated temperatures, and its ability to be used in the high-temperature braze joints used in this application. Titanium-zirconium-molybdenum (TZM) tiles are brazed to the Glidcop AL-15 flange using high-temperature braze alloy. External water-cooling channels are used on the dump flange to prevent water leakage inside the vacuum space. This article discusses the details of the high heat flux dump design, including the computational fluid dynamics (CFD) and structural analyses performed to validate the design to meet the operational requirements of the MPEX device.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Performance Improvement Through Indexing of Turbine Airfoils: Numerical Simulation - Part 2

An experimental/analytical study has been conducted to determine the performance improvements achievable by circumferentially indexing succeeding rows of turbine stator airfoils. A series of tests was conducted to experimentally investigate stator wake clocking effects on the performance of the space shuttle main engine (SSME) alternate turbopump development (ATD) fuel turbine test article (TTA). The results from this study indicate that significant increases in stage efficiency can be attained through application of this airfoil clocking concept. Details of the experiment and its results are documented in part 1 of this paper. In order to gain insight into the mechanisms of the performance improvement, extensive computational fluid dynamics (CFD) simulations were executed. The subject of the present paper is the initial results from the CFD investigation of the configurations and conditions detailed in part 1 of the paper. To characterize the aerodynamic environments in the experimental test series, two-dimensional (2D), time accurate, multistage, viscous analyses were performed at the TTA midspan. Computational analyses for five different circumferential positions of the first stage stator have been completed. Details of the computational procedure and the results are presented. The analytical results verify the experimentally demonstrated performance improvement and are compared with data whenever possible. Predictions of time-averaged turbine efficiencies as well as gas conditions throughout the flow field are presented. An initial understanding of the turbine performance improvement mechanism based on the results from this investigation is described.

Griffin, Lisa W.↗

Air-Loads Prediction of a UH-60A Rotor inside the 40- by 80-Foot Wind Tunnel

The presented research extends the capability of a loose coupling computational fluid dynamics (CFD) and computational structure dynamics (CSD) code to calculate the flow-field around a rotor and test stand mounted inside a wind tunnel. Comparison of predicted air-load results for a full-scale UH-60A rotor recently tested inside the National Full-Scale Aerodynamics Complex (NFAC) 40- by 80-Foot Wind Tunnel at Ames Research Center and in free-air flight are made for three challenging flight data points from the earlier conducted UH-60A Air-loads Program. Overall results show that the extension of the coupled CFD/CSD code to the wind-tunnel environment is generally successful.

Chang, I-Chung↗

The Onset of Nucleate Boiling Prediction in a Mini Channel with the Eulerian Multiphase Flow

The prediction of boiling incipience is a critical issue for a reactor. Vapor lock causes operating instability and undesirable temperature rise, especially for small channels. This study uses the Eulerian multiphase flow boiling model to predict the onset of nucleate boiling in a millimeter scale rectangular channel. Two wall boiling models, namely the RPI (Rensselaer Polytechnic Institute) and the non-equilibrium sub cooled model are applied. Computational fluid dynamics is employed for this study to investigate the impact of operating conditions like liquid flow rate, operating pressure, heat flux, velocity profile, turbulence model, and inlet liquid temperature. Here, the effect of fundamental boiling parameters on the surface wall temperature and vapor volume fraction is also studied. The parameters under consideration are the bubble departure diameter, bubble departure frequency, nucleate site density, quenching time period, and interface heat transfer coefficients. A series of CFD calculations is conducted by varying the considered variables systematically in a wide range of flow conditions covering laminar, transition, and turbulent flows. The sub-component heat fluxes like the liquid convective, the vapor convective, the evaporative, and the quenching terms are monitored to infer the boiling dynamics in the transition region from the single-phase flow to the nucleate boiling zone. It is found that the CFD approach to detecting the boiling incipience point agrees reasonably with available experimental data. However, its limitations, like inaccuracy for lower flow rates and early transition to boiling flow, are also noticed due to the nature of the considered RPI model.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗