Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Theory of nucleation”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

A New Theory of Nucleate Pool Boiling in Arbitrary Gravity

Heat transfer rates specific to nucleate pool boiling under various conditions are determined by the dynamics of vapour bubbles that are originated and grow at nucleation sites of a superheated surface. A new dynamic theory of these bubbles has been recently developed on the basis of the thermodynamics of irreversible processes. In contrast to other existing models based on empirically postulated equations for bubble growth and motion, this theory does not contain unwarrantable assumptions, and both the equations are rigorously derived within the framework of a unified approach. The conclusions of the theory are drastically different from those of the conventional models. The bubbles are shown to detach themselves under combined action of buoyancy and a surface tension force that is proven to add to buoyancy in bubble detachment, but not the other way round as is commonly presumed. The theory ensures a sound understanding of a number of so far unexplained phenomena, such as effect caused by gravity level and surface tension on the bubble growth rate and dependence of the bubble characteristics at detachment on the liquid thermophysical parameters and relevant temperature differences. The theoretical predictions are shown to be in a satisfactory qualitative and quantitative agreement with observations. When being applied to heat transfer at nucleate pool boiling, this bubble dynamic theory offers an opportunity to considerably improve the main formulae that are generally used to correlate experimental findings and to design boiling heat removal in various industrial applications. Moreover, the theory makes possible to pose and study a great deal of new problems of essential impact in practice. Two such problems are considered in detail. One problem concerns the development of a principally novel physical model for the first crisis of boiling. This model allows for evaluating critical boiling heat fluxes under various conditions, and in particular at different gravity levels, with a good agreement with experimental evidence. The other problem bears upon equilibrium shapes of a detached bubble near a heated surface in exceedingly low gravity. In low gravity or in weightlessness, the bubble can remain in the close vicinity of the surface for a long time, and its shape is greatly affected by the Marangoni effect due to both temperature and possible surfactant concentration being nonuniform along the interface. The bubble performs at these conditions like a heat pipe, with evaporation at the bubble lower boundary and condensation at its upper boundary, and ultimately ensures a substantial increase in heat removal as compared with that in normal gravity. Some other problems relevant to nucleate pool and forced convection boiling heat transfer are also discussed.

Buyevich, Y. A.↗
NASA NTRS1995

Grain formation around carbon stars. 1: Stationary outflow models

Asymptotic giant branch (AGB) stars are known to be sites of dust formation and undergo significant mass loss. The outflow is believed to be driven by radiation pressure on grains and momentum coupling between the grains and gas. While the physics of shell dynamics and grain formation are closely coupled, most previous models of circumstellar shells have treated the problem separately. Studies of shell dynamics typically assume the existence of grains needed to drive the outflow, while most grain formation models assume a constant veolcity wind in which grains form. Furthermore, models of grain formation have relied primarily on classical nucleation theory instead of using a more realistic approach based on chemical kinetics. To model grain formation in carbon-rich AGB stars, we have coupled the kinetic equations governing small cluster growth to moment equations which determine the growth of large particles. Phenomenological models assuming stationary outflow are presented to demonstrate the differences between the classical nucleation approach and the kinetic equation method. It is found that classical nucleation theory predicts nucleation at a lower supersaturation ratio than is predicted by the kinetic equations, resulting in significant differences in grain properties. Coagulation of clusters larger than monomers is unimportant for grain formation in high mass-loss models but becomes more important to grain growth in low mass-loss situations. The properties of the dust grains are altered considerably if differential drift velocities are ignored in modeling grain formation. The effect of stellar temperature, stellar luminosity, and different outflow velocities are investigated. The models indicate that changing the stellar temperature while keeping the stellar luminosity constant has little effect on the physical parameters of the dust shell formed. Increasing the stellar luminosity while keeping the stellar temperature constant results in large differences in grain properties. For small outflow velocities, grains form at lower supersaturation ratios and close to the stellar photosphere, resulting in larger but fewer grains. The reverse is true when grains form under high outflow velocities, i.e., they form at higher supersaturation ratios, farther from the star, and are much smaller but at larger quantities.

