Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “parameterized dynamical systems”

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.

231 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

An open-source framework for balancing computational speed and fidelity in production cost models

Studies of bulk power system operations need to incorporate uncertainty and sensitivity analyses, especially around exposure to weather and climate variability and extremes, but this remains a computational modeling challenge. Commercial production cost models (PCMs) have shorter runtimes, but also important limitations (opacity, license restrictions) that do not fully support stochastic simulation. Open-source PCMs represent a potential solution. They allow for multiple, simultaneous runs in high-performance computing environments and offer flexibility in model parameterization. Yet, developers must balance computational speed (i.e. runtime) with model fidelity (i.e. accuracy). In this paper, we present Grid Operations (GO), a framework for instantiating open-source, scale-adaptive PCMs. GO allows users to search across parameter spaces to identify model versions that appropriately balance computational speed and fidelity based on experimental needs and resource limits. Results provide generalizable insights on how to navigate the fidelity and computational speed tradeoff through parameter selection. We show that models with coarser network topologies can accurately mimic market operations, sometimes better than higher-resolution models. It is thus possible to conduct large simulation experiments that characterize operational risks related to climate and weather extremes while maintaining sufficient model accuracy.

42 ENGINEERING↗
DOE OSTI2025

Symmetry breaking as predicted by a phase space Hamiltonian with a spin Coriolis potential

Here, we perform electronic structure calculations for a set of molecules with degenerate spin-dependent ground states ( 3 CH 2 , 2 CH$^{•}_{3}$, 3 O 2 ) going beyond the Born–Oppenheimer approximation and accounting for nuclear motion. According to a phase space approach that parameterizes electronic states (|Φ⟩) and electronic energies (E) by nuclear position and momentum [i.e., |Φ(R, P)⟩ and E(R, P)], we find that the presence of degenerate spin degrees of freedom leads to broken symmetry ground states. More precisely, rather than a single degenerate minimum at (R, P) = (R min , 0), the ground state energy has two minima at (R,P)=(R' min ,±P min ) (where R' min is close to R min ), dramatically contradicting the notion that the total energy of the system can be written in separable form as E = $\frac{P^2}{2M}$ + V el . Although we find that the broken symmetry solutions have small barriers between them for the small molecules, we hypothesize that the barriers should be macroscopically large for metallic solids, thus offering up a new phase-space potential energy surface for simulating the Einstein–de Haas effect.

Berry connection↗
DOE OSTI2020

Assessing impacts of selective logging on water, energy, and carbon budgets and ecosystem dynamics in Amazon forests using the Functionally Assembled Terrestrial Ecosystem Simulator

Abstract. Tropical forest degradation from logging, fire, and fragmentation not only alters carbon stocks and carbon fluxes, but also impacts physical land surface properties such as albedo and roughness length. Such impacts are poorly quantified to date due to difficulties in accessing and maintaining observational infrastructures, as well as the lack of proper modeling tools for capturing the interactions among biophysical properties, ecosystem demography, canopy structure, and biogeochemical cycling in tropical forests. As a first step to address these limitations, we implemented a selective logging module into the Functionally Assembled Terrestrial Ecosystem Simulator (FATES) by mimicking the ecological, biophysical, and biogeochemical processes following a logging event. The model can specify the timing and aerial extent of logging events, splitting the logged forest patch into disturbed and intact patches; determine the survivorship of cohorts in the disturbed patch; and modifying the biomass and necromass (total mass of coarse woody debris and litter) pools following logging. We parameterized the logging module to reproduce a selective logging experiment at the Tapajós National Forest in Brazil and benchmarked model outputs against available field measurements. Our results suggest that the model permits the coexistence of early and late successional functional types and realistically characterizes the seasonality of water and carbon fluxes and stocks, the forest structure and composition, and the ecosystem succession following disturbance. However, the current version of FATES overestimates water stress in the dry season and therefore fails to capture seasonal variation in latent and sensible heat fluxes. Moreover, we observed a bias towards low stem density and leaf area when compared to observations, suggesting that improvements are needed in both carbon allocation and establishment of trees. The effects of logging were assessed by different logging scenarios to represent reduced impact and conventional logging practices, both with high and low logging intensities. The model simulations suggest that in comparison to old-growth forests the logged forests rapidly recover water and energy fluxes in 1 to 3 years. In contrast, the recovery times for carbon stocks, forest structure, and composition are more than 30 years depending on logging practices and intensity. This study lays the foundation to simulate land use change and forest degradation in FATES, which will be an effective tool to directly represent forest management practices and regeneration in the context of Earth system models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Incorporating the Effect of Large‐Scale Vertical Motion on Convection Through Convective Mass Flux Adjustment in E3SMv2

