SEARCH · Engineering Papers

Results for “diffusion length”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,135 records · Page 63

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Validation of the HyRAM+ physics models for use with propane

The Hydrogen Risk Assessment Model Plus (HyRAM+) toolkit combines quantitative risk assessment with simulations of unignited dispersion, ignited turbulent diffusion flames, and indoor accumulation with delayed ignition of fuels. HyRAM+ is differentiated from HyRAM in that it includes models and leak data for other alternate fuels. The models of the physical phenomena need to be validated for each of the fuels in the toolkit. This report shows the validation for propane which is being used as a surrogate for autogas, which is a mixture of propane and butane and used in internal combustion engines in vehicles. For flame length comparisons, five previously published experiments from peer reviewed journals were used to validate our models. The validation looked at flame lengths and flame widths with respect to different leak diameters, mass flow rates, and source pressures. Most of the sources included more than one set of experimental data, which were collected using different methods (CCD cameras, IR visualization etc.). In general, HyRAM+ overpredicts the flame lengths by around 65%. For heat and radiation models, we compared the heat flux and radiation data reported from two different sources to the values calculated by HyRAM+. For higher mass flow rates, the HyRAM+ calculated flame length results gave a better estimate of what is found in the experiments (65% error), but a higher error (85%) is observed between the HyRAM+ calculated lengths and the experimental flame lengthsfor lower mass flows. Some differences can be attributed to outdoor environmental effects (i.e. wind speed) and uncertainties in jet flame shapes. The propane flame trajectory is predicted for a high Reynolds number case with Re = 12,500 and a low Reynolds number case where Re = 2,000. The Re=12,500 case which is momentum dominated matches well with the experimental flame trajectory, but the agreement for the bouancy driven low Reynolds number case is not as good. Dispersion modeling for unignited propane was also analyzed. We compared the mole fraction, mixture fraction, mean velocity, concentration half width, and inverse mass concentration over an axial distance from different credible journals to the values calculated by HyRAM+. The results display good agreement but generally, HyRAM+ predicts a wider profile for mole fraction and mixture fraction experiments. Overall, HyRAM+’s results are reasonable for predicting the flame length, heat flux, flame trajectory, and dispersion for propane and can be used in risk analyses

08 HYDROGEN↗

NASA NTRS2000

Smoke-Point Properties of Nonbuoyant Round Laminar Jet Diffusion Flames

The laminar smoke-point properties of nonbuoyant round laminar jet diffusion flames were studied emphasizing results from long duration (100-230 s) experiments at microgravity carried -out on- orbit in the Space Shuttle Columbia. Experimental conditions included ethylene-and propane-fueled flames burning in still air at an ambient temperature of 300 K, initial jet exit diameters of 1.6 and 2.7 mm, jet exit velocities of 170-1630 mm/s, jet exit Reynolds numbers of 46-172, characteristic flame residence times of 40-302 ms, and luminous flame lengths of 15-63 mm. The onset of laminar smoke-point conditions involved two flame configurations: closed-tip flames with first soot emissions along the flame axis and open-tip flames with first soot emissions from an annular ring about the flame axis. Open-tip flames were observed at large characteristic flame residence times with the onset of soot emissions associated with radiative quenching near the flame tip; nevertheless, unified correlations of laminar smoke-point properties were obtained that included both flame configurations. Flame lengths at laminar smoke-point conditions were well-correlated in terms of a corrected fuel flow rate suggested by a simplified analysis of flame shape. The present steady and nonbuoyant flames emitted soot more readily than earlier tests of nonbuoyant flames at microgravity using ground-based facilities and of buoyant flames at normal gravity due to reduced effects of unsteadiness, flame disturbances and buoyant motion. For example, laminar smoke-point flame lengths from ground-based microgravity measurements were up to 2.3 times longer and from buoyant flame measurements were up to 6.4 times longer than the present measurements at comparable conditions. Finally, present laminar smoke-point flame lengths were roughly inversely proportional to pressure, which is a somewhat slower variation than observed during earlier tests both at microgravity using ground-based facilities and at normal gravity.

