Search NASA⌕ Search

SEARCH · Search NASA

Results for “Salt convection”

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.

90 records · Page 5

Generating a stratocumulus-like cloud top in a convection-cloud chamber

Stratocumulus-topped boundary layers play a crucial role in influencing daily weather and earth energy balance. Entrainment at the stratocumulus cloud top affects the cloud’s lifetime, precipitation, and radiative properties, but our understanding remains limited due to the lack of resolution in both field observations and numerical simulations. A recently proposed convection-cloud chamber with detailed control of sidewall temperatures can provide a unique opportunity to explore this mechanism in a laboratory setting. In this work, we use numerical simulations to demonstrate that this design can produce a cloud top that mimics the entrainment interfacial layer in a stratocumulus cloud. Our results show that a steady-state cloud can be formed by cooling the lower portions of the sidewalls and warming the bottom surface, while a temperature inversion at the cloud top can be generated by keeping the upper sidewalls and top surface warmer than the bottom. The turbulent kinetic energy profile and budget are similar to those found in a convective boundary layer, and inhomogeneous mixing near the cloud top can be observed. These findings significantly enhance the scientific value of constructing the tall convection-cloud chamber.

54 ENVIRONMENTAL SCIENCES↗

The 10 September 2025 M w 4.1 Earthquake in Northeastern Utah, United States: An Archetypal Continental Mantle Event

The 10 September 2025 M w 4.1 earthquake in northeastern Utah, United States, had a focal depth 68 km beneath sea level, which is ∼20–25 km greater than estimates of local crustal thickness, making it a rare example of a continental mantle earthquake (CME). The focal depth is well resolved from arrival-time inversion (nearest station ∼13 km away) and moment tensor inversion of regional waveforms. Similar to other CMEs in the Intermountain West, there were no obvious aftershocks or foreshocks, and the waveforms were enriched in high-frequency energy. Spectral modeling gives a stress drop of ∼80 MPa and a radiation efficiency of ∼0.08, albeit with large uncertainties. The high stress drop and low radiation efficiency are consistent with a dissipative source process such as thermal runaway. Also similar to previous Intermountain West CMEs, the event occurred along the boundary of the Archean Wyoming craton, where pressure–temperature conditions favor ductile deformation. We hypothesize that edge-driven or regional-scale mantle convection produces increased strain rates near the craton boundary that make either conventional brittle failure or thermal runaway feasible at relatively high pressure–temperature conditions. High conductivity inferred around the edge of the craton may suggest that fluids also contribute to CME occurrence.

