Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Time-dependent modeling”

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
DOE OSTI2026

Testing Recommendations for Thermal Failure Analysis

Minimizing ranges of uncertainty in failure QOIs is vital for making informed engineering decisions. To improve our decision-making ability, it is necessary to have a QOI representation that captures more complex failure modes. Capturing the physics of complex time-dependent temperature failure in explicit detail can be computationally and experimentally expensive. To avoid this blooming of complication and costs, an empirical QOI component-level failure/damage model may be sufficient to help identify component-level failure. In this context, failure and damage are abstractions of more complicated processes that lead to the cessation of normal operating behavior, not a representation of mechanical damage, such as void formation. Such a model requires minimal increases in computational cost and experimental burden. The proposed QOI serves as an indicator that something has likely failed, similar to the previous critical failure temperatures or pressures. This QOI will not replace any existing failure QOIs. It will be calculated alongside them, allowing for preservation of legacy style QOI measurements with new context from the new QOI.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Overview of the physics basis for the ARC fusion power plant

Commonwealth Fusion Systems plans to build ARC as the first fusion power plant at a site in Chesterfield County, Virginia, USA by the early 2030s. We present an overview of analysis comprising the physics basis of the ARC V3A design, a high-magnetic-field tokamak with 𝐵 0 =11.4 T, 𝐼 𝑝 =12.0 MA, 𝑅 0 =4.62 m, 𝑎 =1.18 m. ARC V3A is designed to produce 𝑃 𝑓⁡𝑢⁢𝑠 ≈1.13 GW DT fusion power and deliver ⩾400 MW net electric power to the grid. This overview includes quantitative analysis of fundamental issues for design of and operational plasma scenarios for a tokamak power plant, and lays out the design targets and strategic choices for ARC, including empirical fusion performance projections, assessment of H-mode access, ion cyclotron resonance heating simulations, alpha particle physics and time-dependent full-pulse simulations. This is complemented by topical papers on fusion performance and transport, disruption physics, boundary physics and magnetohydrodynamic stability. Critically, these studies identify key model uncertainties and physics risks to be retired through SPARC operation. Due to the modular nature of ARC, early results from SPARC can be incorporated into the design of the first ARC as well as subsequent replacements of the ARC vacuum vessel.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Integration of Online Cross-Section Generation Capability with Depletion and Transient Solvers in Griffin

Griffin is a Multiphysics Object-Oriented Simulation Environment (MOOSE)-based reactor multiphysics analysis application jointly developed by Argonne and Idaho National Laboratories under the DOENE Nuclear Energy Advanced Modeling and Simulation (NEAMS) program. In FY25, an online crosssection generation capability based on the Self-Shielding Application Programming Interface (SSAPI) was demonstrated for TRISO-fueled reactor problems under steady-state conditions. This fiscal year, that capability was extended to support depletion and transient multiphysics calculations, enabling high-fidelity analyses that generate self-shielded cross sections on the fly from the actual evolving composition and temperature states rather than from pre-tabulated libraries. For depletion, a two-way coupling was established in which SSAPI computes compact-averaged self-shielded cross sections that the depletion solver then uses to advance the Bateman equations, with the updated compositions returned to SSAPI at each step; the depletion module was refactored to support both library-based and SSAPI-based cross sections, and additional logic was added to track daughter isotopes and to exclude minor isotopes for efficiency. For transient analysis, the SSAPI multigroup library was extended with the kinetics data required for time-dependent calculations, the Improved Quasi-Static (IQS) scheme was coupled with SSAPI, and several supporting capabilities were implemented, including a self-shielding treatment that lets control rods and drums move within a self-shielded model, which had previously been impossible and had ruled out rod- and drum-movement transients with on-the-fly cross sections altogether, a new mixing scheme for delayed-neutron precursor decay constants, a checkpoint-based restart workflow, and performance improvements such as pointwise cross-section interpolation and the bypassing of unnecessary Dancoff factor calculations. The implemented capabilities were verified against Serpent Monte Carlo solutions. For depletion, a prismatic pin-cell problem based on a Next Generation Nuclear Plant (NGNP) Very High Temperature Reactor benchmark showed excellent agreement, with eigenvalue differences within 200 pcm over the entire burnup range (up to 140 MWD/kgU) and fission-product and actinide inventories agreeing to within 0.8% and 2.5%, respectively; a heat-pipe microreactor assembly problem with a much higher fuel loading confirmed the same behavior and quantified the bias introduced when the multigroup equivalence effect is neglected. For transient analysis, a pin-cell problem with a step reactivity insertion and temperature feedback reproduced the analytically expected asymptotic power and showed close agreement between the direct and IQS solutions, and a two-dimensional microreactor core problem with control-drum rotation exercised the new moving-drum self-shielding treatment and demonstrated successful coupling of the online crosssection generation with both the direct and IQS transient methods. The capability was further exercised on a full-core pebble-bed problem, in which Griffin was coupled with the System Analysis Module (SAM) to simulate load-following operation of the gPBR with the Doppler feedback resolved at the TRISO fuel kernel temperature. These developments in Griffin provide a convenient, high-fidelity approach to cross-section generation for advanced thermal reactors with geometrically complex and highly heterogeneous configurations, including TRISO-fueled prismatic and pebble-bed systems, and support steady-state, depletion, and transient multiphysics calculations. They also enable self-shielded cross sections to be evaluated directly at the actual coupled state of the system, thereby establishing a foundation for high-fidelity, fully coupled multiphysics analysis of advanced reactors

