Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “normal stars”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Multiphysics Design Optimization and Additive Manufacturing of Nuclear Components (Final CRADA Report - Executive Summary)

Westinghouse Electric Company (WEC) actively participated in the advancement of the nuclear fuel and reactor design space and requested the help of Oak Ridge National Laboratory (ORNL) in the creation of a new design tool set. This report details the creation of a collection of software tool sets that are linked together to collectively assist WEC design engineers in developing novel ideas outside the normal scope of traditional nuclear fuel and reactor design formulas. Specifically, Siemens HEEDS, a design space exploration and parametric optimization software, monitored and changed parameters in a collection of softwares to meet the team’s objective. The HEEDS parametric optimization method, SHERPA, was developed to control the Siemens NX CAD platform to adjust the native CAD of a hexahedral spacer grid. This new geometry can be used to execute a topological design optimization by the NX Topology software add-in. The resulting geometry is additively manufacturable. This topological optimization occurred twice—once on the spacer grid’s spring, and once on the dimple geometry. These new geometries were imported by Siemens’ STAR-CCM+, a multiphysics structural and fluid dynamic computational solver in which the spring geometry is deflected to match the rod insertion configuration. Along with the dimple geometry, this new deflected spring was used to complete a hydraulic assessment of a single-unit cell comprising one rod, one spring, and two dimples. The HEEDS SHERPA algorithm ranks the design based on the final mass of the unit cell and the hydraulic pressure drop performance. The ORNL team demonstrated the ability to use this software and provided engineering judgement to apply modern aerospace aerodynamic design. The effort has been focused on thinking outside the conventional design space and redesigning a spacer grid to perform beyond the WEC set objectives. Furthermore, the ORNL team also demonstrated that the HEEDS optimization routine can independently develop a design that meets the WEC design goals. Although these designs were at a low technology readiness level, their demonstration confirmed the team’s capability to create novel advanced nuclear concepts.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Bayesian mixture model approach to quantifying the empirical nuclear saturation point

The equation of state (EOS) in the limit of infinite symmetric nuclear matter exhibits an equilibrium density, $n_0 \approx 0.16 \, \mathrm{fm}^{-3}$, at which the pressure vanishes and the energy per particle attains its minimum, $E_0 \approx -16 \, \mathrm{MeV}$. Although not directly measurable, the nuclear saturation point $(n_0,E_0)$ can be extrapolated by density functional theory (DFT), providing tight constraints for microscopic interactions derived from chiral effective field theory (EFT). However, when considering several DFT predictions for $(n_0,E_0)$ from Skyrme and Relativistic Mean Field (RMF) models together, a discrepancy between these model classes emerges at high confidence levels that each model prediction's uncertainty cannot explain. How can we leverage these DFT constraints to rigorously benchmark nuclear saturation properties of chiral interactions? To address this question, we present a Bayesian mixture model that combines multiple DFT predictions for $(n_0,E_0)$ using an efficient conjugate prior approach. The inferred posterior distribution for the saturation point's mean and covariance matrix follows a Normal-inverse-Wishart class, resulting in posterior predictives in the form of correlated, bivariate $t$-distributions. The DFT uncertainty reports are then used to mix these posteriors using an ordinary Monte Carlo approach. At the 95\% credibility level, we estimate $n_0 \approx 0.157 \pm 0.010 \, \mathrm{fm}^{-3}$ and $E_0 \approx -15.97 \pm 0.40 \, \mathrm{MeV}$ for the marginal (univariate) $t$-distributions. Combined with chiral EFT calculations of the pure neutron matter EOS, we obtain bivariate normal distributions for the nuclear symmetry energy and its slope parameter evaluated at $n_0$: $S_v \approx 32.0 \pm 1.1 \, \mathrm{MeV}$ and $L\approx 52.6\pm 8.1 \, \mathrm{MeV}$ (95\%), respectively. Furthermore, our Bayesian framework is publicly available, so practitioners can readily use and extend our results.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Bayesian estimation of the $S$ factor and thermonuclear reaction rate for 16 O(p, γ) 17 F

