Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “RADIAL DISTRIBUTION”

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Distributions in CFT. Part II. Minkowski space

CFTs in Euclidean signature satisfy well-accepted rules, such as the convergent Euclidean OPE. It is nowadays common to assume that CFT correlators exist and have various properties also in Lorentzian signature. Some of these properties may represent extra assumptions, and it is an open question if they hold for familiar statistical-physics CFTs such as the critical 3d Ising model. Here we consider Wightman 4-point functions of scalar primaries in Lorentzian signature. We derive a minimal set of their properties solely from the Euclidean unitary CFT axioms, without using extra assumptions. We establish all Wightman axioms (temperedness, spectral property, local commutativity, clustering), Lorentzian conformal invariance, and distributional convergence of the s-channel Lorentzian OPE. This is done constructively, by analytically continuing the 4-point functions using the s-channel OPE expansion in the radial cross-ratios ρ, ρ¯. We prove a key fact that |ρ|, |ρ¯| < 1 inside the forward tube, and set bounds on how fast |ρ|, |ρ¯| may tend to 1 when approaching the Minkowski space.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Cluster environment quenches the star formation of low-mass satellite galaxies from the inside-out

ABSTRACT Environment plays a critical role in the star formation history of galaxies. Tidal and hydrodynamical stripping, prominent in cluster environment, can remove the peripheral gas of galaxies and star formation may thus be environmentally suppressed from the outside-in. We revisit the environmental dependence of the radial gradient of specific star formation rate (sSFR) profile. We probe the radial gradient by using the archival spectral indices D4000n and HδA measured from SDSS fibre spectra, to indicate central sSFR and the total sSFR from fitting the spectral energy distribution. Despite the low spatial resolution, the wealth of SDSS data allows to disentangle the dependences on stellar mass, sSFR, and environment. We find that low-mass satellite galaxies in the mass range $9\lt \mathrm{log}\, \mathcal {M}_{\star }/\mathcal {M}_{\odot }\lt 9.8$ on average quench in more inside-out pattern compared to isolated galaxies matched in mass, sSFR, and fibre coverage. This environmental effect is particularly strong for galaxies below the star formation main sequence, and peaks for those in the core of massive clusters where the phase-space diagram reveals clear links between the inside-out quenching and orbital properties. Our results suggest that both tidal and hydrodynamical interactions in cluster environment suppress the star formation of satellites mainly from the inside-out. As accreted gas of low angular momentum from hot gas haloes is an important source for replenishing central gas reservoir, we discuss how gas stripping in clusters may lead to starvation and cause inside-out quenching when the outer star-forming discs are not significantly affected.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Influence of Shear Strength Assumptions on BISON Debonding Simulations

Accurately predicting the thermomechanical response of buffer–IPyC debonding in TRISO fuel particles requires reliable mechanical property inputs for each coating layer, particularly the normal and shear strengths that influence interlayer delamination and stress concentrations. Micro tensile testing of AGR-2 fuel particles provided experimentally measured normal strengths for the buffer, IPyC, and buffer–IPyC interface; however, shear strength was not measured. As a result, BISON simulations of interface debonding must rely on assumed shear strength values, typically estimated as 20–40% of the measured ultimate tensile strength. This study evaluates how these assumed shear strength values influence cohesive zone model (CZM) predictions of buffer–IPyC separation in AGR 2 TRISO particles. Using micro tensile data from three AGR 2 compacts (2 1 3, 5 1 3, and 6 3 3), BISON simulations were performed with multiple shear strength assumptions to quantify their effect on radial and tangential stress evolution, debonding, and gap propagation. The results show that shear strength is a high sensitivity parameter: increasing the assumed shear strength significantly alters the stress distribution at the buffer–IPyC junction, shifts the predicted debonding location, and changes the extent of partial gap formation. While normal strength controls the initiation of interface separation, shear strength strongly influences the mode mixity of the failure process and the resulting stress concentrations transmitted to the IPyC and SiC layers. These findings highlight a critical gap in current TRISO mechanical characterization. Without experimentally measured shear strength, BISON simulations must rely on approximations that introduce uncertainty into predictions of coating layer integrity and fission product barrier performance. Future fuel qualification campaigns should therefore consider measurement of shear strength at the interlayer interfaces to reduce model uncertainty and improve the fidelity of TRISO fuel performance simulations.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Discovery of Extended Tidal Tails around the Globular Cluster Palomar 13

