Search NASA⌕ Search

SEARCH · Search NASA

Results for “HEAT EQUATION”

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

Challenging Common Assumptions of Thick-Wall Chamber Dynamics in Inertial Fusion Systems Using MOOSE-based Multiphysics Simulations

As an increasing number of companies look toward commercial Inertial Fusion Energy (IFE) designs, there is a pressing need to understand the physics of thick-wall chamber gas dynamics. The thick liquid wall approach implements a renewable wall to mitigate the fusion target emissions, thereby reducing the radiation damage rate and significantly extending the lifetime of chamber structures, leading to increased plant availability and reduced waste streams in comparison to dry wall chamber designs. It is necessary, however, to assess the critical performance and safety aspects of these systems. For example, it is crucial to predict (1) where the mass ablated from the liquid walls will vent, which determines the placement of condensing surfaces; (2) debris propagation up the beam lines, which provides essential information for design and protection requirements; (3) peak pressures and impulse on chamber walls, which affect chamber structural design; and (4) momentum transfer to the liquid jets, which constrains the shape and positioning of the jets. In turn, the chamber design and its liquid walls affect shielding requirements, material activation, and tritium fuel cycle. Currently available simulation tools, however, are unable to accurately capture key thick-wall chamber dynamics. Significant assumptions are often made to simplify the system and reduce computational cost and modeling capability needs, but the impact of these assumptions on simulation predictions has not been evaluated. For example, no three-dimensional simulations can be found in the open literature to evaluate gas venting and momentum transfer to the jets with simulations using two-dimensional domains to represent complex three-dimensional geometries. Moreover, limited studies have been dedicated to jet breakup due to both turbulence and neutron heating, and no studies have been found that evaluate how jet breakup can impact shock-jet interaction. Furthermore, effects of radiative heat transfer have rarely been included for the hydrodynamic phase of shock propagation, and integration of proper equations of state in shock dynamics codes has been mostly exploratory. In this study, we use the flexible, high-fidelity Multiphysics Object-Oriented Simulation Environment (MOOSE) to model these complex phenomena and inform design and safety studies. Capabilities to model thick-wall chamber gas dynamics are being developed, and the impact of the assumptions listed above (i.e., two-dimensional vs three-dimensional, absence of jet breakout, no radiative heat transfer, and ideal gas behavior) are being quantified.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

TNet: A Model-Constrained Tikhonov Network Approach for Inverse Problems

Deep learning (DL), in particular deep neural networks, by default is purely data-driven and in general does not require physics. This is the strength of DL but also one of its key limitations when applied to science and engineering problems in which underlying physical properties—such as stability, conservation, and positivity—and accuracy are required. DL methods in their original forms are often not capable of respecting the underlying mathematical models or achieving desired accuracy even in big-data regimes. On the other hand, many data-driven science and engineering problems, such as inverse problems, typically have limited experimental or observational data, and DL would overfit the data in this case. Leveraging information encoded in the underlying mathematical models, we argue, not only compensates for missing information in low data regimes but also provides opportunities to equip DL methods with the underlying physics, hence promoting better generalization. This paper develops a model-constrained DL approach and its variant TNet—a Tikhonov neural network—which are capable of learning not only information hidden in the training data but also in the underlying mathematical models to solve inverse problems governed by partial differential equations in low data regimes. We provide the constructions and some theoretical results for the proposed approaches for both linear and nonlinear inverse problems. Since TNet is designed to learn inverse solutions with Tikhonov regularization, it is interpretable: in fact it recovers Tikhonov solutions for linear cases while potentially approximating Tikhonov solutions for nonlinear inverse problems. We also prove that data randomization can enhance not only the smoothness of the networks but also their generalizations. Comprehensive numerical results confirm the theoretical findings and show that with even as little as 1 training data sample for one-dimensional (1D) deconvolution, 5 for an inverse 2D heat conductivity problem, 100 for inverse initial conditions for a time-dependent 2D Burgers’s equation, and 50 for inverse initial conditions for 2D Navier–Stokes equations, TNet solutions can be as accurate as Tikhonov solutions while being several orders of magnitude faster. Furthermore, this is possible owing to the model-constrained term, replications, and randomization.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Reversible hydrogen storage in multilayer graphane: Lattice dynamics, compressibility, and heat capacity studies

