Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Operator approximation”

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 1,567 records · Page 87

Mitigating Data Center Impact on Grid Stability: A Coordinated Control Strategy Using Verrus StabiliGrid Architecture

Large data centers, which now represent a significant and growing share of the total U.S. grid load, can inadvertently destabilize the electrical grid when they disconnect simultaneously during brief voltage disturbances. The July 10, 2024, Eastern Interconnection incident, in which a sub-100-millisecond transmission fault triggered the cascading loss of approximately 1,500 MW of data center load, illustrates this vulnerability. While commercial battery energy storage systems (BESS) deployed in data centers provide device-level fault ride-through per IEEE 1547, they lack coordination with facility protection logic and uninterruptible power supplies (UPS), limiting their effectiveness as grid-stabilizing assets. This report presents the Verrus StabiliGrid architecture, a coordinated control framework that integrates BESS, UPS, and point-of-interconnection (POI) protection settings to enable data centers to ride through both undervoltage and overvoltage grid contingencies without disconnecting. The four-step strategy encompasses: (1) high-resolution power quality monitoring to detect the grid state during events such as undervoltage, overvoltage, underfrequency, and overfrequency; (2) POI protection settings that allow for extended ride-through and grid-connected operation during grid contingencies; (3) grid state-driven autonomous dispatch of assets to improve grid resilience by reducing power draw during undervoltage or absorbing more power during overvoltage events; and (4) coordinated post-recovery dispatch of data center assets to restore firm load to pre-contingency levels. Validated through controller-hardware-in-the-loop (C-HIL) simulations at the National Laboratory of the Rockies, results show grid import restoration to pre-fault levels within 100 milliseconds of voltage recovery. This work advances the ability of data centers to transition from passive, disturbance-sensitive loads to active participants in grid stability, a capability increasingly required by emerging NERC and ERCOT regulatory frameworks.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Unalakleet Microgrid Optimization for Tribal Community Resilience

The Unalakleet Microgrid Optimization Project aimed to strengthen the reliability and efficiency of the isolated electric power system that serves the Tribal community of Unalakleet, Alaska. The community relies entirely on a local wind-diesel microgrid, consisting of four 475 kW diesel generators and six 100 kW wind turbines, to provide electricity to approximately 745 residents and Tribal facilities. Because Unalakleet is not on a road system and is located nearly 400 miles from the nearest major power grid, maintaining a resilient and efficient local energy system is critical. The scope of this project included upgrading a portion of the transmission line between the wind farm and the power plant to increase voltage and reduce line losses, along with modernizing the Supervisory Control and Data Acquisition (SCADA) system to improve monitoring, control, and data management of the power system. These upgrades were designed to increase wind energy utilization, reduce diesel fuel consumption by tens of thousands of gallons annually, and improve overall grid stability. By allowing more of the community’s electricity to be supplied by local renewable wind resources, the project was designed to lower operating costs, reduce dependence on imported fuel, and strengthen the long-term resilience of the power system. These improvements represent an important step toward the community’s long-term energy vision of expanding renewable generation, incorporating energy storage, and eventually achieving “diesels-off” operation, where essential Tribal loads are powered primarily by local renewable resources.

17 WIND ENERGY↗

Integration of Total-Sky Imager Data with a Physics-Based Smart Persistence Model for Intra-Hour Forecasting of Solar Radiation

