Search NASA⌕ Search

SEARCH · Search NASA

Results for “Computer experiments”

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 829 records · Page 46
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

A comparison of two Stokes ice sheet models applied to the Marine Ice Sheet Model Intercomparison Project for plan view models (MISMIP3d)

Here, we present a comparison of the numerics and simulation results for two "full" Stokes ice sheet models, FELIX-S (Leng et al. 2012) and Elmer/Ice. The models are applied to the Marine Ice Sheet Model Intercomparison Project for plan view models (MISMIP3d). For the diagnostic experiment (P75D) the two models give similar results (< 2 % difference with respect to along-flow velocities) when using identical geometries and computational meshes, which we interpret as an indication of inherent consistencies and similarities between the two models. For the standard (Stnd), P75S, and P75R prognostic experiments, we find that FELIX-S (Elmer/Ice) grounding lines are relatively more retreated (advanced), results that are consistent with minor differences observed in the diagnostic experiment results and that we show to be due to different choices in the implementation of basal boundary conditions in the two models. While we are not able to argue for the relative favorability of either implementation, we do show that these differences decrease with increasing horizontal (i.e., both along- and across-flow) grid resolution and that grounding-line positions for FELIX-S and Elmer/Ice converge to within the estimated truncation error for Elmer/Ice. Stokes model solutions are often treated as an accuracy metric in model intercomparison experiments, but computational cost may not always allow for the use of model resolution within the regime of asymptotic convergence. In this case, we propose that an alternative estimate for the uncertainty in the grounding-line position is the span of grounding-line positions predicted by multiple Stokes models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

An automatic performance model-based scheduling tool for coupled climate system models

The prediction of the climate system is highly depended on the efficient integration of observations and simulations of the Earth. This is regarded as a canonical example of the cyber-physical system. The climate system model, the simulation engine in this cyber-physical system, is one of most challenging applications in scientific computing. It utilizes the multi-physics simulation that couples multiple components, conducts decadal to millennium simulations, and has long been an important application on supercomputers. However, current climate system models suffer from the inefficient task scheduling methods resulting in an intolerable simulation time. Take the Community Earth System Model (CESM), the most widely used climate system model, as an example, one major reason that CESM suffers from bad performances is the huge overhead to rationally distribute processes among the coupled heterogeneous components. According to the report of NCAR, every percent improvement in CESM performance frees up to the equivalent of $250,000 in computing resources in their scientific experiments. To address such challenge, our paper firstly constructs a lightweight and accurate performance model for effectively capturing and predicting the heterogeneous time-to-solution performance of end-to-end CESM components with a given simulation configuration. Then, based on the performance model, we further propose an efficient scheduling strategy based on rectangular packing method to determine the best process layout among different components, and the process numbers assigned to each component. Our evaluations show that we can achieve 58% average run time reductions on CESM comparing to the widely used sequential process layout for a scale of 144 to 480 cores on typical CPU clusters. And we can save 4 million CPU hours when we conduct one standard scientific experiment (a 2870-year simulation), which equals to save $40,089 with a charge of $0.01 per CPU hour. Meanwhile, 26% extra performance improvements could also be gained in our methods comparing to the heuristic branch and bound algorithm with the guidance of the known curve-fitting performance model.

58 GEOSCIENCES↗
DOE OSTI2023

Bentonite buffer under high temperature: laboratory experiments and coupled process modeling

