Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Variational analysis”

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 595 records · Page 33
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Tunable angle-dependent electrochemistry at twisted bilayer graphene with moiré flat bands

Tailoring electron transfer dynamics across solid–liquid interfaces is fundamental to the interconversion of electrical and chemical energy. Stacking atomically thin layers with a small azimuthal misorientation to produce moiré superlattices enables the controlled engineering of electronic band structures and the formation of extremely flat electronic bands. Here, we report a strong twist-angle dependence of heterogeneous charge transfer kinetics at twisted bilayer graphene electrodes with the greatest enhancement observed near the ‘magic angle’ (~1.1°). Furthermore, this effect is driven by the angle-dependent tuning of moiré-derived flat bands that modulate electron transfer processes with the solution-phase redox couple. Combined experimental and computational analysis reveals that the variation in electrochemical activity with moiré angle is controlled by a structural relaxation of the moiré superlattice at twist angles of <2°, and ‘topological defect’ AA stacking regions, where flat bands are localized, produce a large anomalous local electrochemical enhancement that cannot be accounted for by the elevated local density of states alone.

36 MATERIALS SCIENCE↗
DOE OSTI2022

AtomAI framework for deep learning analysis of image and spectroscopy data in electron and scanning probe microscopy

Over the past several decades, electron and scanning probe microscopes have become critical components of condensed matter physics, materials science and chemistry research. At the same time, the infrastructure for establishing a connection between microscopy observations and materials behaviour over a broader parameter space is lacking. In this work, we introduce AtomAI, an open-source software package bridging instrument-specific Python libraries, deep learning and simulation tools into a single ecosystem. AtomAI allows direct applications of deep neural networks for atomic and mesoscopic image segmentation converting image and spectroscopy data into class-based local descriptors for downstream tasks such as statistical and graph analysis. For atomically resolved imaging data, the output is types and positions of atomic species, with an option for subsequent refinement. AtomAI further allows the implementation of a broad range of image and spectrum analysis functions, including invariant variational autoencoders for disentangling structural factors of variation and im2spec type of encoder–decoder models for mapping structure–property relationships. Finally, our framework allows seamless connection to the first principles modelling with a Python interface on the inferred atomic positions.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Applying equity principles leads to higher carbon removal obligations in Canada

Despite net-zero pledges, consensus on national responsibilities for carbon dioxide removal (CDR) strategies is lacking. Here, we use integrated assessment modeling to examine equity-informed estimates of Canada’s remaining carbon budgets, exploring CDR’s role at net-zero and beyond. Gigaton-scale CDR efforts post-2050 are needed to address Canada’s carbon debt under various burden-sharing principles. Cumulative negative emissions (2050-2100) could increase from 7.5 GtCO 2 in the Net-Zero scenario to 20.3 GtCO 2 in equity-informed scenarios. By 2100, a CDR portfolio, including bioenergy with carbon capture and storage, direct air capture, and enhanced weathering could contribute up to ~500 MtCO 2 /year of removals. The projected average CDR growth rates, 2.8%-16%/year, align with the historical adoption rates of ammonia synthesis and biomass consumption in Canada, underscoring the importance of drawing lessons from past successes. Socio-economic and technological sensitivity analysis highlights that, despite variations in the role of individual CDR technologies, CDR remains essential for Canada’s post-net-zero commitments.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Unsupervised machine learning for unbiased chemical classification in X-ray absorption spectroscopy and X-ray emission spectroscopy

Here we report a comprehensive computational study of unsupervised machine learning for extraction of chemically relevant information in X-ray absorption near edge structure (XANES) and in valence-to-core X-ray emission spectra (VtC-XES) for classification of a broad ensemble of sulphorganic molecules. By progressively decreasing the constraining assumptions of the unsupervised machine learning algorithm, moving from principal component analysis (PCA) to a variational autoencoder (VAE) to t-distributed stochastic neighbour embedding (t-SNE), we find improved sensitivity to steadily more refined chemical information. Surprisingly, when embedding the ensemble of spectra in merely two dimensions, t-SNE distinguishes not just oxidation state and general sulphur bonding environment but also the aromaticity of the bonding radical group with 87% accuracy as well as identifying even finer details in electronic structure within aromatic or aliphatic sub-classes. We find that the chemical information in XANES and VtC-XES is very similar in character and content, although they unexpectedly have different sensitivity within a given molecular class. We also discuss likely benefits from further effort with unsupervised machine learning and from the interplay between supervised and unsupervised machine learning for X-ray spectroscopies. Our overall results, i.e., the ability to reliably classify without user bias and to discover unexpected chemical signatures for XANES and VtC-XES, likely generalize to other systems as well as to other one-dimensional chemical spectroscopies.

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