Koper, Keith D. [Univ. of Utah, Salt Lake City, UT↗

NRC Multiphysics Analysis Capability Deployment (FY2021--Part 1)

This report details progress and activities of Idaho National Laboratory (INL) on the Nuclear Regulatory Commission (NRC) project “Development and Modeling Support for Advanced Non-Light Water Reactors.” The tasks completed for this report are as follows: First, Task 1d: The net radiation transfer method was implemented into MOOSE for modeling reactor cavity cooling system geometries. RCCS models for two experiments were created: (1) Natural Convection Shutdown Heat Removal Test Facility (NSTF) experiment R022, and (2) HTTR VCS mockup. For validation, computed temperature distributions were compared to measured temperatures. Next, Task 4c: An algorithm for computing the pebble bed reactor equilibrium core isotopic com-position was developed and an initial version is implemented into the reactor multi-physics code Griffin. Initial results for a simplified axisymmetric pebble bed reactor are presented. Finally, Task 7: generation of a reference plant model for molten salt cooled pebble bed reactor.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Aerosol responses to precipitation along North American air trajectories arriving at Bermuda

North American pollution outflow is ubiquitous over the western North Atlantic Ocean, especially in winter, making this location a suitable natural laboratory for investigating the impact of precipitation on aerosol particles along air mass trajectories. We take advantage of observational data collected at Bermuda to seasonally assess the sensitivity of aerosol mass concentrations and volume size distributions to accumulated precipitation along trajectories (APT). The mass concentration of particulate matter with aerodynamic diameter less than 2.5 µm normalized by the enhancement of carbon monoxide above background (PM 2.5 /ΔCO) at Bermuda was used to estimate the degree of aerosol loss during transport to Bermuda. Results for December–February (DJF) show that most trajectories come from North America and have the highest APTs, resulting in a significant reduction (by 53 %) in PM 2.5 /ΔCO under high-APT conditions (> 13.5 mm) relative to low-APT conditions (< 0.9 mm). Moreover, PM 2.5 /ΔCO was most sensitive to increases in APT up to 5 mm (–0.044 µg m –3 ppbv –1 mm –1 ) and less sensitive to increases in APT over 5 mm. While anthropogenic PM 2.5 constituents (e.g., black carbon, sulfate, organic carbon) decrease with high APT, sea salt, in contrast, was comparable between high- and low-APT conditions owing to enhanced local wind and sea salt emissions in high-APT conditions. The greater sensitivity of the fine-mode volume concentrations (versus coarse mode) to wet scavenging is evident from AErosol RObotic NETwork (AERONET) volume size distribution data. A combination of GEOS-Chem model simulations of the 210 Pb submicron aerosol tracer and its gaseous precursor 222 Rn reveals that (i) surface aerosol particles at Bermuda are most impacted by wet scavenging in winter and spring (due to large-scale precipitation) with a maximum in March, whereas convective scavenging plays a substantial role in summer; and (ii) North American 222 Rn tracer emissions contribute most to surface 210 Pb concentrations at Bermuda in winter (~75 %–80 %), indicating that air masses arriving at Bermuda experience large-scale precipitation scavenging while traveling from North America. A case study flight from the ACTIVATE field campaign on 22 February 2020 reveals a significant reduction in aerosol number and volume concentrations during air mass transport off the US East Coast associated with increased cloud fraction and precipitation. These results highlight the sensitivity of remote marine boundary layer aerosol characteristics to precipitation along trajectories, especially when the air mass source is continental outflow from polluted regions like the US East Coast.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Simulating Droplet-Resolved Haze and Cloud Chemistry Forming Secondary Organic Aerosols in Turbulent Conditions within Laboratory and Cloud Parcels

Most of our existing knowledge of cloud chemistry in regards to forming secondary organic aerosols (SOA) is based on measurements in bulk aqueous solutions. However, SOA reaction kinetics derived from bulk solution measurements might differ from the kinetics in actual cloud droplets, since turbulent mixing and ionic strengths, and ratio of surface area to volume in individual cloud droplets might vary substantially from bulk solutions in the real atmosphere. Three-dimensional models at various scales have been used to simulate aqueous chemistry. However, most of these models do not resolve turbulence down to the smallest length scales of 1 mm and do not simulate cloud chemistry in individual cloud droplets due to large computational costs. Here we incorporate the formation of isoprene epoxydiol SOA (IEPOX-SOA) in individual droplets within a one-dimensional explicit mixing parcel model (EMPM-Chem). We apply EMPM-Chem to simulate turbulence and droplet-resolved IEPOX-SOA formation using a configuration based on the Michigan Tech Pi chamber. We find that the dissolution of IEPOX gases is weighted more towards larger cloud droplets due to their large liquid water content (compared to smaller droplets), while the conversion of dissolved IEPOX to IEPOX-SOA is much greater within smaller deliquesced haze particles due to their higher acidity and ionic strengths compared to cloud droplets. We also find that as droplet residence times increase in the chamber, e.g., due to increasing ammonium bisulfate seed aerosol injection rates and/or increasing heights of the chamber, formation of IEPOX-SOA increases substantially. Thus, our EMPM-Chem model could be used to design future cloud chambers to maximize SOA production from cloud chemistry. We also apply the EMPM-Chem model to simulate how IEPOX-SOA formation evolves in individual cloud droplets within rising cloudy parcels in the atmosphere. We find that as subsaturated air is entrained into and turbulently mixed with the cloud parcel, evaporation causes a reduction in droplet sizes, which leads to corresponding increases in per droplet ionic strength and acidity. Increased droplet acidity in turn greatly accelerates the kinetics of IEPOX-SOA formation. In conclusion, our results provide key insights into single-cloud-droplet chemistry, suggesting that entrainment mixing may be an important process that increases SOA formation in the real atmosphere.

54 ENVIRONMENTAL SCIENCES↗

Updraft Width Modulates Ambient Atmospheric Controls on Convective Cloud Depth

Abstract The depth of convective clouds affects vertical transport of atmospheric constituents, influencing downstream weather and climate. Atmospheric controls on the maximum depth reached by moist convection are investigated with radar‐tracked convective cells tagged with sounding‐derived atmospheric parameters from a field campaign in central Argentina. Regression analyses show that narrow (<12‐km diameter) and wide (>16‐km diameter) cell depths respond to disparate factors, where cell areas are defined using composite reflectivity signatures. Undiluted lifted parcel indices including convective available potential energy (CAPE) and level of neutral buoyancy (LNB) are top predictors of wide cell maximum depth while mid‐tropospheric relative humidity is the top predictor of narrow cell maximum depth. Because narrow cells are more numerous than wide cells, the overall outcome of the full cell population does not strongly correlate with CAPE and LNB conditions. Tracked cells and atmospheric conditions in a simulation with 3‐km grid spacing covering the field campaign produce similar results to those observed. Narrow cells that are relatively deep have a cooler and moister mid‐troposphere with weaker free tropospheric subsidence, while relatively deep wide cells have much warmer and moister lower tropospheric conditions. These atmospheric differences are present 1 hr before cell initiation at both a fixed observing site and variable cell initiation locations. Simulated narrow cell maximum equivalent potential temperature decreases with height at a rate similar to the ambient vertical gradient, causing these cells to fall short of their LNB and supporting the view that entrainment‐driven dilution is a dominant control on their depth.

54 ENVIRONMENTAL SCIENCES↗

3D Cloud Aerosol Precipitation Experiment at kennaook Cape Grim (3D-CAPE-k) Field Campaign Report

The objective of this campaign was to complement the U.S. Department of Energy Atmospheric Radiation Measurement (ARM) User Facility’s Cloud And Precipitation Experiment at kennaook (CAPE-k) in northwest Tasmania by deploying scanning, fixed pointing angle, and in situ aerosol, cloud, and precipitation remote-sensing instruments during the last six months of the experiment. These instruments provided a three-dimensional (3D) context to the second ARM Mobile Facility (AMF2) vertically pointing observations collected at the kennaook Cape Grim (KCG) site and captured a portion of the life cycle of the clouds before and after the passage over the AMF2 vertical column. These 3D cloud measurements will be used in the near future to 1) investigate how the vertical profiles of aerosol and cloud properties are representative of the broader area and how these properties evolve during the portion of the cloud life cycle captured by the scanning instruments, 2) quantify how accurate the cloud fraction derived with assumptions from vertically pointing observations are compared with direct 3D cloud fraction measurements, 3) evaluate cloud fraction and liquid water path in the Australian Community Climate and Earth-System Simulator–Convective-scale (ACCESS-C) forecast model, and 4) evaluate aerosol and cloud products from the European Space Agency/Japanese Aerospace Exploration Agency (ESA/JAXA) Earth Cloud Aerosol Radiation Explorer (EarthCARE) mission using the scanning measurements, offering more chances of exact collocation, thereby complementing the statistical approach we are planning to employ with AMF2 observations for that same purpose.

