Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reproducible science”

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,135 records · Page 63
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Physics or Knob‐Tuning? Tropical Anvil Peak Is Captured by GCMs

Tropical anvil clouds peak near 200 hPa and significantly impact Earth's climate, yet its physical realism in coarse‐resolution General Circulation Models (GCMs) remains debated. We examine anvil cloud formation by performing simulations with a GCM with a hierarchy of cloud fraction schemes ranging from a complex prognostic Tiedtke scheme to a simple binary scheme. All schemes consistently reproduce the anvil peak. The robust anvil peak arises because extremely cold temperatures at the upper troposphere facilitate frequent saturation events, producing clouds that disproportionately influence mean cloud fraction. Sensitivity experiments with enhanced evaporation of cloud condensate unexpectedly show increased anvil coverage, highlighting how slight evaporative moistening reinforces local saturation in cold upper‐tropospheric conditions. These results demonstrate that the tropical anvil cloud peak emerges from fundamental thermodynamic constraints, rather than specific cloud fraction parameterization choices.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

A simple and practical wax-encapsulation method for air-sensitive XAS samples

To facilitate X-ray absorption spectroscopy (XAS) measurements of air-sensitive samples, we present a simple method in which materials are encased in common paraffin wax to protect them from air and moisture. We demonstrate the efficacy of this approach using a highly reducing, air- and moisture-sensitive uranium(III) complex, the tris(amide) U[N(SiMe 3 ) 2 ] 3 (1). When finely dispersed in a boron nitride matrix and subsequently encased in inert paraffin wax, samples of 1 remain stable with no visible or spectroscopic degradation after several days under ambient conditions. The viability of this method for XAS measurements was further evaluated across a series of uranium compounds, ranging from uranyl species to highly air- and moisture-sensitive molecular complexes, at the uranium L 3 -edge. Edge energy determinations were highly reproducible (±0.1 eV between replicates) and, where available, showed excellent agreement with literature values. In conclusion, this low-cost, effective, and versatile method offers a viable solution for XAS studies of air-sensitive compounds and materials.

36 MATERIALS SCIENCE↗
NASA NTRS2022

Perturbations in Stratospheric Aerosol Evolution Due to the Water-Rich Plume of the 2022 Hunga-Tonga Eruption

The January 2022 Hunga Tonga-Hunga Ha'apai (HTHH) volcanic eruption injected a relatively small amount of SO2, but significantly more water into the stratosphere than previously seen in the modern satellite record. Here we show that the large amount of water resulted in large perturbations to stratospheric aerosol evolution. Our Community Earth System Model simulation reproduces the enhanced water vapor observed by the Microwave Limb Sounder at pressure levels between 10 and 50 hPa for three months. Compared with a simulation without a water injection, this additional source of water vapor increases OH, which halves the SO2 lifetime. Subsequent coagulation creates larger sulfate particles that double the stratospheric aerosol optical depth. A seasonal forecast of volcanic plume transport in the southern hemisphere indicates this eruption will greatly enhance the aerosol surface area and water vapor near the polar vortex until at least October 2022, suggesting that there will continue to be an impact of the HTHH eruption on the climate system.

Atmospheric science↗
DOE OSTI2022

UrbanPop: A spatial microsimulation framework for exploring demographic influences on human dynamics

Ensuring the social equity of planning measures in social systems requires an understanding of human dynamics, particularly how individual relationships, activities, and interactions intersect with individual needs. Spatial microsimulation models (SMSMs) support planning for human security goals by representing human dynamics through realistic, georeferenced synthetic populations, that a) provide a complete representation of social systems while b) also protecting individual privacy. In this paper, we present UrbanPop, an open and reproducible SMSM framework for analysis of human dynamics with high spatial, temporal, and demographic resolution. UrbanPop creates synthetic populations of demographically detailed worker and student agents, positioning them first at probable nighttime locations (home), then moving them to probable daytime locations (work/school). Summary aggregations of these populations match the granular detail available at the census block group level in the American Community Survey Summary File (SF), providing realistic approximations of the actual population. UrbanPop users can select particular demographic traits important in their application, resulting in a highly tailored agent population. We first lay out UrbanPop's baseline methodology, including population synthesis, activity modeling, and diagnostics, then demonstrate these capabilities by developing case studies of shifting population distributions and high-risk populations in Knox County, TN during the global COVID-19 pandemic.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Prediction of neutron production and energy spectrum by the inverse kinematic reaction between an incident 7 Li 3+ beam and a proton target in PHITS

