Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “domination number”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Spectroscopic Insight on Neodymium Solvation in Lithium Borohydride-Supported Electrolyte

Borohydride-based electrolytes have recently emerged as promising media for the electrodeposition of electropositive metals, including rare earth (RE) elements. While the presence of supporting alkali metal cations and RE counteranions provides essential electrochemical conductivity for achieving fast metal electrodeposition, interactions between the host ligand and solvated neodymium (Nd) complexes remain unclear. This study provides insights into the coordination structure of a concentrated and directly solvated Nd salt in a lithium borohydride-supported electrolyte. Our spectroscopic results indicate that the RE coordination environment is significantly influenced by the solvation mechanism, which can vary between metathesis and complexation pathways, primarily dictated by stoichiometric factors. Under dilute conditions, nearly complete metathesis of anions leads to a high coordination number for the host ligand (borohydride), consistent with the previously reported solvated Nd speciation in chlorine-free electrolytes. In contrast, concentrated dissolution of the Nd salt in the supported electrolyte is dominated by a complexation pathway featuring a Li-ion-paired complex with a low coordination number of the host ligand. Density functional theory (DFT) calculations indicated that the observed blue shift in the borohydride vibration was the result of an increase in electron density drawn into the terminal B–H interbond region from the hydride as the coordination changed from Li to Nd. In conjunction with DFT results, vibrational analyses allowed correlation of the experimental shifts associated with changes in Nd ligation and coordination spheres, further consolidating the prevalence of highly chloride-coordinated species under concentrated conditions. In conclusion, the outcomes of this work illuminate the distinctive and heterogeneous coordination structures that the electroactive RE species can adopt at high concentrations in lithium borohydride-supported electrolytes, as a key step to comprehend the reported metal electrodeposition performance in these media.

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

Dust-void evolution driven by turbulent dust flux can induce runaway migration of Earth-mass planets

Torques from asymmetric dust structures (so-called dust-void and filamentary structures) formed around low-mass planets embedded in a nonturbulent dust-gas disk can exceed the torques produced by the gas disk component and then go on to dominate the planet’s orbital dynamics. Here, we investigate how these structures (hence the dust torque) change when the effect of turbulent dust diffusion and dust feedback are included, along with the direct implications on the migration of Earth-like planets. Using the FARGO3D code, we performed 2D and 3D multifluid hydrodynamic simulations, focusing on a non-migrating planet with a mass of M p = 1.5 M ⊕ in 2D and on migrating planets with M p ∈ [1.5, 12] M ⊕ in 3D. We varied the δ-dimensionless diffusivity parameter in the range [0, 3 × 10 −3 ] and considered three different Stokes numbers, St = {0.04, 0.26, 0.55}, which are representative of the gas-dominated, the transitional, and the gravity-dominated regimes, respectively. In our 2D models, we find that turbulent diffusion of dust prevents the formation of the dust-void and filamentary structures when δ > 3 × 10 −4 . Otherwise, dust structures survive turbulent diffusion flow. However, dust and total torques become positive only in transitional and gravity-dominated regimes. In our 3D models, we find that the dust-void is drastically modified and the high-density ring-shaped barrier delineating the dust-void disappears if δ ≳ 10 −4 , due to the effect of dust turbulent diffusion along with the back-reaction of the dust. For all values of δ, the filament in front of the planet is replaced by a low-density trench. Remarkably, as we have allowed the planets to migrate, the evolving dust-void can drive either runaway migration or outward (inward) oscillatory-torque migration. Our study thus suggests that low-mass Earth-like planets can undergo runaway migration in dusty disks.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Light nuclei collectivity from s NN = 3 GeV Au+Au collisions at RHIC

In high-energy heavy-ion collisions, partonic collectivity is evidenced by the constituent quark number scaling of elliptic flow anisotropy for identified hadrons. A breaking of this scaling and dominance of baryonic interactions is found for identified hadron collective flow measurements in sNN = 3 GeV Au+Au collisions. In this paper, we report measurements of the first- and second-order azimuthal anisotropic parameters, v1 and v2, of light nuclei (d, t, 3He, 4He) produced in sNN = 3 GeV Au+Au collisions at the STAR experiment. An atomic mass number scaling is found in the measured v1 slopes of light nuclei at mid-rapidity. For the measured v2 magnitude, a strong rapidity dependence is observed. Unlike v2 at higher collision energies, the v2 values at mid-rapidity for all light nuclei are negative and no scaling is observed with the atomic mass number. Calculations by the Jet AA Microscopic Transport Model (JAM), with baryonic mean-field plus nucleon coalescence, are in good agreement with our observations, implying baryonic interactions dominate the collective dynamics in 3 GeV Au+Au collisions at RHIC.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Grain structure and texture selection regimes in metal powder bed fusion