54 ENVIRONMENTAL SCIENCES↗

Interactions Between Aerosols, Meteorology, and Early Convective Cloud Lifecycle as Measured During CACTI (Final Technical Report)

The research supported by this award sought to improve the understanding and forecasting of thunderstorms. We did so by using data from the RELAMPAGO-CACTI project, which deployed a suite of instruments around a mountain range in Argentina that sees thunderstorms erupt over the same general area almost daily. The main portion of our research examined the relationship between the concentration of atmospheric particulates (dust, smoke, etc) on the intensity of thunderstorms. Despite prior studies finding that increased particulate concentration corresponded to more intense storms, we found no effect, or if anything a slightly opposite effect.

54 ENVIRONMENTAL SCIENCES↗

Data-Driven Reduced-Order Modeling of Convective Heat Transfer in Porous Media

This work presents a data-driven Reduced-Order Model (ROM) for parametric convective heat transfer problems in porous media. The intrusive Proper Orthogonal Decomposition aided Reduced-Basis (POD-RB) technique is employed to reduce the porous medium formulation of the incompressible Reynolds-Averaged Navier–Stokes (RANS) equations coupled with heat transfer. Instead of resolving the exact flow configuration with high fidelity, the porous medium formulation solves a homogenized flow in which the fluid-structure interactions are captured via volumetric flow resistances with nonlinear, semi-empirical friction correlations. A supremizer approach is implemented for the stabilization of the reduced fluid dynamics equations. The reduced nonlinear flow resistances are treated using the Discrete Empirical Interpolation Method (DEIM), while the turbulent eddy viscosity and diffusivity are approximated by adopting a Radial Basis Function (RBF) interpolation-based approach. The proposed method is tested using a 2D numerical model of the Molten Salt Fast Reactor (MSFR), which involves the simulation of both clean and porous medium regions in the same domain. For the steady-state example, five model parameters are considered to be uncertain: the magnitude of the pumping force, the external coolant temperature, the heat transfer coefficient, the thermal expansion coefficient, and the Prandtl number. For transient scenarios, on the other hand, the coastdown-time of the pump is the only uncertain parameter. The results indicate that the POD-RB-ROMs are suitable for the reduction of similar problems. The relative L2 errors are below 3.34% for every field of interest for all cases analyzed, while the speedup factors vary between 54 (transient) and 40,000 (steady-state).