We use photometry from the DECam Legacy Survey to detect candidate tidal tails extending ~5° on either side of the Palomar 13 globular cluster. The tails are aligned with the proper motion of Palomar 13 and are consistent with its old, metal-poor stellar population. We identify three RR Lyrae stars (RRLs) that are plausibly associated with the tails, in addition to four previously known in the cluster. From these RRLs, we find that the mean distance to the cluster and tails is 23.6 ± 0.2 kpc and estimate the total (initial) luminosity of the cluster to be , consistent with previous claims that its initial luminosity was higher than its current luminosity. Combined with previously determined proper motion and radial velocity measurements of the cluster, we find that Palomar 13 is on a highly eccentric orbit (e ~ 0.8) with a pericenter of ~9 kpc and an apocenter of ~69 kpc, and a recent pericentric passage of the cluster ~75 Myr ago. We note a prominent linear structure in the interstellar dust map that runs parallel to the candidate tidal features, but conclude that reddening due to dust is unlikely to account for the structure that we observe. If confirmed, the Palomar 13 stellar stream would be one of very few streams with a known progenitor system, making it uniquely powerful for studying the disruption of globular clusters, the formation of the stellar halo, and the distribution of matter within our Galaxy.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Bayesian Cross-matching of High Proper-motion Stars in Gaia DR2 and Photometric Metallicities for ∼1.7 million K and M Dwarfs

We present a Bayesian method to cross-match 5,827,988 high proper-motion Gaia sources (μ > 40 mas yr{sup −1}) to various photometric surveys: Two Micron All Sky Survey, AllWISE data release from the Wide-field Infrared Explorer (WISE) mission, Galaxy Evolution Explorer, Radial Velocity Experiment, Sloan Digital Sky Survey, and Panoramic Survey Telescope and Rapid Response System (Pan-STARRS). To efficiently associate these objects across catalogs, we develop a technique that compares the multidimensional distribution of all sources in the vicinity of each Gaia star to a reference distribution of random field stars obtained by extracting all sources in a region on the sky displaced 2′. This offset preserves the local field stellar density and magnitude distribution, allowing us to characterize the frequency of chance alignments. The resulting catalog with Bayesian probabilities >95% has a marginally higher match rate than current internal Gaia data release 2 (DR2) matches for most catalogs. However, a significant improvement is found with Pan-STARRS, where ∼99.8% of the sample within the Pan-STARRS footprint is recovered, as compared to a low ∼20.8% in Gaia DR2. Using these results, we train a Gaussian process regressor to calibrate two photometric metallicity relationships. For dwarfs of 3500 < T{sub eff} < 5280 K, we use metallicity values of 4378 stars from the Apache Point Observatory Galactic Evolution Experiment and Hejazi et al. to calibrate the relationship, producing results with a 1σ precision of 0.12 dex and few systematic errors. We then indirectly infer the metallicity of 4018 stars with 2850 < T{sub eff} < 3500 K, which are wide companions of primaries whose metallicities are estimated with our first regressor, to produce a relationship with a 1σ precision of 0.21 dex and significant systematic errors. Additional work is needed to better remove unresolved binaries from this second sample to reduce these systematic errors.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

MPACT Theory Manual (V.4.3)

