Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Neutron Multiplication”

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 847 records · Page 47

Neutron-unbound states in 31 Ne

Background: The Island of Inversion near the N = 20 shell gap is home to nuclei with reordered single-particle energy levels compared to the spherical shell model. Studies of 31 Ne have revealed that its ground state has a halo component, characterized by a valence neutron orbiting a deformed 30 Ne core. This lightly-bound nucleus with a separation energy of only S n = 170 keV is expected to have excited states that are neutron unbound. Purpose: The purpose of this experiment was to investigate the low-lying excited states in 31 Ne that decay by the emission of a single neutron. Methods: An 89 MeV/u 33 Mg beam impinged on a segmented Be reaction target. Neutron unbound states in 31 Ne were populated via a two-proton knockout reaction. The 30 Ne fragment and associated neutron from the decay of 31 Ne were detected by the MoNA-LISA-Sweeper experimental setup at the National Superconducting Cyclotron Laboratory. Invariant mass spectroscopy was used to reconstruct the two-body decay energy ( 30 Ne + n). Results: The two-body decay energy spectrum exhibits two features: a low-lying peak at 0.30 ± 0.17 MeV and a broad enhancement at 1.50 ± 0.33 MeV, each fit with an energy-dependent asymmetric Breit-Wigner line shape representing a resonance in the continuum. Accompanying shell model calculations using the FSU interaction within NuShellX, combined with cross-section cal culations using the eikonal reaction theory, indicate that these peaks in the decay energy spectrum are caused by multiple resonant states in 31 Ne. Conclusions: In this study, excited states in 31 Ne were observed for the first time. Transitions from calculated shell model final states in 31 Ne to bound states in 30 Ne are in good agreement with the measured decay energy spectrum. Cross-section calculations for the two-proton knockout populating 31 Ne states as well as spectroscopic factors pertaining to the decay of 31 Ne into 30 Ne are used to examine the results within the context of the shell model expectations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Applications of amorphous track models in radiation biology

The average or amorphous track model uses the response of a system to gamma-rays and the radial distribution of dose about an ion's path to describe survival and other cellular endpoints from proton, heavy ion, and neutron irradiation. This model has been used for over 30 years to successfully fit many radiobiology data sets. We review several extensions of this approach that address objections to the original model, and consider applications of interest in radiobiology and space radiation risk assessment. In the light of present views of important cellular targets, the role of target size as manifested through the relative contributions from ion-kill (intra-track) and gamma-kill (inter-track) remains a critical question in understanding the success of the amorphous track model. Several variations of the amorphous model are discussed, including ones that consider the radial distribution of event-sizes rather than average electron dose, damage clusters rather than multiple targets, and a role for repair or damage processing.

NASA Discipline Radiation Health↗

Phase lags in Cygnus X-1

Coronae of hot electrons are ubiquitous features in models of the inner regions of accreting black holes and nonmagnetic neutron stars. The scattering optical depth of these coronae inferred from observations is often tau approximately 3, so the energy spectrum of the disk in these regions is likely masked by the spectrum of upscattered photons in the corona and some of the disk properties are thus obscured. Observations of the dependence on not available with the energy spectrum. In our picture, the disk emission is modeled as a photon source injected into a Comptonizing corona; Comptonization may also occur in the disk, but for our purposes the disk emission is simply an input to the corona. We show that, contrary to some claims, the functional dependence of lag on Fourier frequency emerges intact from transit through the corona, modulo a multiplicative factor (which may in principle be negative, so that a phase lag can be changed to a phase lead), even if the properties of the corona vary with time. We also show that any frequency dependence of the lag due to variation in the corona itself is only second order in the amplitude of the variation, and cannot exceed the transit time approximately ms of the corona; thus, the lags of up to 0.1 s seen in several black hole candidates come from lags in the emission from the disk. Finally, we predict that plots of the time lag versus Fourier frequency in black hole candidates should have a 'shelf' of constant lag equal to the coronal lag (approximately 1 ms), with the constant being proportional to In (E(sub 2)/E(sub 1)) for the lag between energies E(sub 2) and E(sub 1). The lack of such a shelf in current observations of several galactic black hole candidates constraints the radii of the coronae to be R less than or approximately = 10(exp 8) tau cm.

