Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Cross points”

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 631 records · Page 35

Salt Vessel-Sample Generator Interface

The goal of the salt sampling program at Argonne is to develop and deploy automated molten salt sampling approaches for interfacing relevant unit operations with salt analysis to improve the timeliness and accuracy of sampling-based accountancy measurements. One technology under development is a vacuum sampler loop module. In this system, molten salt is drawn from a process vessel through a heated sampling line and into a charge vessel. Next, one or more samples are split from the volume in the charge vessel and the remaining salt is returned to the process. The vacuum sampler loop module is being developed to interface with other sample processing and analysis modules to enable rapid at-line sample characterization. The main purpose of the vacuum sampling loop module is to enable the collection of samples that are more representative of the bulk salt by replacing traditional point samplers (i.e., dip probes) with a sampling approach that captures a larger cross section of salt. Additionally, the vacuum sampling approach eliminates the risk of dross contamination of samples and avoids the use of moving parts in the salt. Two methods of interfacing the vacuum sampling loop with a precision pneumatic sample generator were investigated in FY21. This report covers the testing of fluidic coupling between the two modules. Two iterations of the fluidically coupled modules were tested. The first iteration system coupled the two modules using a freeze valve to seal the vacuum sampler during filling and to control flow into the pneumatic sample generator. While this integrated system functioned as intended, some changes were implemented to make the system more robust and better suited to remote deployment. Specifically, the system was made to be more modular and active control of the vacuum filling operation was replaced with a passive control mechanism. For passive filling, the salt charge vessel was vented to a small gas tank that was at negative pressure, causing salt to be drawn into the vessel until the force of the fluid head was in equilibrium with the gas pressure. The passive control system performed well and will be used in future iterations. Another change in the second system was a newly configured pneumatic sample generator in which sample ejection occurred through a hole in the reservoir’s stainless-steel side wall instead of through a non-wetted sapphire orifice on the bottom. This alternate configuration may be better suited for near-process deployment because it enables on-line orifice maintenance and an orifice bypass drain back to the process. A third change in the second iteration system was the transition to a two-chamber charge vessel which split off a fraction of the sampled salt as a liquid aliquot. The goal was to create a buffer mechanism that would allow reproducible aliquoting of samples, independently of variability in the charge vessel fill height. While the two-chamber design was functional, the sample size reproducibility was below target values. To improve reproducibility and overcome many of the impediments to remote deployment of the vacuum sampling loop module, separate work was conducted to replace the two-chamber liquid aliquoting mechanism with aliquoting into single-use sample tubes. Solid salt transfer in the sample tubes will replace fluidic coupling for integrating the vacuum sampler with downstream modules. Because the proposed operations can all be executed with simple overhead actuation mechanisms or other existing hot-cell technology, there will be no need for large investments in novel hot cell sample handling technologies using this alternate approach. As such, near-term deployment of the vacuum sampling technology will be achievable. This new approach for automated coupling of sample tubes with down-stream modules is covered in a separate FY21 report, and a remotely operated version of the vacuum sampler loop module with the tube aliquoting feature is planned for FY22.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Optical measurements of the droplet size distribution in the case of fuel atomization in swirl nozzles and planar airblast diffusers

The theoretical principles of optical methods for the determination of the particle sizes of sprays are considered and aspects of the experimental implementation of these principles are discussed. An experimental device for point-intensity measurements makes use of a helium-neon laser. The cross-sectional area of the laser beam is enlarged with the aid of a lens system to the size of the measurement cross-section. The intensity of the laser light scattered by the spray particles is measured as a function of light direction. Approaches which take into account the total energy of the diffractively scattered light are also discussed and an investigation is conducted regarding the measurement error sources. A description is presented of experimental results obtained in studies of a number of fuel nozzle sprays.

Kayser, A.↗

Iterative solution of the multistream electron transport equation. I - Comparison with laboratory beam injection experiments