This manual presents the theory underlying the three-dimensional (3D) whole-core, pin-resolved neutron transport calculation methodologies employed in the MPACT code. MPACT’s primary goal is to provide accurate sub-pin power distributions in a computationally efficient manner. To accomplish this, MPACT offers several different transport method options. The 2D/1D method, in which 3D problems are decomposed into an axial stack of radial “slices,” is currently the most commonly used. In this option, two-dimensional (2D) planar solutions are provided by the method of characteristics (MOC), and axial solutions are provided via one-dimensional (1D) approximate diffusion, or P3 solutions. The radial and axial solutions are coupled by (i) axial and radial transverse leakages, and (ii) a global 3D coarse mesh finite difference (CMFD) solve, which provides both acceleration and stability to the solution iteration scheme. The subsequent chapters of this manual present a range of topics, including MOC, CMFD, axial nodal transport solvers, 2D/1D, self-shielding, depletion, thermal-hydraulics, and transient methods. Development of the underlying theory of the 2D/1D method is presented, diagrams are included to highlight important algorithmic flow, and important concepts are discussed as appropriate. This manual is intended to be self-sufficient, but references to published articles and other materials are included for further reading. MPACT is a relatively new code, with new capabilities and many computational methods that did not exist until recently. This manual is an even newer document, with chapters written by several different code contributors working under time constraints. For these reasons, the manual is not yet complete and finalized. Future versions of this manual will address current deficiencies.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

AGGREGATE: dAta-driven modelinG preservinG contRollable dEr for outaGe mAnagemenT and rEsiliency (Final Report)

The AGGREGATE project team successfully developed and validated various modules for outage management. Brief summaries of each module are provided to showcase their strength for outage management and restoration for a distribution system with a high penetration of connected distribution energy resources (DERs). In recent years, inverter-based DERs have been widely deployed in distribution system. A most of behind-the-meter (BTM) solar power generation is not visible to the utility. The data-driven DER and load estimation modules are using machine learning (ML) and artificial intelligence (AI) to manage this issue, which provides an opportunity for distribution system operators (DSOs) to operate systems and make decisions in real-time for a distribution system with a high penetration of DERs deployed. Also, the estimated DER and true load can be further leveraged in network aggregation and cold-load pick up estimation for reducing the computing complexity and providing for fast restoration. After load demand and DER power generations have been estimated, the information will support topology and state estimation (SE). The topology estimation module demonstrated the viability of mixed integer linear programming (MILP) formulation to estimate the most likely operational radial topology and outage sections using power flow measurements, historical/estimated load and DERs data and smart meter ping measurements. Formulation includes continuous (power flow, load and DERs data) and binary measurements (smart meter ping measurements) in a single formulation. Errors in continuous data and binary data are modeled as normal distribution and Bernoulli distribution, respectively. In the future distribution grid, the power injection from controllable DERs will be essential for efficient and resilient grid operation. However, determining the optimal DER injections and restoration actions is dependent on knowledge of the system states. State estimation (SE), already the cornerstone of transmission energy management systems, will become commonplace in distribution management systems as more measurements become available from deployment of automated metering infrastructure (AMI). Observability analysis is the first step in SE, as it determines the sufficiency of the available measurements for accurately estimating the current system states. A new type of pseudo-measurement called a Correlational Measurement (CM) is introduced in this module, to enhance the observability of the system to enable more accurate SE. CMs encapsulate knowledge of correlation between demand patterns for similar classes of loads as well as injection patterns for same-technology renewable DERs. During grid contingency scenarios, DERs have been traditionally disconnected, without any fault ride-through capabilities. However, with new regulations and better technology, it is feasible for these resources to contribute to the grid’s restoration after an adverse event and hence enhance resilience. The controllability module proposes a two-step restoration scheme for the power system restoration process by leveraging additional degrees of freedom in power electronics interfaced DERs for mitigating voltage problems. In a resilience mode without the utility system, the distribution grid relies on DERs to serve critical load. In such a severe event with multiple faults on the distribution feeders, actuation of various protective devices (PDs) divides the distribution system into electrical islands. The undetected actuated PDs due to fault current contributions from DERs can delay the restoration process, thereby reducing the system resilience. The Advanced Outage Management (AOM) and the Advanced Feeder Restoration (AFR) modules developed in this project provide improved system resilience with multiple DERs. AOM identifies the faulted sections and actuated PDs in a distribution system with DERs by incorporating smart meter data. The most credible outage scenario including fault locations, PD actuations, and fault indicator (FI) failures is identified by a set of binary integer linear programming incorporating hypotheses. The AFR module serves to restore a distribution system with available energy resources taking into consideration the availability of utility sources and DERs. By partitioning the system into islands, critical load will be served with the available generation resources within islands based on the solution of a MILP. When the utility systems become available, the optimal path will be determined by a spanning tree search algorithm that reconnects these islands back to substations and restores the remaining load. The transmission and distribution (T&D) co-simulation module was used to validate the effect of a control action performed on the distribution side assets as it propagates to the transmission side. This ensures that the control action performed results in a feasible operating point on both the transmission and the distribution system. In addition to validation, the team used the T&D co-simulation module to demonstrate how distribution system assets can be used to mitigate issues on the transmission system. Specifically, the team demonstrated that appropriate switching operations on the distribution side can alleviate the line overload condition on the transmission side without causing new operational constraint violations.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

