Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Dynamic equivalent”

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.

501 records · Page 28
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Impedance Scan of Inverter-Based Resources and Diesel Generator for Stability Analysis

Impedance-based methods are widely used for power system stability analysis with inverter-based resources (IBRs), e.g., assessing dynamic interactions between the power grid and an IBR, control interactions between multiple IBRs, and the sub-synchronous oscillation and damping phenomenon. Since it is difficult to get a numerical model 100% matching with the hardware IBR, using the hardware inverter directly to obtain its output impedance has become a prominent approach nowadays. Therefore, this article presents the impedance scan using hardware IBRs, and also a hardware diesel generator as it still stays with the grid before the grid completely goes to renewable. The devices under test (DuTs) for the impedance scan includes two 3-..phi.., 480 V, 60 Hz commercial grid-forming IBRs (one of 250 kVA and another of 125 kVA rating) in series with ..delta..-Y transformers, one 3-..phi.., 480 V, 60 Hz commercial grid-following IBR (of 125 kVA rating), and a 3-..phi.., 480 V, 60 Hz commercial diesel generator (of 187.5 kVA rating). Using voltage signals perturbed with sub-, inter-, and higher harmonic components, and measuring the current response, the positive-sequence impedances are computed via an offline-based post-analysis. Moreover, best-fit transfer functions are estimated that closely resemble the measured data points of the positive-sequence impedances. Based on the observations from various outcomes of the hardware experiments, this article also provides some fundamental insights on the equivalent positive-sequence impedance of a combination of multiple hardware components by comparing the estimated and the empirically computed impedances. A comparative insight on the damping capability of the DuTs using the positive-sequence impedances of the hardware is also discussed.

current measurement↗
DOE OSTI2023

IDAES-PSE 2.0 Release

The Institute for the Design of Advanced Energy Systems (IDAES) Integrated Platform is a versatile computational environment offering extensive process systems engineering (PSE) capabilities for optimizing the design and operation of complex, interacting technologies and systems. IDAES enables users to efficiently search vast, complex design spaces to discover the lowest cost, most environmentally sustainable solutions while supporting the full process modeling lifecycle, from conceptual design to dynamic optimization and control. The extensible, open platform empowers users to create models of novel processes and rapidly develop custom analyses, workflows, and end-user applications. IDAES-PSE 2.0.0 Release Highlights Removal of deprecated features from IDAES v1 Update to Pyomo v6.5 – this required a number of updates to support the new NL solver writer and to address some changes in Pyomo Creation of new testing suite for backward compatibility, model robustness and verification More general implementation of the Helmholtz EoS. This brings some new features like standard property diagrams, choice of mass or mole basis, and new state variable options Standardizing names in Heat Exchanger models (breaking change from v2.0.0a2): Control Volumes named hot_side and cold_side Ports names hot_side_inlet, hot_side_outlet, cold_side_inlet and cold_side_outlet Config Blocks names hot_side_config and cold_side_config Config arguments for user provided names for each side: hot_side_name and cold_side_name. Updating Keras surrogate tool to use v1.1 of OMLT New prototype API for model initialization (idaes.core.initialization) The new API uses "Model Initializer" objects instead of class methods, allowing for the definition of multiple initialization routines for a single model A number of common, model agnostic initialization routines have also been defined, including initialization from data, block-decomposition and a general hierarchical approach equivalent to the existing method for common unit models New metadata for thermophysical properties – valid_range This can be used to record the range of values over which a property value can be trusted, such as the range of experimental data used to regress parameters A number of new utility functions have been added to check for properties with values outside the valid range and to set bounds based on this metadata Updated construction of balance expressions in Control Volumes to remove unneeded terms In the past, unneeded terms were added as a constant 0 term, however they will now be dropped entirely from the expression This was necessary due to more strict unit checking in the new Pyomo solver writer which no longer ignores 0 terms Updates to metadata for thermophysical properties to better define known properties and units of measurement This results in more strict enforcement of standard naming for thermophysical and reaction properties Users can still define custom properties, but these must be done explicitly using the define_custom_properties() method instead of being implicitly created by add_property() Updated convergence tester utility tool to support definition of benchmark files (JSON format) and comparison of performance to benchmarks Set default iteration limit for IPOPT in IDAES config to 200 iterations Update scaling of example models to work with new Pyomo NL solver writer Improve testing of extensions and examples infrastructure to avoid need for downloading files Updated distillation column to centralize common functionality and remove a number of Pyomo warnings