The Neumann iteration method presently used for solving the electron transport equation in which energy, attitude, and pitch angle are independent variables is fast, and can compute numerical point-response-function solutions of the electron transport equation. Because both the inelastic cross sections and angular elastic cross sections of the model are empirically based, the solutions obtained represent a test of compatibility between various sets of cross sections and energy deposition measurements. The use of a numerical quadrature based on analytic phase function forms yields accurate phase function integrals at low computational cost.

Porter, H. S.↗

Setting up a Rayleigh Scattering Based Flow Measuring System in a Large Nozzle Testing Facility

A molecular Rayleigh scattering based air density measurement system has been built in a large nozzle testing facility at NASA Glenn Research Center. The technique depends on the light scattering by gas molecules present in air; no artificial seeding is required. Light from a single mode, continuous wave laser was transmitted to the nozzle facility by optical fiber, and light scattered by gas molecules, at various points along the laser beam, is collected and measured by photon-counting electronics. By placing the laser beam and collection optics on synchronized traversing units, the point measurement technique is made effective for surveying density variation over a cross-section of the nozzle plume. Various difficulties associated with dust particles, stray light, high noise level and vibration are discussed. Finally, a limited amount of data from an underexpanded jet are presented and compared with expected variations to validate the technique.

Panda, Jayanta↗

Androgynous { 10 1 ¯ 2 } twin in zinc

Under ambient conditions, Zn is a hexagonal metal with a large c / a ratio of 1.856. Plastic deformation is predominantly accommodated by basal 〈 a 〉 slip and compression twins on the { 10 1 ¯ 2 } planes. Increasing hydrostatic pressure drastically reduces the c / a ratio of Zn and, when a critical threshold of c / a = 3 at about 10 GPa is crossed, the { 10 1 ¯ 2 } twins are predicted to change from compressive to tensile in nature. What happens at the transition point, when c / a = 3 , remains unknown. Here, we strain-cycle a textured polycrystalline sample of pure Zn at uniform hydrostatic pressures ranging between 2 and 17 GPa, over which the c / a ratio crosses the c / a = 3 compressive-tensile transition for { 10 1 ¯ 2 } twins. During deformation, the state of the sample is monitored through x-ray diffraction to extract texture and internal strain evolution. By comparing the experimental results with the predictions of an elastoviscoplastic polycrystal simulation, we confirm the androgynous nature of { 10 1 ¯ 2 } twin response at low and high pressures. When c / a = 3 , polycrystalline Zn does not display any evidence of twinning and its plastic behavior is controlled by mostly basal and pyramidal 〈 c + a 〉 slip activity, with a very small contribution of prismatic 〈 a 〉 slip. Evidence for the activity of other { 10 1 ¯ n } twinning modes, which have been suggested for Zn under high pressure, are not observed. Published by the American Physical Society 2024

Merkel, Sébastien (ORCID:000000032767581X)↗

Cosmological analysis of three-dimensional BOSS galaxy clustering and Planck CMB lensing cross correlations via Lagrangian perturbation theory

In this work, we present a formalism for jointly fitting pre- and post-reconstruction redshift-space clustering (RSD) and baryon acoustic oscillations (BAO) plus gravitational lensing (of the CMB) that works directly with the observed 2-point statistics. The formalism is based upon (effective) Lagrangian perturbation theory and a Lagrangian bias expansion, which models RSD, BAO and galaxy-lensing cross correlations within a consistent dynamical framework. As an example we present an analysis of clustering measured by the Baryon Oscillation Spectroscopic Survey in combination with CMB lensing measured by Planck. The post-reconstruction BAO strongly constrains the distance-redshift relation, the full-shape redshift-space clustering constrains the matter density and growth rate, and CMB lensing constrains the clustering amplitude. Using only the redshift space data we obtain Ω m = 0.303 ± 0.008, $H_0$ = 69.21 ± 0.78 and $σ_8$ = 0.743 ± 0.043. The addition of lensing information, even when restricted to the Northern Galactic Cap, improves constraints to Ω m = 0.303 ± 0.008, $H_0$ = 69.21 ± 0.77 and $σ_8$ = 0.707 ± 0.035, in tension with CMB and cosmic shear constraints. The combination of Ω m and $H_0$ are consistent with Planck, though their constraints derive mostly from redshift-space clustering. The low $σ_8$ value are driven by cross correlations with CMB lensing in the low redshift bin ($\textit{z}$ ≃ 0.38) and at large angular scales, which show a 20% deficit compared to expectations from galaxy clustering alone. We conduct several systematics tests on the data and find none that could fully explain these tensions.