The 16 O(p, γ) 17 F reaction is the slowest hydrogen-burning process in the CNO mass region. Its thermonuclear rate sensitively impacts predictions of oxygen isotopic ratios in a number of astrophysical sites, including AGB stars. The reaction has been measured several times at low bombarding energies using a variety of techniques. The most recent evaluated experimental rates have a reported uncertainty of about 7.5% below 1 GK. However, the previous rate estimate represents a best guess only and was not based on rigorous statistical methods. We apply a Bayesian model to fit all reliable 16 O(p, γ) 17 F cross section data, and take into account independent contributions of statistical and systematic uncertainties. The nuclear reaction model employed is a single-particle potential model involving a Woods-Saxon potential for generating the radial bound state wave function. The model has three physical parameters, the radius and diffuseness of the Woods-Saxon potential, and the asymptotic normalization coefficients (ANCs) of the final bound state in 17 F. Here, we find that performing the Bayesian S -factor fit using ANCs as scaling parameters has a distinct advantage over adopting spectroscopic factors instead. Based on these results, we present the first statistically rigorous estimation of experimental 16 O(p, γ) 17 F reaction rates, with uncertainties (±4.2%) of about half the previously reported values.

6 ≤ A ≤ 19↗
DOE OSTI2021

Constraining the nonanalytic terms in the isospin-asymmetry expansion of the nuclear equation of state

In this report we examine the properties of the isospin-asymmetry expansion of the nuclear equation of state from chiral two- and three-body forces. We focus on extracting the high-order symmetry energy coefficients that consist of both normal terms (occurring with even powers of the isospin asymmetry) as well as terms involving the logarithm of the isospin asymmetry that are formally nonanalytic around the expansion point of isospin-symmetric nuclear matter. These coefficients are extracted from numerically precise perturbation theory calculations of the equation of state coupled with a new set of finite difference formulas that achieve stability by explicitly removing the effects of higher-order terms in the expansion. We consider contributions to the symmetry energy coefficients from both two- and three-body interactions. It is found that the coefficients of the logarithmic terms are generically larger in magnitude than those of the normal terms from second-order perturbation theory diagrams, but overall the normal terms give larger contributions to the ground state energy. The high-order isospin-asymmetry terms are especially relevant at large densities where they affect the proton fraction in β-equilibrium matter, and in particular we find that at twice saturation density they can reduce the proton fraction by up to 0.02.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

The predicted properties of helium-enriched globular cluster progenitors at high redshift

Globular cluster progenitors may have been detected by Hubble Space Telescope, and are predicted to be observable with James Webb Space Telescope (JWST) and ground-based extremely large telescopes with adaptive optics. This has the potential to elucidate the issue of globular cluster formation and the origins of significantly helium-enriched subpopulations, a problem in Galactic astronomy with no satisfactory theoretical solution. Given this context, we use model stellar tracks and isochrones to investigate the predicted observational properties of helium-enriched stellar populations in globular cluster progenitors. We find that, relative to helium-normal populations, helium-enriched (ΔY = +0.12) stellar populations similar to those inferred in the most massive globular clusters, are expected, modulo some rapid fluctuations in the first ~30 Myr, to be brighter and redder in the rest frame. At fixed age, stellar mass, and metallicity, a helium-enriched population is predicted to converge to being ~0.40 mag brighter at $\lambda \approx 2.0\, {\mu \rm m}$, and to be 0.30-mag redder in the JWST–NIRCam colour (F070W - F200W), and to actually be fainter for $\lambda \lesssim 0.50 \, {\mu \rm m}$. Separately, we find that the time-integrated shift in ionizing radiation is a negligible $\sim \!5{{\ \rm per\ cent}}$, though we show that the Lyman-α escape fraction could end up higher for helium-enriched stars.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Stellar Mass and Stellar Mass-to-light Ratio–Color Relations for Low Surface Brightness Galaxies