A neutron source using the inverse kinematic reaction between lithium and proton, p( 7 Li, n) 7 Be, achieves forward-directed neutrons, potentially enhancing neutron yield in the forward direction. Despite the advantage, no evaluated-cross-section data for this reaction can be used in Monte Carlo simulation codes, such as PHITS. To solve this problem, this study aims to evaluate the applicability of the user-defined cross-section data, Frag data, for p ( 7 Li, n) 7 Be in PHITS. The simulations reproduced collisions between 7 Li 3+ ions and polypropylene targets. The Frag data was edited based on the JENDL-5 by utilizing the two-body collision kinematics. The neutron yield and angular distribution were investigated in the simulation. As a result, the forward neutron convergence with a reasonable neutron yield and energy spectrum was observed. The expected neutron yield in the forward 1-steradian area is 2.46 × 10 10 n/s when lithium-ion energy and current are 16.45 MeV and 0.1 mA.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Alternating-bias assisted annealing of amorphous oxide tunnel junctions

Superconducting quantum bits (qubits) rely on ultra-thin, amorphous oxide tunneling barriers that can have significant inhomogeneities and defects as grown. This can result in relatively large uncertainties and deleterious effects in the circuits, limiting the scalability. Finding a robust solution to the junction reproducibility problem has been a long-standing goal in the field. Here, we demonstrate a transformational technique for controllably tuning the electrical properties of aluminum-oxide tunnel junctions. This is accomplished using a low-voltage, alternating-bias applied individually to the tunnel junctions, with which resistance tuning by more than 70% can be achieved. The data indicates an improvement of coherence and reduction of two-level system defects. Transmission electron microscopy shows that the treated junctions are predominantly amorphous, albeit with a more uniform distribution of alumina coordination across the barrier. This technique is expected to be useful for other devices based on ionic amorphous materials.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Extreme confinement of hydrogen gas within fullerenelike nanoporous carbon

Nanoporous carbons and carbon nanostructures can store hydrogen at cryogenic temperatures but lack the volumetric and gravimetric capacity to be industrially significant. Recent inelastic neutron scattering experiments suggest a highly dense phase of hydrogen at temperatures well above the melting point of solid hydrogen. However, it remains unclear how pore geometry and intermolecular interactions enable these dense phases to exist, with dispersion (van der Waals) or electrostatic/induction suggested to be the key effects in slit and curved pores but their relative contributions have yet to be quantified. In this paper, we perform benchmark electronic structure calculations allowing the interactions between planar and curved aromatic molecules with hydrogen to be accurately determined. Dispersion was found to dominate over electrostatic and inductive effects with some many-body charge transfer (Dobson type-A) effects needed to capture the most highly curved structures. Density functional methods that include type-A many-body effects were found to accurately describe the intermolecular interactions at a fraction of the cost of coupled-cluster simulations and these approaches were used to calculate the energies inside large carbon bowl and slit pores. The interaction energies inside the bowl pores were found to depend on the orientation of the hydrogen molecule. This rotational barrier, modeled as a quantum hindered rotor, could reproduce the peak splitting observed in inelastic neutron scattering experiments, with weak splitting arising from bowl-like fullerene pores and strong splitting from highly confining nanotubelike pores. Increasing the fraction of such curved pores in nanoporous carbons may therefore offer a pathway to enhance their hydrogen-storage capacity. Moreover, the preferential adsorption of ortho hydrogen on nanotubelike pores could enable the storage of high-density hydrogen without the need to remove heat produced during the ortho-para hydrogen conversion.

36 MATERIALS SCIENCE↗
NASA NTRS2001

Chemical Engineering in Space