Urban, D. L.↗

DOE OSTI2023

Elastic Flow Modeling for Hydropower Digital Twins

This report details the elastic, unsteady, one-dimensional flow equations that are used to model flow through a penstock in a hydropower facility. The elastic flow model is accurate even in cases of fast transients, long penstock length, and high gravitation head. The compressibility of the water and elasticity of the pipe walls are explicitly accounted for in the developed models so that the water hammer phenomena can be accurately captured. The elastic flow model is coupled to a mechanistic turbine model to resolve the dynamic feedbacks between the elastic water column and the turbine rotation rate. A finite volume method is employed to solve the governing equations, and the method is shown to have low numerical diffusion in sharp gradient phenomena, such as those encountered when simulating water hammer. The method is applied to single dimensional linear advection benchmark problem and numerical results are compared with analytic results. This is followed on by an application to a full hydropower system with a coupled turbine. The elastic flow model results are compared against an inelastic model.

13 HYDRO ENERGY↗

DOE OSTI2023

In-Situ Spatial Mapping of Hydrogen in Yttrium Hydrides at LANSCE (FY23 Version, Rev. 1)

This report summarizes the development of neutron imaging capabilities and experimental activities performed at the Los Alamos Neutron Science Center (LANSCE) with the main goal of measuring temperature-driven hydrogen diffusion within bulk-yttrium hydride (YH x ) materials. Yttrium hydride is the leading candidate to serve as a solid neutron moderator in microreactor cores, owing to its high density of hydrogen atoms as well as its superior thermal stability compared to all other metal hydrides. The experimental results and technique developments reported herein support the U.S. Department of Energy Office of Nuclear Energy’s (DOE-NE) Microreactor Program under Technology Maturation. In particular, it addresses the critical need to experimentally validate and verify hydrogen-diffusion models of metal hydrides used in high-temperature microreactor designs by means of high-spatial-resolution neutron imaging. These capabilities were designed to apply large temperature gradients across centimeter-sized YH x pellets to simulate conditions faced in the microreactor environment. In principle, neutron imaging, combined with in-situ sample heating, enables near real-time tracking of hydrogen diffusion in YH x on the sub-millimeter scale. In this report, an overview of neutron imaging methodology and technologies are given in the context of recent spatial measures of hydrogen concentrations in similar metal hydrides. Additionally, the commissioning and operation of a custom-built compact dual-zone furnace is given along with details on three in-situ heating measurements of YH x performed over the 2020 to 2022 LANSCE operation cycles. The aims of these experiments ranged from furnace commissioning, determining sample quality, i.e., hydrogen uniformity via neutron computed tomography, and studying the effects of applied temperature-gradients on YH x pellets. Analyses and results from these neutron imaging measurements are given along with outlooks and guidelines for optimal future hydrogen diffusion measurements. Our conclusions are as follows. Image analyses indicate that centimeter-sized yttrium hydride cylindrical pellets exhibit uniform, whole-body hydrogen desorption and absorption without clear temperature dependence as reflected in the image attenuation at the opposing ends of each sample. This suggests that despite the large magnitude in temperature gradients applied by the furnace heating elements, the sample equilibrates to an unknown intermediate temperature. The origin of this result is likely the combination of short sample length (∼1cm) and use of a TZM can for containment where the latter created a thermal short across the sample. Nevertheless, the results from the most recent measurements indicate that neither significant concentration gradients of hydrogen were formed in centimeter-sized samples through the entire temperature range (25 °C to 950 °C) nor any formed due to temperature gradients on the order of 50 °C/cm up to 700 °C/cm. Furthermore, images from the FY2021 and FY2022 measurements indicate that samples of YH x , fabricated from either the direct hydride or powder metallurgy methods, are highly uniform in their hydrogen concentration to within the measurements’ spatial resolutions. The following questions arise from these latest results: 1) What is the intermediate temperature of the pellets in the TZM cans? 2) How quickly does the temperature equilibrate within the sample? and, 3) Do the observed changes in image attenuation follow known pressure-composition-temperature relations of yttrium hydride?

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

