Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ice core”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Timing of iceberg scours and massive ice-rafting events in the subtropical North Atlantic

Abstract High resolution seafloor mapping shows extraordinary evidence that massive (>300 m thick) icebergs once drifted >5,000 km south along the eastern United States, with >700 iceberg scours now identified south of Cape Hatteras. Here we report on sediment cores collected from several buried scours that show multiple plow marks align with Heinrich Event 3 (H3), ~31,000 years ago. Numerical glacial iceberg simulations indicate that the transport of icebergs to these sites occurs during massive, but short-lived, periods of elevated meltwater discharge. Transport of icebergs to the subtropics, away from deep water formation sites, may explain why H3 was associated with only a modest increase in ice-rafting across the subpolar North Atlantic, and implies a complex relationship between freshwater forcing and climate change. Stratigraphy from subbottom data across the scour marks shows there are additional features that are both older and younger, and may align with other periods of elevated meltwater discharge.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Effects from Time Dependence of Ice Nucleus Activity for Contrasting Cloud Types

Abstract The role of time-dependent freezing of ice nucleating particles (INPs) is evaluated with the “Aerosol–Cloud” (AC) model in 1) deep convection observed over Oklahoma during the Midlatitude Continental Convective Cloud Experiment (MC3E), 2) orographic clouds observed over North California during the Atmospheric Radiation Measurement (ARM) Cloud Aerosol Precipitation Experiment (ACAPEX), and 3) supercooled, stratiform clouds over the United Kingdom, observed during the Aerosol Properties, Processes And Influences on the Earth’s climate (APPRAISE) campaign. AC uses the dynamical core of the WRF Model and has hybrid bin–bulk microphysics and a 3D mesoscale domain. AC is validated against coincident aircraft, ground-based, and satellite observations for all three cases. Filtered concentrations of ice (>0.1–0.2 mm) agree with those observed at all sampled levels. AC predicts the INP activity of various types of aerosol particles with an empirical parameterization (EP), which follows a singular approach (no time dependence). Here, the EP is modified to represent time-dependent INP activity by a purely empirical approach, using our published laboratory observations of time-dependent INP activity. In all simulated clouds, the inclusion of time dependence increases the predicted INP activity of mineral dust particles by 0.5–1 order of magnitude. However, there is little impact on the cloud glaciation because the total ice is mostly (80%–90%) from secondary ice production (SIP) at levels warmer than about −36°C. The Hallett–Mossop process and fragmentation in ice–ice collisions together initiate about 70% of the total ice, whereas fragmentation during both raindrop freezing and sublimation contributes <10%. Overall, total ice concentrations and SIP are unaffected by time-dependent INP activity. In the simulated APPRAISE case, the main causes of persistence of long-lived clouds and precipitation are predicted to be SIP in weak embedded convection and reactivation following recirculation of dust particles in supercooled layer cloud.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

System Study: High-Pressure Core Spray 1998-2020

This report presents an unreliability evaluation of the high-pressure core spray (HPCS) at eight U.S. commercial boiling water reactors. Demand, run hour, and failure data from 1998–2020 for selected components were obtained from the Institute of Nuclear Power Operations (INPO) Industry Reporting and Information System (IRIS), formerly the INPO Consolidated Events Database (ICES). The unreliability results are trended for the most recent 10 year period while yearly estimates for system unreliability are provided for the entire active period. No statistically significant increasing or decreasing trends were identified in the HPCS results.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2022

System Study: Reactor Core Isolation Cooling 1998-2020

This report presents an unreliability evaluation of the reactor core isolation cooling (RCIC) system at 31 U.S. commercial boiling water reactors. Demand, run hours, and failure data from calendar years 1998 through 2020 for selected components were obtained from the Institute of Nuclear Power Operations (INPO) Industry Reporting and Information System (IRIS), formerly the INPO Consolidated Events Database (ICES). The unreliability results are trended for the most recent 10-year period while yearly estimates for system unreliability are provided for the entire active period. No statistically significant increasing or decreasing trends were identified in the RCIC results.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2024

Mapping Permafrost Variability and Degradation Using Seismic Surface Waves, Electrical Resistivity, and Temperature Sensing: A Case Study in Arctic Alaska