We estimate the stellar mass for a sample of low surface brightness galaxies (LSBGs) by fitting their multiband spectral energy distributions (SEDs) to the stellar population synthesis model. The derived stellar masses (log M {sub *}/M {sub ⊙}) span from 7.1 to 11.1, with a mean of log M {sub *}/M {sub ⊙} = 8.5, which is lower than that for normal galaxies. The stellar mass-to-light ratio (γ*) in each band varies little with the absolute magnitude but increases with higher M {sub *}. This trend of γ* with M {sub *} is even stronger in bluer bands. In addition, the γ* for our LSBGs slightly declines from the r band to the longer-wavelength bands. The log γ{sub ∗}{sup j} (j = g, r, i, and z) have relatively tight relations with optical colors of g − r and g − i. Compared with several representative γ*–color relations (MLCRs) from the literature, our MLCRs based on LSBG data are consistently among those literature MLCRs previously defined on diverse galaxy samples, and the existing minor differences between the MLCRs are caused by the differences in the SED model ingredients (including initial mass function, star formation history, and stellar population model), line fitting techniques, galaxy samples, and photometric zero-point, rather than the galaxy surface brightness itself, which distinguishes LSBGs from high surface brightness galaxies. Our LSBGs would be very likely to follow those representative MLCRs previously defined in diverse galaxy populations, if those main ingredients were taken into account.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

The high level trigger and express data production at STAR

To meet the demands of the Beam Energy Scan phase-II (BES-II) program, the STAR experiment at the Relativistic Heavy Ion Collider (RHIC) developed a dual real-time framework consisting of a High Level Trigger (HLT) and an Express Data Production system (xProduction). The HLT operates online within the Data Acquisition (DAQ) chain on a dedicated multi-core CPU cluster with the option to offload compute-intensive kernels to Xeon Phi coprocessors. It uses parallelized algorithms, such as the Cellular Automaton (CA) Track Finder, to perform rapid tracking, vertexing, and event filtering. This allows it to select events of interest in real time and provide immediate feedback on detector and beam conditions. In contrast, the xProduction workflow runs concurrently and independently of the DAQ loop. It applies near offline-quality calibration and reconstruction within hours of data collection. The xProduction input is the express data stream, whose content can be enriched by HLT trigger/priority selections under DAQ/HLT resource constraints, and it uses the STAR calibration/conditions framework, incorporating online calibration/QA information when available. This enables early preliminary physics analysis, including the reconstruction of rare signals, such as hyperons and hypernuclei. It also provides collaboration-wide access to analysis-ready datasets. Together, the HLT and xProduction systems form a complementary architecture: the HLT performs online event selection while the xProduction chain delivers high-quality results within a short amount of time. This integrated framework has enabled the prompt reconstruction of the $^5_Λ$ He hypernucleus with high statistical significance and the efficient processing of hundreds of millions of heavy-ion collision events. In conclusion, its demonstrated scalability and robustness establish a model for future high-luminosity experiments requiring both online event filtering and rapid access to analysis-quality data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Star Block Copolymers at Homopolymer Interfaces: Conformation and Compatibilization

Star block copolymers (s-BCPs), comprised of multiple linear diblock copolymers joined at a central point, are shown to segregate to the interface between two immiscible homopolymers that are identical to the blocks of the s-BCPs. The s-BCPs undergo a configurational transition at the interface, with different blocks of copolymers being embedded in their respective homopolymers, thereby bridging the interface and promoting adhesion. A series of 4-arm s-BCPs were synthesized with hydrogenated or deuterated polystyrene (PS/dPS) as the core block and poly(2-vinylpyridine) (P2VP) as the corona block, which was directly placed at the interface between the two homopolymers. Neutron reflectivity (NR) was used to determine the concentration profiles of the PS homopolymer, s-BCP core blocks, and P2VP total segments under equilibrium. The investigation varies the molecular weight (MW) and the total number of s-BCPs at the interface. Self-consistent-field theory (SCFT) was also employed to calculate the concentration profiles of the components at the interface, which were in excellent agreement with experimental results. The NR showed that the interfacial width between the homopolymers increased with the increasing number of s-BCPs at the interface up to a saturation limit. Beyond this limit, additional s-BCPs were released into the corona-miscible phase as unimolecular micelles. For a comparable interlayer thickness of s-BCPs at the interface, lower MW s-BCPs generated a broader interface. SCFT analysis suggested that, at the same packing density, the arms of the low MW s-BCPs align more parallel to the interface, while the arms of high MW s-BCPs adopt a more normal orientation, like their linear BCP counterparts. Furthermore, it was also observed that the core blocks, constrained by the junction points, were oriented more parallel and closer to the interface than the corona blocks. The phase behavior of the polymer blends revealed that s-BCP additives can efficiently reduce the domain size, with the low MW yielding smaller domain sizes due to the greater reduction in the interfacial energy and the high MW arresting phase separation due to their higher binding energy and a jamming of the interfacial assemblies. Asymmetric double cantilever beam (ADCB) tests demonstrated that s-BCPs promoted adhesion more efficiently than their linear BCP counterparts due to stronger binding energy per molecule, suggesting a more efficient compatibilizer for polymer upcycling. In conclusion, the results from these studies provide fundamental insights into the assembly of s-BCPs at homopolymer interfaces, the reduction of domain size, and promotion of adhesion, providing a strategy for the use of s-BCPs as stealth surfactants and universal compatibilizers.

