Search NASA⌕ Search

SEARCH · Search NASA

Results for “coupling analysis”

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 379 records · Page 21

Application of sensitivity analysis in DYMOND/Dakota to fuel cycle transition scenarios

The ability to perform sensitivity analysis has been enabled for the nuclear fuel cycle simulator DYMOND through its coupling with the design and analysis toolkit Dakota. To test and demonstrate these new capabilities, a transition scenario and multi-parameter study were devised. The transition scenario represents a partial transition from the US nuclear fleet to a closed fuel cycle with small modular LWRs and fast reactors fueled by reprocessed used nuclear fuel. Four uncertain parameters in this transition were studied – start date of reprocessing, total reprocessing capacity, the nuclear energy demand growth, and the rate at which the fast reactors are deployed – with respect to their impact on four response metrics. The responses – total natural uranium consumed, maximum annual enrichment capacity required, total disposed mass, and total cost of the nuclear fuel cycle – were chosen based on measures known to be of interest in transition scenarios and to be significantly impacted by the varying parameters. Furthermore, analysis of this study was performed both from the direct sampling and through surrogate models developed in Dakota to calculate the global sensitivity measures Sobol’ indices. This example application of this new capability showed that the most consequential parameter to most metrics was the share of new build capacity that is fast reactors. However, for the cost metric, the scaling factor of the energy demand growth was significant and had synergistic behavior with the fast reactor new build share.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

An assessment of spall failure modes in laser powder bed fusion fabricated stainless steel 316L with low-volume intentional porosity

Here, this paper reports on spall failure and damage modes in Laser Powder Bed Fusion fabricated Stainless Steel 316L (SS316L) with intentional levels of low-volume (1–5 vol. %) porosity and pore sizes of 200, 350, and 500 μm. The fabricated specimens were subjected to uniaxial-strain plate-impact loading at ~4.5 GPa, to initiate incipient spall failure. Analysis of velocimetry profiles measured using multi-probe photon-Doppler velocimetry coupled with post-mortem analysis of soft-recovered samples reveals local suppression of spall failure (termed as spall-dominated) as a function of porosity, as the failure mechanism transitions from spall-centered tensile stress dominated to a pore-centered microstructure-dominated damage mode involving void/crack nucleation and growth at pre-existing pores. The critical porosity level where the suppression of spall failure is first observed, as well as the spall location, is dependent on both the volume fraction and the size of the initially fabricated pores. In samples of 500 μm pore size, the suppression of spall failure is observed with as little as 1 vol. % porosity, while samples with smaller pores (200 μm) still experience spall-centered tensile stress dominated failure with higher levels (5 vol. %) of porosity. In the case of pore-centered microstructure-dominated failure, spall damage can occur but the spall plane is shifted toward the rear free surface, or more generally in areas further away from the region with pores. Highly heterogeneous deformation twinning, shear banding, grain rotation, and cracking are observed in the vicinity of pre-existing pores and expected spall failure sites.

36 MATERIALS SCIENCE↗

New developments and verification of fusion blanket simulation capabilities in the MOOSE framework

Multiphysics modeling capabilities have a crucial role to play in the accelerated deployment of fusion energy. To that end, we developed new multiphysics fusion blanket simulation capabilities in the Multiphysics Object-Oriented Simulation Environment (MOOSE). Firstly, we expanded on the existing capabilities of the previously published work, by coupling 3D tritium transport modeling capabilities using the Tritium Migration Analysis Program, version 8 (TMAP8) to an existing tool including thermal hydraulics, fully three-dimensional (3D) heat transfer, and loosely coupled neutronics analysis. Secondly, we performed a thorough verification of the new capabilities and increased testing code coverage to meet MOOSE’s software quality standards. The MOOSE framework follows a strict software quality assurance plan to be Nuclear Quality Assurance, Level 1 compliant. The new multiphysics fusion blanket simulation capabilities are now held to the same standard. Thirdly, to demonstrate MOOSE’s new fusion blanket modeling capabilities, we performed a fully integrated, multiphysics simulation of a 3D solid ceramic breeder blanket design. This proof-of-concept simulation provides the temperature and tritium distribution across the blanket. In conclusion, the combined efforts towards software quality and the development of multiphysics coupling capabilities provide an effective and reliable framework for modeling solid ceramic fusion blankets using MOOSE.