German, Péter (ORCID:0000000307285283)↗

Uncertainty estimation of bifurcated solutions in the Rayleigh–Bénard problem for advanced nuclear reactors applications

Multiphysics models of nuclear reactors frequently comprise nonlinear systems of equations. The nonlinear nature of these models could lead to solution bifurcations, where a small change in a certain parameter, e.g., the thermophysical properties of the coolant, can lead to a sudden change in the system’s behavior. At the point in parameter space where this happens, called a critical point, the Jacobian matrix of the model’s nonlinear operator becomes singular potentially permitting multiple solutions to coexist. In this paper, we perform uncertainty estimation (UE) in a parameter range that includes bifurcated solutions within the context of Rayleigh–Bénard problem. We perform this analysis assuming uncertain temperature difference, and tilt angle for the iterative solution algorithm with a unit Prandtl number (Pr = 1). Also, we perform this analysis under uncertain thermophysical properties for both FLiBe molten salt and liquid sodium as working fluid. We deploy two approaches to compute statistical moments for the resulting distributions of selected flow-field variables. The first approach is the blind computation of the mean and the standard deviation without any consideration of solution bifurcation, while the second approach utilizes k-means clustering to cluster each branch’s solutions together and compute separate statistical moments for each branch. The statistical distributions are obtained by perturbing the selected parameters about nominal values that correspond to a solution on one of the valid branches, and that solution is used as initial guess for the iterative solution algorithm. We found that perturbation of any parameter when its nominal value is close to its critical point always leads to branch jumping, i.e., the iterations converge to a solution on a branch different from the branch of the initial guess. This produces a statistical ensemble comprised of fundamentally different solutions leading to wrong mean values and uncertainty estimates, whereas clustering provides an efficient way to deal with this type of computation. This work is important for developing Gen IV nuclear systems because many of these systems rely on natural convection for cooling especially in accident conditions.

97 - MATHEMATICS AND COMPUTING↗

Measurements of the thermal conductivity of reference liquids using a modified transient hot-wire needle probe

Many power generation processes feature some type of working fluid that plays a critical role in the transfer of heat throughout the system. To optimize thermal efficiencies, an accurate characterization of the temperature-dependent thermal conductivity of the fluid is needed. Herein, this study develops a time-domain Modified Transient Hot-wire Needle Probe to overcome phenomena that have impaired traditional measurement techniques when the fluid is both corrosive and at an elevated temperature. First, the physical design of the instrument is presented to demonstrate features that limit interference from corrosion, electrical conduction, and convective and radiative heat transfer. Next, a governing equation for 1D (radial) heat conduction is fully developed using the thermal quadrupoles technique. A finite-element study is performed for validation. A sensitivity analysis is done to assess the time-dependent influence of critical parameters in the governing equation for several samples. Finally, thermal conductivity measurements are presented for seven different liquid substances by fitting the thermal model to experimental data across a range of temperatures. Reference correlations are included for comparison. The sample thermal conductivity and ambient temperature of the experimental conditions range from 0.125–16.5$\frac{W}{m·K}$ and 10–400°C, respectively. The total uncertainty, comprised of both the experimental design and fitting process, of all measurements is quantified.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

High-Fidelity and High-Performance Computational Simulations for Rapid Design Optimization of Sulfur Thermal Energy Storage