79 ASTRONOMY AND ASTROPHYSICS↗

Progress in DIII-D towards validating divertor power exhaust predictions

In this study, UEDGE simulations highlight the role of cross-field drifts on the onset of detached conditions, and new calibrated divertor vacuum ultra violet (VUV) spectroscopy is used to challenge the predictions of radiative constituents in these simulations. UEDGE simulations for DIII-D H-mode plasmas with the open divertor with the ion ∇B-drift towards the X-point show a bifurcated onset of the low field side (LFS) divertor detachment, consistent with experimentally observed step-like detachment onset (Jaervinen A.E. et al 2018 Phys. Rev. Lett. 121 075001). The divertor plasma in the simulations exhibits hysteresis in upstream separatrix density between attached and detached solution branches. Reducing the drift magnitude by a factor of 3 eliminates the step-like detachment onset in the simulations, confirming the strong role of drifts in the bifurcated detachment onset. When measured local plasma densities and temperatures are within proximity of predicted values in the simulations, there is no shortfall of the local emission of the dominant resonant radiating lines. However, the simulations systematically predict a factor of two lower total integrated radiated power than measured by the bolometer with the difference lost through radial heat flow out of the computational domain. Even though there is no shortfall in the emission of the dominant lines, a shortfall of total radiated power can be caused by underpredicted spatial extent of the radiation front, indicating a potential upstream or divertor transport physics origin for the radiation shortfall, or shortfall of radiated power in the spectrum between the dominant lines. In addition to the underpredicted spatial extent, in detached conditions, the simulations overpredict the peak radiation and dominant carbon lines near the X-point, which can be alleviated by manually increasing divertor diffusivity in the simulations, highlighting the ad hoc cross-field transport as one of the key limitations of the predictive capability of these divertor fluid codes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Calculation of flow about posts and powerhead model

A three dimensional analysis of the non-uniform flow around the liquid oxygen (LOX) posts in the Space Shuttle Main Engine (SSME) powerhead was performed to determine possible factors contributing to the failure of the posts. Also performed was three dimensional numerical fluid flow analysis of the high pressure fuel turbopump (HPFTP) exhaust system, consisting of the turnaround duct (TAD), two-duct hot gas manifold (HGM), and the Version B transfer ducts. The analysis was conducted in the following manner: (1) modeling the flow around a single and small clusters (2 to 10) of posts; (2) modeling the velocity field in the cross plane; and (3) modeling the entire flow region with a three dimensional network type model. Shear stress functions which will permit viscous analysis without requiring excessive numbers of computational grid points were developed. These wall functions, laminar and turbulent, have been compared to standard Blasius solutions and are directly applicable to the cylinder in cross flow class of problems to which the LOX post problem belongs.

Anderson, P. G.↗

Selecting Charged Current Muon Neutrino Interactions on Argon with the Short-Baseline Near Detector