Abstract. Bentonite buffer in a geological repository will be simultaneously heated from decaying radioactive waste and hydrated from the surrounding host rock, triggering complex and coupled THMC (thermal–hydrological–mechanical–chemical) processes. Understanding the THMC behavior of bentonite-based engineered barrier system (EBS) is key to the evaluation and prediction of its long-term performance. Studies on the THMC process have been focused on conditions under 100 ∘C, as most design concepts impose a thermal limit of 100 ∘C in bentonite. Recently, studies under high-temperature conditions have been conducted to evaluate the possibility of raising the thermal limit and expanding the data/knowledge base to increase the confidence level. In this abstract, we present a series of bench-scale laboratory experiments at high temperatures (up to 200 ∘C) and the corresponding modeling work. Two sets of column tests were conducted, and each set consisted of two test columns: a control column undergoing only hydration (non-heated) and an experiment column experiencing both heating and hydration (heated). During the experiment, frequent X-ray computed tomography (CT) images were collected to provide a 3D visualization of the density distribution and present the spatiotemporal evolution of (1) hydration/dehydration, (2) clay swelling/shrinkage, (3) displacement, and (4) mineral precipitation. The two sets of tests differ with respect to several experimental conditions, such as bentonite type, compacted density and water content, water chemistry, and hydration pressure, but the important difference is that the first set used bentonite powder with a dry density of 1.28 g cm−3, whereas the second set used granulated bentonite (mixture of pellets and powder) with a dry density of 1.45–1.5 g cm−3. In both sets of experiments, a comprehensive post-dismantling characterization of bentonite samples was carried out after the column tests had been running for 1.5 years. Comparing non-heated and heated columns, the temperature gradient led to lower degree of homogenization of bentonite after bentonite became fully saturated; comparing the first and second sets, granulated and powdered bentonite exhibited drastically different hydration behavior. A THM model with a 2D axisymmetric grid system was used to interpret the data from the first set of tests. The model considers the combined impact of saturation, fluid pressure, and porosity change due to swelling/compression on the spatiotemporal distribution of bulk density and movement of the thermocouple modules. Observations from the tests help us understand the early perturbation of bentonite buffer under high temperature, and data from these tests improve the calibration of key constitutive hydrological and mechanical models and, therefore, enhance the modeling capability with respect to calculating the long-term evolution of bentonite buffer.

Zheng, Liange↗
DOE OSTI2017

Thermodynamics and Reaction Mechanisms of Decomposition of the Simplest Protonated Tripeptide, Triglycine: A Guided Ion Beam and Computational Study

We present a thorough characterization of fragmentations observed in threshold collision-induced dissociation (TCID) experiments of protonated triglycine (H{sup +}GGG) with Xe using a guided ion beam tandem mass spectrometer (GIBMS). Kinetic energy-dependent cross-sections for 10 ionic products are observed and analyzed to provide 0 K barriers for six primary products: [b{sub 2}]{sup +}, [y{sub 1} + 2H]{sup +}, [b{sub 3}]{sup +}, CO loss, [y{sub 2} + 2H]{sup +}, and [a{sub 1}]{sup +}; three secondary products: [a{sub 2}]{sup +}, [a{sub 3}]{sup +}, and [y{sub 2} + 2H – CO]{sup +}; and two tertiary products: high energy [y{sub 1} + 2H]{sup +} and [a{sub 2} – CO]{sup +} after accounting for multiple ion-molecule collisions, internal energy of reactant ions, unimolecular decay rates, competition between channels, and sequential dissociations. Relaxed potential energy surface scans performed at the B3LYP-D3/6-311+G(d,p) level of theory are used to identify transition states (TSs) and intermediates of the six primary and one secondary products. Geometry optimizations and single point energy calculations were performed at several levels of theory. These theoretical energies are compared with experimental energies and are found to give reasonably good agreement, in particular for the M06-2X level of theory. This good agreement between experiment and theory validates the reaction mechanisms explored computationally here and elsewhere and allows identification of the product structures formed at threshold energies. The present work presents the first measurement of absolute experimental threshold energies of important sequence ions and non-sequence ions: [y{sub 1} + 2H]{sup +}, [b{sub 3}]{sup +}, CO loss, [a{sub 1}]{sup +}, and [a{sub 3}]{sup +}, and refines those for [b{sub 2}]{sup +} and [y{sub 2} + 2H]{sup +} previously measured. .

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

A 1D–0D–3D coupled model for simulating blood flow and transport processes in breast tissue