IDAES↗
DOE OSTI2021

Influence of Surface and Structural Variations in Donor–Acceptor–Donor Sensitizers on Photoelectrocatalytic Water Splitting

Conjugated organic chromophores composed of linked donor (D) and acceptor (A) moieties have attracted considerable attention for photoelectrochemical applications. In this work, we compare the optoelectronic properties and photoelectrochemical performance of two D–A–D structural isomers with thiophene-X-carboxylic acid (X denotes 3 and 2 positions) derivatives and 2,1,3-benzothiadiazole as the D and A moieties, respectively. 5,5′-(Benzo[c][1,2,5]thiadiazole-4,7-diyl)bis(thiophene-3-carboxylic acid), BTD1, and 5,5′-(benzo[c][1,2,5]thiadiazole-4,7-diyl)bis(thiophene-2-carboxylic acid), BTD2, were employed in the study to understand how structural isomers affect surface attachments within chromophore–catalyst assemblies and their influence on charge-transfer dynamics. Crystal structures revealed that varying the position of the −COOH anchoring group causes the molecules to either contort out of a plane (BTD1) or adopt a near-perfect planar conformation (BTD2). BTD1 and BTD2 were co-loaded with either a water oxidation catalyst, [Ru(2,6-bis(1-methylbenzimidazol-2-yl)pyridine)-(4,4′-((HO) 2 OPCH 2 )2-2,2′-bipyridine)(OH 2 )] 2 , RuCt 2+ , or proton reduction catalyst [Ni(P 2 Ph N 2 C 6 H 4 CH 2 PO 3 H 2 ) 2 ] 2+ , NiCt 2+ , on oxide electrodes to facilitate photodriven water splitting reactions. Emission quenching measurements indicate that both BTD1 and BTD2 inject electrons into n-type SnO 2 |TiO 2 electrodes and holes into p-type NiO semiconductors from their respective excited states at high efficiencies >60%. Photocurrent densities of chromophore–catalyst assemblies obtained using linear sweep voltammetry (LSV) show that BTD2-sensitized photoanodes generate significantly more photocurrent than BTD1-sensitized electrodes; however, both exhibit similar performances at the photocathode. Photoelectrocatyltic measurements demonstrate that both BTD1 and BTD2 performed similarly, generating Faradaic efficiencies of 39 and 38% at the anode or 61 and 79% at the cathode. Transient absorption measurements suggest that the differences between the LSV and photoelectrocatalytic measurements result from the differences in quantum yields of the photogenerated redox equivalents, which is also a reflection of the varying metal oxide surface conformation. Our findings suggest that BTD2 should be investigated further in photocathodic studies since it has the structural advantage of being incorporated into diverse types of chromophore–catalyst assemblies.

Chromophores↗
DOE OSTI2023

SQMBox: Interfacing a semiempirical integral library to modular ab initio electronic structure enables new semiempirical methods

Ab initio and semiempirical electronic structure methods are usually implemented in separate software packages or use entirely different code paths. As a result, it can be time-consuming to transfer an established ab initio electronic structure scheme to a semiempirical Hamiltonian. Here we present an approach to unify ab initio and semiempirical electronic structure code paths based on a separation of the wavefunction ansatz and the needed matrix representations of operators. With this separation, the Hamiltonian can refer to either an ab initio or semiempirical treatment of the resulting integrals. We built a semiempirical integral library and interfaced it to the GPU-accelerated electronic structure code TeraChem. Equivalency between ab initio and semiempirical tight-binding Hamiltonian terms is assigned according to their dependence on the one-electron density matrix. The new library provides semiempirical equivalents of the Hamiltonian matrix and gradient intermediates, corresponding to those provided by the ab initio integral library. This enables the straightforward combination of semiempirical Hamiltonians with the full pre-existing ground and excited state functionality of the ab initio electronic structure code. We demonstrate the capability of this approach by combining the extended tight-binding method GFN1-xTB with both spin-restricted ensemble-referenced Kohn–Sham and complete active space methods. We also present a highly efficient GPU implementation of the semiempirical Mulliken-approximated Fock exchange. The additional computational cost for this term becomes negligible even on consumer-grade GPUs, enabling Mulliken-approximated exchange in tight-binding methods for essentially no additional cost.

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