Miller, M. Coleman↗

Mixed-field GCR Simulations for Radiobiological Research Using Ground Based Accelerators

Space radiation is comprised of a large number of particle types and energies, which have differential ionization power from high energy protons to high charge and energy (HZE) particles and secondary neutrons produced by galactic cosmic rays (GCR). Ground based accelerators such as the NASA Space Radiation Laboratory (NSRL) at Brookhaven National Laboratory (BNL) are used to simulate space radiation for radiobiology research and dosimetry, electronics parts, and shielding testing using mono-energetic beams for single ion species. As a tool to support research on new risk assessment models, we have developed a stochastic model of heavy ion beams and space radiation effects, the GCR Event-based Risk Model computer code (GERMcode). For radiobiological research on mixed-field space radiation, a new GCR simulator at NSRL is proposed. The NSRL-GCR simulator, which implements the rapid switching mode and the higher energy beam extraction to 1.5 GeV/u, can integrate multiple ions into a single simulation to create GCR Z-spectrum in major energy bins. After considering the GCR environment and energy limitations of NSRL, a GCR reference field is proposed after extensive simulation studies using the GERMcode. The GCR reference field is shown to reproduce the Z and LET spectra of GCR behind shielding within 20% accuracy compared to simulated full GCR environments behind shielding. A major challenge for space radiobiology research is to consider chronic GCR exposure of up to 3-years in relation to simulations with cell and animal models of human risks. We discuss possible approaches to map important biological time scales in experimental models using ground-based simulation with extended exposure of up to a few weeks and fractionation approaches at a GCR simulator.

Kim, Myung-Hee Y.↗

Multidimensional Rietveld refinement of high-pressure neutron diffraction data of PbNCN

High-pressure neutron powder diffraction data from PbNCN were collected on the high-pressure diffraction beamline SNAP located at the Spallation Neutron Source (SNS) of Oak Ridge National Laboratory (Tennessee, USA). The diffraction data were analyzed using the novel method of multidimensional (two dimensions for now, potentially more in the future) Rietveld refinement and, for comparison, employing the conventional Rietveld method. To achieve two-dimensional analysis, a detailed description of the SNAP instrument characteristics was created, serving as an instrument parameter file, and then yielding both cell and spatial parameters as refined under pressure for the first time for solid-state cyanamides/carbodiimides. The bulk modulus B 0 = 25.1 (15) GPa and its derivative B ′ 0 = 11.1 (8) were extracted for PbNCN following the Vinet equation of state. Surprisingly, an internal transition was observed beyond 2.0 (2) GPa, resulting from switching the bond multiplicities (and bending direction) of the NCN 2− complex anion. The results were corroborated using electronic structure calculation from first principles, highlighting both local structural and chemical bonding details.

Meinerzhagen, Yannick↗

Geologic Characterization of the South Georgia Rift Basin for Source Proximal CO2 Storage

The project Geologic Characterization of the South Georgia Rift Basin for Source Proximal CO2 Storage is one of 9 site characterization projects that were implemented as part of ARRA (American Recovery and Reinvestment Act). Data from this project was used to improve resolution of data in NATCARB in the area of study. Data related to this study has already been incorporated in NATCARB Atlas. The South Carolina Research Foundation and partners evaluated the feasibility of CCS in the Jurassic/ Triassic (J / TR) saline formations of the buried Mesozoic South Georgia Rift (SGR) Basin that extends from South Carolina into Georgia. The J / TR sequence, based on preliminary assessment of limited geologic and geophysical data, appears to have both the appropriate areal extent and multiple horizons to permanently and safely store CO2 The presence of several igneous rock layers within the sequence may potentially provide adequate seals to prevent upward CO2 migration into the Coastal Plain aquifer systems. Approximately 81 kilometers of 2-D seismic reflection data were collected by Bay Geophysical, Inc. to explore a portion of the SGR located in southern Georgia. The 81 kilometers were divided into two lines approximately 40.5 kilometers each, with Line 1 intersecting Georgia well GGS 3457. Line 2 intersects Line 1 at the southern portion of Line 1 to maximize the extent of coverage away from GGS-3457 (a deep well drilled in the 1980s for oil and gas exploration). This well had a set of usable logs, including gamma and neutron logs that provided promising results related to CO2 storage. Results showed sandstone with porosity values greater than 10 percent and a thickness of 120 meters. The design of the seismic shot was to extrapolate information away from the well and to better define the extent of the SGR and the necessary reservoir and caprock for a successful CO2 injection. A numerical simulation model of CO2 Injection and migration was developed based on the geology log for the GGS-3457 well. The simulation model was used to investigate the feasibility of injecting 30 million metric tons of CO2 into SGR J / TA sediments and integrity of the diabase layers as seals to prevent CO2 migration.

