Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Relative accuracy”

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,045 records · Page 58

Linear point and sound horizon as purely geometric standard rulers

The baryon acoustic oscillations feature (BAO) imprinted in the clustering correlation function is known to furnish us cosmic distance determinations that are independent of the cosmological-background model and the primordial perturbation parameters. These measurements can be accomplished rigorously by means of the purely geometric BAO methods. To date two different purely geometric BAO approaches have been proposed. The first exploits the linear-point standard ruler. The second, called correlation-function model-fitting, exploits the sound-horizon standard ruler. A key difference between them is that, when estimated from clustering data, the linear point makes use of a cosmological-model-independent procedure to extract the ratio of the ruler to the cosmic distance, while the correlation-function model-fitting relies on a phenomenological cosmological model for the correlation function. Nevertheless the two rulers need to be precisely defined independently of any specific observable (e.g., the BAO). We define the linear point and sound horizon and we fully characterize and compare the two rulers’ cosmological-parameter dependence. We find that they are both geometrical (i.e., independent of the primordial cosmological parameters) within the required accuracy, and that they have the same parameter dependence for a wide range of parameter values. We estimate the rulers’ best-fit values and errors, given the cosmological constraints obtained by the Planck Satellite team from their measurements of the cosmic microwave background temperature and polarization anisotropies. We do this for three different cosmological models encompassed by the purely geometric BAO methods. In each case we find that the relative errors of the two rulers coincide and they are insensitive to the assumed cosmological model. Interestingly both the linear point and the sound horizon shift by 0.5σ when we do not fix the spatial geometry to be flat in Λ CDM. Finally, this points toward a sensitivity of the rulers to different cosmological models when they are estimated from the cosmic microwave background.

79 ASTRONOMY AND ASTROPHYSICS↗

Evaluation of Thermal Neutron Scattering Cross Section of Uranium Silicide with Ab Initio Lattice Dynamics

Uranium silicide (U 3 Si 2 ) is a candidate material for the high-density nuclear fuel in commercial light water reactors [1], [2]. Its higher uranium density, 11.3 g-U/cm3, compared to that of uranium dioxide (UO 2 ), 9.7 g-U/cm3, can improve the performance of a nuclear reactor while using low enriched uranium (LEU) and diversify the choice of cladding materials [1]–[3]. It also has a higher thermal conductivity than UO 2 , which can reduce the thermal stress on the material caused by a temperature gradient across the fuel pellet and provide a larger margin for some postulated accidents [1], [2], [4]. Furthermore, compared to U3Si, another high-density fuel candidate, it has better resistance to in-pile swelling due to less irradiation-induced rapid amorphization [1], [3]. Corresponding to its importance in nuclear engineering, many previous studies have reported the properties of U3Si2. Experiments showed that U 3 Si 2 is a paramagnetic (PM) metal, where a slight linear increase in magnetic susceptibility was measured with increasing temperature [5], [6]. In addition, thermodynamic quantities such as thermal expansion coefficient, heat capacity, and thermal conductivity were experimentally determined over a wide temperature range [1], [7], [8]. In several computational studies, ab initio atomistic simulations based on density functional theory (DFT) were performed to calculate various properties including elastic constants, electronic density of states (DOS), and phonon dispersion curves [9]–[12]. Nevertheless, thermal neutron scattering cross sections, which are critical to the prediction of the parameters in reactor physics that are ultimately related to reactor criticality, have not yet been evaluated for U3Si2. The scattering cross section can be calculated from the phonon DOS, or the energy spectrum of lattice vibrations, of the crystalline system [13], [14]. However, there is also no experimental data available for the phonon DOS of U 3 Si 2 . While some computational studies reported the phonon DOS and/or dispersion curves from ab initio simulations [9]–[12], the accuracy cannot be guaranteed because it is unclear whether the spin-polarization behavior of PM U 3 Si 2 was properly described. In the present study, the thermal neutron scattering cross section for U 3 Si 2 is evaluated for the first time by calculating the phonon DOS for U3Si2 from ab initio lattice dynamics (AILD) simulations based on DFT. First, U 3 Si 2 is modeled based on the experimental structure, and AILD simulations are performed on the modeled U3Si2 to optimize the structure. Next, AILD simulations are performed for supercells with atomic displacement to calculate Hellmann-Feynman forces. Based on the calculated forces, partial phonon DOSs for U and Si are obtained, and the thermal neutron scattering law (TSL) for U 3 Si 2 is finally evaluated. To verify the accuracy of the calculations in the present study, the calculation results are compared with experimental data on the structure and heat capacity of U3Si2 [1], [7], [8], [15].