Egan, Michael P.↗
DOE OSTI2022

Nucleation Rate Theory for Coordination Number: Elucidating Water-Mediated Formation of a Zigzag Na 2 SO 4 Morphology

Predicting and controlling nanostructure formation during nucleation can pave the way to synthesizing novel energy materials via crystallization. However, such control over nucleation and crystallization remains challenging due to an inadequate understanding of critical factors that govern evolving atomistic structures and dynamics. Here, in this paper, we utilize coordination number as a reaction coordinate and rate theory to investigate how sodium sulfate, commonly known as a phase-change energy material, nucleates in a supersaturated aqueous solution. In conjunction with ab initio and force field-based molecular dynamics simulation, the rate theoretical analysis reveals that sodium sulfate from an initially dissolved metastable state transits to a heterogeneous mixture of prenucleated clusters and finally to a large cylindrical zigzag morphology. Measurements of Raman spectra and their ab initio modeling confirm that this nucleated morphology contains a few waters for every sulfate. Rate processes such as solvent exchange and desolvation exhibit high sensitivity to the evolving prenucleation/nucleation structures, providing a means to distinguish between critical nucleation precursors. Desolvation and forming the first-shell interionic coordination structure via monomer-by-monomer addition around sulfates are found to explain the formation of large nuclei. Thus, a detailed understanding of the step-by-step structure formation across scales has been achieved. This can be leveraged to predict nucleation-related structures and dynamics and potentially control the synthesis of novel phase-change materials for energy applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2015

Thermodynamic Derivation of the Activation Energy for Ice Nucleation

Cirrus clouds play a key role in the radiative and hydrological balance of the upper troposphere. Their correct representation in atmospheric models requires an understanding of the microscopic processes leading to ice nucleation. A key parameter in the theoretical description of ice nucleation is the activation energy, which controls the flux of water molecules from the bulk of the liquid to the solid during the early stages of ice formation. In most studies it is estimated by direct association with the bulk properties of water, typically viscosity and self-diffusivity. As the environment in the ice-liquid interface may differ from that of the bulk, this approach may introduce bias in calculated nucleation rates. In this work a theoretical model is proposed to describe the transfer of water molecules across the ice-liquid interface. Within this framework the activation energy naturally emerges from the combination of the energy required to break hydrogen bonds in the liquid, i.e., the bulk diffusion process, and the work dissipated from the molecular rearrangement of water molecules within the ice-liquid interface. The new expression is introduced into a generalized form of classical nucleation theory. Even though no nucleation rate measurements are used to fit any of the parameters of the theory the predicted nucleation rate is in good agreement with experimental results, even at temperature as low as 190 K, where it tends to be underestimated by most models. It is shown that the activation energy has a strong dependency on temperature and a weak dependency on water activity. Such dependencies are masked by thermodynamic effects at temperatures typical of homogeneous freezing of cloud droplets; however, they may affect the formation of ice in haze aerosol particles. The new model provides an independent estimation of the activation energy and the homogeneous ice nucleation rate, and it may help to improve the interpretation of experimental results and the development of parameterizations for cloud formation.

Cirrus clouds↗
DOE OSTI2020

Kinetics and Mechanisms of ZnO to ZIF-8 Transformations in Supercritical CO2 Revealed by in situ X-ray Diffraction

ZIF-8 was synthesized in supercritical carbon dioxide (scCO2). In situ powder X-ray diffraction, ex situ microscopy and simulations, provides a encompassing view of the formation of ZIF-8 and intermediary ZnO@ZIF-8 composites in this nontraditional solvent. Time-resolved imaging exposed divergent physicochemical reaction pathways from previous studies of the growth of anisotropic ZIF-8 core@shell structures in traditional solvents. Synthetically relevant physio-chemical properties of scCO2 were integrated into classical nucleation theory, relating interfacial forces with 3D nu-cleation outcomes. The kinetics of crystallization were examined and displayed a characteristic signature of time- and temperature-dependent mechanisms over the extent of the reaction. Lastly, we show that subtle factors, such as the extent of reaction and the size/shape of sacrificial templates can tailor ZIF-8 composition and size, eliciting control over hierarchical porosity in a nonconventional green solvent.