Subsurface processes significantly influence surface dynamics in permafrost regions, necessitating utilizing diverse geophysical methods to reliably constrain permafrost characteristics. This research uses multiple geophysical techniques to explore the spatial variability of permafrost in undisturbed tundra and its degradation in disturbed tundra in Utqiagvik, Alaska. Here, we integrate multiple quantitative techniques, including multichannel analysis of surface waves (MASW), electrical resistivity tomography (ERT), and ground temperature sensing, to study heterogeneity in permafrost’s geophysical characteristics. MASW results reveal active layer shear wave velocities (V s ) between 240 and 370 m/s, and permafrost V s between 450 and 1,700 m/s, typically showing a low-high-low velocity pattern. Additionally, we find an inverse relationship between in situ V s and ground temperature measurements. The V s profiles along with electrical resistivity profiles reveal cryostructures such as cryopeg and ice-rich zones in the permafrost layer. The integrated results of MASW and ERT provide valuable information for characterizing permafrost heterogeneity and cryostructure. Corroboration of these geophysical observations with permafrost core samples’ stratigraphies and salinity measurements further validates these findings. This combination of geophysical and temperature sensing methods along with permafrost core sampling confirms a robust approach for assessing permafrost’s spatial variability in coastal environments. Our results also indicate that civil infrastructure systems such as gravel roads and pile foundations affect permafrost by thickening the active layer, lowering the V s , and reducing heterogeneity. We show how the resulting V s profiles can be used to estimate key parameters for designing buildings in permafrost regions and maintaining existing infrastructure in polar regions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Observations and Machine-Learned Models of Near-Surface Permafrost along the Koyukuk River, Alaska, USA

This dataset contains GeoTIFs (raster) and GeoPackages (vector) that map observations of near-surface permafrost and not-permafrost from a field campaign conducted near the village of Huslia, AK along the Koyukuk River and its floodplain in July 2018. These data were collected as part of a campaign to understand if and how permafrost impacts riverbank erosion. This problem cannot be assessed without knowing where permafrost exists. Permafrost was observed via frost probing (to a maximum depth of one meter), coring (to a maximum depth of two meters) and bank/bar excavations. An additional boat survey was performed wherein expert (Joel Rowland) judgment assessed the presence or absence of distinctive permafrost features (e.g., overhanging tundra mats, thermoerosional niching, ice wedges, active drainage of ice melt from soils). This dataset also contains the input features and results of two machine learning models (random forest and convolutional neural network) that extrapolate the observations to the full floodplain that may be useful for building, testing, or validating other machine-learned permafrost models. Permafrost data are provided as georasters of the same shape and geovectors (polylines/polygons) and are all projected into EPSG:32605. All data can be visualized with a GIS (QGIS, ArcGIS, etc.).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Enrichment of Phosphates, Lead, and Mixed Soil–Organic Particles in INPs at the Southern Great Plains Site

Ice nucleating particles (INPs) are rare particles that initiate primary ice formation, a critical step required for subsequent important cloud microphysical processes that ultimately govern cloud phase and cloud radiative properties. Laboratory studies have found that organic-rich dusts, such as those found in soils, are more efficient INPs compared to mineral dust. However, the atmospheric relevance of these organic-rich dusts are not well understood, particularly in regions with significant agricultural activity. The Agricultural Ice nuclei at the Southern Great Plains field campaign (AGINSGP) was conducted in rural Oklahoma to investigate how soil dusts contribute to INP populations in the Great Plains. We present chemical characterization of ambient and ice crystal residual particles from a single day of sampling, using single particle mass spectrometry (SPMS) and scanning microscopy. Ambient particles were primarily carbonaceous or secondary aerosol, while the fraction of dust particles was higher in the residual particles. We also observed an unusual particle type consisting of a carbonaceous core mixed with dust fragments on the surface, which was found in higher proportion in residuals. Dust particles measured during residual sampling contained greater proportions of phosphate $(^{63}\text{PO}^-_2$ and $^{79}\text{PO}^-_3)$ and lead ( 206 Pb + ). Strong sulfate signals were not seen in the residual dust particles measured by the SPMS, while nitrate was slightly depleted relative to ambient dust. This study shows that organic-rich soils may be important contributors to the ambient INP population in agricultural regions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Classification of Cloud Particle Imagery from Aircraft Platforms Using Convolutional Neural Networks

