Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Effective temperature model”

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 415 records · Page 23

Modeling plasma-induced surface charge effects on CO 2 activation by single atom catalysts supported on reducible and irreducible metal oxides

Abstract The accumulation of negative surface charge on catalytic surfaces in the presence of low-temperature plasma (LTP) could influence catalytic performance. However, it is difficult to disentangle the role of surface charging and other LTP catalytic effects in experiment. Herein, we use density functional theory (DFT) modeling to understand the effect of plasma-induced surface charging on CO 2 activation by atomically dispersed single atom (SA) catalysts on both reducible and irreducible metal oxide supports. We model CO 2 adsorption strength and CO 2 dissociation barriers for Co 1 , Ni 1 , Cu 1 , Rh 1 , Pd 1 , and Ag 1 SAs on both reducible and irreducible supports, namely, CeO 2 (100), TiO 2 (101), and γ -Al 2 O 3 (110), to elucidate trends. We find that accumulated surface charge on the SA increases the CO 2 adsorption strength and decreases the CO 2 dissociation barrier for all studied SA/support combinations. For both charged and uncharged (neutral) systems, SAs on the reducible CeO 2 (100) support generally adsorb CO 2 more weakly compared to when on irreducible supports like γ -Al 2 O 3 (110). SAs on γ -Al 2 O 3 (110) typically have larger barriers for CO 2 dissociation for both charged and uncharged systems compared to TiO 2 (101) and CeO 2 (100). The magnitude of surface charging effects on CO 2 binding energies and dissociation barriers depends sensitively on both the SA and the support. In some cases, the CO 2 activation trends qualitatively change between neutral and charged systems for a fixed SA across different supports. This DFT modeling study demonstrates that surface charging should be considered in strong electric fields because it can have a large effect on molecule adsorption and bond-breaking on catalytic surfaces.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Transformational challenge reactor analysis to inform preconceptual core design decisions: Sensitivity study of transient analysis in a hydride-moderated microreactor

The Transformational Challenge Reactor (TCR) program aims to demonstrate a revolutionary design approach enabled by advanced manufacturing and data analytics in the design of nuclear reactors. This article discusses scoping analyses of preconceptual designs to inform TCR design decisions and the evaluation of sensitivities and uncertainties on postulated transient scenarios. The applicability of the systems codes TRACE and RELAP5-3D to TCR transient analysis are examined, and RELAP5-3D models are used to examine the transient response of two candidate core designs at multiple power levels. Then, the uncertainty quantification code RAVEN is used to quantify the effect of several design parameters on reactivity-initiated accident (RIA) progression at hot zero power (HZP) and hot full power (HFP) as well as to assess the impact of uncertainties in transient parameters for the pressurized loss of forced cooling (PLOFC).When results were compared, TRACE and RELAP5-3D showed good agreement in their ability to predict system behavior, but RELAP5-3D calculations were closer to analytical predictions for the RIA. Furthermore, models for a PLOFC accident in two designs (a UO2 and TRISO core) at multiple power levels showed greater temperature margins for the TRISO core at all power levels. Using this information, along with other scoping analyses and constraints, the TCR design team selected a power level of 3 MWth and a TRISO-based core design. For this design RAVEN was applied to vary RIA parameters in RELAP5-3D models at HZP and HFP to understand the effect on figures of merit such as peak power, fuel and coolant temperature, and energy deposition. This sensitivity study found that the inserted reactivity worth was the most important parameter controlling all figures of merit, but for insertion up to 1.5$ no failure of TCR fuel is anticipated. For constant reactivity insertions, the magnitude of the fuel temperature coefficient was found to have the greatest effect on all figures of merit under most circumstances. These results not only demonstrate the anticipated robust safety of the proposed TCR fuel form but also provide a reference for future metal-hydride moderated systems to understand RIA behavior. In the PLOFC, the impact of heat transfer enhancement due to wavy flow channel effects dominated the variance in peak temperatures, and variations in heat exchanger elevation provided the greatest control on natural circulation flow rate. Additionally, no fuel particle failure is anticipated in the PLOFC.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Modeling direct air carbon capture and storage in a 1.5 °C climate future using historical analogs