NASA NTRS1978

Olivine cooling speedometers

Several kinetic models of zoning in olivines are discussed at length. The effects on predicted cooling rates of various assumptions used in the analyses are evaluated. It is concluded that the models of Walker et al. (1977) and Taylor et al. (1977) both provide underestimates of the cooling rate required to preserve a given compositional profile, and that both models as well as the model of Taylor et al. (1978) can be used to provide order-of-magnitude estimates of cooling rates. A new model is described which considers diffusion in both solid and liquid during crystallization as well as diffusion in the solid after crystallization is complete. The model provides a description of the compositional gradients which develop during crystallization as well as after cooling at various rates. Applied to olivine crystals nucleated at 1272 C in a high-iron analogue to Lunar Composition 15555, the model predicts only slight compositional gradients - in accord with electron beam microprobe measurements on crystals grown isothermally at this temperature.

Onorato, P. I. K.↗

NASA NTRS1998

The Dissociation Energies of He2, HeH, and ArH; A Bond Function Study

The bond energies and bond lengths are determined for He2, HeH, and ArH at the CCSD(T) level using both atom-centered basis sets and those that include bond functions. The addition of bond functions dramatically improves the rate of convergence of the results with respect to the size of the atom-centered basis set; with bond functions, triple zeta atom-centered basis set, outperform quintuple zeta basis sets without bond functions. The addition of bond functions also reduces the number of diffuse functions that must be added to the atom-centered sets. Employing bond functions appear to offer a very cost effective method of computing the interaction between weakly bound systems, especially for He.

Bauschlicher, Charles W., Jr.↗

DOE OSTI2024

A phase-field diffraction model for thermo-hydro-mechanical propagating fractures

This paper introduces a novel diffraction based thermo-hydraulic–mechanical (THM) model for fracture propagation using a phase-field fracture (PFF) approach. The key innovation of the THM-PFF model lies in its integrated treatment of four solution variables—displacements, phase-field, pressure, and temperature—each governed by a combination of conservation of momentum (mechanics problem), a variational inequality (constrained minimization problem), mass conservation (pressure problem), and energy conservation (temperature problem). This leads to a new formulation of a coupled variational inequality system. A major advancement is the development of an extended fixed-stress algorithm, where displacements, phase-field, pressures, and temperatures are solved in a staggered sequence. An important aspect of this work is the global coupling of pressures and temperatures across the domain using diffraction systems, with diffraction coefficients defined by material parameters weighted by the diffusive phase-field variable. To ensure robust local mass conservation, we employ enriched Galerkin finite elements (EG) for both pressure and temperature diffraction equations. By enriching the continuous Galerkin basis functions with discontinuous piecewise constants, EG accurately represents solution and parameter discontinuities while preserving local mass and energy conservation—crucial aspects for THM problems and realistic behavior. Moreover, the use of a predictor–corrector local mesh adaptivity scheme is employed, allowing the model to handle small phase-field length-scale parameters while maintaining high numerical accuracy and reasonable computational cost. Furthermore, these new model and algorithmic developments represent significant advances in the field and have been substantiated through rigorous numerical tests.

Diffraction systems↗

DOE OSTI2024

Prandtl number effects on extreme mixing events in forced stratified turbulence