Additive manufacturing (AM) offers opportunities to produce complex part geometries not possible with conventional processing and in some cases even improve part performance. However, adoption has been slowed by difficulties assessing microstructure variability and there is no straightforward approach to relate processing to grain structure characteristics. In this study, datasets from AdditiveFOAM heat transport simulations of laser powder bed fusion (LPBF) are used to drive ExaCA simulations of grain structure. The GPU utilization of ExaCA and an algorithmic update for modeling melt pool overlap region solidification enabled rapid and parallel simulation across a wider range of process conditions than previously explored with cellular automata-based solidification models. A texture selection angle $θ_s$ is defined based on melt pool overlap geometry, and the range of $θ_s$ over which a commonly observed texture transition occurs in characterized AM builds was well-reproduced by ExaCA simulations over a wide range of melt pool shape, hatch spacing, and layer height. ExaCA simulations with 90 degree rotation of the scan direction on every other layer reproduced a number of trends from the AM literature including grain refinement, the dominance of layers with larger melt pools on the final grain structure, and the weakening or strengthening of texture depending on odd and even layer melt pool overlap geometry. EBSD data from a benchmark AM part is used to validate the simulated mechanism of a layer rotation-induced texture strengthening effect. Importantly, these results expand the understanding of the mechanisms for texture selection in alloys with cubic crystal symmetry and offer an approach to easily evaluate processing conditions. With this new understanding, these modeling tools will enable anticipation of previously unexpected variations in grain structure and target specific microstructures and properties.

36 MATERIALS SCIENCE↗
DOE OSTI2022

High repetition rate exploration of the Biermann battery effect in laser produced plasmas over large spatial regions

Abstract In this paper we present a high repetition rate experimental platform for examining the spatial structure and evolution of Biermann-generated magnetic fields in laser-produced plasmas. We have extended the work of prior experiments, which spanned over millimeter scales, by spatially measuring magnetic fields in multiple planes on centimeter scales over thousands of laser shots. Measurements with magnetic flux probes show azimuthally symmetric magnetic fields that range from 60 G at 0.7 cm from the target to 7 G at 4.2 cm from the target. The expansion rate of the magnetic fields and evolution of current density structures are also mapped and examined. Electron temperature and density of the laser-produced plasma are measured with optical Thomson scattering and used to directly calculate a magnetic Reynolds number of $1.4\times {10}^4$ , confirming that magnetic advection is dominant at $\ge 1.5$ cm from the target surface. The results are compared to FLASH simulations, which show qualitative agreement with the data.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

An–imidophosphorane (An = U–Pu) bond covalency and proton-coupled electron transfer thermodynamics driven by orbital energy matching

A series of mid-actinide (An = U–Pu) tetrahomoleptic complexes supported by highly electron-donating imidophosphorane ligands, NPC ([NP t Bu(pyrr) 2 ] − , where t Bu = C(CH 3 ) 3 ; pyrr = pyrrolidinyl = N(C 4 H 8 )), are systematically investigated computationally and experimentally to elucidate the nature of actinide–ligand (An–L) covalency across the An 3+/4+/5+ oxidation states. Trends in An–L bonding and redox properties for these complexes, together with their protonated counterparts, are examined using orbital-, electron density-, and energy-decomposition-based methods. This integrated approach reveals progressively improved energy matching between α-spin An 5f and N im 2p orbitals with increasing atomic number and oxidation state, becoming particularly pronounced in the ligand-dominant π-bonding orbitals of An 4+ and An 5+ . In contrast to the An 3+ species, the enhanced An 5f π contributions in the higher-valent counterparts drive the increase in An–N im covalency for later An, thereby inverting the covalency trend to U < Np < Pu. Redistribution of electron density towards the An and N im atomic basins due to the growing energy-matching assisted covalency correlates with higher pKa values and increased N im –H bond dissociation free energies in protonated An 4+ complexes. Electron density at Nim in An 4+ shows a linear correlation with the p K a values calculated via the Bordwell equation. Calculations predict a cathodic shift of 0.84–1.00 V in the redox couples upon protonation, a trend validated when experimentally accessible. These findings demonstrate an increasing role of covalency driven by orbital energy matching from U to Pu in tuning the thermodynamic driving force for proton-coupled electron transfer in the An 5+ species.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2021