Bolometer tomography on Wendelstein 7-X for study of radiation asymmetry

The algorithm for bolometer tomography at Wendelstein 7-X (W7-X) has been recently improved using a novel regularization functional, based on relative gradient smoothing (RGS) of the sought radiation profile. It has been validated using radiation patterns provided by 3D modeling under real plasma conditions as phantoms and then applied to bolometer measurements performed during the first divertor operation phase of W7-X. The following results are presented: 1) edge-localized 2D radiation patterns with clearly resolved magnetic island radiation structures, 2) an up-down asymmetry in the impurity radiation that is not captured by the 3D edge plasma transport modeling, 3) reversal of the asymmetry with reversed magnetic field direction. Further analysis reveals a poloidal variation of the emissivity in the outer confined plasma region with a field-direction dependent asymmetry, also supported by the soft X-ray measurements. This asymmetry is considered to be related to asymmetric impurity distributions, driven by the pronounced ion-impurity friction force at the plasma edge where the collisionality of the W7-X plasma is sufficiently high to develop impurity asymmetry as predicted by neoclassical theory of parallel impurity transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Improvements to the Blade Element Momentum Formulation of OpenFAST for Skewed Inflows

In this work, we modify the blade element momentum algorithm of OpenFAST to improve its predictions under large skewed inflow conditions. We use the well-known Glauert's skew correction and introduce continuous extension of the model for high-thrust conditions. We present the rationale behind Glauert's empirical model and discuss the different conventions possible for the axial induction factor. We verify the model against the higher-fidelity lifting-line vortex method and blade-resolved computational fluid dynamics, and we observe that the new implementation enhances the accuracy and reliability of OpenFAST's aerodynamic modeling capabilities in conditions involving large skew angles. For the parametric studies run using the different codes, we find that the power changes with the skew angle as cos 1.7 (θ skew ) and the thrust as cos 0.65 (θ skew ). An analysis of the azimuthal variation of the induced velocities in the rotor plane reveals that current redistribution models used in blade element momentum codes may need to be refined.

17 WIND ENERGY↗
DOE OSTI2023

Combining field and remote sensing data to estimate forest canopy damage and recovery following tropical cyclones across tropical regions

Abstract As tropical forests cycle the most water and carbon, it is crucial to understand the short- and long-term effects of intensifying cyclones on these ecosystems. Soil nutrient status has been shown to moderate forest cyclone responses using field litterfall measurements, but litterfall is one of the multiple cyclone impact metrtics, which may or may not be correlated with one another or with site nutrients. We used remotely sensed vegetation indices to quantify immediate damage and two-year recovery for 42 cases across nine tropical forests in Hawaii, Puerto Rico, Mexico, Australia, and Taiwan affected by 12 cyclones between 2004 and 2017. We tested whether changes in leaf area index (LAI) and enhanced vegetation index (EVI) correlated with changes in litterfall observations and how changes varied with total soil phosphorus (P) concentrations across regions. We compared cyclone-induced changes and recovery of LAI and EVI to litterfall observations compiled in a pantropical meta-analysis. We found large variation in changes in LAI and EVI across forests, with the greatest reductions in LAI (−77%) and EVI (−77%) in Mexico (Jalisco) and Puerto Rico, respectively. LAI ( r = −0.52) and EVI ( r = −0.60) changes correlated with those in litterfall across cases. Post-cyclone data showed recovery of LAI by four months, EVI by two months, and litterfall by ten months. We detected larger changes in LAI and EVI in forests with higher soil P, but these relationships were not significant when accounting for cyclone and site as random effects. Principal component analyses indicated a regional clustering of cases related to their contrasting cyclone regimes, with the frequency and intensity of cyclone events negatively correlated. Overall, remote sensing observations complement but do not substitute for ground observations that reveal cyclone damage and post-cyclone recovery in tropical forests, and soil phosphorus moderates some but not all metrics of stability in response to cyclones.