Sinnwell, Michael A.↗
DOE OSTI2025

How are Heterogeneous Nucleation Rate Observations Influenced by Instrument Resolution?

Experimental measurements of the heterogeneous nucleation rate rely on counting the number of nuclei with time. However, the size of a thermodynamically stable nucleus is often a few nanometers in diameter and is below the resolution of most (in situ) measurement techniques that provide a statistically valid sample. Due to the finite resolution of the instruments and analysis methods, it is challenging to capture the incipient nuclei and the subsequent evolution of nuclei density over time. In this work, we demonstrate the impact of instrument resolution on observed nuclei densities by comparing numerical modeling with experimental results. Further, to achieve this, we implemented heterogeneous nucleation within the pore-scale reactive transport modeling framework using classical nucleation theory (CNT). We compared the modeling results with nucleation rates measured using X-ray nanotomography (XnT) and evaluated how these impact the apparent values of the prefactor and interfacial energy based on CNT and the crystal growth rate. Specifically, we applied a resolution threshold (artificial resolution limit) in the model during nuclei counting to resemble an experimental resolution, ranging from 15 to 500 nm. The findings reveal that the instrument resolution significantly impacts the apparent prefactor and interfacial energy. Both apparent prefactor and interfacial energy decrease with a decrease in the instrument resolution. While deviation in the prefactor due to resolution is anticipated, those in the interfacial energy are unexpected. The approach described here allows one to correct apparent nucleation rates that depend on the instrument’s resolution to derive “intrinsic” CNT parameters for the prefactor and interfacial energy.

47 OTHER INSTRUMENTATION↗
NASA NTRS1983

Homogeneous nucleation kinetics

Homogeneous nucleation kinetics are rederived in a manner fundamentally similar to the approach of classical nucleation theory with the following modifications and improvements. First, the cluster is a parent phase cluster and does not require energization to the parent state. Second, the thermodynamic potential used to describe phase stability is a continuous function along the pathway of phase decomposition. Third, the kinetics of clustering corresponds directly to the diffusional flux of monomers through the cluster distribution and are formally similar to classical theory with the resulting kinetic equation modified by two terms in the preexponential factor. These terms correct for the influence of a supersaturation dependent clustering within the parent phase and for the influence of an asymmetrical cluster concentration as a function of cluster size at the critical cluster size. Fourth, the supersaturation dependence of the nucleation rate is of the same form as that given by classical nucleation theory. This supersaturation dependence must however be interpreted in terms of a size dependent surface tension. Finally, there are two scaling laws which describe supersaturation to either constant nucleation rate or to the thermodynamically determined physical spinodal.

Rasmussen, D. H.↗
NASA NTRS1978

Grain formation in astronomical systems: A critical review of condensation processes

An analysis is presented of the assumption and the applicability of the three theoretical methods for calculating condensations in cosmic clouds where no pre-existing nuclei exist. The three procedures are: thermodynamic equilibrium calculations, nucleation theory, and a kinetic treatment which would take into account the characteristics of each individual collision. Thermodynamics provide detailed results on the composition temperature and composition of the condensate provided the system attains equilibrium. Because of the cosmic abundance mixture of elements, large supersaturations in some cases and low pressures, equilibrium is not expected in astronomical clouds. Nucleation theory, a combination of thermodynamics and kinetics, has the limitations of each scheme. Kinetics, not requiring equilibrium, avoids nearly all the thermodynamics difficulties but requires detailed knowledge of many reactions which thermodynamics avoids. It appears to be the only valid way to treat grain formation in space. A review of experimental studies is given.

Donn, B.↗
NASA NTRS1979

Nucleation of electrically charged droplets