Geometry Optimization↗

Horn Location Sensors (HLS) for LBNF

The Long-Baseline Neutrino Facility (LBNF) will deliver the world's most powerful muon neutrino beam to the Deep Underground Neutrino Experiment (DUNE), initially operating at 1.2 MW and upgradeable to 2.4 MW. Ensuring the accurate direction of this beam is critical for DUNE's precision goals. This talk introduces the Horn Location Sensors (HLS) system, designed to provide precise, relative measurements of the focusing horns, targets, and beam position monitors in the neutrino beamline. The HLS system employs high-precision FSI-based hydrostatic leveling sensors to track vertical motion and tilt, achieving precision on the order of 0.1 mm. Built for minimal maintenance in radioactive environments, the HLS system ensures precise alignment of beamline components during high-power operations. This system is essential for maintaining beam accuracy and enhancing DUNE's scientific performance.

43 PARTICLE ACCELERATORS↗

Comparative Study of Hydrology and Icemelt in Three Nepal River Basins Using the Glacio-Hydrological Degree-Day Model (GDM) and Observations From the Advanced Scatterometer (ASCAT)

An assessment of the water supply and its seasonal and annual changes over the century in the High Mountain Asia (HMA) region is of increasing interest due to its potential impact on one-sixth of the global population. In order to understand the changing hydrology and snow and ice melt, we used remotely sensed Advanced Scatterometer (ASCAT) observations of glacier melt (GM) and a distributed and gridded Glacio-hydrological Degree-day Model (GDM) in three river basins: Tamor, Trishuli and Marsyangdi. The GDM-estimated contribution of snowmelt, icemelt, rainfall and baseflow in river flows is found to be most accurate in the Trishuli River basin, with Nash-Sutcliffe efficiency (NSE) between the estimated and observed discharges of 0.81 and volume differences of -0.5%, and reasonably accurate in the Tamor River basin, with NSE of 0.69 and volume difference of -7.51%. Similarly, NSE of 0.81 and volume difference of 4.64% in Marsyangdi River basin. We find strong similarities in the timing of glacier melting using the GDM and from observations from the ASCAT GM, determining the seasonal start of glacier melting to within 6 days on average. In all basins ASCAT GM observes melting at higher elevations relative to GDM, average of 5,328 m a.s.l. Systematic differences in glacier melting area determined by modeling and satellite observations indicate ASCAT may have suboptimal resolution, view geometry and/or polarimetry for delineating glacier melting at the process-scale in complex topography, especially in the ablation zone. This is the first step in examining the remote sensing products that could potentially be incorporated into hydrologic models to increase the accuracy of the hydrologic flow as well as the ability to estimate river discharge in other basins with limited data.

54 ENVIRONMENTAL SCIENCES↗

Need for advanced research reactors for the next-generation reactor physics, analysis tools, and technology

Full text of publication follows. There is an urgent need for design and deployment of advanced research and test reactors in support of design, licensing and operation of advanced power reactors and education of next generation nuclear workforce. Existing research reactors mostly were designed and constructed decades ago with the main objectives of training operators, performing reactor physics experiments, and educating nuclear engineers and scientists. There are already gaps and significant concern about future capabilities for the existing research reactor facilities to address modern instrumentation and/or flexible environments for: performing reactor physics studies for advanced designs which have significantly different core materials forms and compositions, reactor shapes and size; validation of advanced high-fidelity software; development of machine learning algorithms for enhancement of human-machine collaboration in support of reactor monitoring, operation and safeguards; and, effective education of the next-generation workforce. The authors will focus on the need for advanced research reactors to improve and validate fast and accurate simulation tools for high-fidelity modeling and analysis of nuclear reactors in support of their design, optimization, licensing, operation, and monitoring. In the past, simulation tools were limited to relatively coarse models using approximate methodologies that benefited from two main factors: i) allowance for large margins and tolerances; ii) ability to construct prototype (e.g., zero power) reactors for adjustment of approximate methodologies. The next generation reactors have to be designed mainly by using novel high-fidelity computational tools that are accurate and fast, and therefore can be used for parametric studies and uncertainty quantification. To sufficiently demonstrate the accuracy of these tools, advanced research reactors are needed. The authors argue the need for new computational paradigms such as the MRT (Multistage, Response- function Transport) methodology which has resulted in the development of the novel high-fidelity RAPID (Real-time Analysis for Particle-transport and In-situ Detection) code system. Such code systems have to be robust in modeling any complex system, and should be fast and accurate, henceforth their uncertainties can be quantified at reasonable costs. Again, advanced research reactors are needed for the validation of the fidelity and accuracy of new computational tools. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Evaluation of the beam-induced depolarization of the HJET target at the EIC

