Search NASA⌕ Search

SEARCH · Search NASA

Results for “Graphical model selection”

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 145 records · Page 8

Water Technoeconomic Assessment Pipe-Parity Platform (WaterTAP3)

The Water Technoeconomic Assessment Pipe-Parity Platform (WaterTAP3) was developed under the National Alliance for Water Innovation (NAWI) to facilitate consistent technoeconomic assessments of desalination treatment trains. The WaterTAP3 is an analytically robust modeling tool that can be used to evaluate water technology cost, energy, environmental, and resiliency tradeoffs across different water sources, sectors, and scales. The model simulates steady-state water treatment train performance and costs including flow and constituent mass balance across unit processes, based on source water conditions, configurations of treatment technologies, and system-level techno-economic assumptions. Users can build a new treatment train by connecting any number of unit processes, specific for their context and system, or selecting a train from the treatment train library. The model contains various technical and cost parameter options for a range of treatment processes and a library of influent water quality characteristics for a variety of source waters and case studies. Users can customize water quality parameters to evaluate the technology performance in their context. The model can be set up for different assessment needs including simulation, optimization, and uncertainty and sensitivity analyses. The results from WaterTAP3 can help identify trade-offs among the different system performance metrics, with insight on how particular technologies or systems promote pipe-parity. The flexibility and comprehensive scope of the tool makes it a promising solution to industry-wide water technoeconomic evaluations, leading to more informed water investment decisions and technologies. As a user-friendly, open-source platform, WaterTAP3 can be used by industry, academia, policymakers, planners, and those with or without extensive analytical experience. A publicly available graphical user interface is currently under development.

Miara, Ariel↗

XCT analysis of the influence of melt strategies on defect population in Ti–6Al–4V components manufactured by Selective Electron Beam Melting

Selective Electron Beam Melting (SEBM) is a promising powder bed Additive Manufacturing technique for near-net-shape manufacture of high-value titanium components. However without post-manufacture HIPing the fatigue life of SEBM parts is currently dominated by the presence of porosity. In this study, the size, volume fraction, and spatial distribution of the pores in model samples have been characterised in 3D, using X-ray Computed Tomography, and correlated to the process variables. The average volume fraction of the pores (< 0.2%) was measured to be lower than that usually observed in competing processes, such as selective laser melting, but a strong relationship was found with the different beam strategies used to contour, and infill by hatching, a part section. The majority of pores were found to be small spherical gas pores, concentrated in the infill hatched region; this was attributed to the lower energy density and less focused beam used in the infill strategy allowing less opportunity for gas bubbles to escape the melt pool. Overall, increasing the energy density or focus of the beam was found to correlate strongly to a reduction in the level of gas porosity. Rarer irregular shaped pores were mostly located in the contour region and have been attributed to a lack of fusion between powder particles. - Graphical abstract: Display Omitted - Highlights: • Vast majority of defects detected were small spherical gas pores. • Gas bubbles trapped in the powder granules expand and coalesce in the melt pool. • Pores have been shown not to be randomly distributed. • Larger and deeper melt pools give more opportunity for gas to escape. • Minor changes to melt strategy result in significant reductions in pore population.

36 MATERIALS SCIENCE↗

Measurement of the cross section of top quark-antiquark pair production in association with a W boson in proton-proton collisions at $ \sqrt{s} $ = 13 TeV

The production of a top quark-antiquark pair in association with a W boson ($ \textrm{t}\overline{\textrm{t}}\textrm{W} $) is measured in proton-proton collisions at a center-of-mass energy of 13 TeV. The analyzed data was recorded by the CMS experiment at the CERN LHC and corresponds to an integrated luminosity of 138 fb$^{−1}$. Events with two or three leptons (electrons and muons) and additional jets are selected. In events with two leptons, a multiclass neural network is used to distinguish between the signal and background processes. Events with three leptons are categorized based on the number of jets and of jets originating from b quark hadronization, and the lepton charges. The inclusive $ \textrm{t}\overline{\textrm{t}}\textrm{W} $ production cross section in the full phase space is measured to be 868 ± 40(stat) ± 51(syst) fb. The $ \textrm{t}\overline{\textrm{t}}\textrm{W} ^{+}$ and $ \textrm{t}\overline{\textrm{t}}\textrm{W} ^{−}$ cross sections are also measured as 553 ± 30(stat) ± 30(syst) and 343 ± 26(stat) ± 25(syst) fb, respectively, and the corresponding ratio of the two cross sections is found to be $ 1.61\pm 0.15{\left(\textrm{stat}\right)}_{-0.05}^{+0.07}\left(\textrm{syst}\right) $. The measured cross sections are larger than but consistent with the standard model predictions within two standard deviations, and represent the most precise measurement of these cross sections to date.[graphic not available: see fulltext]

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Synthesis and application of different phthalocyanine molecular sieve catalyst for oxidative desulfurization

