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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Coupled Transport, Reactivity, and Mechanics in Fractured Shale Caprocks

Abstract Shales are low‐permeability caprocks that confine fluid, such as CO 2 , nuclear waste, and hydrogen, in storage formations. Stress‐induced fractures in shale caprocks provide pathways for fluid to leak and potentially contaminate fresh water aquifers. Fractured shales are also increasingly considered as resources for CO 2 sequestration, enhanced geothermal, and unconventional energy recovery. Injecting reactive fluids into shales introduces chemical disequilibrium, causing an onset of a series of dissolution, precipitation, and fines mobilization mechanisms. The reactions have rapid kinetics and significant impact on porosity and permeability; consequently, flow and storage properties of caprocks. While previous research has explored the separate effects of these reactions, this study aims to uncover their simultaneous occurrence and collective influence. This study unveils these highly coupled transport and reactivity mechanisms by tracking and visualizing the reaction‐induced alterations in the matrix, microcracks, and fractures of shales over time. We conducted brine injection experiments sequentially at pH 4 and 2 in a naturally fractured Wolfcamp shale sample while simultaneously imaging the dynamic processes using X‐ray computed tomography (CT). CT images are validated by finer resolution images obtained using micro‐CT and scanning electron microscopy. We also tracked the sample permeability and fluid chemistry using brine permeability and inductively coupled plasma mass spectrometry, respectively. Findings show that fluid primarily flowed through fractures, dissolving reactive minerals and mobilizing fines on fracture surfaces. Dissolution of fracture asperities under confining stress resulted in the closing of fractures. Clogging in narrow fracture pathways, caused by fines accumulation, diverted fluid flow into matrix pores.

58 GEOSCIENCES↗
DOE OSTI2020

GaAs growth rates of 528 μm/h using dynamic-hydride vapor phase epitaxy with a nitrogen carrier gas

Herein we demonstrate record GaAs growth rates approaching 530 μm/h using nitrogen carrier gas and 400 μm/h using hydrogen carrier gas in a dynamic-hydride vapor phase epitaxy reactor. We measured root mean square surface roughness below 1 nm using a 1 μm × 1 μm atomic force microscopy scan for GaAs growth rates up to 483 μm/h using a nitrogen carrier gas and 400 μm/h using a hydrogen carrier gas. We performed computational fluid dynamics modeling to study the effect of the carrier gas choice on the thermal profile within the reactor and how that influences the degree of AsH 3 decomposition. The modeling suggests that the lower thermal conductivity of the nitrogen carrier gas minimizes the amount of AsH 3 that thermally decomposes before reaching the wafer surface and the heavier atomic mass decreases the likelihood that AsH 3 will reach the heated reactor walls, leading to a growth rate enhancement relative to the hydrogen carrier case in a hydride-enhanced growth regime.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Niowave Neutron Source Converter: Lead-Bismuth-Eutectic (LBE) Windowless Target Design and Evaluation

Los Alamos National Laboratory (LANL) is working with Niowave on the design and evaluation of their lead-bismuth-eutectic (LBE) windowless target (i.e., neutron source converter). Niowave plans to use 200 kW electron beam at 40 MeV beam energy to produce neutrons by photonuclear reaction with LBE. Then, the neutrons undergo fission at the surrounding uranium target assembly (UTA) to produce Molybdenum 99 (Mo-99) as a fission product, which eventually decays to Technetium-99 (Tc-99m). Tc-99m is one of important radioisotopes that is used for medical diagnostics. LANL conducted 3D multiphysics analysis for the Niowave neutron converter design and provided design assessment in thermohydraulic aspects. LANL conducted radiation transport calculations using Monte-Carlo N-Particle (MCNP) code with unstructured meshing scheme. The 3D volumetric heating profiles in the LBE and Stainless-Steel (SS) housing were imported into multiphase computational fluid dynamics (CFD) to obtain 3D temperature profiles of LBE and SS through conjugate heat transfer (CHT) analysis. The key findings are: LBE film thickness at the center of the beam is approximately 1.6 cm with a maximum LBE velocity of approximately 1.8 m/s, which is below a 2 m/s limit to avoid erosion issues on supporting structures; Heat deposition in the LBE peaks at ~1 cm depth from the LBE free-surface because of the forward interactions of electron, photon, and neutron with LBE; LBE maximum temperature is ~360 °C, which is below LBE evaporation initiative temperature, ~450 °C; LBE-SS interface temperature is ~350 °C, which is below the safety thermal limit to prevent severe corrosion on SS. The results indicate that Niowave’s neutron converter design satisfies both hydraulic and thermal criteria for safe operation. By virtue of such computational analysis, Niowave can move toward establishing an experimental setup to experimentally test their LBE neutron converter. The following sections describe the detailed work done by LANL.