The aerospace industry has long been perceived as the domain of both physicists and mechanical engineers. This perception has endured even though the primary method of providing the thrust necessary to launch a rocket into space is chemical in nature. The chemical engineering and chemistry personnel behind the systems that provide access to space have labored in the shadows of the physicists and mechanical engineers. As exploration into the cosmos moves farther away from Earth, there is a very distinct need for new chemical processes to help provide the means for advanced space exploration. The state of the art in launch systems uses chemical propulsion systems, primarily liquid hydrogen and liquid oxygen, to provide the energy necessary to achieve orbit. As we move away from Earth, there are additional options for propulsion. Unfortunately, few of these options can compare to the speed or ease of use provided by the chemical propulsion agents. It is with great care and significant cost that gaseous compounds such as hydrogen and oxygen are liquefied and become dense enough to use for rocket fuel. These low-temperature liquids fall within a specialty area known as cryogenics. Cryogenics, the science and art of producing cold operating conditions for use on Earth, in orbit, or on some other nonterrestrial body, has become increasingly important to our ability to travel within our solar system. The production of cryogenic fuels and the long-term storage of these fluids are necessary for travel. As our explorations move farther away from Earth, we need to address how to produce the necessary fuels to make a round-trip. The cost and the size of these expeditions are extreme at best. If we take everything necessary for our survival for the round-trip, we invalidate any chance of travel in the near future. As with the early explorers on Earth, we need to harvest much of our energy and our life support from the celestial bodies. The in situ production of these energy sources is paramount to success. We are currently working on several processes to produce the propellants that would allow us to visit and explore the surface of Mars. The capabilities currently at our disposal for launching and delivering equipment to another planet or satellite dictate that the size and scale of any hardware must be extremely small. The miniaturization of the processes needed to prepare the in situ propellants and life support commodities is a real challenge. Chemical engineers are faced with the prospect of reproducing an entire production facility in miniature so the complex can be lifted into space and delivered to our destination. Another area that does not normally concern chemical engineers is the extreme physical aspects payloads are subjected to with the launch of a spacecraft. Extreme accelerations followed by the sudden loss of nearly all gravitational forces are well outside normal equipment design conditions. If the equipment cannot survive the overall trip, then it obviously will not be able to yield the needed products upon arrival. These launch constraints must be taken into account. Finally, we must consider both the effectiveness and efficiencies of the processes. A facility located on the Moon or Mars will not have an unlimited supply of power or other ancillary utilities. For a Mars expedition, the available electric power is severely limited. The design of both the processes and the equipment must be considered. With these constraints in mind, only the most efficient designs will be viable. Cryogenics, in situ resource utilization, miniaturization, launchability, and power/process efficiencies are only a few of the areas that chemical engineers provide support and expertise for the exploration of space.

Lobmeyer, Dennis A.↗
DOE OSTI2025

Thermal and kinematic properties of ejecta in SN1987A revealed by XRISM

We present an analysis of high-resolution spectra from the shock-heated plasmas in SN 1987A, based on an observation using the Resolve instrument onboard the X-Ray Imaging and Spectroscopy Mission (XRISM). The 1.7–10 keV Resolve spectra are accurately represented by a single-component, plane-parallel shock plasma model, with a temperature of $2.84_{-0.08}^{+0.09}$ keV and an ionization parameter of $2.64_{-0.45}^{+0.58}$ × $10^{11}\,\,{\rm s\,\, cm}^{-3}$. The Resolve spectra are also well reproduced by the 3D magneto-hydrodynamic simulation presented by Orlando et al. (2020, A&A, 636, A22) suggesting substantial contribution from the ejecta. The metal abundances obtained with Resolve align with the Large Magellanic Cloud value, indicating that the X-rays in 2024 originate from “non-metal-rich” shock-heated ejecta and the reverse shock has not reached the inner metal-rich region of ejecta. Doppler widths of the atomic lines from Si, S, and Fe correspond to velocities of 1500–1700 km s$^{-1}$, where the thermal broadening effects in this non-metal-rich plasma are negligible. Therefore, the line broadening seen in Resolve spectra is determined by the large bulk motion of ejecta. For reference, we determined a $90\%$ upper limit on non-thermal emission from a pulsar wind nebula at $4.3 \times 10^{-13}$ erg cm$^{-2}$ s$^{-1}$ in the 2–10 keV range, aligning with NuSTAR findings by Greco et al. (2022, ApJ, 931, 132). Additionally, we searched for the $^{44}$Sc K line feature and found a $1\sigma$ upper limit of $1.0 \times 10^{-6}$ photons cm$^{-2}$ s$^{-1}$, which translates to an initial $^{44}$Ti mass of approximately $2 \times 10^{-4}\, M_{\odot }$, consistent with previous X-ray to soft gamma-ray observations (Boggs et al. 2015, Science, 348, 670; Grebenev et al. 2012, Nature, 490, 373; Leising 2006, ApJ, 651, 1019).

ISM: supernova remnants↗
DOE OSTI2021