Multilayer graphane (hydride of graphite) is a crystalline hydrocarbon of composition CH, which can be synthesized from graphite and molecular hydrogen at pressures above 2GPa [V.E. Antonov et al. Carbon 100 (2016) 465]. Using X-ray diffraction, this compound was tentatively identified as the “graphane II” phase of 3D-graphane predicted by ab initio calculations [X.-D. Wen et al. PNAS 108 (2011) 6833] and consisting of layers of 2D-graphane in the “chair” conformation. When heated in a vacuum, the compound does not form any intermediate hydrocarbons and reversibly decomposes back into graphite and hydrogen at 770–920 K. In the present work, almost single-phase samples of graphite hydride and deuteride were synthesized at 7.4 GPa and 870 K. Their investigation by inelastic neutron scattering supplemented by ab initio calculations gave spectra g(E) of the phonon density of states with a gap of about 15 meV at approx. 100 meV, which is a unique identifier for the chair form of graphane. The equation of state V(P) of the hydride was studied at room temperature and hydrogen pressures up to 53 GPa by synchrotron X-ray diffraction in a diamond anvil cell. Further, the graphane II phase did not react with the surrounding hydrogen and did not undergo any phase transformations upon the compression and after heating to 1500 K at 53GPa. The high thermal and pressure stability of this exotic phase makes it an important part of the C–H system. The obtained g(E) spectra of graphite hydride and deuteride were used to calculate temperature dependences of their heat capacity. Measurements of the heat capacity at temperatures 120–673 K confirmed the good accuracy of these calculations.

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

Mass-conserving implicit–explicit methods for coupled compressible Navier–Stokes equations

Earth system models are composed of coupled components that separately model systems such as the global atmosphere, ocean, and land surface. While these components are well developed, coupling them in a single system can be a significant challenge. Computational efficiency, accuracy, and stability are principal concerns. In this study we focus on these issues. In particular, implicit–explicit (IMEX) tight and loose coupling strategies are explored for handling different time scales. For a simplified model for the air–sea interaction problem, we consider coupled compressible Navier–Stokes equations with an interface condition. Under the rigid-lid assumption, horizontal momentum and heat flux are exchanged through the interface. Several numerical experiments are presented to demonstrate the stability of the coupling schemes. Here, we show both numerically and theoretically that our IMEX coupling methods are mass conservative for a coupled compressible Navier–Stokes system with the rigid-lid condition

42 ENGINEERING↗
DOE OSTI2023

Monsoon Depression Amplification by Horizontal Shear and Humidity Gradients: A Shallow Water Perspective

Transient, synoptic-scale vortices produce a large fraction of total rainfall in most monsoon regions and are often associated with extreme precipitation. However, the mechanism of their amplification remains a topic of active research. For monsoon depressions, which are the most prominent synoptic-scale vortex in the Asian–Australian monsoon, recent work has suggested that meridional gradients in zonal wind in the vortex environment may produce growth through barotropic instability, while meridional gradients in environmental humidity have also been proposed to cause amplification through coupling with precipitating convection. Here, a two-dimensional shallow water model on a sphere with parameterized precipitation is used to examine the relative role played by these two environmental gradients. By systematically varying the meridional moisture gradient and meridional wind shear for both weak, quasi-linear waves and finite-amplitude isolated vortices, we show that rotational winds in the initial vortex are amplified most strongly by meridional shear of the environmental zonal wind, while vortex precipitation rates are most sensitive to environmental moisture gradients. The growth rate in the presence of both gradients is less than the sum of growth rates in the presence of isolated gradients, as the phase relation between moisture and vorticity anomalies becomes distorted with increasing shear. In conclusion, these results suggest that background meridional gradients in both zonal wind and environmental humidity can contribute to the amplification of vortices to monsoon depression strength, but with some degree of decoupling of the dry rotational flow and the moist convection.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