Abstract Recent observational studies suggest that the large‐scale dynamical forcing (vertical motion) plays important roles in deep convection development. In this study we propose a convective mass flux adjustment (MAdj) approach to represent the dynamical effects of large‐scale vertical motion on convection in the Department of Energy's Energy Exascale Earth System Model version 2 (E3SMv2). With MAdj, convection is enhanced (suppressed) when there is large‐scale ascending (descending) motion at the planetary boundary layer top. The coupling of convection with large‐scale circulation significantly improves the simulation of climate variability in E3SMv2 across multiple scales from the diurnal cycle, convectively coupled equatorial waves, to the Madden‐Julian Oscillation (MJO). The standard E3SMv2 tends to simulate overly weak diurnal amplitude of precipitation and overly weak variance of convectively coupled equatorial Kelvin and westward inertio‐gravity (WIG) waves. It also fails to simulate the essential characteristics of the MJO: continuous eastward propagation. With MAdj, the amplitude of diurnal cycle of precipitation is systematically increased and its probability density distribution is much closer to observations. The MAdj can also simulate more realistic eastward propagation of the MJO and much stronger convectively coupled Kelvin and WIG waves. Moreover, the MAdj approach slightly improves the climatology simulations in precipitation, cloud, radiation, circulation, temperature, and moisture fields, with overall root‐mean‐square error (RMSE) of major climatological fields reduced by about 2%. The MAdj approach suppresses excessive grid‐scale precipitation, reducing precipitation wet biases over South China Sea, Philippine Sea, Himalayas, and South Pacific Convergence Zone in western Pacific in summer.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Integrating New Arctic Plant Functional Types in a Land Surface Model Using Above- and Belowground Field Observations: Modeling Archive

This Modeling Archive is in support of an NGEE Arctic publication "Integrating New Arctic Plant Functional Types in a Land Surface Model Using Above- and Belowground Field Observations". We used field observations of biomass and leaf traits (Salmon et al., 2019b) across a gradient of plant communities on the Seward Peninsula in northwest Alaska (Breen et al., 2020; Iversen et al., 2017) to replace the original ELM configuration with nine, arctic-specific PFTs. Original arctic PFTs in the model were (1) broadleaf deciduous boreal shrub and (2) C3 arctic grass. New PFTs were (1) lichen, (2) bryophyte, (3) evergreen dwarf shrub, (4) deciduous dwarf shrub, (5) deciduous low shrub, (6) deciduous low to tall shrub, (7) alder shrub, (8) forb, and (9) graminoid. New PFTs were parameterized and evaluated using site-level trait and biomass measurements (Salmon et al., 2019b), root C:N ratios and rooting depth patterns (Salmon et al., 2019c), and root biomass (Salmon et al., 2020). This archive contains forcing data and model output from ELM simulations conducted at the Kougarok Hillslope site (Kougarok Road Mile Marker 64). Simulations were conducted for the Kougarok Hillslope site using meteorological driving data from the Scenarios Network for Alaska and Arctic Planning (SNAP) downscaled climate projection dataset using NCAR-CCSM forcing (Bieniek et al., 2020 dataset; Walsh et al., 2018). The archive includes three versions of model simulations representing a progression from the original model configuration to an Arctic-specific configuration (see Output section below). Simulation 1 used grid cell data from a global E3SM configuration, including the fractional area of E3SM PFTs assigned to the grid cell containing the Kougarok Hillslope site in global simulations. This simulation did not distinguish between different plant communities on the Kougarok Hillslope, but instead used a single point simulation to represent the entire area. Simulation 2 used default E3SM PFT definitions combined with adjusted depth to bedrock for plant communities with shallow rocky layers to represent the role of abiotic soil factors in driving site differences (based on depth measurements of Iversen et al., 2019a), and adjusted the relative areas of the broadleaf deciduous boreal shrub and C3 arctic grass PFTs to reflect the observed spatial coverage of shrub and non-shrub PFTs across Kougarok Hillslope plant communities. Simulation 3 used the new arctic PFT definitions and parameterizations based on the vegetation types present at the site, along with the adjusted soil depths used in Simulation 2. Included are *.pdf, *.nc (NetCDF), *tar.gz, *.py, and *.bash files. We would like to thank Mary's Igloo Native Corporation for allowing us to perform this research on their land and giving us the opportunity to share this research with their community. Thanks to Shawn Serbin and Alistair Rogers for helpful discussions. Holly Vander Stel, Joanne Childs, and Stan Wullschleger helped with collection of field data. The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research. The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska. Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Manifold learning for coarse-graining atomistic simulations: Application to amorphous solids

