Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “performance gain”

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,171 records · Page 65

Neuroevolution Application to Collaborative and Heuristics-Based Connected and Autonomous Vehicle Cohort Simulation at Uncontrolled Intersection

Artificial intelligence is gaining tremendous attractiveness and showing great success in solving various problems, such as simplifying optimal control derivation. This work focuses on the application of Neuroevolution to the control of Connected and Autonomous Vehicle (CAV) cohorts operating at uncontrolled intersections. The proposed method implementation’s simplicity, thanks to the inclusion of heuristics and effective real-time performance are demonstrated. The resulting architecture achieves nearly ideal operating conditions in keeping the average speeds close to the speed limit. It achieves twice as high mean speed throughput as a controlled intersection, hence enabling lower travel time and mitigating energy inefficiencies from stop-and-go vehicle dynamics. Low deviation from the road speed limit is hence continuously sustained for cohorts of at most 50 m long. This limitation can be mitigated with additional lanes that the cohorts can split into. The concept also allows the testing and implementation of fast-turning lanes by simply replicating and reconnecting the control architecture at each new road crossing, enabling high scalability for complex road network analysis. The controller is also successfully validated within a high-fidelity vehicle dynamic environment, showing its potential for driverless vehicle control in addition to offering a new traffic control simulation model for future autonomous operation studies.

Jacquelin, Frederic (ORCID:0000000183154344)↗

An Assessment of the Influence of Uncertainty in Temporally Evolving Streamflow Forecasts on Riverine Inundation Modeling

Continental-scale river forecasting platforms forecast streamflow at reaches that can be used as boundary conditions to drive a local-scale flood inundation model. Uncertainty accumulated during this process stems not only from any part of the forecasting chain but can also be caused by the daily variations in weather forcing that keeps evolving as the event advances. This work aims to examine the influence of the evolving forecast streamflow on predicting the maximum inundation for extreme floods. A diagnostic case study was made on the basis of a hindcast of Hurricane Matthew striking the eastern U.S. in 2016. The U.S. National Water Model was one-way coupled to a hydrodynamic inundation model through a developed automated workflow. Although the river forcing has significantly mismatched hydrographs versus observations, the simulated peak water surface elevations and maximum extents were validated to be comparable with the observations, which indicates that the inundation model may not be sensitive to the inherited uncertainty from the weather forcing. Moreover, the uncertainty of the forecast streamflow time series caused only one order of magnitude fewer variations in inundation prediction; this dampening effect may become clearer for extreme events with large areas inundated. In addition, the forecast total volume of stream discharge appears to be an important metric for assessing the performance of river forcing for inundation mapping, as a linear correlation between the total volume and the accuracy of the predicted peak water surface elevation and maximum extent was found, with the coefficients of determination all above 0.8. Extra best-practice experience of running similar operational tasks demonstrated the tradeoff between the cost and accuracy gain.

flood, national water model, RIFT↗

Mechanistic Insights and Rational Design of Ca-Doped CeO 2 Catalyst for Acetic Acid Ketonization

Carboxylic acid ketonization has recently gained significant attention to produce biomass-derived hydrocarbon fuels as it not only removes the highly reactive carboxylic functional group but also increases the size of the carbon chain. In this work, Ca-doped CeO 2 -based catalysts were investigated for acetic acid ketonization using a combined experimental and computational approach. Acetic acid conversion was performed across a range of temperatures including higher temperatures relevant to catalytic hot gas filtration (450 °C). Ca addition slightly decreases overall acetic acid ketonization reactivity yet stabilizes the catalyst at the higher temperatures necessary for catalytic hot gas filtration. From density functional theory calculations of the ketonization reaction mechanism, the C–C coupling and water formation steps are identified as two of the most energy-consuming steps on a CeO 2 surface with a proximal oxygen vacancy and the presence of a Ca dopant stabilizes the key intermediates. Calculations predict an optimal structure comprising three Ca ensembles to minimize the reaction free energies for C–C coupling and water formation steps. These findings provide a priori information to guide future experiments for ketonization catalyst design and development.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Design and structural basis of selective 1,4-dihydropyridine inhibitors of the calcium-activated potassium channel K Ca 3.1

