Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Mass diffusivity”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
NASA NTRS1999

Control of Meridional Flow by a Non-Uniform Rotational Magnetic Field

The diffusive mass transfer of species during crystal growth in vertical ampoules is significantly affected by fluid flow in the liquid mother phase (melt). For electrically conductive melts, an elegant way of remotely inducing and controlling this flow is by utilizing a uniform rotational magnetic field (RMF) in the transverse direction. It induces an azimuthal flow which tends to homogenize the thermal and solutal fields. The rotating field also reduces the diffusion boundary layer, stabilizes temperature fluctuations, and promotes better overall crystal growth. For moderate strengths of the applied magnetic field (2-20 m Tesla) with frequencies of up to 400 Hz, the induced secondary meridional flow becomes significant. It typically consists of one roll at the bottom of the liquid column and a second roll (vortex) at the top. The flow along the centerline (ampoule axis) is directed from the growing solid (interface) towards the liquid (melt). In case of convex interfaces (e.g. in floating zone crystal growth) such flow behavior is beneficial since it suppresses diffusion at the center. However, for concave interfaces (e.g. vertical Bridgman crystal growth) such a flow tends to exacerbate the situation in making the interface shape more concave. It would be beneficial to have some control of this meridional flow- for example, a single recirculating cell with controllable direction and flow magnitude will make this technique even more attractive for crystal growth. Such flow control is a possibility if a non-uniform PNE field is utilized for this purpose. Although this idea has been proposed earlier, it has not been conclusively demonstrated so far. In this work, we derive the governing equations for the fluid dynamics for such a system and obtain solutions for a few important cases. Results from parallel experimental measurements of fluid flow in a mercury column subjected to non-uniform RMF will also be presented.

Mazuruk, Konstantin↗
DOE OSTI2023

Direct numerical simulation of a turbulent methane/air flame impinged by a sub-breakdown electric field

Here, the effects of an external electric field on a turbulent methane/air diffusion flame are analyzed in this work using direct numerical simulations. The analyzed configuration consists of a temporally evolving mixing layer of air and a mixture of methane and nitrogen at 1 atm that is impinged by an electric field in the direction normal to the mean mixing plane. The combustion and chemi-ionization reactions involved in the flow are modeled using finite-rate chemistry and a reduced reaction scheme consisting of 26 species and 134 reactions. The mass diffusion and ion-wind effects are modeled using a detailed description of the diffusion coefficients and electric mobilities based on kinetic theory. The presented calculations show that the turbulence generated within the mixing layer is mostly unaffected by the applied electric field for the configuration under exam. In fact, the electric body force is developed away from the mixing region where the flow is uniform. Conversely, the turbulence is able to introduce very high intermittency in the electrically charged species concentration and, consequently, in the electric body force. Such intermittency will constitute a challenge in the future formulation of combustion models that take into account ion-wind effects.

09 BIOMASS FUELS↗
NASA NTRS1991

A model of the evaporation of binary-fuel clusters of drops

A formulation has been developed to describe the evaporation of dense or dilute clusters of binary-fuel drops. The binary fuel is assumed to be made of a solute and a solvent whose volatility is much lower than that of the solute. Convective flow effects, inducing a circulatory motion inside the drops, are taken into account, as well as turbulence external to the cluster volume. Results obtained with this model show that, similar to the conclusions for single isolated drops, the evaporation of the volatile is controlled by liquid mass diffusion when the cluster is dilute. In contrast, when the cluster is dense, the evaporation of the volatile is controlled by surface layer stripping, that is, by the regression rate of the drop, which is in fact controlled by the evaporation rate of the solvent. These conclusions are in agreement with existing experimental observations. Parametric studies show that these conclusions remain valid with changes in ambient temperature, initial slip velocity between drops and gas, initial drop size, initial cluster size, initial liquid mass fraction of the solute, and various combinations of solvent and solute. The implications of these results for computationally intensive combustor calculations are discussed.

Harstad, K.↗
NASA NTRS2000

Effects of G-Jitter on Interfacial Dynamics of Two Miscible Liquids: Application of MIM