The nucleating droplets or clusters in many nucleation environments (various colloidal plasmas in laboratory and technological applications, astrophysical condensation environments, etc.) are likely to be at a finite electric potential. This may be due either to the presence of electrons and ions in the gas phase or to the thermal ionization or photoionization of the droplets. The paper demonstrates that this potential may introduce a nontrivial modification in the conventional nucleation theory. Some results for the typical case of nucleation of water droplets are presented. The general conclusion is that the electric potential makes nucleation harder to achieve, thereby demonstrating the importance of a finite droplet potential in the theory of nucleation.

De, B. R.↗
DOE OSTI2020

The early-stage growth and reversibility of Li electrodeposition in Br-rich electrolytes

Significance Spatial variations in chemical composition and transport properties of the material phases (interphases) formed on reactive metals in liquid electrolytes are thought to be responsible for the propensity of metal battery electrodes to electrodeposit in irregular, nonplanar morphologies. Equilibrium theoretical calculations using joint density functional analysis in vacuum and generic liquid media indicate that in-plane transport at such interphases is enhanced substantially if LiX (X = Br > Cl > F) species predominate. This study employs optical visualization experiments and nucleation theory to experimentally investigate nucleation and early-stage growth dynamics of metallic lithium in electrolytes enriched with LiBr. It is shown that the Li-Br–rich interphases formed profoundly alter the morphology of Li electrodeposits by enhancing Li-ion surface diffusion.

25 ENERGY STORAGE↗
NASA NTRS1983

Transient nucleation in condensed systems

Using classical nucleation theory we consider transient nucleation occurring in a one-component, condensed system under isothermal conditions. We obtain an exact closed-form expression for the time dependent cluster populations. In addition, a more versatile approach is developed: a numerical simulation technique which models directly the reactions by which clusters are produced. This simulation demonstrates the evolution of cluster populations and nucleation rate in the transient regime. Results from the simulation are verified by comparison with exact analytical solutions for the steady state. Experimental methods for measuring transient nucleation are assessed, and it is demonstrated that the observed behavior depends on the method used. The effect of preexisting cluster distributions is studied. Previous analytical and numerical treatments of transient nucleation are compared to the solutions obtained from the simulation. The simple expressions of Kashchiev are shown to give good descriptions of the nucleation behavior.

Kelton, K. F.↗
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↗
DOE OSTI2022

Understanding the Importance of Labile Fe(III) during Fe(II)-Catalyzed Transformation of Metastable Iron Oxyhydroxides

Transformation of metastable Fe(III) oxyhydroxides is a prominent process in natural environments and can be significantly accelerated by the coexisting aqueous Fe(II) (Fe(II) aq ). Recent evidence points to the solution mass transfer of labile Fe(III) (Fe(III) labile ) as the primary intermediate species of general importance. However, a mechanistic aspect that remains unclear is the dependence of phase outcomes on the identity of the metastable Fe(III) oxyhydroxide precursor. Here, we compared the coupled evolution of Fe(II) species, solid phases, and Fe(III) labile throughout the Fe(II)-catalyzed transformation of lepidocrocite (Lp) versus ferrihydrite (Fh) at equal Fe(III) mass loadings with 0.2–1.0 mM Fe(II) aq at pH = 7.0. Similar to Fh, the conversion of Lp to product phases occurs by a dissolution–reprecipitation mechanism mediated by Fe(III) labile that seeds the nucleation of products. Though for Fh we observed a transformation to goethite (Gt), accompanied by the transient emergence and decline of Lp, for initial Lp we observed magnetite (Mt) as the main product. Further, a linear correlation between the formation rate of Mt and the effective supersaturation in terms of Fe(III) labile concentration shows that Fe(II)-induced transformation of Lp into Mt is governed by the classical nucleation theory. When Lp is replaced by equimolar Gt, Mt formation is suppressed by opening a lower barrier pathway to Gt by heterogeneous nucleation and growth on the added Gt seeds. The collective findings add to the mechanistic understanding of factors governing phase selections that impact iron bioavailability, system redox potential, and the fate and transport of coupled elements.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1988

Initial conditions implied by t exp 1/2 solidification of a sphere with capillarity and interfacial kinetics

