Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Core processes”

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,549 records · Page 86

Phylogenomic discovery and engineering of nitrogen fixation into the bioenergy woody crop poplar

Biological nitrogen fixation (BNF) is a key process enabling plants in specific lineages to convert atmospheric dinitrogen (N₂) into bioavailable ammonia through symbioses with diazotrophic microbes. Expanding this capability beyond native nitrogen-fixing clades into non-nodulating crops would reduce synthetic fertilizer use, lowering energy inputs and environmental impacts in agriculture. Supported by DOE Funding Award DE-SC0018247, the NitFix project advanced foundational knowledge required to engineer root-nodule symbioses in new host species. The team generated the most comprehensive phylogenomic analysis to date of all known nodulating lineages, resolving the evolutionary history of nitrogen-fixing symbiosis and identifying core gene suites retained across nodulating taxa. Through multimodal genomics, transcriptomics, and functional analyses in Medicago truncatula and related species, the project mapped regulatory networks underlying nodule organogenesis, bacterial infection, and nitrogen-fixation efficiency. Key discoveries include the identification of conserved signaling modules for rhizobial recognition, transcription factors controlling nodule differentiation, and metabolic pathways integrating fixed nitrogen into plant growth. The project also developed enabling tools—including optimized transformation pipelines, gene-editing workflows, and imaging-based phenotyping—to accelerate engineering efforts in emerging models. Together, these results refine the mechanistic framework of symbiotic nitrogen fixation and highlight transferable components essential for rewiring these traits into non-nodulating crops.

59 BASIC BIOLOGICAL SCIENCES↗

The Galactic Distribution of Phosphorus: A Survey of 163 Disk and Halo Stars*

Phosphorus (P) is a critical element for life on Earth, yet the cosmic production sites of P are relatively uncertain. To understand how P has evolved in the solar neighborhood, we measured abundances for 163 FGK stars over a range of –1.09 < [Fe/H] < 0.47 using observations from the Habitable-zone Planet Finder instrument on the Hobby–Eberly Telescope. Atmospheric parameters were calculated by fitting a combination of astrometry, photometry, and Fe I line equivalent widths. Phosphorus abundances were measured by matching synthetic spectra to a P I feature at 10529.52 Å. Our [P/Fe] ratios show that chemical evolution models generally underpredict P over the observed metallicity range. Additionally, we find that the [P/Fe] differs by ~0.1 dex between thin disk and thick disk stars that were identified with kinematics. The P abundances were compared with α-elements, iron-peak, odd-Z, and s-process elements, and we found that the evolution of P in the disk most strongly resembles that of the α-elements. We also find that molar P/C and N/C ratios for our sample match the scatter seen from other abundance studies. Finally, we measure a [P/Fe] = 0.09 ± 0.1 ratio in one low-α halo star and probable Gaia–Sausage–Enceladus member, an abundance ratio ~0.3–0.5 dex lower than the other Milky Way disk and halo stars at similar metallicities. Overall, we find that P is likely most significantly produced by massive stars in core-collapse supernovae, based on the largest P abundance survey to date.

79 ASTRONOMY AND ASTROPHYSICS↗

Direct numerical simulations of turbulent reacting flows with shock waves and stiff chemistry using many-core/GPU acceleration

Compressible reacting flows may display sharp spatial variation related to shocks, contact discontinuities or reactive zones embedded within relatively smooth regions. The presence of such phenomena emphasizes the relevance of shock-capturing schemes such as the weighted essentially non-oscillatory (WENO) scheme as an essential ingredient of the numerical solver. However, these schemes are complex and have more computational cost than the simple high-order compact or non-compact schemes. In this paper, we present the implementation of a seventh-order, minimally-dissipative mapped WENO (WENO7M) scheme in a newly developed direct numerical simulation (DNS) code called KAUST Adaptive Reactive Flows Solver (KARFS). In order to make efficient use of the computer resources and reduce the solution time, without compromising the resolution requirement, the WENO routines are accelerated via graphics processing unit (GPU) computation. The performance characteristics and scalability of the code are studied using different grid sizes and block decomposition. Furthermore, the performance portability of KARFS is demonstrated on a variety of architectures including NVIDIA Tesla P100 GPUs and NVIDIA Kepler K20X GPUs. In addition, the capability and potential of the newly implemented WENO7M scheme in KARFS to perform DNS of compressible flows is also demonstrated with model problems involving shocks, isotropic turbulence, detonations and flame propagation into a stratified mixture with complex chemical kinetics.

