Search NASA⌕ Search

SEARCH · Search NASA

Results for “conditional generative models”

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 163 records · Page 9

MECHANICAL SHOCK AND VIBRATION ANALYSIS OF SPENT NUCLEAR FUEL CARRIED BY THE ATLAS RAILCAR

Researchers at Pacific Northwest National Laboratory have completed a structural-dynamic analysis of spent nuclear fuel subjected to the mechanical shock and vibration environment that is anticipated during normal conditions of transport in casks carried by the Atlas railcar. The Atlas railcar is a new railcar design that is being developed specifically for the purpose of carrying spent nuclear fuel casks. The analysis used best-estimate railcar dynamics models of the Atlas railcar and considered 17 different spent nuclear fuel transportation cask systems, representing the current fleet of cask options. This work used NUCARS, a specialized railcar dynamics explicit finite element code to calculate railcar dynamic response to prescribed speeds and track configurations. The railcar dynamics models provided cask transient motion for a wide range of speeds and track conditions, generating a relatively large database of potential cask motion. All of the cask motion transients were then applied as loading conditions to LS-DYNA structural-dynamic models of a single fuel rod. The analyses predict that the Equipos Nucleares S.A./U.S. Department of Energy (ENSA/DOE) multimodal transportation test of 2017 provided a relatively stronger vibration environment than is expected from the Atlas railcar. This paper describes the analysis methods, the analysis results, and compares the results of the Atlas transportation analysis to the test results and analyses of the ENSA/DOE multimodal transportation test of 2017.

Klymyshyn, Nicholas A.↗

A superstructure representation, generation, and modeling framework for chemical process synthesis

We present a framework for the efficient representation, generation, and modeling of superstructures for process synthesis. First, we develop a new representation based on three basic elements: units, ports, and conditioning streams. Second, we present four rules based on “minimal” and “feasible” component sets for the generation of simple superstructures containing all feasible embedded processes. Third, in terms of modeling, we develop a modular approach, and formulate models for each basic element. We also present a canonical form of element models using input/output variables and constrained/free variables. The proposed methods provide a coherent framework for superstructure‐based process synthesis, allowing efficient model generation and modification. © 2016 American Institute of Chemical Engineers AIChE J , 62: 3199–3214, 2016

Wu, WenZhao↗

Micromechanics of Fracture Propagation During Multistage Stress Relaxation and Creep in Brittle Rocks

Time-dependent rock deformation caused by the initiation and growth of fractures leads to the weakening of the rock mass. Understanding the fracturing mechanisms involved in the time-dependent behavior in brittle rocks is very important and to achieve this goal, a systematic series of three types of experiments was performed on double-flawed prismatic Barre granite specimens under unconfined compression. The first series aimed to identify the failure mechanism in the short-term failure mode under monotonic loading, the, second series involved multistage relaxation (constant strain) experiments to analyze the damage at different strain levels, and the third series explored the fracture propagation under multistage creep (constant load) experiments. The spatial and temporal evolution of cracking mechanisms were evaluated using the acoustic emission (AE) and two-dimensional digital image correlation (2D-DIC) techniques to observe the whole crack growth process as well as the accumulated inelastic strain at the specified region of interest. Results suggest that in the case of multistage creep experiments, the time to failure was less compared to the multistage relaxation, when loaded above the crack damage threshold (CD) estimated from the monotonic testing. Furthermore, the frequency magnitude distribution of the AE events generated in the three loading conditions followed the Gutenberg Richter model. A relatively lower b-value was obtained for the creep experiments, indicative of high energy AE events and faster crack growth. In addition, the AE and DIC results also revealed high evolution of tensile cracks at-different stages of creep and relaxation compared to shear and mixed-mode cracks.

58 GEOSCIENCES↗

Considering computational speed vs. accuracy: Choosing appropriate mesoscale RVE boundary conditions