In this work, we introduce a generalized machine learning framework to probabilistically parameterize upper-scale models in the form of nonlinear PDEs consistent with a continuum theory, based on coarse-grained atomistic simulation data of mechanical deformation and flow processes. The proposed framework utilizes a hypothesized coarse-graining methodology with manifold learning and surrogate-based optimization techniques. Coarse-grained high-dimensional data describing quantities of interest of the multiscale models are projected onto a nonlinear manifold whose geometric and topological structure is exploited for measuring behavioral discrepancies in the form of manifold distances. A surrogate model is constructed using Gaussian process regression to identify a mapping between stochastic parameters and distances. Derivative-free optimization is employed to adaptively identify a unique set of parameters of the upper-scale model capable of rapidly reproducing the system's behavior while maintaining consistency with coarse-grained atomic-level simulations. The proposed method is applied to learn the parameters of the shear transformation zone (STZ) theory of plasticity that describes plastic deformation in amorphous solids as well as coarse-graining parameters needed to translate between atomistic and continuum representations. We show that the methodology is able to successfully link coarse-grained microscale simulations to macroscale observables and achieve a high-level of parity between the models across scales.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Quantifying and attributing time step sensitivities in present-day climate simulations conducted with EAMv1

Abstract. This study assesses the relative importance of time integration error in present-day climate simulations conducted with the atmosphere component of the Energy Exascale Earth System Model version 1 (EAMv1) at 1∘ horizontal resolution. We show that a factor-of-6 reduction of time step size in all major parts of the model leads to significant changes in the long-term mean climate. Examples of changes in 10-year mean zonal averages include the following: up to 0.5 K of warming in the lower troposphere and cooling in the tropical and subtropical upper troposphere, 1 %–10 % decreases in relative humidity throughout the troposphere, and 10 %–20 % decreases in cloud fraction in the upper troposphere and decreases exceeding 20 % in the subtropical lower troposphere. In terms of the 10-year mean geographical distribution, systematic decreases of 20 %–50 % are seen in total cloud cover and cloud radiative effects in the subtropics. These changes imply that the reduction of temporal truncation errors leads to a notable although unsurprising degradation of agreement between the simulated and observed present-day climate; to regain optimal climate fidelity in the absence of those truncation errors, the model would require retuning. A coarse-grained attribution of the time step sensitivities is carried out by shortening time steps used in various components of EAM or by revising the numerical coupling between some processes. Our analysis leads to the finding that the marked decreases in the subtropical low-cloud fraction and total cloud radiative effect are caused not by the step size used for the collectively subcycled turbulence, shallow convection, and stratiform cloud macrophysics and microphysics parameterizations but rather by the step sizes used outside those subcycles. Further analysis suggests that the coupling frequency between the subcycles and the rest of EAM significantly affects the subtropical marine stratocumulus decks, while deep convection has significant impacts on trade cumulus. The step size of the cloud macrophysics and microphysics subcycle itself appears to have a primary impact on cloud fraction in the upper troposphere and also in the midlatitude near-surface layers. Impacts of step sizes used by the dynamical core and the radiation parameterization appear to be relatively small. These results provide useful clues for future studies aiming at understanding and addressing the root causes of sensitivities to time step sizes and process coupling frequencies in EAM. While this study focuses on EAMv1 and the conclusions are likely model-specific, the presented experimentation strategy has general value for weather and climate model development, as the methodology can help researchers identify and understand sources of time integration error in sophisticated multi-component models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Local carbon reserves are insufficient for phloem terpene induction during drought in Pinus edulis in response to bark beetle–associated fungi