interface↗
DOE OSTI2024

A strainmeter array as the fulcrum of novel observatory sites along the Alto Tiberina Near Fault Observatory

Fault slip is a complex natural phenomenon involving multiple spatiotemporal scales from seconds to days to weeks. To understand the physical and chemical processes responsible for the full fault slip spectrum, a multidisciplinary approach is highly recommended. The Near Fault Observatories (NFOs) aim at providing high-precision and spatiotemporally dense multidisciplinary near-fault data, enabling the generation of new original observations and innovative scientific products. The Alto Tiberina Near Fault Observatory is a permanent monitoring infrastructure established around the Alto Tiberina fault (ATF), a 60 km long low-angle normal fault (mean dip 20°), located along a sector of the Northern Apennines (central Italy) undergoing an extension at a rate of about 3 mm yr –1 . The presence of repeating earthquakes on the ATF and a steep gradient in crustal velocities measured across the ATF by GNSS stations suggest large and deep (5–12 km) portions of the ATF undergoing aseismic creep. Both laboratory and theoretical studies indicate that any given patch of a fault can creep, nucleate slow earthquakes, and host large earthquakes, as also documented in nature for certain ruptures (e.g., Iquique in 2014, Tōhoku in 2011, and Parkfield in 2004). Nonetheless, how a fault patch switches from one mode of slip to another, as well as the interaction between creep, slow slip, and regular earthquakes, is still poorly documented by near-field observation. With the strainmeter array along the Alto Tiberina fault system (STAR) project, we build a series of six geophysical observatory sites consisting of 80–160 m deep vertical boreholes instrumented with strainmeters and seismometers as well as meteorological and GNSS antennas and additional seismometers at the surface. By covering the portions of the ATF that exhibits repeated earthquakes at shallow depth (above 4 km) with these new observatory sites, we aim to collect unique open-access data to answer fundamental questions about the relationship between creep, slow slip, dynamic earthquake rupture, and tectonic faulting.

58 GEOSCIENCES↗
DOE OSTI2024

Numerical simulation of involute-plate research reactor flow behavior using RANS, LES and DNS

This paper investigates the flow behavior of involute-plate research reactors by performing Reynolds-Averaged Navier Stokes simulation (RANS), Large Eddy Simulation (LES) and Direct Numerical Simulation (DNS) of the channel flow between fuel plates. By modeling turbulence with different numerical approaches, this study provides data with three levels of fidelity. For the RANS simulation, three widely used turbulence models, i.e., k-ε, k-ω, Reynolds Stress Turbulence model (RST) are applied by using the commercial CFD code STAR-CCM +. For LES and DNS, the open-source CFD code, Nek5000, is used given its outstanding scalability on High Performance Computer (HPC) and high-order technique. The results from RANS simulations are compared with that from LES and DNS for benchmarking. Both macroscale parameters and turbulence statistics, such as velocity magnitude, lateral velocity and turbulence kinetic energy, are presented and analyzed. The results from RANS simulation achieve good agreement with LES and DNS on velocity and turbulence kinetic energy prediction. The RST turbulence model predicts the most similar flow pattern of lateral velocity as compared to LES and DNS. The Lambda-2 (λ2) criterion with a reasonable threshold is used to demonstrate the instantaneous vortices distribution in the involute channel from both LES and DNS calculation. The DNS simulation captures more detailed turbulence especially near the corner, which explains the discrepancy between LES and DNS results near the corner. The normalized RMS error are defined and calculated to assess the performance of those turbulence models. The RST model captures the anisotropic feature of turbulence, which enable it to outperform other turbulence models for predicting the flow behavior in an involute channel. Although some discrepancies are found between LES and DNS results in the corner, the overall deviations between LES and DNS are found to be small. In conclusion, given that the computational cost of DNS calculation is an order of magnitude higher, using LES data for benchmarking RANS model is a cost-effective approach.