The 1,4-dihydropyridines, drugs with well-established bioavailability and toxicity profiles, have proven efficacy in treating human hypertension, peripheral vascular disorders, and coronary artery disease. Every 1,4-dihydropyridine in clinical use blocks L-type voltage-gated calcium channels. We now report our development, using selective optimization of a side activity (SOSA), of a class of 1,4-dihydropyridines that selectively and potently inhibit the intermediate-conductance calcium-activated K + channel K Ca 3.1, a validated therapeutic target for diseases affecting many organ systems. One of these 1,4-dihydropyridines, DHP-103, blocked K Ca 3.1 with an IC 50 of 6 nM and exhibited exquisite selectivity over calcium channels and a panel of >100 additional molecular targets. Using high-resolution structure determination by cryogenic electron microscopy together with mutagenesis and electrophysiology, we delineated the drug binding pocket for DHP-103 within the water-filled central cavity of the K Ca 3.1 channel pore, where bound drug directly impedes ion permeation. DHP-103 inhibited gain-of-function mutant K Ca 3.1 channels that cause hereditary xerocytosis, suggesting its potential use as a therapeutic for this hemolytic anemia. In a rat model of acute ischemic stroke, the second leading cause of death worldwide, DHP-103 administered 12 h postischemic insult in proof-of-concept studies reduced infarct volume, improved balance beam performance (measure of proprioception) and decreased numbers of activated microglia in infarcted areas. K Ca 3.1-selective 1,4-dihydropyridines hold promise for the many diseases for which K Ca 3.1 has been experimentally confirmed as a therapeutic target.