modeling and simulation↗

Arctic Tipping Points Triggering Global Change (LDRD Final Report)

The Arctic is warming and feedbacks in the coupled Earth system may be driving the Arctic to tipping events that could have critical downstream impacts for the rest of the globe. In this project we have focused on analyzing sea ice variability and loss in the coupled Earth system Summer sea ice loss is happening rapidly and although the loss may be smooth and reversible, it has significant consequences for other Arctic systems as well as geopolitical and economic implications. Accurate seasonal predictions of sea ice minimum extent and long-term estimates of timing for a seasonally ice-free Arctic depend on a better understanding of the factors influencing sea ice dynamics and variation in this strongly coupled system. Under this project we have investigated the most influential factors in accurate predictions of September Arctic sea ice extent using machine learning models trained separately on observational data and on simulation data from five E3SM historical ensembles. Monthly averaged data from June, July, and August for a selection of ice, ocean, and atmosphere variables were used to train a random forest regression model. Gini importance measures were computed for each input feature with the testing data. We found that sea ice volume is most important earlier in the season (June) and sea ice extent became a more important predictor closer to September. Results from this study provide insight into how feature importance changes with forecast length and illustrates differences between observational data and simulated Earth system data. We have additionally performed a global sensitivity analysis (GSA) using a fully coupled ultra- low resolution configuration E3SM. To our knowledge, this is the first global sensitivity analysis involving the fully-coupled E3SM Earth system model. We have found that parameter variations show significant impact on the Arctic climate state and atmospheric parameters related to cloud parameterizations are the most significant. We also find significant interactions between parameters from different components of E3SM. The results of this study provide invaluable insight into the relative importance of various parameters from the sea ice, atmosphere and ocean components of the E3SM (including cross-component parameter interactions) on various Arctic-focused quantities of interest (QOIs).

54 ENVIRONMENTAL SCIENCES↗

Cost-constrained adaptive simulations of transient spray combustion in a gas turbine combustor

Predictive high-fidelity simulations of turbulent spray combustion must capture the combined effects of complex chemistry, multiphase evaporating flow and spray-flame interactions to achieve physical accuracy. Finite-rate chemistry (FRC) combined with a realistic chemical mechanism is a combustion model well-suited for this purpose, but has a high computational cost due to the large number and stiffness of transported chemical species. In contrast, flamelet-based models achieve lower cost by transporting a small number of quantities of reduced stiffness, but assumptions regarding local flame topology, boundary conditions and inter-phase coupling limit their physical accuracy. Recently, the Pareto-efficient combustion (PEC) framework was developed to dynamically assign combustion models based on local cost and accuracy metrics in gas-phase reacting flows. In this work, we extend this PEC framework to spray combustion through the rigorous analysis of the multiphase coupling terms in the governing equations. The derivation shows that spray evaporation causes errors in the prediction of species mass fractions for flamelet-based models due to the sensitivity of the local thermo-chemical state to changes in composition caused by fuel vaporization across combustion regimes present in practical spray combustion devices. Sub-model assignment is formulated as a multiple-choice knapsack problem, where computational cost is directly controlled through the fraction of the domain assigned to the FRC sub-model. The extended PEC formulation is applied to the simulation of a realistic rich-quench-lean gas turbine combustor at steady-state conditions, as well as transient operation resulting in lean blow-out (LBO). Analysis of transient simulations during LBO demonstrates the extended PEC formulation’s capacity to dynamically adapt to changing conditions within the combustor. Transient combustor dynamics are shown to approach convergence with limited increases in computational cost, while retaining substantial computational cost reduction compared to monolithic FRC simulations. Through PEC simulations with increasing fractions of the domain assigned to FRC, monolithic flamelet simulations are shown to over-predict flame stability during LBO. Finally, the extended PEC formulation is thus shown to overcome deficiencies of monolithic models by controlling modeling error for multiphase combustion modeling.

