Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “regional models”

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

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

At least 1,117 records · Page 62

Effect of convective transport in edge and scrape-off-layer plasmas of ADITYA-U tokamak

The 2-D edge plasma fluid transport code, UEDGE has been used to simulate the edge and scrape-off layer regions of circular limiter plasmas of ADITYA-U. The computational mesh defining the limiter geometry of ADITYA-U has been introduced in the UEDGE code through an in-house developed routine. It has successfully been integrated with the UEDGE code and the measured radial profile of electron density (n e ) in the edge region is modelled using the simulations. The plasma discharge for the present study has plasma current I P ~ 130 kA, duration ~ 150 ms, chord-averaged electron density ~ 1.3 x 10 19 m -3 , and central electron temperature ~ 270 eV has been observed. The maximum toroidal magnetic field is ~ 1.2 T. It has been found that an inward convective velocity, v conv ~ 1.5 m/s is required in addition to the constant perpendicular diffusion coefficient, D ⟂ ~ 0.2 m 2 /s to adequately match the measured n e profile in typical discharges of ADITYA-U. The value of D ⟂ ~ 0.2 m 2 /s is found to lie between the estimated neoclassical diffusivity and Bohm diffusivity.

Physics - Plasma physics↗

Numerical discreteness errors in multispecies cosmological N -body simulations

ABSTRACT We present a detailed analysis of numerical discreteness errors in two-species, gravity-only, cosmological simulations using the density power spectrum as a diagnostic probe. In a simple set-up where both species are initialized with the same total matter transfer function, biased growth of power forms on small scales when the solver force resolution is finer than the mean interparticle separation. The artificial bias is more severe when individual density and velocity transfer functions are applied. In particular, significant large-scale offsets in power are measured between simulations with conventional offset grid initial conditions when compared against converged high-resolution results where the force resolution scale is matched to the interparticle separation. These offsets persist even when the cosmology is chosen so that the two particle species have the same mass, indicating that the error is sourced from discreteness in the total matter field as opposed to unequal particle mass. We further investigate two mitigation strategies to address discreteness errors: the frozen potential method and softened interspecies short-range forces. The former evolves particles under the approximately ‘frozen’ total matter potential in linear theory at early times, while the latter filters cross-species gravitational interactions on small scales in low-density regions. By modelling closer to the continuum limit, both mitigation strategies demonstrate considerable reductions in large-scale power spectrum offsets.

79 ASTRONOMY AND ASTROPHYSICS↗

Pressure-induced polyamorphic transition in CaA l 2 O 4 glass

In situ high-pressure ultrasonic velocity measurements of CaAl 2 O 4 glass reveal abrupt irreversible discontinuities in the elastic wave velocities at ~8–10 GPa. Total structure factor and pair distribution functions measured by synchrotron x-ray diffraction show a rapid change in the intermediate range structure attributed to a rearrangement of calcium ions over this narrow pressure region. Atomistic models obtained from molecular dynamics simulations reveal that this intermediate range structure is explained by a transition of Ca–O void radius distribution from a bimodal distribution with peaks at ~2.1 and ~2.4 Å to a single distribution centered at ~2.1 Å. Further, these abrupt structural changes involving the rapid increase in elastic wave velocity are markedly different to the continuous transformations observed in conventional network-forming glasses, such as SiO 2 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Measurements of the 96 Zr(α,n) cross section for astrophysics and applications

The reaction 96 Zr(α,n) plays an important role in ν-driven wind nucleosynthesis in core-collapse supernovae and is a possible avenue for medical isotope production. Here, cross section measurements were performed using the activation technique at the Edwards Accelerator Laboratory. Results were analyzed along with world data on the 96 Zr(α, n) cross section and 96 Zr(α, α) differential cross section using large-scale Hauser-Feshbach calculations. We find general agreement between our data, previous measurements, and a statistical description of the reaction. However, we find a larger cross section at low energies compared to prior experimental results, allowing for a larger astrophysical reaction rate. This may impact results of core-collapse supernova ν-driven wind nucleosynthesis calculations, but does not significantly alter prior conclusions about 99 Mo production for medical physics applications. The results from our large-scale Hauser-Feshbach calculations demonstrate that phenomenological optical potentials may yet be adequate to describe (α, n) reactions of interest for ν-driven wind nucleosynthesis, albeit with regionally-adjusted model parameters.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Automated label flows for excited states of correlation functions in lattice gauge theory