Limiting the rise in global temperature to 1.5 °C will rely, in part, on technologies to remove CO 2 from the atmosphere. However, many carbon dioxide removal (CDR) technologies are in the early stages of development, and there is limited data to inform predictions of their future adoption. Here, we present an approach to model adoption of early-stage technologies such as CDR and apply it to direct air carbon capture and storage (DACCS). Our approach combines empirical data on historical technology analogs and early adoption indicators to model a range of feasible growth pathways. We use these pathways as inputs to an integrated assessment model (the Global Change Analysis Model, GCAM) and evaluate their effects under an emissions policy to limit end-of-century temperature change to 1.5 °C. Adoption varies widely across analogs, which share different strategic similarities with DACCS. If DACCS growth mirrors high-growth analogs (e.g., solar photovoltaics), it can reach up to 4.9 GtCO 2 removal by midcentury, compared to as low as 0.2 GtCO 2 for low-growth analogs (e.g., natural gas pipelines). For these slower growing analogs, unabated fossil fuel generation in 2050 is reduced by 44% compared to high-growth analogs, with implications for energy investments and stranded assets. Residual emissions at the end of the century are also substantially lower (by up to 43% and 34% in transportation and industry) under lower DACCS scenarios. The large variation in growth rates observed for different analogs can also point to policy takeaways for enabling DACCS.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

A Survey of Multizone Airflow Simulation Software

H-Canyon simulates the building and process ventilation to ensure positive contamination control for any modifications as well as evaluating off-normal conditions. Simulations can also be used to identify air quality of personnel occupied areas. The models can be classified into two groups, single and multi-zone models. As the name states, a single zone model assumes that the building can be described by a single, well-mixed zone. Whereas large buildings have structures that must be described as multizone. Multizone airflow models deal with the difficulty of flows in a building by distinguishing the effects of internal flow restrictions. Several models have been developed that take internal partitions into account. CONTAM (currently used) is a multizone airflow program created by NIST to determine airflows including infiltration, exfiltration, wind pressures acting outside the building, buoyancy effects caused by the indoor and outdoor air temperature difference. COMIS is a multizone air infiltration, ventilation, and contaminant transport model that was developed at Lawrence Berkeley National Laboratory. COMIS deals with applied pressure distribution, airflow components, HVAC systems, schedules, and contaminant transport. EnergyPlus is designed by DOE. It supports flexible component level configuration of HVAC models and simulates sub-hourly time steps to control fast dynamic systems. The program features include providing accurate environmental performance data including energy consumption and carbon emissions, and reporting solar gains on surfaces. Heavent is used to model industrial ventilation systems to control airborne contaminants. It's designed based on the American Conference of Governmental Industrial Hygienists (ACGIH) Ventilation manual. DesignBuilder is the most recognized and advanced user interface to EnergyPlus. It gives access to all of the most frequently required simulation capabilities covering building fabric, thermal mass, glazing, shading renewables, and HVAC analysis. The objective of this study is to identify all relevant software available including commercially licensed, freely licensed, and open source, that can be used to simulate large, complex process facility ventilation systems. All discovered software tools are evaluated for salient features and compared to the software tool currently used to simulate the H-Canyon ventilation system. The design criteria of a ventilation system has two basic requirements: Supplying fresh air for the occupants and changing the air in the room sufficiently; therefore, the smells, fumes, and contaminants will be removed. In this study, first, whether the software simulates airflow (infiltration, natural ventilation, and interzonal) or considers them to be inputs is discussed. Second, if indoor contaminant levels are considered; and finally suggestions on how airflow analysis can be more appropriately accounted for when performing airflow simulation of a building. All of the airflow software reviewed for this study can account for constant infiltration and natural ventilation rates that are not affected by changes in indoor and outdoor conditions. Modelling infiltration in EnergyPlus or DesignBuilder are applicable for small buildings and not appropriate for larger buildings like the H-Canyon. The effect of wind on external pressures, infiltration, and natural and mechanical ventilation, can also be calculated using all of these software. These programs have control options to switch between natural to mechanical ventilation modes or both (hybrid). Interzonal airflow can be applied which means one zone can be specified to receive a fixed amount of air from another zone. The following ventilation modeling software is recommended for further testing: CONTAM, COMIS, Heavent. The next step in evaluating the software will consist of creating test facility simulations using evaluation copies of the software.