Relatively strongly stratified turbulent flows tend to self-organise into a ‘layered anisotropic stratified turbulence’ (LAST) regime, characterised by relatively deep and well-mixed density ‘layers’ separated by relatively thin ‘interfaces’ of enhanced density gradient. Understanding the associated mixing dynamics is a central problem in geophysical fluid dynamics. It is challenging to study LAST mixing, as it is associated with Reynolds numbers $Re := UL/\nu \gg 1$ and Froude numbers $Fr :=(2{\rm \pi} U)/(L N) \ll 1$ ( $U$ and $L$ being characteristic velocity and length scales, $\nu$ the kinematic viscosity and $N$ the buoyancy frequency). Since a sufficiently large dynamic range (largely) unaffected by stratification and viscosity is required, it is also necessary for the buoyancy Reynolds number $Re_{b} := \epsilon /(\nu N^{2}) \gg 1$ , where $\epsilon$ is the (appropriately volume-averaged) turbulent kinetic energy dissipation rate. This requirement is exacerbated for oceanically relevant flows, as the Prandtl number $Pr := \nu /\kappa = {O}(10)$ in thermally stratified water (where $\kappa$ is the thermal diffusivity), thus leading (potentially) to even finer density field structures. We report here on four forced fully resolved direct numerical simulations of stratified turbulence at various Froude ( $Fr=0.5, 2$ ) and Prandtl ( $Pr=1, 7$ ) numbers forced so that $Re_{b}=50$ , with resolutions up to $30\,240 \times 30\,240 \times 3780$ . We find that, as $Pr$ increases, emergent ‘interfaces’ become finer and their contribution to bulk mixing characteristics decreases at the expense of the small-scale density structures populating the well-mixed ‘layers’. However, extreme mixing events (as quantified by significantly elevated local destruction rates of buoyancy variance $\chi _0$ ) are always preferentially found in the (statically stable) interfaces, irrespective of the value of $Pr$ .

Mechanics↗

DOE OSTI2023

Fluctuating hydrodynamics and the Rayleigh–Plateau instability

The Rayleigh–Plateau instability occurs when surface tension makes a fluid column become unstable to small perturbations. At nanometer scales, thermal fluctuations are comparable to interfacial energy densities. Consequently, at these scales, thermal fluctuations play a significant role in the dynamics of the instability. These microscopic effects have previously been investigated numerically using particle-based simulations, such as molecular dynamics (MD), and stochastic partial differential equation–based hydrodynamic models, such as stochastic lubrication theory. In this paper, we present an incompressible fluctuating hydrodynamics model with a diffuse-interface formulation for binary fluid mixtures designed for the study of stochastic interfacial phenomena. An efficient numerical algorithm is outlined and validated in numerical simulations of stable equilibrium interfaces. We present results from simulations of the Rayleigh–Plateau instability for long cylinders pinching into droplets for Ohnesorge numbers of Oh = 0.5 and 5.0. Both stochastic and perturbed deterministic simulations are analyzed and ensemble results show significant differences in the temporal evolution of the minimum radius near pinching. Short cylinders, with lengths less than their circumference, were also investigated. As previously observed in MD simulations, we find that thermal fluctuations cause these to pinch in cases where a perturbed cylinder would be stable deterministically. Finally, we show that the fluctuating hydrodynamics model can be applied to study a broader range of surface tension–driven phenomena.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1982

Status of thermoelectronic laser energy conversion, TELEC

A concept known as a thermo-electronic laser energy converter (TELEC), was studied as a method of converting a 10.6 micron CO2 laser beam into electric power. The calculated characteristics of a TELEC seem to be well matched to the requirements of a spacecraft laser energy conversion system. The TELEC is a high power density plasma device which absorbs an intense laser beam by inverse bremsstrahlung with the plasma electrons. In the TELEC process, electromagnetic radiation is absorbed directly in the plasma electrons producing a high electron temperature. The energetic electrons diffuse out of the plasma striking two electrodes which are in contact with the plasma at the boundaries. These two electrodes have different areas: the larger one is designated as the collector, the smaller one is designated as the emitter. The smaller electrode functions as an electron emitter provide continuity of the current. Waste heat is rejected from the collector electrode. An experiment was carried out with a high power laser using a cesium vapor TELEC cell with 30 cm active length. Laser supported plasma were produced in the TELEC device during a number of laser runs over a period of several days. Electric power from the TELEC was observed with currents in the range of several amperes and output potentials of less than 1 volt.

Britt, E. J.↗

NASA NTRS1994

The study of dopant segregation behavior during the growth of GaAs in microgravity