Modeling a material’s microstructure using continuum theories allows for inspection of the relationship between coarse scale and fine scale behaviors. Computational limits generally require selection of a sub-volume from a bulk sample in order to directly model the microstructure. Boundary conditions are applied to the sub-volume to mimic the excluded bulk material. Appropriate selection of boundary conditions helps effectively determine the appropriate spatial scale required of the sub-volume. Applicable boundary conditions include direct displacement, periodic, and uniform traction. While direct displacement and periodic boundary conditions are commonly used, uniform traction boundary conditions have seen limited use due to rigid body stability issues in simulations of compression or shear deformation. A new application of uniform traction boundary conditions was developed through linear constraint equations, similar to approaches employed by direct displacement and periodic boundary conditions, to quench rigid body motions with minimal interference of the relative deformation of the model. These boundary conditions were tested by compressing several synthetically generated periodic microstructures using the finite element method. Evaluating the effective stiffness along the compression axis, the direct displacement boundary condition produced the stiffest response, whereas the uniform traction boundary condition produced the most compliant. Periodic boundary conditions produced the same response for all volumes analyzed and both the direct displacement and uniform traction boundary conditions trended toward the periodic response as the domain volume increased. Computational performance was also evaluated for each boundary condition using implicit and explicit solvers. Direct displacement boundary conditions presented the lowest computational cost of all of the boundary conditions followed by periodic then uniform traction. The computational expense of periodic and uniform traction boundary conditions limited the viable spatial scale and mesh resolutions able to be simulated. Selection of appropriate boundary conditions for specific uses need to be a balance between allowable computational expense and accuracy of the method. Techniques for evaluating which boundary conditions to use are discussed.

42 ENGINEERING↗

Current data are consistent with flat spatial hypersurfaces in the Λ CDM cosmological model but favor more lensing than the model predicts

Here, we study the performance of three pairs of tilted, and a pair of untilited, ΛCDM cosmological models, with three of these four pairs allowing for non-flat spatial hypersurfaces, against cosmic microwave background (CMB) temperature and polarization power spectrum data (P18), measurements of the Planck 2018 lensing potential power spectrum (lensing), and a large compilation of non-CMB data (non-CMB). For the eight models, we measure cosmological parameters and study whether or not pairs of the data sets (as well as subsets of them) are mutually consistent in these models. Half of these models allow the lensing consistency parameter A L , which re-scales the gravitational potential power spectrum, to be an additional free parameter to be determined from data, while the other three have A L = 1 which is the theoretically expected value. The pair of untilted non-flat ΛCDM models are incompatible with P18 data. The tilted spatially-flat models assume the usual primordial spatial inhomogeneity power spectrum that is a power law in wave number. The tilted non-flat models assume either the primordial power spectrum used in the Planck group anal yses [Planck P(q)], that has recently been numerically shown to be a good approximation to what is quantum-mechanically generated from a particular choice of closed inflation model initial conditions, or a recently computed power spectrum [new P(q)] that quantum-mechanically follows from a different set of non-flat inflation model initial conditions. In the tilted non-flat models with A L = 1 we find differences between P18 data and non-CMB data cosmological parameter constraints, which are large enough to rule out the Planck P(q) model at 3σ but not the new P(q) model. No significant differences are found when cosmological parameter constraints obtained with two different data sets are compared within the standard tilted flat ΛCDM model. While both P18 data and non-CMB data separately favor a closed geometry, with spatial curvature density parameter Ω k < 0, when P18+non-CMB data are jointly analyzed the evidence in favor of non-flat hypersurfaces subsides. Differences between P18 data and non-CMB data cosmological constraints subside when A L is allowed to vary. From the most restrictive P18+lensing+non-CMB data combination we get almost model-independent constraints on the cosmological parameters and find that the A L > 1 option is preferred over the Ω k < 0 one, with the A L parameter, for all models, being larger than unity by ~ 2.5σ. According to the deviance information criterion, in the P18+lensing+non-CMB analysis, the varying A L option is on the verge of being strongly favored over the A L = 1 one, which could indicate a problem for the standard tilted flat ΛCDM model. These data are consistent with flat spatial hypersurfaces but more and better data could improve the constraints on Ω k and might alter this conclusion. Error bars on some cosmological parameters are significantly reduced when non-CMB data are used jointly with P18+lensing data. For example, in the tilted flat ΛCDM model for P18+lensing+non-CMB data the Hubble constant H 0 = 68.09 ± 0.38 km s -1 Mpc -1 , which is consistent with that from a median statistics analysis of a large compilation of H 0 measurements as well as with a number of local measurements of the cosmological expansion rate. This H 0 error bar is 31% smaller than that from P18+lensing data alone.