Short-term solar forecasting models based solely on global horizontal irradiance (GHI) measurements are often unable to discriminate the forecasting of the factors affecting GHI from those that can be precisely computed by atmospheric models. Our previous study introduced a Physics-based Smart Persistence model for Intra-hour forecasting of solar radiation (PSPI) that decomposed the forecasting of GHI into the computation of extraterrestrial solar radiation and solar zenith angle and the forecasting of cloud albedo and cloud fraction. The extraterrestrial solar radiation and solar zenith angle were accurately computed by the Solar Position Algorithm (SPA) developed at the National Renewable Energy Laboratory (NREL). A cloud retrieval technique was used to estimate cloud albedo and cloud fraction from surface-based observations of GHI. With the assumption of persistent cloud structures, the cloud albedo and cloud fraction were predicted for future time steps using a two-stream approximation and a 5-minute exponential weighted moving average, respectively. The model evaluation indicated the estimation and forecast of cloud fraction mostly contributed to the uncertainty of the PSPI though it overcame the persistence and smart persistence models in all forecast time horizons between 5 and 60 minutes. This study aims to enhance the PSPI by ingesting surface-based observations of cloud fraction from a total sky imager (TSI). The estimation and forecast of cloud albedo is correspondingly improved by utilizing the cloud fraction observations and thus leads to more accurate GHI forecast. Various time-series analysis methods are also investigated on the forecasting of cloud fraction and cloud albedo for further improving the GHI forecast. These improvements are valuable for many applications, such as forecasting energy use for buildings, grid operations, and ultimately bringing down the cost of solar energy.

14 SOLAR ENERGY↗

HiPACE++: A portable, 3D quasi-static particle-in-cell code

Modeling plasma accelerators is a computationally challenging task and the quasi-static particle-in-cell algorithm is a method of choice in a wide range of situations. In this work, we present the first performance-portable, quasi-static, three-dimensional particle-in-cell code HiPACE++. By decomposing all the computation of a 3D domain in successive 2D transverse operations and choosing appropriate memory management, HiPACE++ demonstrates orders-of-magnitude speedups on modern scientific GPUs over CPU-only implementations. The 2D transverse operations are performed on a single GPU, avoiding time-consuming communications. The longitudinal parallelization is done through temporal domain decomposition, enabling near-optimal strong scaling from 1 to 512 GPUs. HiPACE++ is a modular, open-source code enabling efficient modeling of plasma accelerators from laptops to state-of-the-art supercomputers.

97 MATHEMATICS AND COMPUTING↗

Comparison of multiobjective optimization methods for the $\mathrm{LCLS-II}$ photoinjector

Particle accelerators are among some of the largest science experiments in the world and can consist of thousands of components with a wide variety of input ranges. These systems can easily become unwieldy optimization problems during design and operations studies. Starting in the early 2000s, searching for better beam dynamics configurations became synonymous with heuristic optimization methods in the accelerator physics community. Genetic algorithms and particle swarm optimization are currently the most widely used. These algorithms can take thousands of simulation evaluations to find optimal solutions for one machine prototype. For large facilities such as the Linac Coherent Light Source (LCLS) and others, this equates to a limited exploration of many possible design configurations. In this paper, the LCLS-II photoinjector is optimized with three optimization algorithms. All optimizations were started from both a uniform random and Latin hypercube sample. In all cases, the optimizations started from Latin hypercube samples outperformed optimizations started from uniform samples. All three algorithms were able to optimize the photoinjector, with the model-based methods approximating the Pareto front in fewer simulation evaluations. This work, in combination with previous optimization observations, indicates objective penalties have a strong impact on the efficiency of such methods. In general, we recommend heuristic methods for initial optimizations and model-based methods when information about the objective space is available.

43 PARTICLE ACCELERATORS↗

A fast and accurate physics-informed neural network reduced order model with shallow masked autoencoder

Traditional linear subspace reduced order models (LS-ROMs) are able to accelerate physical simulations in which the intrinsic solution space falls into a subspace with a small dimension, i.e., the solution space has a small Kolmogorov n-width. However, for physical phenomena not of this type, e.g., any advection-dominated flow phenomena such as in traffic flow, atmospheric flows, and air flow over vehicles, a low-dimensional linear subspace poorly approximates the solution. To address cases such as these, we have developed a fast and accurate physics-informed neural network ROM, namely nonlinear manifold ROM (NM-ROM), which can better approximate high-fidelity model solutions with a smaller latent space dimension than the LS-ROMs. Our method takes advantage of the existing numerical methods that are used to solve the corresponding full order models. The efficiency is achieved by developing a hyper-reduction technique in the context of the NM-ROM. Numerical results show that neural networks can learn a more efficient latent space representation on advection-dominated data from 1D and 2D Burgers' equations. A speedup of up to 2.6 for 1D Burgers' and a speedup of 11.7 for 2D Burgers' equations are achieved with an appropriate treatment of the nonlinear terms through a hyper-reduction technique. Lastly, a posteriori error bounds for the NM-ROMs are derived that take account of the hyper-reduced operators.