We designed an experiment to examine the effects of g-jitter on mixing of two miscible liquids using the Microgravity Vibration Isolation Mount (MIM). The global bifurcation of the interface was observed with the MIM operating alternatively to either transmit the g-jitter, isolate from the g-jitter or to provide controlled vibration levels with well defined amplitude and frequency content. With the MIM in isolation mode, the interface remains stationary indicating buoyancy induced convection is negligibly small such that mixing occurs via intrinsic mass diffusion without the masking effect of vibration driven convection. Analytical and computational results are in agreement with the experimental findings. Operation of the MIM in forced mode with conditions typical of g-jitter shows that vibration induced convective flows can excite instability mechanisms such as Kelvin-Helmholtz to generate large amplitude quasi-stationary waves oriented vertically for various cases with Stokes-Reynolds number in the range of 0.003 to 0.5. The two and four mode quasi-stationary waves are also predicted with a mathematical model. Though unplanned, the effect of a primary thruster filing was captured and shown to cause a catastrophic bifurcation, enhancing local mass transport. In light of the findings, experiments planned for the International Space Station should consider the potential effects of g-jitter.

Duval, Walter M. B.↗
NASA NTRS1999

Dynamics of Single Chains of Suspended Ferrofluid Particles

We present an experimental study of the dynamics of isolated chains made of super-paramagnetic particles under the influence of a magnetic field. The motivation of this work is to understand if the chain fluctuations exist and, if it does, how does the fluctuation affect chain aggregation. We find that single chains strongly fluctuate and that the characteristic frequency of their fluctuations is inversely proportional to the magnetic field strength. The higher the field the lower the characteristic frequency of the chain fluctuations. In the high magnetic field limit, chains behave like rigid rods without any internal motions. In this work, we used ferrofluid particles suspended in water. These particles do not have any intrinsic magnetization. Once a magnetic field is applied, a dipole moment is induced in each particle, proportional to the magnetic field. A dipolar magnetic interaction then occurs between particles. If dipole-dipole magnetic energy is higher than the thermal energy, the result is a structure change inside the dipolar fluid. The ratio of these two energies is expressed by a coupling constant lambda as: lambda = (pi(a(exp 3))(chi(exp 2))(mu(sub 0))(H(sub 0))(exp 2))/18kT Where a is the particle radius, mu(sub 0) is the vacuum magnetic permeability, H(sub 0) the applied magnetic field, k the Boltzmann constant and T the absolute temperature. If lambda > 1, magnetic particles form chains along the field direction. The lateral coalescence of several chains may form bigger aggregates especially if the particle volume fraction is high. While many studies and applications deal with the rheological properties and the structural changes of these dipolar fluids, this work focuses on the understanding of the chain dynamics. In order to probe the chain dynamics, we used dynamic light scattering (DLS) in self-beating mode as our experimental technique. The experimental geometry is such that the scattering plane is perpendicular to the magnetic field. Therefore, only motions in this plane are probed. A very dilute sample of a ferrofluid emulsion with a particle volume fraction of 10(exp -5) is used in this experiment. We chose such a low volume fraction to avoid multiple light scattering as well as lateral chain-chain aggregation. DLS measures the dynamic structure factor S(q,t) of the sample (q is the scattering wave vector, t is the time). In the absence of the magnetic field, identical particles of ferrofluid droplets are randomly distributed and S(q,t) reduces to exp(-q(exp 2)2D(sub 0)t). D(sub 0)=(kT/(6(pi)(eta)(a)) is the diffusion coefficient of Brownian particles (where Xi = (6(pi)(eta)(a)) is the Stokes frictional coefficient of a spherical particle in a fluid of viscosity eta). If interactions or polydispersity can not be ignored, an effective diffusion coefficient is introduced. Formally, D(sub eff) is defined as: D(sub eff) = - q(exp -2) partial derivative of (ln(S(q,t)) with respect to time, as t goes to 0. D(sub eff) reduces to D(sub 0) if no interactions and only a few particles size are present. Therefore, we can use DLS to measure particle size. The particle radius was found to be a=0.23 mu m with 7% of polydispersity. In this case, if we vary the scattering angle theta (and so q) we do not have any change in the measured diffusion coefficient: it is q-independent. When a magnetic field is applied, particles aggregate into chains if lambda > 1. We first studied the kinetics of the chain formation when lambda = 406. At a fixed scattering angle, we measured diffusion coefficient D(sub eff) as a function of time. Experimentally, we find that D(sub eff) decreases monotonously with time. Physically, this means that chains are becoming longer and longer. Since we are only sensitive to motions in the scattering plane and since chains have their main axis perpendicular to this plane, the measured diffusion coefficient is the trans-verse diffusion coefficient. We can relate D(sub eff) to the mean number of particles per chain N(t) at a given time and to the diffusion coefficient of an isolated particle D(sub 0) as D(sub eff)=f(N(t))D(sub 0). Since f(N) is known from other recent work, N can be expressed as a function of the time. We found a square root dependency: N(t) proportional to the square root of t. As expected for very low volume fraction, this behavior is characteristic of a diffusion-limited aggregation as suggested by several authors and by our previous work. In this study, we focus on the dependence of the effective diffusion coefficient on the scattering angle and the magnetic field strength. After the magnetic field is applied (lambda = 406) for a long time, typically 6 hours, kinetics of chain formation becomes very slow. Chain size does not vary much over the next hour period. Thus, we can perform different interesting experiments. First, at a fixed magnetic field, we measure the effective diffusion coefficient as a function of the scattering angle (from 5 to 130 deg). Our results show that the measured diffusion coefficient increases linearly with the scattering angle: D(sub eff) proportional to q. If we do the same experiment for different lambda values, D(sub eff) depends on lambda as D(sub eff) proportional to lambda(exp -1/2). We also find for different lambda values that the same asymptotic D(sub eff) value is obtained when q approaches zero. The angle dependency of D(sub eff) suggests that an additional motion exists besides chain drifting. Chain size is constant during experiment, which was verified by measuring the same diffusion coefficient at the beginning and at the end of the angle switching. If chains are rigid, D(sub eff) is independent of q. Therefore, we found that D(sub eff) not only measures the motion of the entire chain but also its internal fluctuations. These internal motions are the fluctuations of the particles in the chain. To understand the q dependency of D(sub eff), let us look at the probing length used. In our study, the characteristic length scale probed is l=2pi/q which is in the range of 0.9<l/a<20. When l is much larger than the particle radius (l/a going to infinity) we are mainly sensitive to the center of mass diffusion of the chain. D(sub eff) is then the diffusion of the entire chain and depends only on N but not on q and lambda. The value D(sub eff) allows us to obtain the number of particles per chain N. In the opposite limit (i.e. l/a < 1), we are sensitive to motions on the size of individual particles. The main contributions to the measured diffusion coefficient come from internal motions of the chains (i.e. particles' fluctuations inside the chain). We investigated also the effect of the chain size on D(sub eff). For the same value of the magnetic interaction lambda and different chain sizes, we found that the slope of D(sub eff) versus q decreases when N increases as a power law D(sub eff) proportional to N(exp -0.7). This is an important behavior because this means that the longer the chain the slower its internal fluctuations besides the whole chain motion.

Cutillas, S.↗
DOE OSTI2022

Matrimid substrates with bicontinuous surface and macrovoids in the bulk: A nearly ideal substrate for composite membranes in CO 2 capture

Scalable fabrication of thin-film composite (TFC) membrane for post-combustion carbon capture is often limited by the availability of high-performance nanoporous substrate. Ideally, the substrate should allow for fast gas transport at the selective layer/substrate interface as well as in the bulk of the substrate. In this study, highly permeable Matrimid substrates were prepared via vapor- and nonsolvent-induced phase separations. The phase separation mechanism of the Matrimid®/N-methyl-2-pyrrolidone casting solution was modulated by the addition of LiCl to control the solution thermodynamic stability, and nanoporous Matrimid substrates were formed with a bicontinuous surface and macrovoids in the bulk. This hierarchically optimized structure resulted in a CO 2 permeance of 2.60 × 10 5 GPU (1 GPU = 1 × 10 –6 cm 3 (STP) cm –2 s –1 cmHg –1 = 3.349 × 10 –10 mol m –2 s –1 Pa –1 ), which was ca. 11 times more permeable than a benchmark polyethersulfone substrate with cellular pores. The benefit of using this new substrate was demonstrated by coating a 170-nm amine-containing polymer to form a TFC facilitated transport membrane. The membrane exhibited a CO 2 permeance of 932 GPU at 57°C, which was 72 GPU higher than the counterpart coated on the benchmark substrate. Meanwhile, the CO 2 /N 2 selectivity was remained at 158. This permeance improvement could be well explained by the resistance-in-series model, where the improved permeance was attributed to the reductions in substrate and lateral diffusion mass transfer resistances. The upper bound (UB) analysis indicates that the Matrimid substrate is nearly ideal for the state-of-the-art polymers for CO 2 /N 2 separation. Here, the substrate improvement could also reduce the parasitic energy associated with the membrane process due to the reduced membrane cost and flue gas compression requirement.

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

Neutron capture of UO 2 targets prepared by spin-coating assisted combustion synthesis

Two uranium dioxide (UO 2 ) targets of (414 ± 23) nm and (1092 ± 93) nm thicknesses were prepared on 6061 aluminum alloy and puratronic grade aluminum backing materials. The targets were deposited with a novel method combining spin coating and solution combustion synthesis (SCS). The target layers consisted of small (3–7 nm) UO 2 grains and uniformly distributed ultra-small (1–3 nm) pores. The prepared targets were tested at the Los Alamos National Laboratory’s LANSCE facility for neutron irradiation damage and suitability for neutron capture experiments. The samples showed no signs of target material loss after the irradiation. However, irradiation caused a significant increase in the grain size (4–10 nm), as well as upward mass diffusion and coalescence of the pores due to the thermal spikes. The magnesium in the aluminum 6061 alloy backing also diffused into the UO 2 layer during neutron irradiation. The structural changes in the target after the irradiation do not affect the data from neutron capture. As a result, the new method can be used more broadly to prepare other actinide targets for nuclear physics experiments.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Versatile TRISO fuel particle modeling in Bison

Tri-structural isotropic (TRISO) fuel particles are a key component in several previous and current reactors as well as in a variety of novel nuclear reactor designs. Interest in TRISO fuel is on the rise, necessitating considerable computer modeling of TRISO fuel behavior in order to support related design and licensing activities. The Bison nuclear fuel performance code, which offers a full set of capabilities for modeling TRISO fuels, makes it easier to explore the various important aspects of TRISO fuel behavior. One key advantage of Bison is its ability to create meshes in 1D, 2D, and 3D. Users can customize these meshes for specific geometries, mesh densities, and use cases. This enables a wide variety of analyses, including thermal, structural, mass diffusion, homogenization, and statistical failure analyses. Furthermore, the meshing capability simplifies analysts’ workflows. The inherent mesh generation capability eliminates the need for separate mesh-generating software and mesh file management. Also, the fact that the meshes are customizable makes it straightforward to automate an investigation over a range of geometric parameters or mesh densities. Here, the present paper highlights the ease with which Bison may be used to create meshes for both simple and relatively complex TRISO fuel particles, and it explores the types of analyses enabled by these meshes.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Simulation and flow physics of a shocked and reshocked high-energy-density mixing layer

This work describes a computational investigation of multimode instability growth and multimaterial mixing induced by multiple shock waves in a high-energy-density (HED) environment, where pressures exceed 1 Mbar. The simulations are based on a series of experiments performed at the National Ignition Facility (NIF) and designed as an HED analogue of non-HED shock-tube studies of the Richtmyer–Meshkov instability and turbulent mixing. A three-dimensional computational modelling framework is presented. It treats many complications absent from canonical non-HED shock-tube flows, including distinct ion and free-electron internal energies, non-ideal equations of state, radiation transport and plasma-state mass diffusivities, viscosities and thermal conductivities. The simulations are tuned to the available NIF data, and traditional statistical quantities of turbulence are analysed. Integrated measures of turbulent kinetic energy and enstrophy both increase by over an order of magnitude due to reshock. Large contributions to enstrophy production during reshock are seen from both the baroclinic source and enstrophy–dilatation terms, highlighting the significance of fluid compressibility in the HED regime. Dimensional analysis reveals that Reynolds numbers and diffusive Péclet numbers in the HED flow are similar to those in a canonical non-HED analogue, but conductive Péclet numbers are much smaller in the HED flow due to efficient thermal conduction by free electrons. It is shown that the mechanism of electron thermal conduction significantly softens local spanwise gradients of both temperature and density, which causes a minor but non-negligible decrease in enstrophy production and small-scale mixing relative to a flow without this mechanism.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Strain-Based Assessment of Shale Caprock during Cyclic Underground Hydrogen Storage

Successful large-scale underground hydrogen storage (UHS) in depleted gas reservoirs depends on the integrity of the overlying caprock to prevent hydrogen loss during cyclic injection and depletion. Prior studies on crushed Marcellus shale, a potential caprock, indicate that cyclic hydrogen injection and depletion induces microstructural changes, increasing porosity and permeability. However, the extent of these changes in intact shale remains unclear. This study presents a strain-based experimental approach to quantify volumetric strain evolution in intact Marcellus shale matrix under unconstrained stress conditions. A quadrant-shaped shale sample without visible fractures underwent eight hydrostatic pore-pressure cycles (injection to 1500 psi and depletion to 500 psi in 250 psi steps). Linear strain gauges measured strain in three orthogonal directions. Results indicate progressive plastic strain accumulation, leading to an ∼12% increase in matrix porosity after eight cycles, with an estimated 19% increase after 30 cycles. This porosity increase follows a logarithmic trend, suggesting a diminishing effect in later cycles. Additionally, permeability and diffusive mass flux are projected to rise by ∼70% over 30 cycles, enhancing hydrogen migration risk. The shale matrix also exhibited mechanical stiffening over successive cycles, limiting large-scale deformation but not preventing porosity enhancement. A new parameter, α, was introduced to characterize shale sensitivity to cyclic loading, aiding UHS caprock assessments. These findings underscore the necessity of incorporating cyclic loading effects in UHS site selection and operational strategies to ensure long-term storage integrity.

08 HYDROGEN↗
DOE OSTI2023

On-Demand Aligned DNA Hydrogel Via Light Scanning

DNA is an anisotropic, water-attracting, and biocompatible material, an ideal building block for hydrogel. The alignment of the anisotropic DNA chains is essential to maximize hydrogel properties, which has been little explored. Here, we present a method to fabricate the anisotropic DNA hydrogel that allows precise control for the polymerization process of photoreactive cationic monomers. Scanning ultraviolet light enables the uniaxial alignment of DNA chains through the polymerization-induced diffusive mass flow using a concentration gradient. In conclusion, while studying anisotropic mechanical properties and orientation recovery according to the DNA chain alignment direction, we demonstrate the potential of directionally controlled DNA hydrogels as smart materials.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A Decentralized Approach for Modeling Organized Convection Based on Thermal Populations on Microgrids

Abstract In this study, a spectral model for convective transport is coupled to a thermal population model on a two‐dimensional horizontal “microgrid,” covering the typical gridbox size of general circulation models. The goal is to explore new ways of representing impacts of spatial organization in cumulus cloud fields. The thermals are considered the smallest building block of convection, with thermal life cycle and movement represented through binomial functions. Thermals interact through two simple rules, reflecting pulsating growth and environmental deformation. Long‐lived thermal clusters thus form on the microgrid, exhibiting scale growth and spacing that represent simple forms of spatial organization and memory. Size distributions of cluster number are diagnosed from the microgrid through an online clustering algorithm, and provided as input to a spectral multiplume eddy‐diffusivity mass flux scheme. This yields a decentralized transport system, in that the thermal clusters acting as independent but interacting nodes that carry information about spatial structure. The main objectives of this study are (a) to seek proof of concept of this approach, and (b) to gain insight into impacts of spatial organization on convective transport. Single‐column model experiments demonstrate satisfactory skill in reproducing two observed cases of continental shallow convection. Metrics expressing self‐organization and spatial organization match well with large‐eddy simulation results. We find that in this coupled system, spatial organization impacts convective transport primarily through the scale break in the size distribution of cluster number. The rooting of saturated plumes in the subcloud mixed layer plays a key role in this process.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Comparative Evaluation of the Ability of the MYNN‐EDMF PBL Scheme in WRF Model to Reproduce Near Surface Wind Speed Over Different Topographical Types

Abstract This study systematically evaluates the performance of the Mellor‐Yamada‐Nakanishi‐Niino‐Eddy‐Diffusion‐Mass‐Flux planetary boundary layer (PBL) scheme within the Weather Research and Forecasting (WRF) model in simulating near‐surface wind speeds across various topographies in New York State (NYS). Simulated wind speeds are compared with in‐situ measurements from 22 surface sites, grouped into six topographic categories: continental plain (CT), lakeside (LS), river valley (RV), Long Island (LI), Block Island (BI), and offshore ocean (OO). A quantitative evaluation based on Relative Euclidean Distance shows that wind speeds at the OO site are the most accurately reproduced, followed by those at LI sites, while the model performs less accurately for the remaining topographic groups. Wind speeds over CT sites tend to be overestimated by approximately 1 m/s, although their diurnal variability (DV) is well captured. In contrast, the model underestimates wind DV at LS, RV, LI, and BI sites, with the largest biases occurring at LI and BI, resulting in underestimated daytime wind speed and/or overestimated nighttime wind speed. The OO winds exhibit minimal diurnal variation, accurately captured by our WRF model. The surface wind diurnal variation is closely linked to PBL development. Among the indicators of PBL development, surface potential temperature biases most strongly correlate with wind speed biases. Our WRF model faces challenges in capturing the distinctions between winds influenced by local circulations and those over continental plains, and the significantly stronger winds at OO compared to BI. Potential causes for these biases are discussed, offering pathways for improving surface wind simulations in future.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Ionic Peltier effect in Li-ion electrolytes

The coupled transport of charge and heat provide fundamental insights into the microscopic thermodynamics and kinetics of materials. We describe a sensitive ac differential resistance bridge that enables measurements of the temperature difference on two sides of a coin cell with a resolution of better than 10 μK. We use this temperature difference metrology to determine the ionic Peltier coefficients of symmetric Li-ion electrochemical cells as a function of Li salt concentration, solvent composition, electrode material, and temperature. The Peltier coefficients Π are negative, i.e., heat flows in the direction opposite to the drift of Li ions in the applied electric field, large, –Π > 30 kJ mol –1 , and increase with increasing temperature at T > 300 K. The Peltier coefficient is approximately constant on time scales that span the characteristic time for mass diffusion across the thickness of the electrolyte, suggesting that heat of transport plays a minor role in comparison to the changes in partial molar entropy of Li at the interface between the electrode and electrolyte. Here, our work demonstrates a new platform for studying the non-equilibrium thermodynamics of electrochemical cells and provides a window into the transport properties of electrochemical materials through measurements of temperature differences and heat currents that complement traditional measurements of voltages and charge currents.

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

Micropatterning of alumina tubular membranes via laser carving for enhanced direct air capture

An approach for preparing and applying micropatterned alumina tubular membranes was developed for improved gas-liquid contact at low air pressure drop in direct air capture application. The paper demonstrated that the laser carving of micropatterns on the outer surface of the alumina tubular membranes could greatly increase their outer surface area and enhance gas-liquid turbulence for a reduced diffusion mass transfer resistance. Six kinds of micropattern configurations were fabricated and studied, including random and regular micropatterns. By introducing a 500-μm solid grid micropattern, the outer surface area doubles and enhances carbon dioxide capture efficiency from 61% to 97% after the membrane was hydrophobically modified with fluoroalkylsilane. Here, the air pressure drop through the tube lumen remained low even when the packing density increased from 382 to 906 m 2 /m 3 . The liquid entry pressure of this micropatterned membrane was the same after testing for 220 h running with the help of periodic drying.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2021

Application of Scalar Filtered Density Function to Turbulent Flows Under Supercritical Condition

Abstract The scalar filtered density function (FDF) methodology is extended and employed for large eddy simulation (LES) of turbulent flows under supercritical condition. To describe real fluid behavior, the extended methodology incorporates the generalized heat and mass diffusion models along with real fluid thermodynamic relations which are derived using the cubic Peng–Robinson equation of state. These models are implemented within the stochastic differential equations comprising the scalar FDF transport. Simulations are conducted of a temporally developing mixing layer under supercritical condition and the results are assessed by comparing with data generated by direct numerical simulation (DNS) of the same layer. The consistency of the proposed FDF methodology is assessed. The LES-FDF predictions are shown to agree favorably with the DNS data and exhibit several key features pertaining to supercritical turbulent flows.

Energy & Fuels↗
DOE OSTI2021

TMAP7

TMAP7 was developed by the Fusion Safety Program at INL primarily as a safety analysis code, to analyze tritium retention and loss in fusion reactor structures and systems during normal operation and accident conditions. TMAP7 solves one-dimensional thermal- and mass-diffusive transport and trapping equations for transport through structures, across the interface between enclosures and structures, and zero dimensional fluid transport between enclosures.

Humrickhouse, PaulW.↗