43 PARTICLE ACCELERATORS↗
NASA NTRS1996

Wetting Hysteresis at the Molecular Scale

The motion of a fluid-fluid-solid contact line on a rough surface is well known to display hysteresis in the contact angle vs. velocity relationship. In order to understand the phenomenon at a fundamental microscopic level, we have conducted molecular dynamics computer simulations of a Wilhelmy plate experiment in which a solid surface is dipped into a liquid bath, and the force-velocity characteristics are measured. We directly observe a systematic variation of force and contact angle with velocity, which is single-valued for the case of an atomically smooth solid surface. In the microscopically rough case, however, we find (as intuitively expected) an open hysteresis loop. Further characterization of the interface dynamics is in progress.

Jin, Wei↗
NASA NTRS2023

CFD Integrated Transition Modeling for High-Speed Flows via Coupled OVERFLOW-LASTRAC Analysis

This work details ongoing efforts at the NASA Langley Research Center to develop and validate a general-use CFD tool that includes built-in predictions of boundary-layer transition in high-speed flows. Existing tools for the coupling of the NASA OVERFLOW structured overset RANS solver and LASTRAC stability analysis code have been extended to capture boundary-layer transition in high-speed flows driven by either Mack’s first-mode or secondmode instabilities and/or crossflow instabilities. The efficacy of this coupled approach has been demonstrated by examining a variety of supersonic and hypersonic test cases for which experimental validation data is available. The intermittency prescription parameters inherited from prior low-speed applications are found to be suitable for some, but not all high-speed flow scenarios. The method is found to be robust for straight cone configurations and improvements necessary to accurately capture transition on geometries with spatially complex amplification factor envelopes with intermediate regions of slow N-factor variation are examined. The automated, coupled analysis is demonstrated for a 3d supersonic test case and found to perform well within the limits of the linear stability analysis on which it relies. Finally, a preliminary investigation of the method’s robustness to non-ideal CFD meshes is conducted.

High-speed flow↗
NASA NTRS2003

Enhancement and Extension of Porosity Model in the FDNS-500 Code to Provide Enhanced Simulations of Rocket Engine Components

In the past, the design of rocket engines has primarily relied on the cold flow/hot fire test, and the empirical correlations developed based on the database from previous designs. However, it is very costly to fabricate and test various hardware designs during the design cycle, whereas the empirical model becomes unreliable in designing the advanced rocket engine where its operating conditions exceed the range of the database. The main goal of the 2nd Generation Reusable Launching Vehicle (GEN-II RLV) is to reduce the cost per payload and to extend the life of the hardware, which poses a great challenge to the rocket engine design. Hence, understanding the flow characteristics in each engine components is thus critical to the engine design. In the last few decades, the methodology of computational fluid dynamics (CFD) has been advanced to be a mature tool of analyzing various engine components. Therefore, it is important for the CFD design tool to be able to properly simulate the hot flow environment near the liquid injector, and thus to accurately predict the heat load to the injector faceplate. However, to date it is still not feasible to conduct CFD simulations of the detailed flowfield with very complicated geometries such as fluid flow and heat transfer in an injector assembly and through a porous plate, which requires gigantic computer memories and power to resolve the detailed geometry. The rigimesh (a sintered metal material), utilized to reduce the heat load to the faceplate, is one of the design concepts for the injector faceplate of the GEN-II RLV. In addition, the injector assembly is designed to distribute propellants into the combustion chamber of the liquid rocket engine. A porosity mode thus becomes a necessity for the CFD code in order to efficiently simulate the flow and heat transfer in these porous media, and maintain good accuracy in describing the flow fields. Currently, the FDNS (Finite Difference Navier-Stakes) code is one of the CFD codes which are most widely used by research engineers at NASA Marshall Space Flight Center (MSFC) to simulate various flow problems related to rocket engines. The objective of this research work during the 10-week summer faculty fellowship program was to 1) debug the framework of the porosity model in the current FDNS code, and 2) validate the porosity model by simulating flows through various porous media such as tube banks and porous plate.