97 MATHEMATICS AND COMPUTING↗

Quantum minimal surfaces from quantum error correction

We show that complementary state-specific reconstruction of logical (bulk) operators is equivalent to the existence of a quantum minimal surface prescription for physical (boundary) entropies. This significantly generalizes both sides of an equivalence previously shown by Harlow[1]; in particular, we do not require the entanglement wedge to be the same for all states in the code space. In developing this theorem, we construct an emergent bulk geometry for general quantum codes, defining ``areas'' associated to arbitrary logical subsystems, and argue that this definition is ``functionally unique.'' We also formalize a definition of bulk reconstruction that we call ``state-specific product unitary’’ reconstruction. This definition captures the quantum error correction (QEC) properties present in holographic codes and has potential independent interest as a very broad generalization of QEC; it includes most traditional versions of QEC as special cases. Our results extend to approximate codes, and even to the ``non-isometric codes'' that seem to describe the interior of a black hole at late times.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Exploring Fabrication of Ultra-low Absorption, Multilayer, Hafnium High Reflectors

Ultra-low absorption, high damage threshold coatings are critical for the advancement of the high power CW laser technologies, and in most cases, are a limiting factor for the output power. Since high power CW laser systems enable extremely high power irradiances, when the light is incident on an HR mirror, part of the energy is absorbed and converted into heat, which can have many detrimental side effects. The worst-case scenarios are optical distortion and laser-induced damage. These optics must operate reliably with irradiances of 100's kW/cm 2 to MW/cm 2 . The highest pulsed laser damage threshold high reflectors are consistently fabricated with hafnia material. However, hafnia is not the material of choice for high power CW laser systems because of very high absorption and scatter properties, factors of >100x higher than tantala-silica, due to its poly crystalline nature. This study focuses on CW laser use applications, where the coatings’ laser damage resistance is largely limited by the inherent absorption of the coating materials. The goal is to study the feasibility of utilizing a “composite” hafnium-silica and a “doped” hafnia-silica target materials to prohibit crystallization. The results from the study were extremely successful. Hafnia high reflectors coatings generated from this study exhibited one micron absorption values of approximately 1 ppm, which are the lowest reported values to date.

36 MATERIALS SCIENCE↗

High-burnup Experiments in Reactivity Initiated Accidents (HERA)

HERA is a joint experimental program (JEEP) operating within the Nuclear Energy Agency’s (NEA’s) framework for irradiation experiments (FIDES). HERA is dedicated to the understanding of light water reactor (LWR) fuel performance at high burnup under reactivity-initiated accidents (RIA). During RIAs the fuel cladding can be breached by one of two mechanisms; pellet-cladding mechanical interaction or high temperature balloon and rupture. Pellet-cladding mechanical interaction (PCMI) failure occurs in the early phase of an RIA transient where the fuel pellet expands rapidly into the cladding prior to any meaningful heat transfer from the pellet to the cladding. This results in high stresses in the cladding resulting in crack propagation through the base Zircaloy metal and axially along the cladding tube. Failure by PCMI is a brittle failure mode and occurs at low cladding strains, approximately 2% hoop strain. If the cladding survives the PCMI phase of the transient, it may still fail due to ballooning and rupture induced by extended exposure to high temperature. Thermal transport from the fuel pin to the coolant is significantly reduced if a transition to film boiling on the cladding surface occurs. This result in both a rapid increase in fuel temperature, which drives fission gas release, and a spike in cladding temperature, which results in loss of mechanical strength. If the internal pressure of the fuel rod is greater than the system pressure the cladding can inelastically deform and eventually rupture. Failure by swelling and rupture occurs at larger strains, generally greater than 5% hoop strain. While PCMI failures are generally seen as more limiting, the dominance of one failure mode over the other can be affected by both the cladding material conditions, and the transient evolution.