The Polarized Atomic Hydrogen Gas Jet Target (HJET) has played a central role in the absolute calibration of proton beam polarization at RHIC and is foreseen as a key component of the hadron polarimetry program at the future Electron–Ion Collider (EIC). The substantially higher beam current, reduced bunch spacing, and shorter bunch length planned for EIC operation motivate a careful reassessment of possible beam-induced depolarization of the jet target. In this paper, the depolarization of ground-state hydrogen atoms caused by the time-dependent magnetic field of the circulating polarized proton beam is quantitatively evaluated. The hydrogen atom is treated as a four-level hyperfine system in a holding magnetic field, and transitions driven by harmonic components of the bunch-induced magnetic field are analyzed using time-dependent quantum-mechanical evolution along atomic trajectories. Numerical tracking of hydrogen atoms through the beam region is performed using nominal EIC beam parameters. It is shown that, for a holding field of 120 mT (as used at RHIC), the resulting depolarization of the jet target at the EIC is negligibly small, ≲ 0.01 %, and well below the level relevant for EIC polarization accuracy requirements. The stability of this result with respect to plausible variations of the EIC proton beam parameters is also evaluated. In addition, possible effects under alternative experimental conditions are also examined.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

High-accuracy absolute magnetometry with application to the Fermilab Muon g-2 experiment

We present details of a high-accuracy absolute scalar magnetometer based on pulsed proton NMR. The B-field magnitude is determined from the precession frequency of proton spins in a cylindrical sample of water after accounting for field perturbations from probe materials, sample shape, and other corrections. Features of the design, testing procedures, and corrections necessary for qualification as an absolute scalar magnetometer are described. The device was tested at B = 1.45 T but can be modified for a range exceeding 1–3 T. The magnetometer was used to calibrate other NMR magnetometers and measure absolute magnetic field magnitudes to an accuracy of 19 parts per billion as part of a measurement of the muon magnetic moment anomaly at Fermilab.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Development and Test Results of Thermocouples Used in the TRISO-Fuel Irradiation Experiment AGR-5/6/7

High-temperature gas reactor irradiation experiments create unique challenges for thermocouple-based temperature measurements. High-temperature industrial thermocouples suffer rapid decalibration due to transmutation of the thermoelements from neutron absorption. For lower-temperature applications, Type K and Type N thermocouples are affected by neutron irradiation only to a limited extent. But until recently, the use of these nickel-based thermocouples was limited when the temperature exceeded 1050°C due to drift related to phenomena other than nuclear irradiation. Certain portions of the AGR-5/6/7 experiment experienced temperatures higher than any of the previous AGR tests, up to 1500°C. Recognizing the limitations of existing thermometry to measure such high temperatures, the sponsor of the AGR-5/6/7 test supported a development and testing program for thermocouples capable of low-drift operation at temperatures above 1100°C. This program included additional development of high-temperature irradiation-resistant thermocouples based on molybdenum/niobium thermoelements, which have been studied at Idaho National Laboratory since circa 2004. A step change in accuracy and long-term stability of this thermocouple type was achieved as part of the AGR-5/6/7 thermometry development program. Additionally, long-term testing (7000+ h) at 1250°C of Type N thermocouples utilizing a customized sheath developed at the University of Cambridge has been completed with excellent low-drift results. The results of this testing as well as testing of the improved high-temperature irradiation-resistant design are reported herein. Both the improved high-temperature irradiation-resistant and the Cambridge Type N thermocouple types were incorporated into the AGR-5/6/7 test, which began irradiation in February 2018 and was completed in July 2020. Finally, a summary of the performance of the thermocouples incorporated into the AGR-5/6/7 test is included herein.

In-pile instrumentation↗

Computational and Experimental Investigation of Chiral and Achiral Two‐Dimensional Organic Lead Bromide Perovskites: Octahedral Distortions and Electronic and Optical Properties