79 ASTRONOMY AND ASTROPHYSICS↗

Optimal Operation of a Plug-In Hybrid Vehicle

In this work, we present a convex optimization control method that has been shown in simulations to increase the fuel efficiency of a plug-in hybrid electric vehicle by over 10%. Using information on energy demand and energy use profiles, the problem is defined to preferentially use battery resources sourced from the grid over petroleum resources. We pose the general nonlinear optimal resource management problem over a predetermined route as a convex optimization problem using a reduced model of the vehicle. This problem is computationally efficient enough to be optimized “on the fly” on the on-board vehicle computer and is thus able to adapt to changing vehicle conditions in real time. Lastly, using this reduced model to generate control inputs for the detailed vehicle simulator autonomie, we record efficiency gains of over 10% as compared to the industry standard charge depleting charge sustaining controller over synthetic mixed urban-suburban routes.

02 PETROLEUM↗

Moving beyond the Aerosol Climatology of WRF-Solar: A Case Study over the North China Plain

Numerical weather prediction (NWP), when accessible, is a crucial input to short-term solar power forecasting. WRF-Solar, the first NWP model specifically designed for solar energy applications, has shown promising predictive capability. Nevertheless, few attempts have been made to investigate its performance under high aerosol loading, which attenuates incoming radiation significantly. The North China Plain is a polluted region due to industrialization, which constitutes a proper testbed for such investigation. Here, in this paper, aerosol direct radiative effect (DRE) on three surface shortwave radiation components (i.e., global, beam, and diffuse) during five heavy pollution episodes is studied within the WRF-Solar framework. Results show that WRF-Solar overestimates instantaneous beam radiation up to 795.3 W m -2 when the aerosol DRE is not considered. Although such overestimation can be partially offset by an underestimation of the diffuse radiation of about 194.5 W m -2 , the overestimation of the global radiation still reaches 160.2 W m -2 . This undesirable bias can be reduced when WRF-Solar is powered by Copernicus Atmosphere Monitoring Service (CAMS) aerosol forecasts, which then translates to accuracy improvements in photovoltaic (PV) power forecasts. This work also compares the forecast performance of the CAMS-powered WRF-Solar with that of the European Centre for Medium-Range Weather Forecasts model. Under high aerosol loading conditions, the irradiance forecast accuracy generated by WRF-Solar increased by 53.2% and the PV power forecast accuracy increased by 6.8%.

54 ENVIRONMENTAL SCIENCES↗

Integrating Carbon Capture, Utilization, & Sequestration into Chemical Pulp Mills