Solitons in lattice field theories via tight-binding supersymmetry

Reflectionless potentials play an important role in constructing exact solutions to classical dynamical systems (such as the Korteweg-de Vries equation), non-perturbative solutions of various large-N field theories (such as the Gross-Neveu model), and closely related solitonic solutions to the Bogoliubov-de Gennes equations in the theory of superconductivity. These solutions rely on the inverse scattering method, which reduces these seemingly unrelated problems to identifying reflectionless potentials of an auxiliary one-dimensional quantum scattering problem. There are several ways of constructing these potentials, one of which is quantum mechanical supersymmetry (SUSY). In this paper, motivated by recent experimental platforms, we generalize this framework to develop a theory of lattice solitons. We first briefly review the classical inverse scattering method in the continuum limit, focusing on the Korteweg-de Vries (KdV) equation and SU(N) Gross-Neveu model in the large N limit. We then generalize this methodology to lattice versions of interacting field theories. Our analysis hinges on the use of trace identities, which are relations connecting the potential of an equation of motion to the scattering data. For a discrete Schrödinger operator, such trace identities had been known as far back as Toda; however, we derive a new set of identities for the discrete Dirac operator. We then use these identities in a lattice Gross-Neveu and chiral Gross-Neveu (Nambu-Jona-Lasinio) model to show that lattice solitons correspond to reflectionless potentials associated with the discrete scattering problem. These models are of significance as they are equivalent to a mean-field theory of a lattice superconductor. To explicitly construct these solitons, we generalize supersymmetric quantum mechanics to tight-binding models. We show that a matrix transformation exists that maps a tight-binding model to an isospectral one which shares the same structure and scattering properties. The corresponding soliton solutions have both modulated hopping and onsite potential, the former of which has no analogue in the continuum limit. We explicitly compute both topological and non-topological soliton solutions as well as bound state spectra in the aforementioned models.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Groundwater and Surface Water Flow (GSFLOW) model files to explore bedrock circulation depth and porosity in Copper Creek, Colorado

This data package contains integrated hydrological model input and output files for Copper Creek, Colorado (24 km2), a tributary of the East River located in the headwaters of the Upper Colorado River Basin. The model code is the U.S. Geological Survey (USGS) Groundwater and Surface Water Flow (GSFLOW) model. The model contains a 100-m grid resolution and a daily timestep. The land surface model is dynamically linked to a three-dimensional groundwater flow model that allows for streamflow gaining and losing conditions. The groundwater model contains 12 model layers and extends 400 m below land surface. The original Copper Creek model was modified to contain geologic layers representing saprolite, shallow bedrock, and deep bedrock. Endmember depth versus hydraulic conductivity relationships and porosity values for fractured crystalline rock are simulated. For the shallow case, median flow depths occur in the shallow saprolite at depths <8 m, while the deep case promotes a median groundwater flow depth of 100 m. With this modeling framework we compare streamflow response to a plausible worst-case drought lasting up to five years. Streamflow metrics of analysis include average streamflow, fraction of stream network that is dry, no-flow duration, average groundwater flow to streams and time to recovery following the drought. Results and implications are presented in a paper submitted to Geophysical Research Letters titled, "The role of bedrock circulation depth and porosity in mountain streamflow response to prolonged drought" by Rosemary WH. Carroll, Andrew H. Manning and Kenneth H Williams. A Readme.txt file provides instructions on how to download all model files and execute each model scenario. In addition to the GSFLOW output/prms/copper_drought.csv file containing daily basin water stores and fluxes (refer to GSFLOW manual) and the output/prms/copper_drought_statvar.dat file with output defined in the gsflow3.control file (refer to GSFLOW Manual), output files also include spatially distributed daily values of total evapotranspiration, canopy evaporation, precipitation, snowfall, infiltration, snow water equivalent, potential evapotranspiration, recharge, sublimation, soil moisture, contributing interflow, water table elevations, changes in groundwater storage, groundwater evapotranspiration, interbasin groundwater flow (limited to the alluvium below the stream outlet), and surface-groundwater exchanges within the river system.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Canonical Transformation of Potential Model Hamiltonian Mechanics to Geometrical Form I