A computational investigation is presented, in conjunction with synthesis and experimental characterization, into the structural, electronic, and optical properties of layered two-dimensional organic lead bromide perovskites. Materials based on the chiral (R/S)-4-fluoro-α-methylbenzylammonium (R/S-FMBA), which have been shown to lead to bright room-temperature circularly polarized luminescence, are contrasted with the similar achiral 4-fluorobenzylammonium (FBA). Using density functional theory (DFT) with van der Waals (vdW) corrections, relaxed structures (compared with X-ray diffraction, XRD) and optical absorption spectra (compared with experiments) are studied, as well as band structure and orbital character of transitions. A Python code is developed and provided to calculate octahedral distortions and compare DFT and XRD results, finding that vdW corrections are important for accuracy and that DFT overestimates octahedral tilt angles. (FMBA) 2 PbBr 4 shows among the largest tilt angle differences (often termed Δ β ) reported, 14°–15°, indicating strong inversion symmetry-breaking, which enables its chiral emission. A large resulting Dresselhaus spin-splitting effect is found. The lowest-energy optical transitions involve the perovskite only and are polarized within the layer. This work furthers understanding of structure-property relations with applications to optoelectronics and spintronics.

UV/vis spectroscopy↗

Impact of uranium oxide (UO 2 ) fuel with molybdenum (Mo) inserts on pressurized water reactor performance and safety

This work investigates nuclear reactor performance and safety characteristics of UO 2 with high thermal conductivity Mo insert structures by using multiphysics modeling techniques. Additionally, the purpose of this study is to use scoping analyses to quantify the impact of using Mo inserts from neutronic and heat transfer standpoints. Attention is given to reactor performance parameters, such as cycle length, maximum fuel temperature, temperature gradients in the fuel, and stored energy in the fuel. The finite-element code BISON and the Monte Carlo particle transport code Serpent were used to perform sensitivity analyses on the Mo insert geometry to optimize the insert design and inform larger scale modeling that required the homogenization of the UO 2 and Mo. Although BISON is often used as a fuel performance analysis tool, it is used in this context for heat transfer analysis only. Fuel performance optimization is outside the scope of the current study, but would be important for future work focused on this concept. The results showed that the insert had little impact on neutronic performance and that homogenizing the UO 2 and Mo was acceptable for reactor physics calculations. Reactivity temperature coefficients calculated using homogeneous UO 2 -Mo were shown to be relatively similar to UO 2 , but higher Mo content and 235 U enrichment can reduce the worth of soluble boron and control rods. The effect of insert geometry on heat transfer was much greater, and an approximately 15–20% difference in maximum fuel temperature was predicted between the best and worst performing heat transfer geometries. Furthermore, thermal conductivity calibration based on the finite element analysis results was performed to improve the accuracy of temperature predictions in reactor analysis models that homogenized the UO 2 -Mo fuel. Compared with UO 2 in a pressurized water reactor (PWR), the optimized UO 2 -Mo design increased the margin to fuel melt by 13–32% across the fuel cycle, but it requires the 235U enrichment to exceed 5% to match the cycle length of conventional UO 2 .

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Analysis of the EMCCD point-source response using x-rays

Electron Multiplying Charge Coupled Devices, EMCCD are used as x-ray detectors. The NSLS-II Soft Inelastic x-ray Scattering (SIX) beam line has two EMCCDs for x-ray detection in the spectrometer arm. The spectrometer with high resolving power disperses x-rays vertically. The x-ray vertical position on the sensor plane is related to its energy. This allows for very accurate x-ray energy measurements through x-ray coordinates. X-rays interact with silicon and create a number of electron-hole pairs proportional to the x-ray energy. Electrons drift and diffuse toward pixel gates and are collected there. The diffused electrons form a charge cloud distributed over several neighboring pixels. This charge sharing enables coordinate measurements with accuracy better than the pixel pitch. The charge distribution shape has to be taken into account to achieve ultimate accuracy in coordinate measurements. In this paper, we present a method of the charge distribution shape analysis and demonstrate its applications. The drift and diffusion of electrons ftom the point of generation to pixel gates results in the bell-shaped electron cloud usually approximated by Gaussian shape. The number of electrons collected under a pixel is proportional to the shape function integral. These electron packets get transferred to the sense node of the output amplifier. The transfer process could introduce distortions to the original charge distribution. For example, during transfers, electrons in the packet could be exposed to traps if they are present in the sensor. The trapping and later the release processes distort the apparent shape of the charge distribution. Therefore, deviations of the charge distribution shape ftom the originally symmetrical form can indicate the presence of trap centers in the sensor and can be used for sensor diagnostics.