Abstract A vast amount of ice crystal imagery exists from a variety of field campaign initiatives that can be utilized for cloud microphysical research. Here, nine convolutional neural networks are used to classify particles into nine regimes on over 10 million images from the Cloud Particle Imager probe, including liquid and frozen states and particles with evidence of riming. A transfer learning approach proves that the Visual Geometry Group (VGG-16) network best classifies imagery with respect to multiple performance metrics. Classification accuracies on a validation dataset reach 97% and surpass traditional automated classification. Furthermore, after initial model training and preprocessing, 10 000 images can be classified in approximately 35 s using 20 central processing unit cores and two graphics processing units, which reaches real-time classification capabilities. Statistical analysis of the classified images indicates that a large portion (57%) of the dataset is unusable, meaning the images are too blurry or represent indistinguishable small fragments. In addition, 19% of the dataset is classified as liquid drops. After removal of fragments, blurry images, and cloud drops, 38% of the remaining ice particles are largely intersecting the image border (≥10% cutoff) and therefore are considered unusable because of the inability to properly classify and dimensionalize. After this filtering, an unprecedented database of 1 560 364 images across all campaigns is available for parameter extraction and bulk statistics on specific particle types in a wide variety of storm systems, which can act to improve the current state of microphysical parameterizations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

FROSch Preconditioners for Land Ice Simulations of Greenland and Antarctica

Numerical simulations of Greenland and Antarctic ice sheets involve the solution of large-scale highly nonlinear systems of equations on complex shallow geometries. This work is concerned with the construction of Schwarz preconditioners for the solution of the associated tangent problems, which are challenging for solvers mainly because of the strong anisotropy of the meshes and wildly changing boundary conditions that can lead to poorly constrained problems on large portions of the domain. In this study, two-level generalized Dryja--Smith--Widlund (GDSW)--type Schwarz preconditioners are applied to different land ice problems, i.e., a velocity problem, a temperature problem, as well as the coupling of the former two problems. We employ the message passing interface (MPI)--parallel implementation of multilevel Schwarz preconditioners provided by the package FROSch (fast and robust Schwarz) from the Trilinos library. The strength of the proposed preconditioner is that it yields out-of-the-box scalable and robust preconditioners for the single physics problems. To the best of our knowledge, this is the first time two-level Schwarz preconditioners have been applied to the ice sheet problem and a scalable preconditioner has been used for the coupled problem. The preconditioner for the coupled problem differs from previous monolithic GDSW preconditioners in the sense that decoupled extension operators are used to compute the values in the interior of the subdomains. Additionally, several approaches for improving the performance, such as reuse strategies and shared memory OpenMP parallelization, are explored. In our numerical study we target both uniform meshes of varying resolution for the Antarctic ice sheet as well as nonuniform meshes for the Greenland ice sheet. We present several weak and strong scaling studies confirming the robustness of the approach and the parallel scalability of the FROSch implementation. Among the highlights of the numerical results are a weak scaling study for up to 32,K processor cores (8,K MPI ranks and 4 OpenMP threads) and 566,M degrees of freedom for the velocity problem as well as a strong scaling study for up to 4,K processor cores (and MPI ranks) and 68,M degrees of freedom for the coupled problem.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Spherical compression of an applied magnetic field in inertial confinement fusion

Applying an external magnetic field to laser-driven inertial confinement fusion implosions is a promising approach for enhancing fusion yield. The field is compressed with the plasma, producing a magnetized hotspot that anisotropically suppresses thermal losses and traps alpha particles, making performance sensitive to the compressed field orientation. Here, we derive a simple, readily applicable analytic model that enables rapid evaluation of the compressed field topology and show that ablation into the hotspot amplifies the central field, while the ablated ice near the hotspot edge develops a decaying, radially bent field, with a discontinuity in the field direction. The radially bent field renders thermal insulation at the hotspot edge negligible and largely independent of the applied field strength, whereas insulation in the hotspot core still depends strongly on the applied field. Applying the model to non-axial initial field configurations, we find that an initially applied mirror field provides the greatest suppression, followed by the standard axial field.

Physics - Plasma physics↗
DOE OSTI2021

FROSch Preconditioners for Land Ice Simulations of Greenland and Antarctica