An investigation into the segregation behavior of selenium doped gallium arsenide during directional solidification in the microgravity environment was conducted using the Crystal Growth Furnace (CGF) aboard the first United States Microgravity Laboratory (USML-1). The two crystals grown were 1.5 cm in diameter and 16.5 cm in length with an initial melt length of 14 cm. Two translation periods were executed, the first at 2.5 microns/s and after a specified time, which was different between the two experiments, the translation rate was doubled to 5.0 microns/s. The translation was then stopped and the remaining sample melt was solidified using a gradient freeze technique in the first sample and a rapid solidification in the second experiment. Measurement of the selenium dopant distribution, using quantitative infrared transmission imaging, indicates that the first sample initially achieved diffusion controlled growth as desired. However, after about 1 cm of growth, the segregation behavior was driven from a diffusion controlled growth regime to a complete mixing regime. Measurements in the second flight sample indicated that the growth was always in a complete mixing regime. In both experiments, voids in the center line of the crystal, indicative of bubble entrapment, were found to correlate with the position in the crystal when the translation rates were doubled.

Matthiesen, David H.↗

NASA NTRS1993

A diffusion-limited aggregation model for the evolution of drainage networks

We propose a modified diffusion-limited aggregation (DLA) model for the evolution of fluvial drainage networks. Random walkers are introduced randomly on a grid, and each two-dimensional random walk proceeds until the walker finds a drainage network on which to accrete. This model for headward growth of drainage networks generates drainage patterns remarkably similar to actual drainages. The model also predicts statistical features which agree with actual networks, including the frequency-order (bifurcation) ratio (R(sub b) = 3.98) and the stream length-order (R(sub r) = 2.09). Using the definition of network fractal dimension D = log R(sub b)/log R(sub r), we find that our DLA model gives D = 1.87, near the observed range of D approximately equal to 1.80 - 1.85.

Masek, Jeffrey G.↗

NASA NTRS1997

Evolution of Structure in the Intergalactic Medium and the Nature of the LY-Alpha Forest

We have performed a detailed statistical study of the evolution of structure in a photoionized intergalactic medium (IGM) using analytical simulations to extend the calculation into the mildly nonlinear density regime found to prevail at z = 3. Our work is based on a simple fundamental conjecture: that the probability distribution function of the density of baryonic diffuse matter in the universe is described by a lognormal (LN) random field. The LN distribution has several attractive features and follows plausibly from the assumption of initial linear Gaussian density and velocity fluctuations at arbitrarily early times. Starting with a suitably normalized power spectrum of primordial fluctuations in a universe dominated by cold dark matter (CDM), we compute the behavior of the baryonic matter, which moves slowly toward minima in the dark matter potential on scales larger than the Jeans length. We have computed two models that succeed in matching observations. One is a nonstandard CDM model with OMEGA = 1, h = 0.5, and GAMMA = 0.3, and the other is a low-density flat model with a cosmological constant (LCDM), with OMEGA = 0.4, OMEGA(sub LAMBDA) = 0.6, and h = 0.65. In both models, the variance of the density distribution function grows with time, reaching unity at about z = 4, where the simulation yields spectra that closely resemble the Ly-alpha forest absorption seen in the spectra of high-z quasars. The calculations also successfully predict the observed properties of the Ly-alpha forest clouds and their evolution from z = 4 down to at least z = 2, assuming a constant intensity for the metagalactic UV background over this redshift range. However, in our model the forest is not due to discrete clouds, but rather to fluctuations in a continuous intergalactic medium. At z = 3; typical clouds with measured neutral hydrogen column densities N(sub H I) = 10(exp 13.3), 10(exp 13.5), and 10(exp 11.5) /sq cm correspond to fluctuations with mean total densities approximately 10, 1, and 0.1 times the universal mean baryon density. Perhaps surprisingly, fluctuations whose amplitudes are less than or equal to the mean density still appear as "clouds" because in our model more than 70% of the volume of the IGM at z = 3 is filled with gas at densities below the mean value.

Bi, Hongguang↗

NASA NTRS2006

Numerical Modeling of Flow Control in a Boundary-Layer-Ingesting Offset Inlet Diffuser at Transonic Mach Numbers