Predicting nonresonant pressure-driven MHD modes in equilibria with low magnetic shear

Nonresonant internal modes can be difficult to anticipate as there is no resonant surface in the plasma. However, equilibria that are unstable to multiple nonresonant magnetohydrodynamic (MHD) modes may be more prone to global loss of confinement since these instabilities generate spatially extended linear displacements, potentially enhancing magnetic field line chaos via nonlinear interactions. Here, we successfully predict the unstable nonresonant pressure-driven modes for equilibria with zero shear in the plasma core, irrational q on axis, and a central pressure gradient, which is consistent with pre-crash profiles in sawtoothing tokamak plasmas in the large-aspect-ratio limit. A criterion for identifying nonresonant modes most likely to be unstable is developed from the convergents of the continued fraction representation of q 0 . A higher-order analysis of the standard Energy Principle reveals the conditions under which these modes are expected to dominate. Linear growth rate spectra, as a function of toroidal mode number (up to n = 30), calculated using the initial-value extended-MHD code, M3D-C1, recover the characteristic dependence observed for ideal infernal modes. Nonresonant modes have also been invoked in some ideal sawtooth crash models. Finally, this work provides a mechanism to predict the mode numbers of infernal modes and, potentially, the width of some post-sawtooth-crash profiles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Instantons in chiral magnets

In this work, we exhaustively construct instanton solutions and elucidate their properties in one-dimensional antiferromagnetic chiral magnets based on the $\textit{O}(3)$ nonlinear sigma model description of spin chains with the Dzyaloshinskii-Moriya (DM) interaction. By introducing an easy-axis potential and a staggered magnetic field, we obtain a phase diagram consisting of ground-state phases with two points (or one point) in the easy-axis dominant cases, a helical modulation at a fixed latitude of the sphere, and a tricritical point allowing helical modulations at an arbitrary latitude. We find that instantons (or skyrmions in two-dimensional Euclidean space) appear as composite solitons in different fashions in these phases: temporal domain walls or wall-antiwall pairs (bions) in the easy-axis dominant cases, dislocations (or phase slips called merons) with fractional instanton numbers in the helical state, and isolated instantons and calorons living on the top of the helical modulation at the tricritical point. We also show that the models with DM interaction and an easy-plane potential can be mapped into those without them, providing a useful tool to investigate the model with the DM interaction.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

A pion-argon cross section measurement in the ProtoDUNE-SP experiment with cosmogenic muon

Neutrinos are tiny mysterious fundamental particles with small cross sections. Through neutrino physics, scientists across the world are trying to answer many intriguing questions about nature such as the dominance of matter over antimatter, CP violation in the lepton sector, number of supernovas in the early universe, etc. Detection of neutrinos requires massive particle detectors and intense neutrino beam owing to their small cross section. Deep Underground Neutrino Experiment (DUNE) is a next-generation neutrino experiment that is planned to start taking data beginning in 2026. DUNE will consist of 4 massive detectors, the first of which will be using single-phase liquid argon time projection chamber (LArTPC) technology. The ProtoDUNE-SP experiment is a prototype of the DUNE built at the CERN neutrino platform and uses the same detector technology that will be used in DUNE first module. The ProtoDUNE-SP experiment collected months of test beam and cosmic ray data beginning in September 2018. It was built to provide a testbed for the installation of detector parts for DUNE, showing long-term stability of the detector, understanding detector response for different test beam particles (including protons, pions, electrons, kaons, muons), and measurement of hadron-argon cross sections. When a particle passes through LArTPC electron-ion pairs are produced. To reconstruct the position and energy of a particle passing through the medium knowledge of ionization electron drift velocity is essential. The electron drift velocity is distorted by an excess positive charge built up in the detector, known as space charge. This study discusses a novel technique for measuring the ionization electron drift velocity using cosmic-ray muons. The technique uses tracks that travel the entire drift distance of the TPC for drift velocity determination. Secondly, the study discusses a method for converting the charge deposited into energy. The method is carried out in two step s. In th e first step detector response for energetic cosmic ray muons crossing the entire the TPC is used to make the charge deposition uniform throughout the TPC, and in the second step stopping cosmic-ray muons are used for determining the energy scale. Finally, the study discusses a pion-argon cross section measurement based on reweighting of Monte Carlo simulations using J. Calcutt's Geant4Reweight framework. Neutrinos cannot be directly detected; they are identified based on the interaction products. Pions are a common interaction product in a neutrino interaction. For precise modeling of neutrino event generators, it is essential to understand the pion-argon interaction. Pion-argon cross section measurement serves as an important input for neutrino interaction models. The results of the pion-argon total reaction cross section using the Geant4 reweighting technique are found to be in good agreement with Geant4 predictions. The many studies carried out in the ProtoDUNE-SP experi ment wil l be useful for current and future neutrino experiments using LArTPC technology including ICARUS, MicroBooNE, DUNE

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