33 ADVANCED PROPULSION SYSTEMS↗

Enhanced sensitivity to trace 238 U impurity of sapphire via coincidence neutron activation analysis

Sapphire has mechanical and electrical properties that are advantageous for the construction of internal components of radiation detectors such as time projection chambers and bolometers. However, it has proved difficult to assess its 232 Th and 238 U content down to the picogram per gram level. Here, this work reports an experimental verification of a computational study that demonstrates 𝛾⁢𝛾 coincidence counting, coupled with neutron activation analysis (NAA), can reach ppt sensitivities. Combining results from 𝛾⁢𝛾 coincidence counting with those of earlier single-𝛾 counting based NAA shows that a sample of Saint Gobain sapphire has 232 Th and 238 U concentrations of <0.26 ppt and <2.3 ppt, respectively; the best constraints on the radiopurity of sapphire.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Climatology of Cloud‐Land‐Surface Coupling Across Different ARM Sites

Land-atmosphere interactions play a critical role in the evolution and formation of low-level clouds. The different states of coupling between low-level clouds and the surface are uncertain, primarily over continental regions, where complex thermodynamics complicates their investigation. This study uses observations from the Atmospheric Radiation Measurement User Facility to explore cloud-surface coupling and perform a climatological analysis of this interaction in five countries across three continents. The results reveal consistent coupling thresholds and average percentages across the five sites, with coupled clouds accounting for 66% of the cases and decoupled clouds for 34%. Thermodynamic and dynamic evaluations show distinct differences between coupled and decoupled clouds. Coupled clouds are characterized by humid environments, in which vertical motions connect the surface and lower atmosphere to the cloud base, conditions that favor the formation of boundary layer clouds. Decoupled clouds prefer to occur in a drier and colder environment with vertical motions inside the boundary layer being detached from the cloud base, under which boundary layer clouds are hard to form. Coupled clouds peak during warmer hours and seasons, and vice versa for decoupled clouds. This study underscores the complexity of cloud-land-surface interactions and paves the way for further investigations into cloud formation and evolution under different atmospheric environments.

54 ENVIRONMENTAL SCIENCES↗

Exhaled breath condensate profiles of U.S. Navy divers following prolonged hyperbaric oxygen (HBO) and nitrogen-oxygen (Nitrox) chamber exposures

Prolonged exposure to hyperbaric hyperoxia can lead to pulmonary oxygen toxicity (PO 2 tox). PO 2 tox is a mission limiting factor for special operations forces divers using closed-circuit rebreathing apparatus and a potential side effect for patients undergoing hyperbaric oxygen (HBO) treatment. In this study, we aim to determine if there is a specific breath profile of compounds in exhaled breath condensate (EBC) that is indicative of the early stages of pulmonary hyperoxic stress/PO 2 tox. Using a double-blind, randomized 'sham' controlled, cross-over design 14 U.S. Navy trained diver volunteers breathed two different gas mixtures at an ambient pressure of 2 ATA (33 fsw, 10 msw) for 6.5 h. One test gas consisted of 100% O 2 (HBO) and the other was a gas mixture containing 30.6% O 2 with the balance N 2 (Nitrox). The high O 2 stress dive (HBO) and low O 2 stress dive (Nitrox) were separated by at least seven days and were conducted dry and at rest inside a hyperbaric chamber. EBC samples were taken immediately before and after each dive and subsequently underwent a targeted and untargeted metabolomics analysis using liquid chromatography coupled to mass spectrometry (LC-MS). Following the HBO dive, 10 out of 14 subjects reported symptoms of the early stages of PO 2 tox and one subject terminated the dive early due to severe symptoms of PO 2 tox. No symptoms of PO 2 tox were reported following the nitrox dive. A partial least-squares discriminant analysis of the normalized (relative to pre-dive) untargeted data gave good classification abilities between the HBO and nitrox EBC with an AUC of 0.99 (±2%) and sensitivity and specificity of 0.93 (±10%) and 0.94 (±10%), respectively. Furthermore, the resulting classifications identified specific biomarkers that included human metabolites and lipids and their derivatives from different metabolic pathways that may explain metabolomic changes resulting from prolonged HBO exposure.