97 MATHEMATICS AND COMPUTING↗

Optical Evidence of Compositional Fractioning between Plasma‐Condensed and Melt Pool Matter

Control of multicomponent alloys during welding is challenging because it lacks a real‐time understanding of composition. The optical emissions of plasma formed during laser‐induced metal welding correlate with the composition of particles ejected from the melt pool. Plasma emissions observed in this study contain large iron, manganese, chrome, and copper signatures, which match the composition of emitted particles. Particles recovered closest to the melt pool exhibit a core–shell morphology that is composed of iron‐manganese‐chrome intermetallic cores within copper shells. Particles collected farther from the melt pool, do not share this core–shell morphology, though similar elemental compositions are observed. The correlation between plasma optical emissions and particle composition can be used to predict the composition of the melt pool, allowing for real‐time welding and sintering control.

36 MATERIALS SCIENCE↗

TRACER Perspectives on Gulf-Breeze and Bay-Breeze Circulations and Coastal Convection

Abstract This study explores gulf-breeze circulations (GBCs) and bay-breeze circulations (BBCs) in Houston–Galveston, investigating their characteristics, large-scale weather influences, and impacts on surface properties, boundary layer updrafts, and convective clouds. The results are derived from a combination of datasets, including satellite observations, ground-based measurements, and reanalysis datasets, using machine learning, changepoint detection method, and Lagrangian cell tracking. We find that anticyclonic synoptic patterns during the summer months (June–September) favor GBC/BBC formation and the associated convective cloud development, representing 74% of cases. The main Tracking Aerosol Convection Interactions Experiment (TRACER) site located close to the Galveston Bay is influenced by both GBC and BBC, with nearly half of the cases showing evident BBC features. The site experiences early frontal passages ranging from 1040 to 1630 local time (LT), with 1300 LT being the most frequent. These fronts are stronger than those observed at the ancillary site which is located further inland from the Galveston Bay, including larger changes in surface temperature, moisture, and wind speed. Furthermore, these fronts trigger boundary layer updrafts, likely promoting isolated convective precipitating cores that are short lived (average convective lifetime of 63 min) and slow moving (average propagation speed of 5 m s −1 ), primarily within 20–40 km from the coast.

54 ENVIRONMENTAL SCIENCES↗

Drop size distribution broadening mechanisms in a bin microphysics Eulerian model

In this study, processes that broaden drop size distributions (DSDs) in Eulerian models with two-moment bin microphysics are analyzed. Numerous tests are performed to isolate the effects of different physical mechanisms that broaden DSDs in two- and three-dimensional Weather Research and Forecasting model simulations of an idealized ice-free cumulus cloud. Additionally, sensitivity of these effects to modifying horizontal and vertical model grid spacings is also examined. As expected, collision-coalescence is a key process broadening the modeled DSDs. In-cloud droplet activation also contributes substantially to DSD broadening, whereas evaporation has only a minor effect and sedimentation has little effect. Cloud dilution (mixing of cloud-free and cloudy air) also broadens the DSDs considerably, whether or not it is accompanied by evaporation. This mechanism involves the reduction of droplet concentration from dilution along the cloud’s lateral edges, leading to locally high supersaturation and enhanced drop growth when this air is subsequently lifted in the updraft. DSD broadening ensues when the DSDs are mixed with those from the cloud core. Decreasing the horizontal and vertical model grid spacings from 100 m to 30 m has limited impact on the DSDs. However, when these physical broadening mechanisms (in-cloud activation, collision-coalescence, dilution, etc.) are turned off, there is a reduction of DSD width by up to ~20-50% when the vertical grid spacing is decreased from 100 m to 30 m, consistent with effects of artificial broadening from vertical numerical diffusion. Nonetheless, this artificial numerical broadening appears to be relatively unimportant overall for DSD broadening when physically based broadening mechanisms in the model are included for this cumulus case.