DNS↗
DOE OSTI2020

Compact galaxies and the size–mass galaxy distribution from a colour-selected sample at 0.04 < z < 0.15 supplemented by ugrizYJHK photometric redshifts

ABSTRACT The size–mass galaxy distribution is a key diagnostic for galaxy evolution. Massive compact galaxies are potential surviving relics of a high-redshift phase of star formation. Some of these could be nearly unresolved in Sloan Digital Sky Survey (SDSS) imaging and thus not included in galaxy samples. To overcome this, a sample was selected from the combination of SDSS and UKIRT Infrared Deep Sky Survey (UKIDSS) photometry to r < 17.8. This was done using colour–colour selection, and then by obtaining accurate photometric redshifts (photo-z) using scaled flux matching (SFM). Compared to spectroscopic redshifts (spec-z), SFM obtained a 1σ scatter of 0.0125 with only 0.3 per cent outliers (|Δln (1 + z)| > 0.06). A sample of 163 186 galaxies was obtained with 0.04 < z < 0.15 over $2300\, {\rm deg}^2$ using a combination of spec-z and photo-z. Following Barro et al. log Σ1.5 = log M* − 1.5log r50, maj was used to define compactness. The spectroscopic completeness was 76 per cent for compact galaxies (log Σ1.5 > 10.5) compared to 92 per cent for normal-sized galaxies. This difference is primarily attributed to SDSS ‘fibre collisions’ and not the completeness of the main galaxy sample selection. Using environmental overdensities, this confirms that compact quiescent galaxies are significantly more likely to be found in high-density environments compared to normal-sized galaxies. By comparison with a high-redshift sample from 3D-HST, log Σ1.5 distribution functions show significant evolution, with this being a compelling way to compare with simulations such as EAGLE. The number density of compact quiescent galaxies drops by a factor of about 30 from z ∼ 2 to log (n/Mpc−3) = − 5.3 ± 0.4 in the SDSS–UKIDSS sample. The uncertainty is dominated by the steep cut off in log Σ1.5, which is demonstrated conclusively using this complete sample.

Baldry, Ivan K.↗
DOE OSTI2025

Leptodermic corrections to the TOV equations and nuclear astrophysics within the effective surface approximation

The macroscopic model for a neutron star (NS) as a liquid drop at the equilibrium is used to extend the Tolman-Oppenheimer-Volkoff (TOV) equations taking into account the gradient terms responsible for the system surface. The parameters of the Schwarzschild metric in the spherical case are found with these surface corrections to the known leading (zero) order of the leptodermic approximation a/R << 1, where a is the NS effective-surface (ES) thickness, and R is the effective NS radius. The energy density $\mathscr{E}$ is considered in a general form including the functions of the particle number density and of its gradient terms. The macroscopic gravitational component $Φ$(ρ) of the energy density is taken into account in the simplest form as expansion in powers of $ρ$ – $\overline{ρ}$, where $\overline{ρ}$ is the saturation density, up to second order, in terms of its contributions to the separation particle energy and incompressibility. Density distributions ρ across the NS ES in the normal direction to the ES, which are derived in the simple analytical form at the same leading approximation, was used for the derivation of the modified TOV (MTOV) equations by accounting for their NS surface corrections. As a result, the MTOV equations are analytically solved at first order and the results are compared with the standard TOV approach of the zero order.