2-D seismic↗

Foundations of Molecular 'Isotomics'

The naturally occurring rare isotopes are versions of common elements, such as hydrogen, carbon and oxygen, that contain a larger than usual number of neutrons in their atomic nuclei and therefore are higher in mass than the common atoms of that element. Isotopes exist for most elements and are found in most natural and synthetic materials, but are uneven in their distribution because chemical and physical processes are isotope-selective (e.g., a chemical reaction may proceed more rapidly for one isotope than for another). For this reason, abundances of isotopes in a material of interest can provide a record, or ‘signature’ of various features of that material’s origin and history. These signatures have been used in the geo, life, chemical and physical sciences in a wide variety of ways over close to 8 decades. However, many such applications struggle to reach unique interpretations of isotopic data because multiple factors combine to control a given sample’s overall isotopic content. That is, the factors controlling isotopic content are too numerous and complex to fully constrain from a simple measurement of a material’s isotope abundances. However, the distribution of isotopes within materials, at molecular scales potentially provides a vastly larger number and diversity of constraints on the chemical and physical processes that comprise a material’s history. The rare isotopes may be concentrated into one atomic position in a molecule relative to another, some proportion of molecules in a sample may contain two or more rare isotopes, and those multiply-isotope-substituted forms of molecules may also have uneven distributions of those isotopes across individual atomic sites. For these reasons, even small, seemingly simple molecules, such as sugars, amino acids or drug compounds, actually exist in a vast number of isotopically unique forms (often millions or more), and each one of those forms is in some sense an independent ‘vote’ on that sample’s history. This project has focused on opening this rich archive of information by enabling the creation of routinely and widely applicable ways of measuring and interpreting isotopic structures of molecules. This work has included the development of core technologies and analytical methods, advancing fundamental understanding of the physical and chemical properties of isotopic versions of molecules, and conducting proof of concept studies of illustrative geochemical, cosmochemical and forensic problems in order to show how these technologies, methods and principles come together to solve problems in new ways. A key to the success of this project was the adaptation of ‘Fourier transform mass spectrometry’ (FTMS) to the task of precisely measuring proportions of the rare, naturally occurring isotopic forms of molecules. FTMS is a highly specialized form of mass spectrometry that traps ions within magnetic or electrostatic cavities and, effectively, ‘listens’ (through registering of subtle electrical signals) to the harmonic signals they make while rapidly orbiting within those cavities. These signals have periods that are a function of their mass and strength (or ‘loudness’) that is proportional to their abundances. Thus, these signals constrain relative amounts of molecules that differ in their mass due to various isotopic substitutions. This technology has been essential to the identification of organic molecules in the life, chemical and environmental sciences for over 4 decades, but generally has lacked the control, stability and precision to meaningfully measure rare isotope forms of molecules. This project’s most fundamental contribution has been to modify FTMS, both in terms of hardware and methods, to enable such measurements. The raw data of molecular isotopic structure is tremendously voluminous and complex, so another important activity of this project has been developing the theoretical and data-science tools needed to interpret the data generated by this new form of isotopic measurement. A particularly challenging part of this task has been predicting molecular isotopic structure, as only through the comparison of measurements with predictions can we make progress on hypothesis driven research questions. We have attacked this this prediction task through a combination of first-principles chemical-physics models of the effects of isotope substitution on molecule properties and data-science models that permit us to generalize that chemical physics to cases that have not yet been studied by detailed chemical physics theory. The proof of concept applications we have pursued over the course of this study include biological reactions of amino acids and other biomolecules, non-biological synthesis of organic molecules in extra-terrestrial settings such as meteorites, petroleum geoscience questions concerning the origin and evolution of natural gas, oil and kerogen compounds, and forensic questions such as the sourcing of chemical weapons. The successes of these applications have laid the groundwork for the next phase of this field’s development, which will include larger scale and more ambitious studies of molecular isotopic structure as a means of diagnosing human diseases, such as cancer, and reconstructing detailed interpretations of the origin and evolution of organic molecules in modern and geological environments.