54 ENVIRONMENTAL SCIENCES↗

Non-invasive Pipe Pressure Monitoring for Safeguards

Our project explores a new solution to complement current safeguards toolbox by independently, and noninvasively, determining the pressure inside uranium hexafluoride (UF 6 )-carrying Gas Centrifuge Enrichment Plant (GCEP) pipes. Proof-of-principle, non-invasive pipe pressure monitoring tests were conducted through measurements of hoop and axial strain using an optical-interference dilatometry technique based on fiber Bragg gratings (FBGs) inscribed in the core of telecommunication-type 125 µm diameter optical fibers. The fibers were affixed to commercial stainless steel and aluminum tubes and pipes of the type expected to be found at GCEP’s, each measuring four inches in diameter. When the internal air pressure in the pipes changes their external diameter and length follow, and the pitch in the FBG shrinks and/or expands accordingly. In turn, these changes result in a shift of the light wavelength reflected by the FBGs. Several FBGs physically attached to the pipe, as well as control FBGs not attached, were interrogated during the process using state-of-the-art commercial interrogators that operate on swept-wavelength laser technology. Strains measured were contrasted against calculations based on the known elastic and mechanical properties of the pipes. Good agreement is found, and strategies are proposed for further improvement of the experimental resolution as well as future implementation of this novel approach for GCEP safeguards.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

High-Resolution Frequency-Domain Spectroscopic and Modeling Studies of Photosystem I (PSI), PSI Mutants and PSI Supercomplexes

Photosystem I (PSI) is one of the two main pigment–protein complexes where the primary steps of oxygenic photosynthesis take place. This review describes low-temperature frequency-domain experiments (absorption, emission, circular dichroism, resonant and non-resonant hole-burned spectra) and modeling efforts reported for PSI in recent years. In particular, we focus on the spectral hole-burning studies, which are not as common in photosynthesis research as the time-domain spectroscopies. Experimental and modeling data obtained for trimeric cyanobacterial Photosystem I (PSI 3 ), PSI 3 mutants, and PSI 3 –IsiA 18 supercomplexes are analyzed to provide a more comprehensive understanding of their excitonic structure and excitation energy transfer (EET) processes. Detailed information on the excitonic structure of photosynthetic complexes is essential to determine the structure–function relationship. We will focus on the so-called “red antenna states” of cyanobacterial PSI, as these states play an important role in photochemical processes and EET pathways. The high-resolution data and modeling studies presented here provide additional information on the energetics of the lowest energy states and their chlorophyll (Chl) compositions, as well as the EET pathways and how they are altered by mutations. We present evidence that the low-energy traps observed in PSI are excitonically coupled states with significant charge-transfer (CT) character. The analysis presented for various optical spectra of PSI 3 and PSI 3 -IsiA 18 supercomplexes allowed us to make inferences about EET from the IsiA 18 ring to the PSI 3 core and demonstrate that the number of entry points varies between sample preparations studied by different groups. In our most recent samples, there most likely are three entry points for EET from the IsiA 18 ring per the PSI core monomer, with two of these entry points likely being located next to each other. Therefore, there are nine entry points from the IsiA 18 ring to the PSI 3 trimer. We anticipate that the data discussed below will stimulate further research in this area, providing even more insight into the structure-based models of these important cyanobacterial photosystems.

59 BASIC BIOLOGICAL SCIENCES↗

Chemical Analysis of the Ultrafaint Dwarf Galaxy Grus II. Signature of High-mass Stellar Nucleosynthesis

We present a detailed abundance analysis of the three brightest member stars at the top of the giant branch of the ultrafaint dwarf (UFD) galaxy Grus II. All stars exhibit a higher than expected [Mg/Ca] ratio compared to metal-poor stars in other UFD galaxies and in the Milky Way (MW) halo. Nucleosynthesis in high-mass ( 20 M ⊙) core-collapse supernovae has been shown to create this signature. The abundances of this small sample (three) stars suggests the chemical enrichment of Grus II could have occurred through substantial high-mass stellar evolution, and is consistent with the framework of a top-heavy initial mass function. However, with only three stars it cannot be ruled out that the abundance pattern is the result of a stochastic chemical enrichment at early times in the galaxy. The most metal-rich of the three stars also possesses a small enhancement in rapid neutron-capture (r-process) elements. The abundance pattern of the r-process elements in this star matches the scaled r-process pattern of the solar system and r-process enhanced stars in other dwarf galaxies and in the MW halo, hinting at a common origin for these elements across a range of environments. All current proposed astrophysical sites of r-process element production are associated with high-mass stars, thus the possible top-heavy initial mass function of Grus II would increase the likelihood of any of these events occurring. The time delay between the α and r-process element enrichment of the galaxy favors a neutron star merger as the origin of the r-process elements in Grus II.