42 ENGINEERING↗

Volatilization Behavior of Supported Au Nanoparticle Arrays under H 2 at High Temperature

In this work, the fabrication and experimental study of model Au nanoparticle arrays (nanodots) was used to allow a direct comparison of reactive evaporation at elevated temperatures and nonambient H 2 -containing gas atmospheres with theoretical models. Strong temperature and H 2 -concentration dependencies were confirmed as expected, and the observed dependence in particle diameter with time was effectively described using a reactive evaporation and gas phase diffusion based quantitative model.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Effect of boundary conditions on fluctuations of the Bose condensate of interacting atoms

For a gas of trapped interacting atoms forming a Bose condensate, we demonstrate a possibility of a significant influence of the boundary conditions on the statistics of the number of particles in the fundamental mode of the system. The analysis is carried out in the Bogoliubov – Popov low-temperature approximation for model homogeneous cubic traps with periodic or combined periodic and zero boundary conditions. It is shown that the effect of the boundary conditions does not weaken even in the region of parameters corresponding to the relatively strong interaction in the Thomas – Fermi asymptotic case, and does not disappear when proceeding to the thermodynamic limit. (paper)

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

SATS Enhanced Capabilities and Demonstration of Improved Ramp Rates for In-Cell Testing

This report presents the development of the second-generation Severe Accident Test Station, referred to as SATS 2.0. SATS 2.0 represents a major advancement in addressing critical research needs for the United States nuclear industry, particularly in relation to fuel fragmentation, relocation, and dispersal (FFRD) concerns associated with the nuclear industry’s goal to extend fuel burnup beyond a peak rod average of 62 GWd/tU. The central objective of SATS 2.0 was the deployment a new furnace capable of achieving heating rates up to 100°C/s. This achievement addresses critical capability gaps, enabling assessment of prioritized research objectives established by the Electric Power Research Institute’s (EPRI) Collaborative Research on Advanced Fuel Technologies (CRAFT) working group. These research objectives are directly related to high-burnup loss-of-coolant accident (LOCA) conditions and the effects of prolonged exposure to elevated temperatures, both of which are crucial to enhancing nuclear safety and efficiency. The SATS 2.0 system builds upon prior experience with the original SATS system, which focused on evaluating accident-tolerant fuel (ATF) cladding concepts during accident conditions. However, the new system not only expands upon prior core capabilities by achieving higher heating rates with a better furnace but also plans to incorporate novel auxiliary systems to create a versatile platform for addressing current research needs. For example, a system was developed for the quantification and characterization of fission gas released during high-temperature transients. Additionally, in situ measurement capabilities were developed, such as digital image correlation which enabled the real-time capture of strain related to balloon and burst events and fiber-optic sensors that allowed for high-fidelity characterization of temperature gradients. Demonstration tests of the new 12-lamp furnace achieved heating rates up to 120°C/s. Notably, SATS 2.0 demonstrated its capacity to simulate LOCA burst tests at different pressures, yielding burst data that align with historical empirical models. Moreover, the system exhibited its capability to simulate complex conditions observed in anticipated operational occurrences, while effectively mitigating temperature overshoots. These accomplishments mark significant progress toward overcoming FFRD challenges and advancing the United States nuclear industry's safety basis and technical capabilities for extended burnup.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Uncertainty Quantification of Calculated Temperatures for the AGR 5/6/7 Experiment