Industrial process heating (IPH) accounts for approximately 70% of US manufacturing energy use and is primarily produced by fossil fuel combustion. Approximately 1500 TWht (approximately 60%) of IPH demand is in the temperature range of 100-300. Industrial applications in this temperature range include drying, hydrothermal processing, thermal enhanced oil recovery, food and beverage, bioethanol production, etc. Cost-effective thermal energy storage (TES) that increases the utilization of waste and renewable heat (solar, geothermal, etc.) could provide significant energy savings and reliable heat sources, decrease emissions, and increase US manufacturing competitiveness through reductions in fuel consumption. TES development has historically been dominated by technologies suitable for deployment with concentrating solar power (CSP). State-of-the-art thermal storage deployed commercially with power tower CSP plants uses a 60%/40% NaNO3/KNO3 molten salt and operates between temperatures of approximately 280 degrees Celsius and 570 degrees Celsius using a two-tank configuration. However, these nitrate salts are unsuitable for operation outside of this temperature range due to a high freezing point of approximately 220 degrees Celsius, and limits on high-temperature salt stability and corrosion resistance of containment alloys. Other materials being investigated for TES include those based on: (1) sensible energy storage (various molten salt compositions, inert solid particles, rocks or pebble beds, sulfur, water, concrete, graphite, etc.), (2) latent energy storage in materials that undergo solid-liquid phase change at relevant temperatures (organic materials for low-temperature applications, inorganic salts and/or metals for high-temperature applications), or (3) thermochemical energy storage (hydrides, hydroxides, carbonates, metal oxides, etc.). The application temperature and challenges pertaining to storage material and/or containment cost, energy density, long-term thermal and cyclic stability, and charge/discharge heat transfer effectiveness drive material selection for a given IPH or electricity generation application. Sulfur is a cheap commodity at $80/ton compared to $1100 - 1300/ton for conventional salts. When using a metric of storage cost per kWh, sulfur costs around 2-3 $/kWh. Previous sulfur TES development focused on high temperature (>600 degrees) concentrated solar power applications with sulfur encapsulated in pipes and flow of gaseous HTF (air) in the shell side. However, for lower-temperature IPH applications in the range of approximately 100-300 degrees Celsius Element 16 adopted a compact and scalable TES design with molten sulfur in the shell and HTF pipes submerged in the molten sulfur bath. The low-cost molten sulfur TES for dispatchable IPH has deployment potential for broad applications. The spatial and temporal evolution of the HTF and sulfur temperature is critical to the TES system performance, and thus detailed modeling can improve understanding of the performance and facilitate design improvements. Using high performance computing and computational fluid dynamics (CFD) a low-cost molten sulfur thermal energy storage (TES) system for industrial process heating (IPH) applications was developed. The unique challenges in CFD modeling of sulfur TES are the sharp property changes of sulfur relevant to the working temperatures. Above 159, liquid sulfur undergoes polymerization, and the viscosity of sulfur rapidly increases by several orders of magnitude between 159 degrees Celsius and 188 degrees Celsius, followed by a decrease in viscosity beyond 188 degrees Celsius due to thermal bound dissociation. In addition, various concentrations of H2S impurities can also modify sulfur viscosity. This numerical challenge is especially relevant to transient simulation of the sulfur TES charging and discharging processes as the extreme property variations limit the applicability of traditional heat transfer correlations. Transient CFD simulations including the temperature-dependent sulfur properties and geometric complexity of the TES design were used to predict the effect of natural convection during charging and discharging on the heat transfer process, sulfur temperature uniformity, charge/discharge rates, and performance of the storage devices. The CFD model was validated with experimental results for a full charge and discharge cycle. The work will show 3D and 2D simulation comparisons aimed to facilitate rapid design iterations and a machine learning based design optimization approach.

CFD↗

Prototype Modeling for a Light-Trapping Planar-Cavity Enclosed Particle Solar Receiver

