Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “feature impact”

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,243 records · Page 69
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Modeling Analysis of Ball-Milling Process for Battery-Electrode Synthesis

The mechanical alloying process is a promising method for synthesizing electrode materials for batteries owing to its benefits such as the ability to produce nanostructured, high-performing electrode alloys, no adverse effects on the solid electrolyte for solid-state batteries, stable production of thick electrodes, simple processing steps, and low processing costs. It is gaining intensive attention in the battery industry as one of the best methods to replace the conventional wet-slurry-solvent method, and its application is rapidly increasing these days. However, the operation is currently conducted purely based on trial-and-error methods without fully utilizing the features of its functions. Here, this may be attributed to a lack of understanding of the effect of operating parameters on the alloying process and final products. Surprisingly, there is a scarcity of the literature conducting fundamental research to comprehend the underlying physics of the entire mechanical alloying process, resulting in a significant knowledge gap. To address this knowledge gap, extensive research was conducted. The existing literature on mechanical alloying was reviewed to comprehend the current state of understanding and to discuss the direction for future research. Mathematical expressions were developed to create physics-based models capable of capturing the entire mechanical alloying process, including milling kinetics and defect-enhanced phase evolution. These methods were then applied to investigate the impact of operating parameters such as milling frequency, initial mole ratio of the alloyed materials, density of grinding balls, and energy required for the powders to become amorphous (i.e., the amorphization energy threshold). This research aimed not only to comprehend the direct effects of these operating parameters but also to unveil the physics underlying the ball-milling process. The results of our study can serve as crucial information for the battery industry in designing or operating the ball-milling process.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Detailed simulations of the first deuterium-tritium-filled double shell implosions on the National Ignition Facility

The first indirectly driven, liquid DT-filled double shell inertial confinement fusion (ICF) implosions have recently been successfully performed on the National Ignition Facility (NIF). Double shells are a class of alternative designs that use a low-Z outer shell to compress a foam cushion that accelerates a high-Z inner shell to efficiently compress a liquid DT core. Double shells are challenging to fabricate, field, and model. Important engineering features enabling double shell fabrication include a fill-tube penetrating all shells and a carefully designed and very narrow (few μm) step-joint in the ablator. Due to the higher density materials involved, high Atwood number instabilities are also important at many material interfaces. In this paper, numerical simulations of double shell implosions using the Los Alamos National Laboratory multi-physics radiation-hydrodynamics code xRAGE will be discussed. An extensive effort has been under way for several years to develop the code capabilities for ICF simulations in a common modeling framework to allow ease of simulation setup and standardization of the computational methodology. This paper will present a wide range of simulation results capturing, quantifying, and comparing the impact of all these degradation mechanisms on implosion performance. Brief comparisons with recent experimental results and suggestions for future improvements will also be discussed. Our results suggest that capsule surface roughness and the step-joint gap have the largest impact on implosion performance. Initial experimental data may suggest that the sensitivity to the step-joint gap could provide the dominant explanation for DT-filled double shell experiments that have been fielded on NIF thus far.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Extracting dark-matter velocities from halo masses: A reconstruction conjecture

Increasing attention has recently focused on nontraditional dark-matter production mechanisms which result in primordial dark-matter velocity distributions with highly nonthermal shapes. Here, in this paper, we undertake an assessment of how the detailed shape of a general dark-matter velocity distribution impacts structure formation in the nonlinear regime. In particular, we investigate the impact on the halo-mass and subhalo-mass functions, as well as on astrophysical observables such as satellite and cluster-number counts. We find that many of the standard expectations no longer hold in situations in which this velocity distribution takes a highly nontrivial, even multimodal shape. For example, we find that the nominal free-streaming scale alone becomes insufficient to characterize the effect of free-streaming on structure formation. In addition, we propose a simple one-line conjecture which can be used to “reconstruct” the primordial dark-matter velocity distribution directly from the shape of the halo-mass function. Although our conjecture is completely heuristic, we show that it successfully reproduces the salient features of the underlying dark-matter velocity distribution even for nontrivial distributions which are highly nonthermal and/or multimodal, such as might occur for nonminimal dark sectors. Moreover, since our approach relies only on the halo-mass function, our conjecture provides a method of probing dark-matter properties even for scenarios in which the dark and visible sectors interact only gravitationally.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Blockchain for Fault-Tolerant Grid Operations