DC electric field generation and distribution in magnetized plasmas

Very large DC and AC electric fields cannot be sustained between conducting electrodes because of volume gas breakdown and/or surface field emission. However, very large potential fields are now routinely generated in plasma structures, such as laser generated wake in unmagnetized plasmas. In magnetized plasmas, large DC fields can also be sustained and controlled perpendicular to the magnetic field, but the metallic end plates limiting the plasma, terminating the magnetic field lines, and usually providing the voltage drop feed between the field lines impose severe restrictions on the maximum field. However, it is shown that very large radial DC voltage drops can be sustained by injecting waves of predetermined frequencies and wave vectors, traveling along the azimuthal direction of an axially magnetized plasma cylinder, or by injecting fast neutral particles beams along this azimuthal direction. The large conductivity along the magnetic field lines and the small conductivity between the field lines then distribute this voltage drop. Furthermore, the global power balance and control parameters of wave and beam generated large DC electric fields in magnetized plasmas are identified, described, and analyzed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Modeling of the high column density systems in the Lyman-Alpha forest

The Lyman-α forests observed in the spectra of high-redshift quasars can be used as a tracer of the cosmological matter density to study baryon acoustic oscillations (BAO) and the Alcock-Paczynski effect. Extraction of cosmological information from these studies requires modeling of the forest correlations. While the models depend most importantly on the bias parameters of the intergalactic medium (IGM), they also depend on the numbers and characteristics of high-column-density systems (HCDs) ranging from Lyman-limit systems with column densities log N HI /1cm -2 > 17 to damped Lyman-α systems (DLAs) with log N HI /1cm -2 > 20.2. These HCDs introduce broad damped absorption characteristic of a Voigt profile. Consequently they imprint a component on the power spectrum whose modes in the radial direction are suppressed, leading to a scale-dependent bias. Using mock data sets of known HCD content, we test a model that describes this effect in terms of the distribution of column densities of HCDs, the Fourier transforms of their Voigt profiles and the bias of the halos containing the HCDs. Our results show that this physically well-motivated model describes the effects of HCDs with an accuracy comparable to that of the ad-hoc models used in published forest analyses. We also discuss the problems of applying the model to real data, where the HCD content and their bias is uncertain.

Lyman alpha forest↗
DOE OSTI2026

Localized Interactions in Neutrino Simulations

The Deep Underground Neutrino Experiment (DUNE) requires precise modeling of neutrino--nucleus interactions to achieve its neutrino-oscillation measurement goals. GENIE, the Monte Carlo event generator used by DUNE, exhibits a known discrepancy with MicroBooNE measurements of transverse kinematic imbalance (TKI): the data display a larger high-TKI tail while maintaining a peak similar to that predicted by the baseline GENIE model. Previous variations of final-state interaction (FSI) strength affected both the peak and the tail and therefore did not resolve the discrepancy. This work investigates whether inconsistencies between local and global treatments of intranuclear physics contribute to the observed mismodeling. The GENIE FSI routines were modified to use the particle position when generating nucleons, thereby introducing a local-density treatment in the hA and hN intranuclear models for their 2018 and 2025 implementations. Meson-exchange-current (MEC) localization was also tested for the hN 2018 model. Comparisons of the intranuclear scattering-center momentum and its radial dependence confirm that the localization was implemented as intended. Localization produces modest changes in the proton momentum spectrum, primarily at low momentum, but only a small change in the TKI distribution between approximately \SI{0.25}{\giga\electronvolt} and \SI{0.40}{\giga\electronvolt}. These changes are insufficient to account for the discrepancy with MicroBooNE data. Although a localized treatment improves the internal consistency of the GENIE model, the origin of the TKI discrepancy remains unresolved.