This report documents the quantification of uncertainty of the calculated temperature data for the Advanced Gas Reactor (AGR) 5/6/7 fuel irradiation experiment conducted in the Advanced Test Reactor at Idaho National Laboratory in support of the Advanced Reactor Technologies? research and development program. Recognizing uncertainties inherent in physics and thermal simulations of the AGR 5/6/7 capsules, the results of the numerical simulations are used in combination with statistical analysis methods to improve qualification of measured data. The calculated fuel temperatures for AGR tests are also used for validation of the fission product transport and fuel performance simulation models. These crucial roles of the calculated fuel temperatures in ensuring achievement of the AGR experimental program objectives require accurate determination of the model temperature uncertainties. This report covers temperature uncertainty results for each of the five AGR 5/6/7 capsules. To quantify the uncertainty of calculated temperatures determined using the ABAQUS finite element heat transfer code, this study identifies and analyzes model parameters of potential importance to the calculated temperatures of fuel compacts and thermocouples. The selection of input parameters for uncertainty quantification is based on the ranking of their influences upon temperature predictions. Thus, selected input parameters include those with high sensitivity and those with the largest uncertainty. Propagation of model parameter uncertainty and sensitivity is then used to quantify the overall uncertainty of calculated temperatures. Measurement uncertainty, analysis of modeling assumptions, and expert judgment are used as the basis to quantify the uncertainty range for selected input parameters. The input uncertainties are dynamic, accounting for the effect of unplanned events and changes in thermal properties of capsule components over extended exposure to high temperatures and fast neutron irradiation. The sensitivity analysis performed in this work went beyond the traditional local sensitivity. Using experimental design, analysis of pairwise interactions of model parameters was performed to establish sufficiency of the time dependent first order (linear) expansion terms in constructing the temperature response surface. To achieve completeness, uncertainty propagation made use of pairwise noise correlations of model parameters. Furthermore, using an interpolation scheme over the input parameter domain, the analysis obtains time dependent sensitivity over the test campaign duration. This allows computation of uncertainty for the calculated fuel temperatures and the calculated graphite temperatures at thermocouple locations during the entire irradiation period.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Toward efficient elastocaloric systems: Predicting material thermal properties with high fidelity

A critical need to accurately model thermal behaviors of materials that exhibit strong elastocaloric effects, including predicting the effects themselves at varying stresses and temperatures, has been addressed using a simple and versatile physics-based approach. The key factor leading to the high precision is approximating the underlying elastic phase transition as a smooth modification of lattice entropy between coexisting phases. Once the phase transformation entropy is modeled to match experimentally measured strain as a function of temperature and applied stress, estimating the heat capacity, entropy, and isothermal entropy and adiabatic temperature changes in temperature-stress coordinates becomes straightforward. This approach provides insight into how thermal properties of elastocaloric materials vary through the transition based on strain measurements that are simple to perform and interpret. In addition to aiding in the rapid evaluation of new and existing elastocaloric materials, this advance is expected to prove invaluable for accurate heat transfer modeling aimed at designing efficient regenerative elastocaloric cooling devices.

36 MATERIALS SCIENCE↗

A mortar thermomechanical contact computational framework for nuclear fuel performance simulation

Nuclear fuel performance simulations involve the modeling of complex physical phenomena, ranging from fission gas release to fuel swelling and other temperature-induced effects. For light-water reactors (LWRs), swelling of the fuel and the pressure it imposes on the clad when they come into contact causes permanent clad deformation. Accurately characterizing the fuel-cladding interaction, which involves multiple physics, is essential to accurately simulate the fuel/cladding system. Thermomechanical modeling of this problem using a variationally consistent enforcement (e.g., a mortar approach) has been shown to improve the quality of results and facilitate convergence. Here, we present a general multiphysics computational framework for solving nuclear fuel problems using a mortar approach in BISON, a nuclear fuel performance code. In this study analyses show that using the mortar approach, which enables variationally consistent constraint enforcement, improves the quality of results as compared to the more commonly used node-on-face enforcement for representative LWR nuclear fuel simulations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Estimation of secondary organic aerosol viscosity from explicit modeling of gas-phase oxidation of isoprene and α -pinene