Cesar, Jaime↗

The spectra of accretion disks and their application to low-mass X-ray binaries

The influence of external radiation from the boundary layer and/or from a neutron star surface on the emitted spectrum of an accretion disk is studied. The surface properties of the accretion disk are reviewed in the context of the alpha disk model, and the spectral regimes associated with each are discussed. It is found that the spectrum differs significantly from that obtained in the blackbody assumption, especially in cases where the disk is irradiated at fluxes near the Eddington limit. Generally, the spectrum is found to peak at lower energies than the blackbody assumption would indicate and is considerably flatter for energies between about 0.1 and 1 keV. For an irradiated disk where the color temperature characterizing the external radiation is greater than about 1 keV, the spectrum emitted from the disk is harder than that of multiple blackbodies. The emitted spectrum from the disk is sensitive to the inclination angle of the system and the location of the inner disk edge.

Taam, Ronald E.↗

Exploring the NΛ–NΣ coupled system with high precision correlation techniques at the LHC

The interaction of Λ and Σ hyperons (Y) with nucleons (N) is strongly influenced by the coupled-channel dynamics. Due to the small mass difference of the NΛ and NΣ systems, the sizable coupling strength of the NΣ↔NΛ processes constitutes a crucial element in the determination of the NΛ interaction. In this letter we present the most precise measurements on the interaction of pΛ pairs, from zero relative momentum up to the opening of the NΣ channel. The correlation function in the relative momentum space for pΛ ⊕$\bar{p}$$\bar{Λ}$ pairs measured in high-multiplicity triggered pp collisions at $\sqrt{s}$=13 TeV at the LHC is reported. The opening of the inelastic NΣ channels is visible in the extracted correlation function as a cusp-like structure occurring at relative momentum k*=289 MeV/c. This represents the first direct experimental observation of the NΣ↔NΛ coupled channel in the pΛ system. The correlation function is compared with recent chiral effective field theory calculations, based on different strengths of the NΣ↔NΛ transition potential. A weaker coupling, as possibly supported by the present measurement, would require a more repulsive three-body NNΛ interaction for a proper description of the Λ in-medium properties, which has implications on the nuclear equation of state and for the presence of hyperons inside neutron stars.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Multiple valence-space and core-excited band terminations in Ho and Er isotopes

High-spin structures of holmium (𝑍 = 67) and erbium (𝑍 = 68) isotopes are presented near mass 160, a region where collective and single-particle modes of excitation compete in the generation of angular momentum at the highest experimentally attainable spins in atomic nuclei. The level schemes of the 156,157,158 Ho isotopes have been significantly extended, reaching spins in the 40⁢ℏ to 55⁢ℏ range, through multiple experiments utilizing the Gammasphere spectrometer. In particular, several new noncollective band-terminating states are reported for these nuclei. The nature of these states, including both valence-space terminations and core-excited oblateconfigurations, is discussed in terms of proton particle-hole (p-h) excitations across the semimagic 𝑍 = 64 shell closure. Cranked-Nilsson–Strutinsky (CNS) calculations have been used to interpret the nature of these noncollective states. The corresponding terminating states in 154–159 Er are also included, providing a systematic description of band termination in this mass region. In addition, three weak rotational structures have been assigned to 157 Ho at the highest spins. These collective bands are proposed to be based on a triaxial strongly deformed nuclear shape, involving neutron p-h excitations across the 𝑁 = 82 shell gap.

150 ≤ A ≤ 189↗

Quantum Choreography of the Nucleus: Rotations, Vibrations, and Emergent Structure

