Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Thermal fluid physics”

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 145 records · Page 8

Fluid and gyrokinetic turbulence in open field-line, helical plasmas

Two-fluid Braginskii codes have simulated open-field line turbulence for over a decade, and only recently, it has become possible to study these systems with continuum gyrokinetic codes. This work presents a first-of-its-kind comparison between fluid and (long-wavelength) gyrokinetic models in open field-lines, using the GDB and Gkeyll codes to simulate interchange turbulence in the Helimak device at the University of Texas [T. N. Bernard et al., Phys. Plasmas 26, 042301 (2019)]. Partial agreement is attained in a number of diagnostic channels when the GDB sources and sheath boundary conditions (BCs) are selected carefully, especially the heat-flux BCs that can drastically alter the temperature. The radial profile of the fluctuation levels is qualitatively similar and quantitatively comparable on the low-field side, although statistics such as moments of the probability density function and the high-frequency spectrum show greater differences. This comparison indicates areas for future improvement in both simulations, such as sheath BCs, and improvements in GDB like particle conservation and spatially varying thermal conductivity, in order to achieve better fluid-gyrokinetic agreement and increase fidelity when simulating experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Analysis of convection heat transfer on multiscale rough superhydrophobic and liquid infused surfaces

Multiscale rough superhydrophobic or slippery liquid infused porous surfaces have gained much interest in recent years for their improved transport phenomena properties. While there have been several studies on drag reduction and condensation on non-wetting surfaces, convection heat transfer that is important in many thermal and thermochemical applications has not been addressed systematically. Here, this article utilizes a fractal description of rough surface topographies to develop analytical models for the Nusselt number and the thermal hydraulic factor for fluid flow and heat transfer inside a cylinder with non-wetting surfaces. For air-infused superhydrophobic surfaces, the model considers the dynamic stability of the air/fluid interface in the asperities. Using the analytical formulations and the stability criteria, systematic studies are presented on the effects of the fractal surface parameters, cylinder radius and Reynolds number on the convective heat transfer characteristics, from which surface texture design maps are developed for maximizing the convection heat transfer. It is shown that multiscale non-wetting surfaces are most effective in the range of lower Reynolds number and small cylinder radius for achieving the best convective heat transfer and thermal hydraulic performance. Applying the models to actual non-wetting surface topographies fabricated using electrodeposition and chemical etching, it is shown that contrary to prevailing notion, superhydrophobicity, characterized by the highest contact angles, does not always lead to the maximum convective heat transfer performance, and that under certain fluid flow conditions, hydrophobic surfaces may offer a greater thermal performance.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Heat transfer in directly-irradiated high-temperature solid–gas flows laden with polydisperse particles

Heat transfer in directly-irradiated high-temperature solid–gas flows laden with polydisperse particles is investigated using a novel transient three-dimensional computational fluid dynamics model. The model couples particle–gas hydrodynamics of solid–gas flows laden with polydisperse particles, radiative heat transfer in non-grey absorbing, emitting and anisotropically-scattering multi-component participating media, conduction heat transfer in the gas phase, and interfacial convection heat transfer. The multiphase particle-in-cell method is used to predict high-fidelity solid–gas flow characteristics, such as the local discrete particle size distribution, with increased computational efficiency by combining the advantages of both Eulerian and Lagrangian methods. The multi-component radiative transfer model is implemented using an advanced collision-based Monte Carlo ray-tracing method. The number of the prescribed discrete particle components is found to be the key parameter affecting the computational accuracy and efficiency, which primarily depends on the size distribution of the particles. For the model particle–gas flow featuring free-falling Gamma-distributed ceramic particles exposed to concentrated solar irradiation, the particle volume fraction, radiative, fluid flow and thermal characteristics appear to converge with the increasing number of the discrete particle components. Five particle components are sufficient to obtain physically meaningful results. A further increase in the number of the particle components only slightly increases the accuracy of the numerical predictions at the expense of a rapidly increasing computational time. For five particle components, the particle vertical velocity at the receiver exit for particles with the diameter of 43.4 μm is 57% of that for the particles with the diameter of 202.8 μm. The temperatures of these two particle components increase from the initial ambient values by factors of 2 and 1.2, respectively, during the simulation time. In conclusion, the model developed allows for increased fidelity of particle–gas flow simulations with significant radiative effects.