Bloom, Dellena E. (ORCID:0000000205981747)↗
DOE OSTI2026

Glauber-theory analysis of nuclear reactions on a 12 C target with variational Monte Carlo wave functions

The application of Glauber theory has been playing an increasingly important role with the study of unstable or exotic nuclei. Its adaptation to medium and high-energy nucleus-nucleus collisions is severely limited because one has to evaluate the matrix elements of multiple-scattering operators. The extraction of physical observables has been done using ‘approximate’ Glauber theory whose validity is hard to evaluate. Here, we perform a full calculation of the matrix elements using Monte Carlo integration and analyze the elastic differential cross sections and the total reaction cross sections for p+¹²C, ⁴,⁶He+¹²C, and ¹²C+¹²C collisions. We use the variational Monte Carlo wave functions for ⁴,⁶He and ¹²C obtained by using realistic two- and three-nucleon potentials. We demonstrate the performance of the Glauber-theory calculations by comparing with available experimental data. We further discuss the accuracy of the conventional approximate methods in the light of the cumulant expansion for Glauber’s phase-shift function.

Horiuchi, W. [Osaka Metropolitan University (Japan↗
DOE OSTI2022

Detailed analysis of excited-state systematics in a lattice QCD calculation of 𝑔 𝐴

Excited state contamination remains one of the most challenging sources of systematic uncertainty to control in lattice QCD calculations of nucleon matrix elements and form factors: early time separations are contaminated by excited states and late times suffer from an exponentially bad signal-to-noise problem. High-statistics calculations at large time separations ≳ 1 fm are commonly used to combat these issues. In this work, focusing on g A , we explore the alternative strategy of utilizing a large number of relatively low-statistics calculations at short to medium time separations (0.2–1 fm), combined with a multistate analysis. On an ensemble with a pion mass of approximately 310 MeV and a lattice spacing of approximately 0.09 fm, we find this provides a more robust and economical method of quantifying and controlling the excited state systematic uncertainty. A quantitative separation of various types of excited states enables the identification of the transition matrix elements as the dominant contamination. The excited state contamination of the Feynman-Hellmann correlation function is found to reduce to the 1% level at approximately 1 fm while, for the more standard three-point functions, this does not occur until after 2 fm. Critical to our findings is the use of a global minimization, rather than fixing the spectrum from the two-point functions and using them as input to the three-point analysis. We find that the ground state parameters determined in such a global analysis are stable against variations in the excited state model, the number of excited states, and the truncation of early-time or late-time numerical data.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Evaluation of Converter Performance Considering Static and Dynamic Device Part-to-Part Variability

This paper presents a methodology to incorporate and analyze the impact of semiconductor device part-to-part variation on power converter performance. By integrating extensive static and dynamic device characterization data with an automated compact model generation process that reflects manufacturing variability, device models with inherent variability features are utilized in converter simulations for a comprehensive assessment of performance impacts. The traditional converter performance evaluation process typically yields fixed efficiency values, often dismissing the inherent part-to-part variability caused by the manufacturing process of semiconductor devices. To address this limitation, a large population of devices was characterized to capture variations in static parameters-such as transfer, output, and capacitance characteristics-as well as dynamic behaviors, including switching losses. This data-driven approach enables the development of individual compact models, which were then integrated into converter simulations to evaluate efficiency ranges rather than single point estimated values. The converter simulation results show that part-to-part component variation can lead to significant efficiency deviations, exceeding several percentage points in high-power conversion applications. By offering a more accurate representation of converter behavior under real-world manufacturing conditions, this methodology enables designers to anticipate performance variability, improving the robustness of power converter designs.

device characterization↗
DOE OSTI2023

Assisted migration in a warmer and drier climate: less climate buffering capacity, less facilitation and more fires at temperate latitudes?

Assisted tree migration has been proposed as a conceptual solution to mitigate lags in biotic responses to anthropogenic climate change. The rationale behind this concept is that tree species currently growing under warmer and drier climates will be more resistant and resilient to the new climatic conditions than tree species naturally growing in currently wetter and colder climates. However, we hypothesize that, by being more stress‐tolerant to warmer and drier conditions, translocated species should exhibit different functional attributes, which could induce important ecological and societal costs and overcome the desired benefits of maintaining wood production and other ecosystem services. We used principal component analysis (PCA) to analyze variation in seven traits of 106 tree and tall shrub species from contrasting latitudinal distributions in western North America and Europe to predict the potential functional changes of forest ecosystems due to the translocation of tree species from low to high latitudes. We show that species from both continents differed primarily by their position on the leaf economy spectrum (LES) and their size traits. Even though, in Europe, differences in LES were significantly correlated to species southern latitudinal positions, in both continents differences in size traits were significantly correlated to latitude. These results suggest that assisted migration by translocating more conservative species of shorter stature in currently cooler climates should decrease the buffering capacity of forest canopies, decrease facilitation for understory species, and increase wildfire risks, whose effects have the potential to accelerate climate warming through negative atmospheric feedback processes. As an alternative solution to assisted migration that may accelerate rather than mitigate climate change, we recommend that foresters gradually diversify the vertical structure and layering of the existing forest canopy to maintain a sustainable water cycle and energy balance between the soil, the tree and the atmosphere without increasing the wildfire risk.

Environmental Sciences & Ecology↗
DOE OSTI2022

Estimating SHmax azimuth with P sources and vertical geophones: Use P-P reflection amplitudes or use SV-P reflection times?

We compared two methods for extracting the azimuth of maximum horizontal stress (SHmax) from 3D land-based seismic data generated by a P source and recorded with vertical geophones. In the first method, we used the direct-SV mode that is produced by all land-based P sources. P sources generate SV illumination that radiates in all azimuth directions from a source station and creates SV-P reflections that are recorded by vertical geophones. Unless stratigraphy has steep dip, SV-P raypaths recorded by vertical geophones are the reverse of P-SV raypaths recorded by horizontal geophones. Thus, SV-P data provide the same S-wave sensitivity to stress fields as popular P-SV data do. In the second method, we retrieved P-P reflections and then performed an amplitude-variation-with-azimuth (AVA) analysis of the amplitude-gradient behavior of P-P reflection wavelets. We did this analysis in narrow azimuth corridors to determine the gradient of reflection-wavelet amplitudes as a function of azimuth. This P-P AVA amplitude-gradient method has been of great interest in the reflection seismology community since it was introduced in the late 1990s. Each of these methods, AVA analysis of the gradient of P-P reflection amplitudes and azimuth-dependent arrival times of SV-P reflections, can be used to determine the azimuth of SHmax stress. We compare the results of the two methods with ground truth measurements of SHmax azimuth at a CO 2 sequestration site in the Michigan Basin. SHmax azimuths were determined from P-P and SV-P data at three major boundaries at depths of approximately 3500 ft (1067 m), 5500 ft (1676 m), and 7500 ft (2286 m). Two estimates of SHmax azimuth (one using SV-P data and one using P-P data) were made at each stacking bin inside a 24 mi 2 (62 km 2 ) image space. The result was approximately 98,000 estimates of SHmax azimuth across each of these three boundaries for each of these two prediction strategies. Histogram displays of PP AVA gradient estimates had peaks at correct azimuths of SHmax at all three depths, but the spread of the distributions widened with depth and split into two peaks at the deepest boundary. In contrast, each histogram of SHmax azimuth predicted by azimuth-dependent SV-P traveltimes had a single, definitive peak that was positioned at the correct SHmax azimuth at all three boundary depths.

Geochemistry & Geophysics↗
DOE OSTI2022

Classification and Localization of Fracture-Hit Events in Low-Frequency Distributed Acoustic Sensing Strain Rate with Convolutional Neural Networks

Summary Distributed acoustic sensing (DAS) has been used in the oil and gas industry as an advanced technology for surveillance and diagnostics. Operators use DAS to monitor hydraulic fracturing activities, examine well stimulation efficacy, and estimate complex fracture system geometries. Particularly, low-frequency DAS can detect geomechanical events such as fracture hits because hydraulic fractures propagate and create strain rate variations in the rock. Analysis of DAS data today is mostly done post-job and subject to interpretation methods. However, the continuous and dense data stream generated live by DAS poses the opportunity for more efficient and accurate real-time data-driven analysis. The objective of this study is to develop a machine learning-based workflow that can identify and locate fracture-hit events in simulated strain rate responses correlated with low-frequency DAS data. In this paper, “fracture hit” refers to a hydraulic fracture originating from a stimulated well intersecting an offset well. We start with building a single fracture propagation model to produce strain rate patterns observed at a hypothetical monitoring well. This model is used to generate two sets of strain rate responses with one set containing fracture-hit events. The labeled synthetic data are then used to train a custom convolutional neural network (CNN) model for identifying the presence of fracture-hit events. The same model is trained again for locating the event with the output layer of the model replaced with linear units. We achieved near-perfect predictions for both event classification and localization. These promising results prove the feasibility of using CNN for real-time event detection from fiber-optic sensing data. Additionally, we use edge detection techniques to recognize fracture-hit event patterns in strain rate images. The fracture-hit location can be identified using recognized pixels in the image. The accuracy of edge detection-based location identification is also plausible, but edge detection is dependent on the assumption of pattern shape and image quality, hence it is less robust compared to CNN models. This comparison further supports the need for CNN applications in image-based real-time fiber-optic sensing event detection.

Engineering↗
DOE OSTI2019

North American Winter Dipole: Observed and Simulated Changes in Circulations

In recent years, a pair of large-scale circulation patterns consisting of an anomalous ridge over northwestern North America and trough over northeastern North America was found to accompany extreme winter weather events such as the 2013–2015 California drought and eastern U.S. cold outbreaks. Referred to as the North American winter dipole (NAWD), previous studies have found both a marked natural variability and a warming-induced amplification trend in the NAWD. In this study, we utilized multiple global reanalysis datasets and existing climate model simulations to examine the variability of the winter planetary wave patterns over North America and to better understand how it is likely to change in the future. We compared between pre- and post-1980 periods to identify changes to the circulation variations based on empirical analysis. It was found that the leading pattern of the winter planetary waves has changed, from the Pacific–North America (PNA) mode to a spatially shifted mode such as NAWD. Further, the potential influence of global warming on NAWD was examined using multiple climate model simulations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Harvesting Reactor Pressure Vessel Beltline Material from the Decommissioned Zion Nuclear Power Plant Unit 1

The decommissioning of the Zion Nuclear Power Plant (NPP) provided a unique opportunity to harvest and study service-aged reactor pressure vessel (RPV) beltline materials. This work, conducted through the U.S. Department of Energy’s Light Water Reactor Sustainability (LWRS) Program, aims to improve the understanding of radiation-induced embrittlement to support extended nuclear plant operations. Material segments containing the Linde 80 flux, wire heat 72105 (WF-70) beltline weld and the A533B Heat B7835-1 base metal, obtained from the intermediate shell region with a peak fluence of 0.7 × 10 19 n/cm 2 (E > 1.0 MeV), were extracted, cut into blocks, and machined into test specimens for mechanical and microstructural characterization. The segmentation process involved oxy-propane torch-cutting, followed by precision machining using wire saws and electrical discharge machining (EDM). A chemical composition analysis confirmed the expected variations in alloying elements, with copper levels being notably higher in the weld metal. The harvested specimens enable a detailed evaluation of through-wall embrittlement gradients, a comparison with the existing surveillance data, and the validation of predictive embrittlement models. This study provides critical data for assessing long-term reactor vessel integrity, informing aging-management strategies, and supporting regulatory decisions to extend the life of nuclear plants. This article is a revised and expanded version of a paper entitled, “Current Status of the Characterization of RPV Materials Harvested from the Decommissioned Zion Unit 1 Nuclear Power Plant”, PVP2017-65090, which was accepted and presented at the ASME 2017 Pressure Vessels and Piping Conference, Waikoloa, HI, USA, 16–20 July 2017.

harvesting beltline material↗
DOE OSTI2024

Fire–precipitation interactions amplify the quasi-biennial variability in fires over southern Mexico and Central America

Fires have great ecological, social, and economic impact. However, fire prediction and management remain challenges due to a limited understanding of their roles in the Earth system. Fires over southern Mexico and Central America (SMCA) are a good example of this, greatly impacting local air quality and regional climate. Here we report that the spring peak (April–May) of fire activities in this region has a distinct quasi-biennial signal based on multiple satellite datasets measuring different fire characteristics. The variability is initially driven by quasi-biennial variations in precipitation. Composite analysis indicates that strong fire years correspond to suppressed ascending motion and weakened precipitation over the SMCA. The anomalous precipitation over the SMCA is further found to be mostly related to the East Pacific–North Pacific (EP-NP) pattern 2 months prior to the fire season. The positive phase of the EP-NP leads to enhanced precipitation over the eastern US but suppressed precipitation over the SMCA, similar to the spatial pattern of precipitation differences between strong and weak fire years. Meanwhile, the quasi-biennial signals in precipitation and fires appear to be amplified by their interactions through a positive feedback loop at short timescales. Model simulations show that in strong fire years, more aerosol particles are released and transported downstream over the Gulf of Mexico and the eastern US, where suspended light-absorbing aerosols warm the atmosphere and cause the ascending motion of the air aloft. Subsequently, a compensating downward motion is formed over the region of the fire source and ultimately suppresses precipitation and intensifies fires. Statistical analysis shows the different durations of the two-way interaction, where the fire suppression effect of precipitation lasts for more than 20 d, while fire leads to a decrease in precipitation at shorter timescales (3–5 d). This study demonstrates the importance of fire–climate interactions in shaping the fire activities on an interannual scale and highlights how precipitation–fire interactions at short timescales contribute to the interannual variability in both fire and precipitation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Prognostic simulations of mixed-phase clouds with model AC-1D v1.0: the impact of aerosol types and freezing parameterizations on ice crystal budgets

Mixed-phase clouds at high latitudes contribute to the uncertainty in predicting cloud feedbacks and climate sensitivity, mainly due to the complexity of microphysical processes that influence the partitioning between the supercooled liquid and ice phases, and hence, cloud radiative effects on regional scales. Particularly in Arctic mixed-phase clouds, the activation of ice-nucleating particles (INPs) from various aerosol populations remains a leading source of uncertainty. We developed an aerosol-cloud one-dimensional (AC-1D) model, which provides a novel framework to prognostically treat INP and ice crystal budgets while explicitly accounting for polydisperse and multicomponent aerosol that activate INPs following different freezing parameterizations. The AC-1D model is informed by large-eddy simulations to probe the impact of INP representation on predicted ice crystal number concentrations (N i ) and ice crystal budgets in mixed-phase Arctic stratus. We apply three immersion freezing (IMF) parameterizations, two time-independent (singular) and one time-dependent (classical nucleation theory), to predict the evolution of the INP reservoir and resulting ice crystal budget from polydisperse mineral dust, organic (humic-like substances), and sea spray aerosol particle size distributions. Our analysis focuses on how variations in aerosol number concentration and cloud system parameters such as cloud cooling rate, cloud-top entrainment rate, and ice crystal fall speed influence the INP reservoir and ice crystal budgets. Furthermore, this study investigates the competitive ice nucleation dynamics in mixed aerosol environments and provides a process-level quantification of the INP budget terms, which directly controls ice crystal budgets. For all studied case scenarios, the aerosol types and associated particle size distributions significantly impact INP and N i , and the choice between a time-dependent and a singular freezing description yields orders-of-magnitude differences in the predicted INP and N i over the 10 h simulation time, reflecting typical cloud lifetimes. Our results show that the influence of cloud cooling, INP entrainment, and sedimentation varies significantly depending on the chosen freezing parameterization. These findings underscore the critical need for robust IMF parameterizations and precise cloud system observations to enhance the accuracy of models in predicting mixed-phase cloud structure and evolution.

Arctic clouds↗