This paper will investigate the validation of a NASA developed, Reynolds-averaged Navier-Stokes (RANS) flow solver, OVERFLOW, for a boundary-layer-ingesting (BLI) offset (S-shaped) inlet in transonic flow with passive and active flow control devices as well as the baseline case. Numerical simulations are compared to wind tunnel results of a BLI inlet conducted at the NASA Langley 0.3-Meter Transonic Cryogenic Tunnel. Comparisons of inlet flow distortion, pressure recovery, and inlet wall pressures are performed. The numerical simulations are compared to the BLI inlet data at a freestream Mach number of 0.85 and a Reynolds number of approximately 2 million based on the length of the fan-face diameter. The numerical simulations with and without wind tunnel walls are performed, quantifying effects of the tunnel walls on the BLI inlet flow measurements. The wind tunnel test evaluated several different combinations of jet locations and mass flow rates as well as a vortex generator (VG) vane case. The numerical simulations will be performed on a single jet configuration for varying actuator mass flow rates at a fix inlet mass flow condition. Validation of the numerical simulations for the VG vane case will also be performed for varying inlet mass flow rates. Overall, the numerical simulations were able to predict the baseline circumferential flow distortion, DPCPavg, very well for comparisons made within the designed operating range of the BLI inlet. However the CFD simulations did predict a total pressure recovery that was 0.01 lower than the experiment. Numerical simulations of the baseline inlet flow also showed good agreement with the experimental inlet centerline surface pressures. The vane case showed that the CFD predicted the correct trends in the circumferential distortion for varying inlet mass flow but had a distortion level that was nearly twice as large as the experiment. Comparison to circumferential distortion measurements for a 15 deg clocked 40 probe rake indicated that the circumferential distortion levels are very sensitive to the symmetry of the flow and that a miss alignment of the vanes in the experiment could have resulted in this difference. The numerical simulations of the BLI inlet with jets showed good agreement with the circumferential inlet distortion levels for a range of jet actuator mass flow ratios at a fixed inlet mass flow rate. The CFD simulations for the jet case also predicted an average total pressure recovery that was 0.01 lower than the experiment as was seen in the baseline. Comparison of the flow features the jet case revealed that the CFD predicted a much larger vortex at the engine fan-face when compare to the experiment.

Allan Brian G.↗

DOE OSTI2025

Integral X-ray diffuse scattering for studying irradiation-induced dislocation loops in single crystals

Numerical diffuse scattering cross-section calculations are used to establish a rigorous basis for determining the concentration and size distribution of dislocation loops in irradiated single crystals from integral X-ray diffuse scattering (XRDS) measurements. Differential XRDS intensities for prismatic {111} type dislocation loops are numerically calculated as a function of loop radius R and wavevectors q relative to Bragg reflections in tungsten. The results show the well known 1/q 2 Huang scattering form at small q that transitions to a ∼1/q 4 dependence associated with the Stokes–Wilson approximation for q ≳ 1/R. More importantly, they show further that the 1/q 4 falloff is not the asymptotic large-q form of the diffuse scattering for small loops (R < 200 Å) as has often been assumed. Rather, for loop sizes as small as R ≃ 5 Å with strong curvature, the calculations show definitively that the scattering transitions to a robust 1/q 5 falloff at larger q that arises due to the local strains near the dislocation core defining the circumference of the dislocation loops. The presence of this 1/q 5 asymptotic form for both small and large loops is experimentally confirmed using an integral XRDS measurement around the 110 reflection on self-ion-irradiated tungsten combined with numerically calculated integral XRDS cross-sections. Accordingly, the historical two-region theoretical treatment of the cross-sections for integral XRDS is extended to a three-region model that has direct sensitivity to the (first-moment) dislocation line lengths of dislocation loops. In conclusion, these developments enable the use of both numerical and analytically modeled cross-sections to make accurate integral XRDS determinations of dislocation loop sizes and concentrations using modest-intensity laboratory X-ray sources.

42 ENGINEERING↗

DOE OSTI2023

High-speed mapping of surface charge dynamics using sparse scanning Kelvin probe force microscopy