Abstract. Secondary organic aerosols (SOA) are major components of atmospheric fine particulate matter, affecting climate and air quality. Mounting evidence exists that SOA can adopt glassy and viscous semisolid states, impacting formation and partitioning of SOA. In this study, we apply the GECKO-A (Generator of Explicit Chemistry and Kinetics of Organics in the Atmosphere) model to conduct explicit chemical modeling of isoprene photooxidation and α-pinene ozonolysis and their subsequent SOA formation. The detailed gas-phase chemical schemes from GECKO-A are implemented into a box model and coupled to our recently developed glass transition temperature parameterizations, allowing us to predict SOA viscosity. The effects of chemical composition, relative humidity, mass loadings and mass accommodation on particle viscosity are investigated in comparison with measurements of SOA viscosity. The simulated viscosity of isoprene SOA agrees well with viscosity measurements as a function of relative humidity, while the model underestimates viscosity of α-pinene SOA by a few orders of magnitude. This difference may be due to missing processes in the model, including autoxidation and particle-phase reactions, leading to the formation of high-molar-mass compounds that would increase particle viscosity. Additional simulations imply that kinetic limitations of bulk diffusion and reduction in mass accommodation coefficient may play a role in enhancing particle viscosity by suppressing condensation of semi-volatile compounds. The developed model is a useful tool for analysis and investigation of the interplay among gas-phase reactions, particle chemical composition and SOA phase state.

54 ENVIRONMENTAL SCIENCES↗

Accumulation-Type Ohmic van der Waals Contacts to Nearly Intrinsic WSe 2 Nanosheet-Based Channels: Implications for Field-Effect Transistors

We report the fabrication of ohmic van der Waals (vdW) contacts to nearly intrinsic WSe 2 nanosheet-based channels in field-effect transistors (FETs) using degenerately p-doped MoS2 (p + -MoS2) as a contact metal. We demonstrate that accumulation-type ohmic contacts and the high device performance are achievable without electrostatically gating the drain/source contact regions despite the nearly intrinsic nature of WSe 2 . Back-gated WSe 2 FETs with p + -MoS 2 bottom contacts (which screen the back-gate electric field in the drain/source regions) exhibit linear output characteristics, a high on/off ratio of 10 8 , and a high two-terminal field-effect mobility up to ~200 cm2 V –1 s –1 at room temperature. Our theoretical modeling reveals that the p + -MoS 2 /WSe 2 vdW junction behaves like a metal/semiconductor ohmic contact signified by a vanishingly thin space-charge region of ~1 nm on the p + -MoS 2 side and a substantial accumulation layer of free holes on the WSe 2 side, which is further verified by additional temperature-dependent and dual-gated measurements of WSe 2 FETs. We attribute the formation of accumulation-type ohmic contacts free of a Schottky barrier to the near absence of Fermi-level pinning at the vdW interface and the work function of p + -MoS 2 being larger than the ionization energy of WSe 2 . This study represents an important step toward achieving low-resistance ohmic contacts to two-dimensional (2D) semiconductors by eliminating the Fermi-level pinning effects, which is expected to have significant implications for next-generation 2D semiconductor-based nanoelectronics.

36 MATERIALS SCIENCE↗

THE IMPACT OF OXIDATION-INDUCED DEGRADATION ON MATERIALS USED IN HYDROGEN-FIRED MICROTURBINES