(HERA)↗

High-burnup Experiments in Reactivity Initiated Accidents (HERA)

HERA is a joint experimental program (JEEP) operating within the Nuclear Energy Agency’s (NEA’s) framework for irradiation experiments (FIDES). HERA is dedicated to the understanding of light water reactor (LWR) fuel performance at high burnup under reactivity-initiated accidents (RIA). During RIAs the fuel cladding can be breached by one of two mechanisms; pellet-cladding mechanical interaction or high temperature balloon and rupture. Pellet-cladding mechanical interaction (PCMI) failure occurs in the early phase of an RIA transient where the fuel pellet expands rapidly into the cladding prior to any meaningful heat transfer from the pellet to the cladding. This results in high stresses in the cladding resulting in crack propagation through the base Zircaloy metal and axially along the cladding tube. Failure by PCMI is a brittle failure mode and occurs at low cladding strains, approximately 2% hoop strain. If the cladding survives the PCMI phase of the transient, it may still fail due to ballooning and rupture induced by extended exposure to high temperature. Thermal transport from the fuel pin to the coolant is significantly reduced if a transition to film boiling on the cladding surface occurs. This result in both a rapid increase in fuel temperature, which drives fission gas release, and a spike in cladding temperature, which results in loss of mechanical strength. If the internal pressure of the fuel rod is greater than the system pressure the cladding can inelastically deform and eventually rupture. Failure by swelling and rupture occurs at larger strains, generally greater than 5% hoop strain. While PCMI failures are generally seen as more limiting, the dominance of one failure mode over the other can be affected by both the cladding material conditions, and the transient evolution.

(HERA)↗

Time dependent supervisory control update with FARM using rolling window

This report describes improvements to the Feasible Actuator Range Modifier (FARM) component of the RAVEN-based HYBRID framework for analysis of Integrated Energy Systems (IES). FARM supports the HERON plug-in that solves the power dispatch problem. The solution to the dispatch problem involves economically optimal dispatches that satisfy limits on production variables and their rates of variation (explicit constraints) as well as process variables tied to the service life of equipment (implicit constraints). FARM serves to validate or confirm that a HERON solution for explicit constraints also satisfies the implicit constraints. FARM-alpha was released by Argonne National Laboratory in January 2021 followed by FARM-Beta in January 2022 with the latter providing increased flexibility for the user. In this report, FARM-Gamma, the latest version of the code, is described. The major improvement is the implementation of a system identification algorithm based on the Dynamic Mode Decomposition with Control (DMDc) coupled with a “Rolling Window” scheme that allows obtaining linear time-varying state-space models. This feature equips FARM with the most accurate approximation of system dynamics, and it relieves the user from the burden of performing an exhaustive off-line characterization of the dynamics. FARM-Gamma capabilities are assessed by solving the power dispatch problem for a representative IES unit. The simulation times corresponding to the different releases are estimated and compared. These values capture the increasing computational burden of the successively higher fidelity state-space models adopted by FARM-Alpha, FARM-Beta and FARM-Gamma. The code implementation provides significant flexibility, i.e., the user can always select the most suitable version of FARM according to the problem to be solved and the available computational resources. It is anticipated that FARM will play a role in addressing several future IES applications. We outline how it can support the coordinated management and safe operation of a nuclear plant coupled to industrial processes to produce hydrogen and synfuels.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Discretization Writeup for Grey Flux-Limited Radiation Diffusion

This report documents the time and space discretizations for grey flux-limited diffusion applied to the thermal radiative transfer (TRT) equations. We begin with a description of the physics being solved before moving into the diffusion approximation. Once we have the TRT system, we show a finite-volume-inspired discretization from Jim Morel (Texas A&M University, NUEN 627 class notes, lecture 8). As systems become hotter they emit more photons in the form of blackbody radiation. Because average photon energy of the blackbody source is proportional to the temperature of the system, we call these thermal photons or thermal radiation. As material temperatures increase, increasing fractions of the total energy in the system go into the radiation field. In addition, radiation can deposit energy and momentum non-locally, making it an important phenomenon for heating and impulse. In order to accurately study systems at high temperatures, we wish to add the physics of thermal radiation to our hydrodynamic system. In practice, coupling radiation and hydrodynamics is often done by operator-splitting each timestep into two consecutive, non-overlapping phases: (1) update the hydrodynamics for a fixed radiation state (2) update the radiation and internal energy for an otherwise fixed hydrodynamic state. Because of this clean separation of physics updates, in this report we show only the latter phase, which involves solely the TRT equations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Development and Experimental Optimization of High-Temperature Modeling Tools and Methods for Concentrated Solar Power Particle - Systems