Ong, Seow Theng [Lee Kong Chian School of Medicine↗

Deliquescence of Eutectic LiCl-KCl Diluted with NaCl for Interim Waste Salt Storage

Molten eutectic LiCl-KCl salt is a widely used electrolyte for electrorefining uranium from spent nuclear fuel. Due to the hygroscopic nature of this salt, such operations must be performed under controlled atmospheric conditions, and waste salts require careful storage to avoid deliquescence and corrosion of container materials. This study investigated a potential processing path for reducing the degree of deliquescence through dilution to varying extents with NaCl. Here, the hydration behavior of LiCl-KCl salts diluted with NaCl was evaluated in terms of mass gain due to water absorption, degree of deliquescence (including first appearances of standing water), and evidence of corrosion to stainless steel containers in a humid air environment (40°C, 20% relative humidity). In this humid air environment, pure eutectic LiCl-KCl exhibited a 50 mass % increase due to water absorption and showed evidence of standing water after 24 h. Waste salt diluted with NaCl required loadings of 89 mass % NaCl in order to prevent deliquescence and exhibited a 3 mass % increase due to water absorption. After periodic observation for 48 h, standing water was observed near all ingots with the exception of the 89 mass % NaCl samples. Dilution with 89% NaCl was also found to reduce evidence of corrosion when stored in stainless steel crucibles. While dilution with NaCl greatly decreases steady-state hydration, the storage volume is increased ~10× through this procedure.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Modeling Molten Salt Spreading and Heat Transfer using MELTSPREAD – An Uncertainty Analysis

This report summarizes spreading and heat transfer results for molten salt pouring onto a flat stainless steel substrate that were calculated using MELTSPREAD to gain insight into bulk salt behavior during a molten salt spill accident. MELTSPREAD was developed at Argonne to model the one-dimensional flow and freezing behavior of molten corium and is being applied to model molten salts for the first time as part of this work. An uncertainty analysis of thermophysical properties and spill conditions was performed using MELTSPREAD to determine which properties and conditions have the greatest impact on model outcome. The model was run with and without the inclusion of decay heat. Eutectic FLiNaK was used as the salt composition because it is well-characterized and appropriate for model development. Small salt volumes were used for the initial model runs described in this report so that the model results can be compared to the results from ongoing experiments on molten salt spreading and heat transfer currently being conducted at Argonne at a benchtop scale. The spreading of a small volume of FLiNaK (50 mL) was found to be limited by the balance between the molten salt surface tension and gravity and not by freezing for the pour conditions and assumed mechanisms of heat transfer in the model. Heat transfer from the salt was highly inefficient to due to its low thermal conductivity, high heat capacity, and high heat of fusion and varying these properties had no effect on the spreading behavior. Changing the viscosity and density of the salt slightly affected the spreading behavior. Setting the salt surface tension to approximately zero to simulate the salt wetting the substrate had the greatest effect on the spreading behavior and this increased both the spreading velocity and overall spreading area. Including decay heat in the model had a significant effect on the long-term cooling behavior of the salt but had no effect on the overall spreading behavior. Future effort should involve modeling a larger volume of molten salt to represent a reactor case spill scenario because the spreading behavior of large volumes of salt may differ from the small volumes that were the focus of this document.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Capsule defects measurement [Slides]

Recent high yield shots a NIF showed that capsule defects are limiting shot performance. Voids in the capsule are not well understood in how they degrade yield performance. Recent simulations show that voids can reduce compressibility without any diagnostic signatures but reduced yield. In extreme cases loss of compression is expected. Up to 10,000 voids might be necessary to measure the different effects. Experiments using capsules with engineered defects can help guide simulations and improve our understanding of the physics involved. The knowledge gained can be applied to other problems.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Modeling and Simulation of Fuel Dispersal During the Loss-of-Coolant Accident

This document is the compilation of the milestone portion to a larger end of project NEUP report. The executive summary of the modeling portion is provided below: In the event of cladding rupture during a postulated LOCA in a pressurized water reactor, fuel particles, along with fission gases, can be expelled into the reactor core from the fractured fuel rod, a phenomenon referred to as fuel dispersal. The initial stage of fuel dispersal is strongly influenced by the high-pressure ejection of fuel fragments, the size and geometry of the ruptured cladding, and the depressurization history of the fuel rod during the postulated LOCA transient. Depending on the location of the burst orifice relative to the quench front, the dispersal event represents an intricate three-phase flow and heat transfer phenomenon, where high-temperature fuel particles carried by the fission gases interact with the coolant within the narrow subchannels of the fuel assemblies, inducing localized phase change. Given the unique multiphysics nature of this phenomena, the current study develops a dedicated computational framework to predict the mass distribution and cooling of dispersing fuel particles, facilitating post-accident assessment and management of the fuel assemblies. Considering the scale of nuclear reactor applications, a continuum three-fluid model is proposed for simulating the transport of solids within the reactor core. With high-temperature fuel fragments within the liquid media, nucleation sites inducing phase changes are dispersed within the flow domain. Coupled with the fact that the transient dispersal event occurs on different time scales than other three-phase flow applications, this study derives a time-averaged three-fluid flow model without losing generality. The assumptions regarding the continuum treatment of the solid phase and the modeling of fuel dispersal behavior are incorporated to simplify the governing equations and derive applicable closure relations. The computational validation of the model was conducted using adiabatic experimental results obtained from ongoing research at Oregon State University, focusing on characterizing fuel dispersal behavior during simulated LOCA conditions. Settlement characteristics of the solids, quantified by the probability distribution of equivalent particles, closely matched the probability density functions reported in experimental studies. The transport of fuel particles within a scaled 5 × 5 lattice of a pressurized-water reactor rod bundle geometry was modeled through a two-fluid Eulerian framework. The required boundary conditions were evaluated from the fuel performance code BISON in a postulated large-break LOCA scenario. The modeling framework considered solid fuel particles as granular matter, interacting with the gaseous dry steam phase and fission gases through the governing interfacial momentum exchange between the participating fluids. The simulation results provided the volume fraction of the solids obtained at the bottom surface of the enclosing tank geometry. Postulated LOCA leading to fuel dispersal phenomena involves the strong coupling between fuel thermomechanics, cladding deformation, thermal-hydraulics, and fuel particle transport. Incorporation of such a strong coupling in numerical simulation is performed by coupling the multiphysics solvers. In the case of fuel dispersal, a strong coupled simulation can be performed by coupling the BISON code for fuel performance, the TRACE code for system-level thermal hydraulics, and fuel particle transport in Multiphysics Object-Oriented Simulation Environment (MOOSE). For such intricate infrastructure, the MOOSE Framework eases the data transfer between codes. The recent version of MOOSE has incorporated the Navier-Stokes module for the fluid flow. An exploratory exercise was done to gain familiarity with finite volume capabilities in the MOOSE framework to incorporate the Spalart-Allmaras (SA) turbulence model. New finite-volume and auxiliary kernels were introduced to assemble the SA transport equation, compute turbulent viscosity, and evaluate wall distance and diagnostic turbulence terms, fully integrated with existing Navier-Stokes modules. A turbulent lid-driven cavity at a Reynolds number of approximately 10,000 is used for verification. MOOSE shows the robust solver convergence and produces the turbulent features. But it underpredicts the velocity profile and turbulent quantities, emphasizing the need to develop improved SA near-wall treatments (e.g., low-Re corrections or wall functions) as a key direction for future work.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Numerical Analysis of Regular Material Point Method and its Application to Multiphase Flows

The material point method (MPM) is gaining wide popularity in engineering research to model and simulate complex multiphase flow dynamics. The method relies on solving the governing equations of motion and transport in a Lagrangian framework using particles also known as material points. The fluid and kinematic properties are stored on the material points while the spatial gradient calculation and temporal integration are performed on a background grid. This Lagrangian framework allows for large deformations, easy integration of constitutive models, and direct import of complex geometries as particles. However, despite their increasing popularity, very few studies have addressed the issues of numerical resolution and stability of MPM techniques. The presence of additional factors such as the number of material points-per-cell, the location of the material points, the CFL-like condition used in time update, and the grid shape functions also increase the complexity of the error analysis when compared to other finite element methods. In this presentation, we analyze the various forms of error incurred in the application of MPM to continuum mechanics and multiphase flows. The effect of the previously mentioned factors on the error dynamics is studied. The application of these principles to canonical and industrial problems is also presented.

high pressure reverse osmosis↗

Resilient cooling through geothermal district energy system

Decarbonization and resilience to heat waves have recently become high priorities for building and district energy systems. Geothermal coupled district heating and cooling systems that operate a water loop near ground temperature gain increasing adoption to support decarbonization. In these systems, vapor-compression machines, distributed in the energy transfer stations, lift the temperature up or down to the needs of the particular building. In principle, these systems can provide low-power, free cooling from the geothermal bore field during heat waves when electricity is often scarce. However, the performance of such a resilience operation mode and its implication on the energy system configuration and the sizing of the bore field and HVAC equipment is not yet understood. Consequently, we are assessing their resilience, power use and design implications under a scenario of a heat wave on five working days during which chillers are switched off to reduce electrical consumption. Our analysis is based on high-fidelity, coupled dynamic models of district energy, building-side HVAC and actual control logic, with whole building energy simulation used to assess thermal conditions in a 2004 vintage multi-zone office building in Chicago, IL. The results show that relying only on waterside economizer cooling, the indoor thermal conditions can be maintained in a tolerable range for the majority of the building zones with half the electrical energy compared to standard chiller operation. Thermal comfort in the hottest zones can be further improved by oversizing the cooling coil. However, the waterside economizer has significant implications on the system configuration and sizing: The geothermal bore field needs to be sized about 30% larger than the upper limit of the range observed for conventional geothermal systems. Nevertheless, if a central chiller plant is added, the bore field can be downsized to the typical design range. The latter configuration still allows compressor-less cooling during the heat wave with peak power reduced by 60% compared to the standard design and chiller operation.

15 GEOTHERMAL ENERGY↗

High-density stable glasses formed on soft substrates

Enabled by surface-mediated equilibration, physical vapour deposition can create high-density stable glasses comparable with liquid-quenched glasses aged for millions of years. Deposition is often performed at various rates and temperatures on rigid substrates to control the glass properties. Here, in this work, we demonstrate that on soft, rubbery substrates, surface-mediated equilibration is enhanced up to 170 nm away from the interface, forming stable glasses with densities up to 2.5% higher than liquid-quenched glasses within 2.5 h of deposition. Gaining similar properties on rigid substrates would require 10 million times slower deposition, taking ~3,000 years. Controlling the modulus of the rubbery substrate provides control over the glass structure and density at constant deposition conditions. These results underscore the significance of substrate elasticity in manipulating the properties of the mobile surface layer and thus the glass structure and properties, allowing access to deeper states of the energy landscape without prohibitively slow deposition rates.

36 MATERIALS SCIENCE↗

Performance of microchannel plate based detectors for <25 keV x rays: Monte Carlo simulations and comparisons with experimental results

We present the results of Monte Carlo simulations of the microchannel plate (MCP) response to x rays in the 250 eV to 25 keV energy range as a function of both x-ray energy and impact angle and their comparisons with the experimental results from the X8A beamline at the National Synchrotron Light Source at Brookhaven National Laboratory. Incoming x rays interact with the lead glass of the microchannel plate, producing photoelectrons. Transport of the photoelectrons is neglected in this model, and it is assumed that photoelectrons deposit all their energy at the point they are created. This deposition leads to the generation of many secondary electrons, some fraction of which diffuse to the MCP pore surface where they can initiate secondary electron cascades in the pore under an external voltage bias. X-ray penetration through multiple MCP pore walls is increasingly important above 5 keV, and the effect of this penetration on MCP performance is studied. In agreement with past measurements, we find that the dependence of MCP sensitivity with angle relative to the pore bias changes from a cotangent dependence to angular independence and then proceeds to a secant dependence as the x-ray energy increases. We also find that with the increasing x-ray energy, the MCP gain sensitivity as a function of bias voltage decreases. The simulations also demonstrate that for x rays incident normal to the MCP surface, spatial resolution shows little dependence on the x-ray energy but degrades with the increasing x-ray energy as the angle of incidence relative to the surface normal increases. This agrees with experimental measurements. Simulation studies have also been completed for MCPs gated with a subnanosecond voltage pulse. In this work, we find that the optical gate profile width increases as the x-ray energy is increased above 5 keV, a consequence of increased x-ray penetration at energies >5 keV. Simulations of the pulsed dynamic range show that the dynamic range varies between ~100 and 1000 depending on x-ray energy and peak voltage.

47 OTHER INSTRUMENTATION↗

Synchrotron light source X-ray detection with Low-Gain Avalanche Diodes

Abstract Low Gain Avalanche Diodes (LGADs) represent the state-of-the-art in timing measurements and will instrument the future Timing Detectors of ATLAS and CMS for the High-Luminosity LHC. While initially conceived as a sensor for charged particles, the intrinsic gain of LGADs makes it possible to detect low-energy X-rays with good energy resolution and excellent time resolution (tens of picoseconds). Using the Stanford Synchrotron Radiation Lightsource (SSRL) at SLAC, several LGADs designs were characterized with energies from 5 to 70 keV. The SSRL provides 10 ps pulsed X-ray bunches separated by 2 ns intervals with an energy dispersion (ΔE/E) of 10 -4 . LGADs from Hamamatsu Photonics (HPK) and Brookhaven National Laboratory (BNL) with different thicknesses ranging from 20 μm to 50 μm and different gain layer designs were read out using fast amplification boards and digitized with a high bandwidth and high sampling rate oscilloscope. PIN devices from HPK and AC-LGADs from BNL were characterized as well. A systematic and detailed characterization of the devices' energy linearity, resolution, and time resolution as a function of X-ray energy was performed for different biasing voltages at room temperature and are reported in this work. The charge collection and multiplication mechanism were simulated using Geant4 and TCAD Sentaurus, providing an important handle for interpreting the data.

Instruments & Instrumentation↗

Direct-drive double-shell implosion: A platform for burning-plasma physics studies

Double-shell ignition designs have been investigated with the indirect-drive inertial confinement fusion (ICF) scheme in both simulations and experiments in which the inner-shell kinetic energy was limited to ~10 to 15 kJ, even driven by megajoule-class lasers such as the National Ignition Facility. Since direct-drive ICF can couple more energy to the imploding shells, we have performed a detailed study on direct-drive double-shell (D 3 S) implosions with state-of-the-art physics models implemented in radiation-hydrodynamic codes (LILAC and DRACO), including nonlocal thermal transport, cross-beam energy transfer (CBET), and first-principles-based material properties. To mitigate classical unstable interfaces, we have proposed the use of a tungsten/beryllium–mixed inner shell with gradient-density layers that can be made by magnetron sputtering. In our D 3 S designs, a 70-μm-thick beryllium outer shell is driven symmetrically by a high-adiabat (α ≥ 10), 1.9-MJ laser pulse to a peak velocity of ~240 km/s. Upon spherical impact, the outer shell transfers ~30 to 40 kJ of kinetic energy to the inner shell filled with deuterium–tritium gas or liquid, giving neutron-yield energies of ~6 MJ in 1-D simulations. Two-dimensional high-mode DRACO simulations indicated that such high-adiabat D 3 S implosions are not susceptible to laser imprint, but the long-wavelength perturbations from the laser port configuration along with CBET can be detrimental to the target performance. Yet, neutron yields of ~0.3- to 1.0-MJ energies can still be obtained from our high-mode DRACO simulations. The robust α-particle bootstrap is readily reached, which could provide a viable platform for burning-plasma physics studies. Once CBET mitigation and/or more laser energy becomes available, we believe that breakeven or moderate energy gain might be feasible with the proposed D 3 S scheme.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Building and Sustaining a Community Resource for Best Practices in Scientific Software: The Story of BSSw.io

The development of scientific software—a cornerstone of long-term collaboration and scientific progress—parallels the development of other types of software but still poses distinct challenges, especially in high-performance computing. Although web searches yield numerous resources on software engineering, there is still a scarcity specifically for scientific software development. Here, this article introduces the Better Scientific Software site (https://bssw.io), a platform that hosts a community of researchers, developers, and practitioners who share their experiences and insights on scientific software development. Since 2017, this collaborative hub has gained traction within the scientific computing community, attracting a growing number of readers and contributors eager to share ideas and elevate their software development practices. In sharing the BSSw.io site’s story, we hope to encourage further growth of the BSSw.io community through both readership and contributors, with a long-term goal of fostering culture change by increasing emphasis on best practices in scientific software.

97 MATHEMATICS AND COMPUTING↗

Efficient Training of Deep Neural Operator Networks via Randomized Sampling

Neural operators (NOs) employ deep neural networks to learn the mappings between infinitedimensional function spaces. Deep operator network (DeepONet), a popular NO architecture, has demonstrated success in the real-time prediction of complex dynamics across various scientific and engineering applications. In this work, we introduce a random sampling technique to be adopted during the training of DeepONet, aimed at improving the generalization ability of the model, while significantly reducing the computational time. The proposed approach targets the trunk network of the DeepONet model that outputs the basis functions corresponding to the spatiotemporal locations of the bounded domain on which the physical system is defined. While constructing the loss function, DeepONet training traditionally considers a uniform grid of spatiotemporal points at which all the output functions are evaluated for each iteration. This approach leads to a larger batch size, resulting in poor generalization and increased memory demands, due to the limitations of the stochastic gradient descent (SGD) optimizer. The proposed random sampling over the inputs of the trunk net mitigates these challenges, improving generalization and reducing the memory requirements during training, resulting in significant computational gains. We validate our hypothesis through three benchmark examples, demonstrating substantial reductions in training time while achieving comparable or lower overall test errors relative to the traditional training approach. Our results indicate that incorporating randomization in the trunk network inputs during training enhances the efficiency and robustness of DeepONet, offering a promising avenue for improving the framework’s performance in modeling complex physical systems.

Karumuri, Sharmila [Department of Civil & Systems ↗

MOSCATO Development and Integration in Fiscal Year 2024

MOSCATO (Molten Salt Chemistry and Transport) is a multiphysics code that provides high-fidelity, coupled simulations of fluid flow, heat transfer, mass transfer, chemistry, electrochemical phenomena, and alloy evolution for molten salt equipment. In FY24, significant developments were made to the code package, enhancing its capabilities in many aspects. The improvements and advancements can be summarized as follows: 1. Implementation of tritium transport capabilities and validation with experimental data: To enable modeling of tritium and other fission gases within MSRs, we implemented gas transport within MOSCATO via inclusion of couple mass transport equations within the salt and structural alloys. Comparisons to experimental data from literature showed good agreement with respect to tritium release rates. 2. Preliminary implementation of two-phase flow models in MOSCATO: To model tritium and other gases above their solubility limits, we implemented preliminary two-phase flow models within MOSCATO to account for bubble transport. The first model adopted was the Level-Set approach, which can handle the high void fraction regime, but with a requirement for high mesh resolution thus high computational expense. In this report, we present a verification of the Level-Set method using a simple benchmark case. We also performed a demonstration of the code as applied to an experimental case involving cover gas flow through salt in an experimental vessel. The second model adopted was the Eulerian-Eulerian dispersed flow model, which is computationally cheaper but limited to low void fraction regimes, such as bubbly flow. Validation and verification have not yet been performed for the Eulerian-Eulerian approach, but a preliminary implementation was completed. 3. Validation with static corrosion experiments: Static corrosion experimental data for stainless steel coupons within molten salts was used to further validate the corrosion model in MOSCATO. To do so, we leveraged the existing models in MOSCATO and simulated the sample mass loss and mass gain phenomena. Several ion species, including Cr 2+ , Fe 2+ and H + , were simulated in salt using the PNP solver, while Cr 0 and Fe 0 were simulated with a diffusion solver in stainless steel. The mass loss of the samples was compared with experimental data, and good agreement was achieved. These combined activities served to further expand the capabilities of MOSCATO and make it more generally applicable to the full range of phenomena that can control chemistry and corrosion in molten salt reactors.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

A Spectroscopic Classification Survey to Search for New ρ Puppis Stars

The ρ Puppis stars are mid-F-type stars that show peculiar chemical-abundance patterns similar to those of the Am stars. Typically they exhibit overabundances of iron-peak elements such as Fe and Ni and s- and r-process elements such as Sr and Eu, and underabundances of certain other elements including He, Ca, and Sc. It has been proposed that these stars are evolving Am stars passing through the short-lived phase that occurs between the re-establishment of convection and the consequent erasure of their chemical peculiarities. In this paper we suggest a second hypothesis: these stars may have acquired their peculiar abundance patterns in a fashion similar to the Barium stars, i.e., they may have gained the s-process element enhancements via mass transfer from a once asymptotic giant-branch companion star, now turned white dwarf. This study will detail our efforts to investigate the ρ Puppis stars with regard to these two hypotheses, with a view to understanding these stars and their significance in terms of stellar evolution. We have performed a spectral classification survey, and a detailed chemical-abundance analysis of selected ρ Puppis stars. This paper gives an overview of the ρ Puppis stars and describes the spectral classification survey which has increased the number of ρ Puppis stars currently known from 6 to 49. A future paper will describe the chemical-abundance analysis which should provide insight into the nature of these stars.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