59 BASIC BIOLOGICAL SCIENCES↗

Optical Microscopy-Guided Laser Ablation Electrospray Ionization Ion Mobility Mass Spectrometry: Ambient Single Cell Metabolomics with Increased Confidence in Molecular Identification

Single cell analysis is a field of increasing interest as new tools are continually being developed to understand intercellular differences within large cell populations. Laser-ablation electrospray ionization mass spectrometry (LAESI-MS) is an emerging technique for single cell metabolomics. Over the years, it has been validated that this ionization technique is advantageous for probing the molecular content of individual cells in situ. Here, we report the integration of a microscope into the optical train of the LAESI source to allow for visually informed ambient in situ single cell analysis. Additionally, we have coupled this ‘LAESI microscope’ to a drift-tube ion mobility mass spectrometer to enable separation of isobaric species and allow for the determination of ion collision cross sections in conjunction with accurate mass measurements. This combined information helps provide higher confidence for structural assignment of molecules ablated from single cells. Here, we show that this system enables the analysis of the metabolite content of Allium cepa epidermal cells with high confidence structural identification together with their spatial locations within a tissue.

59 BASIC BIOLOGICAL SCIENCES↗

Geometric phase in coupled cluster theory

It has been well-established that the topography around conical intersections between excited electronic states is incorrectly described by coupled cluster and many other single reference theories (the intersections are “defective”). Despite this, we show both analytically and numerically that the geometric phase effect (GPE) is correctly reproduced upon traversing a path around a defective excited-state conical intersection (CI) in coupled cluster theory. The theoretical analysis is carried out by using a non-Hermitian generalization of the linear vibronic coupling approach. Interestingly, the approach qualitatively explains the characteristic (incorrect) shape of the defective CIs and CI seams. Moreover, the validity of the approach and the presence of the GPE indicate that defective CIs are local (and not global) artifacts. This implies that a sufficiently accurate coupled cluster method could predict nuclear dynamics, including geometric phase effects, as long as the nuclear wavepacket never gets too close to the conical intersections.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

ANTARES: a new IRSN computational chain for safety analysis

Since some years, the IRSN has developed a new project to simulate accurately the interaction between neutronic and thermal hydraulic phenomena. This project called ANTARES (Advanced Neutronics and Thermal-hydraulic for the Analysis of Reactor Safety) couples the thermal hydraulic code CATHARE-3 with the 3D neutronic nodal code PARCS. The purpose of this project is to take the advantage of the complementary of these 2 codes, by the knowledge of the 3D core power distribution in CATHARE-3 (computed by PARCS) and improving the thermal hydraulic boundaries conditions in PARCS. In this way, the IRSN diversifies and enlarges its available tools to perform safety analysis with improved accuracy. The current technique usually adopted in France for the safety demonstrations for Reactivity Insertion Accident (RIA) is the so-called 'conservative' approach, which consists of reducing all the feedback (Doppler and moderator effects) and in modifying some physical quantities in such a way to increase a power peak. For this reason, these facilities ('penalties') have been implemented in ANTARES. In this paper we will give two examples of accidental transients that can be simulated with ANTARES: a Main Steam Line Break (MSLB) and an inadvertent boron dilution event. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Energy transport analysis of NSTX plasmas with the TGLF turbulent and NEO neoclassical transport models