The uncertainties associated with neutrino-nucleus cross section modelling need to be reduced in order to measure neutrino oscillation parameters, the mass hierarchy and CP violation to a high degree of precision. SBND, a 112 ton liquid argon time projection chamber, will be the near detector of the short-baseline neutrino program at Fermilab. Once data taking begins in 2022, it will provide flux constraints for sterile neutrino searches and produce world leading neutrino-argon cross sections with up to seven million neutrino events in three years. This thesis describes the simulation, calibration and reconstruction tools developed to measure charged current muon neutrino interactions.Backgrounds from cosmic ray muons, neutrino interactions outside of the detector and different neutrino interaction channels were identified and removed as part of this process.The signal events were selected with purity of 90% and an efficiency of 67%, providing a strong starting point for meas urements of exclusive channels. The expected statistical and systematic uncertainties on rate predictions and cross section measurements at the full 6.6×10$^{20}$ proton on target exposure were evaluated.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

ICESat Observations of Inland Surface Water Stage, Slope, and Extent: a New Method for Hydrologic Monitoring

River discharge and changes in lake, reservoir and wetland water storage are critical terms in the global surface water balance, yet they are poorly observed globally and the prospects for adequate observations from in-situ networks are poor (Alsdorf et al., 2003). The NASA-sponsored Surface Water Working Group has established a framework for advancing satellite observations of river discharge and water storage changes which focuses on obtaining measurements of water surface height (stage), slope, and extent. Satellite laser altimetry, which can achieve centimeter-level elevation precision for single, small laser footprints, provides a method to obtain these inland water parameters and contribute to global water balance monitoring. Since its launch in January, 2003 the Ice, Cloud, and land Elevation Satellite (ICESat), a NASA Earth Observing System mission, has achieved over 540 million laser pulse observations of ice sheet, ocean surface, land topography, and inland water elevations and cloud and aerosol height distributions. By recording the laser backscatter from 80 m diameter footprints spaced 175 m along track, ICESat acquires globally-distributed elevation profiles, using a 1064 nm laser altimeter channel, and cloud and aerosol profiles, using a 532 nm atmospheric lidar channel. The ICESat mission has demonstrated the following laser altimeter capabilities relevant to observations of inland water: (1) elevation measurements with a precision of 2 to 3 cm for flat surfaces, suitable for detecting river surface slopes along long river reaches or between multiple crossings of a meandering river channel, (2) from the laser backscatter waveform, detection of water surface elevations beneath vegetation canopies, suitable for measuring water stage in flooded forests, (3) single pulse absolute elevation accuracy of about 50 cm (1 sigma) for 1 degree sloped surfaces, with calibration work in progress indicating that a final accuracy of about 12 cm (1 sigma) will be achieved for clear atmosphere conditions, suitable for detection of stage changes through time, (4) ability to precisely point the spacecraft so as to position the laser profile on the Earth the surface with a cross-track accuracy of 50 m (1 sigma), enabling small water bodies and specific locations to be targeted and re-observed through time, (5) adequate signal levels from specular water surfaces up to 5 degrees off-nadir, enabling complete global access to any location on the Earth's surface from the ICESat repeat orbit by off-nadir pointing, and (6) day and night operation with successful laser ranging to the Earth's surface through thin to moderate cloud cover, enabling more frequent measurements than can be achieved by passive optical sensors. Here we illustrate these capabilities by showing ICESat observations through time for selected river and lake locations.

Harding, David J.↗

Correction of downwash in wind tunnels of circular and elliptic sections

The downwash velocity distribution behind the wing was determined for the free jet and for the closed tunnel of both circular and elliptic cross sections. The wing was placed at the center of the tunnel. The theory makes it possible to determine the downwash at any point in the jet. The computations were performed for points in the plane determined by the jet axis and the center-of-pressure line of the wing. The downwash proved to be proportional to the wing lift and inversely proportional to the cross-sectional area of the tunnel.

Lotz, Irmgard↗

Ab Initio Evaluation of Plutonium Dioxide $S(α,β)$ and Thermal Neutron Cross Sections

