Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “stability theory”

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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Electropolishing-induced topographic defects in niobium: insights and implications for superconducting radio frequency applications

Electropolishing (EP) is the premier surface preparation method for high-Q, high-gradient superconducting RF cavities made of Nb. This leaves behind an apparently smooth surface, yet the achievable peak magnetic fields fall well below the superheating field of Nb, in most cases. Here, in this work, the ultimate surface finish of EP was investigated by studying its effect on highly polished Nb samples. EP introduces high slope angle sloped-steps at grain boundaries. The magnetic field enhancement and superheating field suppression factors associated with such a geometry are calculated in the London theory. Despite the by-eye smoothness of electropolished Nb, such defects compromise the stability of the low-loss Meissner state, likely limiting the achievable peak accelerating fields in superconducting RF cavities. Finally, the impact of surface roughness on impurity diffusion is investigated which can link surface roughness to the effectiveness of heat treatments like low-temperature baking or nitrogen infusion in the vortex nucleation or hydride hypotheses. Surface roughness tends to decrease the effective dose of impurities as a result of the expansion of impurities into regions with greater internal angle. The effective dose of impurities can be protected by minimizing slope angles and step heights, ensuring uniformity.

Hryhorenko, Oleksandr [Thomas Jefferson National A↗
DOE OSTI2022

Fate of multiparticle resonances: From Q-balls to 3 He droplets

We study a system of N nonrelativistic particles which form a near-threshold resonance. Assuming no subset of these particles can form a bound state, the resonance can only decay through an “explosion” into N particles. We find that the decay width of the resonance scales as E Δ –5/2 in the limit when the energy E of the resonance goes to zero, where Δ is the ground-state energy of a system of N particles in a spherical harmonic trap with unit frequency. Here, the formula remains valid when some pairs of final particles have zero-energy s-wave resonance, but the Efimov effect is not present. In the limit of large N, we show that the final particles follow a Maxwell-Boltzmann distribution if they are bosons and a semicirclelike law if they are fermions. We expect our general result to be applicable to various systems that exist in nature. In particular, we argue that metastable 3 He droplets exist with the lifetime varying over many orders of magnitude ranging from a fraction of a nanosecond to values greatly exceeding the age of the Universe.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Towards Robustness Guarantees for Feedback-Based Optimization

Feedback-based online optimization algorithms have gained traction in recent years because of their simple implementation, their ability to reject disturbances in real time, and their increased robustness to model mismatch. While the robustness properties have been observed both in simulation and experimental results, the theoretical analysis in the literature is mostly limited to nominal conditions. In this work, we propose a framework to systematically assess the robust stability of feedback-based online optimization algorithms. We leverage tools from monotone operator theory, variational inequalities and classical robust control to obtain tractable numerical tests that guarantee robust convergence properties of online algorithms in feedback with a physical system, even in the presence of disturbances and model uncertainty. The results are illustrated via an academic example and a case study of a power distribution system.

approximation algorithms↗
DOE OSTI2022

Effects of experimental nitrogen deposition on soil organic carbon storage in Southern California drylands

Atmospheric nitrogen (N) deposition is enriching soils with N across biomes. Soil N enrichment can increase plant productivity and affect microbial activity, thereby increasing soil organic carbon (SOC), but such responses vary across biomes. Drylands cover ~45% of Earth's land area and store ~33% of global SOC contained in the top 1 m of soil. Nitrogen fertilization could, therefore, disproportionately impact carbon (C) cycling, yet whether dryland SOC storage increases with N remains unclear. To understand how N enrichment may change SOC storage, we separated SOC into plant–derived, particulate organic C (POC), and largely microbially derived, mineral–associated organic C (MAOC) at four N deposition experimental sites in Southern California. Theory suggests that N enrichment increases the efficiency by which microbes build MAOC (C stabilization efficiency) if soil pH stays constant. But if soils acidify, a common response to N enrichment, then microbial biomass and enzymatic organic matter decay may decrease, increasing POC but not MAOC. We found that N enrichment had no effect on C fractions except for a decrease in MAOC at one site. Specifically, despite reported increases in plant biomass in three sites and decreases in microbial biomass and extracellular enzyme activities in two sites that acidified, POC did not increase. Furthermore, microbial C use and stabilization efficiency increased in a non–acidified site, but without increasing MAOC. Instead, MAOC decreased by 16% at one of the sites that acidified, likely because it lost 47% of the exchangeable calcium (Ca) relative to controls. Indeed, MAOC was strongly and positively affected by Ca, which directly and, through its positive effect on microbial biomass, explained 58% of variation in MAOC. Finally, long–term effects of N fertilization on dryland SOC storage appear abiotic in nature, such that drylands where Ca–stabilization of SOC is prevalent and soils acidify, are most at risk for significant C loss.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Thermodynamics and kinetics of H adsorption and intercalation for graphene on 6 H -SiC(0001) from first-principles calculations

Previous experimental observations for H intercalation under graphene on SiC surfaces motivate the clarification of configuration stabilities and kinetic processes related to intercalation. From first-principles density-functional-theory calculations, we analyze H adsorption and intercalation for graphene on a 6H-SiC(0001) surface, where the system includes two single-atom-thick graphene layers: the top-layer graphene (TLG) and the underling buffer-layer graphene (BLG) above the terminal Si layer. Our chemical potential analysis shows that in the low-H coverage regime (described by a single H atom within a sufficiently large supercell), intercalation into the gallery between TLG and BLG or into the gallery underneath BLG is more favorable thermodynamically than adsorption on top of TLG. However, intercalation into the gallery between TLG and BLG is most favorable. We obtain energy barriers of about 1.3 and 2.3 eV for a H atom diffusing on and under TLG, respectively. From an additional analysis of the energy landscape in the vicinity of a step on the TLG, we assess how readily one guest H atom on the TLG terrace can directly penetrate the TLG into the gallery between TLG and BLG versus crossing a TLG step to access the gallery. Finally, we also perform density functional theory calculations for higher H coverages revealing a shift in favorability to intercalation of H underneath BLG and characterizing the variation with H coverage in interlayer spacings.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Machine learning pipeline to predict defect behavior in metallic alloy systems

The interaction between defect and solute atoms is critical to the thermodynamic and kinetic behavior of metallic alloys under exposure to high-energy radiation, causing irradiation damage in materials. Radiation can generate non-equilibrium concentrations of point defects such as vacancies and interstitials. The excess point defects not only accelerate diffusional processes such as precipitation that cause radiation embrittlement, but also change the pathway of phase transformations, including nucleation processes. Understanding these defect behaviors is complicated by the challenge and complexity of addressing each possible local and discrete distribution of environments and chemical interactions around targeted defects-solute or solute-solute complexes. To resolve the challenge, machine learning regression techniques have emerged as powerful tools that can train and construct an energy model to accurately describe the chemical interactions of solutes and defects. In Fiscal Year 2022, the work focused on the workflow development and demonstration using machine learning regression, density functional theory, cluster expansion, and Monte Carlo simulation to predict the effects of ternary solute elements (e.g., aluminum and molybdenum) and point defects on the Cr-rich $\alpha^{\prime}$ precipitation in multicomponent FeCr model alloys. The computational outcomes include the prediction of the ternary phase diagram, vacancy formation energy for different compositions, and the effect of vacancies on the nucleation of Cr-rich clusters. The simulations predict a pronounced change of Cr solubility in bcc Fe by the addition of Al and the rejection of Al atoms from $\alpha^{\prime}$ precipitates. Additionally, the simulations show the formation of Cr-vacancy clusters as the initial nuclei for stable nucleation and growth of $\alpha^{\prime}$ particles. The results demonstrate important outcomes and applications of using machine learning pipeline to study model or commercial alloys with multicomponent solute species and point defects.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Electropolishing-induced topographic defects in niobium: Insights and Implications for SRF

Electropolishing (EP) is the premier surface preparation method for high-Q, high-gradient superconducting RF cavities made of Nb. This leaves behind an apparently smooth surface, yet the achievable peak magnetic fields fall well below the superheating field of Nb, in most cases. In this work, the ultimate surface finish of EP was investigated by studying its effect on highly polished Nb samples. EP introduces high slope angle sloped-steps at grain boundaries. The magnetic field enhancement and superheating field suppression factors associated with such a geometry are calculated in the London theory. Despite the by-eye smoothness of electropolished Nb, such defects compromise the stability of the low-loss Meissner state, likely limiting the achievable peak accelerating fields in superconducting RF cavities. Finally, the impact of surface roughness on impurity diffusion is investigated which can link surface roughness to the effectiveness of heat treatments like low-temperature baking or nitrogen infusion in the vortex nucleation or hydride hypotheses. Surface roughness tends to decrease the effective dose of impurities as a result of the expansion of impurities into regions with greater internal angle. The effective dose of impurities can be protected by minimizing slope angles and step heights, ensuring uniformity.

Lechner, Eric [Thomas Jefferson National Accelerat↗
DOE OSTI2026

Transition Metal Taggants in UO 2 from First Principles

The incorporation of transition metals into nuclear fuel has gained attention both to improve fuel properties and as a possible nuclear forensics tool. Recent experimental studies by Ulrich et al. and Adorno Lopes et al. have investigated Ni and Fe as candidate transition metal dopants for potential nuclear forensics purposes and found that there is minimal alteration to key UO 2 fuel properties. In the present work, we performed density functional theory (DFT) investigations into possible defect structures for Ni and Fe incorporation. The dynamic stability of these defect structures was validated by calculating the phonon density of states. We found that transition metal incorporation likely occurs via substitution at U sites in the fluorite crystal structure with a nearby O vacancy for charge-balancing, which agrees with experimentally proposed structures. Additionally, this defect structure does not cause long-range alterations to the crystal parameters—an important consideration for use as a nuclear fuel taggant. Future needed work involves computational investigation of additional defect concentrations using larger supercells, additional transition metal charge states, and thermal effects. Such investigations can be carried forward into sintering models for a more complete understanding of the suitability of Ni and Fe as fuel taggants.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

First-principles calculations of metal surfaces. II. Properties of low-index platinum surfaces toward understanding electron emission

The stability of low-index platinum surfaces and their electronic properties is investigated with density functional theory, toward the goal of understanding the surface structure and electron emission, and identifying precursors to electrical breakdown, on nonideal platinum surfaces. Propensity for electron emission can be related to a local work function, which, in turn, is intimately dependent on the local surface structure. The (1×N) missing row reconstruction of the Pt(110) surface is systematically examined. The (1×3) missing row reconstruction is found to be the lowest in energy, with the (1×2) and (1×4) slightly less stable. In the limit of large (1×N) with wider (111) nanoterraces, the energy accurately approaches the asymptotic limit of the infinite Pt(111) surface. Furthermore, this suggests a local energetic stability of narrow (111) nanoterraces on free Pt surfaces that could be a common structural feature in the complex surface morphologies, leading to work functions consistent with those on thermally grown Pt substrates.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Pellet triggering of edge localized modes in low collisionality pedestals at DIII-D

Edge localized modes (ELMs) are triggered using deuterium pellets injected into plasmas with ITER-relevant low collisionality pedestals, and the resulting peak ELM energy fluence is reduced by approximately 25%–50% relative to natural ELMs destabilized at similar pedestal pressures. Cryogenically frozen deuterium pellets are injected from the low-field side of the DIII-D tokamak at frequencies lower than the natural ELM frequency, and heat flux is measured by infrared cameras. Ideal MHD pedestal stability calculations show that without pellet injection, these low collisionality pedestals were limited by their current density (peeling-limited) rather than their pressure gradient (ballooning-limited). ELM triggering success correlates strongly with pellet mass, consistent with the theory that a large pressure perturbation is required to trigger an ELM in low collisionality discharges that are far from the ballooning stability boundary. For sufficiently large pellets, both instantaneous and time-integrated ELM energy deposition measured by infrared cameras is reduced with respect to naturally occurring ELMs at the inner strike point, which is the position where it is largest for natural ELMs. Energy fluence at the outer strike point is less effected. Cameras observing both heat flux and D-alpha emission often find significant toroidally asymmetric striations in the outboard far scrape-off layer resulting from ELMs that are triggered by pellets. Toroidal asymmetries at the inner strike point are similar between natural and pellet-triggered ELMs, indicating that the reduction in peak heat flux and total fluence at that location is robust for the conditions reported here.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

New insights on wind turbine wakes from large‐eddy simulation: Wake contraction, dual nature, and temperature effects

Abstract Large‐eddy simulation (LES) has been adopted to study wind turbine wakes because it can capture fine‐scale details of turbulent wind flows and interactions with wind turbines. Here, we use the LES version of the Weather Research and Forecasting (WRF) model with an actuator disk model to gain insights on several wake effects that have been traditionally difficult to measure. The first finding is that the wake has a “dual nature,” meaning that the wind speed deficit behaves differently from the added turbulent kinetic energy (TKE) and the two are not co‐located in space. For example, the wind speed deficit peaks at hub height and reaches the ground within 8D (D is the rotor diameter), but added TKE peaks near the rotor tip and generally remains aloft. Second, temperature changes near the ground are driven by the added TKE in the rotor area and by atmospheric stability. The combination of these two factors determines the sign and intensity of the vertical heat flux divergence below the rotor, with convergence and warming associated with stable conditions and weak divergence and modest cooling with unstable conditions. Third, wakes do not expand indefinitely, as suggested by similarity theory applied to the wind speed deficit, but eventually stop expanding and actually contract, at different rates depending on atmospheric stability. The implication of these findings is that, in order to study wakes, it is not sufficient to focus on wind speed deficit alone, because TKE is also important and yet behaves differently from the wind speed deficit.

17 WIND ENERGY↗
DOE OSTI2021

Exomoons in Systems with a Strong Perturber: Applications to α Cen AB

The presence of a stellar companion can place constraints on occurrence and orbital evolution of satellites orbiting exoplanets, i.e., exomoons. In this work we revise earlier orbital stability limits for retrograde orbits in the case of a three-body system consisting of a star, planet, and satellite. The revised limit reads a{sub sat}{sup crit}≈0.668(1-1.236e{sub p}) for e {sub p} ≤ 0.8 in units of the Hill Radius and represents the lower critical orbit as a function of the planetary eccentricity e {sub p}. A similar formula is determined for exomoons hosted by planets in binary star systems, where e {sub p} is replaced with the components of free and forced eccentricity from secular orbit evolution theory. By exploring the dynamics of putative exomoons in α Centauri AB we find that the outer stability limit can be much less than half the Hill Radius due to oscillations in the planetary orbital eccentricity caused by the gravitational interaction with the binary star. We show, furthermore, how the resulting truncation of the outer stability limit can affect the outward tidal migration and potential observability of exomoons through transit-timing variations (TTVs). Typical TTV (rms) amplitudes induced by exomoons in binary systems are ≲10 minutes and appear more likely for planets orbiting the less massive stellar component.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Machine Learning for Improving Surface-Layer-Flux Estimates

Abstract Flows in the atmospheric boundary layer are turbulent, characterized by a large Reynolds number, the existence of a roughness sublayer and the absence of a well-defined viscous layer. Exchanges with the surface are therefore dominated by turbulent fluxes. In numerical models for atmospheric flows, turbulent fluxes must be specified at the surface; however, surface fluxes are not known a priori and therefore must be parametrized. Atmospheric flow models, including global circulation, limited area models, and large-eddy simulation, employ Monin–Obukhov similarity theory (MOST) to parametrize surface fluxes. The MOST approach is a semi-empirical formulation that accounts for atmospheric stability effects through universal stability functions. The stability functions are determined based on limited observations using simple regression as a function of the non-dimensional stability parameter representing a ratio of distance from the surface and the Obukhov length scale (Obukhov in Trudy Inst Theor Geofiz AN SSSR 1:95–115, 1946), $$z/L$$ z / L . However, simple regression cannot capture the relationship between governing parameters and surface-layer structure under the wide range of conditions to which MOST is commonly applied. We therefore develop, train, and test two machine-learning models, an artificial neural network (ANN) and random forest (RF), to estimate surface fluxes of momentum, sensible heat, and moisture based on surface and near-surface observations. To train and test these machine-learning algorithms, we use several years of observations from the Cabauw mast in the Netherlands and from the National Oceanic and Atmospheric Administration’s Field Research Division tower in Idaho. The RF and ANN models outperform MOST. Even when we train the RF and ANN on one set of data and apply them to the second set, they provide more accurate estimates of all of the fluxes compared to MOST. Estimates of sensible heat and moisture fluxes are significantly improved, and model interpretability techniques highlight the logical physical relationships we expect in surface-layer processes.

Meteorology & Atmospheric Sciences↗
DOE OSTI2025

Multi-fidelity learning for interatomic potentials: low-level forces and high-level energies are all you need

The promise of machine learning interatomic potentials (MLIPs) has led to an abundance of public quantum mechanical (QM) training datasets. The quality of an MLIP is directly limited by the accuracy of the energies and atomic forces in the training dataset. Unfortunately, most of these datasets are computed with relatively low-accuracy QM methods, e.g. density functional theory with a moderate basis set. Due to the increased computational cost of more accurate QM methods, e.g. coupled-cluster theory with a complete basis set (CBS) extrapolation, most high-accuracy datasets are much smaller and often do not contain atomic forces. The lack of high-accuracy atomic forces is quite troubling, as training with force data greatly improves the stability and quality of the MLIP compared to training to energy alone. Because most datasets are computed with a unique level of theory, traditional single-fidelity (SF) learning is not capable of leveraging the vast amounts of published QM data. In this study, we apply multi-fidelity learning (MFL) to train an MLIP to multiple QM datasets of different levels of accuracy, i.e. levels of fidelity. Specifically, we perform three test cases to demonstrate that MFL with both low-level forces and high-level energies yields an extremely accurate MLIP—far more accurate than a SF MLIP trained solely to high-level energies and almost as accurate as a SF MLIP trained directly to high-level energies and forces. Therefore, MFL greatly alleviates the need for generating large and expensive datasets containing high-accuracy atomic forces and allows for more effective training to existing high-accuracy energy-only datasets. Indeed, low-accuracy atomic forces and high-accuracy energies are all that are needed to achieve a high-accuracy MLIP with MFL.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Uptake of Pb and the Formation of Mixed (Ba,Pb)SO 4 Monolayers on Barite During Cyclic Exposure to Lead-Containing Sulfuric Acid

Barite (BaSO 4 ) is a common additive in lead-acid batteries, where it acts as a nucleating agent to promote the reversible formation and dissolution of PbSO 4 during battery cycling. However, little is known about the molecular-scale mechanisms the nucleation and cyclic evolution of PbSO 4 over a battery’s lifetime. In this study, we explore the responses of a barite (001) surface to cycles of high and low lead concentrations in the presence of 100 mM sulfuric acid using in-situ atomic force microscopy and high-resolution X-ray reflectivity. Here we find that PbSO 4 epitaxial films readily nucleate on the barite surface, even from solutions that are undersaturated relative to bulk PbSO 4 . Despite this, barite (001) proves to be an ineffective nucleator of bulk PbSO 4 , as multilayer growth is suppressed even in highly supersaturated solutions. Instead, we find evidence that Pb 2+ ions can directly exchange with Ba 2+ to create mixed (Ba,Pb)SO 4 surfaces. These chemically mixed surfaces do not host PbSO4 monolayers as readily as pristine barite, and the original reactivity is not regained until a fresh surface is reestablished by aggressive etching. Our results can be partly explained by traditional models of Stranski-Krastanov (S-K) growth, in which monolayer films are stabilized by a reduction in surface energy, but multilayer growth is inhibited by epitaxial strain. Complementary density functional theory calculations confirmed the basic energetic-terms of traditional S-K models, but also showed evidence form more complex, thickness-dependent, energetics than would be predicted from the S-K models. The experimental results can be better understood by extending the S-K model to consider the formation of mixed surfaces and films, which can have reduced strain and interfacial energies relative to pure films, while also being stabilized by entropy of mixing. These insights into non-stoichiometric heteroepitaxy will enable better predictions of how barite affects PbSO 4 nucleation in battery environments.

25 ENERGY STORAGE↗
DOE OSTI2020

Two-Dimensional Conductive Ni-HAB as a Catalyst for the Electrochemical Oxygen Reduction Reaction

Catalytic systems whose properties can be systematically tuned via changes in synthesis conditions are highly desirable for the next-generation catalyst design and optimization. Herein, we present a two-dimensional (2D) conductive metal–organic framework consisting of M-N 4 units (M = Ni, Cu) and a hexaaminobenzene (HAB) linker as a catalyst for the oxygen reduction reaction. By varying synthetic conditions, we prepared two Ni-HAB catalysts with different crystallinities, resulting in catalytic systems with different electric conductivities, electrochemical activity, and stability. We show that crystallinity has a positive impact on conductivity and demonstrate that this improved crystallinity/conductivity improves the catalytic performance of our model system. Additionally, density functional theory simulations were performed to probe the origin of M-HAB’s catalytic activity, and they suggest that M-HAB’s organic linker acts as the active site with the role of the metal being to modulate the linker sites’ binding strength.

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

Multi-modal characterization of the B2 phase in the Ta-Re binary system

The energy and transportation industries demand materials that retain their mechanical property at high temperatures. Refractory complex concentrated alloys (RCCAs) with a BCC + B2 microstructure offer a potential solution, where maintaining the high temperature mechanical properties can be achieved by precipitation strengthening. This depends on the B2 phase in RCCAs being thermodynamically stable with a high solvus temperature. Recently, we predicted the high temperature stability of the B2 structure in the Ta-Re binary system, using density functional theory. Here, in this study, we provide experimental evidence for the existence of this phase for the first time, using a Ta 65 Re 35 alloy. Despite Ta-Re binary phase diagrams predicting a single-phase BCC microstructure for Ta 65 Re 35 , we show that a high Z nanoscale secondary phase appears after heat treatment at 1550 °C and 1100 °C. Scanning transmission electron microscopy (STEM) revealed that this phase has a cubic structure and is equiatomic TaRe though B2 superlattice reflections were absent in fast Fourier transforms (FFT) and diffraction patterns (DPs). DP simulations indicate that the B2 TaRe superlattice reflections are up to two orders of magnitude weaker than their fundamental reflections, making their detection challenging via electron microscopy. Neutron diffraction confirmed the second phase had a B2 structure. This study identified a previously unobserved high temperature stable B2 phase in the Ta-Re system, enabling the development of new high temperature BCC + B2 RCCAs.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Theory of In-Cloud Activation of Aerosols and Microphysical Quasi-Equilibrium in a Deep Updraft

The microphysical quasi-equilibrium in an ascending adiabatic parcel is elucidated by an analytical theory in 0D with drastic simplifications. The theory predicts how in-cloud activation is most likely to be triggered by the onset of precipitation during sufficient ascent, with the ascent only needing to approach almost twice the cloud-base updraft speed aloft. The precipitation and cloud condensate mass fields are coupled in a closed system in a simplified version of the model. The initial state of no precipitation is unstable with respect to a perturbation. In the 2D phase space of both mass fields, there is a neutral line. Unstable growth of precipitation mass drives the system to cross the line into a regime of stability. An attractor is then approached consisting of precipitation mass balanced by accretion of cloud mass and fallout. The cloud-particle number concentration also approaches a stable equilibrium governed by the balance between in-cloud activation aloft from the inexorably increasing supersaturation during vertical acceleration and accretion of cloud droplets by precipitation. Dimensionless numbers characterizing the microphysical equilibria and their stability are derived mathematically, including a condensation–precipitation efficiency and an in-cloud activation efficiency. The theory explains common observations of the orders of magnitude of liquid water content in convective and stratiform clouds. Finally, sensitivity tests of the numerically integrated theoretical equations are documented with respect to variations in cloud condensation nucleus (CCN) aerosols and updraft speed. In conclusion, it is shown that this theory of in-cloud activation, with an increasing supersaturation during ascent, applies to both ice-only and liquid-only cloud.

54 ENVIRONMENTAL SCIENCES↗