Extracting excited states from lattice gauge theory correlation functions can be achieved through $χ^2$ minimization fits or algebraic approaches such as the variational method and Prony's method. Performing any kind of error analysis often leads to overlapping confidence regions of model parameters, even when the spectrum is not particularly dense. To correctly estimate errors, one must beware of mislabeling the states. In this work, we provide an algorithm that we call automated label flows which consistently and systematically identifies a deterministic labeling of states. This is a black-box approach in the sense that it gives a sensible set of labels without user guidance. Further, as an example, we pair one black box method with another, analyzing a lattice correlation function from real data using automated label flows in the context of Prony's method.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

First Measurement of Missing Energy due to Nuclear Effects in Monoenergetic Neutrino Charged-Current Interactions

We present the first measurement of the missing energy due to nuclear effects in monoenergetic, muon neutrino charged-current interactions on carbon, originating from 𝐾 + → 𝜇 + ⁢𝜈 𝜇 decay at rest (𝐸 𝜈 𝜇 = 235.5 MeV), performed with the J-PARC Sterile Neutrino Search at the J-PARC Spallation Neutron Source liquid scintillator based experiment. Toward characterizing the neutrino interaction, ostensibly 𝜈 𝜇 ⁢𝑛 → 𝜇 − ⁢𝑝 or 𝜈 𝜇 ⁢ 12 C → 𝜇 − ⁢ 12 N, we define the missing energy as the energy transferred to the nucleus (𝜔) minus the kinetic energy of the outgoing proton(s), 𝐸 𝑚 ≡ 𝜔−∑ 𝑇 𝑝 , and relate this to visible energy in the detector, 𝐸 𝑚 = 𝐸 𝜈 𝜇 ⁡(235.5 MeV) − 𝑚 𝜇⁡ (105.7 MeV) + [𝑚 𝑛 − 𝑚 𝑝⁡ (1.3 MeV)] − 𝐸 vis . The missing energy, which is naively expected to be zero in the absence of nuclear effects (e.g., nucleon separation energy, Fermi momenta, and final-state interactions), is uniquely sensitive to many aspects of the interaction, and has previously been inaccessible with neutrinos. The shape-only, differential cross section measurement reported, based on a (77 ± 3)% pure double-coincidence kaon decay-at-rest signal (621 total events), provides detailed insight into neutrino-nucleus interactions, allowing even the nuclear orbital shell of the struck nucleon to be inferred. The measurement provides an important benchmark for models and event generators at hundreds of MeV neutrino energies, characterized by the difficult-to-model transition region between neutrino-nucleus and neutrino-nucleon scattering, and relevant for applications in nuclear physics, neutrino oscillation measurements, and Type-II supernova studies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Record High 2022 September-Mean Temperature in Western North America

Human-induced warming is estimated to have increased occurrence probability (magnitude) of the record-breaking September 2022 heat event in western North America by 6–67 times (0.6–1 K) by E3SMv2 and even higher by coupled regional refined model (RRM) simulations.

54 ENVIRONMENTAL SCIENCES↗

Forest Breeze–Cold Pool Interactions Drive Convective Organization over Heterogeneous Vegetation

Heterogeneous landscapes can influence the development of convection through the generation of thermally driven mesoscale circulations. To assess the impacts of these circulations and their interaction with sea breezes, we simulated convection in an idealized coastal environment using the Regional Atmospheric Modeling System (RAMS). We compared simulations with striped patterns of surface vegetation to those of uniform vegetation to identify the importance of vegetation heterogeneity in impacting convective development. Under dry soil conditions representative of those during the Tracking Aerosol Convection Interactions Experiment (TRACER) and Experiment of Sea Breeze Convection, Aerosols, Precipitation, and Environment (ESCAPE) campaigns in June 2022, we found that these vegetation-induced circulations, referred to in the literature as “forest breezes,” are more important than the sea breeze in determining the location of convection initiation. Convection and precipitation are also found to be favored over forests and suppressed over pasture and suburban landscapes as a result of greater surface sensible heat flux over the forest. Our findings also indicate that forest breezes are important for initiating convection along the boundaries of the forest, but that cold pools may play a key role in propagating the forest breezes toward the center of the forest stripe. In our simulations, the collisions of these breezes in the center of the forest stripe lead to uplift and strong convection there; however, a different width of the forest stripe would alter when the forest breezes collide or whether they collide at all. The presence of these cold pools may therefore impact the “ideal stripe width,” the width of each vegetation stripe which maximizes domain-wide precipitation.