A novel, open-source radiative modeling toolset was developed to extend the functionality of particle-based modeling software (e.g. discrete element method (DEM)) to environmental conditions relevant to concentrated solar power applications. This toolset was optimized for deployment on desktop workstations instead of high-performance computing systems, to render such tools more accessible to the research community. Both particle-based modeling and radiative exchange modeling are computationally expensive and often require specialized programming expertise, making these methods cumbersome to use. Recent developments in DEM software by DCS Computing have greatly reduced these challenges, providing a graphical-user-interface based platform and modeling optimization for desktop workstations, HPCs, and cloud computing. The University of Dayton leveraged the experience of DCS Computing in developing a user-friendly, open-source radiative heat transfer expansion for DEM modeling. The University of Dayton DEM+ radiative modeling toolset was developed using a combination of fundamental experimental measurements, modeling, and simplified flow experiments over a range of temperatures and flow conditions. The toolset provides researchers with access to multiple radiative models including an accelerated Monte-Carlo Ray Tracing (application agnostic, highly computationally expensive), an expanded database of distance-based approximations (application limited, computationally light), and a weighted blending of the two methods capable of achieving over 90% reduction in computation time with equivalent accuracy compared to Monte-Carlo Ray Tracing. Through a graphical user interface, users can customize the radiative models to match their desired accuracy and available computational resources, improving access to particle based modeling for the research community. Ceramic sintered bauxite proppants were used in modeling and experimentally as a baseline. Both the radiative heat transfer and flow properties for particulate systems were investigated at elevated temperatures up to 800 °C. The major accomplishments for this work include a verified, open-source radiative modeling toolset to be distributed amongst the research community and the fabrication of three small-scale test facilities to investigate particle behavior and tune DEM flow properties for operation up to 800 °C. The findings have been shared with the research community via conference modeling workshops, deployment of the tools in DCS Computing Aspherix®, and open-source access to the developed radiative modeling tool. The development of next-generation CSP facilities and thermal energy storage systems based on ceramic particles requires providing access to computationally efficient and accurate modeling tools. Particles will experience a wide range of environments (20-800 °C) and handling conditions (dilute curtains or dense packing), requiring specially designed and optimized equipment. Optimizing solid particle physics models and establishing best-practices for particle modeling in CSP environments will assist researchers with designing optimized equipment, accelerating the deployment of more economically-competitive CSP facilities.

14 SOLAR ENERGY↗

Speeding Up Hartree–Fock in JuliaChem with Density Fitting

In this work, the density fitting (DF) approximation is added to the restricted Hartree–Fock (RHF) implementation in the JuliaChem computational chemistry code. Utilizing a DF algorithm that uses symmetry and integral screening, a significant reduction in time to compute the Fock matrix is achieved. The symmetry and screening DF-RHF techniques were adapted to be performed on graphics processing units (GPUs), which are well suited to perform the matrix multiplications that comprise the bulk of the Fock build time in DF-RHF. The JuliaChem DF-RHF GPU algorithm employs a novel approach that automatically switches between two DF-RHF algorithms depending on the number of basis functions in the calculation. The JuliaChem GPU DF-RHF implementation demonstrates up to 2× speedup for Fock build times compared to the existing best-in-class GPU DF-RHF implementation by operating directly on screened intermediate matrices. Due to the high portability of the Julia language code, the JuliaChem CPU and GPU DF-RHF implementations could be benchmarked on a variety of CPU and GPU architectures from multiple hardware vendors.