79 ASTRONOMY AND ASTROPHYSICS↗

Investigation of a bacteriochlorin-containing pentad array for panchromatic light-harvesting and charge separation

A new pentad array designed to exhibit panchromatic absorption and charge separation has been synthesized and characterized. The array is composed of a triad panchromatic absorber (a bis(perylene-monoimide)-porphyrin) to which are appended an electron acceptor (perylene-diimide) and an electron donor/hole acceptor (bacteriochlorin) in a crossbar arrangement. The motivation for incorporation of the bacteriochlorin versus a free-base or zinc chlorin utilized in prior constructs was to facilitate hole transfer to this terminal unit and thereby achieve a higher yield of charge separation across the array. The intense S 0 → S 1 (Q y ) band of the bacteriochlorin also enhances absorption in the near-infrared spectral region. Due to synthetic constraints, a phenylethyne linker was used to join the bacteriochlorin to the core porphyrin of the panchromatic triad rather than the diphenylethyne linker employed for the prior chlorin-containing pentads. Static and time-resolved photophysical studies reveal enhanced excited-state quenching for the pentad in benzonitrile and dimethyl sulfoxide compared to the prior chlorin-containing analogues. Success was only partial, however, as a long-lived charge separated state was not observed despite the improved energetics for the final ground-state hole/electron-shift reaction. The apparent reason is more facile competing charge-recombination due to the shorter bacteriochlorin – porphyrin linker that increases electronic coupling for this process. Furthermore, the studies highlight design criteria for balancing panchromatic absorption and long-lived charge separation in molecular architectures for solar-energy conversion.

14 SOLAR ENERGY↗

Atomic-level mechanisms of short-circuit diffusion in materials

This paper reviews the recent progress in understanding the atomic mechanisms of short-circuit diffusion along materials interfaces, such as grain and interphase boundaries, as well as lattice and interfacial dislocations/disconnections. Recent atomistic computer simulations have shown that short-circuit diffusion is dominated by collective atomic rearrangements in the form of strings and rings of mobile atoms. The process is dynamically heterogeneous in space and time and has many features in common with atomic dynamics in supercooled glass-forming liquids. We discuss examples of grain boundary, interphase boundary, and dislocation diffusion in metals and alloys, including the solute effect on the diffusion rates and mechanisms. Interphase boundaries are exemplified by Al–Si interfaces with diverse orientation relationships and atomic structures. The hierarchy of short-circuit diffusion paths in materials is reviewed by comparing the rates of grain boundary, interphase boundary, and dislocation diffusion. Future directions in the field of short-circuit diffusion in defect core regions are discussed.

36 MATERIALS SCIENCE↗

Digital Miniature Cathode Ray Magnetometer

In this study, we introduce the concept and construction of an innovative Digital Miniature Cathode Ray Magnetometer designed for the precise detection of magnetic fields. This device addresses several limitations inherent to magnetic probes such as D.C. offset, nonlinearity, temperature drift, sensor aging, and the need for frequent recalibration, while capable of operating in a wide range of magnetic fields. The core principle of this device involves the utilization of a charged particle beam as the sensitivity medium. The system leverages the interaction of an electron beam with a scintillator material, which then emits visible light that is captured by an imager. The emitted scintillation light is captured by a CMOS sensor. This sensor not only records the scintillation light but also accurately determines the position of the electron beam, providing invaluable spatial information crucial for magnetic field mapping. The key innovation lies in the combination of electron beam projection, CMOS imager scintillation-based detection, and digital image signal processing. By employing this synergy, the magnetometer achieves remarkable accuracy, sensitivity and dynamic range. The precise position registration enabled by the CMOS sensor further enhances the device’s utility in capturing complex magnetic field patterns, allowing for 2D field mapping. In this work, the optimization of the probe’s performance is tailored for applications related to the characterization of insertion devices in light sources, including undulators.