System Design of a 2.0 MWth Sodium/Molten Salt Pilot System

Nitrate molten salt concentrating solar power (CSP) systems are currently deployed globally and are considered state-of the art heat transfer fluids (HTFs) for present day high-temperature operation. Although slightly higher limits may be possible with molten salt, to fully realize SunShot efficiency goals of $15/kWhth HTFs and an LCOE of 6¢/kWh, HTF technologies working at higher temperatures (e.g., 650 °C to 750 °C) will require an alternative to molten salts, such as with alkali metal systems. This investigation explores the development of a 2.0 MWth sodium receiver system that employs a sodium receiver as the HTF, as well as with a ternary chloride (20%NaCl/40%MgCl/40%KCl by mol wt.%) salt as a thermal energy storage (TES) medium to facilitate a 6-hr. storage duration. A sodium-to-salt heat exchanger model as well as a salt-to-sCO2 primary heat exchanger model are employed and evaluated in this investigation. A thermodynamic system design model was developed using Engineering Equation Solver (EES) where state properties were calculated at inlets and outlets along both hot and cold legs of the pilot-scale plant. This investigation assesses receiver performance as well as system efficiency studies for the pump and system operational ranges. Results found that high efficiency sodium receivers were found to have higher heat transfer coefficients and required far less spreading of incident flux. The system performance model results suggest that for a pump speed of 2400 RPM, respective hot and cold pump TDH values were determined to be 260.1–307 ft. and 260.1–307 ft for pump flow rates of 90–120 GPM.

41 EE - Solar Energy Technologies Office (EE-4S)↗
DOE OSTI2025

Equipment List Comparing Balance of Plant Containing a Heat Pump against a Reference Electricity Generating Plant

Approximately two-thirds of U.S. energy consumption in the industrial and transportation sectors relies on fossil fuels. These sectors require high-quality heat, i.e., thermal energy at very high temperatures, for molecular transformation processes. The Integrated Energy Systems (IES) program aims to assess the economic potential of utilizing nuclear-grade heat from Advanced Reactors (ARs) to meet the high-quality heat demands. By having industrial processes (IPs) supplied with nuclear-generated heat, manufacturers could benefit from more stable and potentially lower energy costs, reducing reliance on volatile fossil fuel markets. The main outcome of the FY24 research was the thermodynamic assessment and gap analysis of steam generation for IP applications. The study completed in June 2024 demonstrated that the required steam temperatures could be achieved by integrating a heat pump into an AR power plant. After identifying the thermal demands of target IPs, multiple balance of plant configurations for the Xe-100 reactor by X-energy, incorporating a heat pump, were analyzed. Their technical feasibilities were assessed, including the design of suitable axial compressors for these applications. Comparative performance analysis showed that thermal efficiency alone is insufficient to evaluate system suitability. To address this, a new indicator (heat factor) was introduced in the report released in September 2024 to quantify the low-quality thermal power needed to produce one unit of high-quality heat for the IP. Results showed that integrating heat pumps into Rankine cycles enables higher steam temperatures, though at the expense of increased thermal energy input. This report builds upon and completes the foundational work previously undertaken. It focuses on the design of two Balance of Plant (BOP) configurations, both based on Rankine energy conversion cycles: “Case 1”, which involves electricity generation only, and “Case 2”, which combines electricity and high-temperature heat generation for industrial use. For each configuration, a comprehensive equipment list was developed, detailing all major components such as turbines, compressors, heat exchangers, pumps, and control systems. These lists will serve as the basis for future comparative cost analyses, with the goal of assessing the number and type of components required to integrate a heat pump into the Rankine cycle and to establish a heat transport system capable of delivering thermal energy from the nuclear plant to an industrial facility. Using the constitutive equations presented in the June 2024 milestone, the operating conditions of all BOP components for both “Case 1” and “Case 2” were evaluated. These parameters—such as temperature, pressure, mass flow rate, and steam quality—served as the basis for calculating the associated thermal and mechanical power flows. The net power required from the heat pump to raise the steam temperature to the target level was also determined. The thermodynamic performance of the configuration was then assessed using the heat factor metric. The key outcome of this analysis is a comparative table that presents the equipment that was used in the “Case 1” and “Case 2” configurations. This study offers a preliminary comparison of the two designs, providing insight into the impact of integrating a heat pump in terms of component requirements and thermal efficiency. This equipment list, together with the evaluated operating conditions, also serves as a foundation for the economic analyses scheduled for the current fiscal year.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