Abstract Nuclei are complex many-body quantum systems where interactions of the neutrons and protons via the strong, the weak, and the electromagnetic forces lead to the emergence of simple patterns of energy states that have been described by various theoretical approaches. One of the goals of all the theoretical models is the development of a universal theory that can be applied across the entire chart of nuclides. Significant progress has been made by experiments as well as the increasing sophistication of models, but a universal theory has yet to be established. A recent review of nuclei in the Z = 50–82 region of the chart of nuclides has analyzed all the available compiled data from several decades of studies towards a clarification of the low-lying structure of nuclei. Other reviews have reported and explained the emergence of multiple different shapes in nuclei at somewhat higher excitation energies than the ground state. Some have challenged the interpretation of the first excited K π = 0 + band as a vibration of ground state. This work attempts to provide a guide to determining the nature of the first excited K π = 0 + band in nuclei by the combined use of nuclear lifetimes, energy level evolutions, dynamic moments of inertia, and intrinsic quadrupole moments extracted from transition probabilities. The result is that for a subset of the nuclei in this region, the K π = 0 + band is consistent with the traditional 𝛽-vibration description of an oscillation built on the ground state.

CHFB+5DCH↗

Probability Density Function for the spatial and intensity distribution of neutron-induced defects in Silicon

The ability to model semiconductor device degradation under neutron irradiation depends upon having a robust modeling capability for the neutron-induced collision cascades as well as a means to analytically fit the resulting probability distributions of defect production and ionizing energy deposition for purposes of extrapolation to low-probability, high-consequence scenarios. In this paper, the widely-utilized binary collision approximation codes MARLOWE and SRIM are deployed in conjunction with a critical examination of their parameterizations as benchmarked against higher-fidelity molecular dynamics simulations. A simple 3-parameter form described by the Generalized Logistic Distribution is shown to be a good fit to Frenkel pair and ionization intensity distributions in bulk silicon. The BCA codes are then applied to simulate cascades in 5 nm layers of a representative gate-all-around nanosheet transistor, where joint probability distributions of threshold levels of damage to multiple layers are evaluated.

36 MATERIALS SCIENCE↗

Pulse profile modelling of thermonuclear burst oscillations – II. Handling variability

ABSTRACT Pulse profile modelling is a relativistic ray-tracing technique that can be used to infer masses, radii, and geometric parameters of neutron stars. In a previous study, we looked at the performance of this technique when applied to thermonuclear burst oscillations from accreting neutron stars. That study showed that ignoring the variability associated with burst oscillation sources resulted in significant biases in the inferred mass and radius, particularly for the high count rates that are nominally required to obtain meaningful constraints. In this follow-on study, we show that the bias can be mitigated by slicing the bursts into shorter segments where variability can be neglected, and jointly fitting the segments. Using this approach, the systematic uncertainties on the mass and radius are brought within the range of the statistical uncertainty. With about 106 source counts, this yields uncertainties of approximately 10 per cent for both the mass and radius. However, this modelling strategy requires substantial computational resources. We also confirm that the posterior distributions of the mass and radius obtained from multiple bursts of the same source can be merged to produce outcomes comparable to that of a single burst with an equivalent total number of counts.

Astronomy & Astrophysics↗

Carbon, nitrogen, magnesium, silicon, and titanium isotopic compositions of single interstellar silicon carbide grains from the Murchison carbonaceous chondrite