This work presents a study of plasma transport at low aspect ratio on the National Spherical Torus Experiment tokamak, where the turbulent and neoclassical energy fluxes calculated by the quasilinear Trapped Gyro Landau Fluid (TGLF) model and the multi species drift-kinetic Neoclassical solver (NEO) are validated against experimental data. The turbulent energy transport of two plasma discharges, one in the L-mode confinement regime and another in the H-mode regime, is dominated by electrostatic drift-wave instabilities, while the ion heat transport has a significant neoclassical contribution. The data analysis workflow is described in detail to understand how the variations of mapping and fitting of experimental data affect the power balance solution and subsequent flux-matching plasma profile predictions with the TGYRO solver. On average, the predicted plasma profiles are consistent with experimental data. However, the solutions are sensitive to various input parameters, including boundary conditions, and the electron-ion coupling. Linear gyrokinetic stability analysis demonstrates close agreement of the real frequencies of unstable modes between TGLF and CGYRO gyrokinetic simulations, but higher growth rates are predicted by TGLF, especially for the H-mode case. Estimates of the low-k modes' contributions to the total flux are consistent with linear stability analysis and the E × B suppression of turbulence in TGLF simulations with the SAT1 saturation model, while the SAT2 saturation model over-predicts the low-k modes' contribution in the H-mode case. Moreover, the results with SAT1 model are consistent with power balance analysis, which indicates only neoclassical ion energy fluxes inside ρ < 0.4 in the L-mode case and $\rho \unicode{x2A7D} 0.7$ in the H-mode case. The presence of multi-scale turbulence and ion-scale driven zonal flow mixing effects are also observed in TGLF scans of the electron turbulent heat flux over a range of temperature gradients and the electron-ion temperature ratio, which could explain the strong model sensitivity to variations of input parameters.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Two-step neutronics calculations with Shift and Griffin for advanced reactor systems

This research develops the initial coupling of the Shift Monte Carlo (MC) code and the Griffin reactor physics code for reactor analysis of non–light-water reactor systems. The novelty of this work is twofold. It is the first application of Shift to produce the multigroup cross sections needed for Griffin as applied to a non–light-water reactor system; and, the first investigation and analysis of characteristics of the Empire microreactor benchmark that should be considered for steady state and transient reactor physics calculations. This application uses the previously developed two-step neutronics analysis workflow to demonstrate this initial coupling. Here, we outline the two-step neutronics analysis workflow in which the Shift MC code is used to generate the multigroup cross sections and fluxes needed by the Griffin deterministic solver. Details on how these multigroup cross sections are generated using MC tallies are given, as well as the practicalities and limitations of the two-step neutronics workflow. The Empire microreactor benchmark was used to investigate and validate this coupling. Results using this benchmark show good agreement between Griffin calculations using Serpent-generated cross sections and Shift-generated cross sections. Analysis of the characteristics of this Empire benchmark show larger eigenvalue differences between heterogeneous and pin–homogenized solutions compared to those of traditional light-water reactor (LWR) designs, thus requiring super homogenization factor corrections for accurate eigenvalue and power distribution predictions.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Magnetic structure and spin dynamics of the quasi-two-dimensional antiferromagnet Zn-doped copper pyrovanadate

We report magnetic properties of the antiferromagnet Zn x Cu 2-x V 2 O 7 (ZnCVO) with x ≈ 0.06 have been thoroughly investigated on powder and single-crystal samples. The crystal structure determination using powder x-ray and neutron diffraction confirms that our ZnCVO samples are isostructural with β-Cu 2 V 2 O 7 (β-CVO) with a small deviation in the lattice parameters. Macroscopic magnetic property measurements also confirm the similarity between the two compounds. The Cu 2+ spins were found to align along the crystallographic c axis, antiparallel to their nearest neighbors connected by the leading exchange interaction J 1 . Spin dynamics reveals a typical symmetric spin-wave dispersion with strong interactions in the bc plane and weak interplane coupling. The exchange interaction analysis indicates that the spin network of ZnCVO is topologically consistent with the previous DFT prediction but the values of leading exchange interactions are contradictory. Furthermore, rather than the predicted 2D honeycomb structure, the spin network in ZnCVO could be better described by the anisotropic 2D spin network composed of J 1 , J 5 , and J 6 interactions, four bonds per one spin site, coupled by weak interplane interactions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Geothermal Play-Fairway Analysis of Washington State Prospects: Final Report