Midwest US Croplands Determine Model Divergence in North American Carbon Fluxes

Large uncertainties in North American terrestrial carbon fluxes hinder regional climate projections. Terrestrial biosphere models (TBMs), the essential tools for understanding continental-scale carbon cycle, diverge on whether temperate forests or croplands dominate carbon uptake in North America. Evidence from novel photosynthetic proxies, such as those based on chlorophyll fluorescence, has cast doubt on the “weak cropland, strong forest” carbon uptake patterns simulated by most TBMs. However, no systematic evaluation of TBMs has yet been attempted to pin down space-time patterns that are most consistent with regional CO 2 observational constraints. Here, we leverage atmospheric CO 2 observations and satellite-observed photosynthetic proxies to understand emergent space-time patterns in North American carbon fluxes from a large suite of TBMs and data-driven models. To do so, we evaluate how well the atmospheric signals resulting from carbon flux estimates reproduce the space-time variability in atmospheric CO 2 , as is observed by a network of continuous-monitoring towers over North America. Models with gross or net carbon fluxes that are consistent with the observed CO 2 variability share a salient feature of growing-season carbon uptake in Midwest US croplands. Conversely, the remaining models place most growing-season uptake in boreal or temperate forests. Differences in model explanatory power depend mainly on the simulated annual cycles of cropland uptake—especially, the timing of peak uptake—rather than the distribution of annual mean fluxes across biomes. Our results suggest that improved model representation of cropland phenology is crucial to robust, policy-relevant estimation of North American carbon exchange.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Rainfall manipulation experiments as simulated by terrestrial biosphere models: Where do we stand?

Changes in rainfall amounts and patterns have been observed and are expected to continue in the near future with potentially significant ecological and societal consequences. Modelling vegetation responses to changes in rainfall is thus crucial to project water and carbon cycles in the future. In this study, we present the results of a new model-data intercomparison project, where we tested the ability of 10 terrestrial biosphere models to reproduce the observed sensitivity of ecosystem productivity to rainfall changes at 10 sites across the globe, in nine of which, rainfall exclusion and/or irrigation experiments had been performed. The key results are as follows: (a) Inter-model variation is generally large and model agreement varies with timescales. In severely water-limited sites, models only agree on the interannual variability of evapotranspiration and to a smaller extent on gross primary productivity. In more mesic sites, model agreement for both water and carbon fluxes is typically higher on fine (daily–monthly) timescales and reduces on longer (seasonal–annual) scales. (b) Models on average overestimate the relationship between ecosystem productivity and mean rainfall amounts across sites (in space) and have a low capacity in reproducing the temporal (interannual) sensitivity of vegetation productivity to annual rainfall at a given site, even though observation uncertainty is comparable to inter-model variability. (c) Most models reproduced the sign of the observed patterns in productivity changes in rainfall manipulation experiments but had a low capacity in reproducing the observed magnitude of productivity changes. Models better reproduced the observed productivity responses due to rainfall exclusion than addition. (d) All models attribute ecosystem productivity changes to the intensity of vegetation stress and peak leaf area, whereas the impact of the change in growing season length is negligible. Finally, the relative contribution of the peak leaf area and vegetation stress intensity was highly variable among models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Can we achieve atmospheric chemical environments in the laboratory? An integrated model-measurement approach to chamber SOA studies

Secondary organic aerosol (SOA), atmospheric particulate matter formed from low-volatility products of volatile organic compound (VOC) oxidation, affects both air quality and climate. Current 3D models, however, cannot reproduce the observed variability in atmospheric organic aerosol. Because many SOA model descriptions are derived from environmental chamber experiments, our ability to represent atmospheric conditions in chambers directly affects our ability to assess the air quality and climate impacts of SOA. Here, we develop an approach that leverages global modeling and detailed mechanisms to design chamber experiments that mimic the atmospheric chemistry of organic peroxy radicals (RO 2 ), a key intermediate in VOC oxidation. Drawing on decades of laboratory experiments, we develop a framework for quantitatively describing RO 2 chemistry and show that no previous experimental approaches to studying SOA formation have accessed the relevant atmospheric RO 2 fate distribution. We show proof-of-concept experiments that demonstrate how SOA experiments can access a range of atmospheric chemical environments and propose several directions for future studies.

Science & Technology - Other Topics↗
DOE OSTI2021

A crystal plasticity model with an atomistically informed description of grain boundary sliding for improved predictions of deformation fields