The U.S. pulp and paper industry presents a unique and largely untapped opportunity for large- scale carbon dioxide removal (CDR). Unlike most industrial sectors, pulp mills rely heavily on biomass, meaning that much of their carbon emissions originate from atmospheric CO₂ that was recently captured by plants. If this biogenic CO₂ can be captured and permanently stored, pulp mills can be transformed from carbon emitters into net carbon removal facilities. This project was motivated by that opportunity and aimed to develop and evaluate integrated, low-cost strategies for capturing, utilizing, and sequestering CO₂ within existing chemical pulping operations. The scope of this work focused on four complementary innovations designed to integrate seamlessly into kraft pulp mill infrastructure: (1) in situ CO₂ capture within the recovery cycle, (2) oxy-fuel retrofitting of the rotary lime kiln to produce a high-purity CO₂ stream, (3) ex situ CO₂ capture and mineralization using pulp mill residues (dregs, grits, and lime mud), and (4) beneficial reuse of these residues as mineral carbonate fertilizers. The project combined process modeling, laboratory experimentation, life cycle assessment (LCA), and field trials to evaluate the technical feasibility, economic viability, and environmental impact of these approaches. The results demonstrate that pulp mills can serve as effective platforms for carbon removal when equipped with integrated carbon capture systems. Process modeling showed that combining sodium spiking with oxy-fuel calcination significantly enhances CO₂ capture efficiency while reducing costs by up to 31% compared to conventional configurations. Experimental work further revealed that calcination behavior in high-CO₂ environments differs substantially from traditional systems, leading to the development of a new kinetic model that predicts reaction rates under these conditions. This model provides essential design guidance for next-generation decarbonized lime kilns. In parallel, the project demonstrated that alkaline mineral residues generated during pulping operations can be repurposed as a sustainable alternative to agricultural lime. Across a wide range of soils in the southeastern United States, these materials performed equivalently to commercial lime in adjusting soil pH while offering lower greenhouse gas emissions and reduced cost. Field and greenhouse studies confirmed that crop and tree growth responses were comparable, supporting their viability as a drop-in replacement. This co-product pathway provides a practical utilization strategy that offsets costs and improves overall system economics. A major contribution of this project is the first comprehensive life cycle assessment of carbon removal in pulp and paper systems across multiple system boundaries. Results show that retrofitted mills can achieve carbon removal efficiencies ranging from 12% to 92%, depending on how the system is defined. This finding highlights a critical issue in carbon accounting: reported performance is highly sensitive to methodological choices. By explicitly quantifying these differences, this work provides valuable guidance for policymakers, carbon registries, and project developers working to standardize carbon removal metrics. From a commercialization perspective, the technologies investigated in this project are well- aligned with existing industrial infrastructure, minimizing the need for entirely new facilities. 3 DE-EE0009413 Industry engagement throughout the project—including collaboration with pulp and paper companies, equipment manufacturers, and carbon removal developers—has accelerated the transition from research to deployment. Notably, a commercial developer is actively pursuing carbon capture projects at pulp mills in the southeastern United States and has cited this research as a contributing foundation. The emergence of voluntary carbon markets and long-term offtake agreements further strengthens the business case for implementation. The broader public benefits of this work are significant. By enabling large-scale carbon removal using existing industrial systems, this approach offers a near-term pathway to reduce atmospheric CO₂ concentrations while supporting domestic manufacturing and rural economies. The reuse of industrial residues as fertilizers reduces reliance on mined materials, lowers costs for farmers, and decreases environmental impacts associated with conventional lime production. In addition, the project has supported workforce development by training graduate students and researchers in carbon capture technologies, helping to build capacity in a critical area of national interest. In conclusion, this project demonstrates that integrated carbon capture, utilization, and sequestration in pulp mills is both technically feasible and economically promising. By combining process innovation, experimental validation, and systems-level analysis, the work advances the understanding of how biomass-based industries can contribute to climate mitigation. The findings provide a strong foundation for commercial deployment and offer a scalable solution for transforming a major U.S. industry into a source of durable carbon removal.

09 BIOMASS FUELS↗

FUEL PERFORMANCE MODELING OF FECRAL CLADDING FAILURE DURING LOSS OF COOLANT CONDITIONS