Using the methods of symplectic geometry, we establish the existence of a canonical transformation from potential model Hamiltonians of standard form in a Euclidean space to an equivalent geometrical form on a manifold, where the corresponding motions are along geodesic curves. The advantage of this representation is that it admits the computation of geodesic deviation as a test for local stability, shown in recent previous studies to be a very effective criterion for the stability of the orbits generated by the potential model Hamiltonian. We describe here an algorithm for finding the generating function for the canonical transformation and describe some of the properties of this mapping under local diffeomorphisms. We give a convergence proof for this algorithm for the one-dimensional case, and provide a precise geometric formulation of geodesic deviation which relates the stability of the motion in the geometric form to that of the Hamiltonian standard form. We apply our methods to a simple one-dimensional harmonic oscillator and conclude with a discussion of the relation of bounded domains in the two representations for which Morse theory would be applicable.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Kinetic Monte Carlo simulations of aging in δ -Pu

We have developed a first-passage kinetic Monte Carlo approach for materials aging to investigate the sensitivity of void swelling to model parameters, including helium bubble density and size distribution. In addition to explicitly accounting for the spatial distribution of individual point defects, bubbles, and voids, our approach can simulate total doses equivalent to 100 years of natural aging on statistically representative volumes of materials. This technique enables us to study the effects on swelling and radiation damage evolution due to temperature and dose rate (as altered in artificially aged experiments), differences in effective interaction radii between vacancies and interstitials, and varying defect diffusion activation energies, while providing more detailed information than previous rate-equation based approaches. In conclusion, our results indicate that spatial effects that are not modeled in mean-field rate theories could play a significant role in void swelling initiation and growth for certain regimes of model parameters.

Actinides↗
DOE OSTI2024

Performance Testing of a Moving-Bed Gasifier Using Coal, Biomass, and Waste Plastic Blends with Washed and Unwashed Legacy Coals and Other Waste Fuels to Generate White Hydrogen