The paper explores the initial conditions implied by t exp 1/2 growth of a spherical crystal solidifying from a pure, undercooled melt, including the effects of both capillarity and interface kinetics, and relates the findings to initial conditions that would be expected on the basis of classical nucleation theory. For crystal sizes near the nucleation radius, the calculated temperature profiles show a cold region ahead of the advancing interface that is even more undercooled than the undercooled bath. This cold region acts as a local heat sink that compensates for the reduced growth speed that would otherwise result from capillarity and kinetics, leading to precisely the same t exp 1/2 growth law that would have been obtained had both capillarity and kinetics been neglected.

Sekerka, R. F.↗
DOE OSTI2020

Extreme Antiscaling Performance of Slippery Omniphobic Covalently Attached Liquids

Scale formation presents an enormous cost to the global economy. Classical nucleation theory dictates that to reduce the heterogeneous nucleation of scale, the surface should have low surface energy and be as smooth as possible. Past approaches have focused on lowering surface energy via the use of hydrophobic coatings and have created atomically smooth interfaces to eliminate nucleation sites, or both, via the infusion of low-surface-energy lubricants into rough superhydrophobic substrates. Although lubricant-based surfaces are promising candidates for antiscaling, lubricant drainage inhibits their utilization. Here, we develop methodologies to deposit slippery omniphobic covalently attached liquids (SOCAL) on arbitrary substrates. Similar to lubricant-based surfaces, SOCAL has ultralow roughness and surface energy, enabling low nucleation rates and eliminating the need to replenish the lubricant. To enable SOCAL coating on metals, we investigated the surface chemistry required to ensure high-quality functionalization as measured by ultralow contact angle hysteresis (<3°). Using a multilayer deposition approach, we first electrophoretically deposit (EPD) silicon dioxide (SiO 2 ) as an intermediate layer between the metallic substrate and SOCAL. The necessity of EPD SiO 2 is to smooth (<10 nm roughness) as well as to enable the proper surface chemistry for SOCAL bonding. To characterize antiscaling performance, we utilized calcium sulfate (CaSO 4 ) scale tests, showing a 20× reduction in scale deposition rate than untreated metallic substrates. Descaling tests revealed that SOCAL dramatically decreases scale adhesion, resulting in rapid removal of scale buildup. Furthermore, our work not only demonstrates a robust methodology for depositing antiscaling SOCAL coatings on metals but also develops design guidelines for the creation of antifouling coatings for alternate applications such as biofouling and high-temperature coking.

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

Dust Formation and Destruction in Core Collapse Supernovae Explosions [Dissertation]

Dust is a crucial component of the Interstellar Medium (ISM) and an important driver in many astrophysical processes. Dust formed from evolved stars increases the metallicity of the ISM, enriching it with heavier elements. Dust makes up the source material of collapsing clouds in the ISM leading to stellar and galactic formation. The stabilization of these systems are aided by dust through absorption and scattering of excess energy. As stellar systems evolve and undergo supernova and mass loss events, enriched dust is produced, further enhancing the metallicity of the ISM. However, despite the importance of dust, the processes that produce dust are not well understood. Core-Collapse Supernovae are one source of dust. Dying massive stars explode, sending their constituent material into the Interstellar Medium, which as the ejecta expands and cools, material condenses and nucleates forming dust grains. I investigate the formation, survivability, and properties of dust. In Chapter 2, dust grain nucleation in CCSNe is discussed. A revised formulation of Kinetic Nucleation Theory (KNT) is used to track the nucleation and growth of dust grains in the ejecta. The affects of the progenitor system on the amount and type of dust is explored, showing a strong dependence on explosion energy. I discuss the code and methods used to model dust nucleation and destruction in Chapter 3. The sputtering and erosion of dust grains in CCSNe shocks are explored in Chapter 4. The erosion and survival of this nucleated dust is modeled through thermal and non-thermal sputtering by the surrounding gas and shocks. The optical properties and scattering of light by surviving dust grains is modeled by Mie scattering in Chapter 5. A discussion of observations of dust in variable star systems and supernovae is outlined in Chapter 6.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Model form and sensitivity analysis of CALPHAD-based nucleation models in b-stabilized Ti alloys