Crystal plasticity (CP) is a powerful meso-scale technique for deformation modeling in polycrystalline materials. Crystal plasticity models typically do not include an explicit description of grain boundary (GB) sliding, which could lead to inaccurate predictions of strain distributions, especially in the vicinity of GBs. In the present study, a CP model is developed that includes GB sliding as an additional deformation mechanism and models the interaction between slip and GB sliding. Atomistic simulations are used to formulate the constitutive model for GB sliding and its interaction with incident slip. The deformation fields and structure of the GBs are obtained directly from experimental characterization and are faithfully reproduced as bicrystal systems for molecular dynamics simulations. The GBs are modeled as random-type (non-coincidence, asymmetrical) boundaries as observed from experimental data. The underlying atomistic-scale mechanism for pure sliding of random GBs was found to be analogous to fluid-flow. The interaction between incident slip and a sliding GB caused local increases in stress concentration, which further led to a momentary increase in the local sliding rate. The displacement profiles at the sliding GBs computed from the CP model with sliding accommodation are 38% more accurate than the baseline CP model as quantified by the mean squared error between the simulations and experiment. Here these results help improve the sophistication and accuracy of deformation modeling by including physics-based descriptions for GB sliding and its interaction with slip, eventually leading to more reliable predictions of micromechanical quantities.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A critical look at the prediction of the temperature field around a laser-induced melt pool on metallic substrates

The study of microstructure evolution in additive manufacturing of metals would be aided by knowing the thermal history. Since temperature measurements beneath the surface are difficult, estimates are obtained from computational thermo-mechanical models calibrated against traces left in the sample revealed after etching, such as the trace of the melt pool boundary. Here we examine the question of how reliable thermal histories computed from a model that reproduces the melt pool trace are. To this end, we perform experiments in which one of two different laser beams moves with constant velocity and power over a substrate of 17-4PH SS or Ti-6Al-4V, with low enough power to avoid generating a keyhole. We find that thermal histories appear to be reliably computed provided that (a) the power density distribution of the laser beam over the substrate is well characterized, and (b) convective heat transport effects are accounted for. Poor control of the laser beam leads to potentially multiple three-dimensional melt pool shapes compatible with the melt pool trace, and therefore to multiple potential thermal histories. Ignoring convective effects leads to results that are inconsistent with experiments, even for the mild melt pools here.

42 ENGINEERING↗
DOE OSTI2026

Modeling the behavior of concentrated aqueous HNO 3 using machine learning interatomic potentials

We develop two multi-defect machine learning interatomic potentials (MLIPs) trained at the BLYP-D2 and PBE-D3 density functional theories using the DeepMD-kit, allowing for the investigation of structural and thermodynamic properties of nitric acid over a wide range of concentrations via molecular dynamics (MD) simulations. We directly compute the degree of dissociation, α, and pK a from MD simulations, revealing that HNO 3 behaves as a weaker acid at higher concentrations, noting that our standard-state pK a value is in excellent agreement with the experimental one. In general, good agreement is observed with experimental results such as α and density outside the training dataset, with only modest deviations at low-to-medium concentrations. We benchmark our custom multi-defect DeepMD MLIPs against foundational models MACE-MP0 and MACE-OFF23. The foundation models capture some aspects of HNO 3 /NO 3 − solvation in concentrated nitric acid but show noticeable density errors and miss subtle structural features relevant to spectroscopy, whereas the bespoke DeepMD MLIPs yield more compact solvation shells, reproduce density-concentration trends, and run ∼12–15× faster than MACE-MP0. Although classical FFs are still more efficient and match experimental densities better, they lack chemical reactivity and thus cannot predict α or pK a , underscoring the need for system-specific reactive MLIPs beyond universal MLIPs.