54 ENVIRONMENTAL SCIENCES↗

2014 Southern Nevada Household Travel Survey

The 2014 Southern Nevada Household Travel Survey collected information from residents in the Las Vegas area to update the regional travel demand model and assess travel behavior. The Regional Transportation Commission of Southern Nevada contracted with Westat to conduct the survey. The survey was conducted in two phases—from March to May 2014 and from August to October 2014. Participants provided demographic and travel data (via a travel log). A 10% subsample (1,694 participants) was randomly selected to take part in a wearable global positioning system (GPS) technology-based component of the study, the purpose of which was to assess the level of trip under-reporting in the self-reported travel logs. Participants in the GPS portion of the study were instructed to record trips in their travel logs on the first day only, while passively recording their travel for three full days.

1Hz data↗

Urban Parameters Maricopa County 100m Grid Spacing

132 urban parameters for Maricopa County are provided in WRF-readable binary form at 100m grid spacing with index. These parameters allow regional weather modelers to account for the contribution of built structures in modeling the urban heat island. This data is also provided in CSV format. Data are generated using 2023 Model America Version 1 data (https://www.osti.gov/biblio/1774134) and will be used to represent the urban morphology in Maricopa County to the WRF model for experiments conducted for the Southwest Urban Corridor Integrated Field Laboratory (SW-IFL). Three files associated with these parameters are included in this repository: index, binary file 00001-01549.00001-01796, Maricopa_Parameters_03072024.csv . To use the binary file with the Weather Research and Forecasting (WRF) model, the binary file and the index file must be placed in their own directory in WRF_GEOG and accessed in the same way the National Urban Database and Access Portal Tool data (NUDAPT44) would be accessed. The name of the binary file must be in the format: 00001-02715.00001-01803. If, upon download, the binary filename has additional or different characters, please change them to match this format.

54 ENVIRONMENTAL SCIENCES↗

A User’s Guide to the PLTEMP/ANL Code

PLTEMP/ANL V4.3 is a program that obtains a steady-state flow and temperature solution for a nuclear reactor core, or for a single fuel assembly. It is based on an evolutionary sequence of codes originally used for plate temperatures, hence “PLTEMP”, developed at Argonne National Laboratory over several decades. Fueled and non-fueled regions are modeled. Each fuel assembly consists of one or more plates or tubes separated by coolant channels. The fuel plates may have one to five layers of different materials, each with heat generation. The width of a fuel plate may be divided into multiple longitudinal stripes, each with its own axial power shape. Depending upon the heat transfer method selected, the temperature solution is effectively 2-dimensional or 3-dimensional (Appendix XV). It begins with a solution across all coolant channels and fuel plates or tubes within a given fuel assembly, at the entrance to the assembly. The temperature solution is repeated for each axial node along the length of the fuel assembly. The geometry may be either slab or radial, corresponding to fuel assemblies made of a series of flat (or slightly curved) plates, or of nested tubes. A variety of thermal-hydraulic correlations are available to determine safety margins such as onset-of-nucleate boiling ratio(ONBR), departure from nucleate boiling ratio (DNBR), and onset of flow instability ratio (OFIR). Coolant properties for either light or heavy water are obtained from FORTRAN functions rather than from tables. The code is intended for thermal-hydraulic analysis of research reactor performance in the sub-cooled boiling regime. Both turbulent and laminar flow regimes can be modeled. Options to calculate both forced flow and natural circulation are available. A general search capability is available (Appendix XII) to greatly reduce the reactor analyst’s time.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

A User’s Guide to the PLTEMP/ANL Code (V.4.3)

PLTEMP/ANL V4.3 is a program that obtains a steady-state flow and temperature solution for a nuclear reactor core, or for a single fuel assembly. It is based on an evolutionary sequence of codes originally used for plate temperatures, hence “PLTEMP”, developed at Argonne National Laboratory over several decades. Fueled and non-fueled regions are modeled. Each fuel assembly consists of one or more plates or tubes separated by coolant channels. The fuel plates may have one to five layers of different materials, each with heat generation. The width of a fuel plate may be divided into multiple longitudinal stripes, each with its own axial power shape. Depending upon the heat transfer method selected, the temperature solution is effectively 2-dimensional or 3-dimensional (Appendix XV). It should be noted that for a 3-dimensional solution the axial power shape is modeled without accounting for axial heat conduction. It begins with a solution across all coolant channels and fuel plates or tubes within a given fuel assembly, at the entrance to the assembly. The temperature solution is repeated for each axial node along the length of the fuel assembly. The geometry may be either slab or radial, corresponding to fuel assemblies made of a series of flat (or slightly curved) plates, or of nested tubes. A variety of thermal-hydraulic correlations are available to determine safety margins such as onset-of-nucleate boiling ratio(ONBR), departure from nucleate boiling ratio (DNBR), and onset of flow instability ratio (OFIR). Coolant properties for either light or heavy water are obtained from FORTRAN functions rather than from tables. The code is intended for thermal-hydraulic analysis of research reactor performance in the sub-cooled boiling regime. Both turbulent and laminar flow regimes can be modeled. Options to calculate both forced flow and natural circulation are available. A general search capability is available (Appendix XII) to greatly reduce the reactor analyst’s time.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

A Tutorial for Using an Open-Source Solver for the Regional Energy Deployment System (ReEDS) Model

ReEDS is a publicly available model developed at NREL that can be used to analyze the potential evolution of the U.S. electric power system into the future. ReEDS is formulated as a linear program, written in GAMS, and solved using a linear programming solver (e.g., CPLEX). This tutorial provides context and understanding for the implications of using an open-source solver for ReEDS.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Assessment of MOOSE-Based Tools for Calculating Radial Core Expansion

Radial core expansion in liquid-metal cooled fast reactor systems is a well-known phenomenon that produces strong reactivity feedback effects. An inherently safe reactor design takes advantage of negative reactivity feedback in accident scenarios by utilizing a core restraint system which produces a bowed shape that allows for radial expansion of the fuel regions. Detailed modeling and simulation of core radial expansion itself as well as subsequent reactivity feedback is a challenging task involving contact of many fuel assembly elements and physics feedback from neutronics, thermal hydraulics, and thermal mechanical response. A variety of physics codes have been developed to model aspects of radial core expansion but in general invoke geometrical or physics approximations. No code system currently exists which tightly and robustly couples these physics with enough detail to fully resolve the complex core radial expansion reactivity feedback effects. The future availability of such a code system is of vital importance to fully understanding the reactivity feedback effects that occur due to radial expansion, and consequently to optimizing the design of the core restraint system. A high-fidelity code will also be used to benchmark existing lower fidelity, faster running models to understand their benefits, limitations, and range of applications. A code development path using MOOSE-based tools is proposed in order to leverage the detailed geometry capabilities and natural tight coupling and robustness of MOOSE-based applications for modeling this complex phenomena. While simulation of the full phenomenon involves several physics, an assessment has been initiated on the capabilities and readiness of the currently available Tensor Mechanics module within MOOSE for calculation of the structural mechanical responses which occur within the reactor core. This report focuses on modeling the force-deformation response which mimics the physics of a duct contact deformation, as well as differential thermal expansion which produces a thermal bowed shaped for the fuel assemblies in a core. Simple examples were initially performed such as simple supported beam bending under load. The complexity of examples was progressively increased to better mimic the duct behavior by including a differential thermal example and inclusion of hexagonal cross-sections in the geometry. Further assessment of the structural mechanical response simulation capability is still required for modeling duct contact interactions and irradiation creep and swelling. Companion thermal hydraulic and neutronics assessments will also be required; these activities are planned for future years. Finally, integration of the multiple physics components through MOOSE is required to predict the core radial expansion and subsequent feedback effects.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Investigating Benefits and Challenges of Converting Retiring Coal Plants into Nuclear Plants

A coal-to-nuclear (C2N) transition means siting a nuclear reactor at the site of a recently retired coal power plant. Three overarching questions from the C2N transition guide this research: where in the United States are retired coal facilities located and what factors make a site feasible for transition; what factors of technology, cost, and project timeline drive investor economics over such a decision; and how will C2N impact local communities? The study team evaluated the siting characteristics of recently retired plants and those operating coal-fired power plant sites run by a utility or an independent power producer utilizing publicly available data to screen U.S. coal power plant sites to nuclear-feasible locations. After screening all retired coal sites to a set of 157 potential candidates and screening operating sites to a set of 237 candidates, the study team estimates that 80% of retired and operating coal power plant sites that were evaluated have the basic characteristics needed to be considered amenable to host an advanced nuclear reactor. For the recently retired plant sites evaluated, this represents a capacity potential of 64.8 GWe to be backfit at 125 sites. For the operating plant sites evaluated, this represents a capacity potential of 198.5 GWe to be backfit at 190 sites. This report evaluates a case study for the detailed impacts and potential outcomes from a C2N transition. Based on the nuclear technology choices and sizes evaluated to replace a large coal plant of 1,200 MWe generation capacity at the case study site, nuclear overnight costs of capital could decrease by 15% to 35% when compared to a greenfield construction project, through the reuse of infrastructure from the coal facility. Nuclear replacement designs can have a lower capacity size because nuclear power plants run at higher capacity factors than coal power plants. In the case study replacing coal capacity with 924 GWe of nuclear capacity, the study team found regional economic activity could increase by as much as $\$275 million$ and add 650 new, permanent jobs to the region of analysis. The evaluated site choice in the report is hypothetical for analysis purposes only and based on available data and documented assumptions. Consequently, the findings only inform at a general level. A community, investor, or other interested stakeholder can use these results to set up a detailed, in-depth analysis for a specific application of interest, such as evaluating a C2N transition of a specific coal power plant and a specific nuclear technology design. The report was subjected to independent peer reviews by experts in systems engineering and regional economic modeling to evaluate analysis and assumptions.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Progress on Demonstration of a MOOSE-Based Coupled Capability for Hot Channel Factors in Fast Reactors

Hot channel factors (HCFs) are computed values that account for the impact on predicted peak fuel, cladding, and coolant temperatures due to uncertainties in the as-built reactor’s material properties and geometry as well as uncertainties due to modeling approximations. Reduction in computed HCF values via reduction or elimination of modeling approximations may translate to significant economic savings if the reactor power can be raised due to the extra temperature margin gained. While limited historical datasets exist for sodium-cooled fast reactors (SFRs), there are no available HCF data for lead-cooled fast reactors (LFRs) outside of work generated previously within NEAMS. The computation of HCFs involves insights from reactor physics, thermal fluids and heat conduction calculations to determine how the peak temperatures respond to various uncertainties in the design. Due to the significant advantages for multi-physics coupling offered by the MOOSE framework, Griffin (MOOSE-based reactor physics code), MOOSE Heat Conduction Module, and Cardinal (MOOSE-wrapped multi-physics application which includes the NekRS thermal fluids code) are being coupled together using the MOOSE MultiApp System to develop a highfidelity multi-physics modeling capability for HCF simulations. This high-fidelity coupling workflow may also be beneficial for other fast reactor applications in the future. In previous work, Griffin and NekRS were individually assessed to ensure the necessary capabilities were in place. This work describes initial efforts to couple the codes (including folding in the MOOSE Heat Conduction Module) and determining the workflow for the perturbed calculations which will leverage the Stochastic Tools Module (STM). To our knowledge, this is the first coupling of Griffin and NekRS as well as the first exploratory use of Stochastic Tools Module for Cardinal. In this report, the neutronics code Griffin, the heat conduction solver in MOOSE, and the MOOSE-wrapped application containing NekRS (Cardinal) are linked together to demonstrate the coupled capability. Griffin and Cardinal are linked dynamically by specifying shared libraries. Different coupling hierarchies are tested for selecting the most appropriate coupling strategy. A coupling scheme is selected based on the efficiency of calculation and ease of data communication. Multiple tests are performed to choose suitable mesh structure, model configurations, scheme setup and boundary conditions to avoid loss of energy due to data interpolation between different modules or weak imposition of fluxes in finite element codes. Computational experiments are performed to study the tolerance control of each type of iteration to avoid false convergence. The coupled capability is demonstrated in both single pin and 7-pin models based on LFR materials and geometry. The study finds that the use of too large a time step size in the heat conduction module can lead to temperature oscillation even though the heat conduction equation does not have a time-derivative kernel, but only the time-dependent boundary condition. A 7-pin model without duct region achieved good convergence in the coupled calculation while a 7 pin model with duct region experienced data communication issues which need to be resolved.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Managing Increased Electric Vehicle Shares on Bulk Power Systems

The convergence of transportation electrification, other load growth, and the integration of many power sector technologies presents a complex planning problem requiring realistic, region-specific modeling and analysis to ensure the power system can support new loads quickly and affordably.

33 ADVANCED PROPULSION SYSTEMS↗