Hayes, John J. [Ames Laboratory, and Iowa State Un↗

All-perovskite tandem solar cells with improved grain surface passivation

All-perovskite tandem solar cells hold the promise of surpassing the efficiency limits of single-junction solar cells; however, until now, the best-performing all-perovskite tandem solar cells have exhibited lower certified efficiency than have single-junction perovskite solar cells. A thick mixed Pb–Sn narrow-bandgap subcell is needed to achieve high photocurrent density in tandem solar cells, yet this is challenging owing to the short carrier diffusion length within Pb–Sn perovskites. Here we develop ammonium-cation-passivated Pb–Sn perovskites with long diffusion lengths, enabling subcells that have an absorber thickness of approximately 1.2 μm. Molecular dynamics simulations indicate that widely used phenethylammonium cations are only partially adsorbed on the surface defective sites at perovskite crystallization temperatures. The passivator adsorption is predicted to be enhanced using 4-trifluoromethyl-phenylammonium (CF3-PA), which exhibits a stronger perovskite surface-passivator interaction than does phenethylammonium. By adding a small amount of CF3-PA into the precursor solution, we increase the carrier diffusion length within Pb–Sn perovskites twofold, to over 5 μm, and increase the efficiency of Pb–Sn perovskite solar cells to over 22%. Here, we report a certified efficiency of 26.4% in all-perovskite tandem solar cells, which exceeds that of the best-performing single-junction perovskite solar cells. Encapsulated tandem devices retain more than 90% of their initial performance after 600 h of operation at the maximum power point under 1 Sun illumination in ambient conditions.

14 SOLAR ENERGY↗

Non-Hermitian Quantum Mechanics Approach for Extracting and Emulating Continuum Physics Based on Bound-State-like Calculations

Here, this Letter introduces a unified emulation framework for studying continuum physics in finite quantum systems. Using a reduced basis method, we construct powerful emulators for the inhomogeneous Schrödinger equation that operate in a combined parameter space of complex energy (𝐸) and other inputs (𝜽). Within the space, the emulators simultaneously perform analytical continuation in 𝐸—extracting continuum physics from numerically simpler bound-state-like calculations—and interpolate this entire process across 𝜽. This yields a small, non-Hermitian system whose properties (e.g., resonances and scattering observables) can be rapidly predicted for any 𝜽. Crucially, the complex-𝐸 emulation provides a pathway to compute continuum observables for complex systems where advanced bound-state methods exist but direct continuum calculations are yet to be developed, while the 𝜽 emulation enables rapid parameter-space exploration and can be adapted to accelerate other existing continuum calculations. Demonstrations with two- and three-body systems highlight the method’s effectiveness and suggest its connection to (near-)optimal rational approximation. This Letter presents the key results, with further details reserved for a companion paper.

ab initio calculations↗

Grid Value Analysis of Medium-Voltage Back-to-Back Converter on DER Hosting Enhancement

This paper presents an analysis of the value that can be realized by medium-voltage back-to-back (MVB2B) converters in terms of the increased utilization rate of distributed energy resources (DERs) and the improvement in operational conditions. A systematic, transferrable, and scalable methodology has been designed to analyze and quantify the increased DER value from three perspectives: 1) curtailment reduction of the DER generation, 2) size reduction of the energy storage needed to otherwise realize DER hosting levels, and 3) hosting capacity improvement of the DERs compared to the base distribution circuit capability. In the case study, the proposed methodology is applied to two utility distribution systems for analysis and quantification of the grid value of the MVB2B converter, installed in the distribution circuit, and provided to the solar photovoltaic (PV) DERs. Here, the analysis results demonstrate that the MVB2B converter can deliver significant value to the PV hosting enhancement of two adjacent distribution systems when they are connected by the MVB2B converter. Based on this case study, this paper analyzes and summarizes the approximate realized grid value of the MVB2B converter for distribution systems dominated by different shares of customer classes.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Ice Cryo-Encapsulation Balloon (Project ICEBall) Field Campaign Report