Turqueti, Marcos↗

Dynamic Distribution of High-Rate Data Processing from CERN to Remote HPC Data Centers

The prompt reconstruction of the data recorded from the Large Hadron Collider (LHC) detectors has always been addressed by dedicated resources at the CERN Tier-0. Such workloads come in spikes due to the nature of the operation of the accelerator and in special high load occasions experiments have commissioned methods to distribute (spill-over) a fraction of the load to sites outside CERN. The present work demonstrates a new way of supporting the Tier-0 environment by provisioning resources elastically for such spilled-over workflows onto the Piz Daint Supercomputer at CSCS. Furthermore, this is implemented using containers, tuning the existing batch scheduler and reinforcing the scratch file system, while still using standard Grid middleware. ATLAS, CMS and CSCS have jointly run selected prompt data reconstruction on up to several thousand cores on Piz Daint into a shared environment, thereby probing the viability of the CSCS high performance computer site as on demand extension of the CERN Tier-0, which could play a role in addressing the future LHC computing challenges for the high luminosity LHC.

97 MATHEMATICS AND COMPUTING↗

Do not forget the electrons: Extending moderately-sized nuclear networks for multidimensional hydrodynamic codes

Context.Nuclear networks are widely used coupled with hydrodynamical simulations of explosive scenarios to account for the change of nuclear species and energy generation rate due to nuclear reactions. In this way, there is a feedback mechanism between the hydrodynamical state and the nuclear processes. Unfortunately, the timescale of nuclear reactions is orders of magnitude smaller than the dynamical timescale that drives hydrodynamical simulations. Therefore, these nuclear networks are usually very small, reduced in most cases to a dozen elements, especially when simulations are carried out in more than one dimension. Aims.We present here an extended nuclear network, with 90 species, designed for being coupled with hydrodynamic simulations, which includes neutrons, protons, electrons, positrons, and the corresponding neutrino and anti-neutrino emission. This network is also coupled with temperature, making it extremely robust and, together with its size, unique of its kind. The inclusion of electron captures on free protons makes the network very appropriate for multidimensional studies of Type Ia supernova explosions, especially when the exploding object is a massive white dwarf. Methods.We perform several tests that are relevant to simulate explosive scenarios, such as Type Ia supernovae and core-collapse supernovae. We compare the results of the 90 nuclei network with a standardα-chain network with 14 elements to evaluate the differences in the energy generation rate. We also evaluate the relevance of including the electrons in the network in terms of generated yields and how it affects the pressure of a degenerate fluid such as that of white dwarfs. The results obtained with the 90-nuclei network have been verified with a much larger 2000-nuclei network built from REACLIB (WinNet), in terms of nuclear energy generation rate, pressure, and produced yields. Results.The results obtained with the proposed medium-sized network compare fairly well, to a few percent, with those computed withWinNetin scenarios reproducing the gross physical conditions of current Type Ia supernova explosion models. In those cases where the carbon and oxygen fuel ignites at high density, the high-temperature plateau typical of the nuclear statistical equilibrium regime is well defined and stable, allowing large integration time steps. We show that the inclusion of electron captures on free protons substantially improves the estimation of the electron fraction of the mixture. Therefore, the pressure is better determined than in networks where electron captures are excluded, which will ultimately lead to more reliable hydrodynamic models. Explosive combustion of helium at low density, occurring near the surface layer of a white dwarf, is also better described with the proposed network, which gives nuclear energy generation rates much closer toWinNetthan typical reduced alpha networks. Conclusions.A nuclear network withN= 90 species, including electrons, aimed at multidimensional calculations of supernova explosions is described and verified. The proposed network is suitable for the study of Type Ia supernova explosions because it provides better values of pressure and electron abundance than other existing networks with smaller or even a similar size but without including electron capture processes.

Astronomy & Astrophysics↗