Dinpajooh, Mohammadhasan [Pacific Northwest Nation↗
DOE OSTI2025

Multistage nucleation pathway in LiF molten salt mirrors the crystal–melt interface structure

Despite over a century of studies, fundamental questions remain about the processes governing crystal nucleation from melts or solutions. Research over the past three decades has presented mounting evidence for kinetic pathways of crystal nucleation that are more complex than envisioned by the simplest forms of classical theory. Such observations have been presented for colloidal and elemental systems with covalent and metallic bonding. Despite the technological and geochemical importance of molten salts, similar studies for these ionically bonded systems are currently lacking. Here we develop a machine learning interatomic potential for a model ionic system: LiF. The potential features quantum-level accuracy for both liquid and multiple solid polymorphs over wide temperature and pressure ranges and accurately reproduces experimentally measured properties. Thanks to the efficiency of the potential, which enables microsecond-scale molecular dynamics simulations, induction times for nucleation of LiF solids from their melts are computed over a range of undercoolings. With the aid of a set of robust local order parameters established here, the simulations reveal that homogeneous crystal nucleation in undercooled melts preferentially initiates from liquid regions showing slow dynamics and high bond orientational order simultaneously, and the second-shell order of both precritical nuclei and the surface of postcritical nuclei is dominated by hexagonal close packing and body-centered cubic local structure, even though the nucleus core is dominated by face-centered cubic structure corresponding to the stable rocksalt crystal structure. Finally, we establish a connection between the crystallization pathway and the equilibrium crystal-melt interface structure.

Applied Physical Sciences↗
DOE OSTI2022

Assessment of Bacterial Inoculant Delivery Methods for Cereal Crops

Despite growing evidence that plant growth-promoting bacteria can be used to improve crop vigor, a comparison of the different methods of delivery to determine which is optimal has not been published. An optimal inoculation method ensures that the inoculant colonizes the host plant so that its potential for plant growth-promotion is fully evaluated. The objective of this study was to compare the efficacy of three seed coating methods, seedling priming, and soil drench for delivering three bacterial inoculants to the sorghum rhizosphere and root endosphere. The methods were compared across multiple time points under axenic conditions and colonization efficiency was determined by quantitative polymerase chain reaction (qPCR). Two seed coating methods were also assessed in the field to test the reproducibility of the greenhouse results under non-sterile conditions. In the greenhouse seed coating methods were more successful in delivering the Gram-positive inoculant (Terrabacter sp.) while better colonization from the Gram-negative bacteria (Chitinophaga pinensis and Caulobacter rhizosphaerae) was observed with seedling priming and soil drench. This suggested that Gram-positive bacteria may be more suitable for the seed coating methods possibly because of their thick peptidoglycan cell wall. We also demonstrated that prolonged seed coating for 12 h could effectively enhance the colonization of C. pinensis, an endophytic bacterium, but not the rhizosphere colonizing C. rhizosphaerae. In the field only a small amount of inoculant was detected in the rhizosphere. This comparison demonstrates the importance of using the appropriate inoculation method for testing different types of bacteria for their plant growth-promotion potential.

Terrabacter↗
DOE OSTI2022

MPAS-Seaice (v1.0.0): sea-ice dynamics on unstructured Voronoi meshes

Abstract. We present MPAS-Seaice, a sea-ice model which uses the Model for Prediction Across Scales (MPAS) framework and spherical centroidal Voronoi tessellation (SCVT) unstructured meshes. As well as SCVT meshes, MPAS-Seaice can run on the traditional quadrilateral grids used by sea-ice models such as CICE. The MPAS-Seaice velocity solver uses the elastic–viscous–plastic (EVP) rheology and the variational discretization of the internal stress divergence operator used by CICE, but adapted for the polygonal cells of MPAS meshes, or alternatively an integral (“finite-volume”) formulation of the stress divergence operator. An incremental remapping advection scheme is used for mass and tracer transport. We validate these formulations with idealized test cases, both planar and on the sphere. The variational scheme displays lower errors than the finite-volume formulation for the strain rate operator but higher errors for the stress divergence operator. The variational stress divergence operator displays increased errors around the pentagonal cells of a quasi-uniform mesh, which is ameliorated with an alternate formulation for the operator. MPAS-Seaice shares the sophisticated column physics and biogeochemistry of CICE and when used with quadrilateral meshes can reproduce the results of CICE. We have used global simulations with realistic forcing to validate MPAS-Seaice against similar simulations with CICE and against observations. We find very similar results compared to CICE, with differences explained by minor differences in implementation such as with interpolation between the primary and dual meshes at coastlines. We have assessed the computational performance of the model, which, because it is unstructured, runs with 70 % of the throughput of CICE for a comparison quadrilateral simulation. The SCVT meshes used by MPAS-Seaice allow removal of equatorial model cells and flexibility in domain decomposition, improving model performance. MPAS-Seaice is the current sea-ice component of the Energy Exascale Earth System Model (E3SM).

58 GEOSCIENCES↗