In this work, we present mixed dimensional models for simulating blood flow and transport processes in breast tissue and the vascular tree supplying it. These processes are considered, to start from the aortic inlet to the capillaries and tissue of the breast. Large variations in biophysical properties and flow conditions exist in this system necessitating the use of different flow models for different geometries and flow regimes. In total, we consider four different model types. First, a system of 1D nonlinear hyperbolic partial differential equations (PDEs) is considered to simulate blood flow in larger arteries with highly elastic vessel walls. Second, we assign 1D linearized hyperbolic PDEs to model the smaller arteries with stiffer vessel walls. The third model type consists of ODE systems (0D models). It is used to model the arterioles and peripheral circulation. Finally, homogenized 3D porous media models are considered to simulate flow and transport in capillaries and tissue within the breast volume. Sink terms are used to account for the influence of the venous and lymphatic systems. Combining the four model types, we obtain two different 1D–0D–3D coupled models for simulating blood flow and transport processes: The first model results in a fully coupled 1D–0D–3D model covering the complete path from the aorta to the breast combining a generic arterial network with a patient specific breast network and geometry. The second model is a reduced one based on the separation of the generic and patient specific parts. The information from a calibrated fully coupled model is used as inflow condition for the patient specific sub-model allowing a significant computational cost reduction. Several numerical experiments are conducted to calibrate the generic model parameters and to demonstrate realistic flow simulations compared to existing data on blood flow in the human breast and vascular system. Moreover, we use two different breast vasculature and tissue data sets to illustrate the robustness of our reduced sub-model approach.

multiscale blood flow models↗
DOE OSTI2019

Scalable Algorithms for MPI Intergroup Allgather and Allgatherv

MPI intergroup collective communication defines message transfer patterns between two disjoint groups of MPI processes. Such patterns occur in coupled applications, and in modern scientific application workflows, mostly with large data sizes. However, current implementations in MPI production libraries adopt the “root gathering algorithm”, which does not achieve optimal communication transfer time. In this paper, we propose algorithms for the intergroup Allgather and Allgatherv communication operations under single-port communication constraints. Here we implement the new algorithms using MPI point-to-point and standard intra-communicator collective communication functions. We evaluate their performance on the Cori supercomputer at NERSC. Using message sizes per compute node ranging from 64KBytes to 8MBytes, our experiments show significant performance improvements of up to 23.67 times on 256 compute nodes compared with the implementations of production MPI libraries.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Metal phosphides as potential thermoelectric materials

We report there still exists a crucial need for new thermoelectric materials to efficiently recover waste heat as electrical energy. Although metal phosphides are stable and can exhibit excellent electronic properties, they have traditionally been overlooked as thermoelectrics due to expectations of displaying high thermal conductivity. Based on high-throughput computational screening of the electronic properties of over 48000 inorganic compounds, we find that several metal phosphides offer considerable promise as thermoelectric materials, with excellent potential electronic properties (e.g. due to multiple valley degeneracy). In addition to the electronic band structure, the phonon dispersion curves of various metal phosphides were computed indicating low-frequency acoustic modes that could lead to low thermal conductivity. Several metal phosphides exhibit promising thermoelectric properties. The computed electronic and thermal properties were compared to experiments to test the reliability of the calculations indicating that the predicted thermoelectric properties are semi-quantitative. As a complete experimental study of the thermoelectric properties in MPs, cubic-NiP2 was synthesized and the low predicted lattice thermal conductivity (~1.2 W m -1 K -1 at 700 K) was confirmed. The computed Seebeck coefficient is in agreement with experiments over a range of temperatures and the phononic dispersion curve of c-NiP 2 is consistent with the experimental heat capacity. Finally, the predicted high thermoelectric performance in several metal phosphides and the low thermal conductivity measured in NiP 2 encourage further investigations of thermoelectric properties of metal phosphides.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Denoising and Extension of Response Functions in the Time Domain

Response functions of quantum systems, such as electron Green’s functions, magnetic, or charge susceptibilities, describe the response of a system to an external perturbation. They are the central objects of interest in field theories and quantum computing and measured directly in experiment. Further, response functions are intrinsically causal. In equilibrium and steady-state systems, they correspond to a positive spectral function in the frequency domain. Since response functions define an inner product on a Hilbert space and thereby induce a positive definite function, the properties of this function can be used to reduce noise in measured data and, in equilibrium and steady state, to construct positive definite extensions for data known on finite time intervals, which are then guaranteed to correspond to positive spectra.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Soft pions and transport near the chiral critical point