Distribution system fault/failure has a direct impact on customers. Distribution systems are not designed with the contingency of one or more elements failing as the transmission system, however, distribution systems are more likely to have faults in comparison to transmission and generation systems. Most of the distribution systems have a radial design with the protection system assumption of unidirectional power flow. The radial design results in the disconnection of customers for any given component failure. The rapid deployment of distributed energy resources conflicts with the assumption of unidirectional power flow. To maintain/improve the distribution system fault tolerance for a grid with high penetration of distributed energy resources requires improvement. Blockchain can add value to improve fault-tolerant grid operations. That can be achieved using blockchain’s core features of distributed consensus-based decision-making process and immutability. In the process of preparing for an event or system restoration, there is a need for reliable data and system situational awareness. The system situational awareness enables accurate planning for possible scenarios and system restoration after an event, an area where blockchain’s distributed architecture can help. The value of blockchain for a fault-tolerant grid operation is elaborated with three use cases: 1) fault location, isolation, and service restoration; 2) topology identification; and 3) data configuration. These use cases demonstrate how blockchain-based architecture can facilitate bridging the vulnerabilities of the existing implementation.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2026

Prairie State Generating Company 2D Seismic Interpretation

The objectives of the Illinois Storage Corridor (ISC) project are to accelerate commercial deployment of carbon capture utilization and storage at two individual sites and receive approvals for Underground Injection Control (UIC) Class VI permits for construction at each site (ISC Project Narrative, 2020). As part of this project, and as part of the subsurface geologic characterization, two-dimensional (2D) seismic data was acquired at both sites. This report summarizes the findings from two phases of 2D acquisition and seismic interpretation at the Prairie State Generating Company site near Marissa, Illinois. The seismic data indicates the presence of three faults that completely transect the storage and confining units. The initial 2021 2D seismic acquisition revealed the presence of a feature two miles east of the Lively Grove #1 characterization well. The feature is a polyphase fault with a component of strike-slip motion, forming a small positive flower structure. The second phase of 2D seismic acquisition in 2022 constrained the maximum extent of the fault, indicating it has a relatively limited length. Other 2022 seismic lines revealed the presence of two other faults that transect the storage and confining units in the southwest and northeast portions of the 2D seismic acquisition area. The 2021 and 2022 2D seismic surveys have identified three specific locations within the project area that may have an elevated risk of out of zone CO 2 migration due to faulting. Further technical work will be needed to quantify this risk and determine how and if this risk will impact the placement of CO 2 injection wells across the project area.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2020

Use of Remote Sensing and In-Situ Observations to Develop and Evaluate Improved Representations of Convection and Clouds for the ACME Model

The overachieving goal of the whole CMDV-MCS project is to improve understanding of warm season continental convection and to develop treatments of convection and microphysics capable of representing mesoscale convective systems (MCSs) features in large-scale models. Our tasks for this project contributing to the overachieving goal include: (1) Improve the ice nucleation formulation for MG2 and P3 cloud microphysics schemes; (2) Improve the treatment of subgrid dynamics and thermodynamics driving the ice nucleation in E3SM; and (3) Test the performance of improved ice microphysics in E3SM with observation data. In this project, we have (1) Improved the ice nucleation parameterization for MG2 and P3 in E3SM by implementing two advanced empirical parameterizations with connection to aerosols. The two deterministic heterogeneous ice nucleation parameterizations (i.e., DeMott et al., 2015; Niemand et al., 2012) were merged with the MG2 and P3 cloud microphysics schemes in E3SM. Long-term simulations were conducted to examine the impacts of these new parameterizations on simulated cloud properties; (2) Improved the treatment of subgrid dynamics and thermodynamics driving the ice nucleation in E3SM. We evaluated the double Gaussian PDF of vertical velocity simulated by the Cloud Layers Unified By Binormals (CLUBB) and the sub-column vertical velocity sampled from the Subgrid Importance Latin Hypercube Sampler (SILHS) in E3SM. We introduced the vertical velocity variance induced by topographic gravity waves for ice nucleation and droplet activation; and (3) Tested the performance of improved ice microphysics in E3SM with observation data. We tested the new treatments of ice nucleation in the single column model (SCM) mode for the stratiform mixed-phase clouds observed during 9-10 October 2004 in the DOE ARM Mixed-Phase Arctic Cloud Experiment (M-PACE) and for the convective clouds observed on 20 May 2011 in the Midlatitude Continental Convective Clouds Experiment (MC3E). Modeled ice nucleating particles (INPs) concentrations were compared against observations collected around the globe.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Relative molecular orientation can impact the onset of plasticity in molecular crystals