Stomatal closure during drought inhibits carbon uptake and may reduce a tree's defensive capacity. Limited carbon availability during drought may increase a tree's mortality risk, particularly if drought constrains trees' capacity to rapidly produce defenses during biotic attack. We parameterized a new model of conifer defense using physiological data on carbon reserves and chemical defenses before and after a simulated bark beetle attack in mature Pinus edulis under experimental drought. Attack was simulated using inoculations with a consistent bluestain fungus (Ophiostoma sp.) of Ips confusus, the main bark beetle colonizing this tree, to induce a defensive response. Trees with more carbon reserves produced more defenses but measured phloem carbon reserves only accounted for c. 23% of the induced defensive response. Our model predicted universal mortality if local reserves alone supported defense production, suggesting substantial remobilization and transport of stored resin or carbon reserves to the inoculation site. Our results show that de novo terpene synthesis represents only a fraction of the total measured phloem terpenes in P. edulis following fungal inoculation. Without direct attribution of phloem terpene concentrations to available carbon, many studies may be overestimating the scale and importance of de novo terpene synthesis in a tree's induced defense response.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

SNICAR-ADv4: a physically based radiative transfer model to represent the spectral albedo of glacier ice

Abstract. Accurate modeling of cryospheric surface albedo is essential for our understanding of climate change as snow and ice surfaces regulate the global radiative budget and sea-level through their albedo and mass balance. Although significant progress has been made using physical principles to represent the dynamic albedo of snow, models of glacier ice albedo tend to be heavily parameterized and not explicitly connected with physical properties that govern albedo, such as the number and size of air bubbles, specific surface area (SSA), presence of abiotic and biotic light absorbing constituents (LACs), and characteristics of any overlying snow. Here, we introduce SNICAR-ADv4, an extension of the multi-layer two-stream delta-Eddington radiative transfer model with the adding–doubling solver that has been previously applied to represent snow and sea-ice spectral albedo. SNICAR-ADv4 treats spectrally resolved Fresnel reflectance and transmittance between overlying snow and higher-density glacier ice, scattering by air bubbles of varying sizes, and numerous types of LACs. SNICAR-ADv4 simulates a wide range of clean snow and ice broadband albedo (BBA), ranging from 0.88 for (30 µm) fine-grain snow to 0.03 for bare and bubble-free ice under direct light. Our results indicate that representing ice with a density of 650 kg m−3 as snow with no refractive Fresnel layer, as done previously, generally overestimates the BBA by an average of 0.058. However, because most naturally occurring ice surfaces are roughened “white ice”, we recommend modeling a thin snow layer over bare ice simulations. We find optimal agreement with measurements by representing cryospheric media with densities less than 650 kg m−3 as snow and larger-density media as bubbly ice with a Fresnel layer. SNICAR-ADv4 also simulates the non-linear albedo impacts from LACs with changing ice SSA, with peak impact per unit mass of LACs near SSAs of 0.1–0.01 m2 kg−1. For bare, bubble-free ice, LACs actually increase the albedo. SNICAR-ADv4 represents smooth transitions between snow, firn, and ice surfaces and accurately reproduces measured spectral albedos of a variety of glacier surfaces. This work paves the way for adapting SNICAR-ADv4 to be used in land ice model components of Earth system models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Modeling Interfaces to Support Low-Level Waste Disposal System Performance Assessments - 20366