Seven hundred and twenty SiC grains from the Murchison CM2 chondrite, ranging in size from 1 to 10 micrometers, were analyzed by ion microprobe mass spectrometry for their C-isotopic compositions. Subsets of the grains were also analyzed for N (450 grains), Si (183 grains), Mg (179 grains), and Ti (28 grains) isotopes. These results are compared with previous measurements on 41 larger SiC grains (up to 15 x 26 micrometers) from a different sample of Murchison analyzed by Virag et al. (1992) and Ireland, Zinner, & Amari (1991a). All grains of the present study are isotopically anomalous with C-12/C-13 ratios ranging from 0.022 to 28.4 x solar, N-14/N-15 ratios from 0.046 to 30 x solar, Si-29/Si-28 from 0.54 to 1.20 x solar, Si-30/Si-28 from 0.42 to 1.14 x solar, Ti-49/Ti-48 from 0.96 to 1.95 x solar, and Ti-50/Ti-48 from 0.94 to 1.39 x solar. Many grains have large Mg-26 excesses from the decay of Al-26 with inferred Al-26/Al-27 ratios ranging up to 0.61, or 12,200 x the ratio of 5 x 10(exp -5) inferred for the early solar system. Several groups can be distinguished among the SiC grains. Most of the grains have C-13 and N-14 excesses, and their Si isotopic compositions (mostly excesses in Si-29 and Si-30) plot close to a slope 1.34 line on a Delta Si-29/Si-28 versus Delta Si-30/Si-28 three-isotope plot. Grains with small C-12/C-13 ratios (less than 10) tend to have smaller or no N-14 excesses and high Al-26/Al-27 ratios (up to 0.01). Grains with C-12/C-13 greater than 150 fall into two groups: grains X have N-15 excesses and Si-29 and Si-30 deficits and the highest (0.1 to 0.6) Al-26/Al-27 ratios; grains Y have N-14 excesses and plot on a slope 0.35 line on a Si three-isotope plot. In addition, large SiC grains of the Virag et al. (1992) study fall into three-distinct clusters according to their C-, Si-, and Ti-isotopic compositions. The isotopic diversity of the grains and the clustering of their isotopic compositions imply distinct and multiple stellar sources. The C- and N-isotopic compositions of most grains are consistent with H-burning in the CNO cycle. These and s-process Kr, Xe, Ba, and Nd suggest asymptotic giant branch (AGB) or Wolf-Rayet stars as likely sources for the grains, but existing models of nucleosynthesis in these stellar sites fail to account in detail for all the observed isotopic compositions. Special problems are posed by grains with C-12/C-13 less than 10 and almost normal and heavy N-isotopic compositions. Also the Si- and Ti-isotopic compositions, with excesses in Si-29 and Si-30 relative to Si-28 and excesses in all Ti isotopes relative to Ti-48, do not precisely conform with the compositions predicted for slow neutron capture. Additional theoretical efforts are needed to achieve an understanding of the isotopic composition of the SiC grains and their stellar sources.

Hoppe, Peter↗

Ferroelectric Self-Polarization Controlled Magnetic Stratification and Magnetic Coupling in Ultrathin La 0.67 Sr 0.33 MnO 3 Films

Multiferroic oxide heterostructures consisting of ferromagnetic and ferroelectric components hold the promise for nonvolatile magnetic control via ferroelectric polarization, advantageous for the low-dissipation spintronics. Modern understanding of the magnetoelectric coupling in these systems involves structural, orbital, and magnetic reconstructions at interfaces. Previous works have long proposed polarization-dependent interfacial magnetic structures; however, direct evidence is still missing, which requires advanced characterization tools with near-atomic-scale spatial resolutions. Here, extensive polarized neutron reflectometry (PNR) studies have determined the magnetic depth profiles of PbZr 0.2 Ti 0.8 O 3 /La 0.67 Sr 0.33 MnO 3 (PZT/LSMO) bilayers with opposite self-polarizations. When the LSMO is 2–3 nm thick, the bilayers show two magnetic transitions on cooling. However, temperature-dependent magnetization is different below the lower-temperature transition for opposite polarizations. PNR finds that the LSMO splits into two magnetic sublayers, but the inter-sublayer magnetic couplings are of opposite signs for the two polarizations. Near-edge X-ray absorption spectroscopy further shows contrasts in both the Mn valences and the Mn–O bond anisotropy between the two polarizations. Finally, this work completes the puzzle for the magnetoelectric coupling model at the PZT/LSMO interface, showing a synergic interplay among multiple degrees of freedom toward emergent functionalities at complex oxide interfaces.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

VERA-Grizzly Ex-Core Calculations: Watts Bar Unit 1 Cycles 1-2