Abstract Creating or moving dislocations is the first step to dissipating mechanical energy via plastic deformation under contact loading. In molecular crystals there is both a lattice that defines crystal orientation and a relative orientation of the basis of the molecules. We define a normalization parameter which relates strain at yield, the hardness of the bulk crystal, and a distance parameter analogous to a Burgers vector that nominally predicts the relative ease of initiating plasticity in this broad class of materials. Analyzing the yield behavior of 10 different molecular crystals of varying space groups shows the inter-molecular orientation predicts the experimentally observed applied stress needed to nucleate dislocations. When molecules are oriented ‘parallel’ relative to one another the normalized maximum shear stress at the onset of plasticity is on the order of 3–5 times lower than when molecules within the crystal are ‘anti-parallel’, and molecules with a more equiaxed shape fall in between these bounds. This provides an initial indication of a structural feature which predicts the relative ease of initiating plasticity during contact loading in molecular crystals.

36 MATERIALS SCIENCE↗
DOE OSTI2021

The pervasive impact of critical fluctuations in liquid–liquid extraction organic phases

Liquid-liquid extraction is an essential chemical separation technique where polar solutes are extracted from an aqueous phase into a nonpolar organic solvent by amphiphilic extractant molecules. A fundamental limitation to the efficiency of this important technology is third phase formation, wherein the organic phase splits upon sufficient loading of polar solutes. The nanoscale drivers of phase splitting are challenging to understand in the complex hierarchically structured organic phases. In this study, we demonstrate that the organic phase structure and phase behavior are fundamentally connected in a way than can be understood with critical phenomena theory. For a series of binary mixtures of trialkyl phosphate extractants with linear alkane diluents, we combine small angle x-ray scattering and molecular dynamics simulations to demonstrate how the organic phase mesostructure over a wide range of compositions is dominated by critical concentration fluctuations associated with the critical point of the third phase formation phase transition. These findings reconcile many longstanding inconsistencies in the literature where small angle scattering features, also consistent with such critical fluctuations, were interpreted as reverse micellar-like particles. Altogether, this study shows how the organic phase mesostructure and phase behavior are intrinsically linked, deepening our understanding of both and providing a new framework for using molecular structure and thermodynamic variables to control mesostructure and phase behavior in liquid-liquid extraction.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

ADAM: A web platform for graph-based modeling and optimization of supply chains

Modeling and optimization are essential tasks that arise in the analysis and design of supply chains (SCs). SC models are essential for understanding emergent behavior such as transactions between participants, inherent value of products exchanged, as well as impact of externalities (e.g., policy and climate) and of constraints. Unfortunately, most users of SC models have limited expertise in mathematical optimization, and this hinders the adoption of advanced decision-making tools. Here, in this work, we present ADAM, a web platform that enables the modeling and optimization of SCs. ADAM facilitates modeling by leveraging intuitive and compact graph-based abstractions that allow the user to express dependencies between locations, products, and participants. ADAM model objects serve as repositories of experimental, technology, and socio-economic data; moreover, the graph abstractions facilitate the organization and exchange of models and provides a natural framework for education and outreach. Here, we discuss the graph abstractions and software design principles behind ADAM, its key functional features and workflows, and application examples.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Enhanced Detection of Primary Biological Aerosol Particles Using Machine Learning and Single-Particle Measurement