The diurnal susceptibility of subtropical clouds to aerosols

The diurnal susceptibility of clouds and their radiative properties to aerosols is examined during their Lagrangian transition from subtropical stratocumulus to shallow cumulus regimes. Using large-eddy simulations, we analyze the six-day evolution of an air mass along a 3800 km observed trajectory from the coast of Peru toward the equator. Pristine and polluted scenarios are simulated with forcing imposed from weather reanalysis. The polluted scenario exhibits stronger diurnal variations in cloud water, cloud fraction, and albedo, with enhanced nighttime entrainment and suppressed precipitation. The overall response of cloud properties and outgoing shortwave radiation to droplet number concentration follows a distinct diurnal pattern: strong positive cloud adjustments dominate at night and in the morning, while weak negative adjustments prevail in the afternoon. This cycle is driven by the competition between precipitation suppression, which enhances cloud water and coverage, and entrainment drying, which depletes them. In polluted conditions, enhanced entrainment leads to a deeper and more decoupled boundary layer that cannot be sustained by surface fluxes in the afternoon, resulting in negative cloud adjustments. The enhanced entrainment rate under polluted conditions is caused by the reduced sedimentation of cloud and precipitation water from the entrainment zone. While the Twomey effect dominates the diurnal average albedo response, the diurnal variation in the competing cloud adjustments lead to a near-neutral net susceptibility in the afternoon, highlighting the critical role of diurnally varying processes in aerosol-cloud interactions.

Geosciences↗
DOE OSTI2024

Stress-controlled medium-amplitude oscillatory shear (MAOStress) of PVA–Borax

We report the first-ever complete measurement of MAOStress material functions, which reveal that stress can be more fundamental than strain or strain rate for understanding linearity limits as a function of Deborah number. The material used is a canonical viscoelastic liquid with a single dominant relaxation time: polyvinyl alcohol (PVA) polymer solution cross-linked with tetrahydroborate (Borax) solution. We outline experimental limit lines and their dependence on geometry and test conditions. These MAOStress measurements enable us to observe the frequency dependence of the weakly nonlinear deviation as a function of stress amplitude. The observed features of MAOStress material functions are distinctly simpler than MAOStrain, where the frequency dependence is much more dramatic. The strain-stiffening transient network model was used to derive a model-informed normalization of the nonlinear material functions that accounts for their scaling with linear material properties. Moreover, we compare the frequency dependence of the critical stress, strain, and strain-rate for the linearity limit, which are rigorously computed from the MAOStress and MAOStrain material functions. While critical strain and strain-rate change by orders of magnitude throughout the Deborah number range, critical stress changes by a factor of about 2, showing that stress is a more fundamental measure of nonlinearity strength. This work extends the experimental accessibility of the weakly nonlinear regime to stress-controlled instruments and deformations, which reveal material physics beyond linear viscoelasticity but at conditions that are accessible to theory and detailed simulation.

Mechanics↗
DOE OSTI2024

Inverter Intensive Hybrid Power Plant Modeling with Small-Signal Stability Augmentation through Flexible Operation Mode Transition