The Ice Cryo-Encapsulation Balloon (ICEBall) field campaign was designed to sample the ice crystals that compose high-altitude cirrus with a passive device. The campaign made use of a new instrument, ICEBall, which is a balloon-borne ice crystal sampling system. The ice crystal sounding system is capable of measuring ice crystal concentration, temperature, atmospheric pressure, ice crystal habit, aerosol particle morphology, and residual composition. The 3-kg instrument is carried upwards at 5 m s-1 by a high-altitude balloon. The instrument can be cut down from the balloon at any altitude up to 20 km, and the apparatus returns to the surface by parachute. Basic measurements such as temperature and pressure are recorded onboard, and high-frequency Global Positioning System (GPS) records altitude and latitude/longitude. Ice crystal concentrations are measured through the use of a high-resolution video camera mounted on the device. Ice crystals are collected through an open aperture leading to insulated collection chambers cooled with dry ice. Upon exiting the top of the cloud system, the chamber aperture is closed, and the ~1 mm3 sample cell is magnetically sealed and isolated at -78 °C, ensuring that ice particles do not sublimate or grow after collection. Once the crystals are returned to the surface, they are double-sealed and immersed at liquid nitrogen temperature in “dry-cryo shippers” before being transported back to the laboratory. Dr. Magee’s laboratory at The College of New Jersey contains a cryo-stage scanning electron microscope (SEM), which was used to interrogate the crystals and aerosol particles. The main purpose of this pilot field campaign was to provide an unprecedented level of detail on the crystal habits and ice surface complexity in mid-latitude cirrus, which may help resolve issues associated with habit identification and classification in cirrus. The ICEBall campaign was originally scheduled to run from March 28 to April 18 of 2021 at the U.S. Department of Energy Atmospheric Radiation Measurement (ARM) Southern Great Plains observatory. The COVID-19 pandemic intervened and caused us to shift the dates of the experiment to October 16-November 6 of 2021. This period is also climatologically favorable for cirrus. The approximately six-month gap between our original field campaign dates and the actual dates afforded us the opportunity to build two new ICEBall payload instruments. These instruments were tested during an August 2021 trip to The College of New Jersey. During this field testing phase, we decided to launch the ICEBall payload upstream from the ARM SGP site with the goal of landing in the vicinity of the site. Our goal was to sample the ice crystals before the cirrus were advected over the remote-sensing instruments at the SGP site. The team assembled for the field campaign consisted of the Principle Investigator (PI) and Co-Principle Investigator (Co-PI) (Drs. Harrington and Magee), The Pennsylvania State University research scientist Dr. Alfred Moyle, and two graduate students (Ms. Marley Majetic and Gwenore Pokrifka). The team operated out of a house rented in Enid, Oklahoma. We successfully sampled seven cirrus cloud systems during the three-week field campaign (October 21, 23-26, 31, and November 1). This was a much higher success rate than either of the PIs anticipated (our goal was closer to sampling three or four cases). The balloon was typically launched from oil pads or farm fields northwest of Enid and the payload was typically retrieved somewhat north of the SGP site. We never landed directly at the SGP site, and so did not need regular access to the SGP facilities. Our greatest concern going into the field campaign was the longer-term storage of crystals in the -196°C cryo dry-shipper dewars and the subsequent transport across the country. We had tested storage and transport prior to the field campaign, but we had never stored crystals for a few weeks nor had we transported the dewars over long distances. To our great relief, the storage and transport worked flawlessly and we were able to image a large number of crystals from six of the seven cases. Working with the staff at the ARM SGP office was excellent. They not only helped us find the sources we needed for helium, liquid nitrogen, and other materials, but also helped with contacts within the Federal Aviation Administration (FAA) and Vance Air Force Base. One goal of our field project was to tie the in situ measurements of ice crystal habits to the radar signatures derived from the Ka-band ARM Zenith-pointing Radar (KAZR). Unfortunately KAZR was down for the duration of our experiment. However, the Ka-band Scanning ARM Cloud Radar (KASACR) was put into vertically pointing mode during the ICEBall campaign and those data, along with Doppler lidar measurements, have proved very useful.

54 ENVIRONMENTAL SCIENCES↗