Accurately identifying primary biological aerosol particles (PBAPs) using analytical techniques poses inherent challenges due to their resemblance to other atmospheric carbonaceous particles. Here, we present a study of an enhanced method for detecting PBAPs by combining single-particle measurement with advanced supervised machine learning (SML) techniques. We analyzed ambient particles from a variety of environments and lab-generated standards, focusing on chemical composition for traditional rule-based and clustering approaches and incorporating morphological features into the SML approaches, neural networks and XGBoost, for improved accuracy. This study demonstrates that SML methods outperform traditional methods in quantifying PBAPs, achieving significant improvements in precision, recall, F1-score, and accuracy, leading to an increased number of detected PBAPs by at least 19%. The adaptability of the proposed XGBoost-based SML model is showcased in comparison to traditional methods in categorizing PBAPs for blind data sets from different geographical locations. Two field case studies were investigated, over agricultural land and Amazonia rain forest, representing relatively low and high concentrations of PBAPs, respectively, where XGBoost consistently detected up to 3.5 times more PBAPs than traditional methods. Precise detection of PBAPs in the atmosphere could significantly improve the prediction of climatic impacts by them.

42 ENGINEERING↗
DOE OSTI2024

Use of Satellite, Surface Observations and Numerical Weather Prediction Model Data to Improve Cloud Base Height and Cloud Base Vertical Velocity Estimation

Cloud base height (CBH) and cloud base vertical velocity (CBVV) are important variables that impact the overall climate in a region as they influence the formulation, longevity, and evolution of clouds. Retrieval of both parameters have long used ground instrumentation (e.g., Doppler lidar (DL), ground base radar); however, retrieving CBH from satellites is particularly challenging given that space-based instruments only observe cloud tops. In this manuscript, CBH is retrieved using a multi-linear regression equation, while CBVV used a random forests model. Both retrievals combine satellite and numerical weather prediction data. The satellite data used are the Visible Infrared Imaging Radiometer Suite imagery, while measurements of CBH and CBVV include DL and radiosonde data at the Southern Great Plains (SGP) Atmospheric Radiation Measurement observatory. Data from 83 summer days (May-August) in 2018–2021 featuring cumulus clouds forced by solar heating were examined and used to train the models, with years 2022–2023 used for validation. Various spatial domains were defined with one large (2.4° longitude by 2.0° latitude) SGP domain being split into smaller sections (smallest being 0.99° and 0.61° longitude and latitude respectably). CBH and CBVV values obtained from the DL as compared to the models show root mean square errors between 150 and 200 m, with CBVV values between 0.45 and 1 ms -1 . Finally, it was found that the CBH formulation performs well over all domains, while the CBVV retrievals become less accurate due to more turbulence being introduced into the observations as the number of DL stations decreases in the smaller domains.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Impact of molten salt inflow on the temperature distribution in thermal energy storage tanks at startup for central receiver concentrating solar power plants

Concentrating Solar Power (CSP) systems with molten salt thermal energy storage (TES) tanks are one of the most promising, renewable-based energy conversion technologies for larger-scale power generation. The TES tank is one of the most critical components in CSP plants due to its high-temperature operation (up to 565 °C), daily thermal cycling, and intermittent solar radiation conditions. The plant startup is one of the most challenging operation conditions that could lead to damaging thermal gradients due to low salt inventory levels. In this study an analytical model for the sparger ring was developed and integrated with a detailed computational fluid dynamics model of a commercial-scaled molten salt tank. The integrated model allows an accurate representation of the tank operation to evaluate the effect of molten salt inflow on the mixing process. The tank filling process during plant startup was analyzed considering sparger rings with variations in design features, including inlet orifice configurations, number of orifices, direction of the inlets, and orifice diameter. The results demonstrated that higher temperature gradients are obtained in the tank floor during the plant startup. A sparger ring configuration with a predetermined orifice inlet inclination (30°, 45° and 60°) leads to significant temperature differences in the floor, between 57 °C and 62 °C, but better homogeneity in the temperature of the salt inventory. Lower salt inflow velocities result in a more homogeneous floor temperature, with maximum temperature differences under 38 °C. The sparger ring configuration with 52 orifices of 1-in. diameter and vertical flow showed better homogeneity in the temperature differences as a function of the salt level and lower temperature gradients in the tank floor. Because large temperature gradients in the tank's floor have been identified as one of the main factors contributing to tank failures, assessing various sparger ring design features is fundamental to determining proper inflow conditions that lead to low-temperature gradients and reducing failure susceptibility.