M{sub 2}(PcAN){sub 2} (M=Fe, Co, Ni, Cu, Zn and Mn) anchored onto W-HZSM-5 (M{sub 2}(PcAN){sub 2}–W-HZSM-5) or the M{sub 2}(PcTN){sub 2} doping W-HZSM-5 (M{sub 2}(PcTN){sub 2}/W-HZSM-5) were prepared and their catalytic performances were tested for oxidative desulfurization in the presence of oxygen. Thiophene (T), benzothiophene (BT), and dibenzothiophene (DBT) were considered as sulfur compounds. Among zeolite-based catalysts, the Cu{sub 2}(PcAN){sub 2}–W-HZSM-5 and Cu{sub 2}(PcTN){sub 2}/W-HZSM-5 showed superior desulfurization performance and the activity of selectivity followed the order: T>BT>DBT. The effects of phthalocyanine concentration were studied by UV–Vis and calcination temperature was obtained by TG-DSC for Cu{sub 2}(PcTN){sub 2}/W-HZSM-5. Catalysts were characterized by EA, IR, XRD, SEM, TEM, ICP, and N{sub 2} adsorption. Reaction time, temperature and the amount of catalyst were investigated as the important parameters for optimization of the reaction. Furthermore, a possible process of oxidative desulfurization and the reaction products were proposed. - Graphical abstract: The ODS reaction schematic shows the reaction mechanism of ultra-deep desulfurization. The sulfur compounds are oxidized to their corresponding sulfoxides or sulfones through the use of oxygen and catalysts. The reaction process of ultra-deep desulfurization. - Highlights: • A kind of novel catalyst for deep desulfurization was synthesized. • Cu{sub 2}(PcAN){sub 2}–W-HZSM-5 exhibits excellent catalytic performance for desulfurization. • The reaction conditions that affect desulfurization efficiency are investigated. • The reaction process of model sulfur compounds is proposed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

pop-cosmos : redshifts and physical properties of KiDS-1000 galaxies

ABSTRACT Principled Bayesian inference of galaxy properties has not previously been performed for wide-area weak-lensing surveys with millions of sources. We address this gap by applying the pop-cosmos generative model to perform spectral energy distribution (SED) fitting for 4 million KiDS (Kilo-Degree Survey)-1000 galaxies. Calibrated on deep COSMOS2020 photometric data, pop-cosmos specifies a physically motivated prior over the galaxy population up to $z \simeq 6$ in stellar population synthesis (SPS) parameter space. Using the Speculator SPS emulator with GPU (graphics processing unit)-accelerated Markov Chain Monte Carlo sampling, we perform full posterior inference at 8.2 GPU seconds per galaxy, obtaining joint constraints on galaxy redshifts and physical properties. We validate photometric redshifts against $\sim \!185\,\!000$ KiDS galaxies cross-matched to Dark Energy Spectroscopic Instrument Data Release 1 spectroscopic samples, achieving low bias ($2\times 10^{-3}$), scatter ($\sigma _{\mathrm{MAD}}=0.03$), and outlier fraction (3.2 per cent) for the Bright Galaxy Survey, with comparable performance (bias $3\times 10^{-2}$, $\sigma _{\mathrm{MAD}}=0.05$, 1.0 per cent outliers) for luminous red galaxies (LRGs). Within the LRG sample, we identify massive, dusty, star-forming contaminants at $z \simeq 0.4$ satisfying standard colour selections for quenched populations. We infer trends in stellar mass, star formation, metallicity, and dust across five tomographic redshift bins consistent with established scaling relations. Using specific star formation rate constraints, we identify $\sim$7 per cent of KiDS-1000 galaxies as quenched, versus 37 per cent implied by conservative colour cuts. This enables the construction of weak-lensing samples defined by physical properties while mitigating intrinsic alignment systematics and preserving statistical power. Our analysis validates pop-cosmos out of sample, establishing it as a scalable approach for galaxy evolution and cosmological analyses with photometric surveys.