In low-level waste (LLW) disposal sites, interfaces between cementitious materials, used as waste forms and/or engineered barriers, and the surrounding soil or backfill material are often encountered. Reactions across these interfaces may lead to chemical and structural alteration of the cementitious and backfill materials that may ultimately affect long-term performance. For example, the ingress of carbon dioxide from soil gas into a waste tank concrete shell at the Hanford Site or Savannah River Site (SRS) or a low-level waste disposal vault at SRS leads to carbonation of cement hydration products (e.g. portlandite, calcium-silicate hydrate or C-S-H, and ettringite). The result of carbonation is a decrease in pH in the cement paste portion of the tank wall or vault concrete structure that could ultimately lead to de-passivation of the embedded structural steel and, eventually, to cracking. Cracking can lead to increased ingress of water into the structure, corresponding increased release of constituents of concern, and to subsequent increased transport of these constituents to the surrounding environment. The mobility of trace constituents (e.g., radionuclides of concern) in waste forms may increase in response to changes in pH, pore structure, and mineralogical gradients. Depletion of portlandite results in subsequent decalcification of C-S-H that leads to changes in the cement strength and may lead to structural failure. At the boundary between waste forms and concrete vaults, the migration of sulfate ions from the salt waste form (e.g., the saltstone waste form used for LLW disposal at the Savannah River Site in Saltstone Disposal Units, or SDUs) into the barrier or vault concrete has been predicted to result in expansive mineral phase formation (or 'sulfate attack'). Expansive mineral formation could result in cracking and potential loss of structural integrity in the vault concrete. In the performance assessment (PA) for the SRS SDUs, prediction of both the carbonation and sulfate attack fronts are critical to assessing the long-term performance of these LLW disposal vaults. A general purpose geochemical reactive transport model has been developed in LeachXS/ORCHESTRA to evaluate typical interfaces for LLW disposal environments, including SDUs at SRS. The model may be used to predict the long-term performance of interfaces between cementitious materials and saltstone or backfill materials with respect to primary phases by considering mass transport, cement chemistry, geochemical speciation, and multi-ionic diffusion over interfaces between different materials. The present study compares the performance and interactions between vault concrete and saltstone waste form for a representative SDU scenario. The materials were characterized using pH-dependence leaching test (EPA 1313) and semi-dynamic transport tests (EPA 1315) with test results used to parameterize the interface model. The development of a credible inter model provides an important basis for supporting future PAs where carbonation and sulfate attack are important aging and degradation mechanisms. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2021

Characterization of oceanic post-cold frontal clouds and their model representation. Final report

Low-level clouds are ubiquitous over the oceans, and their impact on the Earth’ radiative budget is important. However, general circulation models (GCMs) that are used to represent our climate and its evolution still experience problems to represent their correct amount and radiative impact. In the extratropical latitudes (~30-60 N/S), cloud properties are strongly modulated by the occurrence of anti-cyclones (high pressure systems) and cyclones (low pressure systems). While deep clouds dominate in cyclones near the center and along the cold and warm fronts, low-level clouds populate the colder regions of the cyclones, and more specifically the region in the wake of the cold fronts, the post-cold frontal region (PCF). The GCMs mentioned above were found to underestimate the amount of these PCF clouds, and this caused errors in their representation of the amount of solar radiation at the surface of the southern hemisphere oceans, causing an excess in absorption and errors in long-term predictions. To better understand which processes are ill-represented in the GCMs to cause such issues, we used ground-based observations from the Eastern North Atlantic (ENA) ARM site and a high resolution regional model (called WRF) to examine the properties of PCF clouds. With the WRF model, we simulated the passage of a cold front at the ENA site using theoretically distinct representations of convection and the physics of the boundary layer. While these physical representations had little impact on the timing, structure and circulation of the cold front, we found a high sensitivity to how the large scale information was fed into the model. Errors rapidly developed if the incoming flow information was injected too far from the site, but 1000 km was found to be an optimal distance between the outer boundary and the ENA site for a realistic cold front passage. With this model, we then conducted a series of experiment where convection and boundary layer schemes were changed, and revealed that 1) the clouds were more sensitive to the convection than boundary layer schemes and 2) the interaction between the two has a large impact. Furthermore, changing the representation of the physics in the model had greater impacts on PCF clouds than changing the initial conditions. As large-scale models increase their spatial resolution, convection will eventually be resolved but it remains that their representation of the physics of the boundary layer will matter for the correct representation of shallow (low-level) clouds. Observations were used to find which large-scale information had a clear relationship with PCF cloud and precipitation properties, in view of helping with GCM parameterizations and evaluation. Both lower troposphere stability and surface forcing were found to be important for PCF clouds that experience more vigorous dynamics than their tropical or subtropical counterparts. For cloud base and top heights, a simple measure of the change in potential temperature between the surface and the 800 hPa was found to be a good predictor. This relationship was found to hold as well when using data from campaigns recently conducted in the southern ocean south of Tasmania. A recent GCM was found to succeed in representing similar relations. For cloud depth, water content, and different characteristics of precipitation both the temperature contrast between the surface and near-air and the near-surface wind speed had a strong relationship. To summarize, through observational constraints and regional model experiments, the project has provided guidance to evaluate GCMs and tools to improve their representation of low-level clouds in the extratropics. It has also highlighted differences in behavior between these clouds and those typically found in the subtropics, thus opening new avenues of research for shallow clouds in general.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Development and assessment of hierarchical multi-reward reinforcement learning based potential for silicene with state-of-the-art models