14 SOLAR ENERGY↗
DOE OSTI2021

Deep learning for Alzheimer's disease: Mapping large-scale histological tau protein for neuroimaging biomarker validation

Abnormal tau inclusions are hallmarks of Alzheimer's disease and predictors of clinical decline. Several tau PET tracers are available for neurodegenerative disease research, opening avenues for molecular diagnosis in vivo. However, few have been approved for clinical use. Understanding the neurobiological basis of PET signal validation remains problematic because it requires a large-scale, voxel-to-voxel correlation between PET and (immuno) histological signals. Large dimensionality of whole human brains, tissue deformation impacting co-registration, and computing requirements to process terabytes of information preclude proper validation. We developed a computational pipeline to identify and segment particles of interest in billion-pixel digital pathology images to generate quantitative, 3D density maps. The proposed convolutional neural network for immunohistochemistry samples, IHCNet, is at the pipeline's core. We have successfully processed and immunostained over 500 slides from two whole human brains with three phospho-tau antibodies (AT100, AT8, and MC1), spanning several terabytes of images. Our artificial neural network estimated tau inclusion from brain images, which performs with ROC AUC of 0.87, 0.85, and 0.91 for AT100, AT8, and MC1, respectively. Introspection studies further assessed the ability of our trained model to learn tau-related features. We present an end-to-end pipeline to create terabytes-large 3D tau inclusion density maps co-registered to MRI as a means to facilitate validation of PET tracers.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Holistic Marine Energy Resource Assessments: A Wave and Offshore Wind Perspective of Metocean Conditions

Offshore wind and wave energy resources are possibly the last significantly untapped renewable energy resource, and could play a significant role in mitigating the worst impacts of climate change via the generation of renewable electricity. However, offshore wind and wave energy resource assessments have been conducted as independent parallel processes, with little interaction with regard to best practices, lessons learnt, or opportunities to create compatible methodologies for future utilization by the broader marine energy sector. Based on the latest technical specifications from the International Electrotechnical Commission, and the highest fidelity publicly available datasets, the offshore wind and wave conditions at the PacWave site off Oregon, USA were quantified. The results clearly show a significant untapped energy resource an annual average wave energy flux of ~35kW/m and a mean average wind speed of 7.8m/s. The offshore wind and wave energy resources both show significant seasonal variation, with offshore wind also featuring a consistent daily profile during summer. Finally, opportunities and challenges associated with developing a holistic assessment of offshore marine energy resources were discussed and recommendations provided.

Robertson, Bryson↗
DOE OSTI2025

Electric dipole excitations near the neutron separation energies in 96 Mo

Electric dipole strength near the neutron separation energy significantly impacts nuclear structure properties and astrophysical scenarios. These excitations are complex in nature and may involve the so-called pygmy dipole resonance (PDR). Transition densities play a crucial role in understanding the nature of nuclear excited states, including collective excitations, as well as in constructing transition potentials in DWBA or coupled-channels equations. In this work, we focus on electric dipole excitations in spherical molybdenum isotopes, particularly 96 Mo, employing fully consistent Hartree-Fock-Bogoliubov (HFB) and Quasiparticle Random Phase Approximation (QRPA) methods. We analyze the dipole strength near the neutron separation energy, which represents the threshold for neutron capture processes, and examine the isospin characteristics of PDR states through transition density calculations. Examination of proton and neutron transition densities reveals distinctive features of each dipole state, indicating their isoscalar and isovector nature. We observe that the primary component in the enhanced low-energy region exhibits isovector character. The PDR displays a mixture of isoscalar and isovector nature, distinguishing it from the isovector giant dipole resonance (IVGDR). These findings lay the groundwork for future investigations into the role of transition densities in reaction models and for their application to inelastic scattering calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