The objective of this effort, primarily funded by the United States Department of Energy (DOE), and led by the Electric Power Research Institute, Inc. (EPRI), with support by Hamilton Maurer International (HMI) and Sotacarbo S.p.A. (Sotacarbo), has been to qualify coal, biomass, and plastic waste blends based on performance testing of selected fuel pellet compositions in a pilot-scale updraft moving-bed (UDMB) gasifier. The testing provided relevant data to advance the commercial-scale design of the moving-bed gasifier to be able to successfully use these feedstocks to produce hydrogen. In particular, the effects of waste plastics on feedstock development (i.e., blending and pelletizing) and the resulting products (i.e., syngas compositions, organic condensate production, and ash characteristics) are the focus. The gasifier used for testing is HMI’s moving-bed gasifier, which has been proven capable of gasifying nearly all coal ranks. It has also shown the ability in prior testing work to gasify wood chips (biomass). However, mixtures of these fuels with plastic wastes have not been prepared and gasified together. The three feedstocks were densified and pelletized by California Pellet Mill (CPM) to meet the feedstock size required by Sotacarbo’s 30mm ID UDMB gasifier, under contract to HMI. The technical tasks and results from this two-year research project included: (1) Feed Procurement and Preparation: Nine different tri-fuel pellets were prepared from varying compositions of fresh mined PRB coal, corn stover biomass, and car fluff waste plastics. Tri-fuel pellets were produced by CPM and shipped to Sotacarbo’s test facility in Carbonia, Sardinia, Italy. (2) Test Plan Development: A test plan was created to define the test runs to be performed. The test plan detailed the different UDMB gasification tests to be performed in Sotacarbo’s 12-inch ID pilot scale gasifier, the process monitoring instrumentation used, and the extractive samples recovered for analysis of the total gasification process mass and energy balance. (3) Gasifier Testing: Nine different gasification runs were performed in the pilot-scale gasifier at Sotacarbo using nine different fuel feedstock compositions generated from varying mixtures of PRB coal, biomass, and plastic wastes. The testing generated performance data on gasification reaction efficiency and performance, yielding relevant data for models used to scale up the gasifier design. This task also included work to refurbish and reassemble the pilot gasifier at Sotacarbo and perform a baseline 100% PRB coal run. (4) Data Analysis and Reporting: Review of the data, determination of figures of merit, and interpretation of the results are reported in the project’s final report, published in March 2024. The results show that all tri-fuel pellets gasified well and maintained structural integrity throughout the gasification process. The syngas generated can be shifted to hydrogen by using commercial syngas shifting technologies. (5) High Fidelity computational fluid dynamics (CFD) Simulation: The National Energy Technology Laboratory (NETL) team performed CFD simulations of the UDMB gasifier for two of the tri-fuel pellets gasified in Sotacarbo’s pilot scale gasifier. The kinetic mechanisms for the pyrolysis of each constituent, PRB coal, corn stover biomass, and waste plastics are based on thermogravimetric analysis performed by Sotacarbo. The gasification model was validated by comparing the predicted syngas composition at the exit of the gasifier with the measured syngas composition. In addition, the reactor’s measured internal temperature profile agreed well with the predicted internal reactor temperature profile. These results validate that the model can be used to predict the performance of the updraft moving bed gasifier for different feedstocks and operating conditions. This paper summarizes the results of the completed work in which the pelletizing procedure was validated to ensure the viability of the tri-fuel pellets for the gasification runs performed at Sotacarbo’s 30 mm UDMB gasifier. The gasification performance data from this series of nine runs will enable modeling of a full-scale HMI industrial scale gasifier supporting both combined heat and power, and Hydrogen production from coal (both fresh mined and legacy) combined with various biomass and waste plastics. Additionally, plans and progress on a follow-up project, being executed by the same project team, will be presented. In this project, a total of twenty (20) different feedstocks are being prepared from varying compositions of biomass (both woody biomass and corn stover) with a mixture of legacy coal waste, plastic waste, and refuse-derived fuel (RDF). The testing will provide information on gasification reaction efficiency/performance, yielding relevant data for models used to scale up the gasifier design to 50 megawatt electric (MWe) (equivalent hydrogen production). Tests will also be performed on a bench-scale fluidized-bed gasifier for comparison purposes. The results of this testing will be used to specify the range of feedstock blends that can be successfully gasified as well as quantify gasifier outputs based on specific blends.

08 HYDROGEN↗
DOE OSTI2025

Calcium is associated with specific soil organic carbon decomposition products at Blodgett Forest Research Center, Georgetown, California as analysed with scanning transmission X-ray microscopy carbon near-edge X-ray absorption fine structure spectroscopy