Accident tolerant materials are defined as those that provide significantly increased response time in the event of an accident while providing similar or improved performance as the conventional UO2/Zircaloy-4 fuel rods during normal operation. One avenue for materials to demonstrate enhanced accident tolerance is to have improved reaction (e.g., oxidation) kinetics with steam, resulting in slower hydrogen (or other combustible gases) generation rate, while maintaining acceptable thermo-mechanical properties, fuel-clad interactions, and fission product behavior. Due to their improved reaction kinetics with steam compared to Zircaloy-4, iron-chromium-aluminum (FeCrAl) alloys are being considered as an accident tolerant cladding material. In order to determine whether FeCrAl alloys are indeed accident tolerant materials, comparisons between conventional zirconium-based claddings under accident conditions is required. To this end, we have developed a burst (failure) criterion for FeCrAl using existing experimental data for first-generation alloys. This model was incorporated into the BISON fuel performance code to perform comparative studies with Zircaloy-4 cladding tubes under non-oxidizing and oxidizing LOCA conditions. The separate effects LOCA conditions considered were adopted from the PUZRY and REBEKA test series. Due to the fact that FeCrAl alloys have a larger thermal absorption cross-section and higher stiffness than Zircaloy-4, cladding thickness will be smaller for FeCrAl than for Zircaloy-4. Three different cladding thicknesses were analyzed in this work. The results indicate that for both non-oxidizing (PUZRY) and oxidizing (REBEKA) LOCA conditions, FeCrAl has similar burst behavior (i.e., time to burst and pressure at burst) as Zircaloy-4. It should be noticed that this study used data for first-generation FeCrAl alloys to develop the burst criterion. Although no burst data?currently exists for C35M (a second-generation FeCrAl alloy under development at Oak Ridge National Laboratory), improved burst behavior is anticipated. As burst data for C35M becomes available the burst criterion will be updated and the simulations revisited to determine the accident tolerance of second-generation FeCrAl alloys.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Fuel performance modeling of FeCrAl cladding failure during loss of coolant conditions

Accident tolerant materials are defined as those that provide significantly increased response time in the event of an accident while providing similar or improved performance as the conventional UO2/Zircaloy-4 fuel rods during normal operation. One avenue for materials to demonstrate enhanced accident tolerance is to have improved reaction (e.g., oxidation) kinetics with steam, resulting in slower hydrogen (or other combustible gases) generation rate, while maintaining acceptable thermo-mechanical properties, fuel-clad interactions, and fission product behavior. Due to their improved reaction kinetics with steam compared to Zircaloy-4, iron-chromium-aluminum (FeCrAl) alloys are being considered as an accident tolerant cladding material. In order to determine whether FeCrAl alloys are indeed accident tolerant materials, comparisons between conventional zirconium-based claddings under accident conditions is required. To this end, we have developed a burst (failure) criterion for FeCrAl using existing experimental data for first-generation alloys. This model was incorporated into the BISON fuel performance code to perform comparative studies with Zircaloy-4 cladding tubes under non-oxidizing and oxidizing LOCA conditions. The separate effects LOCA conditions considered were adopted from the PUZRY and REBEKA test series. Due to the fact that FeCrAl alloys have a larger thermal absorption cross-section and higher stiffness than Zircaloy-4, cladding thickness will be smaller for FeCrAl than for Zircaloy-4. Three different cladding thicknesses were analyzed in this work. The results indicate that for both non-oxidizing (PUZRY) and oxidizing (REBEKA) LOCA conditions, FeCrAl has similar burst behavior (i.e., time to burst and pressure at burst) as Zircaloy-4. It should be noticed that this study used data for first-generation FeCrAl alloys to develop the burst criterion. Although no burst data?currently exists for C35M (a second-generation FeCrAl alloy under development at Oak Ridge National Laboratory), improved burst behavior is anticipated. As burst data for C35M becomes available the burst criterion will be updated and the simulations revisited to determine the accident tolerance of second-generation FeCrAl alloys.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Probabilistic wave parameters for WEC survival analysis at US Navy's Wave Energy Test Site in Hawaii.

The persistent and diversified wave conditions make the Navy's Wave Energy Test Site (WETS) at Marine Corps Base Hawaii, a suitable place for WEC testing. Hurricane waves, despite their infrequent occurrence, pose a challenge to the survivability of the devices. We utilized 252 hurricane scenarios within 200 km of WETS from downscaling simulations of a climate model and the third generation spectral model WaveWatch III to produce a range of severe wave conditions at WETS. The computed wave height presented as a function of annual exceedance probability enables a risk-based survival analysis of WECs. This study provides a proof-of-concept of the probabilistic approach and a baseline for analysis of hurricane wave conditions under the present and future climate projections.

16 TIDAL AND WAVE POWER↗

Development of a Dynamic Model and Control System for Load-Following Studies of Supercritical Pulverized Coal Power Plants