A Verification of Flux Sensitivity Estimates Using the MCNP Tally Perturbation Tool

Nuclear data is commonly used in applications such as nuclear nonproliferation, safeguards, and criticality safety. More specifically, nuclear data is used in predictive simulation codes like the Monte-Carlo N-Particle (MCNP ® ) transport code, Serpent, and similar radiation transport codes. The improvement of nuclear data enables more precise and accurate simulations, which result in higher fidelity designs and reduced operational/procedural costs. Therefore, the improvement of nuclear data is of paramount importance across the nuclear community. Nuclear data is improved and validated through integral benchmark experiments. The design of benchmark experiments is an extensive process; therefore, these experiments are often optimized on multiple characteristics, including sensitivity to the nuclear data, during the design process. Sensitivity is a measure of how much a quantity changes due to changes in independent variables such as experimental configuration. An experimental design that has a larger sensitivity to the nuclear data of interest will have a larger impact on the accuracy and precision of the validated data. Past integral benchmark experiments have primarily used the effective multiplication factor ($k_{eff}$) as the predominant measured quantity; however, experiments designed with other quantities in mind would be able to optimize on validating different areas of the nuclear data. A primary goal of the EUCLID project is to design, constrain, and reduce compensating errors in experiments focused on quantities other than $k_{eff}$ to better validate nuclear data across the board. Currently, there is a capability in MCNP to easily calculate the sensitivity of $k_{eff}$ to specific nuclear data of numerous reactions types and isotopes (KSEN card); however, the sensitivity of other quantities must be estimated in more strenuous manners. For example, the perturbation feature (PERT card) of MCNP can be used to estimate first-order sensitivities of some response in fixed source simulations. A recent announcement revealed that the first- and second-order perturbation features in previous releases of MCNP contained a bug. It was identified that particles were being scored into the wrong energy bin. The bug is in the most recent public release (MCNP6.2); however, a patch has been added to the most up to date version (MCNP6.2.2) that has not been released publicly. A direct comparison of the PERT card results for an F4 (neutron flux averaged over a cell) tally before and after the patch are shown in figure 1. All simulations used in the sensitivity estimates in this report were performed with MCNP6.2.2. This work verifies the patched MCNP perturbation tool by comparing first order sensitivities made using the PERT card to estimates made using manual perturbation of the compact ENDF (ACE) files.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Exploration of Rayleigh-Taylor seeding mechanisms in laser- and pulsed-power-driven inertial fusion

The ubiquitous Rayleigh-Taylor instability (RTI) has been observed in fluids and plasmas in a wide range of parameter regimes from astrophysical to laboratory plasmas. In high-energy-density (HED) laboratory plasmas, such as laser-driven inertial confinement fusion (ICF) plasmas, the RTI can have a detrimental impact on achieving fusion ignition through the generation of hydrodynamic mix. Hence, an understanding of the seeding mechanisms that produce the RTI and identifying ways to mitigate hydrodynamic mix is of critical importance. Appropriately aligned magnetic fields have been known to stabilize short-wavelength RTI. A number of numerical and experimental studies have demonstrated the benefits of using imposed magnetic fields in laser-driven ICF to achieve higher ion temeratures, higher neutron yields, and a relaxation of the ignition criteria. This work addresses seeding mechanisms in laser-driven implosions to understand critical early-stage physics that ultimately leads to substantial growth of the RTI along with mechanisms for mitigation of this growth. Surface perturbations due to machining tolerances and single-feature seeds (for example due to fill tubes) can produce substantial RTI growth in the ignition-relevant high-convergence ratio targets for laser- driven ICF implosions. RTI growth from single-feature seeding in laser-based implosions has the potential to be mitigated through appropriately aligned externally applied magnetic fields and the goal of this work has been to quantify that numerically and experimentally. Resistive-magnetohydrodynamic (MHD) simu- lations are used to study the seeding and evolution of the RTI leveraging previous support from the DOE HEDLP program. Furthermore, experimental data has been obtained and applied towards code validation of unmagnetized and magnetized evolution of single-feature seeded RTI growth.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Developing And Scaling an OpenFOAM Model to Study Turbulent Flow in a HFIR Coolant Channel