Plutonium Dioxide (PuO 2 ) is an advanced fuel-cycle nuclear fuel material for thermal and fast reactors, with a melting point at 3261.8 K and thermal conductivity above 3 Wm -1 K -1 for temperatures below 1500 K. Plutonium from spent UO 2 is reprocessed to form chemically stable mixed oxide (MOX) fuel. In contrast to metallic fuels, the low PuO 2 conductivity results in high thermal gradients between fuel center and surface, resulting in structural annealing and efficient fission product migration over burnup. Crystal binding affects the doppler broadening of epithermal resonances for nuclear fuels and scattering behavior of thermal neutrons. Currently, there is no ENDF/B Thermal Scattering Law, i.e., $S(α,β)$, evaluation for PuO 2 , only free atom cross sections for Plutonium and Oxygen, respectively, which neglect lattice contributions to total cross section. Herein, ab initio lattice dynamics (AILD) techniques are employed to calculate the phonon density of states (DOS) using spin-orbit-coupling density functional theory (DFT) to predict the PuO2 paramagnetic ground-state structure. These results highlight the evaluation of $S(α,β)$ for Plutonium and Oxygen, respectively in PuO 2 and consequential generation of thermal neutron scattering cross sections for high fidelity criticality safety analysis and reactor calculations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Theoretical studies of the radar properties of the icy Galilean moons of Jupiter

The icy Galilean satellites of Jupiter - Europa, Ganymede, and Callisto - have unusual radar scattering properties compared with those of the terrestrial planets or Earth's Moon. There are three main features of the data that distinguish these targets: (1) the radar cross-section normalized by the geometrical cross-section is an order of magnitude larger than that of any terrestrial planet; (2) the reflected power is almost evenly distributed between two orthogonal polarizations with more power being returned in the same circular polarization as was transmitted whereas virtually all of the power returned from the terrestrial planets is contained in the opposite circular polarization to the one that was transmitted; and (3) the echo power spectra have a broad shape indicating a nearly uniformly radar-bright surface in contrast to the spectra from the terrestrial planets that contain a strong quasi-specular component from the vicinity of the sub-radar point and very little reflected power from the rest of the surface. The normalized radar cross-sections decrease as the areal water ice coverage decreases from Europa to Ganymede to Callisto. Recently, radar echoes from the polar caps of Mars and Mercury, and from Saturn's satellite Titan imply similarly strong cross-sections and have classically unexpected polarization properties and it is also thought that this is due to the presence of ice on the surface. A model called the radar glory model is analyzed and it is shown that the main features of the radar echoes calculated from this model agree well with the observations from all three icy Galilean satellites. This model involves long radar paths in the ice below the surface and special structures in which the refractive index decreases abruptly at a hemispherical boundary. It is not known whether such structures exist or how they could be created, but possible scenarios can be imagined such as the formation of an impact crater followed by deposition of a frost layer followed by a resurfacing event in which a layer of solid ice is placed above the layer of frost. Regardless of the exact geophysical processes required to create such structures, the superior ability of this model to account for all of the important observations with very few adjustable parameters and with no ad hoc assumptions is a compelling argument in support of at least the electromagnetic model. The key features of the electromagnetic model are multiple subsurface scattering events, total internal reflection, and a low degree of randomness imposed on a deterministic geometry that strongly favors backscattering.

Eshleman, Von R.↗

Asymmetric response of Amazon forest water and energy fluxes to wet and dry hydrological extremes reveals onset of a local drought-induced tipping point

