SEARCH · Engineering Papers

Results for “coupling”

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.

SourceResearch connectionFrom yearThrough year

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

At least 865 records · Page 48

All sources NASA NTRS DOE OSTI

DOE OSTI2017

Novel phases in strongly coupled four-fermion theories

Here, we study a lattice model comprising four massless reduced staggered fermions in four dimensions coupled through an SU(4)-invariant four-fermion interaction. We present both theoretical arguments and numerical evidence that no bilinear fermion condensates are present for any value of the four-fermi coupling, in contrast to earlier studies of Higgs{Yukawa models with different exact lattice symmetries. At strong coupling we observe the formation of a four-fermion condensate and a mass gap in spite of the absence of bilinear condensates. Unlike those previously studied systems we do not find a ferromagnetic phase separating this strong-coupling phase from the massless weak-coupling phase. Instead we observe long-range correlations in a narrow region of the coupling, still with vanishing bilinear condensates. While our numerical results come from relatively small lattice volumes that call for caution in drawing conclusions, if this novel phase structure is verified by future investigations employing larger volumes it may offer the possibility for new continuum limits for strongly interacting fermions in four dimensions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2018

Preheating after multifield inflation with nonminimal couplings. III. Dynamical spacetime results

This paper concludes our semianalytic study of preheating in inflationary models comprised of multiple scalar fields coupled nonminimally to gravity. Using the covariant framework of paper I in this series, we extend here the rigid-spacetime results of paper II by considering both the expansion of the Universe during preheating, as well as the effect of the coupled metric perturbations on particle production. The adiabatic and isocurvature perturbations are governed by different effective masses that scale differently with the nonminimal couplings and evolve differently in time. The effective mass for the adiabatic modes is dominated by contributions from the coupled metric perturbations immediately after inflation. The metric perturbations contribute an oscillating tachyonic term that enhances an early period of significant particle production for the adiabatic modes, which ceases on a time scale governed by the nonminimal couplings ${{\xi}}_{I}$. The effective mass of the isocurvature perturbations, on the other hand, is dominated by contributions from the fields' potential and from the curvature of the field-space manifold (in the Einstein frame), the balance between which shifts on a time scale governed by ${{\xi}}_{I}$. As in papers I and II, we identify distinct behavior depending on whether the nonminimal couplings are small [${{\xi}}_{I}{\lesssim}\mathcal{O}(1)$], intermediate [${{\xi}}_{I}{\sim}\mathcal{O}(1{-}10)$], or large (${{\xi}}_{I}{\ge}100$).

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2018

Quantum Field Theory of Nematic Transitions in Spin-Orbit-Coupled Spin-1 Polar Bosons

We theoretically study an ultracold gas of spin-1 polar bosons in a one-dimensional continuum, which are subject to linear and quadratic Zeeman fields and a Raman induced spin orbit coupling. Concentrating on the regime in which the background fields can be treated perturbatively, we analytically solve the model in its low-energy sector; i.e., we characterize the relevant phases and the quantum phase transitions between them. Depending on the sign of the effective quadratic Zeeman field ε, two superfluid phases with distinct nematic order appear. In addition, we uncover a spin-disordered superfluid phase at strong coupling. We employ a combination of renormalization group calculations and duality transformations to access the nature of the phase transitions. At ε=0, a line of spin-charge separated pairs of Luttinger liquids divides the two nematic phases, and the transition to the spin-disordered state at strong coupling is of the Berezinskii-Kosterlitz-Thouless type. In contrast, at ε≠0, the quantum critical theory separating nematic and strong coupling spin-disordered phases contains a Luttinger liquid in the charge sector that is coupled to a Majorana fermion in the spin sector (i.e., the critical theory at finite ε maps to a quantum critical Ising model that is coupled to the charge Luttinger liquid). Because of an emergent Lorentz symmetry, both have the same logarithmically diverging velocity. We discuss the experimental signatures of our findings that are relevant to ongoing experiments in ultracold atomic gases of 23 Na.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2024

Comparing Exciton-Polariton Couplings in Two-Dimensional MoS2 Cavities

Exciton-polaritons (EPs) are an emerging approach to achieve strong light-matter coupling and robust quantum phenomena, with application to areas such as photonic neuromorphic computing for instance reservoir computing. 2D transition metal dichalcogenides (TMDs) are promising candidates to host EPs due to their potential for strong coupling as evidenced by large Rabi splitting energies, which in the case of few-layer MoS2 has been demonstrated to reach 293 meV when coupled to the C exciton, placing it in the ultrastrong coupling regime, at room temperature. Here we examine how the Rabi splitting of CVD-grown monolayer MoS2 in a cavity can be optimized by tuning two variables: cavity geometry (thickness of silver and dielectric spacers) and the targeted exciton for coupling (A at about 1.9 eV, B at about 2.06 eV, or C at about 2.9 eV). With the cavity geometry, there is a tradeoff between optical field confinement and transmission. With the targeted exciton, there may be a tradeoff between Rabi splitting and lifetime. Transfer matrix simulations are used to inform how thick the cavity layers need to be, and at what angle measurements should be performed, to achieve Rabi splitting at a certain exciton energy. Experimentally, the Rabi splitting is measured by placing the cavity at a specified angle and measuring the optical reflectance and transmittance. By fitting the data, we find Rabi splittings of 124 meV, 194 meV, and 288 meV for the A, B, and C excitons, respectively, which places them in the strong, and on the edge of the ultrastrong, coupling regime.

exciton-polaritons↗

DOE OSTI2022

Evaluation of coupled wind-wave model simulations of offshore winds in the Mid-Atlantic Bight using lidar-equipped buoys.

From 2014 to 2017, two Department of Energy buoys equipped with Doppler lidar were deployed off the U.S. East Coast to provide long term measurements of hub-height wind speed in the marine environment. In this study, we performed simulations of selected cases from the deployment using a 5-km configuration of the Weather Research and Forecasting (WRF) model, to see if simulated hub height speeds could produce closer agreement with the observations than existing reanalysis products. For each case we performed two additional simulations: one in which marine surface roughness height was one-way coupled to forecast wave parameters from a standalone WaveWatch III (WW3) simulation, and another in which WRF and WW3 were two-way coupled using the Coupled-Ocean-Atmosphere-Wave-Sediment-Transport (COAWST) framework. It was found that all the 5-km WRF simulations improved 90-m wind speed statistics for the tropical cyclone case of 08 May 2015 and the cold frontal case of 25 Mar 2016, but not the nor-easter of 18 Jan 2016. The impact of wave coupling on buoy-level (4 m) wind speed was modest and case dependent, but when present, the impact was typically seen at 90 m as well, being as large as 10% in stable conditions. One-way wave coupling consistently reduced wind speeds, improving biases for 25 Mar 2016 but worsening them for 08 May 2015. Two-way wave coupling mitigated these negative biases, improved wave field representation and statistics, and mostly improved 4-m wind field correlation coefficients, at least at the VA buoy, largely due to greater self-consistency between wind and wave fields.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Exceptional coupling in photonic anisotropic metamaterials for extremely low waveguide crosstalk

Electromagnetic coupling is ubiquitous in photonic systems and transfers optical signals from one device to the other, creating crosstalk between devices. While this allows the functionality of some photonic components such as couplers, it limits the integration density of photonic chips, and many approaches have been proposed to reduce the crosstalk. However, due to the wave nature of light, complete elimination of crosstalk between closely spaced, identical waveguides is believed to be impossible and has not been observed experimentally. Here we show an exceptional coupling that can completely suppresses the crosstalk utilizing highly anisotropic photonic metamaterials. The anisotropic dielectric perturbations in the metamaterial mutually cancel the couplings from different field components, resulting in an infinitely long coupling length. We demonstrate the extreme suppression of crosstalk via exceptional coupling on a silicon-on-insulator platform, which is compatible with a complementary metal-oxide-semiconductor process. The idea of exceptional coupling with anisotropic metamaterials can be applied to many other electromagnetic devices, and it could drastically increase the integration density of photonic chips.

Mia, Md Borhan (ORCID:0000000320851331)↗

DOE OSTI2024

Exploring the impact of vibrational cavity coupling strength on ultrafast CN + c -C 6 H 12 reaction dynamics

Abstract Molecular polaritons, hybrid light-matter states resulting from strong cavity coupling of optical transitions, may provide a new route to guide chemical reactions. However, demonstrations of cavity-modified reactivity in clean benchmark systems are still needed to clarify the mechanisms and scope of polariton chemistry. Here, we use transient absorption to observe the ultrafast dynamics of CN radicals interacting with a cyclohexane ( c -C 6 H 12 ) and chloroform (CHCl 3 ) solvent mixture under vibrational strong coupling of a C–H stretching mode of c- C 6 H 12 . By modulating the c -C 6 H 12 :CHCl 3 ratio, we explore how solvent complexation and hydrogen (H)-abstraction processes proceed under collective cavity coupling strengths ranging from 55 to 85 cm −1 . Reaction rates remain unchanged for all extracavity, on-resonance, and off-resonance cavity coupling conditions, regardless of coupling strength. These results suggest that insufficient vibrational cavity coupling strength may not be the determining factor for the negligible cavity effects observed previously in H-abstraction reactions of CN with CHCl 3 .

Chen, Liying↗

DOE OSTI2024

Machine Learning-Based Identification of the Interface Regions for Coupling Local and Nonlocal Models

Local-nonlocal coupling approaches provide a means to combine the computational efficiency of local models and the accuracy of nonlocal models. However, the coupling process can be challenging, requiring expertise to identify the interface between local and nonlocal regions. Here, this study introduces a machine learning-based approach to automatically detect the regions in which the local and nonlocal models should be used. The method uses loading functions evaluated at grid points to decide the model selection at those points. Training of the networks is based on datasets provided by classes of loading functions for which reference coupling configurations are computed using accurate coupled solutions, where accuracy is measured in terms of the relative error between the solution to the coupling approach and the solution to the nonlocal model. We study two approaches that vary in data structure. The first, the full-domain input data approach, uses the entire load vector and outputs a complete label vector, performing a global classification. The second, a window-based approach, processes loads into windows and addresses the problem as a node-wise classification where each window's central point is classified individually. The classification problems are solved via deep learning algorithms based on convolutional neural networks. The performance of these approaches is studied on one-dimensional numerical examples using F1-scores and accuracy metrics. Notably, the windowing approach achieves an accuracy of 0.96 and an F1-score of 0.97, highlighting its potential to automate coupling processes effectively and enhance computational efficiency in material science applications.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2019

Coupled Hydrodynamic-Geomechanical Modelling of CO2-WAG Field Development at Farnsworth Unit: A Case Study

Abstract The SWP project is located in a mature waterflood undergoing conversion to CO2-WAG operations at Farnsworth, Texas, USA. Utilized CO2 is anthropogenic, sourced from a fertilizer and an ethanol plant. Major project goals are optimizing the storage/production balance, ensuring storage permanence, and developing best practices for CCUS. This paper provides a review of work performed toward development of a 3D coupled Mechanical Earth Model (MEM) for use in assessment of caprock integrity, fault reactivation potential, and evaluation of stress dependent permeability in reservoir forecasting. Mechanical property estimates computed from geophysical logs at selected wellbores were integrated with 3D seismic elastic inversion products to create a 3D "static" mechanical property model sharing the same geological framework as the existing reservoir simulation model including 3 major faults. Stresses in the MEM were initialized from wellbore stress estimates and reservoir simulation pore pressures. One way and two way coupled simulations were performed using a compositional hydrodynamic flow model and geomechanical solvers. Coupled simulations were performed on history matched primary, secondary (waterflood), and tertiary (CO2 WAG) recovery periods, as well as an optimized WAG prediction period. These simulations suggest that the field has been operating at conditions which are not conducive to either caprock failure or fault reactivation. Two way coupled simulations were performed in which permeability was periodically updated as a function of volumetric strain using the Kozeny-Carmen porosity-permeability relationship. These simulations illustrate the importance of frequent permeability updating when recovery scenarios result in large pressure changes such as in field re-pressurization through waterflood after a long primary depletion recovery period. Conversely, production forecasting results are less sensitive to permeability update frequency when pressure cycles are short and shallow as in WAG cycles. This paper describes initial work on development of a mechanical earth model for use in assessment of geomechanical risks associated with CCUS operations at FWU. The emphasis of this work is on integration of available geomechanical data for creation of the static mechanical property model. Preliminary coupled hydro-mechanical simulations are presented to illustrate some of the key diagnostic output from coupled simulations which will be used in later work for in depth evaluation of specific risk factors such as induced seismicity and caprock integrity.

02 PETROLEUM↗

DOE OSTI2014

Simulating High Flux Isotope Reactor Core Thermal-Hydraulics via Interdimensional Model Coupling

A coupled interdimensional model is presented for the simulation of the thermal-hydraulic characteristics of the High Flux Isotope Reactor core at Oak Ridge National Laboratory. The model consists of two domains a solid involute fuel plate and the surrounding liquid coolant channel. The fuel plate is modeled explicitly in three-dimensions. The coolant channel is approximated as a twodimensional slice oriented perpendicular to the fuel plate s surface. The two dimensionally-inconsistent domains are linked to one another via interdimensional model coupling mechanisms. The coupled model is presented as a simplified alternative to a fully explicit, fully three-dimensional model. Involute geometries were constructed in SolidWorks. Derivations of the involute construction equations are presented. Geometries were then imported into COMSOL Multiphysics for simulation and modeling. Both models are described in detail so as to highlight their respective attributes in the 3D model, the pursuit of an accurate, reliable, and complete solution; in the coupled model, the intent to simplify the modeling domain as much as possible without affecting significant alterations to the solution. The coupled model was created with the goal of permitting larger portions of the reactor core to be modeled at once without a significant sacrifice to solution integrity. As such, particular care is given to validating incorporated model simplifications. To the greatest extent possible, the decrease in solution time as well as computational cost are quantified versus the effects such gains have on the solution quality. A variant of the coupled model which sufficiently balances these three solution characteristics is presented alongside the more comprehensive 3D model for comparison and validation.

COMSOL HFIR LEU↗

DOE OSTI2017

Warthog: Coupling Status Update

The Warthog code was developed to couple codes that are developed in both the Multi-Physics Object-Oriented Simulation Environment (MOOSE) from Idaho National Laboratory (INL) and SHARP from Argonne National Laboratory (ANL). The initial phase of this work, focused on coupling the neutronics code PROTEUS with the fuel performance code BISON. The main technical challenge involves mapping the power density solution determined by PROTEUS to the fuel in BISON. This presents a challenge since PROTEUS uses the MOAB mesh format, but BISON, like all other MOOSE codes, uses the libMesh format. When coupling the different codes, one must consider that Warthog is a light-weight MOOSE-based program that uses the Data Transfer Kit (DTK) to transfer data between the various mesh types. Users set up inputs for the codes they want to run, and then Warthog transfers the data between them. Currently Warthog supports XSProc from SCALE or the Sub-Group Application Programming Interface (SGAPI) in PROTEUS for generating cross sections. It supports arbitrary geometries using PROTEUS and BISON. DTK will transfer power densities and temperatures between the codes where the domains overlap. In the past fiscal year (FY), much work has gone into demonstrating two-way coupling for simple pin cells of various materials. XSProc was used to calculate the cross sections, which were then passed to PROTEUS in an external file. PROTEUS calculates the fission/power density, and Warthog uses DTK to pass this information to BISON, where it is used as the heat source. BISON then calculates the temperature profile of the pin cell and sends it back to XSProc to obtain the temperature corrected cross sections. This process is repeated until the convergence criteria (tolerance on BISON solve, or number of time steps) is reached. Models have been constructed and run for both uranium oxide and uranium silicide fuels. These models demonstrate a clear difference in power shape that is not accounted for in a stand-alone BISON run. Future work involves improving the user interface (UI), likely through integration with the Nuclear Energy Advanced Modeling and Simulation (NEAMS) Workbench. Furthermore, automating the input creation would ease the user experience. The next priority is to continue coupling the work with other codes in the SHARP package. Efforts on other projects include work to couple the Nek5000 thermo-hydraulics code to MOOSE, but this is in the preliminary stages.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2017

Progress on Fabrication of Planar Diffusion Couples with Representative TRISO PyC/SiC Microstructure

Release of fission products from tristructural-isotropic (TRISO) coated particle fuel limits the fuel’s operational lifetime and creates potential safety and maintenance concerns. A need for diffusion analysis in representative TRISO layers exists to provide fuel performance models with high fidelity data to improve fuel performance and efficiency. An effort has been initiated to better understand fission product transport in, and release from, quality TRISO fuel by investigating diffusion couples with representative pyrocarbon (PyC) and silicon carbide (SiC). Here planar PyC/SiC diffusion couples are being developed with representative PyC/SiC layers using a fluidized bed chemical vapor deposition (FBCVD) system identical to those used to produce laboratory-scale TRISO fuel for the Advanced Gas Reactor Fuel Qualification and Development Program’s (AGR) first fuel irradiation. The diffusivity of silver, the silver and palladium system, europium, and strontium in the PyC/SiC will be studied at elevated temperatures and under high temperature neutron irradiation. The study also includes a comparative study of PyC/SiC diffusion couples with varying TRISO layer properties to understand the influence of SiC microstructure (grain size) and the PyC/SiC interface on fission product transport. The first step in accomplishing these goals is the development of the planar diffusion couples. The diffusion couple construction consists of multiple steps which includes fabrication of the primary PyC/SiC structures with targeted layer properties, introduction of fission product species and seal coating to create an isolated system. Coating development has shown planar PyC/SiC diffusion couples with similar properties to AGR TRISO fuel can be produced. A summary of the coating development process, characterization methods, and status are presented.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2018

GMD Coupling to Power Systems and Disturbance Mitigation

Presentation includes slides on Geomagnetic Disturbance: Ground Fields; Geomagnetic Disturbance: Coupling to Bulk Electric System; Geomagnetic Disturbance: Transformers; GMD Assessment Workflow (TPL-007-1); FERC order 830; Goals; SuperMag (1 min data) Nov. 20-21, 2003 Storm (DST = -422); Spherical Harmonics; Spherical Harmonics Nov. 20-21, 2003 Storm (DST = -422); DST vs HN0,0; Fluctuations vs. DST; Fluctuations; Conclusions and Next Steps; GMD Assessment Workflow (TPL-007-1); EMP E3 Coupling to Texas 2000 Bus Model; E3 Coupling Comparison (total GIC) Varying Ground Zero; E3 Coupling Comparison (total MVAR) Varying Ground Zero; E3 Coupling Comparison (GIC) at Peak Ground Zero; E3 Coupling Comparison (GIC) at Peak Ground Zero; and Conclusion.

58 GEOSCIENCES↗

DOE OSTI2019

Selection of glasses in support of the Sludge Batch 9 variability study for coupled operation with the Salt Waste Processing Facility

Prior to the processing of Sludge Batch 9 (SB9) with Frit 803 at the Defense Waste Processing Facility (DWPF), assessments and experimental glass work were completed to support the qualification of SB9 coupled processing with streams from the Actinide Removal Process (ARP) and Modular Caustic Side Solvent Extraction Unit (MCU). This previous work did not address coupled processing with streams from the Salt Waste Processing Facility (SWPF) because the Product Composition Control System (PCCS) models limited the TiO 2 concentration in glass to 2 weight percent (wt.%). These high activity streams from SWPF include monosodium titanate (MST) and sludge solids from the Sludge Solids Receipt Tank as well as Cs-containing strip effluent. The incorporation of these SWPF streams is expected to reach TiO2 concentrations in glass greater than 2 wt.% based on an MST addition of 0.4 g MST/L of salt solution. If double strike operation is required, then the MST addition is increased to 0.8 g MST/L of salt solution. Therefore, revisions to the models contained in PCCS were completed to allow the evaluation of glasses containing greater than 2 wt.% TiO 2 , but less than 6 wt.%. To support the implementation of the revised PCCS for SB9 – Frit 803 processing at DWPF prior to SWPF startup, the applicability of the new durability models for the SB9 variability study glasses was verified. Since the initial SB9 qualification for ARP-MCU coupled operation, additional Tank 40 analyses have been completed, which served as the basis for an updated projection of the SB9 blend composition in Tank 40. To identify a frit optimized for SB9 coupled operation with SWPF, the Savannah River National Laboratory (SRNL) performed a paper-study by utilizing Measurement Acceptance Region (MAR) assessments with DWPF PCCS. Based on these MAR assessment results, Frit 625 was recommended for processing both single and double MST strike scenarios. Savannah River Remediation (SRR) requested that SRNL perform a variability study with both Frit 803 and Frit 625 for SB9 coupled operation with SWPF. The objective of the variability study is to demonstrate that the Product Consistency Test (PCT) responses of the glasses are acceptable relative to the chemical durability of the Environmental Assessment (EA) reference glass and predictable by the current DWPF PCCS models for durability. This report documents the development of a test matrix that will support the variability study. To provide operational flexibility prior to SWPF startup, the predictability of glasses within the glass region for Frit 625-SB9 coupled operation with ARP-MCU was also confirmed.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2020

Coupling CTH to Linear Acoustic Propagation across an Air-Earth Interface

The interface between the Earth and the atmosphere forms a strong contrast in material properties. As such, numerical issues can arise when simulating an elastic wavefield across such a boundary when using a numerical simulation scheme. This is exacerbated when two different simulation codes are coupled straddling that interface. In this report we document how we implement the coupling of CTH, a nonlinear shock physics code, to a linearized elastic/acoustic wave propagation algorithm, axiElasti, across the air-earth interface. We first qualitatively verify that this stable coupling between the two algorithms produces expected results with no visible effects of the coupling interface. We then verify the coupling interface quantitatively by checking consistency with results from previous work and with coupled acoustic-elastic seismo-acoustic source inversions in three earth materials.

58 GEOSCIENCES↗

DOE OSTI2020

Coupling of CTF and TRANSFORM using the Functional Mockup Interface

An in-memory coupling between the sub-channel thermal hydraulics code COBRA-TF (CTF), which is included in the Virtual Environment for Reactor Applications (VERA), and the systems code Transient Simulation Framework of Reconfigurable Models (TRANSFORM) was developed in this work. Data exchange is accomplished by using the Functional Mock-Up Interface (FMI), an open standard for coupling dynamic systems models together. The FMI-based coupling necessitated the development of a novel FORTRAN wrapper for communicating with Functional Mock-Up Units. This wrapper will facilitate future FMI-based code couplings with other FORTRAN-based programs. CTF-TRANSFORM is exercised on a simplified Molten Salt Reactor Experiment (MSRE) model in steady-state and transient configurations. The coupled CTF-TRANSFORM model is shown to predict core temperature deltas similar to those available in historical MSRE operational data and is shown to be robust to fast power transients.The coupling paves the way to performing sensitivity and design studies impossible with VERA/CTF alone, in which secondary and tertiary loop operating conditions and design parameters may be perturbed and their impact on the core operating conditions studied.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

Coupled Ocean Wave-Atmosphere Models for Offshore Wind Energy Applications

Coupling ocean wave models to mesoscale atmospheric models is necessary to represent the effect of waves on wind turbine hub-height winds. In this report we provide a review of the most widely used ocean waves models and the phased-averaged spectral wave modeling paradigm that they are based on. Methodologies used to couple these wave models to mesoscale atmospheric models are described along with details of existing coupled modeling systems. We summarize impacts on offshore wind resource assessments that have been investigated with such coupled modeling systems to-date. Specifically, coupling in the North and Baltic Seas was shown to have a small negative effect on the offshore wind energy resource at weak-to-moderate wind speeds, with little impact at higher wind speeds. Finally, limitations of these existing coupled modeling systems and impact assessments are discussed, including examples of and the potential for improved parameterizations of the air-sea fluxes, and the use of fine-scaled simulations using phased-resolved wave models.

17 WIND ENERGY↗

DOE OSTI2022

Initial Investigation of Coupled Electron-Photon Transport with SCEPTRE and RAPTURE (Rev.1.1)

The SCEPTRE1 code models coupled photon and electron transport in material media using discrete ordinates methods. The code can use different solvers in 1D and 3D in parallel, promising users the ability to run problems that are significant both in physics and computational scale. One capability that is perhaps unique at this time in the NNSA complex is SCEPTRE’s coupled photon– electron transport physics in a discrete ordinates approach. For example, incident x–rays can generate a significant electron population, whose transport can be simulated within the same code. In this context, fully–coupled means that photons can produce electrons and electrons can produce photons; uncoupled or no–coupling means photons do not produce electrons. Partial–coupling, an intermediate option in which photons can produce electrons but electrons do not produce photons, is not considered in the present work. The MCNP code can treat this physics with a Monte Carlo approach. MCNP should be viewed as complementary to SCEPTRE in both computational method (Monte Carlo versus discrete ordinates) as well as data libraries. SCEPTRE and MCNP should provide important checks on each other, especially from complementary modeling frameworks of a S n (discrete ordinates) mesh based method versus Monte Carlo particle based method. In order to learn more about SCEPTRE, and how it might inform engineering work at the Laboratory, we used it to carry out a set of initial studies. This was a series of 1D simulations of photon beams hitting Aluminum slabs. Our studies mainly used SCEPTRE directly, and made little use of the closely associated RAPTURE code (more on this below). The goals of these studies were first, to familiarize ourselves with operating the code, and second, to begin to assess the importance of coupled electron–photon transport. To support these studies, we began to develop a set of scripts to initialize inputs, run tools, process outputs, and plot results. An initial version of these scripts has been committed to the Lightning workflow repository. This memo relates the results of those studies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