47 OTHER INSTRUMENTATION↗

Assessment of RELAP5-3D Code with Molten Salt Heat Transfer Experiments

The accurate thermal-hydraulic assessment of thermal storage systems is crucial for addressing the integrity and performance of the storage system design. This is particularly important for thermal storage systems using molten salts as a heat transport and storage medium, as the chemical corrosion and erosion by the high-temperature molten salts can cause serious damage to the structural components of the storage system. To understand these effects for TerraPower’s Natrium® Demonstration Reactor design, the system thermal-hydraulic analysis code, RELAP5-3D, is utilized to analyze the thermal storage system. This study evaluates two heat transfer correlations implemented in the RELAP5-3D code—the Dittus-Boelter and Gnielinski correlations—using selected benchmark cases from two molten salt heat transfer experiments conducted at Xi’an Jiao Tong University and the German Aerospace Center. The Nusselt numbers calculated using the Gnielinski correlation agree with the experimental data within ±5% of the relative deviations at 300°C and 400°C. With the increasing Reynolds numbers, the Dittus-Boelter correlation underestimates the Nusselt numbers by up to 22% compared to the Gnielinski correlation. Comparison of the RELAP5-3D calculation results with experimental data at a high temperature of 550°C revealed that the thermo-chemical characteristics of solar salt at high temperatures surpassing the decomposition temperatures of nitrate salts reduce the accuracy of heat transfer model of the RELAP5-3D code. Specifically, it has been observed that the temperatures of the wall surface and liquid near the wall can surpass the decomposition temperatures, even though the bulk temperature of molten salts remains significantly lower than these decomposition temperatures. In summary, this study recommends the use of the Gnielinski correlation for the heat transfer analysis of molten salt, rather than the Dittus-Boelter correlation.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗

Experimental testing of additively manufactured embedded fiber optic smart devices for clean energy applications

Abstract An additively manufactured prototype smart device was created to investigate in-flow temperature distributions using embedded high-definition fiber optic sensors within a component for clean energy systems. The devices were created using Ultrasonic Additive Manufacturing to create a unique embedded pathway within a flow conditioner for the high-definition fiber optic sensors to be placed within. The fibers used allowed for temperature measurements to be taken every 0.65 mm along the fiber. The high-resolution fibers were thermally calibrated enable the 2D reconstruction of the temperature profile in the flow path of the structure. This is due to the temperature-related strain response of the material and strain measurements of the fibers. Hot airflow testing of these devices showed the ability to identify localized temperature differences in the flow. The observed strain response within the smart device had higher strain concentrations in the thicker web regions than in the thinner web regions. These higher strain regions resulted in higher uncertainties for the temperature responses. Further calibration is needed to improve the accuracy of the smart devices, specifically within the inner web structures of a flow straightening component. Further investigations of the devices within flow showed the temperature sensing to be independent of the effects of flow velocity. The devices were able to distinguish temperature differences within single and two-phase flow and showed local sensitivity to the temperature changes with the identification of hot and cold spots. The presented results showed the viability of the smart device for obtaining detailed temperature distributions using common industrial components. Eventually, the goal for these smart devices will be to withstand higher temperature and pressure environments such as those experienced in nuclear, fusion, and concentrated solar energy systems.

Donlan, Connor F. (ORCID:0000000223317882)↗

Field-theoretic analysis of hadronization using soft drop jet mass

One of the greatest challenges in quantum chromodynamics is understanding the hadronization mechanism, which is also crucial for carrying out precision physics with jet substructure. In this paper, we bring together recent advancements in our understanding of nonperturbative structure of the soft drop jet mass based on field theory, with precise perturbative calculations at next-to-next-to-leading logarithmic accuracy of its multidifferential variants. This allows for a model-independent analysis of power corrections associated with hadronization in a systematic manner. We test and calibrate hadronization models and their interplay with parton showers by comparing our universality predictions with various event generators for quark and gluon initiated jets in both lepton-lepton and hadron-hadron collisions. Our findings reveal that hadronization models perform better for quark jets relative to gluon jets. Our results provide a valuable toolbox for precision studies with the soft drop jet mass and pave the way for future analyses using real-world collider data. The stringent constraints derived in our framework are useful for improving the modeling of hadronization and its interplay with parton showers in next-generation event generators.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Freely jointed chain models with extensible links