97 MATHEMATICS AND COMPUTING↗

Coupled radiation transport and plasma physics for ICF simulations [Slides]

Describing high-temperature plasmas requires inclusion of several nonlinear coupled physics, including ion kinetics - for each isotope species present in the plasma, electron kinetics, thermal radiative transfer (TRT), electromagnetic description, and nuclear reactions (fusion burn). Simulating plasmas in ICF hohlraums aims to solve kinetic ions (particles), fluid electrons, energetic electrons, electromagnetics, kinetic radiation (MG DP solver), and laser-plasma interactions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

SINDA/FLUINT and Thermal Desktop Multi-Node Settled and Unsettled Propellant Tank Modeling of Zero Boil Off Test

Cryogenic propellant storage tank self-pressurization involves complex physical phenomena which are usually analytically modelled via complex multidimensional CFD (Computational Fluid Dynamics) codes. Unfortunately these codes, even when modelling axisymmetric domains, may take weeks or longer to obtain transient pressure and temperature information for relatively short periods of time (several seconds to several hours). Propellant tank storage end-to-end mission simulations can last a duration of days to weeks to months. Multi-node modelling of propellant tanks is a viable alternative to traditional CFD modelling and presents the advantage of greatly reduced run times on the order of hours and days compared to the weeks or longer for CFD codes. A multi-node model represents the fluid within the storage tank, as well as the storage tank itself, as a fluid-thermal network. This type of setup is not necessarily geometrically based. This can be accomplished using a commercial generalized fluid-thermal network code, such as SINDA/FLUINT (SF). The advantage of using a fluid-thermal network code like SF lies in its extensive ability to model the external environment of the storage tank through the graphical user interface, Thermal Desktop (TD). The total heat load into the tank may be a function of heaters and a complex radiative environment as well. Thermal Desktop may be used to address the detailed radiative environment of the tank as well as building a geometrically accurate depiction of the storage tank itself.

Thermal Desktop↗

Peeling-ballooning modes in spherical tokamaks: Multi-branch instabilities and effects beyond ideal MHD

A number of important physics effects on the stability of relatively high-n (n is the toroidal mode number) peeling-ballooning modes (PBMs) are investigated based on an equilibrium reconstructed from a NSTX discharge, utilizing extended magnetohydrodynamic (MHD) eigenvalue solvers. For a given toroidal mode number n, multiple branches of instabilities are computed, with the total number of unstable branches roughly linearly scaling with n. Most of the unstable branches are located in the plasma core region, but edge-localized branches, i.e., PBMs, are also identified at higher n-numbers. For the single-fluid-wise most unstable PBM with n = 19⁠, stabilizing/destabilizing effects due to various physics beyond ideal MHD are systematically investigated. Plasma toroidal flow is found to be weakly stabilizing. Local flow shear is generally stabilizing as well, with the degree of stabilization depending on the initial growth rate (without flow shear) of the mode. The plasma resistivity can strongly destabilize the PBM within the single-fluid framework. Anisotropic thermal transport, strong parallel sound wave damping, as well as two-fluid effects are all stabilizing to the mode. In particular, diamagnetic stabilization (within the two-fluid model) is found to be very strong for this mode.

Linear stability analysis↗

Atmospheric cloud physics laboratory project study

Engineering studies were performed for the Zero-G Cloud Physics Experiment liquid cooling and air pressure control systems. A total of four concepts for the liquid cooling system was evaluated, two of which were found to closely approach the systems requirements. Thermal insulation requirements, system hardware, and control sensor locations were established. The reservoir sizes and initial temperatures were defined as well as system power requirements. In the study of the pressure control system, fluid analyses by the Atmospheric Cloud Physics Laboratory were performed to determine flow characteristics of various orifice sizes, vacuum pump adequacy, and control systems performance. System parameters predicted in these analyses as a function of time include the following for various orifice sizes: (1) chamber and vacuum pump mass flow rates, (2) the number of valve openings or closures, (3) the maximum cloud chamber pressure deviation from the allowable, and (4) cloud chamber and accumulator pressure.