During the expansion of a heavy ion collision, the system passes close to the $\textit{O}$(4) critical point of QCD, and thus the fluctuations of the order parameter $(σ, \overrightarrow{π}$) are expected to be enhanced. In this work, our goal is to compute how these enhanced fluctuations modify the transport coefficients of QCD near the pseudocritical point. We also make a phenomenological estimate for how chiral fluctuations could effect the momentum spectrum of soft pions. We first formulate the appropriate stochastic hydrodynamic equations close to the $\textit{O}$(4) critical point. Then, working in mean field, we determine the correlation functions of the stress tensor and the currents which result from this stochastic real-time theory, and use these correlation functions to determine the scaling behavior of the transport coefficients. The hydrodynamic theory also describes the propagation of pion waves, fixing the scaling behavior of the dispersion curve of soft pions. We present scaling functions for the shear viscosity and the charge conductivities near the pseudocritical point and estimate the absolute magnitude of the critical fluctuations to these parameters and the bulk viscosity. Using the calculated pion dispersion curve, we estimate the expected critical enhancement of soft pion yields, and this estimate provides a plausible explanation for the excess seen in experiment relative to ordinary hydrodynamic computations. Our results motivate further phenomenological and numerical work on the implications of chiral symmetry on real-time properties of thermal QCD near the pseudocritical point.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2016

High-order shock-fitted detonation propagation in high explosives

In this paper, a highly accurate numerical shock and material interface fitting scheme composed of fifth-order spatial and third- or fifth-order temporal discretizations is applied to the two-dimensional reactive Euler equations in both slab and axisymmetric geometries. High rates of convergence are not typically possible with shock-capturing methods as the Taylor series analysis breaks down in the vicinity of discontinuities. Furthermore, for typical high explosive (HE) simulations, the effects of material interfaces at the charge boundary can also cause significant computational errors. Fitting a computational boundary to both the shock front and material interface (i.e. streamline) alleviates the computational errors associated with captured shocks and thus opens up the possibility of high rates of convergence for multi-dimensional shock and detonation flows. Several verification tests, including a Sedov blast wave, a Zel'dovich–von Neumann–Döring (ZND) detonation wave, and Taylor–Maccoll supersonic flow over a cone, are utilized to demonstrate high rates of convergence to nontrivial shock and reaction flows. Comparisons to previously published shock-capturing multi-dimensional detonations in a polytropic fluid with a constant adiabatic exponent (PF-CAE) are made, demonstrating significantly lower computational error for the present shock and material interface fitting method. For an error on the order of 10 m / s , which is similar to that observed in experiments, shock-fitting offers a computational savings on the order of 1000. In addition, the behavior of the detonation phase speed is examined for several slab widths to evaluate the detonation performance of PBX 9501 while utilizing the Wescott–Stewart–Davis (WSD) model, which is commonly used in HE modeling. It is found that the thickness effect curve resulting from this equation of state and reaction model using published values is dramatically more steep than observed in recent experiments. Utilizing the present fitting strategy, in conjunction with a nonlinear optimizer, a new set of reaction rate parameters improves the correlation of the model to experimental results. In conclusion, this new model is tested against two dimensional slabs as a validation test.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2020

Phase Equilibria and Thermochemistry of Advanced Fuels: Modeling Burnup Behavior (Final Report)