Traditional energy production plants are increasingly forced to cycle their load and operate under low-load conditions in response to growth in intermittent renewable generation. A plant-wide dynamic model of a supercritical pulverized coal (SCPC) power plant has been developed in the Aspen Plus Dynamics® (APD) software environment and the impact of advanced control strategies on the transient responses of the key variables to load-following operation and disturbances can be studied. Models of various key unit operations, such as the steam turbine, are developed in Aspen Custom Modeler® (ACM) and integrated in the APD environment. A coordinated control system (CCS) is developed above the regulatory control layer. Three control configurations are evaluated for the control of the main steam; the reheat steam temperature is also controlled. For studying servo control performance of the CCS, the load is decreased from 100% to 40% at a ramp rate of 3% load per min. The impact of a disturbance due to a change in the coal feed composition is also studied. The CCS is found to yield satisfactory performance for both servo control and disturbance rejection.

20 FOSSIL-FUELED POWER PLANTS↗

Improving Seasonal Forecast Using Probabilistic Deep Learning

The path toward realizing the potential of seasonal forecasting and its socioeconomic benefits relies on improving general circulation model (GCM) based dynamical forecast systems. To improve dynamical seasonal forecasts, it is crucial to set up forecast benchmarks, and clarify forecast limitations posed by model initialization errors, formulation deficiencies, and internal climate variability. With huge costs in generating large forecast ensembles, and limited observations for forecast verification, the seasonal forecast benchmarking and diagnosing task proves challenging. Here, we develop a probabilistic deep learning-based statistical forecast methodology, drawing on a wealth of climate simulations to enhance seasonal forecast capability and forecast diagnosis. By explicitly modeling the internal climate variability and GCM formulation differences, the proposed Conditional Generative Forecasting (CGF) methodology enables bypassing crucial barriers in dynamical forecast, and offers a top-down viewpoint to examine how complicated GCMs encode the seasonal predictability information. We apply the CGF methodology for global seasonal forecast of precipitation and 2 m air temperature, based on a unique data set consisting 52,201 years of climate simulation. Results show that the CGF methodology can faithfully represent the seasonal predictability information encoded in GCMs. We successfully apply this learned relationship in real-world seasonal forecast, achieving competitive performance compared to dynamical forecasts. Using this CGF as benchmark, we reveal the impact of insufficient forecast spread sampling that limits the skill of the considered dynamical forecast system. Finally, we introduce different strategies for composing ensembles using the CGF methodology, highlighting the potential for leveraging the strengths of multiple GCMs to achieve advantgeous seasonal forecast.

54 ENVIRONMENTAL SCIENCES↗

Stochastic Model Predictive Control With Gaussian Wind Direction Preview for Wake Steering

This article addresses the problem of wake steering control for wind farms that explicitly consider the tradeoff between farm-level power generation and yaw duty cycle under variable and uncertain wind conditions. A novel stochastic model predictive control (MPC) algorithm is presented, which utilizes a stochastic model of the freestream wind field components in a receding horizon framework to compute optimal yaw set points that maximize the expected value of the farm power while constraining the yaw actuation. Different configurations of the algorithm are evaluated using a steady-state wind farm simulator. The proposed stochastic MPC algorithm can plan control actions over a future prediction horizon based on probabilistic estimates of the incoming wind magnitude and direction.

17 WIND ENERGY↗

Ensemble cell-wide kinetic modeling of anaerobic organisms to support fuels and chemicals production

The goal of this project is to understand the metabolism and underlying regulatory networks of two anaerobic microorganisms, Moorella thermoacetica and Clostridium thermocellum by making use of Ensemble Modeling (EM) paradigm through integration of multiple omic information (transcriptomic, proteomic, metabolomic & fluxomic). Kinetic models have the potential to provide a mechanistic description of cellular metabolism by combining several layers of regulatory events into an integrated framework. The foundation of the kinetic model is a stoichiometric network description. To this end, we have developed and curated genome-scale metabolic (GSM) models of the organisms that will serve as scaffold for building the kinetic models. For M. thermoacetica, the GSM model was constructed and tested by the Stephanopoulos group with experimental growth rate. The GSM elucidated the energy conservation mechanism under autotrophic growth conditions. For C. thermocellum, the Maranas group has developed the second-generation GSM model through integration of experimental data. The Maranas group has also constructed a core metabolic kinetic model of the organism which includes the central metabolism and all relevant pathways for biofuels production. The stoichiometric and kinetic model predictions were tested extensively using experimental data for wild-type and 20 mutant strains. We have outlined below in detail the progress related to specific individual aims. Below are listed publications of research carried out based on this research grant support.