Halder, Anik [Institute of Astronomy and Kavli Ins↗

VENTSAR V3.0: A Python GUI for Estimating Contaminant Concentrations on or Near Buildings Due to Building Effects and Plume Rise and For Calculating Inhalation and Plume Shine Doses

VENTSAR: Originated as VENTAX (Smith and Weber 1983) on the IBM Mainframe at the Savannah River Site (SRS) as a Fortran program. Updated to VENTSAR XL V1.0 (Simpkins 1997) as a spreadsheet version using macros created in Microsoft Excel. Later updated to VENTSAR XL V2.0 (Dixon 2018) as a spreadsheet executing a modified version of the original VENTAX Fortran code. Estimates contaminant concentrations on or near a building from a release at a nearby location. Calculates concentrations for a given meteorological exceedance probability or for a given stability and wind speed combination. Can model a single building with or without a penthouse on top or a ground location from either a stack or ground release. Plume rise can be considered. Contaminant releases can be chemical or radioactive with downwind concentrations determined at user-specified distances. Wind passing over and around buildings creates a complicated dispersion pattern. Air-intake vents may be located on building roofs or near the ground downwind of a release source. An estimation of pollutant concentrations on or near a structure is important in determining expected pollutant levels. Meteorological data are selected based on a specific area of the site. Fortran-based VENTAX able to make fast, complex mathematical calculations. Not user-friendly VENTSAR XL V1.0 no longer supported due to obsolete Excel macros. VENTSAR XL V2.0 an interim solution using an Excel spreadsheet to interface with the modified VENTAX Fortran code. Requires user to have Microsoft Excel. Not an intuitive interface for someone unfamiliar with VENTSAR. Does not perform dose calculations. VENTSAR V3.0 goal to combine the strengths of previous versions of VENTSAR into a single, powerful yet user-friendly program with a Graphical User Interface (GUI). Not dependent upon a specific program or plug-ins. Self-contained package to be used on any computer running Microsoft Windows. User-friendly, intuitive GUI created in Python for parameter input. Executes same modified VENTAX Fortran code as VENTSAR XL V2.0. Outputs formatted text file of completed calculations for easy review and dissemination. Performs inhalation and plume shine dose calculations for up to 11 user-selected nuclides from a nuclide dose factor library of almost 500 nuclides. 12 test cases were created for verification of V1.0 and V2.0. Same test cases run in V3.0 to verify correct performance. V3.0 produced similar results to V1.0. Confirmed GUI did not alter calculations in any way. Comparisons of the test case results confirm the V3.0 GUI does not influence the VENTSAR calculations. Simply passes same input parameters to VENTAX Fortran code for execution. Provides end user an easy-to-use, intuitive tool to quickly make building effect and plume rise calculations, as well as, inhalation and plume shine dose calculations. No experience with or working knowledge of Fortran, Excel, command line, or macros required. Self-contained package allows VENTSAR be deployed to any Windows computer without the need for specific programs or plug-ins.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

TEM study of the (SbS){sub 1+δ}(NbS{sub 2}){sub n}, (n=1, 2, 3; δ~1.14, 1.20) misfit layer phases

In the Sb–Nb–S system four new misfit layer phases have been found and carefully investigated via Transmission Electron Microscopy (TEM). Their structures are of composite modulated structure type with stoichiometries that can be formulated as (SbS){sub 1+δ}(NbS{sub 2}){sub n}; for n=1, δ~1.14 and 1.19; for n=2, δ~1.18 and for n=3, δ~1.19. Selected Area Electron Diffraction (SAED) patterns show an almost commensurate fit between the pseudo-tetragonal (SbS) and the pseudo-orthohexagonal (NbS{sub 2}){sub n} subcells along the misfit direction a, with 3(SbS)≈5(NbS{sub 2}), being b the same for both sub-lattices and c the stacking direction. For n=1, a commensurate phase with 4a{sub SbS}=7a{sub NbS2} has also been observed. In addition to the characteristic misfit and associated modulation of the two sub-structures, a second modulation is also present which appears to be primarily associated with the (SbS) sub-structure of both the n=1 and n=2 phases. High Resolution Transmission Electron Microscopy (HRTEM) images show ordered stacking sequences between the (SbS) and (NbS{sub 2}){sub n} lamellae for each of the four phases, however, disordered intergrowths were also occasionally found. Most of the crystals showed different kinds of twinning defects on quite a fine scale. Many crystals showed curled up edges. In some cases the lamellar crystals were entirely folded giving rise to similar diffraction patterns as found for cylindrical crystals. - Graphical abstract: Idealized structure models of the first three members of the homologous series (SbS){sub 1+δ}(NbS{sub 2}){sub n}. - Highlights: • Transmission Electron Microscopy study of misfit layer sulfides (SbS){sub 1+δ}(NbS{sub 2}){sub n}. • The structures consist of a (SbS) layer interleaved between n (NbS{sub 2}) layers. • Two different members n=1, one n=2 and one n=3 have been studied. • Twinning, intergrowths and different modulations in the (SbS) substructure.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Ground and excited state gradients with end-to-end differentiable semiempirical quantum chemistry

Accurate and efficient gradients of molecular energy with respect to nuclear degrees of freedom are essential for geometry optimization and molecular dynamics, including simulations that go beyond the Born–Oppenheimer regime. A common approach involves deriving analytical formulas for new electronic structure methods, which is often conceptually difficult and requires tedious coding. Here, we implement analytical, semi-numerical, and automatic differentiation (AD)-based gradient pathways for semiempirical Hamiltonian models in the PYSEQM software package, leveraging both graphics processing unit (GPU) and central processing unit (CPU) architectures. We further extend these capabilities to excited states calculated using the configuration interaction singles and time-dependent Hartree–Fock ansätze. We benchmark wall time, peak memory usage, and accuracy across three molecular families of varying chemical complexity, including systems of up to a thousand atoms. For ground-state simulations, analytical and AD gradients achieve near-identical GPU runtimes, while semi-numerical gradients are slower on GPU but remain competitive on CPU. For excited states, both analytical and custom AD approaches using implicit differentiation show similar performance and low memory requirements, whereas gradients with full AD are memory-limited. AD gradients match analytical ones in accuracy across all tested systems, aided by a quaternion-based diatomic frame rotation for two-center quantities that ensures smooth energy surfaces. Overall, automatic differentiation emerges as a practical alternative to analytical gradients in semiempirical quantum chemistry, offering high accuracy while allowing seamless integration in AI-driven workflows and popular packages, such as PyTorch and JAX. Our results provide actionable guidance for selecting optimal gradient strategies in large-scale ground- and excited-state molecular dynamics simulations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Acute phytotoxicity of seven metals alone and in mixture: Are Italian soil threshold concentrations suitable for plant protection?

Metals can pollute soils in both urban and rural areas with severe impacts on the health of humans, plants and animals living there. Information on metal toxicity is therefore important for ecotoxicology. This study investigated the phytotoxicity of different metals frequently found as pollutants in soils: arsenic, cadmium, chromium, lead, mercury, nickel and zinc. Cucumber (Cucumis sativus), sorghum (Sorghum saccharatum) and cress (Lepidium sativum) seeds were used as models for other plants used in human nutrition such as cereals, rice, fruits and vegetables. The 72-h germination rate and root elongations were selected as short-term ecotoxicological endpoints in seeds exposed to single metals and mixtures. Metals were spiked onto OECD standard soils in concentrations comparable to current Italian contamination threshold concentrations for residential and commercial soils. Arsenic, chromium, mercury and nickel were the most toxic metals in our experimental conditions, particularly to cress seeds (5.172, 152 and 255.4 mg/kg as 72 h IC50 for arsenic, mercury and nickel respectively). Italian limits were acceptable for plant protection only for exposure to each metal alone but not for the mixtures containing all the metals concentrations expected by their respective legislative threshold. The effects of the mixture were class-specific: trends were comparable in dicots but different in monocots. The response induced by the mixture at high concentrations differed from that theoretically obtainable by summing the effects of the individual metals. This might be due to partial antagonism of the metals in soil or to the formation of complexes between the metals, which reduce the bioavailability of the pollutants for plants. - Graphical abstract: Metals investigated: Arsenic, Cadmium, Chromium, Lead, Mercury, Nickel and Zinc. - Highlights: • The short-term phytotoxicity of seven metals was investigated with 3 higher plants. • Italian limits for arsenic and nickel in residential soils can damage plants. • Root elongation is heavily affected by metals. • Italian limits are not suitable for plant protection when metals are in mixture.

54 ENVIRONMENTAL SCIENCES↗

An Applied Strategy for Using Empirical and Hybrid Models in Online Monitoring

The monitoring of plant equipment for failure prediction is one of the key contributors to operation and maintenance (O&M) costs for a nuclear power plant (NPP) because O&M monitoring depends on labor-intensive activities that are required to meet high equipment reliability standards. These activities rely primarily on humans for information gathering, condition diagnosis, and predictive analysis. Online monitoring aims to automate these activities by relying on sensors to replace human information gathering and machine learning to replace human analysis and decision making. To facilitate automated monitoring, a systematic strategy for anomaly detection is needed to optimally use the available sensor data, empirical models, and physics-supported models. This strategy is essential to provide credible reasoning on why and when an empirical (i.e., purely data-driven) versus hybrid (i.e., physics-supported) approach should be used and to determine the ideal mix of these two approaches for a defined anomaly detection scope. The extant methods usually adopt an ad hoc trial-and-error approach that, in addition to being time-consuming and costly, is also highly subjective; it is impacted by the background and the skill set of the personnel making the decisions. Thus, such an approach cannot guarantee an optimum outcome. This represents the motivation of the current research effort, which is focused on devising a scientifically supported strategy for the optimum selection of anomaly detection methods. This report presents a detailed assessment of the main anomaly detection techniques within the empirical or hybrid method streams. Empirical methods include pattern, statistical, and causal inference. Hybrid methods include the use of physics models to train and test data methods, reduce data dimensionality, reduce data-model complexity, augment data, and reduce empirical uncertainty; hybrid methods also include the use of data to tune physics models. The listed techniques within these two streams represent the vast majority of techniques performed for anomaly detection. Using the techniques as outcomes, a strategy was developed to enable a systematic decision-making process to lead to one of these techniques. The strategy is driven by key decision points related to data relevance, simple modeling feasibility, data inference, physics-modeling value, data dimensionality, physics knowledge, method of validation, performance, data availability and suitability for training and testing, cause-effect, entropy inference, and model fitting. Each of these decision points in the strategy is explained in detail in this report with examples, along with the scientific basis behind the decisions and outcomes in common and simplified terminology. The strategy is developed for use by any NPP staff with basic engineering or science knowledge. A user-friendly graphical state flow diagram was also developed as a visual presentation of the strategy. The strategy was tested and demonstrated through two pilot projects for the application of anomaly detection at an NPP. Each pilot had two use cases: an initial case in which certain decisions were made that resulted in one or more empirical techniques and a revised use case where one or more key decisions were modified resulting in using a set of hybrid methods.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Alkali metal – yttrium borohydrides: The link between coordination of small and large rare-earth

The system Li–A–Y–BH{sub 4} (A=K, Rb, Cs) is found to contain five new compounds and four further ones known from previous work on the homoleptic borohydrides. Crystal structures have been solved and refined from synchrotron X-ray powder diffraction, thermal stability of new compounds have been investigated and ionic conductivity measured for selected samples. Significant coordination flexibility for Y{sup 3+} is revealed, which allows the formation of both octahedral frameworks and tetrahedral complex anions with the tetrahydroborate anion BH{sub 4} both as a linker and terminal ligand. Bi- and trimetallic cubic double-perovskites c-A{sub 3}Y(BH{sub 4}){sub 6} or c-A{sub 2}LiY(BH{sub 4}){sub 6} (A=Rb, Cs) form in all the investigated systems, with the exception of the Li–K–Y system. The compounds with the stoichiometry AY(BH{sub 4}){sub 4} crystallize in all investigated systems with a great variety of structure types which find their analog amongst metal oxides. In-situ formation of a new borohydride – closo-borane is observed during decomposition of all double perovskites. - Graphical abstract: The system Li–A–Y–BH{sub 4} (A=K, Rb, Cs) is found to contain five novel compounds and four further ones previously reported. Significant coordination flexibility of Y{sup 3+} is revealed, which can be employed to form both octahedral frameworks and tetrahedral complex anions, very different structural topologies. Versatility is also manifested in three different simultaneously occurring coordination modes of borohydrides for one metal cation, as proposed by DFT optimization of the monoclinic KY(BH{sub 4}){sub 4} structural model observed by powder diffraction. - Highlights: • The system Li-A-Y-BH{sub 4} (A=K, Rb, Cs) contains nine compounds in total. • Y{sup 3+} forms octahedral frameworks and tetrahedral complex anions. • Bi- and trimetallic double-perovskites crystallize in most systems. • Various AY(BH{sub 4}){sub 4} crystallize with structure types analogous to metal oxides. • Double-perovskites decompose and form a novel borohydride-closo-borane.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

TOPAS-nBio: An Extension to the TOPAS Simulation Toolkit for Cellular and Sub-cellular Radiobiology

The TOPAS Monte Carlo (MC) system is used in radiation therapy and medical imaging research, having played a significant role in making Monte Carlo simulations widely available for proton therapy related research. While TOPAS provides detailed simulations of patient scale properties, the fundamental unit of the biological response to radiation is a cell. Thus, our goal was to develop TOPAS-nBio, an extension of TOPAS dedicated to advance understanding of radiobiological effects at the (sub-)cellular, (i.e., the cellular and sub-cellular) scale. TOPAS-nBio was designed as a set of open source classes that extends TOPAS to model radiobiological experiments. TOPAS-nBio is based on and extends Geant4-DNA, which extends the Geant4 toolkit, the basis of TOPAS, to include very low-energy interactions of particles down to vibrational energies, explicitly simulates every particle interaction (i.e., without using condensed histories) and propagates radiolysis products. To further facilitate the use of TOPAS-nBio, a graphical user interface was developed. TOPAS-nBio offers full track-structure Monte Carlo simulations, integration of chemical reactions within the first millisecond, an extensive catalogue of specialized cell geometries as well as sub-cellular structures such as DNA and mitochondria, and interfaces to mechanistic models of DNA repair kinetics. We compared TOPAS-nBio simulations to measured and published data of energy deposition patterns and chemical reaction rates (G values). Our simulations agreed well within the experimental uncertainties. Additionally, we expanded the chemical reactions and species provided in Geant4-DNA and developed a new method based on independent reaction times (IRT), including a total of 72 reactions classified into 6 types between neutral and charged species. Chemical stage simulations using IRT were a factor of 145 faster than with step-by-step tracking. Finally, we applied the geometric/chemical modeling to obtain initial yields of double-strand breaks (DSBs) in DNA fibers for proton irradiations of 3 and 50 MeV and compared the effect of including chemical reactions on the number and complexity of DSB induction. Over half of the DSBs were found to include chemical reactions with approximately 5% of DSBs caused only by chemical reactions. In conclusion, the TOPAS-nBio extension to the TOPAS MC application offers access to accurate and detailed multiscale simulations, from a macroscopic description of the radiation field to microscopic description of biological outcome for selected cells. TOPAS-nBio offers detailed physics and chemistry simulations of radiobiological experiments on cells simulating the initially induced damage and links to models of DNA repair kinetics.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Thermal and mechanical stabilization process of the organic fraction of the municipal solid waste

Graphical abstract: Display Omitted - Highlights: • A domestic scale prototype for the pre-treatment of OFMSW has been tested. • Two grinding techniques are compared and thermopress is used for the drying stage. • Increasing temperature up to 170 °C reduces energy consumption of the drying stage. • In the range 5–10 bar a reduction of 97% of the initial volume is obtained. • In most cases energy recovery from the dried waste matches energy consumption. - Abstract: In the present study a thermo-mechanical treatment for the disposal of the Organic Fraction of Municipal Solid Waste (OFMSW) at apartment or condominium scale is proposed. The process presents several advantages allowing to perform a significant volume and moisture reduction of the produced waste at domestic scale thus producing a material with an increased storability and improved characteristics (e.g. calorific value) that make it available for further alternative uses. The assessment of the applicability of the proposed waste pretreatment in a new scheme of waste management system requires several research steps involving different competences and application scales. In this context, a preliminary study is needed targeting to the evaluation and minimization of the energy consumption associated to the process. To this aim, in the present paper, two configurations of a domestic appliance prototype have been presented and the effect of some operating variables has been investigated in order to select the proper configuration and the best set of operating conditions capable to minimize the duration and the energy consumption of the process. The performances of the prototype have been also tested on three model mixtures representing a possible daily domestic waste and compared with an existing commercially available appliance. The results obtained show that a daily application of the process is feasible given the short treatment time required and the energy consumption comparable to the one of the common domestic appliances. Finally, the evaluation of the energy recovered in the final product per unit weight of raw material shows that in most cases it is comparable to the energy required from the treatment.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Multi-node and Multi-core Performance Studies of a Monte Carlo code RMC

Over the past 15 years the architecture of high-performance computing systems has become increasingly heterogeneous with individual computing cores being subdivided into nodes and even auxiliary co-processing units such as graphics processing units or the Intel many integrated core co-processors. In fact, future computer platforms is moving toward a larger number of nodes and cores with lower memory available per node. Today, several production Monte Carlo codes use some types of hybrid parallelism strategy that makes use of shared-memory parallelism within a computing node and message-passing parallelism between nodes, which enables the codes to utilize available computing resources as much as possible, both among and within computing nodes on a cluster. The so-called hybrid parallelization is also possible with Reactor Monte Carlo Code RMC by activating both MPI and OpenMP parallel algorithms together. On HPC' which typically consists of several or several-ten thousand nodes equipping 16-32 computing intra-cores per node, simulation of RMC can be carried out in the hybrid parallelization with the MPI inter-node parallelization and the OpenMP intra-node parallelization. RMC has adopted a hybrid message-passing and shared-memory approach to parallelism, which enables the code to utilize available computing resources, both among and within computing nodes on a cluster. Shared-memory parallelism is especially valuable for large-model simulations because it enables multiple compute cores within a node to share common model information such as the geometry description, material, cross section and tally data, thus reducing the memory usage per processor. The calculation on the Tianhe-2 supercomputer took 6144 cores, resulting in a rate of 173,314,000 neutrons/minute for BEAVRS benchmark. While pin-wise full core H-M benchmark burnup calculations with millions of depletion regions was selected using MPI/OpenMP parallelism, the speed of 6,144 cores reach 1,843,200 burnable regions per minute. Keeping up with the continuing trend for increases in the number of computing resources within modern HPC systems is a significant challenge for teams that develop and maintain high-performance MC radiation transport solvers. order to address this challenge, RMC has adopted a hybrid message-passing and shared-memory approach to parallelism, which enables the code to utilize available computing resources, both among and within computing nodes on a cluster. Shared-memory parallelism is especially valuable for large-model simulations because it enables multiple compute cores within a node to share common model information such as the geometry description, material and composition definitions, and cross section data, thus reducing the memory usage per processor. The parallel performance has been studied by using two complex and representative benchmarks, H-M and BEAVRS benchmark. All simulations were run on Intel Xeon E5-2692 12 core, 2.2 GHz processors at MilkyWay-2 Super Computer with up to 12,288 available computing cores. While employing MPI parallelism, excellent parallel efficiency is achieved even above 12,000 processors and RMC reaches the parallel efficiency of 96% on 12,288 cores for the pin-by-pin computation of the BEAVRS benchmark. While BEAVRS criticality benchmark was selected using MPI/OpenMP parallelism, the speed of 6,144 cores reach 173,314,000 particles per minute. While pin-wise full core H-M benchmark burnup calculations with millions of depletion regions was selected using MPI/OpenMP parallelism, the speed of 6,144 cores reach 1,843,200 burnable regions per minute. Results from these studies demonstrate that RMC, using the hybrid parallel model described in this paper, is able to handle pin-wise full core burnup calculations with millions of depletion regions and scales well through thousands of processors. (authors)

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Real-Time Drilling Optimization System for Improved Overall Rate of Penetration and Reduced Cost Per Foot in Geothermal Drilling

The key to success in geothermal drilling is economic feasibility, and a major cost in the development of geothermal resources is the actual drilling of the wells. In this project, a real-time drilling optimization system for geothermal drilling was developed. The system couples three individual components while drilling. The first component is a drill stem vibration analysis model, the second is Mechanical Specific Energy (MSE) analyses, and the third is a detailed PDC Rate of Penetration (ROP) drill bit model for optimum RPM and WOB combinations. The benefit of the coupled system is that the range of WOB and RPM could be selected to avoid drill stem vibrations. Secondly, MSE is used as an efficiency measure and the detailed PDC drill bit model ensures the drill bit does not endure temperatures that exceed the temperature at which the PDC cutters experience accelerated wear. The new detailed PDC bit model is based on rock/bit interaction that physically tracks the PDC cutter wear flats as the bit drills ahead giving the capability to calculate the temperature being generated underneath the worn cutters to better advise on operational parameters to avoid accelerated cutter wear and failure and to ensure that operational parameters are applied so that overall ROP is maximized. By combining the drill stem vibrations and the detailed PDC bit cutter wear and “safe” non-accelerated cutter wear temperature and optimum ranges of operating parameters, it results in higher ROP and lower cost drilling. Single cutter PDC testing performed in different lithologies at Sandia was utilized to verify the PDC cutter forces and depth of cut for new and worn cutters. Based on single cutter PDC temperature modeling, verification using single cutter data from the testing done by National Oilwell Varco (NOV) was performed. Sandia’s Hard-Rock Drilling Facility (HRDF) was utilized to test different drill bit configurations with different cutter designs and wear status with different induced modes of vibration to obtain the critical bit RPM/WOB ranges resulting in ineffective drilling and low ROP. The collected test data were further used to verify and calibrate the full hole PDC ROP model that was developed based on single cutter interaction data. A full coupled drill stem vibration model was formulated and verified with geothermal field data from the Chocolate Mountain Aerial Gunnery Range (CMAGR). A graphical user interface (GUI) was developed using Tkinter library in the computer programming language Python, which integrates all the developed models in one system. The developed system consists mainly of the PDC ROP model, PDC bit wear model, PDC cutter temperature model, Mechanical Specific Energy (MSE) model, and drillstring vibration model integrated into one system. The developed system can be used for both, post well analysis and real-time optimization using different criteria such as ROP maximization or MSE minimization. The software uses Differential Evolution Algorithm (DEA) to find optimum values for operational parameters based on last foot drilled while avoiding the drillstring vibration and cutter temperature critical operating parameters.

15 GEOTHERMAL ENERGY↗

Radiological Performance Modeling: Platform Comparison

Zapadnoe Uranium Mill Tailings Site, Pridneprovsky Chemical Plant Ukraine: Two exposure scenarios were utilized to develop performance assessment models in GoldSim and NORMALYSA. The first scenario is an onsite worker exposed to radioactivity through external ground irradiation as well as radon and radioactive dust inhalation. The second being an offsite Resident 300 years after institutional control exposed through the ingestion of garden vegetables irrigated by contaminated ground water, ingestion of fish from Konoplyanka River, and agricultural activities in contaminated garden. NORMALYSA utilizes a modular interface to connect extensive pre-created modules to develop contaminant transport and risk assessment models for safety assessments and performance assessments at NORM and legacy sites. Developed by Facilia AB, the program was tested with the creation of the Zapadnoe model under the International Atomic Energy Agency MODARIA Il Working Group 1 (WG1). Default values are provided for input parameters where site specific data and distributions lack. Probabilistic simulations and inputs are limited to parameters with assigned Gaussian probability functions. Similar results for transient worker inhalation dose were found in comparison with the Zapadnoe GoldSim model. However, the NORMALYSA model has results larger by an order of magnitude for external radiation in the transient worker and for aquatic food ingestion, inhalation, and garden food ingestion to the resident. Performance Assessment Modeling: Radiological performance models promote good management and support in remediation and environmental decision making by optimizing radioactive waste disposal site performance. a performance assessment model of the Zapadnoe uranium mill tailings site and a model for calculating authorized limits at Los Alamos National Laboratory are created to demonstrate and highlight the respective capabilities of three modeling platforms: GoldSim, NORMALYSA, and RESRAD-Onsite. GoldSim is a general-purpose modeling platform with the ability to simulate site-specific FT radionuclide transport and risk assessment using a graphical interface to create and link elements that act as indivisible building blocks of the model. The Zapadnoe performance assessment model created in GoldSim is composed of input parameters collected from site specific parameters from WG1 as well as standard assumptions seen in the Generic PA Model from Neptune and Company. The multivariate Zapadnoe simulation contains both deterministic and probabilistic simulation capabilities with the ability to quantitatively represent the uncertainty and associated risks of the site via the Monte Carlo simulation method. Model parameters, documentation, and assumptions are built directly into the model creating a transparent, reproducible, and more purposeful site specific model. To demonstrate this, an ingestion of water pathway was included in addition to the exposure pathways included in the NORMALYSA model. GoldSim offers a wide variety of probability distributions to choose from as well as many tools to visualize the uncertainty of the models and explore which pathways and radionuclides dominate dose. In contrast to the NORMALYSA model, the GoldSim model uses dynamic simulation to calculate tritium concentrations. Authorized limits support the evaluation of radiological health effects from radionuclides present in bulk contaminated soil in order to assist in the process of clearing land for release under different uses. Calculations of authorized limits in RESRAD-Onsite and GoldSim models support the evaluation for 15 radionuclides using a residential exposure scenario, indoor and outdoor inhalation of dust and vapor phase radionuclides (tritium), incidental soil ingestion, and ingestion grown fruits and Vegetables. RESRAD-Onsite, developed by Argonne National Laboratory is a dedicated computer code for calculating radiation doses and cancer risks to an individual located on top of radioactively contaminated soils. It can also be used to provide guidelines for deriving radionuclide soil levels corresponding to a specific dose criterion. The menu driven user interface facilitates selection of radionuclides and soil concentrations, as well as defines parameter values related to materials, environmental factors, and exposure pathways. The program contains radionuclide-specific data such as dose conversion factors and risk coefficients, as well as soil-water partition coefficients and transfer factors for biotic media. The GoldSim model successfully emulates the RESRAD-Onsite code using identical radiological decay and ingrowth factors, radiation dose conversion factors, and other input parameters. DE EEDE ETS Each model supports dynamic modeling of radionuclide ingrowth and decay over time, and integrate this capability with the transport of radionuclides in the environment via processes such as advection and/or diffusion in fluids, or within porous solid media, The applicability of each program is heavily dependent on project context. In the context of developing Authorized Limits, RESRAD-Onsite may be a more desirable option as it was developed for this particular application. For simple disposal or legacy sites or sites that do not need site-specific model structure, NORMALYSA may be a desirable modeling platform, GoldSim can be used to model any problem for which RESRAD or NORMALYSA is applicable, but its strength lies in its flexibility, analytical, and probabilistic capabilities.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Linear and nonlinear instability in vertical counter-current laminar gas-liquid flows

We consider the genesis and dynamics of interfacial instability in vertical gas-liquid flows, using as a model the two-dimensional channel flow of a thin falling film sheared by counter-current gas. The methodology is linear stability theory (Orr-Sommerfeld analysis) together with direct numerical simulation of the two-phase flow in the case of nonlinear disturbances. We investigate the influence of two main flow parameters on the interfacial dynamics, namely the film thickness and pressure drop applied to drive the gas stream. To make contact with existing studies in the literature, the effect of various density contrasts is also examined. Energy budget analyses based on the Orr-Sommerfeld theory reveal various coexisting unstable modes (interfacial, shear, internal) in the case of high density contrasts, which results in mode coalescence and mode competition, but only one dynamically relevant unstable interfacial mode for low density contrast. A study of absolute and convective instability for low density contrast shows that the system is absolutely unstable for all but two narrow regions of the investigated parameter space. Direct numerical simulations of the same system (low density contrast) show that linear theory holds up remarkably well upon the onset of large-amplitude waves as well as the existence of weakly nonlinear waves. For high density contrasts, corresponding more closely to an air-water-type system, linear stability theory is also successful at determining the most-dominant features in the interfacial wave dynamics at early-to-intermediate times. Nevertheless, the short waves selected by the linear theory undergo secondary instability and the wave train is no longer regular but rather exhibits chaotic motion. The same linear stability theory predicts when the direction of travel of the waves changes — from downwards to upwards. We outline the practical implications of this change in terms of loading and flooding. The change in direction of the wave propagation is represented graphically in terms of a flow map based on the liquid and gas flow rates and the prediction carries over to the nonlinear regime with only a small deviation.

42 ENGINEERING↗

An Overview of the Molten Salt Thermal Properties Database–Thermophysical, Version 4.0 (MSTDB-TP V.4.0)

A central repository of thermophysical and thermochemical properties of molten salt compositions of relevance to molten salt reactors (MSRs) is vital in supporting the broad community of MSR developers, who are at various stages of developing and deploying their reactor designs. In general, these MSR designs differ significantly from developer to developer (e.g., with respect to the hardness of the neutron spectra, level of fissile loading, target multicomponent temperatures and power levels, and moderating capabilities). Therefore, the fuel and coolant salts being considered vary greatly: they may be chlorides or fluorides, they utilize different actinides at different ratios, and the cations in the melt are selected based on perceived advantages and disadvantages. Considering the general need for thermal properties, and the vastness of the array of potential candidate salt mixtures, the Molten Salt Thermal Properties Database (MSTDB) was initiated in 2018 with the goal of providing thermophysical and thermochemical characterization of key molten salt compounds and mixtures across their temperature and compositional domains. The MSTDB is thus divided into the thermophysical arm (MSTDB-TP) and the thermochemical arm (MSTDB-TC). The MSTDB is an effort funded by the Department of Energy, Office of Nuclear Energy (DOE-NE) Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, and the MSR Campaign. This report provides an overview of the MSTDB-TP v4.0 in terms of the data contained within, the state of the tools used to access the data, the availability of predictive models that leverage the raw data in the database, the preliminary status of developmental efforts that are currently underway, and an account of future goals for MSTDB-TP. The primary goal for the update from MSTDB-TP v.3.1 to v4.0 was the incorporation of surface tension data into the database; this property is important for thermal hydraulics modeling and species transport in other tools that have been developed under the NEAMS program. A breakdown of the surface tension data that have been added into MSTDB-TP v4.0 is provided herein, and the manner in which the quality of the data has been assessed is also documented. For MSTDB-TP v4.0, newly published thermophysical property data—primarily from collaborative experimental efforts under the MSR Campaign—have been incorporated into the database, and the resulting expansion is documented here. Because of the size to which MSTDB-TP has grown, the raw data format has now been recast into JavaScript Object Notation (JSON) format for easier connection with the MSTDB-TP application programming interface (API). Saline; the pre-existing comma-separated value (CSV) format has been deprecated but is still maintained, accessible, and up to date. As a final effort in packaging the MSTDB-TP v4.0 update, the graphical user interface (GUI) for MSTDB has been updated to allow full accessibility to the density and viscosity predictive models, which are based on Redlich-Kister expansions of MSTDB-TP raw data. Some other major aspects of this report, in terms of preliminary and future work, include: (1) documentation of the formalism and preliminary testing of a kinetic theory model that may act as a predictive model for thermal conductivity; (2) documentation of the candidate predictive models that may be considered in the future for surface tension, making use of the surface tension data now in MSTDB-TP v4.0; (3) a preliminary account of a data collection process that will enable the filling of additional gaps within MSTDB-TP, namely with data which have been collected computationally (e.g., through ab initio molecular dynamics).

22 GENERAL STUDIES OF NUCLEAR REACTORS↗