Improving the understanding of how computational fluid dynamics (CFD) direct numerical simulations (DNS) of flows in the High Flux Isotope Reactor (HFIR) perform when run in parallel using the high performance computing (HPC) platform Summit at the Oak Ridge Leadership Computing Facility (OLCF) is of particular importance to boost the computational tools used to support HFIR conversion to low enriched fuel (LEU). Evaluation of scaling performance was driven by the increasing importance of graphics processing unit (GPU) usage in HPC, which is becoming the standard for modern supercomputers such as Summit. The desired results are to obtain a strong positive correlation between the computational resources dedicated to a problem and the relative speed-up of the simulation in comparison to a benchmark. This capability will allow substantially improvement in HFIR flow analytical capabilities, specifically when predicting turbulence properties at high Reynolds numbers. The study leverages previous simulation results performed with code PHASTA (finite element) on HPC platforms Cori (NERSC) and Theta (ALCF) [1] with computing options provided in the computing platform OpenFOAM (finite volume) at OLCF. Transitioning from PHASTA to OpenFOAM will (1) eliminate dependence on third-party software for mesh generation and manipulation, (2) reduce resource needs by employing modern architectures, and (3) build expertise for future modeling of HFIR-specific problems like heat transfer in involute geometry, entrance effects, flow structure in channel corners, and so on—all important issues when defining the available thermal margins in the transition to LEU. CPUs and GPUs differ significantly in their architecture and utilization, as discussed in the literature [2]. The most important differences are in the approach to computations and their memory. A single GPU contains a large quantity of cores, enabling it to perform with a much higher throughput than a CPU, but execution requires a different approach. GPU codes execute instructions using the Single-Instruction Multiple-Thread (SIMT) approach in which a single instruction is used for groups of threads called warps. A warp typically consists of 32 threads which must execute the same set of instructions, although on separate threads. Alternately, a CPU has far fewer cores that are much more flexible in their operation, excelling at quickly performing more complex serial computations. This is why GPUs have greater throughput when properly utilized. The second important difference is seen when comparing their memory spaces. Limited memory allocations and CPU–GPU communications cause a significant bottleneck in GPU-accelerated programs. Further study was required to properly take advantage of GPU resources. A comprehensive analysis of code performance and the model-specific features of turbulence constitutes the core of this work. In this study, a DNS simulation of HFIR channel turbulence was performed with the finite volume CFD code OpenFOAM v2112 and CUDA v11.0 on Red Hat Enterprise Linux v8.2. The OpenFOAM installation had AMGx integrated to enable GPU acceleration and utilizes the PETSc4FOAM library. The computational resources and the problem size were scaled on CPU and CPU + GPU architectures to gain a better understanding of the performance of a DNS problem on modern computing hardware. The study aimed to analyze the scaling of the code exclusively on CPUs and then to examine the scaling of the codes with GPU acceleration enabled. Scaling studies included CPU and GPU acceleration on a mesh of varying resolution to analyze the impact of problem size relative to computational resources. In the course of preparing the GPU configuration on Summit, mainly using the AMGX solvers, difficulties were encountered stemming from constant changes resulting from extensive ongoing development activities and the changing environment. This resulted in the inability to complete the GPU portion of the work. The code was compiled and tested, but production runs to assess acceleration were not performed because the used discretional compute time allocation expired as year-end approached. The Summit HPC platform is scheduled for decommissioning in 2024, making it unattractive for future use with Nvidia-based GPUs. Therefore, the work will be moved onto NERSC machines in FY24. An application was prepared and submitted, and sufficient node-hours were awarded to continue the research in the next calendar year. This report summarizes work performed thus far, which mostly focused on CPU OpenFOAM computing.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