The [M 6 (S 2 C 2 Ph 2 ) 6 ] (M = Ni, Pd, Pt) Series: Multielectron Reservoirs That Sustain Ligand-Based Oxidations and Metal-Based Reductions

A complete [M 6 (S 2 C 2 R 2 ) 6 ] series (M = Ni (1), Pd (2), Pt (3); R = Ph), the rarest variety among homoleptic dithiolene transition-metal compounds, has been prepared by reaction between [M(S 2 C 2 Ph 2 ) 2 ] and a M 0 source. The platinum member of this set is the first of its type. Diffraction-quality crystals, grown with high reproducibility by evaporation from PhNO 2 solutions, reveal fully reduced [Ph 2 C 2 S 2 ] 2– dianions and an octahedral M 6 core that is reduced to C 2 symmetry by the fusion of a mononuclear D 2h [M(S 2 C 2 Ph 2 ) 2 ] fragment upon a C 4 -symmetric base. The [Ni 6 (S 2 C 2 Ph 2 ) 2 ] 1– monoanion, prepared by Cp* 2 Co reduction, shows only modest structural differences from its neutral counterpart. In CH 2 Cl 2 , 1 and 2 can undergo two reductions and an oxidation, while 3 sustains two reductions and two oxidations. In benzonitrile, 1 sustains three reversible oxidations at potentials that are shifted appreciably to less positive values. The cathodic processes are shown by density functional theory (DFT) calculations to involve an MO largely of metal–sulfur composition that has contributions throughout the C 4 -symmetric pentametallic base of the assembly, while the oxidations are largely ligand-based and confined to the monometallic [M(S 2 C 2 Ph 2 ) 2 ] cap. The absorption spectra are marked by multiple overlapping bands that produce a continuous, tapering absorption profile of unresolved shoulders and swells.

Ligands↗

Verification and Validation of the STAT7 Code

This document presents the verification and validation of the STAT7 Version 1.1 code, which has been developed to perform the steady-state thermal hydraulics (TH) safety calculations required to support conversion of the Massachusetts Institute of Technology Research Reactor from highly enriched uranium fuel to low- enriched uranium fuel. The key capabilities of the code, including statistical processing, water property generation, and TH solution, are verified independently. While the first two capabilities are checked using statistical tests and NIST data tables, each component of the TH solution is verified using simple geometry test problems. The comparison between results obtained by the code and by hand calculations on these problems shows that each of these components works correctly. Additionally, a comprehensive comparison between STAT7 Version 1.1 and PLTEMP/ANL Version 4.3 Revision 162 calculations is performed for the validation. Good agreement is achieved between the results obtained with the two codes for all TH parameters belonging to the fuel core region.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Reduced Order Models Generation for HTGRs Pebble Shuffling Procedure Optimization Studies

This report provides an initial study for producing reduced-order models (ROMs) of pebble-bed high temperature gas reactor (HTGR) models for the purposes of design optimization. As an initial study, this work is meant to be exploratory---identifying useful workflows and methods for ROM generation---and not meant to be a catch-all analysis of HTGR ROM generation and usage for optimization. This report summarizes three tasks performed in Fiscal Year 2022: 1) the creation of HTGR model, 2) the sensitivity analysis of model design parameters, and 3) an introduction to ROM generation techniques. The representative HTGR model created in this work is a multiphysics equilibrium-core using the BlueCRAB (comprehensive reactor analysis bundle) reactor analysis application, coupling four physical phenomena: neutronics, streamline depletion, porous flow thermal hydraulics, and pebble heat conduction. Part of the model creation was identifying some design parameters and quantities of interest that are relevant in an optimization analysis and adjustable in the model. The sensitivity analysis utilized a polynomial chaos expansion methodology to compute global sensitivity metrics. This analysis showed that thermal hydraulics parameters and quantities of interest had a relatively small impact on simulation results. Finally, the ROM generation work involved exploring three different ROM methodologies: polynomial regression, a Gaussian process, and artificial neural networks. Using a cross-validation technique to characterize ROM performance, the Gaussian process and single-layer artificial neural networks showed the most promising results. Overall, this study was insightful and the lessons learned will be invaluable for the eventual development of an HTGR design optimization workflow.

97 MATHEMATICS AND COMPUTING↗

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↗