The critical structures that comprise light-water reactor (LWR) nuclear power plants are subjected to operating environments that can challenge their integrity. Structures in close proximity to the reactor core, such as the reactor pressure vessel (RPV) and the biological shield wall, are subjected to high levels of radiation emanating from the core, as well as elevated temperatures. As the US fleet of operating LWRs ages, the effects of these operating environments on the integrity of these structures must be considered to ensure their continued safe operation. Extending the lifetime of commercial reactors and maintaining the aging reactor fleet require accurate prediction of the exposure of ex-core components to neutron and photon radiation. In particular, concrete degradation studies must be performed to evaluate the safety and long-term operation of reactors with lifetime extensions. The concrete reactor bioshield is important for providing radiological protection during operation and must last for the entire lifetime of the reactor. Recent interest in lifetime extensions furthers the need to accurately simulate concrete material degradation in the reactor bioshield. As a result of this need, the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program has funded this study to couple its tools, Virtual Environment for Reactor Applications (VERA) and Grizzly. VERA allows users to set up models to calculate time-dependent and fully coupled solutions (with thermal feedback) for ex-core quantities of interest such as vessel and coupon fluence and detector responses for multiple statepoints and cycles. Grizzly is a finite-element application based on the Multiphysics Object Oriented Simulation Environment (MOOSE) framework that is used to enable aging materials calculations. This report highlights the work performed to calculate the fluence in the vessel and concrete for Watts Bar Nuclear Plant Unit 1 (WBN1) Cycles 1 and 2. The fluences obtained from VERA were successfully transferred to Grizzly using a Python script. Four simulations were run with Grizzly: (1) the Mazars model with the initial Young’s modulus being the instantaneous modulus, (2) the Mazars model with the initial Young’s modulus being the delayed modulus, (3) the Mazars model with the initial Young’s modulus being the delayed modulus with the addition of the effects of micro-damage caused by irradiation, and (4) the Mazars model with the initial Young’s modulus being the instantaneous modulus, and with the addition of micro-damage and creep. Details regarding the methods used to obtain the fluence and the statistical errors associated with the VERA Monte Carlo Shift calculations are discussed in greater detail in this report. The results obtained from the four Grizzly models are also presented in this report.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Ordering Unstructured Meshes for Sparse Matrix Computations on Leading Parallel Systems

The ability of computers to solve hitherto intractable problems and simulate complex processes using mathematical models makes them an indispensable part of modern science and engineering. Computer simulations of large-scale realistic applications usually require solving a set of non-linear partial differential equations (PDES) over a finite region. For example, one thrust area in the DOE Grand Challenge projects is to design future accelerators such as the SpaHation Neutron Source (SNS). Our colleagues at SLAC need to model complex RFQ cavities with large aspect ratios. Unstructured grids are currently used to resolve the small features in a large computational domain; dynamic mesh adaptation will be added in the future for additional efficiency. The PDEs for electromagnetics are discretized by the FEM method, which leads to a generalized eigenvalue problem Kx = AMx, where K and M are the stiffness and mass matrices, and are very sparse. In a typical cavity model, the number of degrees of freedom is about one million. For such large eigenproblems, direct solution techniques quickly reach the memory limits. Instead, the most widely-used methods are Krylov subspace methods, such as Lanczos or Jacobi-Davidson. In all the Krylov-based algorithms, sparse matrix-vector multiplication (SPMV) must be performed repeatedly. Therefore, the efficiency of SPMV usually determines the eigensolver speed. SPMV is also one of the most heavily used kernels in large-scale numerical simulations.

Oliker, Leonid↗

CFD Evaluation of Pressure Change Along Coolant Passages in Sodium-Cooled Fast Reactor with Nek5000

To support the design efforts of advanced sodium-cooled fast reactors (SFRs), a series of computational fluid dynamics (CFD) simulations are performed to investigate the pressure change along various flow passages in the proposed SFR system. The simulations are carried out with the state-of-the-art spectral element flow solver, Nek5000. Two specific case studies are presented in this paper: the flow exiting the axial neutron reflector channels and the flow entering the fuel pin bundle. Due to the high Reynolds numbers expected, a Reynolds-averaged Navier-Stokes (RANS) approach is necessary to model the turbulence. A newly developed regularized k-omega RANS model is adopted in the related CFD calculations. The first case study explores the effect of Reynolds number on the pressure change when flow exits the reflector channels. The pressure change in this case has two major contributors: the change due to wall friction and the Bernoulli effect. It is noted that the nondimensional pressure loss follows a log-linear trend up to Re = 10 5 , and then the trend is flattened. In the second case study, the advanced NekNek coupling capability is tested where an integral domain can be divided into multiple subdomains with coupling interfaces, which would greatly ease the meshing process of complex engineering geometries and potentially save computational resources. The preliminary results obtained so far confirm the consistency between the NekNek results and those produced by regular Nek5000 simulation. The presented work demonstrates the readiness and flexibility of the related CFD techniques, which is part of the broader effort to leverage cutting-edge CFD to inform the advanced nuclear reactor designs.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