Cheng, Gary↗
NASA NTRS2007

Computational Analyses of Offset Stream Nozzles for Noise Reduction

The Wind computational fluid dynamics code was used to perform a series of simulations on two offset stream nozzle concepts for jet noise reduction. The first concept used an S-duct to direct the secondary stream to the lower side of the nozzle. The second concept used vanes to turn the secondary flow downward. The analyses were completed in preparation of tests conducted in the NASA Glenn Research Center Aeroacoustic Propulsion Laboratory. The offset stream nozzles demonstrated good performance and reduced the amount of turbulence on the lower side of the jet plume. The computer analyses proved instrumental in guiding the development of the final test configurations and giving insight into the flow mechanics of offset stream nozzles. The computational predictions were compared with flowfield results from the jet rig testing and showed excellent agreement.

Dippold, Vance, III↗
NASA NTRS2016

Space Launch System Booster Separation Aerodynamic Testing in the NASA Langley Unitary Plan Wind Tunnel

A wind-tunnel investigation of a 0.009 scale model of the Space Launch System (SLS) was conducted in the NASA Langley Unitary Plan Wind Tunnel to characterize the aerodynamics of the core and solid rocket boosters (SRBs) during booster separation. High-pressure air was used to simulate plumes from the booster separation motors (BSMs) located on the nose and aft skirt of the SRBs. Force and moment data were acquired on the core and SRBs. These data were used to corroborate computational fluid dynamics (CFD) calculations that were used in developing a booster separation database. The SRBs could be remotely positioned in the x-, y-, and z-direction relative to the core. Data were acquired continuously while the SRBs were moved in the axial direction. The primary parameters varied during the test were: core pitch angle; SRB pitch and yaw angles; SRB nose x-, y-, and z-position relative to the core; and BSM plenum pressure. The test was conducted at a free-stream Mach number of 4.25 and a unit Reynolds number of 1.5 million per foot.

Wilcox, Floyd J., Jr.↗
NASA NTRS2011

Sonic Boom Computations for a Mach 1.6 Cruise Low Boom Configuration and Comparisons with Wind Tunnel Data

Accurate analysis of sonic boom pressure signatures using computational fluid dynamics techniques remains quite challenging. Although CFD shows accurate predictions of flow around complex configurations, generating grids that can resolve the sonic boom signature far away from the body is a challenge. The test case chosen for this study corresponds to an experimental wind-tunnel test that was conducted to measure the sonic boom pressure signature of a low boom configuration designed by Gulfstream Aerospace Corporation. Two widely used NASA codes, USM3D and AERO, are examined for their ability to accurately capture sonic boom signature. Numerical simulations are conducted for a free-stream Mach number of 1.6, angle of attack of 0.3 and Reynolds number of 3.85x10(exp 6) based on model reference length. Flow around the low boom configuration in free air and inside the Langley Unitary plan wind tunnel are computed. Results from the numerical simulations are compared with wind tunnel data. The effects of viscous and turbulence modeling along with tunnel walls on the computed sonic boom signature are presented and discussed.

Elmiligui, Alaa A.↗
NASA NTRS1991

Theoretical and computational dynamics of a compressible flow

An introduction to the theoretical and computational fluid dynamics of a compressible fluid is presented. The general topics addressed include: thermodynamics and physical properties of compressible fluids; 1D flow of an inviscid compressible fluid; shock waves; fundamental equations of the dynamics of a compressible inviscid non-heat-conducting and radiating fluid, method of small perturbations, linearized theory; 2D subsonic steady potential flow; hodograph and rheograph methods, exact solutions of 2D insentropic steady flow equations, 2D steady transonic and hypersonic flows, method of characteristics, linearized theory of 3D potential flow, nonlinear theory of 3D compressibe flow, anisentropic (rotational) flow of inviscid compressible fluid, electromagnetogasdynamics, multiphase flows, flows of a compressible fluid with transport phenomena.

Pai, Shih-I↗
NASA NTRS2011

Continuing Validation of Computational Fluid Dynamics for Supersonic Retropropulsion

A large step in the validation of Computational Fluid Dynamics (CFD) for Supersonic Retropropulsion (SRP) is shown through the comparison of three Navier-Stokes solvers (DPLR, FUN3D, and OVERFLOW) and wind tunnel test results. The test was designed specifically for CFD validation and was conducted in the Langley supersonic 4 x4 Unitary Plan Wind Tunnel and includes variations in the number of nozzles, Mach and Reynolds numbers, thrust coefficient, and angles of orientation. Code-to-code and code-to-test comparisons are encouraging and possible error sources are discussed.