A Probabilistic Scheme for Semilinear Nonlocal Diffusion Equations with Volume Constraints

This work presents a probabilistic scheme for solving semilinear nonlocal diffusion equations with volume constraints and integrable kernels. The nonlocal model of interest is defined by a time-dependent semilinear partial integro-differential equation (PIDE), in which the integro-differential operator consists of both local convection-diffusion and nonlocal diffusion operators. Here, our numerical scheme is based on the direct approximation of the nonlinear Feynman–Kac formula that establishes a link between nonlinear PIDEs and stochastic differential equations. The exploitation of the Feynman–Kac representation avoids solving dense linear systems arising from nonlocal operators. Compared with existing stochastic approaches, our method can achieve first-order convergence after balancing the temporal and spatial discretization errors, which is a significant improvement of existing probabilistic/stochastic methods for nonlocal diffusion problems. Error analysis of our numerical scheme is established. The effectiveness of our approach is shown in two numerical examples. The first example considers a three-dimensional nonlocal diffusion equation to numerically verify the error analysis results. The second example presents a physics problem motivated by the study of heat transport in magnetically confined fusion plasmas.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Kinetic simulation of a 50 mTorr capacitively coupled argon discharge over a range of frequencies and comparison to experiments

The effect of driving frequency in the range of 13.56–73 MHz on electron energy distribution and electron heating modes in a 50 mTorr capacitively coupled argon plasma discharge is studied using 1D-3V particle-in-cell simulations. Calculated electron energy probability functions exhibit three distinct “temperatures” for low-, mid-, and high-energy electrons at all the studied driving frequencies. When compared to published experimental data, the calculated probability functions show a reasonable agreement for the energy range resolved in the measurements (about 2–10 eV). Discrepancies due to limitations in experimental energy resolution outside this range lead to differences between computational and experimental values of the electron number density determined from the distribution functions, and the predicted effective electron temperature is within 25% of experimental values. The impedance of the discharge is interpreted in terms of a homogeneous equivalent circuit model, and the driving frequency dependence of the inferred combined sheath thickness is found to obey a known, theoretically derived, power law. The average power transferred from the field to the electrons (electron heating) is computed, and a region of negative heating near the sheath edge, particularly at higher driving frequencies, is identified. Analysis of the electron momentum equation shows that electron inertia, which on temporal averaging would be zero in a linear regime, is responsible for negative values of power deposition near the sheath edge at high driving frequencies due to the highly nonlinear behavior of the discharge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Moment-Fourier approach to ion parallel fluid closures and transport for a toroidally confined plasma