Schultz, W. E.↗

Composition-transferable machine learning potential for LiCl-KCl molten salts validated by high-energy x-ray diffraction

Unraveling the liquid structure of multicomponent molten salts is challenging due to the difficulty in conducting and interpreting high-temperature diffraction experiments. Here, motivated by this challenge, we developed composition-transferable Gaussian approximation potential (GAP) for molten LiCl-KCl. A DFT-SCAN accurate GAP is active-learned from only ~1100 training configurations drawn from 10 unique mixture compositions enriched with metadynamics. The GAP-computed structures show strong agreement across high-energy x-ray diffraction experiments, including for a eutectic not explicitly included in model training, thereby opening the possibility of composition discovery.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Physical phenomena and the microgravity response

The living biological cell is not a sack of Newtonian fluid containing systems of chemical reactions at equilibrium. It is a kinetically driven system, not a thermodynamically driven system. While the cell as a whole might be considered isothermal, at the scale of individual macromolecular events there is heat generated, and presumably sharp thermal gradients exist at the submicron level. Basic physical phenomena to be considered when exploring the cell's response to inertial acceleration include particle sedimentation, solutal convection, motility electrokinetics, cytoskeletal work, and hydrostatic pressure. Protein crystal growth experiments, for example, illustrate the profound effects of convection currents on macromolecular assembly. Reaction kinetics in the cell vary all the way from diffusion-limited to life-time limited. Transport processes vary from free diffusion, to facilitated and active transmembrane transport, to contractile-protein-driven motility, to crystalline immobilization. At least four physical states of matter exist in the cell: aqueous, non-aqueous, immiscible-aqueous, and solid. Levels of order vary from crystalline to free solution. The relative volumes of these states profoundly influence the cell's response to inertial acceleration. Such subcellular phenomena as stretch-receptor activation, microtubule re-assembly, synaptic junction formation, chemotactic receptor activation, and statolith sedimentation were studied recently with respect to both their basic mechanisms and their responsiveness to inertial acceleration. From such studies a widespread role of cytoskeletal organization is becoming apparent.

Todd, Paul↗

Modeling Cyclic Phase Change and Energy Storage in Solar Heat Receivers

Numerical results pertaining to cyclic melting and freezing of an encapsulated phase change material (PCM), integrated into a solar heat receiver, have been reported. The cyclic nature of the present melt/freeze problem is relevant to latent heat thermal energy storage (LHTES) systems used to power solar Brayton engines in microgravity environments. Specifically, a physical and numerical model of the solar heat receiver component of NASA Lewis Research Center's Ground Test Demonstration (GTD) project was developed and results compared with available experimental data. Multi-conjugate effects such as the convective fluid flow of a low-Prandtl-number fluid, coupled with thermal conduction in the phase change material, containment tube and working fluid conduit were accounted for in the model. A single-band thermal radiation model was also included to quantify reradiative energy exchange inside the receiver and losses through the aperture. The eutectic LiF-CaF2 was used as the phase change material (PCM) and a mixture of He/Xe was used as the working fluid coolant. A modified version of the computer code HOTTube was used to generate results for comparisons with GTD data for both the subcooled and two-phase regimes. While qualitative trends were in close agreement for the balanced orbit modes, excellent quantitative agreement was observed for steady-state modes.

Hall, Carsie A., III↗

Space Nuclear Power Autonomous Control Algorithm and Control Element Test Bed