Schauerhamer, Daniel Guy↗
NASA NTRS2014

Afterbody Heat Flux Measurements in the NASA Ames HFFAF Ballistic Range

In order to measure afterbody heat fluxes over a model in the ballistic range, the required modifications to a proven technique for measuring forebody heat fluxes are described. This involves the use of an extended helium gas plume to remove the glowing wake and the use of special high conductivity, high temperature capable graphite-filled plastic for the afterbody. The models and test conditions are described. Data in the form of plots of the surface temperature of the models are presented. Finally, experimental and computational fluid dynamic (CFD) heat flux data for forebody and afterbody heat fluxes are presented and compared. Data are presented for a 45 degree sphere-cone (with a projecting rear stud) at 2.70 km/s and for a sphere at 4.76 km/s. Both models were launched into 76 Torr of CO2 gas. The experimental forebody heat fluxes were within 1.5% of the CFD values. The experimental afterbody heat fluxes were within 1% of the CFD values for the sphere, but only 51% of the CFD values for the sphere-cone.

Bogdanoff, David W.↗
DOE OSTI2024

RELAP5-3D in the OECD-NEA HTGR Thermal Hydraulics Benchmark

The OECD-NEA HTGR Thermal Hydraulics benchmark uses data from the High Temperature Test Facility for a series of code-to-code and code-to-data exercises aimed at improving the state of knowledge on existing thermal hydraulics modeling and simulation tools for prismatic HTGR applications. This benchmark includes problems representing hot gas mixing in the lower plenum, the depressurized conduction cooldown accident and the pressurized conduction cooldown accident. This presentation discusses the benchmark and RELAP5-3D's role in that benchmark, including as a tool for predicting behavior in the core and as a tool for providing boundary conditions to computational fluid dynamics analysis in the lower plenum.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

CFD Simulation of the Dosing Behavior within the Atomic Layer Deposition Feeding System

The effective operation of atomic layer deposition (ALD) feeding system is the premise of realizing specific ALD processes. In the present work, a detailed computational fluid dynamics (CFD) model of the feeding system has been developed and validated, which accounts for the roles of ALD valves and manifolds. A numerical simulation of the compressible fluid flow and heat/mass transfer within the feeding system was conducted. The dosing amounts and the spatiotemporal distributions of the precursors can be accurately predicted using the CFD model, as validated by experimental results. Different precursors, operating conditions, and structures of the feeding system were simulated and analyzed to examine the operating flexibility of the feeding system. The simulation results can be adopted as the upstream boundary conditions for simulations of the ALD process in the reaction chamber. The substrate-scale simulation indicates that the effect of the feeding system on the film deposition is highly related to the surface kinetics of ALD. The present work can serve as a guide for the development and optimization of different ALD-based processes via proper operation and even the design of the feeding system.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Pressure Response Optimization of an Eddy Current-Driven Flyer Plate Valve for the ITER Shattered Pellet Injection System

One technique for mitigating disruptions in a tokamak is shattered pellet injection (SPI). SPI is a process in which a large solid pellet consisting of deuterium, neon, or argon is desublimated in a pipe gun barsrel and launched downstream. Pellets are shattered just before entering the plasma by an impact with an angled tube. Injection of these materials into the plasma radiates stored thermal energy, limits current decay rates, suppresses the generation of runaway electrons, and dissipates runaway electrons if necessary. A critical element of the SPI system is a fast-acting valve that releases high-pressure gas to dislodge and accelerate pellets directly, or indirectly via a mechanical punch. A prototype valve sized for the ITER SPI system has been designed and fabricated. A pulsed high-voltage power supply energizes the valve’s internal magnetic coil, which induces eddy currents in the adjacent flyer plate resulting in a repulsive force between the flyer plate and the coil. The flyer plate action lifts a valve seat, allowing high-pressure gas to flow from the valve plenum to the downstream (breech) location of the pellet or mechanical punch. All of the valve’s internal components are designed to operate in ITER-level static background magnetic fields. Here, a study was conducted to optimize the downstream pressure response for a range of valve sizes and operating pressures. In particular, the study analyzes the breech pressure response associated with varying plenum pressures as well as varying breech volumes. A computational fluid dynamics simulation was built in STAR-CCM+ and validated against data from laboratory experiments. The resulting simulation outputs, in the form of downstream responses for a variety of initial plenum pressures and breech volumes, will be used as a complement to experimental data to ensure the pressure pulse is suitable for pellet survivability. These data, combined with studies on pellet shear strength and shock response, will be applied to optimization of overall operating parameters of the ITER SPI system.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Ring Model Development and Validation for Prismatic HTGR Core Thermal-hydraulics and Safety Analysis