A general method of solving the drift kinetic equation is developed for an axisymmetric magnetic field. Expanding a distribution function in general moments, a set of ordinary differential equations is obtained. Successively expanding the moments and magnetic-field involved quantities in Fourier series, a set of linear algebraic equations is obtained. The set of full (Maxwellian and non-Maxwellian) moment equations is solved to express the first-order density, temperature, and flow velocity in terms of radial gradients of the zeroth-order pressure and temperature. Closure relations that connect parallel heat flux density and viscosity to the radial gradients and parallel gradients of temperature and flow velocity are also obtained by solving the non-Maxwellian moment equations. The closure relations combined with the linearized fluid equations reproduce the same solution obtained directly from the full moment equations. Furthermore, the method can be generalized to derive closures and transport for an electron-ion plasma and a multi-ion plasma in a general magnetic field.

neoclassical transport↗
DOE OSTI2022

Dimensionless parameters for cloudy Rayleigh-Bénard convection: Supersaturation, Damköhler, and Nusselt numbers

In steady-state Rayleigh-Bénard convection, heat is transported by turbulent thermal convection from the bottom, hot surface to the top, cold surface, leading to a height-independent sensible heat flux. When water vapor is present and cloud formation occurs, there is also an additional latent heat flux. Heat transport in cloudy Rayleigh-Bénard convection depends on turbulent flow as well as the microphysical state of the clouds: specifically, whether substantial supersaturations exist and whether cloud liquid water is removed through sedimentation/precipitation. In this article we bridge between the Rayleigh-Bénard convection literature and the atmospheric literature. We express the governing equations for cloudy convection in dimensionless form, thereby explicitly identifying the governing parameters relevant to the cloudy case, including Schmidt, Damköhler, supersaturation, and sedimentation numbers. We further connect to the atmospheric literature by obtaining a Nusselt number (dimensionless heat flux) for a cloud-convection system, directly from the conservation equations for temperature and water vapor. This flux has the same form as that identified by Zhang et al. [L. Zhang, K. L. Chong, and K.-Q. Xia, J. Fluid Mech. 874, 1041 (2019)] for convection with water vapor, but is extended to the cloudy case. For equal thermal and water vapor diffusivities, the flux corresponds to the widely used atmospheric quantities equivalent temperature and moist static energy. Using large eddy simulation (LES) of an idealized cloudy Rayleigh-Bénard convection system with fixed boundary conditions, we find that the equivalent heat flux (Nusselt number) is only weakly dependent on the microphysical details of the system, such as liquid water mixing ratio and cloud droplet number concentration. Finally, from the results, we show the vertical profiles of sensible and latent heat fluxes depend on the liquid water content, whereas the equivalent heat flux remains a constant throughout the height of the chamber.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Quantitative assessment of eddy viscosity rans models for turbulent mixed convection in a differentially heated plane channel

Turbulent mixed convection between two vertical, infinite parallel plates at different temperatures is studied using various two-equation turbulence models. The numerical simulations are performed at a turbulent Reynolds number of Re τ = 150 and a Grashof number of Gr = 9.6 × 10 5 . Comparisons are made against the highly trusted DNS results. Consistent with the DNS approach, the current simulations are performed using constant properties and the Boussinesq approximation to predict the influence of buoyancy. Previous studies have provided assessments of two-equation turbulence models for various scenarios, but often rely on a qualitative “eye” test in order to determine the most appropriate model to predict a given flow. This study aims to provide a new form of quantitative assessment that accounts for both the physics captured by the turbulence model as well as the magnitude of the system response quantities (SRQ) using a modified symmetric mean absolute percent error (SMAPE) method. This method is designed to be approachable to researchers at any level and can be applied to system response quantities from multiple research fields. Uncertainty quantification is also performed to determine the discretization error for each turbulence model. Recommendations are made as to which turbulence models best capture the physics – hydrodynamically and thermally – using both local and global validation metrics. Lastly, a sensitivity analysis is performed on the damping functions used in the most accurate models. This underpins the potential of model developments and adjustments most worth pursuing for buoyant flows. Finally, this framework provides a more physics-based comparative analysis of the selected turbulence models.

42 ENGINEERING↗
DOE OSTI2021

Thermochemical reduction modeling in a high-temperature moving-bed reactor for energy storage: 1D model