Achieving the goal of developing advanced fuel concepts that meet the DOE objectives of being robust, demonstrating high performance, and are more tolerant of accident conditions than current fuel systems will require a thorough understanding of the thermophysical and thermochemical properties of the constituent materials. Non-oxide fuel systems are being explored under the Advanced Fuels Program that hold significant promise for improved performance and accident tolerance, including the uranium silicide-based system considered in the current work. Prospective cladding materials currently considered that contribute to improved accident tolerance include silicon carbide composites and ferritic alloys (Fe-Cr-Al base compositions). Thus, the effort developed thermochemical models and values, supported with targeted experiments, to evaluate the ferritic alloy and silicon carbide composite cladding systems in contrast to current zirconium alloy cladding. The developed detailed understanding will serve to aid in relatively early screening of candidate systems to avoid wasted effort, guide development of new fuel forms, and to provide a basis for predicting and modeling fuel performance. Major deliverables for the project included: Thermochemical assessment and models of phases in the U-Si and U-Si-N systems; thermochemical evaluation supported by experimental measurements of fuel-cladding interactions of silicide fuel with baseline zirconium, ceramic composite, and ferritic alloy cladding; thermochemical assessment and models of phases supported by experimental measurements for silicide fuel with key fission products provided in a dataset and reported in refereed publications. Within the project a significantly refined U-Si phase diagram was developed and reported that now includes homogeneity ranges for key phases, such at the U 3 Si 2 proposed fuel phase, and settles issues with regard to uncertainty in the stability of some phases. Computational efforts together with key experiments has determined phase formation in interactions between U 3 Si 2 and Zircaloy-4 cladding material, a ferritic FeCrAlY alloy of interest as an advanced cladding material, and silicon carbide, also of interest as a fiber-reinforced composite cladding. As expected, very significant reactions occur between U3Si2 and Zircaloy-4, with much less interaction at higher temperatures for the ferritic alloy, and finally interactions with SiC only in the region of contact. A potentially major issue is the stability of U 3 Si 2 fuel that has undergone significant burnup. The result is the loss of the uranium metal, liberating silicon, and the formation of concomitant fission product elements that either dissolve in the U 3 Si 2 phase or form independent, and possibly silicide phases. A combination of experimental determinations of phase formation of U 3 Si 2 reacted with representative fission products yttrium, gadolinium, cerium, zirconium, and molybdenum and first principles calculations has helped understand the fuel chemistry. The behavior of the U 3 Si 2 phase and the partitioning of silicon to possible fission product phases with burnup was thus determined, with significant dissolution in U 3 Si 2 of cerium, gadolinium, zirconium, and plutonium predicted along with independent phase formation of a U-Mo-Si ternary phase, yttrium silicide, and elemental selenium. It can be concluded that at significant burnup there will be a very minor amount of the U 3 Si 2 fuel phase that will decompose to a lower silicide or a uranium alloy as silicon preferentially forms a secondary phase.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Computation of conventional and alternative jet fuel sensitivity to lean blowout

Large Eddy Simulations (LES) are performed to compute the sensitivity of a conventional (A-2) and an alternate bio-jet (C-1) fuel to Lean Blowout (LBO). A realistic aviation gas turbine engine combustor configuration is considered. Reliable experimental LBO data and OH* chemiluminescence data for the conventional and alternate jet fuel in the combustor configuration have recently become available. The present work utilizes a highly automated, on-the-fly meshing strategy, along with adaptive mesh refinement, to demonstrate the feasibility of capturing the realistic combustion processes. A Lagrangian framework, with initial conditions specified using measurements of spray statistics, is used to model the fuel spray. Newly developed compact reaction mechanisms based on fuel surrogates are validated for the A-2 and the C-1 fuels. The compact reaction mechanisms are implemented using a detailed finite rate chemistry solver. Spray statistics computed by the present LES simulations compare well with available measurements at stable flame conditions near the lean blowout limit. The computed shape of the stable flame as represented by line integrated OH concentrations compares well with the experimental OH* chemiluminescence data. Lean blowout is reached by gradually decreasing the fuel flow rate in the computations, similar to that in the experiments. The results of the LES simulations effectively capture the fuel composition effects and estimate the sensitivity of the LBO limits to the fuel type. The computed trends in LBO limits agree within engineering accuracy with the experimental results for conventional and alternative aviation fuels. The methodology for predicting the fuel composition effects on the lean blowout limits in a fully resolved realistic, complex combustor is established for the first time.

Adaptive mesh refinement↗
DOE OSTI2023