Because of the complex core geometry, prismatic high temperature gas-cooled reactors (prismatic HTGRs) often exhibit complex thermal fluid behaviors during both normal operating and transient conditions. Most HTGR designs rely on passive safety system for decay heat removal, such as the reactor cavity cooling system (RCCS). During postulated accidents like Pressurized Conduction Cooldown (PCC) event, the decay heat is first radially transferred from the core region to the reactor vessel outer surface, then to the RCCS cooling panels. The peak fuel temperature is controlled by heat transfer mechanisms with two distinctive characteristic length scales, i.e., the core-wise effective heat conduction and local heat conduction in the fuel pellet scale. As both length scales are essential to determine the fuel temperature, from the modeling perspective, computer codes must be able to capture heat transfer in both scales. This is challenging for both computational fluid dynamics (CFD) tools because of extremely large amount of computation resources required, and for system analysis codes because of the challenge to model the complex core geometry. Under the support of DOE-NE’s Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, efforts have been pursued to support HTGR technology development and its modeling and simulation needs. There is a particular need for advanced modeling and simulation tools to predict thermal-fluid behavior during safety-related transients. In our previous studies, the ring model was adopted in SAM and further developed to simulate the normal operating condition and a PCC event using the MHTGR-350 design of General Atomics as the reference design. This current work represents a continuation of these previous efforts, and the focus is to critically review and examine simplifications and assumptions made to develop the ring model, and to perform code validation using experimental data from an integral-effect test facility, the High Temperature Test Facility (HTTF) at the Oregon State University. In this study, the test PG-27 from the HTTF test suite was selected for code validation purpose. The test PG-27 is a transient test designed to simulate the PCC event of the MHTGR design. Very good agreements between SAM prediction and experimental measurements were found in both coolant and solid structure temperatures during the transient.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
NASA NTRS2015

Unsteady Three-Dimensional Simulation of a Shear Coaxial GO2/GH2 Rocket Injector with RANS and Hybrid-RAN-LES/DES Using Flamelet Models

Historically, the analysis and design of liquid rocket engines (LREs) has relied on full‐scale testing and one‐dimensional empirical tools. The testing is extremely expensive and the one‐dimensional tools are not designed to capture the highly complex, and multi‐dimensional features that are inherent to LREs. Recent advances in computational fluid dynamics (CFD) tools have made it possible to predict liquid rocket engine performance, stability, to assess the effect of complex flow features, and to evaluate injector‐driven thermal environments, to mitigate the cost of testing. Extensive efforts to verify and validate these CFD tools have been conducted, to provide confidence for using them during the design cycle. Previous validation efforts have documented comparisons of predicted heat flux thermal environments with test data for a single element gaseous oxygen (GO2) and gaseous hydrogen (GH2) injector. The most notable validation effort was a comprehensive validation effort conducted by Tucker et al. [1], in which a number of different groups modeled a GO2/GH2 single element configuration by Pal et al [2]. The tools used for this validation comparison employed a range of algorithms, from both steady and unsteady Reynolds Averaged Navier‐Stokes (U/RANS) calculations, large‐eddy simulations (LES), detached eddy simulations (DES), and various combinations. A more recent effort by Thakur et al. [3] focused on using a state‐of‐the‐art CFD simulation tool, Loci/STREAM, on a two‐dimensional grid. Loci/STREAM was chosen because it has a unique, very efficient flamelet parameterization of combustion reactions that are too computationally expensive to simulate with conventional finite‐rate chemistry calculations. The current effort focuses on further advancement of validation efforts, again using the Loci/STREAM tool with the flamelet parameterization, but this time with a three‐dimensional grid. Comparisons to the Pal et al. heat flux data will be made for both RANS and Hybrid RANSLES/ Detached Eddy simulations (DES). Computation costs will be reported, along with comparison of accuracy and cost to much less expensive two‐dimensional RANS simulations of the same geometry.

Westra, Doug G.↗