Hydrogen-fueled microturbines are being considered as part of the future green microgrid. However, the use of hydrogen as a fuel presents new challenges for selection and development of suitable high temperature materials for hydrogen combustion. The burning of hydrogen is expected to result in higher operating temperatures and higher than typically observed water vapor contents in exhaust gases than after burning natural gas. In the present work, foil specimens of various Fe- and Ni-based alloys were oxidized in air + 10vol.% H2O and air + 60vol.% H2O for up to 5,000 h at 700 °C to simulate the exhaust atmosphere of natural gas and hydrogen-fueled microturbines. The impact of alloy composition and water vapor content on the oxidation/ volatilization induced loss of wall thickness was experimentally evaluated. Enhanced external oxidation and volatilization of Cr¬2O3 and Ti-doped Cr2O3 scales was observed in 60 % H2O compared to 10% H2O. No significant impact of the higher water vapor content was observed on Al2O3 scales formed on Fe-based alumina forming alloys. Lifetime modeling was employed to predict the combined effects of water vapor content, gas flow rates, temperature and alloy composition on the oxidation-induced lifetime of the investigated materials

Romedenne, Marie↗

Towards understanding the magnetic properties of the breathing pyrochlore compound Ba3Yb2Zn5O11through single-crystal studies

Abstract Ba 3 Yb 2 Zn 5 O 11 is exceptional among breathing pyrochlore compounds for being in the nearly-decoupled limit where inter-tetrahedron interactions are weak, hosting isolated clusters or molecular magnet-like tetrahedra of magnetic ytterbium (Yb 3+ ) ions. In this work, we present the study carried out on single-crystal samples of the breathing pyrochlore Ba 3 Yb 2 Zn 5 O 11 , using a variety of magnetometry and neutron scattering techniques along with theoretical modeling. We employ inelastic neutron scattering to investigate the magnetic dynamics as a function of applied field (with respect to both magnitude and direction) down to a temperature of 70 mK, where inelastic scattering reveals dispersionless bands of excitations as found in earlier powder sample studies, in good agreement with a single-tetrahedron model. However, diffuse neutron scattering at zero field and dc-susceptibility at finite field exhibit features suggesting the presence of excitations at low-energy that are not captured by the single tetrahedron model. Analysis of the local structure down to 2 K via pair distribution function analysis finds no evidence of structural disorder. We conclude that effects beyond the single tetrahedron model are important in describing the low-energy, low-temperature physics of Ba 3 Yb 2 Zn 5 O 11 , but their nature remains undetermined.

36 MATERIALS SCIENCE↗

Simulation of a TRISO MiniFuel irradiation experiment with data-informed uncertainty quantification

An irradiation experiment using tristructural isotropic (TRISO) fuel particles and the miniature fuel (MiniFuel) irradiation vehicle was performed in Oak Ridge National Laboratory’s High Flux Isotope Reactor (HFIR) to support development of the Kairos Power fluoride salt–cooled, high-temperature reactor (KP-FHR). Here, this paper describes modeling predictions of temperatures and fuel burnup for the as-built experiment. An uncertainty quantification (UQ) analysis was performed to determine the effect of TRISO particle volume and position on the temperature predictions at various fuel heat generation rates (HGRs). This UQ study utilized fuel kernel position and volume measurements previously collected using X-ray computed tomography (XCT) techniques and Monte Carlo sampling methods to generate fuel compact cases that were then analyzed using a finite element thermal model. The UQ analysis indicated that uncertainty in calculated temperatures caused by varying TRISO particle arrangement is relatively small, even at high fuel HGR. Final predictions of particle temperatures throughout the irradiation are shown to be relevant to KP-FHR normal and off-normal operating conditions and to previous TRISO irradiation experiments. The combination of XCT with UQ analyses will inform post-irradiation examination (PIE) of the irradiated fuel compacts, and these analyses can be used to develop fuel performance models for coated particle fuel forms. Both PIE of separate-effects irradiation data and enhanced fuel performance modeling support accelerated qualification of TRISO fuels for a broad range of advanced reactor applications. The novel approach demonstrated here of measuring TRISO particle configurations with XCT methods and generating representative fuel compacts for finite element modeling and UQ analysis could be leveraged by the broader particle fuel community in the development of other TRISO fuel experiments in which these variables may have a significant impact on key outcomes.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Integrating plant physiology into simulation of fire behavior and effects