The design of robust and efficient high-temperature thermochemical reactors and determination of operating conditions are critical steps toward enabling high-efficiency long-duration solar energy storage. This work presents a computational model for the thermal reduction of a metal oxide material (Mg-Mn-O) up to 1450 °C and the coupled complex transport phenomena in a novel tubular thermal reactor design that features the capability for a high extent-of-reduction (high energy storage density) and inherent heat recuperation. A one-dimensional model coupling counter-current gas–solid flow, two-phase heat transfer, thermochemical redox reactions, and species transport in a moving-bed reactor is developed. Simplified versions of the model are validated with published results in the literature for packed beds with both inert and reactive particles; the fully coupled model is also validated with experimental measurements of a moving-bed reactor in terms of local temperatures and oxygen release at the exit. Detailed comparisons on the effects of different boundary conditions in the reaction zone (prescribed wall temperature vs. heat flux conditions) and formulations based on a simple uniform flow assumption vs. plug flow using Ergun equation for gas flow are investigated. The results are compared with experimental measurements, and for all cases, the energy flow components in the reactor system and the thermal to chemical conversion efficiency and overall system efficiency are computed. Finally, the predicted high thermal-to-chemical efficiency ~95% and system efficiency ~30% agree with experimental measurements.

25 ENERGY STORAGE↗
DOE OSTI2022

ThAu 2 – , ThAu 2 O – , and ThAuOH – anions: Photoelectron spectroscopic and theoretical characterization

The thorium–gold negative ions ThAu 2 – , ThAu 2 O – , and ThAuOH – have been observed and experimentally characterized by anion photoelectron spectroscopy. These experiments are accompanied by extensive ab initio electronic structure calculations using a relativistic composite methodology based primarily on coupled cluster singles and doubles with perturbative triples calculations. The theoretical electron affinities (EAs) at 0 K agree with the experimental adiabatic EAs to within 0.02 eV for all species. Two separate isomers were located in the calculations for ThAuOH – , and detachment from both of these appears to be present in the photoelectron spectrum. Excited electronic states of the neutral molecules are reported at the equation of motion-coupled cluster singles and doubles level of theory. Atomization energies and heats of formation are also calculated for each neutral species and have expected uncertainties of 3 and 4 kcal/mol, respectively. The σ bonds between Th and Au are determined by natural bond orbital analysis to consist of predominately sd hybrids on Th bonding with the Au 6s orbital. In order to investigate the correspondence between the bonding in Th–Au and Th–F molecules, a limited number of calculations were also carried out on most of the F-analogs of this study. These results demonstrate that Au does behave like F in these cases, although the Th–F σ bonds are much more ionic compared to Th–Au. This results in an EA for ThF 2 that is 10 kcal/mol smaller than that of ThAu 2 . The EA values for the Th(IV) species, i.e., ThX 2 O and ThXOH, only differed, however, by 3–4 kcal/mol.

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

Large-scale molecular dynamics simulations of bubble collapse in water: Effects of system size, water model, and nitrogen