The Washington State Geothermal Play-Fairway Analysis overcomes the exploration challenges posed by dense vegetation, glacial deposits, and extreme precipitation. The geothermal play-fairways we target are locations where heat, permeability, and saturated porosity are present in sufficient volume to provide adequate heat exchange at depths accessible by modern drilling technology. The three study areas lie along the Cascade Range magmatic arc and are near Mount Baker, Mount St. Helens, and the Wind River Valley. The seven-year project is divided into three phases. In Phase 1 we build on a previous statewide assessment of geothermal resources and develop an initial modeling approach. The results are a series of favorability, uncertainty, and risk maps for three targeted study areas. Based on these initial results, we collect new geologic and geophysical data to further refine our modeling and reduce exploration uncertainty in Phase 2. We improve the modeling method to handle the new data and update the favorability, uncertainty, and risk maps. We also update the conceptual geothermal resource models. In Phase 3 we validate our modeling approach by drilling two temperature-gradient holes and collecting and analyzing core, image logs, and new geochemistry. Our modeling approach improves on an earlier statewide method through a more-rigorous and detailed assessment of heat and permeability. Permeability potential is assessed through geomechanical modeling of the deformation that can generate and maintain reservoir porosity and permeability. Metrics to inform heat potential include temperature-gradient wells, which are sparse in Washington; proximity of Quaternary volcanic vents and young intrusive rock; spring temperature; and reservoir temperature inferred from geothermometry. We weight the individual components using an expert-guided approach known as the Analytical Hierarchy Process. During Phase 2 we also develop a fluid-filled fracture model, and an infrastructure model that helps to delineate areas which are more favorable for geothermal development based on proximity to transmission lines, elevation, land ownership and use restrictions, and availability of process water. New geologic and geophysical data is collected during Phase 2 in each of our three main study areas. At Mount Baker and north of Mount St. Helens we conduct 1:24,000-scale geologic mapping and lidar analysis to better constrain the location and character of surface faults; detailed mapping in the Wind River Valley was completed just prior to the start of this project. Ages of intrusive rocks are determined with 40 Ar/ 39 Ar geochronology, though all of our samples are Miocene or older. We collect ground based gravity observations (a total of 1,580 new stations) in all of our study areas and ground-based magnetic lines (a total of 93 km) at Mount Baker. These data are combined with existing gravity and aeromagnetic data and used to constrain fault locations and geometry. Two to three cross sections are constructed at each study area using the mapped surface geology and forward-modeling of the gravity and magnetic data; these cross sections form the basis for our updated conceptual models. We collect magnetotelluric surveys at Mount Baker and Mount St. Helens and these data are inverted to form a resistivity model from the surface to about 10 km depth; each model shows conductive zones that can be interpreted as upwelling geothermal fluids. At Mount St. Helens we deploy a passive seismic array and use the newly detected events to refine the location of the Saint Helens seismic zone. We also employ ambient-noise tomography to develop a detailed seismic-velocity model for the study area and use this model to help constrain our cross sections and conceptual model. Based on the new data collected during Phase 2—and our updated models—we develop a campaign of temperature-gradient holes and core analysis to validate our modeling in Phase 3. Drill hole MB76-31 is located near Little Park Creek, 11 km west-southwest of the summit of Mount Baker, and is 1,471 ft deep. About 410 ft of core from the lower portion of the hole—and image logs from ~175 ft below ground surface to the bottom—are collected and analyzed. Water samples are collected and processed for geothermometry. Drill hole MSH17-24 is located along upper Schultz Creek, 16 km north-northeast of Mount St. Helens and has core from 470 ft to the bottom at 1,053 ft. We did not collect image logs due to borehole stability concerns, but water samples are collected and analyzed for geothermometry. Repeat temperature-gradient measurements are made at both sites and thermal conductivity is measured from core samples. At MB76-31, the equilibrated temperature gradient of 64°C/km and calculated heat flow of 141–159 mW/m 2 is more than twice the regional average. Detailed mapping and analysis of the core, coupled with correlation to the image logs, indicates a history of permeability generation consistent with our predictions of high permeability. Because the site has high favorability in the Phase 2 model, we consider the results a positive validation of the modeling. At site MSH17-24, the equilibrated temperature gradient of ~15°C/km and calculated heat flow of 41–43 mW/m 2 are similar to regional. Geochemical analysis of the water samples indicates a meteoric source without any geothermal component. Detailed outcrop-based mapping of fault exposures near the drill site and analysis of image logs from nearby boreholes indicates a history of permeability generation consistent with our predictions. Because the site has low favorability in the Phase 2 model, we consider the results a positive validation of the modeling. Together, the two sites provide a reasonably positive validation of the Phase 2 modeling and should encourage future use of this modeling approach.