This data is from the paper calcium is associated with specific soil organic carbon decomposition products, published in SOIL. DOI: https://doi.org/10.5194/soil-11-381-2025, 2025.This file contains CSVs with spectral data and bulk soil data and there is no specific program required to open this data. The data includes Scanning transmission X-ray microscopy carbon near-edge X-ray absorption fine structure spectroscopy. data from the measurement of samples from the Whole-soil Warming project, run by the Belowground Biogeochemistry team at Blodgett Forest Research Center, Georgetown, California run by the University of California, Berkeley. It also includes bulk soil chemical properties. The University of California's Blodgett Forest Research Station (Forest) is situated in the Sierra Nevada foothills (1370 m a.s.l.) near Georgetown, California. The samples were collected from here: 38.912013, -120.661469, https://maps.app.goo.gl/291bCJ1zVqUhgktz6. The Forest soils were characterised as Alfisols, which are equivalent to Dystric Cambisols (IUSS Working Group WRB, 2015), and formed in granitic parent materials, in a temperate climate, under thinned, mixed-coniferous forest (Fig. S3; Gaudinski et al., 2009). With these analyses we aimed to answer the question, is calcium associated with a specific type of organic matter enriched in aromatic and phenolic carbon at the microscale in samples from Blodgett Forest Research Center? and how does this specific type of carbon respond to experiments targetted at removing and adding calcium to the soils, specifically cation exchange and incubation after calcium addition? Abstract from the paper can be found below: Calcium (Ca) may contribute to the preservation of soil organic carbon (SOC) in more ecosystems than previously thought. Here we provide evidence that Ca is co-located with SOC compounds that are enriched in aromatic and phenolic groups, across different acidic soil-types and locations with different ecosystem properties, differing in terms of climate, parent material, soil type, and vegetation. In turn, this co-localised fraction of Ca-SOC is removed through cation-exchange, and the association is then only re-established during decomposition in the presence of Ca (Ca addition incubation). Thus, highlighting a causative link between decomposition and the co-location of Ca with a characteristic fraction of SOC. Decomposition increases the relative proportion of negatively charged functional groups, which can increase the propensity for the association between SOC and Ca, and in turn, this association inhibits dissolved organic carbon export or further decomposition. We propose that this mechanism could be driven by Ca hotspots on the microscale shifting local decomposition processes and thereby explaining the colocation of Ca with SOC of a specific composition across different acidic soil environments. Incorporating this biogeochemical process into Earth System Models could improve our understanding, predictions, and management of carbon dynamics in soils, and account for their response to Ca-rich amendments.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Toward Computation-Guided Design of Tunable Organic-Inorganic CdS Quantum Dot Binary Superlattices

Combining the advantages of structural programmability in sequence-defined biomimetic molecules and the controllable packing geometry in nanoparticle superlattices, we demonstrate a self-assembled organic-inorganic superlattice whose structure can be altered with the slightest change in the sequence of the organic counterpart. Here, oleate-coated CdS quantum dots (QDs) form a square-packed superlattice with a 1:1 molar equivalence of a di-block amphiphilic peptoid (Nbrpe6Dig) in chloroform. In contrast, no apparent structure is observed in the organic solvent alone. Based on theoretical evidence, we show that the assembly is a binary superlattice where both the CdS QDs and the peptoids serve as building blocks and further predict a correlation between the superlattice structure and the peptoid sequence. The computationally guided prediction is validated by experiments where superlattice transformation is observed with modified peptoids. The mechanism identified in our work inspires new ways to control and tune organic-inorganic hybrid nanomaterial self-assembly.

Qi, Xin↗
DOE OSTI2020

Experimental investigation of turbulent wake flows in a helically wrapped rod bundle in presence of localized blockages

In nuclear sodium fast reactors, bundles of rods are tightly packed into a triangular lattice, enclosed in a hexagonal duct, and each pin is spirally wrapped with a thin wire. Flow blockages can potentially impact the local flow characteristics and heat transfer mechanisms in the bundle due to its small subchannel size. The effects of the blockage on the flow structures and heat transfer mechanisms are important aspects that require an accurate investigation. In this study, the flow-field characteristics in the vicinity of a blockage located in the exterior subchannel of rod bundles with helically wrapped wires were experimentally investigated. The velocity fields in the exterior subchannel were acquired by applying matched-index-of-refraction and time-resolved particle image velocimetry (TR-PIV) techniques for Reynolds numbers of Re 1 = 4000 and Re 2 = 17 000, i.e., equivalent to Rew 1 = 19 600 and Rew 2 = 83 200, respectively, based on the blockage width. The results from the TR-PIV measurements revealed an arch-shaped vortex with a large flow recirculation and a pair of counter-rotating vortices in the wake region downstream of the blockage, which is commonly observed in the wake flow of bluff bodies. The relative lateral distance and angle between the two vortices decreased when the Reynolds numbers increased. Profiles of maximum turbulence intensity along the shear layers illustrated the transition process including the growth, peak, and decay along the flow direction. From the spectral analysis of the turbulent velocities extracted at points along the shear layer, the Strouhal numbers ( St ) representing the vortex shedding frequency were found to be St = 0.25 and St = 0.56 for the left and right shear layers, respectively. Characteristics of shear layers generated by the blockage in the exterior subchannel were investigated via the two-point cross correlation of fluctuating velocities. The patiotemporal cross correlations of turbulent velocities, computed at points in the region where the left shear layer exhibited rolling effects and vortex breakdowns, were considerably wider and longer. The convection velocity U c was estimated to be ~0.82 U m to 0.93 U m . Proper orthogonal decomposition (POD) analysis was applied to the instantaneous velocity fields to extract the statistically dominant flow structures. It was found that POD modes 2–3 and 4–5 formed the pair modes when the corresponding POD temporal coefficients depicted sinusoidal shapes and exhibited nearly circular orbits in the phase space. Spectral analysis of the POD temporal coefficients confirmed the vortex shedding frequencies detected in the analysis of turbulent velocities.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