Molecular dynamics simulations in the microcanonical ensemble are performed to study the collapse of a bubble in liquid water using the single-site mW and the four-site TIP4P/2005 water models. To study system size effects, simulations for pure water systems are performed using periodically replicated simulation boxes with linear dimensions, L, ranging from 32 to 512 nm with the largest systems containing 8.7 × 10 6 and 4.5 × 10 9 molecules for the TIP4P/2005 and mW water models, respectively. The computationally more efficient mW water model allows us to reach converging behavior when the bubble dynamics results are plotted in reduced units, and the limiting behavior can be obtained through linear extrapolation in L –1 . Qualitative differences are observed between simulations with the mW and TIP4P/2005 water models, but they can be explained by the models’ differences in predicted viscosity and surface tension. Although bubble collapse occurs on time scales of only hundreds of picoseconds, the system sizes used here are sufficiently large to obtain bubble dynamics consistent with the Rayleigh–Plesset equation when using the models’ thermophysical properties as input. For the conditions explored here, extreme heating of the interfacial water molecules near the time of collapse is observed for the larger mW water systems (but the model underpredicts the viscosity), whereas heating is less pronounced for the TIP4P/2005 water systems because its larger viscosity contribution slows the collapse dynamics. The presence of nitrogen within the bubble only starts to affect bubble dynamics near the very end of the initial collapse, leading to an incomplete collapse and strong rebound for the mW water model. Although nitrogen is non-condensable at 300 K, it becomes highly compressed and reaches a liquid-like density near the collapse point. We find that the dissolution of nitrogen is much slower than the movement of the collapsing water front, and the re-expansion of the dense nitrogen droplet gives rise to bubble rebound. The incompatibility of the collapse and dissolution time scales should be considered for continuum-scale modeling of bubble dynamics. Finally, we also confirm that the diffusion coefficient for dissolved nitrogen is insensitive to pressure as the liquid transitions from a compressed to a stretched state.

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

Bohm Criterion of Plasma Sheaths away from Asymptotic Limits

The plasma exit flow speed at the sheath entrance is constrained by the Bohm criterion. The so-called Bohm speed regulates the plasma particle and power exhaust fluxes to the wall, and it is commonly deployed as a boundary condition to exclude the sheath region in quasineutral plasma modeling. In this work, the Bohm criterion analysis is performed in the intermediate plasma regime away from the previously known limiting cases of adiabatic laws and the asymptotic limit of infinitesimal Debye length in a finite-size system, using the transport equations of an anisotropic plasma. The resulting Bohm speed has explicit dependence on local plasma heat flux, temperature isotropization, and thermal force. Comparison with kinetic simulations demonstrates its accuracy over the plasma-sheath transition region in which quasineutrality is weakly perturbed and the Bohm criterion applies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

AEROICESTUDY-Colorado State University Ice Spectrometer

The dataset comprises measures, using Colorado State University's Ice Spectrometer (IS, an immersion freezing device with a range from 0 °C down to -26 to -29 °C) of atmospheric ice nucleating particle (INP) concentrations taken at ARM's Southern Great Plains (SGP) Site Guest Instrument Facility (GIF), in Lamont, Oklahoma, during the Aerosol-Ice Formation Closure Pilot Study (AEROICESTUDY). A filter sampler was mounted approximately 6 m above ground level atop the GIF rooftop platform. Single-use filter units open to the atmosphere fitted with pre-cleaned and pre-sterilized, 47-mm diameter Nuclepore polycarbonate filters were used. Filters were typically drawn for either a full day period (9AM-5PM CDT) or broken up into morning (9AM-12PM CDT) and afternoon (1PM-5PM CDT) periods. Additionally, 4 overnight filters were taken (6PM-8AM CDT). The volume of air filtered was typically between 3000 and 8000 standard liters for daily samples. After collection, filters were stored at -20 °C until processed. Initial processing to obtain spectra of INP number concentration active via the immersion freezing mechanism versus temperature was conducted using CSU's IS instrument (McCluskey et al., 2018). For measurement of INPs, collected aerosol particles were re-suspended in 8 mL of 0.1 µm-filtered deionized (DI) water. Aliquots of each suspension, and serial dilutions, were dispensed into trays which were fit into aluminum blocks in the IS. Samples were cooled at 0.33 °C per min and the freezing temperatures of wells recorded automatically. Cumulative INP concentrations were determined by first calculating the INPs per mL of suspension based on Vali (1971) and then converting to concentration per standard liter of air using the proportion of the total liquid sample dispensed and the air sample volumes. Two-tailed, 95% confidence intervals for binomial sampling were based on Equation 2 from Agresti and Coull (1998). Aliquots of suspensions from selected samples were also heat treated (95 °C for 20 min) to denature and deactiv

54 ENVIRONMENTAL SCIENCES↗