Concentrating solar thermal (CST) systems present a promising avenue for affordable and reliable energy production. Solar receivers are key components that determine the efficiency and longevity of these systems. Particle-based solar receivers have emerged as a compelling alternative to traditional technologies, offering several advantages that address limitations in current CST systems. This is especially true as next-generation CST technologies target applications including electricity generation, thermochemical processes, and industrial process heat, many of which necessitate higher operating temperatures than current commercial molten salt systems. Molten-salt thermal energy storage (TES) systems, commonly used in CSP, face challenges related to freezing and corrosion. Particle-based TES systems, in contrast, do not experience these issues, as particles are stable at high temperatures, exceeding 1000 degrees Celsius. This capability allows for a wider range of applications, including those requiring higher temperatures for industrial processes and efficient electricity generation. A novel innovation in particle-based solar receiver technology is the light-trapping planar cavity receiver (LTPCR) configuration developed by NREL. The LTPCR design consists of small cavity-like structures using opaque planar surfaces, enabling efficient capture and absorption of solar energy. A high incident flux concentration at the cavity aperture is absorbed on the receiver walls, and subsequently transferred to particles on the inside of cavities. The particles flow through the system, forming a fluidized bed inside of the receiver panels, effectively capturing the absorbed solar heat. Air is used as a fluidizing medium in this process to enhance particle heat transfer and mixing. The effectiveness of this design lies in its ability to manage solar flux conditions and ensure high solar-to-thermal receiver efficiency. A 100-kW prototype is currently being tested at the King Saud University in Saudi Arabia to assess the receiver performance. A range of modeling analyses for the optical, thermal, and mechanical effects were conducted to assess the performance of the receiver under on-sun conditions. The solar flux resulting from the KSU heliostat field was modeled using NREL SolTrace software and produced up to 600 kW/m2 at the receiver aperture. The solar flux absorbed on the receiver walls was then used within a computational fluid dynamics (CFD) model to predict wall temperature distributions along with radiation and convection loss. A two-phase CFD model was developed for the fluidized bed of silica sand inside the receiver panels to predict local wall-to-particle heat transfer coefficients, particle temperature distributions, and outlet temperature of the particles. We have also conducted analyses to understand the thermomechanical behavior of these innovative enclosed light-trapping solar receivers optimized for particle heating. We used finite element analysis (FEA) to predict the receiver's performance using temperature distributions obtained from CFD and based on the resulting stress profiles, evaluated creep-fatigue damage with a goal of achieving a 30-year service life. Analysis showed a significant impact of the particle-to-wall heat transfer coefficients (HTCs) on receiver performance, with higher HTCs resulting in reduced stress and increased lifespan. For instance, when using Inconel 740H, increasing the HTC from 800 W/m2 K to 1400 W/m2 K increased the creep life from 4,000 hours to over 100,000 hours. This highlights the importance of understanding and optimizing heat transfer in the design of high-efficiency receivers.

14 SOLAR ENERGY↗

A Model Intercomparison Study of Aerosol‐Cloud‐Turbulence Interactions in a Cloud Chamber: 1. Model Results

This study presents the first model intercomparison of aerosol‐cloud‐turbulence interactions in a controlled cloudy Rayleigh‐Bénard Convection chamber environment, utilizing the Pi Chamber at Michigan Technological University. We analyzed simulated cloud chamber‐averaged statistics of microphysics and thermodynamics in a warm‐phase, cloudy environment under steady‐state conditions at varying aerosol injection rates. Simulation results from seven distinct models (DNS, LES, and a 1D turbulence model) were compared. Our findings demonstrate that while all models qualitatively capture observed trends in droplet number concentration, mean radius, and droplet size distributions at both high and low aerosol injection rates, significant quantitative differences were observed. Notably, droplet number concentrations varied by over two orders of magnitude between models for the same injection rates, indicating sensitivities to the model treatments in droplet activation and removal and wall fluxes. Furthermore, inconsistencies in vertical relative humidity profiles and in achieving steady‐state liquid water content suggest the need for further investigation into the mechanisms driving these variations. Despite these discrepancies, the models generally reproduced consistent power‐law relationships between the microphysical variables. This model intercomparison underscores the importance of controlled cloud chamber experiments for validating and improving cloud microphysical parameterizations. Recommendations for future modeling studies are also highlighted, including constraining wall conditions and processes, investigating droplet/aerosol removal (including sidewall losses), and conducting simplified experiments to isolate specific processes contributing to model divergence and reduce model uncertainties.

54 ENVIRONMENTAL SCIENCES↗

Volumetrically Absorbing Thermal Insulator (VATI) for High-temperature Receivers