The R-Process Alliance: Hunting for gold in the near-UV spectrum of 2MASS J05383296–5904280

Over the past few years, the R-Process Alliance (RPA) has successfully carried out a search for stars that are highly enhanced in elements produced via the rapid neutron-capture (r-) process. In particular, the RPA has identified a number of relatively bright, highly r-process-enhanced (r-II) stars, suitable for observations with the Hubble Space Telescope (HST), facilitating abundance derivation of elements such as gold (Au) and cadmium (Cd). This paper presents the detailed abundances derived for the metal-poor ([Fe/H] = −2.55) highly r-process-enhanced ([Eu/Fe] = +1.29) r-II star 2MASS J05383296–5904280. One-dimensional local thermodynamic equilibrium (LTE) elemental abundances were derived via equivalent width and spectral synthesis using high-resolution high signal-to-noise near-UV HST/STIS and optical Magellan/MIKE spectra. Abundances were determined for 43 elements, including 26 neutron-capture elements. In particular, abundances of the rarely studied elements Nb, Mo, Cd, Lu, Os, Pt, and Au are derived from the HST spectrum. These results, combined with RPA near-UV observations of two additional r-II stars, increase the number of Cd abundances derived for r-process-enriched stars from seven to ten and Au abundances from four to seven. A large star-to-star scatter is detected for both of these elements, highlighting the need for more detections enabling further investigations, specifically into possible non-LTE effects.

chemically peculiar stars↗
DOE OSTI2025

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↗
DOE OSTI2025

The R -Process Alliance: Exploring the cosmic scatter among ten r -process sites with stellar abundances

Context. The astrophysical origin of the rapid neutron-capture process (r-process), responsible for producing roughly half of the elements heavier than iron, remains uncertain. Detailed chemical signatures from the oldest, most metal-poor stars, which act as fossil records of the earliest nucleosynthesis events, can be used to identify the dominant r-process sites. Aims. We present a homogeneous chemical abundance analysis of ten r-process element-enhanced stars. These old and metal-poor stars are strongly enriched in r-process elements with minimal contamination from other nucleosynthetic sources. By focusing on this chemically pure sample, we aim to investigate intrinsic variations in the r-process abundance patterns and explore their implications for the nature and potential diversity of r-process sites. Methods. We performed a detailed chemical abundance analysis of high-resolution, high-signal-to-noise spectra. For each star, we inspected over 1400 individual absorption lines using a combination of equivalent width measurements and spectral synthesis. The analysis was conducted under the assumption of 1D local thermodynamic equilibrium and employing the MOOG radiative transfer code. Results. We derived abundances for 54 chemical species, including 29 neutron-capture (n-capture) elements, covering the full mass range of the r-process abundance pattern. A kinematic analysis reveals that stars likely originated from ten kinematically distinct systems. Based on this assumption, we used the sample to probe the maximum variation expected from ten independent r-process nucleosynthesis events and computed the intrinsic dispersion of each element relative to Zr and Eu for the light and heavy r-process elements, respectively. This exercise resulted in a remarkably low cosmic scatter across the ten r-process sites enriching these stars; for the rare earth and third peak elements, for example, we find σ [La/Eu] = 0.08 and σ [Os/Eu] = 0.11 dex, while the scatter between light and heavy elements, σ [Zr/Eu] , is slightly higher at 0.18 dex. Conclusions. The elemental abundance patterns across the ten independent r-process sites show remarkably small cosmic dispersions. This minimal dispersion suggests a high degree of uniformity in r-process yields across diverse astrophysical environments.

Astronomy and AstroPhysics↗