Understanding the effects of intensification of Amazon basin hydrological cycling—manifest as increasingly frequent floods and droughts—on water and energy cycles of tropical forests is essential to meeting the challenge of predicting ecosystem responses to climate change, including forest “tipping points”. Here, we investigated the impacts of hydrological extremes on forest function using 12+ years of observations (between 2001–2020) of water and energy fluxes from eddy covariance, along with associated ecological dynamics from biometry, at the Tapajós National Forest. Measurements encompass the strong 2015–2016 El Niño drought and La Niña 2008–2009 wet events. We found that the forest responded strongly to El Niño-Southern Oscillation (ENSO): Drought reduced water availability for evapotranspiration ( ET ) leading to large increases in sensible heat fluxes ( H ). Partitioning ET by an approach that assumes transpiration ( T ) is proportional to photosynthesis, we found that water stress-induced reductions in canopy conductance ( G s ) drove T declines partly compensated by higher evaporation ( E ). By contrast, the abnormally wet La Niña period gave higher T and lower E , with little change in seasonal ET . Both El Niño-Southern Oscillation (ENSO) events resulted in changes in forest structure, manifested as lower wet-season leaf area index. However, only during El Niño 2015–2016, we observed a breakdown in the strong meteorological control of transpiration fluxes (via energy availability and atmospheric demand) because of slowing vegetation functions (via shutdown of G s and significant leaf shedding). Drought-reduced T and G s , higher H and E , amplified by feedbacks with higher temperatures and vapor pressure deficits, signaled that forest function had crossed a threshold, from which it recovered slowly, with delay, post-drought. In conclusion, identifying such tipping point onsets (beyond which future irreversible processes may occur) at local scale is crucial for predicting basin-scale threshold-crossing changes in forest energy and water cycling, leading to slow-down in forest function, potentially resulting in Amazon forests shifting into alternate degraded states.

54 ENVIRONMENTAL SCIENCES↗

Initial Mobility Analysis for ORNL VA-EDH Synthetic Populations

Travel burdens are a major barrier to healthcare access among US Veteran patient populations, particularly those residing in rural areas. Spatial accessibility to points of care for US Veteran populations is commonly assessed in two ways. The first approach uses open data from the US Census to represent collective travel burdens, for example the distance between population-weighted census tract centroids and VHA points of care. The second approach uses restricted-access VHA patient data to measure travel costs (e.g., distance, time) for accessing points of care with respect to geolocated patient addresses and real or approximated transportation networks. While the advantage of the open data approach lies in its reproducibility, it has notable limitations in its tendency to infer individual travel behavior from aggregate population characteristics, a problem known as ecological fallacy. Conversely, while the patient data approach is able to account for individual travel behavior, its ability to account for localized access disparities (e.g., a neighborhood with exceptionally high transportation costs) and patient demographics is limited as protecting individual patient data requires their storage in closed systems with limited capacity for adequately modeling real-world travel patterns or for supplementing patient attributes. Additionally, the patient data approach cannot account for veterans who are not enrolled in the VHA system but who may be eligible for care. These challenges limit the ability to perform “what if” analyses on the effects of place-specific interventions on veteran populations with high access barriers to healthcare. To address these challenges, we explore the application of realistic synthetic populations to examine travel burdens and spatial accessibility issues among veteran patient populations. Synthetic populations provide a virtual, individually-resolved and cross-sectional representation of the veteran patient population that enables investigation of spatial access to points of care in ways in which aggregate data and patient data do not. First, synthetic populations allow one to directly assess how individuals access points of care, from synthesized residential locations to outpatient facilities on real-world transportation networks. Modeling access to points of care at the individual scale addresses the ecological fallacy problem associated with using aggregated census data to represent veteran populations and patterns of movement. Second, synthetic populations provide a means of completely representing an area’s veteran population using only publicly available, anonymized census microdata from the American Community Survey (ACS) to ensure the privacy of real-world individuals. Generating synthetic populations from the ACS also expands descriptive characteristics beyond what patient data typically offers to include socio-demographic, economic, housing, and mobility attributes. More detailed profiles of both VHA patient populations and veterans not enrolled in the VA system will provide a comprehensive picture of groups that may benefit from interventions or outreach. As an initial exercise for using synthetic populations to measure veteran travel burdens to VA care, we apply Oak Ridge National Laboratory’s (ORNL) UrbanPop capability to generate a series of synthetic VHA patient populations for 9 Veterans Integrated Services Networks (VISN) market areas in 9 Census Divisions across the continental United States, which are listed in Table 1. We use UrbanPop to produce synthetic populations for the VISN markets selected for each US Census Division, then assign VA outpatient clinic destinations to synthetic VHA patients based on travel about each VISN market’s road network. To demonstrate using the synthetic populations to evaluate healthcare travel burdens, we compare the time-based impedance between simulated home locations and VA outpatient clinics in each VISN market. We then perform validation exercises on the synthetic populations with respect to neighborhood (block group) demographic composition as well as patient mobility, comparing aggregate origin-destination statistics for the synthetic population to outpatient visits available in restricted patient data from the VA’s Corporate Data Warehouse (CDW) database.