This project seeks to exploit volumetric absorption of concentrated solar irradiation in a thin, high-temperature, thermally conductive medium. High efficiency thermal conversion is achieved by volumetric absorption of both incoming solar irradiation (short wavelength) and emitted irradiation (long wavelength), demonstrating a volumetrically absorbing thermally insulating (VATI) effect. The thermalized energy is conducted to the back wall, enclosing the working fluid (molten salt or supercritical CO 2 ). Reliance on conductive thermal transport requires a thermally conductive medium, in the absence of which large temperature gradients drive losses due to emission. The project’s outcomes are important to realize a cost-effective approach to reduce optical and thermal losses from CSP receivers at high temperatures (720°C). High-temperature stable and commercially available porous SiC structures (open-cell foams) were explored as a volumetrically absorbing and thermally insulating layer. To the best of our knowledge, the current state-of-the-art for industrially deployed coatings is Pyromark. However, Pyromark suffers degradation at 700+°C and diurnal temperature changes. This necessitates periodic recoating, and the downtime increases the overall levelized cost of energy (LCOE) production. On the other hand, contemporary research activities have generated significant advances in the development of selective emitter coatings, which present challenges with costs, scalability, and stability. The pursued approach alleviates these concerns by developing a receiver that utilizes the inherent structure of high-temperature stable porous materials to enable robust and cost-effective receivers which require no periodic maintenance downtimes. The overall goal of the project is to experimentally demonstrate a figure of merit (FOM) of 0.92 at a temperature of 720°C and solar irradiation of 1000x concentration with porous receivers. The relevant crystallographic (phase) optical and thermal properties of porous SiC were first characterized. Second, the 3-D geometry of the porous SiC was analyzed using micro-X-ray tomography and converted to CAD data using image processing analysis. Utilizing this 3-D geometry and relevant optical/thermal properties, Monte Carlo- Ray Tracing (MCRT) analysis was performed to extract important parameters governing solar-thermal energy conversion such as extinction coefficient (β, 1/m), scattering albedo (ω) and the scattering phase function (Φ). These properties were then integrated into an in-house radiative and conduction transport model to solve for temperature and transport fluxes characterizing the solar-thermal energy conversion. This model was utilized to predict the FOM for various SiC porous geometries and to identify the highest possible FOM. Testing of the optimized porous structures will be accomplished with a custom-built high-accuracy (< ±4%) FOM measurement test-stand and a 1000x solar concentrator. When neglecting convective losses and resistance at the open boundary and the back wall, the optimized SiC foam leads to a FOM of 0.84. This FOM does not surpass the performance of Pyromark 2500. Yet, conversely to Pyromark 2500, SiC is stable at high temperatures and does not degrade over time. Also, it is possible to boost the FOM of SiC by engineering its effective thermal conductivity and its scattering albedo. A FOM of 0.92 is predicted for an optimized foam by accounting for convective losses and resistance at the open boundary and the back wall. This largely exceeds the FOM of Pyromark 2500 (~0.86) predicted by neglecting convective losses and resistance. Therefore, further engineering of the foam could lead to unprecedented FOMs.

14 SOLAR ENERGY↗

NRC Multiphysics Analysis Capability Deployment FY21: Part 3

This report details the progress and activities of Idaho National Laboratory (INL) on the Nuclear Regulatory Commission (NRC) project “Development and Modeling Support for Advanced Non-Light Water Reactors.” The deliverables completed for this report are: Deliverable 1c: the capability to model gas mixtures was added to Pronghorn. A test problem mimicking the conditions achieved in a depressurized loss of forced cooling (DLOFC) event was solved with both RELAP-5 and Pronghorn. Pronghorn employed a finite vol ume method with the Kurganov-Tadmor discretization. The comparison between the mass fraction spatial profiles computed with RELAP-5 and Pronghorn clearly shows the presence of numerical artifacts (i.e., overly diffusive behavior at low Mach numbers). We confirmed that the problem disappears at higher Mach numbers. We recommend future work on the implementation of a low Mach finite volume formulation to better treat low Mach number problems. Deliverable 2a: we demonstrated two approaches to model the radiation/conduction/natural convection heat transfer across a stagnant gas for the PBMR-400 design using Pronghorn. The first approach is based on the net radiation method, which relies on the computation of view factors with the Multiphysics Object-Oriented Simulation Environment (MOOSE) ray tracing capability. The second method is a traditional thermal resistance approach. The test problems include both 2D and 3D geometries. In all cases, the results show very good agreement during a DLOFC transient. This confirms that the faster thermal resistance method produces solutions that are equivalent to the net radiation method for this geometry. Deliverable 3d: we demonstrated the use of the advection kernel for the delayed neutron precursor equation in Griffin with a 2D MSFR model. The results appear physical but further verification is recommended. We also recommend the addition of conjugate heat transfer to compute the temperatures and model the thermomechanic behavior of the reflectors and other structures. Significant memory and performance issues were encountered in the 3D axisymmetric model. Future work is recommended in this area. Task 8g: this task allows multidimensional MOOSE applications to be coupled to system codes (RELAP-7 and SAM). We implemented a faster multiphysics iteration coupling algorithm, which provides an overall 6× acceleration of the 3D-1D coupling of the core multidi- mensional fluid flow solver and the 1D primary and secondary loop model.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