Unraveling local dynamic charge processes is vital for progress in diverse fields, from microelectronics to energy storage. This relies on the ability to map charge carrier motion across multiple length- and timescales and understanding how these processes interact with the inherent material heterogeneities. Towards addressing this challenge, we introduce high-speed sparse scanning Kelvin probe force microscopy, which combines sparse scanning and image reconstruction. This approach is shown to enable sub-second imaging (>3 frames per second) of nanoscale charge dynamics, representing several orders of magnitude improvement over traditional Kelvin probe force microscopy imaging rates. Bridging this improved spatiotemporal resolution with macroscale device measurements, we successfully visualize electrochemically mediated diffusion of mobile surface ions on a LaAlO 3 /SrTiO 3 planar device. Such processes are known to impact band-alignment and charge-transfer dynamics at these heterointerfaces. Furthermore, we monitor the diffusion of oxygen vacancies at the single grain level in polycrystalline TiO 2 . Through temperature-dependent measurements, we identify a charge diffusion activation energy of 0.18 eV, in good agreement with previously reported values and confirmed by DFT calculations. Together, these findings highlight the effectiveness and versatility of our method in understanding ionic charge carrier motion in microelectronics or nanoscale material systems.

25 ENERGY STORAGE↗

DOE OSTI2023

Dynamic Subgrid Turbulence Modeling for Shallow Cumulus Convection Simulations beyond LES Resolutions

A scale-dependent dynamic Smagorinsky model is implemented in the Met Office/NERC Cloud (MONC) model using two averaging flavors, along Lagrangian pathlines and local moving averages. The dynamic approaches were compared against the conventional Smagorinsky–Lilly scheme in simulating the diurnal cycle of shallow cumulus convection. The simulations spanned from the LES to the near-gray-zone and gray-zone resolutions and revealed the adaptability of the dynamic model across the scales and different stability regimes. The dynamic model can produce a scale- and stability-dependent profile of the subfilter turbulence length scale across the chosen resolution range. At gray-zone resolutions the adaptive length scales can better represent the early precloud boundary layer leading to temperature and moisture profiles closer to the LES compared to the standard Smagorinsky. As a result, the initialization and general representation of the cloud field in the dynamic model is in good agreement with the LES. In contrast, the standard Smagorinsky produces a less well-mixed boundary layer, which fails to ventilate moisture from the boundary layer, resulting in the delayed spinup of the cloud layer. Moreover, strong downgradient diffusion controls the turbulent transport of scalars in the cloud layer. However, the dynamic approaches rely on the resolved field to account for nonlocal transports, leading to overenergetic structures when the boundary layer is fully developed and the Lagrangian model is used. Introducing the local averaging version of the model or adopting a new Lagrangian time scale provides stronger dissipation without significantly affecting model behavior.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1979

A cesium TELEC experiment at Lewis Research Center

The thermoelectronic laser energy converter (TELEC), was studied as a method of converting a 10.6 mm CO2 laser beam into electric power. The calculated characteristics of a TELEC seem to be well matched to the requirements of a spacecraft laser energy conversion system. The TELEC is a high power density plasma device which absorbs an intense laser beam by inverse bremsstrahlung with the plasma electrons. In the TELEC process, electromagnetic radiation is absorbed directly in the plasma electrons producing a high electron temperature. The energetic electrons diffuse out of the plasma striking two electrodes which are in contact with the plasma at the boundaries. These two electrodes have different areas: the larger one is designated as the collector, the smaller one is designated as the emitter. The smaller electrode functions as an electron emitter to provide continuity of the current. Waste heat is rejected from the collector electrode. An experiment was carried out with a high power laser using a cesium vapor TELEC cell with 30 cm active length. Laser supported plasma was produced in the TELEC device during a number of laser runs over a period of several days. Electric power from the TELEC was observed with currents in the range of several amperes and output potentials of less than 1 volt. The magnitudes of these electric outputs were smaller than anticipated but consistent with the power levels of the laser during this experiment.

Britt, E. J.↗