Bulla, Braden [Unlisted, US, IL]↗
DOE OSTI2025

Optimal Network Reconfiguration and Scheduling With Hardware-in-the-Loop Validation for Improved Microgrid Resilience

With the increased occurrence of various major extreme weather events, power outages and prompt power system restorations have recently drawn more attention to the resilience and recovery of power systems. From the perspective of a more resilient power delivery at the distribution grid, system restoration using network topology reconfiguration together with optimal scheduling of distributed energy resources are adopted in this paper. The proposed optimization model aims at minimizing the total load shedding cost and other operational costs, in which linearized topological constraints borrowed from graph theory and linearized DistFlow models are respectively used to maintain the radial network topology and power flow balance after system contingencies. To demonstrate the applicability of the proposed strategy, a real-world case study of a networked three-microgrid system in Adjuntas, Puerto Rico, is used with the consideration of different independent/interconnected microgrid scenarios, contingencies, and fairness settings. Furthermore, hardware-in-the-loop testing is conducted for the same three-microgrid network, where the closely matched results with the simulated ones have validated the effectiveness of the proposed restoration strategy, which is now ready to move one step forward towards field deployment. Finally, to test the proposed restoration strategy in a larger networked system, the modified IEEE-33 bus test distribution system is considered, and the results show a more resilient power delivery for critical loads under three and four line outages.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2026

Signatures of a Tidally Induced Spiral Arm at the Anticenter of the Milky Way and a Kinematically Extended Anticenter Stream Using DESI Data Release 2

Using the Dark Energy Spectroscopic Instrument (DESI) Milky Way Survey, we examine the six-dimensional space of the anticenter region of the Milky Way stellar disk (150° < Galactic longitude < 220°) using 61,883 main-sequence turnoff stars. We focus on two well-known stellar overdensities in the anticenter: the Monoceros Ring (MRi) and Anticenter Stream (ACS). We find that the MRi overdensity has kinematic signatures consistent with a tidally induced spiral arm, a type of dynamic spiral arm created by an interaction with a satellite galaxy, most likely the Sagittarius dwarf spheroidal galaxy (Sgr). We use the kinematics of the MRi to calculate the two most recent passage times of Sgr, finding 0.25 ± 0.09 Gyr and 1.10 ± 0.23 Gyr from the present day. We validate that the ACS is kinematically decoupled from the MRi because they are moving in opposite radial and vertical directions. We find that the kinematics associated with the ACS extends beyond our defined overdensity. The features we see in the ACS region are likely part of a broader distribution of stars with the same kinematic signature as detected in other places, like the vertical wave in the outer disk and phase spiral.