97 MATHEMATICS AND COMPUTING↗

Method for Calculation of Laminar Heat Transfer in Air Flow Around Cylinders of Arbitrary Cross Section (including Large Temperature Differences and Transpiration Cooling)

The solution of heat-transfer problems has become vital for many aeronautical applications. The shapes of objects to be cooled can often be approximated by cylinders of various cross sections with flow normal to the axis as, for instance heat transfer on gas-turbine blades and on air foils heated for deicing purposes. A laminar region always exists near the stagnation point of such objects. A method previously presented by E. R. G. Eckert permits the calculation of local heat transfer around the periphery of cylinders of arbitrary cross section in the laminar region for flow of a fluid with constant property values with an accuracy sufficient for engineering purposes. The method is based on exact solutions of the boundary-layer equations for incompressible wedge-type flow and on the postulate that at any point on the cylinder the boundary-layer growth is the same as that on a wedge with comparable flow conditions. This method is extended herein to take into account the influence of large temperature differences between the cylinder wall and the flow as well as the influence of transpiration cooling when the same medium as the outside flow is used as coolant.

Eckert, E R↗

Griffin Software Development Plan

Griffin is a MOOSE-based reactor physics application for advanced reactor multiphysics modeling and simulation. The application is developed in a consistent multiphysics environment with strong software quality assurance. Griffin inherited most of the capabilities of MAMMOTH/Rattlesnake and is adopting the capabilities from PROTEUS that are needed in the code. The toolset includes a variety of deterministic radiation transport solvers for fixed source, k-eigenvalue, ad-joint, and subcritical multiplication, as well as transient solvers for point-kinetics, improved quasi-static, and spatial dynamics. The code contains the cross-section preparation capabilities applicable to fast and thermal reactors, including TRISO-fueled reactors. Core management capabilities include core performance, fuel depletion and shuffling, equilibrium core calculation, pebble-bed reactor run-in and equilibrium core, molten-salt reactor delayed neutron precursor drift, and control rod and drum movement with cusping correction. This software development plan presents the current and future capabilities and features in Griffin for the design and analysis of non-light-water reactor systems in steady-state and transient conditions.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Eigenstate entanglement scaling for critical interacting spin chains

With increasing subsystem size and energy, bipartite entanglement entropies of energy eigenstates cross over from the groundstate scaling to a volume law. In previous work, we pointed out that, when strong or weak eigenstate thermalization (ETH) applies, the entanglement entropies of all or, respectively, almost all eigenstates follow a single crossover function. The crossover functions are determined by the subsystem entropy of thermal states and assume universal scaling forms in quantum-critical regimes. This was demonstrated by field-theoretical arguments and the analysis of large systems of non-interacting fermions and bosons. Here, we substantiate such scaling properties for integrable and non-integrable interacting spin-1/2 chains at criticality using exact diagonalization. In particular, we analyze XXZ and transverse-field Ising models with and without next-nearest-neighbor interactions. Indeed, the crossover of thermal subsystem entropies can be described by a universal scaling function following from conformal field theory. Furthermore, we analyze the validity of ETH for entanglement in these models. Even for the relatively small system sizes that can be simulated, the distributions of eigenstate entanglement entropies are sharply peaked around the subsystem entropies of the corresponding thermal ensembles.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