Numerical simulations of Greenland and Antarctic ice sheets involve the solution of large-scale highly nonlinear systems of equations on complex shallow geometries. This work is concerned with the construction of Schwarz preconditioners for the solution of the associated tangent problems, which are challenging for solvers mainly because of the strong anisotropy of the meshes and wildly changing boundary conditions that can lead to poorly constrained problems on large portions of the domain. Here, two-level GDSW (Generalized Dryja–Smith–Widlund) type Schwarz preconditioners are applied to different land ice problems, i.e., a velocity problem, a temperature problem, as well as the coupling of the former two problems. We employ the MPI-parallel implementation of multi-level Schwarz preconditioners provided by the package FROSch (Fast and Robust Schwarz)from the Trilinos library. The strength of the proposed preconditioner is that it yields out-of-the-box scalable and robust preconditioners for the single physics problems. To our knowledge, this is the first time two-level Schwarz preconditioners are applied to the ice sheet problem and a scalable preconditioner has been used for the coupled problem. The pre-conditioner for the coupled problem differs from previous monolithic GDSW preconditioners in the sense that decoupled extension operators are used to compute the values in the interior of the sub-domains. Several approaches for improving the performance, such as reuse strategies and shared memory OpenMP parallelization, are explored as well. In our numerical study we target both uniform meshes of varying resolution for the Antarctic ice sheet as well as non uniform meshes for the Greenland ice sheet are considered. We present several weak and strong scaling studies confirming the robustness of the approach and the parallel scalability of the FROSch implementation. Among the highlights of the numerical results are a weak scaling study for up to 32 K processor cores (8 K MPI-ranks and 4 OpenMP threads) and 566 M degrees of freedom for the velocity problem as well as a strong scaling study for up to 4 K processor cores (and MPI-ranks) and 68 M degrees of freedom for the coupled problem.

58 GEOSCIENCES↗
DOE OSTI2021

Superionic iron oxide–hydroxide in Earth’s deep mantle

Water ice becomes a superionic phase under the high pressure and temperature conditions of deep planetary interiors of ice planets such as Neptune and Uranus, which affects interior structures and generates magnetic fields. The solid Earth, however, contains only hydrous minerals with a negligible amount of ice. Here we combine high pressure and temperature electrical conductivity experiments, Raman spectroscopy and first-principles simulations to investigate the state of hydrogen in the pyrite-type FeO 2 H x (x ≤ 1), which is a potential H-bearing phase near the core–mantle boundary. We find that when the pressure increases beyond 73 GPa at room temperature, symmetric hydroxyl bonds are softened and the H + (or proton) becomes diffusive within the vicinity of its crystallographic site. Increasing temperature under pressure, the diffusivity of hydrogen is extended beyond the individual unit cell to cover the entire solid, and the electrical conductivity soars, indicating a transition to the superionic state, which is characterized by freely moving protons and a solid FeO 2 lattice. Furthermore, the highly diffusive hydrogen provides fresh transport mechanisms for charge and mass, which dictate the geophysical behaviours of electrical conductivity and magnetism, as well as geochemical processes of redox, hydrogen circulation and hydrogen isotopic mixing in Earth’s deep mantle.

58 GEOSCIENCES↗
DOE OSTI2022

Miscibility of rock and ice in the interiors of water worlds

Abstract Super-Earths and sub-Neptunes are the most common planet types in our galaxy. A subset of these planets is predicted to be water worlds, bodies that are rich in water and poor in hydrogen gas. The interior structures of water worlds have been assumed to consist of water surrounding a rocky mantle and iron core. In small planets, water and rock form distinct layers with limited incorporation of water into silicate phases, but these materials may interact differently during the growth and evolution of water worlds due to greater interior pressures and temperatures. Here, we use density functional molecular dynamics (DFT-MD) simulations to study the miscibility and interactions of enstatite (MgSiO 3 ), a major end-member silicate phase, and water (H 2 O) at extreme conditions in water world interiors. We explore pressures ranging from 30 to 120 GPa and temperatures from 500 to 8000 K. Our results demonstrate that enstatite and water are miscible in all proportions if the temperature exceeds the melting point of MgSiO 3 . Furthermore, we performed smoothed particle hydrodynamics simulations to demonstrate that the conditions necessary for rock-water miscibility are reached during giant impacts between water-rich bodies of 0.7–4.7 Earth masses. Our simulations lead to water worlds that include a mixed layer of rock and water.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Big Sky Regional Carbon Sequestration Partnership (Phase III Final Scientific/Technical Report)