Hybrid power plants (HPPs) prompt the penetration of inverter-based renewable energy sources (RES) in transmission systems; however, given their low-inertia nature, HPPs are dominated by power electronic inverters, so there are inevitable challenges in system stability when increasing numbers of HPPs are integrated into the modern power grids. To boost the penetration level of HPPs without jeopardizing system stability, it is desirable to equip them with operational characteristics (i.e., grid-forming capabilities) that are comparable to those of conventional power plants dominated by synchronous generators (SGs). In this paper, a holistic model of inverter intensive HPP is derived and a bi-level hierarchical control is developed to allow HPPs to flexibly switch among the designed operation modes (i.e., P-Q, P-V, and isochronous modes). Compared to SGs, which have limited controllability, the operation mode of each HPP could vary as requested. Such flexible mode transitions could be integrated into the secondary plant-level control and be leveraged as an additional control degree to further augment system stability. Further, the system small-signal stability margin is quantified with varying HPP operation modes. More importantly, modal analysis is thereby conducted to quantify the impacts of mode transition on the system oscillatory modes. The effectiveness of the proposed HPP bi-level hierarchical control is verified using extensive case studies based on the simplified real-world island power grid, and the results validate that the system small-signal stability margin can be enhanced with the additional degree of control flexibility enabled by HPP operation mode transition. The real-time hardware-in-the-loop (HIL) results are also provided to verify the proposed method.

grid-forming control↗
DOE OSTI2022

Cannibalism's lingering imprint on the matter power spectrum

The early universe may have contained internally thermalized dark sectors that were decoupled from the Standard Model. In such scenarios, the relic dark thermal bath, composed of the lightest particle in the dark sector, can give rise to an epoch of early matter domination prior to Big Bang Nucleosynthesis, which has a potentially observable impact on the smallest dark matter structures. This lightest dark particle can easily and generically have number-changing self-interactions that give rise to “cannibal” behavior. We consider cosmologies where an initially sub-dominant cannibal species comes to temporarily drive the expansion of the universe, and we provide a simple map between the particle properties of the cannibal species and the key features of the enhanced dark matter perturbation growth in such cosmologies. Finally, we further demonstrate that cannibal self-interactions can determine the small-scale cutoff in the matter power spectrum even when the cannibal self-interactions freeze out prior to cannibal domination.

cosmological perturbation theory↗
DOE OSTI2021

Real-time visualization of particle evolution during reactive flux-assisted processing of aluminum melts

Here, a multi-modal, multi-scale correlative tomography investigation of Al-TiC metal matrix composites processed via flux-assisted reaction synthesis is reported. Synchrotron X-ray microradiography is utilized to visualize the reaction and particle evolution in real-time. Changes in particle diameter and areal number density suggest that the process is nucleation- rather than growth-dominated. At 950 °C, the bulk of the reaction takes place in a relatively short time span of less than 600 s. The microstructure is imaged at higher resolution in 2D (scanning electron microscopy) and 3D (synchrotron X-ray nanotomography), revealing the formation of carbide particles with a hexagonal platelet morphology. We propose that the morphology arises due to the incorporation of Si impurities during the experiment. It is expected that the correlative tomography workflow and analysis may guide future metal matrix composite (MMC) processing strategies.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Assessing Metal Ion Assignment Accuracy in Protein Data Bank Models via Elemental Spectroscopy

Accurate representation of metal ions in macromolecular structures is critical for chemical interpretation, computational modeling, and machine-learning methods that rely on Protein Data Bank (PDB) entries. However, the elemental identity of metals modeled in crystallographic structures is often inferred indirectly and rarely validated experimentally. Here, we combine Particle Induced X-ray Emission (PIXE) and X-ray Fluorescence Spectroscopy (XRFS) to determine the elemental composition of protein samples used to generate 70 deposited metalloprotein crystal structures. By analyzing the original protein material employed for crystallization, but before the addition of crystallization buffer solutions, we assess whether the modeled metal ions in deposited structures are consistent with experimentally detectable elemental content. We find that in a majority of cases, the metals modeled in the corresponding PDB entries are inconsistent with the metals present in the protein samples before crystallization, or that additional metals are present but not represented in the structural models. Spectroscopic results were integrated with automated crystallographic validation metrics, including real-space Z-difference (RSZD) analysis and systematic rerefinement, to evaluate atomic-number mismatch at metal sites. PIXE and XRFS show strong agreement for dominant elemental signals and provide complementary, scalable approaches for identifying suspect metal assignments. This work does not address physiological or functional metalation but instead highlights a widespread data integrity issue in deposited macromolecular structures, PDB-wide. These results establish an experimentally corroborated link between elemental identity and crystallographic validation metrics, enabling the large-scale detection of chemically inconsistent annotations in structural databases used for computational modeling and machine learning.