Accurate prediction of α-phase nucleation and growth in β-stabilized titanium alloys is crucial for designing heat treatments to optimize mechanical properties in additively manufactured lightweight components. Ideally, predictions of nucleation and growth would incorporate both top-down observations of past experimental heat treatments and bottom-up modeling of phase transformations; however, the appropriate method of combining these information sources is not self-evident. Combining top-down and bottom-up information requires a unified form of model that can connect between spatiotemporal scales, as well as sets of fitting parameters that can be identified by each data source. The selection of which parameters to fit to which data source can be made based on expert opinion, or by performing a sensitivity analysis. In solid-solid nucleation, direct observation of the nucleation and growth process is challenging. Most data on the heat treatment-controlled phase transformations are not in-situ. To predict the process and outcome of the nucleation, growth and coarsening of precipitates, theoretical models of the nucleation pathway are used to bridge the gap. Many sources of uncertainty affect the modeling of this nucleation process. It can be influenced by small variations in the thermomechanical processing history, chemical composition, and initial microstructure. If molecular dynamics (MD) simulations are used to determine thermodynamic quantities and inform CALPHAD modeling, additional uncertainty can be introduced and accounted for using Bayesian methods. Top-down uncertainties require additional steps to quantify. The influence of nucleation model form on the sensitivity of predictions to input parameters and physical conditions is the focus of this study. Classical nucleation theory (CNT) allows modeling to formulate the nucleation as homogeneous or, more commonly, heterogeneous. Non-classical nucleation models are also increasingly explored as a means of reconciling top-down and bottom-up data. In this study, the sensitivity of the intragranular nucleation of α in a β-annealed, slow-cooled aging (BASCA) heat treatment of β-stabilized Ti5553 alloy is explored using CNT and both heterogeneous and homogeneous assumptions. The Kampmann-Wagner Numerical model of precipitate nucleation and growth is employed. Using open-source tools (pyCalphad and thermodynamic modeling of TiMo as a surrogate system, a sensitivity analysis is performed to measure variations in key parameters, including chemical driving force, interfacial energy, and diffusivity, as they relate to predictions of precipitate number density. The inclusion of top-down and bottom-up data in selection of nucleation model form is discussed.

Rodriguez Negron, A. M.↗
DOE OSTI2023

A 1D Model for Nucleation of Ice From Aerosol Particles: An Application to a Mixed‐Phase Arctic Stratus Cloud Layer

Abstract Mixed‐phase clouds (MPCs) have been identified as significant contributors to uncertainties in climate projections, attributable to model representation of processes controlling the formation and loss of supercooled water droplets and ice particles from the atmosphere. Arctic MPCs are commonly widespread and long‐lived, with sustained ice crystal formation processes that challenge current understanding. This study examines the ice‐nucleating particle (INP) reservoir dynamics governing immersion‐mode heterogeneous freezing in an observed case of Arctic MPCs using a simplified 1D aerosol‐cloud model. The model setup includes prescribed dynamical forcings and thermodynamic profiles, and represents INPs as multicomponent and polydisperse particle size distributions. Diagnostic and prognostic approaches to immersion freezing parameterization are compared, including time‐independent (singular) number‐ and surface area‐based descriptions and a time‐dependent description following classical nucleation theory (CNT). The choice of freezing parameterization defines the size of the INP reservoir. The CNT‐based description yields an orders of magnitude larger INP reservoir than the singular parameterizations, which is the dominant factor for sustained ice crystal formation. The efficiency of the freezing process and cloud cooling are of secondary importance. A diagnostic treatment neglecting INP loss is only accurate when the INP reservoir size is large and INP depletion weak. Since a larger INP reservoir sustains ice crystal formation substantially longer, and ice water path scales with ice crystal concentrations for the conditions considered, resolving the source of differences in INP reservoir dynamics due to model implementation is a high priority for advancing climate model physics.

54 ENVIRONMENTAL SCIENCES↗