The Big Sky Carbon Sequestration Partnership (BSCSP) pursued a Phase III demonstration project at Kevin Dome in north central Montana. Kevin Dome covers approximately 700 square miles and is a naturally occurring CO 2 reservoir that is flanked by oil and gas fields. The carbon dioxide (CO 2 ) is in the upper Devonian Duperow (carbonate) formation and does not reach the spill point of the dome; therefore, the dome has potential as a CO 2 sequestration reservoir, a CO 2 supply, or as both if anthropogenic sources and enhanced oil recovery (EOR) operations are associated with the dome. Kevin Dome could potentially act as a buffer to continue accepting anthropogenic CO 2 when EOR flooding operations are interrupted or completed. The project objective was to produce one million tonnes of CO 2 from the gas cap of Kevin Dome, pipe it laterally, inject, and re-store it in the brine leg of the same formation to test the hub / buffer storage concept. This was to be accomplished by drilling up to five production wells, building a short pipeline and compression facilities, and drilling an injection well and several monitoring wells. BSCSP commenced outreach and site characterization activities including acquiring baseline data for near-surface insurance monitoring, acquiring 3-dimensional, 9-component surface seismic over the project area, drilling two test wells (one in the production area and one in the injection area), coring key intervals, and performing comprehensive logging. Well tests of those wells revealed two barriers to the project. The production (Danielson 33-17) well was perforated in multiple zones but failed to produce any significant CO 2 . This was despite being drilled in the near vicinity of a historic well that had produced 3700 MCF per day in a drill stem test. Modeling indicated that this was likely due to a phase change during production causing a temperature drop resulting in hydrate and/or water ice formation that clogged the formation. Tests of the injection zone (Wallewein 22-1) well indicated total dissolved solids (TDS) slightly below the EPA required 10,000 parts per million (ppm) for a Class VI underground injection control permit. While the project was initiated before Class VI rules were promulgated, and this was an experimental project (seemingly qualified for a Class V permit), the Environmental Protection Agency (EPA) indicated that injection would require a Class VI permit. The low salinity result was unexpected as contours plotted based on regional formation water quality data indicated an expected TDS above 20,000 ppm, and wells between the recharge zone and the Wallewein well tested above 10,000 ppm. Faced with the inability to obtain an injection permit, the demonstration project could not proceed. However, the project had generated valuable samples and data on a large natural analog including 32 sq. mi. of 3-D, 9-C seismic, 430 ft. of carbonate core covering seven different depositional environments taken from areas with, and without the presence of CO 2 , 30 ft. of core of two caprocks, a tight carbonate and an anhydrite, a full set of modern logs on both wells, and well tests. DOE decided to re-scope the project around completing studies utilizing this data. This report covers both the initial scope and the re-scope (Task / Section numbers preceded with an R). While the report covers a wide range of project activities, highlights of this work include: Development of a geostatic model using neural nets to match well logs to facies and using multi-waveform seismic to inform reservoir heterogeneity; Unique mechanical testing of permeability – stress relationship in two caprock materials; Development of full waveform inversion to generate a high resolution velocity model; Model development for dual permeability (fracture and matrix) systems to better account for matrix-matrix interactions; Joint seismic wave inversion (including the first quadr-joint inversion) exhibiting better imaging of a challenging reservoir zone in stiff rock; Core flow and core flood results on a reactive carbonate; and Innovative laboratory measurements of seismic response of fractured core as a function of fluid fill.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Free-spray characteristics and spray-wall interactions of methanol on a gasoline direct injector under flash-boiling and non-flash-boiling conditions