Lambert, Mika [University of California, Santa Cru↗
DOE OSTI2024

The XRISM first-light observation: Velocity structure and thermal properties of the supernova remnant N 132D

Abstract We present an initial analysis of the X-Ray Imaging and Spectroscopy Mission (XRISM) first-light observation of the supernova remnant (SNR) N 132D in the Large Magellanic Cloud. The Resolve microcalorimeter has obtained the first high-resolution spectrum in the 1.6–10 keV band, which contains K-shell emission lines of Si, S, Ar, Ca, and Fe. We find that the Si and S lines are relatively narrow, with a broadening represented by a Gaussian-like velocity dispersion of $\sigma _v \sim 450$ km s$^{-1}$. However, the Fe He$\alpha$ lines are substantially broadened with $\sigma _v \sim 1670$ km s$^{-1}$. This broadening can be explained by a combination of the thermal Doppler effect due to the high ion temperature and the kinematic Doppler effect due to the SNR expansion. Assuming that the Fe He$\alpha$ emission originates predominantly from the supernova ejecta, we estimate the reverse shock velocity at the time when the bulk of the Fe ejecta were shock heated to be $-1000 \lesssim V_{\rm rs}$ (km s$^{-1}$) $\lesssim 3300$ (in the observer frame). We also find that Fe Ly$\alpha$ emission is redshifted with a bulk velocity of $\sim 890$ km s$^{-1}$, substantially larger than the radial velocity of the local interstellar medium surrounding N 132D. These results demonstrate that high-resolution X-ray spectroscopy is capable of providing constraints on the evolutionary stage, geometry, and velocity distribution of SNRs.

individual objects↗
DOE OSTI2026

Studies of Localization Effects on Transverse Kinematic Imbalance in the GENIE Generator

The forthcoming Deep Underground Neutrino Experiment (DUNE) requires precise modeling of neutrino--nucleus interactions to achieve its neutrino-oscillation measurement goals. GENIE, the Monte Carlo event generator used by DUNE and other Fermilab-based experiments as a central value model for neutrino-nucleus cross sections, exhibits known discrepancies in comparison to MicroBooNE cross section data, especially when considering neutrino scattering events' transverse kinematic imbalance (TKI). For quasielastic-like measurements, the data show a larger transverse missing momentum tail while maintaining a peak similar to that predicted by the baseline GENIE model. Previous attempts at a solution include variations of final-state interaction (FSI) strengths, incorrectly decreasing both the peak and tail, thus leaving the discrepancy unresolved. This work investigates whether inconsistencies between local and global treatments of intranuclear physics contribute to the observed mismodeling. In this work, GENIE FSI routines have been modified to directly utilize intranuclear particle positions, thereby introducing a newly localized momentum treatment in correlation with the nominal density in both the hA and hN intranuclear cascade models for their 2018 and 2025 variants. Meson-exchange-current (MEC) localization was also tested for the hN 2018 model. Comparisons of the intranuclear nucleons' ("scattering center") momentum and its radial dependence validate the localization implementation. However, while such consistent localization produces modest changes in the proton momentum spectrum at low-to-medium momenta, only small changes are observed in the TKI distributions. These changes are insufficient to account for the discrepancy with MicroBooNE data. Although a localized treatment improves the internal consistency of the GENIE model, the origin of the TKI discrepancy remains unresolved by such a solution.

Bulla, Braden [Ctr. Coll., Danville; Fermilab]↗
DOE OSTI2022

Dark and luminous mass components of Omega Centauri from stellar kinematics

Here we combine proper motion data from Gaia EDR 3 and HST with line-of-sight velocity data to study the stellar kinematics of the ω Cen globular cluster. Using a steady-state, axisymmetric dynamical model, we measure the distribution of both the dark and luminous mass components. Assuming both Gaussian and Navarro–Frenk–White mass profiles, depending on the data set, we measure an integrated mass of ≲10 6 M ⊙ within the ω Cen half-light radius for a dark component that is distinct from the luminous stellar component. For the HST and radial velocity data, models with a non-luminous mass component are strongly statistically preferred relative to a stellar mass-only model with a constant mass-to-light ratio. While a compact core of stellar remnants may account for a dynamical mass up to ~5 x 10 5 M ⊙ , they likely cannot explain the higher end of the range. This leaves open the possibility that this non-luminous dynamical mass component comprises non-baryonic dark matter. In comparison to the dark matter distributions around dwarf spheroidal galaxies, the ω Cen dark mass component is much more centrally concentrated. Interpreting the non-luminous mass distribution as particle dark matter, we use these results to obtain the J-factor, which sets the sensitivity to the annihilation cross-section. For the data sets considered, the range of median J-factors is ~10 22 -10 24 GeV 2 cm -5 , which is larger than that obtained for any dwarf spheroidal galaxy.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Measurement of Transport Properties of Woody Biomass Feedstock Particles Before and After Pyrolysis by Numerical Analysis of X-Ray Tomographic Reconstructions

Lignocellulosic biomass has a complex, species-specific microstructure that governs heat and mass transport during conversion processes. A quantitative understanding of the evolution of pore size and structure is critical to optimize conversion processes for biofuel and bio-based chemical production. Further, improving our understanding of the microstructure of biochar coproduct will accelerate development of its myriad applications. This work quantitatively compares the microstructural features and the anisotropic permeabilities of two woody feedstocks, red oak and Douglas fir, using X-ray computed tomography (XCT) before and after the feedstocks are subjected to pyrolysis. Quantitative analysis of the three-dimensional (3D) reconstructions allows for direct calculations of void fractions, pore size distributions and tortuosity factors. Next, 3D images are imported into an immersed boundary based finite volume solver to simulate gas flow through the porous structure and to directly calculate the principal permeabilities along longitudinal, radial, and tangential directions. The permeabilities of native biomass are seen to differ by three to four orders of magnitude in the different principal directions, but we find that this anisotropy is substantially reduced in the biochar formed during pyrolysis. The quantitative transport properties reported here enhance the ability of pyrolysis simulations to account for feedstock-specific effects and thereby provide a useful touchstone for the biorefining community.

09 BIOMASS FUELS↗
DOE OSTI2020

BD+14°3061: A Luminous Yellow Post-AGB Star in the Galactic Halo

I report the discovery that the ninth-mag Galactic-halo star BD+14°3061 is a member of the rare class of luminous metal-poor “yellow post-AGB” stars. Its Gaia DR2 parallax implies an absolute magnitude of M V = -3.44 ± 0.27, and it is a very high-velocity star moving in a retrograde Galactic orbit. BD+14°3061 is a field analog of the half-dozen yellow PAGB stars known in Galactic globular clusters, which have closely similar absolute magnitudes. These objects are the visually brightest members of old stellar populations; their apparently narrow luminosity function makes them potentially useful as Population II standard candles. The spectral-energy distribution of BD+14°3061 out to 22 μm shows no evidence for circumstellar dust. The star is a low-amplitude semiregular pulsating variable, with typical periods of 30–32 days. A radial-velocity study suggests that it is a spectroscopic binary with a period of 429.6 days, making it similar to known binary yellow PAGB stars such as HD 46703 and BD+39°4926.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Transient Modeling and Simulation of a Generic Stable Salt Reactor

A SAM system-level model of a generic stable salt reactor has been developed to investigate thermal-hydraulic behavior and safety performance under steady and transient conditions. The model integrates information generated from a reactor physics analysis using PROTEUS and PERSENT, and a computation fluid dynamics (CFD) analysis using STAR-CCM+. A loose, iterative coupling scheme between PROTEUS and SAM is implemented to calculate the equilibrium power and temperature distributions in the steady-state critical core condition. The converged steady-state model is then used in PERSENT to calculate the four reactivity feedback temperature coefficients (Doppler, fuel density, coolant density, and core radial expansion) and kinetic parameters that are needed in SAM to model the temperature feedback effects in transient simulations. Within the fully enclosed liquid fuel pins, natural convection is the dominant heat transfer mechanism. The STAR-CCM+ model of the fuel pin considers conjugate heat transfer from the liquid fuel salt to the pin cladding and external reactor coolant. The CFD results of the axial and radial temperature profiles are used to empirically determine an effective fuel salt thermal conductivity in the SAM fuel pin model so that the temperatures predicted by the SAM model match as closely as possible the CFD results. In the central region of the fuel pin, the effective thermal conductivity is as high as similar to 60 times the physical fuel salt thermal conductivity. The whole-plant SAM model is then used to simulate an unprotected station blackout transient. The results of this simulation showed that the large negative fuel axial expansion reactivity feedback reduces fission power to similar to 2.4% nominal power. The core is cooled by natural circulation, which removes heat in the core to the emergency heat removal system, and ultimately, to the ambient. However, peak fuel salt and cladding temperatures can potentially reach as high as 1500 K, albeit briefly, if the shutdown mechanism fails to operate.

stable salt reactor; transient simulations; system↗