Nuclear thermal rockets are currently NASA’s preferred option for use in a manned mission to Mars in the 2040s. The communication delay between an Earth ground station and a spacecraft heading toward Mars can be up to 20 min. Therefore, controlling the nuclear rocket engine would require either a full-time reactor operator on the mission or an autonomous control system for the reactor. The latter idea of making space nuclear reactors fully autonomous has drawn more interest from stakeholders, but such an autonomous control system must be rigorously tested and validated before it is certified for human use. The cost of a full ground test for a space nuclear reactor is tremendous, so a nonnuclear mock reactor test bed was created to test and validate control elements and control algorithms for space nuclear reactors. The test bed consists of control element hardware that inputs physical measurement data into a reactor emulator to produce the reactor’s performance under steady-state, transient, and fault conditions. The control element hardware consists of six full-sized control drums equipped with servo drives and motors and is instrumented with optical encoders, resolvers, and torque sensors for drum movement characterization. In addition to the drums, a two-phase flow loop was designed and built to mimic the valves and turbomachinery associated with the propellant flow through a nuclear thermal rocket engine; components such as pressure sensors, flow meters, thermocouples, and tachometers are instrumented throughout the loop to characterize the fluid flow, valve, and turbomachinery behavior of the system. The data from the physical hardware (e.g., drum position, propellant flow rates) are input to a nuclear reactor simulator to determine the actual nuclear reactor parameters, and the data are sent back to a control algorithm to complete the control loop. The ability to conduct numerous tests of the control systems and autonomous algorithms can help validate the instrumentation and control aspects for a space nuclear reactor for every possible fault situation.

Wilson, Brandon↗

Manipulation of particles by weak forces

Quantitative relations between various force fields and their effects on the motion of particles of various sizes and physical characteristics were studied. The forces considered were those derived from light, heat, microwaves, electric interactions, magnetic interactions, particulate interactions, and sound. A physical understanding is given of the forces considered as well as formulae which express how the size of the force depends on the physical and electrical properties of the particle. The drift velocity in a viscous fluid is evaluated as a function of initial acceleration and the effects of thermal random motion are considered. A means of selectively sorting or moving particles by choosing a force system and/or environment such that the particle of interest reacts uniquely was developed. The forces considered and a demonstration of how the initial acceleration, drift velocity, and ultimate particle density distribution is affected by particle, input, and environmental parameters are tabulated.

Adler, M. S.↗

Optimization of contoured hypersonic scramjet inlets with a least-squares parabolized Navier-Stokes procedure

A new optimization procedure, in which a parabolized Navier-Stokes solver is coupled with a non-linear least-squares optimization algorithm, is applied to the design of a Mach 14, laminar two-dimensional hypersonic subscale flight inlet with an internal contraction ratio of 15:1 and a length-to-throat half-height ratio of 150:1. An automated numerical search of multiple geometric wall contours, which are defined by polynomical splines, results in an optimal geometry that yields the maximum total-pressure recovery for the compression process. Optimal inlet geometry is obtained for both inviscid and viscous flows, with the assumption that the gas is either calorically or thermally perfect. The analysis with a calorically perfect gas results in an optimized inviscid inlet design that is defined by two cubic splines and yields a mass-weighted total-pressure recovery of 0.787, which is a 23% improvement compared with the optimized shock-canceled two-ramp inlet design. Similarly, the design procedure obtains the optimized contour for a viscous calorically perfect gas to yield a mass-weighted total-pressure recovery value of 0.749. Additionally, an optimized contour for a viscous thermally perfect gas is obtained to yield a mass-weighted total-pressure recovery value of 0.768. The design methodology incorporates both complex fluid dynamic physics and optimal search techniques without an excessive compromise of computational speed; hence, this methodology is a practical technique that is applicable to optimal inlet design procedures.

Korte, J. J.↗

Propagation of Input Uncertainties in Numerical Simulations of Laser Powder Bed Fusion

Laser powder bed fusion has the potential of redefining state-of-the-art processing and production methods, but defect formation and inconsistent build quality have limited the implementation of this process. Numerical models are widely used to study this process and predict the formation of these defects. Presently, the uncertainties of model input parameters and thermophysical properties used by these numerical simulations have not been investigated. In the present study, the uncertainty in these input parameters and material properties are quantified for laser powder bed fusion, with and without a simulated powder bed, to determine their influence on the predictive accuracy of an experimentally validated numerical model. Accounting for all possible sources of uncertainty quickly becomes computationally expensive on account of the curse of dimensionality. Uncertainty in laser absorption, solid, and liquid specific heat of the metal were found to have the largest effect on model prediction reliability with or without the use of a powder bed. Results also illustrate that accounting for these three uncertain parameters still captures the majority of model prediction uncertainty. Furthermore, the methodology of this study may be used to understand the uncertainty in as-built microstructure through propagation to microstructure prediction models, or applied under processing conditions where high Péclet numbers are observed and the thermal convection and fluid flow within the molten pool are substantial.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Comprehensive Review of DPF-based Flash Neutron Radiography Viability