Methanol is considered a promising alternative fuel for internal combustion engines (ICEs) due to its high-octane number, fast laminar flame speed, and elevated latent heat of vaporization, all of which support higher compression ratios and improved thermal efficiency. However, its substantial latent heat of vaporization also poses cold-start challenges, such as misfire and fuel film deposition. This study aims to investigate methanol spray morphology and spray-wall interaction using the Spray M injector from the Engine Combustion Network within a constant-pressure flow vessel. A recently developed unified numerical framework capable of modeling both flash and non-flash boiling sprays is validated against experimental liquid volume fraction data acquired via 3-D computed tomography. Here, the results reveal that flash boiling significantly alters the spray morphology, leading to smaller droplets and spray collapse due to enhanced air-entrainment-induced turbulence. Quantitative agreement between experiments and simulations confirms this behavior. Coupled 0-D equilibrium and 3-D computational fluid dynamics analyses show that flash boiling accelerates evaporation and reduces fuel residence time, while non-flash conditions maintain a persistent liquid core more susceptible to wall wetting. Wall temperature diagnostics reveal that spray collapse alters heat transfer patterns by shifting cooling effects. Mixture fraction analysis indicates that evaporation is primarily governed by shear-layer turbulence, though deviations from adiabatic equilibrium mixing emerge under low-turbulence conditions. Finally, increasing fuel, ambient, and wall temperatures reduces wall wetting and film thickness, mitigating cold-start risks. These findings enhance the understanding of methanol sprays’ behavior and support its adoption as a viable, alternative fuel for ICEs.

Engine Combustion Network↗
DOE OSTI2021

High yields in direct-drive inertial confinement fusion using thin-ice DT liner targets

In this work, we report simulation results that suggest record neutron yields above 3 x 10 14 can be achieved in laser-driven fusion on OMEGA by imploding large-outer-diameter targets (~1010μm) containing ice layers as thin as 34μm. Historically, the performance of thin and fast capsules is thought to be highly degraded by hydrodynamic instabilities. We present pulse designs that push implosion velocities above 650 km/s while maintaining enough ablative stabilization and low enough in-flight aspect ratios to survive the acceleration phase. We show these targets also possess improved stability during deceleration stemming from augmented core temperatures and low fuel convergence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Earth's water reservoirs in a changing climate

Progress towards achieving a quantitative understanding of the exchanges of water between Earth's main water reservoirs is reviewed with emphasis on advances accrued from the latest advances in Earth Observation from space. These exchanges of water between the reservoirs are a result of processes that are at the core of important physical Earth-system feedbacks, which fundamentally control the response of Earth's climate to the greenhouse gas forcing it is now experiencing, and are therefore vital to understanding the future evolution of Earth's climate. The changing nature of global mean sea level (GMSL) is the context for discussion of these exchanges. Different sources of satellite observations that are used to quantify ice mass loss and water storage over continents, how water can be tracked to its source using water isotope information and how the waters in different reservoirs influence the fluxes of water between reservoirs are described. The profound influence of Earth's hydrological cycle, including human influences on it, on the rate of GMSL rise is emphasized. The many intricate ways water cycle processes influence water exchanges between reservoirs and thus sea-level rise, including disproportionate influences by the tiniest water reservoirs, are emphasized.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Refactoring the elastic–viscous–plastic solver from the sea ice model CICE v6.5.1 for improved performance

This study focuses on the performance of the elastic–viscous–plastic (EVP) dynamical solver within the sea ice model, CICE v6.5.1. The study has been conducted in two steps. First, the standard EVP solver was extracted from CICE for experiments with refactored versions, which are used for performance testing. Second, one refactored version was integrated and tested in the full CICE model to demonstrate that the new algorithms do not significantly impact the physical results. The study reveals two dominant bottlenecks, namely (1) the number of Message Parsing Interface (MPI) and Open Multi-Processing (OpenMP) synchronization points required for halo exchanges during each time step combined with the irregular domain of active sea ice points and (2) the lack of single-instruction, multiple-data (SIMD) code generation. The standard EVP solver has been refactored based on two generic patterns. The first pattern exposes how general finite differences on masked multi-dimensional arrays can be expressed in order to produce significantly better code generation by changing the memory access pattern from random access to direct access. The second pattern takes an alternative approach to handle static grid properties. The measured single-core performance improvement is more than a factor of 5 compared to the standard implementation. The refactored implementation of strong scales on the Intel® Xeon® Scalable Processors series node until the available bandwidth of the node is used. For the Intel® Xeon® CPU Max series, there is sufficient bandwidth to allow the strong scaling to continue for all the cores on the node, resulting in a single-node improvement factor of 35 over the standard implementation. This study also demonstrates improved performance on GPU processors.

58 GEOSCIENCES↗