A Concept of a Convection–Cloud Chamber to Study Aerosol–Cloud–Drizzle Interactions

Understanding and quantifying the full chain of processes from aerosol activation to drizzle formation, and the associated feedbacks to the aerosol chemical and physical properties, all within a turbulent cloud are some of the toughest challenges in atmospheric chemistry and physics and are keys to the cloud–precipitation puzzle. This paper describes a concept for a new type of research facility consisting of a cloud chamber plus associated instrumentation and computational models, to explore aerosol–cloud interactions and processing, cloud optical properties, entrainment–cloud interactions, and quantitative assessment of drizzle onset. The envisioned design is for a 3 m × 3 m × 9 m chamber, such that the height is sufficient to achieve long lifetimes for aerosol processing and for significant drizzle growth by collision and coalescence. A suite of computational tools for simulating microphysical properties in the chamber provides a digital twin for designing the chamber and a range of example experiments. Theory and test results from novel remote sensing systems for exploring chemical and physical interactions and evolution of aerosols, cloud droplets, and drizzle within turbulent clouds are described. Testing of technology needed for the operation of a large-volume chamber, including aerosol generation methods and novel materials for water vapor boundary conditions, is described. Simulations suggest that spatially uniform turbulence and microphysical properties can be sustained in a steady state, with reasonable aerosol and water vapor fluxes, and that substantial drizzle can be produced through collision and coalescence of cloud droplets. Remaining challenges for more detailed engineering design and a discussion of possible first-light experiments are described.

54 ENVIRONMENTAL SCIENCES↗

Mesoscale Convective Systems Tracking Method Intercomparison (MCSMIP): Application to DYAMOND Global km‐Scale Simulations

Abstract Global kilometer‐scale models represent the future of Earth system modeling, enabling explicit simulation of organized convective storms and their associated extreme weather. Here, we comprehensively evaluate tropical mesoscale convective system (MCS) characteristics in the DYAMOND (DYnamics of the atmospheric general circulation modeled on non‐hydrostatic domains) simulations for both summer and winter phases. Using 10 different feature trackers applied to simulations and satellite observations, we assess MCS frequency, precipitation, and other key characteristics. Substantial differences (a factor of 2–3) arise among trackers in observed MCS frequency and their precipitation contribution, but model‐observation differences in MCS statistics are more consistent across trackers. DYAMOND models are generally skillful in simulating tropical mean MCS frequency, with multi‐model mean biases ranging from −2%–8% over land and −8%–8% over ocean (summer vs. winter). However, most DYAMOND models underestimate MCS precipitation amount (23%) and their contribution to total precipitation (17%). Biases in precipitation contributions are generally smaller over land (13%) than over ocean (21%), with moderate inter‐model variability. While models better simulate MCS diurnal cycles and cloud shield characteristics, they overestimate MCS precipitation intensity and underestimate stratiform rain contributions (up to a factor of 2), particularly over land, albeit observational uncertainties exist. Additionally, models exhibit a wide range of precipitable water in the tropics compared to reanalysis and satellite observations, with many models showing exaggerated sensitivity of MCS precipitation intensity to precipitable water. The MCS metrics developed here provide process‐oriented diagnostics to guide future model development.

54 ENVIRONMENTAL SCIENCES↗