Crystallization↗
DOE OSTI2022

Anomalous Second Harmonic Generation from Atomically Thin MnBi 2 Te 4

MnBi 2 Te 4 is a van der Waals topological insulator with intrinsic intralayer ferromagnetic exchange and A-type antiferromagnetic (AFM) interlayer coupling. Theoretically, it belongs to a class of structurally centrosymmetric crystals whose layered antiferromagnetic order breaks inversion symmetry for even layer numbers, making optical second harmonic generation (SHG) an ideal probe of the coupling between the crystal and magnetic structures. Here, we perform magnetic field and temperature dependent SHG measurements on MnBi 2 Te 4 flakes ranging from bulk to monolayer thickness. We find that the dominant SHG signal from MnBi 2 Te 4 is unexpectedly unrelated to both magnetic state and layer number. As a result, we suggest that surface SHG is the likely source of the observed strong SHG, whose symmetry matches that of the MnBi 2 Te 4 -vacuum interface. Our results highlight the importance of considering the surface contribution to inversion symmetry breaking in van der Waals centrosymmetric magnets.

2D magnets↗
DOE OSTI2021

Onset of Plasmoid Reconnection during Magnetorotational Instability

The evolution of current sheets in accretion flows undergoing magnetorotational instability (MRI) is examined through two- and three-dimensional numerical modeling of the resistive MHD equations in global cylindrical geometry. With an initial uniform magnetic field aligned in the vertical (z) direction, MRI produces radially extended toroidal (azimuthal) current sheets. In both 2D and 3D when axisymmetric modes dominate, these current sheets attract each other and merge in the poloidal (rz) plane, driving magnetic reconnection when the Lundquist number S > 3 × 10 2 , making it a possible source of plasmoids (closed magnetic loops) in accretion disks. At high Lundquist numbers in the 2D regime, starting at S = 5 × 10 3 , self-consistent MRI-generated current sheets become thin and subject to plasmoid instability, and therefore spontaneous magnetic reconnection. When nonaxisymmetric 3D modes dominate, turbulence makes the azimuthal current sheets more unstable and stretch vertically. Toroidally extended vertical current sheets in the inner region, as well as larger 3D magnetic islands in the outer regions of the disks are also formed. These findings have strong ramifications for astrophysical disks as potential sources of plasmoids that could cause local heating, particle acceleration, and high energy EM radiation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Impacts of secondary ice production on Arctic mixed-phase clouds based on ARM observations and CAM6 single-column model simulations

For decades, measured ice crystal number concentrations have been found to be orders of magnitude higher than measured ice-nucleating particle number concentrations in moderately cold clouds. This observed discrepancy reveals the existence of secondary ice production (SIP) in addition to the primary ice nucleation. However, the importance of SIP relative to primary ice nucleation remains highly unclear. Furthermore, most weather and climate models do not represent SIP processes well, leading to large biases in simulated cloud properties. This study demonstrates a first attempt to represent different SIP mechanisms (frozen raindrop shattering, ice–ice collisional breakup, and rime splintering) in a global climate model (GCM). The model is run in the single column mode to facilitate comparisons with the Department of Energy (DOE)'s Atmospheric Radiation Measurement (ARM) Mixed-Phase Arctic Cloud Experiment (M-PACE) observations. We show the important role of SIP in four types of clouds during M-PACE (i.e., multilayer, single-layer stratus, transition, and frontal clouds), with the maximum enhancement in ice crystal number concentrations up to 4 orders of magnitude in moderately supercooled clouds. We reveal that SIP is the dominant source of ice crystals near the cloud base for the long-lived Arctic single-layer mixed-phase clouds. The model with SIP improves the occurrence and phase partitioning of the mixed-phase clouds, reverses the vertical distribution pattern of ice number concentrations, and provides a better agreement with observations. The findings of this study highlight the importance of considering SIP in GCMs.

54 ENVIRONMENTAL SCIENCES↗