09 BIOMASS FUELS↗

A 30-yr high-resolution weather research and forecasting model downscaling data over California and Nevada

This dataset presents a 30-year high resolution meteorological dataset obtained using the WRF model (Advanced version Research WRF version 4.4). We used WRF and European Centre for Medium-Range Weather Forecasts Reanalysis v5 as initial and boundary conditions to generate gridded meteorological variables. A large number of surface weather stations was used for model validation. A multi-physics analysis was first developed to identify a good physics suite extended from 6 November 00 UTC to 10 November 23 UTC, 2018, which included the Camp Fire in northern California. Based on the best physics suite, the downscaling dataset extends from 1 December to 28 February, 1990–2021 and the horizontal domain has 1.5 km grid spacing covering the entire states of California and Nevada in the United States. Comparisons between hourly surface observations and WRF simulations of air temperature, relative humidity and wind speeds show mean absolute errors on the order of (1.6-2.0 C), (10 %) and 1.2–1.5 m s -1 , respectively.

54 ENVIRONMENTAL SCIENCES↗

Modeling of Aircraft Measurements of Ice Crystal Concentration in the Arctic and a Parameterization for Mixed-Phase Cloud

Abstract In this study, two parameterizations of ice nucleation rate on dust particles are used in a parcel model to simulate aircraft measurements of ice crystal number concentration Ni in the Arctic. The parcel model has detailed microphysics for droplet and ice nucleation, growth, and evaporation with prescribed vertical air velocities. Three dynamic regimes are considered, including large-scale ascent, cloud-top generating cells, and their combination. With observed meteorological conditions and aerosol concentrations, the parcel model predicts the number concentrations of size-resolved ice crystals, which may be compared to aircraft measurements. Model results show rapid changes with height/time in relative humidity, Ni, and thermodynamic phase partitioning, which is not resolved in current climate and weather forecasting models. Parameterizations for ice number and nucleation rate in mixed-phase stratus clouds are thus developed based on the parcel model results to represent the time-integrated effect of some microphysical processes in large-scale models.

54 ENVIRONMENTAL SCIENCES↗

Evaluating the Practicability of Multiphase CFD for BWR Fuel Assembly Analysis

Multiphase computational fluid dynamics (M-CFD) models of a boiling water reactor (BWR) fuel assembly have been generated in the commercial CFD code, STAR-CCM+ Ver.11.04.010-R8. The Eulerian multiphase approach was employed, in combination with a baseline ‘Zero Closure model’ for two-phase flow boiling, applicable to high void fraction regimes. The capability and practicability of MCFD for BWR fuel assembly were evaluated against the international OECD/NRC BWR Full-size Finemesh Bundle Test (BFBT) benchmark data. The benchmark selected, BFBT-4101-61, is a steady-state full power case, with high void fraction conditions. Geometrical modeling and mesh construction of BWR assemblies introduce challenges in generating high quality computational meshes, while retaining a low cell count to allow acceptable runtimes. In order to evaluate the scalability of the CFD code, parallel CFD simulations were carried out with the High Performance Computing (HPC) system of the Idaho National Laboratory. M-CFD results display good agreement with the experimental data for both local void fraction distribution and exit quality. The parallel speedup factor is linearly proportional to the number of processors up to 500 cores, and retains over 92% efficiency with 1000 cores. The execution time of the simulation decreased exponentially with increasing number of processors while the cost of the simulation increased linearly. The proposed baseline closure, in combination with the parallel performance of the STAR-CCM+ code demonstrated promising capabilities for analyzing BWR fuel assemblies.

97 MATHEMATICS AND COMPUTING↗