15 GEOTHERMAL ENERGY↗

NRC Multiphysics Analysis Capability Deployment (FY2021 - Part 2)

This report details the progress and activities of Idaho National Laboratory (INL) on the Nuclear Regulatory Commission (NRC) project “Development and Modeling Support for Advanced Non-Light Water Reactors.” Task 4b was completed for this report. INL developed a sample problem showing how to use Serpent 2 to calculate macroscopic cross sections for use in Griffin for a typical sodium fast reactor Unprotected Loss of Flow (ULOF) transient. The complicating factor is the relatively large axial gradient of the coolant (and hence fuel) temperature and its evolution during the transient. A 3D Griffin model of the Advanced Burner Test Reactor (ABTR) is coupled to the System Analysis Module (SAM) to perform the ULOF transient. The model includes various explicit feedback mechanisms for fast reactors, including Doppler, radial expansion from the displacement of the support plate, and axial expansion from the displacement of the fuel pins. The results for the ULOF transient are consistent with published values. We observed small effects in the transient results that arise from the superhomogenization equivalence correction of uniform and nonuniform temperature data sets, but they diminish as more dominant reactivity mechanisms are added to the model. Potential improvements to the neutronics model include: adding the differential control rod and sodium feedback mechanisms, implementing cross section re-homogenization based on the partial volume of the various materials, and enhancing the kinetics parameters. Furthermore, the fluids model only includes four subchannels to represent the entire active core due to performance limitations in SAM when attempting to model one subchannel per assembly.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Scalar and tensor resonances in $J/\psi $ radiative decays

We perform a systematic analysis of the J/ψ → γπ 0 π 0 and → γ$K$$^{0}_{S}$$K$$^{0}_{S}$ partial waves measured by BESIII. We use a large set of amplitude parametrizations to reduce the model bias. We determine the physical properties of seven scalar and tensor resonances in the 1–2.5GeV mass range. These include the well known f 0 (1500) and f 0 (1710), that are considered to be the primary glueball candidates. The hierarchy of resonance couplings determined from this analysis favors the latter as the one with the largest glueball component.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Spatially Resolved Near Field Spectroscopy of Vibrational Polaritons at the Small N Limit

Vibrational polaritons, which have been primarily studied in Fabry–Pérot cavities with a large number of molecules (N ~ 10 6 –10 10 ) coupled to the resonator mode, exhibit various experimentally observed effects on chemical reactions. However, the exact mechanism is elusively understood from the theoretical side, as the large number of molecules involved in an experimental strong coupling condition cannot be represented completely in simulations. This discrepancy between theory and experiment arises from computational descriptions of polariton systems typically being limited to only a few molecules, thus failing to represent the experimental conditions adequately. To address this mismatch, we used surface phonon polariton (SPhP) resonators as an alternative platform for vibrational strong coupling. SPhPs exhibit strong electromagnetic confinement on the surface and thus allow for coupling to a small number of molecules. As a result, this platform can enhance nonlinearity and slow down relaxation to the dark modes. In this study, we fabricated a pillar-shaped quartz resonator and then coated it with a thin layer of cobalt phthalocyanine (CoPc). By employing scattering-type scanning near-field optical microscopy (s-SNOM), we spatially investigated the dependency of vibrational strong coupling on the spatially varying electromagnetic field strength and demonstrated strong coupling with 38,000 molecules only–reaching to the small N limit. Through s-SNOM analysis, we found that strong coupling was observed primarily on the edge of the quartz pillar and the apex of the s-SNOM tip, where the maximum field enhancement occurs. In contrast, a weak resonance signal and lack of coupling were observed closer to the center of the pillar. This work demonstrates the importance of spatially resolved polariton systems in nanophotonic platforms and lays a foundation to explore polariton chemistry and chemical dynamics at the small N limit–one step closer to reconcile with high-level quantum calculations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