We develop a new interatomic force field for Silicene, a 2D material with a buckled hexagonal lattice structure with high polymorphism. We introduce new parameterizations of a Tersoff model using a hierarchical multi-reward reinforcement learning (RL) methodology coupled with a continuous Monte Carlo Tree Search optimization. Our model significantly outperforms existing methods by enhancing the accuracy of predictions for the structural and thermodynamic properties of seven silicene polymorphs-including structure, energy, equation of state, elasticity, and phonon dispersion-when compared to established models. We further make a comprehensive comparison of the various models in predicting the mechanical and thermal properties of silicene. We trace the origin of the improved performance to the description of the angular dependence in the bond-order term, suggesting that modifying the angular terms in short-range models is essential to capture the structural diversity in low dimensional systems.

2D materials↗
DOE OSTI2023

Terrain-Influenced Winds and Fire-Fire Interactions in Wildland Fire Simulations [Dissertation]

Ensemble-based approaches to prescribed fire planning cannot be supported by computational fluid dynamics based models like FIRETEC and the Wildland-Urban Interface Fire Dynamics Simulator (WFDS) because they are too computationally expensive and cannot leverage large eddy simulation approaches like CAWFE and WRF-SFIRE because they have too coarse of resolution. QUIC-Fire was developed to fill this gap but it cannot currently address complex terrain, that is typical for instance in the Western United States. This dissertation describes a variety of improvements made to QUIC-Fire and its various incorporated algorithms in an effort to make it a viable tool in simulating wildland and prescribed fires on terrain. The modifications made to QUIC-Fire are described in three chapters. The first chapter describes the extension of the diagnostic wind model QUIC-URB, the wind engine of QUIC-Fire, to a terrain-following coordinate system. The terraininfluenced winds it generates are analyzed and compared. In particular, this chapter presents the mathematical derivation of the wind solver leading to a linear system of equations that are solved through the successive over-relaxation method. The model is validated against a standard test used in previous works (the Askervein Hill) and against a new dataset from measurements in the Socorro Mountains, New Mexico. The terrain-following implementation captures the correct phenomenology for the isolated Askervein Hill, with a wind speed up at the top of the hill. The model agrees well with measurements on the upwind side of the peak, but overestimates speed-up on the downwind side of the hill. This is due to the inability of the model to generate flow separation and wake-eddy dynamics. On a common laptop, the divergence-free wind field is obtained in 6 s, making the solver appealing for coupled fire-atmosphere simulations. The Socorro Mountain is highly complex, with many cliff faces, peaks, and valleys. Although the model captures the magnitude and direction of inlet and outlet areas of the domain, it performs rather poorly in the valley region and in the regions near the steep cliffs. Hence, the model shows good agreement with data in areas of open sloped terrain but lacks in areas where flow separation and thermally driven effects may be present (neither effect is addressed in this work). In the second chapter the implementation of the terrain-following version of QUIC-URB into QUIC-Fire, and the necessary changes needed to include terrain are described. No changes to the underlying fire spread algorithm are made other than what is required to correctly account for the inclusion of terrain. Previously published FIRETEC results that use five different topographies that share the same centerline profile are compared to simulation results from the modified QUIC-Fire that use the same topographies and fuels. QUIC-Fire results show overall similar behaviors in terms of how the topographies affect fire shapes and trends in spread rates. Due to the terrain-following version of QUIC-URB being unable to generate flow separations at the crest of hills, fire spread rates in these regions across all topographies are over-predicted when compared to FIRETEC. Lateral fire growth shows similar trends with FIRETEC between topographies but does not capture the increase in spread due to a diagonal interface between grassland and forested fuel region of the domain. These results suggest that there are three algorithms within QUIC-Fire that could use improvement: how flame tilt angle is accounted for, the incorporation of non-local drag effects, and the inclusion of the wake-eddy parameterizations that are used in QUIC-URB. Lastly, the third chapter describes a modification to the initial guess used for the calculation of the QUIC-URB mass-conserved wind solution during fire simulations. The modification is aimed at improving fire-fire interactions in QUIC-Fire simulations. The modification consists of using the solution from the previous timestep as the starting point for the calculation of the solution for the next timestep. Fire-fire interactions is greatly improved by the change but a new source of error is introduced. Due to how plumes are modelled in QUIC-Fire the new solution contains errors where gaps in the plume structure are present. However, these errors are mostly limited to the upper atmosphere, where they do not affect fire behavior at the surface, and their magnitude isn’t significant enough to discount the amount of new fire phenomenology now captured in QUIC-Fire with the change.