MJOLNIR (MegaJOule Neutron Imaging Radiography) is a Dense Plasma Focus (DPF) being developed by LLNL as a prototype to assess the viability of a DPF as a flash neutron radiography source. To date, MJOLNIR has discharged up to 1.3 MJ of stored energy into a deuterium plasma load and achieved yields up to 1.2 × 10 12 neutrons per discharge. The MJOLNIR pulse length already meets the preliminary requirements for flash neutron radiography and we demonstrate in this report plausible paths forward to meet the yield requirement in a deuterium plasma. MJOLNIR’s neutron spot size has been characterized to be a factor of two greater than radiography requirements and presents the greatest challenge. A path forward to sufficiently shrinking spot size is outlined in this report, using a combination of gas dopants and a smaller anode implosion radius. A DPF is a plasma device with coaxial electrodes whose discharge ends with a stagnated hot and dense plasma column on-axis, at the tip of the central anode. Inside the MJOLNIR DPF plasma column, both thermal and beam target processes generate neutrons. To understand the underlying physics of neutron generation inside the DPF, we model the experiment using a combination of kinetic, fluid, and reduced-order models we have been developing for over ten years at LLNL. Our numerical tools are also pivotal to charting the path forward because they allow us to numerically test out which modifications would improve the DPF’s characteristic output.

42 ENGINEERING↗

Physics Based Model for Cryogenic Chilldown and Loading. Part I: Algorithm

We report the progress in the development of the physics based model for cryogenic chilldown and loading. The chilldown and loading is model as fully separated non-equilibrium two-phase flow of cryogenic fluid thermally coupled to the pipe walls. The solution follow closely nearly-implicit and semi-implicit algorithms developed for autonomous control of thermal-hydraulic systems developed by Idaho National Laboratory. A special attention is paid to the treatment of instabilities. The model is applied to the analysis of chilldown in rapid loading system developed at NASA-Kennedy Space Center. The nontrivial characteristic feature of the analyzed chilldown regime is its active control by dump valves. The numerical predictions are in reasonable agreement with the experimental time traces. The obtained results pave the way to the development of autonomous loading operation on the ground and space.

heat transfer↗

Modeling Validation and Its Design Relevance for NTP Fuel Elements

The development of Nuclear Thermal Propulsion engines depends heavily on modeling and simulation of a number of disciplines—neutronic / fluid / thermal / structural are only a start. Further, rockets involve extreme conditions, particularly high material temperatures in the nuclear fuel. This paper considers validation of the underling physical models and codes and how their accuracy influences several design challenges. Several issues are identified, where models, material data, could enable design solutions. In particular, thermal and CTE mismatch stresses are a concern, and stress predictions would benefit from improved mechanical property data at high temperatures. Nuclear inter-element heating causes propellant flow maldistribution and performance risks, but a simple model provides insights. High temperature chemical diffusion should be modeled to understand fuel element mass loss.

Nuclear Propulsion↗

Modeling Validation and Its Design Relevance for NTP Fuel Elements

The development of Nuclear Thermal Propulsion engines depends heavily on modeling and simulation of a number of disciplines—neutronic / fluid / thermal / structural are only a start. Further, rockets involve extreme conditions, particularly high material temperatures in the nuclear fuel. This paper considers validation of the underling physical models and codes and how their accuracy influences several design challenges. Several issues are identified, where models, material data, could enable design solutions. In particular, thermal and CTE mismatch stresses are a concern, and stress predictions would benefit from improved mechanical property data at high temperatures. Nuclear inter-element heating causes propellant flow maldistribution and performance risks, but a simple model provides insights. High temperature chemical diffusion should be modeled to understand fuel element mass loss.

Nuclear Propulsion↗