We report analytical relations for the mechanical response of single polymer chains are valuable for modeling purposes, on both the molecular and the continuum scale. These relations can be obtained using statistical thermodynamics and an idealized single-chain model, such as the freely jointed chain model. To include bond stretching, the rigid links in the freely jointed chain model can be made extensible, but this almost always renders the model analytically intractable. Here, an asymptotically correct statistical thermodynamic theory is used to develop analytic approximations for the single-chain mechanical response of this model. The accuracy of these approximations is demonstrated using several link potential energy functions. This approach can be applied to other single-chain models, and to molecular stretching in general.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Dynamic Phasor Modeling of Various Multipulse Rectifiers and a VSI Fed by 18-Pulse Asymmetrical Autotransformer Rectifier Unit for Fast Transient Analysis

To study the fast modelling of dynamic characteristics caused by power electronic switching, the dynamic phasor (DP) theory based on the time-varying Fourier decomposition and frequency shifting is firstly applied to the modelling of various multipulse rectifiers. The DP model for symmetrical 12-pulse phase-shifting reactor rectifier unit (PSR-RU) is derived with model order reduction, relations between ac and dc terminals, and Taylor series expansion. The DP model of asymmetrical 18-pulse autotransformer rectifier unit (AT-RU) is proposed based on the switching functions expressed in the DP domain. Meanwhile, the DP model of voltage source inverter (VSI) fed by asymmetrical 18-pulse AT-RU is built with the harmonic state-space (HSS) equations. Under both balanced and unbalanced conditions, the good calculation accuracy and rapid simulation speed of the developed DP models are validated by the detailed time-domain (TD) simulation.

42 ENGINEERING↗

Polimaster PM1401K-3 Test Report

Extensive testing, including measurements on large quantities of special nuclear material, were taken to assess the functionality of the Polimaster PM1401K-3 to other detectors used in the field. The Polimaster can detect alpha, beta, gamma, and neutron radiation. However, efficiency and accuracy has been sacrificed to accomplish this. Compared to the IdentiFINDER 2, SRM, and FH-40 alpha probe, the Polimaster underperforms in almost all categories regarding radiation detection. The Helium-3 neutron detector in the Polimaster, due to a larger volume, outperformed the IdentiFINDER 2 for neutron measurements. The alpha search mode showed higher efficiency than the FH-40 alpha probe, but uses a more complicated process. The user interface is relatively easy to navigate, clear, and comparable to other detectors. However, the detector has added complications when trying to detect alpha or beta radiation or save a spectrum.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

A Comparison of GPU-Accelerated Multiphase CFD Solvers on the Polaris Supercomputer: Part 1

This report is in support of the Innovative and Novel Computational Impact on Theory and Experiment (INCITE) program sponsored by the U.S. Department of Energy (USDOE). With INCITE-level resources, one project, titled BubblyFlow, was granted computational resources for the 2025 calendar year on the Polaris supercomputer at the Argonne Leadership Computing Facility (ALCF). The project aims to conduct simulations to understand the fundamental characteristics of turbulent bubbly flow phenomena in nature. Staff at the ALCF and Argonne’s Computational Science division, along with collaborators at the City College of New York and University of Illinois at Chicago, helped a summer student to assess the accuracy and performance of two high performance computing (HPC) codes. Both codes, ImExLBM and FluTAS, are fundamentally different in their mathematical and numerical modeling. However, both may be used to solve the same physical problem. The collaboration sought to better understand the differences between both codes in terms of accuracy and efficiency. This would ultimately help the BubblyFlow project better utilize resources and establish a knowledge-base of code capabilities in future simulation campaigns. We compare ImExLBM and FluTAS, two high-performance multiphase computational fluid dynamics (CFD) solvers, in terms of physical fidelity, time-to-solution, and parallel efficiency. We validate ImExLBM (Implicit-Explicit Lattice Boltzmann Method) against a canonical benchmark and assess it’s performance relative to FluTAS (Fluid Transport Accelerated Solver), a well-established open-source CFD code.

97 MATHEMATICS AND COMPUTING↗