58 GEOSCIENCES↗
DOE OSTI2020

Predicting secondary organic aerosol phase state and viscosity and its effect on multiphase chemistry in a regional-scale air quality model

Atmospheric aerosols are a significant public health hazard and have substantial impacts on the climate. Secondary organic aerosols (SOAs) have been shown to phase separate into a highly viscous organic outer layer surrounding an aqueous core. This phase separation can decrease the partitioning of semi-volatile and low-volatile species to the organic phase and alter the extent of acid-catalyzed reactions in the aqueous core. A new algorithm that can determine SOA phase separation based on their glass transition temperature (T g ), oxygen to carbon (O:C) ratio and organic mass to sulfate ratio, and meteorological conditions was implemented into the Community Multiscale Air Quality Modeling (CMAQ) system version 5.2.1 and was used to simulate the conditions in the continental United States for the summer of 2013. SOA formed at the ground/surface level was predicted to be phase separated with core–shell morphology, i.e., aqueous inorganic core surrounded by organic coating 65.4 % of the time during the 2013 Southern Oxidant and Aerosol Study (SOAS) on average in the isoprene-rich southeastern United States. Our estimate is in proximity to the previously reported ~70 % in literature. The phase states of organic coatings switched between semi-solid and liquid states, depending on the environmental conditions. The semi-solid shell occurring with lower aerosol liquid water content (western United States and at higher altitudes) has a viscosity that was predicted to be 10 2 –10 12 Pa s, which resulted in organic mass being decreased due to diffusion limitation. Organic aerosol was primarily liquid where aerosol liquid water was dominant (eastern United States and at the surface), with a viscosity <10 2 Pa s. Phase separation while in a liquid phase state, i.e., liquid–liquid phase separation (LLPS), also reduces reactive uptake rates relative to homogeneous internally mixed liquid morphology but was lower than aerosols with a thick viscous organic shell. The sensitivity cases performed with different phase-separation parameterization and dissolution rate of isoprene epoxydiol (IEPOX) into the particle phase in CMAQ can have varying impact on fine particulate matter (PM 2.5 ) organic mass, in terms of bias and error compared to field data collected during the 2013 SOAS. This highlights the need to better constrain the parameters that govern phase state and morphology of SOA, as well as expand mechanistic representation of multiphase chemistry for non-IEPOX SOA formation in models aided by novel experimental insights.

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

The GFDL‐CM4X Climate Model Hierarchy, Part II: Case Studies

This paper is Part II of a two‐part paper that documents the Climate Model version 4X (CM4X) hierarchy of coupled climate models developed at the Geophysical Fluid Dynamics Laboratory. Part I of this paper is presented in Griffies et al. (2025a, https://doi.org/10.1029/2024MS004861 ). Here we present a suite of case studies that examine ocean and sea ice features that are targeted for further research, which include sea level, eastern boundary upwelling, Arctic and Southern Ocean sea ice, Southern Ocean circulation, and North Atlantic circulation. The case studies are based on experiments that follow the protocol of version 6 from the Coupled Model Intercomparison Project. The analysis reveals a systematic improvement in the simulation fidelity of CM4X relative to its CM4.0 predecessor, as well as an improvement when refining the ocean/sea ice horizontal grid spacing from the 0.25° of CM4X‐p25 to the 0.125° of CM4X‐p125. Even so, there remain many outstanding biases, thus pointing to the need for further grid refinements, enhancements to numerical methods, and/or advances in parameterizations, each of which target long‐standing model biases and limitations.

54 ENVIRONMENTAL SCIENCES↗