Code Coverage Status of ARC Code-DIF3D

The Argonne Reactor Code (ARC) software system supports users in their fast reactor design goals by providing neutronic, thermal-hydraulic, and structural analysis capabilities. DIF3D plays a pivotal role in the ARC system as the primary homogenized assembly neutronic calculation methodology for fast reactor problems. Over its 40 years history, ARC software usage with DIF3D has been applied to numerous fast and thermal spectrum reactor analysis projects with good to excellent comparison against experiments. With continued improvement of computation resources, many of the geometry modeling capabilities in DIF3D that were primarily used in low order schemes are not really needed anymore. Today, the diffusion and transport capabilities of DIF3D-VARIANT are primarily used in the reactor design process with some scattered usage of DIF3D-FD and DIF3D-Nodal. In recent work, the DIF3D software verification was completed for DIF3D-FD and DIF3D-VARIANT on the geometry options used in the Versatile Test Reactor project. While we can be confident that these capabilities of DIF3D are well used and thus trusted, it does not demonstrate that all possible input options of DIF3D are actually working, but just those that were tested as part of VTR are and that they are correct. Thus, the purpose of the present work is to identify a set of test problems for DIF3D and assess the code coverage of DIF3D for those test problems. The goal is to document what parts of the existing DIF3D code are touched by the set of test problems and which are not. Because the verification work done on DIF3D-VARIANT and DIF3D-FD was focused on the most common uses of DIF3D for fast reactor analysis, the code coverage assessment of those capabilities is the highest priority. This will ensure that nothing is being missed by the existing verification test problems that DIF3D relies upon. The DIF3D-Nodal capability will also be inspected for code coverage as part of this work to further ensure that regular regression testing of DIF3D will trap any likely errors the end user might experience with the DIF3D software. The code coverage analysis of DIF3D was performed with the Code Coverage Tool of the Intel Fortran compiler which requires modifications to the compilation of DIF3D. The detailed coverage tables are given for each submodule of DIF3D separately, and for the submodules which are primarily developed for DIF3D, most of the source files could be at least partially touched. Most of the uncovered parts/files could be easily ignored, because they are either for error message and debugging output or obviously not needed by DIF3D. Out of the entire source codes of DIF3D, only a few uncovered modules deserve further investigation.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Code Coverage Status of the ARC Code DIF3D

The Argonne Reactor Code (ARC) software system supports users in their fast reactor design goals by providing neutronic, thermal-hydraulic, and structural analysis capabilities. DIF3D plays a pivotal role in the ARC system as the primary homogenized assembly neutronic calculation methodology for fast reactor problems. Over its 40 years history, ARC software usage with DIF3D has been applied to numerous fast and thermal spectrum reactor analysis projects with good to excellent comparison against experiments. With continued improvement of computation resources, many of the geometry modeling capabilities in DIF3D that were primarily used in low order schemes are not really needed anymore. Today, the diffusion and transport capabilities of DIF3D-VARIANT are primarily used in the reactor design process with some scattered usage of DIF3D-FD and DIF3D-Nodal. In recent work, the DIF3D software verification was completed for DIF3D-FD and DIF3D-VARIANT on the geometry options used in the Versatile Test Reactor project. While we can be confident that these capabilities of DIF3D are well used and thus trusted, it does not demonstrate that all possible input options of DIF3D are actually working, but just those that were tested as part of VTR are and that they are correct. Thus, the purpose of the present work is to identify a set of test problems for DIF3D and assess the code coverage of DIF3D for those test problems. The goal is to document what parts of the existing DIF3D code are touched by the set of test problems and which are not. Because the verification work done on DIF3D-VARIANT and DIF3D-FD was focused on the most common uses of DIF3D for fast reactor analysis, the code coverage assessment of those capabilities is the highest priority. This will ensure that nothing is being missed by the existing verification test problems that DIF3D relies upon. The DIF3D-Nodal capability will also be inspected for code coverage as part of this work to further ensure that regular regression testing of DIF3D will trap any likely errors the end user might experience with the DIF3D software. The code coverage analysis of DIF3D was performed with the Code Coverage Tool of the Intel Fortran compiler which requires modifications to the compilation of DIF3D. The detailed coverage tables are given for each submodule of DIF3D separately, and for the submodules which are primarily developed for DIF3D, most of the source files could be at least partially touched. Most of the uncovered parts/files could be easily ignored, because they are either for error message and debugging output or obviously not needed by DIF3D. Out of the entire source codes of DIF3D, only a few uncovered modules deserve further investigation.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2017