Park, H.↗
DOE OSTI2025

Egg yolk as a model for gelation: From rheometry to flow physics

Egg yolks are an excellent model for studying sol-gel transitions, particularly the power law viscoelasticity that defines the critical point of gelation. However, prior studies lack comprehensive datasets and fail to visualize flow behavior linked to temperature and time-dependent linear and nonlinear rheology. Here, we present a detailed dataset characterizing egg yolk viscoelasticity across temperature, time, and forcing amplitude using oscillatory shear, step strain, step stress, and constant high strain rate. Novel protorheology visualizations link rheological properties with observable flow behavior. Our findings highlight the nuanced determination of the critical gel point, emphasizing observation timescale dependencies. We compare methods to identify critical temperatures for gelation, including power law viscoelasticity, moduli crossover, diverging zero-shear viscosity, and emerging equilibrium elastic modulus, while visualizing flow consequences near these transitions. Egg yolk is an accessible, realistic, and nontoxic material relevant to the physicist and the chef alike, making it ideal for understanding the rheology of critical gels. By integrating protorheology photos and videos with rigorous rheometric data, we deepen the understanding of critical gels, with broader impacts for studying other materials with sol-gel transitions.

Marsh, Maxwell C. [Department of Mechanical Scienc↗
DOE OSTI2026

New correlated numerical methods for attosecond molecular single and double ionization

This project resulted in the development of ASTRA, a new molecular ionization code capable of describing single and double ionization processes in polyatomic molecules with correlated electrons and time-dependent fields. Key achievements include: (i) the development of the PRISM hybrid-integral library enabling core-ionization and high-angular-momentum convergence; (ii) implementation of time-dependent Schrödinger equation solvers for pump–probe spectroscopies; (iii) validation against experimental and theoretical benchmarks for atoms and molecules; and (iv) extension to two-photon double ionization using the finite-pulse virtual-sequential model. These developments establish ASTRA as a versatile platform for attosecond molecular dynamics.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2026

A Full-Induction Magnetohydrodynamics Solver for Liquid Metal Fusion Blankets in Vertex-CFD

Multiphysics modeling of liquid metal fusion blankets, which produce tritium and convert energy of neutrons created via fusion reactions into heat, is crucial for predicting performance, ensuring structural integrity, and optimizing energy production. While traditional blanket modeling of liquid metal flows during normal steady operating conditions commonly employs the inductionless approximation of the magnetohydrodynamics (MHD) equations, transient scenarios, when the plasma-confining magnetic field varies on millisecond time scales, require a full-induction MHD approach that dynamically evolves the magnetic field via the time-dependent induction equation. This paper presents the formulation, implementation, and initial verification of a full-induction MHD solver integrated within the open-source Vertex-CFD framework, which aims to achieve tight multiphysics coupling, a flexible software design enabling easy extension and addition of physics models, and performance portability across computing platforms. The solver utilizes finite element spatial discretization, implicit Runge–Kutta time integration, and an inexact Newton method to solve the resulting discrete nonlinear system, leveraging Trilinos packages for efficient computation. Verification against selected benchmark problems demonstrates accuracy and robustness of the solver. Furthermore, when the solver is applied to an idealized blanket model in 2.5D and full 3D, results obtained with Vertex-CFD are in good agreement with recently published quasi-2D simulations. These findings establish a computational foundation for future simulations of transient MHD phenomena in liquid metal blankets with Vertex-CFD, and open avenues for future extensions and performance optimizations.

Endeve, Eirik [ORNL] (ORCID:0000000312519507)↗
DOE OSTI2026

Muon-induced fission of actinide nuclei

A negative muon captured by an actinide cascades down through the muonic atomic levels; deeply bound transitions can proceed via inverse internal conversion, depositing the muonic transition energy directly into the nucleus and, when the deposited energy exceeds the fission barrier, inducing prompt fission. Because the muon mean lifetime exceeds the saddle-to-scission timescale by orders of magnitude, the muon can survive the entire fission event as a 1⁢𝑠 spectator and ultimately attach to one or both of the emerging fragments. Its postscission attachment probability to the light fragment, 𝑃 𝐿 , can be used as a direct electromagnetic probe of fission dynamics on a timescale of 10 −21 s. In previous work, we introduced a three-dimensional lattice solution of the time-dependent Dirac equation coupled to the electromagnetic field generated by a fissioning nucleus and reported 𝑃 𝐿 for several actinides at a single dissipation strength. In this work, we extend that framework to a systematic survey of 232 Th , 238 U , and 240 Pu and implement a more realistic fission model which incorporates dynamic pairing correlations. We find that 𝑃 𝐿 falls steeply with the fragment charge asymmetry, a robust structural fingerprint of the fissioning system, while its dependence on nuclear dissipation is secondary and sensitive to the phenomenological friction prescription. These results establish 𝑃 𝐿 as a clean electromagnetic probe of fragment charge asymmetry and motivate a self-consistent, coordinate- and time-dependent treatment of nuclear dissipation as the natural next step.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Total absorption spectroscopy of two isomers in 70 Cu influencing nucleosynthesis signatures

Isomers have long been known to be important for astrophysical nucleosynthesis processes, yet they are often neglected in network calculations due to computational limitations or lack of data. "Astromers" are astrophysically metastable nuclear states that can greatly impact nucleosynthesis pathways. In this work we show that astromers further impact the time-dependent electromagnetic signal during heavy element nucleosynthesis. In an experiment performed at the National Superconducting Cyclotron Laboratory, three 𝛽-decaying states of 70 Cu (6 − ground state, and two isomeric states: 3 − and 1 + ) were produced. 𝛽-feeding values were extracted from experimental spectra and compared to shell-model and QRPA+PVC calculations. Here, average 𝛾-ray energies from the 𝛽-decay events were incorporated into simulations of heavy element nucleosynthesis and were found to exhibit different energy release profiles over time, which may impact, in aggregate, time-dependent observable signals.

59 ≤ A ≤ 89↗
NASA NTRSDate not supplied

Battery Pack Shape Optimization using Transient Heat Conduction Coupled with Cell-Discharge Analysis

Battery electric systems exhibit significant time-dependence, especially when evaluated in the context of an aircraft mission profile with continually changing power demands. Additionally, when evaluating battery-powered aircraft concepts, it is important to accurately compute the temperature of the batteries and properly characterize the thermal response of the system. The temperature of the batteries has a significant impact on cell performance, in addition to safety considerations of maintaining battery temperatures below their operating limit. Because of these considerations, battery models for preliminary design and optimization of aircraft should include the capability to accurately compute the temperature distribution within the battery pack. Furthermore, battery pack designs should be as light-weight as possible to maximize the pack energy density, while also considering battery temperature limits. Here, we demonstrate a simultaneous trajectory and shape optimization of a battery pack concept, using a transient heat transfer finite element model coupled with a time-varying cell-discharge battery model to provide this capability. Including the transient finite element problem in the loop enables accurate temperatures that can be passed back to the cell discharge model, while the cell discharge model can supply the finite element model with time-varying heat boundary conditions to the finite element problem, further benefiting the fidelity of the thermal response of the batteries. We first demonstrate the coupling capability between the battery cell-discharge model and the transient finite-element heat transfer through an optimization which computes the optimal current profile for the battery pack while ensuring the battery temperatures remain below their operational limit. We then build on this optimization by adding shape optimization to the problem, which allows us to consider a composite objective function which also minimizes the mass of the battery pack, while also producing an optimal current discharge profile.

Optimization↗
DOE OSTI2026

Trigonometric continuous-variable gates and hybrid quantum simulations of the sine-Gordon model

Hybrid qubit-qumode quantum computing platforms provide a natural setting for simulating interacting bosonic quantum field theories. However, existing continuous-variable gate constructions rely predominantly on polynomial functions of canonical quadratures. In this work, we introduce a complementary universality paradigm based on trigonometric continuous-variable gates, which enable a Fourier-like representation of bosonic operators and are particularly well suited for periodic and non-perturbative interactions. We present an ancilla-based framework for implementing trigonometric gates with arguments given by arbitrary Hermitian functions of qumode quadratures. The protocol yields unitary gates deterministically, and non-unitary gates through probabilistic post-selection. As a concrete application, we develop a hybrid qubit-qumode quantum simulation of the lattice sine-Gordon model. Using these gates, we prepare ground states via quantum imaginary-time evolution, simulate real-time dynamics, compute time-dependent vertex two-point correlation functions, and extract quantum kink profiles under topological boundary conditions. Our results demonstrate that trigonometric continuous-variable gates provide a physically natural framework for simulating interacting field theories on near-term hybrid quantum hardware, while establishing a parallel route to universality beyond polynomial gate constructions. We expect that the trigonometric gates introduced here to find broader applications, including quantum simulations of condensed matter systems, quantum chemistry, and biological models.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Dynamics and observational signatures of core-collapse supernovae with central engines: hydrodynamics simulations with Monte Carlo post-processing

A long-lived central engine embedded in expanding supernova ejecta can alter the dynamics and observational signatures of the event, producing an unusually luminous, energetic, and/or rapidly evolving transient. We use 2D hydrodynamics simulations to study the effect of a central energy source, varying the amount, rate, and isotropy of the energy deposition. We post-process the results with a time-dependent Monte Carlo radiation transport code to extract observational signatures. The engine excavates a bubble at the centre of the ejecta, which becomes Rayleigh–Taylor unstable. Sufficiently powerful engines are able to break through the edge of the bubble and accelerate, shred, and compositionally mix the entire ejecta. The breakout of the engine-driven wind occurs at distinct rupture points, and the outflowing high-velocity gas may eventually give rise to radio emission. The dynamical impact of the engine leads to faster rising optical light curves, with photon escape facilitated by the faster expansion of the ejecta and the opening of low-density channels. For models with strong engines, the spectra are initially hot and featureless, but later evolve to resemble those of broad-line Ic supernovae. Under certain conditions, line emission from ionized, low-velocity material near the centre of the ejecta may be able to escape and produce narrow emission similar to that seen in interacting supernovae. We discuss how variability in the engine energy reservoir and injection rate could give rise to a heterogeneous set of events spanning multiple observational classes, including the fast blue optical transients, broad-line Ic supernovae, and superluminous supernovae.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Diabatization with Electrostatic Embedding for Studying Photophysics in Organic Molecular Crystals

Highly emissive organic molecular crystals find applications in several areas, such as organic electronics, solar cells, and sensors. Understanding the excited-state mechanisms underlying these applications is essential for optimizing and controlling them effectively. Exciton models coupled with nonadiabatic dynamics, particularly quantum dynamics, provide crucial insights into photochemical and photophysical processes in molecular crystals. Nevertheless, there remains a lack of general tools and automated workflows to facilitate such simulations. In this paper, we present a computational strategy to investigate the photoactivated dynamics of organic molecular crystals, bridging methodologies traditionally used for molecular systems and materials science, with a particular focus on the interplay between local excitations and charge transfer (CT) processes. We have implemented an interface between the fromage and Overdia programs, enabling the construction of vibronic Hamiltonians for molecular crystals within an excited-state ONIOM(QM:QM′) framework, incorporating long-range electrostatics through a RESP-based Ewald summation. Fragment-based diabatization provides a route to quantum dynamics simulations in weak-to-intermediate coupling regimes. The method was applied to the photophysics of dibenzo[g,p]chrysene (DBC) crystals using time-dependent DFT. The fromage/ Overdia interface was employed to compute the couplings of local excitations and CT states for 18 unique DBC dimers in the crystal and to quantify the influence of electrostatic embedding, which was found to be modest (10−20%). Simulations on π-stacked dimers reproduced the small red shift observed experimentally from solution to crystal, attributed to electronic interactions among fixed monomers rather than crystal electrostatics. Quantum dynamics simulations revealed ultrafast population transfer from bright local excitations to CT states. This approach establishes a robust framework linking molecular and solid-state excited-state dynamics, with potential applications for studying excitations, defects, and impurities in molecular crystals.

Crystals↗
DOE OSTI2026

Processing-Dependent Structure and Poroelasticity of Nafion in Liquid Water

Ionomers act as the solid polymer electrolyte membrane in many modern electrochemical devices, yet the role of their nanostructure in modulating the poroelastic response remains poorly understood, especially in liquid water, where few techniques can measure simultaneous transport-mechanical properties. Poroelastic Relaxation Indentation (PRI) is uniquely suited for measuring time-dependent transport-mechanical properties of porous solids, specifically hydraulic diffusivity, elastic modulus, Poisson’s ratio, and intrinsic permeability, for porous solids. While ionomers such as Nafion are not porous in the typical sense, Nafion has a nanophase-segregated structure that, when fully swollen in liquid water, behaves as a poroelastic solid with a coupled mechanical-transport response. Using a poroelastic framework, we investigate how casting and pretreatment of Nafion membranes alter their poroelastic response in liquid environments. We characterize both extruded and dispersion-cast Nafion membranes pretreated in water at 25 or 100 °C and relate the mechanical-transport properties to the ionomer structure via hydrophilic and intercrystalline domain spacing measured using Small-Angle X-ray Scattering (SAXS). At 25 °C, dispersion-cast membranes exhibit distinctly lower hydraulic diffusivity and intrinsic permeability than extruded membranes but with comparable mechanical properties. Pretreatment at 100 °C increases hydrophilic domain spacing, improving transport but at the expense of mechanical stiffness. Dispersion-cast membranes respond more strongly to pretreatment than extruded membranes. Using the Carman-Kozeny pore network model and the hydrophilic domain-spacing, we estimate the pore radius and tortuosity to show how pretreatment reduces structure-related tortuosity differences between dispersion-cast and extruded membranes. Here, in this work, we show that nanophase-segregated materials such as Nafion can be rigorously characterized using poroelasticity, resulting in physically meaningful transport-mechanical properties. Coupling PRI with SAXS provides insights into the nanostructural features that govern the coupled mechanical-transport response. By establishing PRI for a nanophase-segregated material, this approach opens avenues for this technique’s application in other hydrated polymeric materials not typically considered to be poroelastic.

Shen, Margaret [University of California, Berkeley↗
DOE OSTI2026

TANTE: Time-adaptive operator learning via neural Taylor expansion

Operator learning for time-dependent partial differential equations (PDEs) has seen rapid progress in recent years, enabling efficient approximation of complex spatiotemporal dynamics. However, most existing methods rely on fixed time step sizes during rollout, which limits their ability to adapt to varying temporal complexity and often leads to error accumulation. In this work, we propose the Time-Adaptive Transformer with Neural Taylor Expansion (TANTE), a novel operator-learning framework that produces continuous-time predictions with adaptive step sizes. TANTE predicts future states by performing a Taylor expansion at the current state, where neural networks learn both the higher-order temporal derivatives and the local radius of convergence. This allows the model to dynamically adjust its rollout based on the local behavior of the solution, thereby reducing cumulative error and improving computational efficiency. We demonstrate the effectiveness of TANTE across a wide range of PDE benchmarks, achieving superior accuracy and adaptability compared to fixed-step baselines, delivering accuracy gains of 60-80 % and speed-ups of 30-40 % at inference time.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Relative Source Time Functions, Spectral Ratios, and Near‐Source Spallation in the Source Physics Experiment Phase I Chemical Explosions

The Source Physics Experiment (SPE) Phase I was composed of six chemical explosions at the Nevada National Security Sites (NNSS) between 2011 and 2016. The experiment provided a robust set of dense, local to regional distance observations suitable for characterizing shallow chemical explosions located within the same borehole. We investigate the time-dependent source processes of each chemical explosion using Relative Source Time Function (RSTF) estimation, derived from five linear seismometer profiles located between 100 and 2000 m of the shot location. The RSTF estimate provide a detailed characterization of each chemical explosion's time-history. Subsequent modeling efforts suggest these measurements may be useful for precise characterization of explosion processes and spallation. RSTF estimation allows us to expand our understanding of the SPE Phase I chemical explosion series. Application of these techniques contributes to our understanding of explosion source physics and facilitates future applications to additional seismic source experiments and characterization of explosion phenomenology.

58 GEOSCIENCES↗
DOE OSTI2026

Quantum effects on the dynamics and properties of soft materials

The quantum effects of nuclear and electronic motion play an important role in the structure, dynamics, and function of soft materials, yet they are difficult to capture with conventional classical simulations or static electronic–structure methods. In this work several complementary approaches for treating quantum effects in polymeric and soft–matter systems are demonstrated, with a focus being on the hydrogen-bonded networks, ion and charge transport, and photoactive chromophores. The proton transfer, tunneling, and isotope effects are captured within the reduced-dimensionality models by implementing grid-based nuclear quantum dynamics in terms of the discrete variable and Fourier bases. The nuclear quantum dynamics is extended to larger systems by employing the quantum trajectories and quantum–thermal bath schemes combined with on-the-fly electronic structure, enabling the description of high-dimensional polymeric environments at feasible cost. The dynamics in the electronic degrees of freedom, simulating the optical response in large chromophores such as chlorophylls, is performed using the real-time time-dependent density functional theory implemented in the real-space multigrid (RMG) code. These approaches are demonstrated on case studies of the proton and hydroxide transport in hydrated polymer membranes, charge transfer in conjugated polymers, and the optical spectra of chlorophyll chromophores relevant to polymerized chlorophyll materials and chlorophyll–polymer hybrids. The reviewed methods and applications highlight practical routes of including quantum effects in simulations of soft functional materials.

Garashchuk, Sophya [Univ. of South Carolina, Colum↗
DOE OSTI2026

Fermionic mean-field dynamics for spin systems beyond free fermions

We introduce the fermionized time-dependent Hartree–Fock (fTDHF), a real-time quantum dynamics method for spin-1/2 Hamiltonians following their mapping to fermions via the Jordan-Wigner transformation. fTDHF is formally equivalent to exact dynamics in the case of free fermions, and can efficiently handle non-local string operators arising from long-range interactions via transition matrix elements between non-orthogonal Slater determinants. We show that the fTDHF method can be implemented on a classical computer with a cost that scales polynomially with system size, and linearly with the time steps. We benchmark fTDHF against exact dynamics on three separate spin-1/2 models, representing adiabatic preparation of states with long-range correlations, disorder-driven observation of many-body localization, and particle production in the Schwinger model. For each of these systems, fTDHF is shown to reproduce the qualitative dynamics generated by the exact evolutions, while maintaining a simple physical picture due to its mean-field nature.

Dutta, Rishab↗
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↗