Magner, A. G. [Institute for Nuclear Research, Kyi↗
DOE OSTI2020

Blind Modeling Validation Exercises Using the Horizontal Dry Cask Simulator

The U.S. Department of Energy (DOE) established a need to understand the thermal-hydraulic properties of dry storage systems for commercial spent nuclear fuel (SNF) in response to a shift towards the storage of high-burnup (HBU) fuel (> 45 gigawatt days per metric ton of uranium, or GWd/MTU). This shift raises concerns regarding cladding integrity, which faces increased risk at the higher temperatures within spent fuel assemblies present within HBU fuel compared to low-burnup fuel (≤ 45 GWd/MTU). A dry cask simulator (DCS) was built at Sandia National Laboratories (SNL) in Albuquerque, New Mexico to produce validation-quality data that can be used to test the accuracy of the modeling used to predict cladding temperatures. These temperatures are critical to evaluating cladding integrity throughout the storage cycle of commercial spent nuclear fuel. A model validation exercise was previously carried out for the DCS in a vertical configuration. Lessons learned during the previous validation exercise have been applied to a new, blind study using a horizontal dry cask simulator (HDCS). Three modeling institutions – the Nuclear Regulatory Commission (NRC), Pacific Northwest National Laboratory (PNNL), and Empresa Nacional del Uranio, S.A., S.M.E. (ENUSA) – were granted access to the input parameters from the DCS Handbook, SAND2017-13058R, and results from a limited data set from the horizontal BWR dry cask simulator tests reported in the HDCS update report, SAND2019-11688R. With this information, each institution was tasked to calculate peak cladding temperatures and air mass flow rates for ten HDCS test cases. Axial as well as vertical and horizontal transverse temperature profiles were also calculated. These calculations were done using modeling codes (ANSYS/Fluent, STAR-CCM+, or COBRA-SFS), each with their own unique combination of modeling assumptions and boundary conditions. For this validation study, the ten test cases of the horizontal dry cask simulator were defined by three independent variables – fuel assembly decay heat (0.5 kW, 1 kW, 2.5 W, and 5 kW), internal backfill pressure (100 kPa and 800 kPa), and backfill gas (helium and air). The plots provided in Chapter 3 of this report show the axial, vertical, and horizontal temperature profiles obtained from the dry cask simulator experiments in the horizontal configuration and the corresponding models used to describe the thermal-hydraulic behavior of this system. The tables provided in Chapter 3 illustrate the closeness of fit of the model data to the experiment data through root mean square (RMS) calculations of the error in peak cladding temperatures (PCTs), PCT axial locations, axial temperature profiles, vertical and horizontal temperature profiles at two different axial locations, and air mass flow rates for the ten test cases, normalized by the experimental results. The model results are assigned arbitrary model numbers to retain anonymity. Due to the relatively flat axial temperature profiles, small temperature gradients resulted in large deviations of all models’ PCT axial location from the experimental PCT axial location. When the PCT axial location error is excluded in the calculation of the combined RMS of the normalized errors that considers PCT, the temperature profiles, and the air mass flow rates, the model data fits the experimental data to within 5%. When the vault information is excluded, the model data fits the experimental data to within 2.5%. An error analysis was developed further for one model, using the model and experimental uncertainties in each validation parameter to calculate validation uncertainties. The uncertainties for each parameter were used to define quantifiable validation criteria. For this analysis, the model was considered validated for a given comparison metric if the normalized error in that metric divided by the validation uncertainty was less than or equal to 1. When considering the combined RMS of the normalized errors of all metrics divided by their validation uncertainties, the model was found to have satisfied the criterion for model validation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Quasars with Proper Motions and the Link to Double and Multiple AGNs

Gaia used a large sample of photometrically selected active galactic nuclei (AGNs) and quasars to remove the residual spin of its global proper motion system in order to achieve a maximally inertial reference frame. A small fraction of these reference objects have statistically significant astrometric proper motions in Gaia EDR3. We compile a source sample of 105,593 high-fidelity AGNs with accurate spectroscopically determined redshifts above 0.5 from the SDSS and normalized proper motions below 4. The rate of genuinely perturbed proper motions is at least 0.17%. A smaller high completeness sample of 152 quasars with excess proper motions at a confidence level of 0.9995 is examined in detail. Pan-STARRS images and Gaia-resolved pairs reveal that 29% of the sample are either double sources or gravitationally lensed quasars. An Anderson–Darling test on parameters of a smaller high-reliability sample and their statistical controls reveals 17 significant factors that favor multiplicity and multi-source structure as the main cause of perturbed astrometry. Using a nearest-neighbor distance statistical analysis and counts of close companions in Gaia on a much larger initial sample of AGNs, an excess of closely separated sources in Gaia is detected. At least 0.33% of all optical quasars are genuinely double or multiply imaged. We provide a list of 44 candidate double or multiple AGNs and four previously known gravitational lenses. Many proper motion quasars may be more closely separated, unresolved doubles exhibiting the variability imposed motion effect, and a smaller fraction may be chance alignments with foreground stars causing weak gravitational lensing.

Astrometric microlensing effect↗
DOE OSTI2022

Preliminary Thermal-Hydraulic Analysis of the UTA–2 Subcritical Linear Accelerator–Driven System

The National Nuclear Security Administration’s mission includes establishing a reliable supply of 99 Mo without highly enriched uranium. Oak Ridge National Laboratory (ORNL) supports this objective through collaborative research and development with industrial partners. Niowave Inc., a current partner, is currently designing a subcritical linear accelerator-driven system (ADS) and preparing for the US Nuclear Regulatory Commission’s licensing process. Niowave’s system has the potential to efficiently supply medical radioisotopes. The technology includes a superconducting electron accelerator and a pile of both natural and low-enriched uranium (LEU) targets. The process fissions uranium and many valuable isotopes can then be extracted from the targets. Niowave is currently iterating through conceptual and detailed design processes for several system sizes. This report discusses UTA–2, which is at the demonstration stage. UTA–2 will validate numerical modeling results with experimental measurements before progressing to the detailed design of UTA–3, the commercial-sized ADS. The thermal-hydraulic behavior of the UTA–2 core design was numerically investigated using STAR-CCM+, as described in this report. STAR-CCM+, a state-of-the-art computational fluid dynamics (CFD) software that was commercially developed by Siemens, has an extensive user base and a set of validation studies. It is also compliant with the American Society of Mechanical Engineers’ Nuclear Quality Assurance 1 standard. Three cases are investigated in this report. Case 1 quantified the temperature field within the UTA–2 assembly and water tank using only conduction as the method of thermal energy transport. This simplified approach was overly conservative and yielded wetted cladding temperatures above the coolant saturation temperature. In this case, the maximum temperature of the wetted cladding surface of the highest power rod exceeded the saturation temperature of water by 63.8°C. Because of the overly conservative approach taken in Case 1 and its negative subcooled margin, buoyancy-driven natural circulation flow physics were implemented in Case 2. Adding coolant motion significantly distributed the thermal energy of the system through convective heat transfer. This relatively small amount of convective heat transfer significantly reduced system temperatures and increased the subcooled margin from –63.8 to 61.3°C. This margin confirmed that no boiling was expected during normal operation of UTA–2 at 230 W. Case 3 had no additional physics models but considered an overpower event during which the power of each LEU and natural uranium rod was at its respective peak values. This resulted in a study with the total assembly power equal to 176% of the nominal power of 230 W considered in Cases 1 and 2. The resulting natural circulation flows were slightly enhanced. The subcooled margin decreased slightly to 46.3°C, which is still a significant margin to local boiling of the water within the tank. This margin confirmed that no boiling was expected during an abnormal operation of UTA–2 at 406 W.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2020

3D turbulent reconnection: Theory, tests, and astrophysical implications

Magnetic reconnection, topological changes in magnetic fields, is a fundamental process in magnetized plasmas. It is associated with energy release in regions of magnetic field annihilation, but this is only one facet of this process. Astrophysical fluid flows normally have very large Reynolds numbers and are expected to be turbulent, in agreement with observations. In strong turbulence, magnetic field lines constantly reconnect everywhere and on all scales, thus making magnetic reconnection an intrinsic part of the turbulent cascade. We note in particular that this is inconsistent with the usual practice of magnetic field lines as persistent dynamical elements. A number of theoretical, numerical, and observational studies starting with the paper done by Lazarian and Vishniac [Astrophys. J. 517, 700–718 (1999)] proposed that 3D turbulence makes magnetic reconnection fast and that magnetic reconnection and turbulence are intrinsically connected. In particular, we discuss the dramatic violation of the textbook concept of magnetic flux-freezing in the presence of turbulence. We demonstrate that in the presence of turbulence, the plasma effects are subdominant to turbulence as far as the magnetic reconnection is concerned. The latter fact justifies a magnetohydrodynamiclike treatment of magnetic reconnection on all scales much larger than the relevant plasma scales. We discuss the numerical and observational evidence supporting the turbulent reconnection model. In particular, we demonstrate that the tearing reconnection is suppressed in 3D, and unlike the 2D settings, 3D reconnection induces turbulence that makes magnetic reconnection independent of resistivity. We show that turbulent reconnection dramatically affects key astrophysical processes, e.g., star formation, turbulent dynamo, and acceleration of cosmic rays. We provide criticism of the concept of “reconnection-mediated turbulence” and explain why turbulent reconnection is very different from enhanced turbulent resistivity and hyper-resistivity and why the latter have fatal conceptual flaws.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

The morphology of the X-ray afterglows and of the jetted GeV emission in long GRBs

ABSTRACT We recall evidence that long gamma-ray bursts (GRBs) have binary progenitors and give new examples. Binary-driven hypernovae (BdHNe) consist of a carbon–oxygen core (COcore) and a neutron star (NS) companion. For binary periods ∼5 min, the COcore collapse originates the subclass BdHN I characterized by (1) an outstanding supernova (SN; the ‘SN-rise’); (2) a black hole (BH), born from the NS collapse by SN matter accretion, leading to a GeV emission with luminosity $L_{\rm GeV} = A_{\rm GeV}\, t^{-\alpha _{\rm GeV}}$, observed only in some cases; and (3) a new NS (νNS), born from the SN, originating from the X-ray afterglow with $L_\mathrm{ X} = A_{\rm X}\, t^{-\alpha _{\rm X}}$, observed in all BdHN I. We record 378 sources and present for four prototype GRBs 130427A, 160509A, 180720B, and 190114C: (1) spectra, luminosities, SN-rise duration; (2) AX, αX = 1.48 ± 0.32, and (3) the νNS spin time evolution. We infer (i) AGeV, αGeV = 1.19 ± 0.04 and (ii) the BdHN I morphology from time-resolved spectral analysis, three-dimensional simulations, and the GeV emission presence/absence in 54 sources within the Fermi-Large Area Telescope boresight angle. For 25 sources, we give the integrated and time-varying GeV emission, 29 sources have no GeV emission detected and show X/gamma-ray flares previously inferred as observed along the binary plane. The 25/54 ratio implies the GeV radiation is emitted within a cone of half-opening angle ≈60° from the normal to the orbital plane. We deduce BH masses of 2.3–8.9 M⊙ and spin of 0.27–0.87 by explaining the GeV emission from the BH rotational energy extraction, while their time evolution validates the BH mass–energy formula.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Measurements of the z > 5 Lyman-α forest flux autocorrelation functions from the extended XQR-30 data set

We present the first observational measurements of the Lyman-α (Ly α) forest flux autocorrelation functions in ten redshift bins from 5.1 ≤ z ≤ 6.0. We use a sample of 35 quasar sightlines at z > 5.7 from the extended XQR-30 data set; these data have signal-to-noise ratios of >20 per spectral pixel. We carefully account for systematic errors in continuum reconstruction, instrumentation, and contamination by damped Ly α systems. With these measurements, we introduce software tools to generate autocorrelation function measurements from any simulation. Our measurements of the smallest bin of the autocorrelation function increase with redshift when normalizing by the mean flux, $\langle{F}\rangle$. This increase may come from decreasing $\langle{F}\rangle$ or increasing mean free path of hydrogen-ionizing photons, λmfp. Recent work has shown that the autocorrelation function from simulations at z > 5 is sensitive to λmfp, a quantity that contains vital information on the ending of reionization. For an initial comparison, we show our autocorrelation measurements with simulation models for recently measured λmfp values and find good agreements. Further work in modelling and understanding the covariance matrices of the data is necessary to get robust measurements of λmfp from this data.

79 ASTRONOMY AND ASTROPHYSICS↗