XRootD popularity on hadoop clusters

Performance data and metadata of the computing operations at the CMS experiment are collected through a distributed monitoring infrastructure, currently relying on a traditional Oracle database system. This paper shows how to harness Big Data architectures in order to improve the throughput and the efficiency of such monitoring. A large set of operational data - user activities, job submissions, resources, file transfers, site efficiencies, software releases, network traffic, machine logs - is being injected into a readily available Hadoop cluster, via several data streamers. The collected metadata is further organized running fast arbitrary queries; this offers the ability to test several Map&Reduce-based frameworks and measure the system speed-up when compared to the original database infrastructure. By leveraging a quality Hadoop data store and enabling an analytics framework on top, it is possible to design a mining platform to predict dataset popularity and discover patterns and correlations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Federated Access from DOE Labs to Distributed Storage in the EIC Era of Computing

The Electron Ion Collider (EIC) collaboration and future experiment is a unique scientific ecosystem within Nuclear Physics as the experiment starts right off as a crosscollaboration between Brookhaven National Lab (BNL) & Jefferson Lab (JLab). As a result, this muti-lab computing model tries at best to provide services accessible from anywhere by anyone who is part of the collaboration. While the computing model for the EIC is not finalized, it is anticipated that the computational and storage resources will be made accessible to a wide range of collaborators across the world. The use of federated ID seems to be a critical element to the strategy of providing such services, allowing seamless access to each lab site computing resources. However, providing Federated access to a Federated storage is not a trivial matter and has its share of technical challenges. In this contribution, we focus on the steps we took towards the deployment of a distributed object storage system that integrates with Amazon S3 and Federated ID. We will first cover for and explain the first stage storage solutions provided to the EIC during the detector design phase. Our initial test deployment consisted of Lustre storage using MinIO, hence providing an S3 interface. High Availability load balancers were added later to provide the initial scalability it lacked. Performance of that system will be shown. While this embryonic solution worked well, it had many limitations. Looking ahead, the Ceph object storage is considered a top-of-the-line solution in the storage community - since the Ceph Object Gateway is compatible with the Amazon S3 API out of the box, our next phase will use a native S3 storage. Our Ceph deployment will consist of erasure coded storage nodes to maximize storage potential along with multiple Ceph Object Gateways for redundant access. We will compare performance of our next stage implementations. Finally, we will present how to leverage OpenID Connect with the Ceph Object Gateway’s to enable Federated ID access. We hope this contribution will serve the community needs as we move forward with cross-lab collaborations and the need for Federated ID access to distributed compute facilities.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

On the integration of molecular dynamics, data science, and experiments for studying solvent effects on catalysis

Computational workflows that combine molecular dynamics (MD) simulations and emerging data-centric (DC) methods can accelerate the screening and analysis of solvent systems experimentally and computationally. Here, MD simulations provide atomic positions and velocities of reactant, solvent, and catalyst materials that can be manipulated into data representations that in turn can be used by DC techniques to conduct predictive modeling, feature extraction, and experimental design. For liquid-phase catalytic applications, emerging DC techniques such as Convolutional and Graph Neural Networks (CNN/GNN), Topological Data Analysis (TDA), and Active Learning (AL) can leverage MD and experimental data to quickly predict solvent effects on reaction outcomes. For instance, in recent studies, 3D solvent environments obtained with MD have been exploited by CNNs to predict experimental reaction rates for homogeneous acid-catalyzed lignocellulosic processes. In this perspective, we discuss basic principles of DC methods and how these can be combined with MD to enable high-throughput screening of solvent selection for diverse catalysis applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