Summary Wildfires are a global crisis, but current fire models fail to capture vegetation response to changing climate. With drought and elevated temperature increasing the importance of vegetation dynamics to fire behavior, and the advent of next generation models capable of capturing increasingly complex physical processes, we provide a renewed focus on representation of woody vegetation in fire models. Currently, the most advanced representations of fire behavior and biophysical fire effects are found in distinct classes of fine‐scale models and do not capture variation in live fuel (i.e. living plant) properties. We demonstrate that plant water and carbon dynamics, which influence combustion and heat transfer into the plant and often dictate plant survival, provide the mechanistic linkage between fire behavior and effects. Our conceptual framework linking remotely sensed estimates of plant water and carbon to fine‐scale models of fire behavior and effects could be a critical first step toward improving the fidelity of the coarse scale models that are now relied upon for global fire forecasting. This process‐based approach will be essential to capturing the influence of physiological responses to drought and warming on live fuel conditions, strengthening the science needed to guide fire managers in an uncertain future.

54 ENVIRONMENTAL SCIENCES↗

Ten-moment fluid model for low-temperature magnetized plasmas

In this paper, a one-dimensional 10-moment multi-fluid plasma model is developed and applied to low-temperature magnetized plasmas. The 10-moment model solves for six anisotropic pressure terms, in addition to density and three components of fluid momentum, which allows the model to capture finite kinetic effects. The results are benchmarked against a 5-moment model, which assumes that the gas constituents follow a Maxwellian velocity distribution function (VDF), and a particle-in-cell Monte Carlo collision model that allows for arbitrary non-Maxwellian VDFs. The models are compared in a one-dimensional, low-temperature, partially magnetized plasma test case. The 10-moment results accurately reproduce the anisotropic temperature profile in low-temperature magnetized plasmas, where shear gradients exist due to the E×B drift. We discuss the mechanisms by which the anisotropic pressure can be generated in low-temperature magnetized plasmas. In addition, the importance of a self-consistent heat flux closure to the 10-moment model is studied, showing consistency with other models only when the assumptions of the underlying model are met. The 10-moment model allows for study of electron inertia effects and non-Maxwellian VDFs without the need for kinetic methods that are more computationally expensive.

Kuldinow, Derek Amur (ORCID:0000000319730196)↗

A Monte Carlo Approach to Approximating the Effects of Pore Geometry on the Phase Behavior of Soil Freezing

Abstract Freezing in porous media is associated with a host of dynamic phenomena that stem from the presence and mobility of premelted liquid at subzero temperatures. Accurate assessments of the progressive liquid‐ice phase transition is required for predictive models of frost damage, glacier‐till coupling, and many other cold regions processes, as well as for evaluating the capacity for water storage in near‐surface extraterrestrial environments. We use a Monte Carlo approach to sample the pore space in a synthetic 3D packing of poly‐dispersed spherical particles and evaluate local geometrical constraints that allow us to assess changes in the relative proportions of pore fluid and ice. By approximating the phase boundary geometry in fine‐grained pores while considering both the curvature of the liquid‐ice interface and wetting interactions with matrix particles, our model predicts changes in phase equilibrium in granular media over a broad temperature range, where present accounting for the colligative effects of chloride and perchlorate solutes. In addition to formulating the constitutive behavior needed to better understand properties and processes in frozen soils, our results also provide insight into other aspects of phase equilibria in porous media, including the formation of methane hydrates in permafrost and marine sediments, and the partitioning between liquid water and vapor in the vadose zone.

54 ENVIRONMENTAL SCIENCES↗