Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “D CODES”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

CEFR simulation using diffusion code system RAST-F

In this study, the RAST-F code system, which is based on the two-step approach that couples a multi-group cross-section generation Monte-Carlo (MC) code and a multi-group nodal diffusion solver, was used for the neutronic simulation of the CEFR start-up experiments. The numerical solution of the RAST-F system was verified against the full core MC solution MCS at all control rods fully inserted and withdrawn states. The RAST-F solution of the selected experimental simulations was compared against the measurement data. A good agreement between RAST-F and MCS solutions was observed with less than 120 pcm discrepancies and 1.2% root-mean-square error in terms of k{sub eff} and power distribution, respectively. Meanwhile, the RAST-F result agreed well with the experimental data within two-sigma of experimental uncertainty. The good agreement of these results indicates that RAST-F can be used to neutronic steady-state simulations for small core-size SFR, which was challenged to deterministic code systems. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2024

NDE of Novel Materials (CFRP) [Slides]

CFRP composites have not been used for nuclear safety-related applications until recently. In 2019, the ASME Boiler and Pressure Vessel Code Committee approved a new Code Case N-871 for internal repairs of Class 2 and 3 safety-related piping using CFRP for Service Levels A, B, C, and D for a service life of 50 years. The NRC did not review N-871 for inclusion in the Code Case Regulatory Guides. This project focuses on evaluating the capabilities and limitations of Nondestructive Evaluation (NDE) methods for examining the Carbon Fiber Reinforced Polymer (CFRP) repairs in commercial nuclear power plants (NPPs). This Report summarizes the work done on this project until the end of FY24.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Intrinsic Toroidal Rotation Driven by Turbulent and Neoclassical Processes in Tokamak Plasmas from Global Gyrokinetic Simulations

Gyrokinetic tokamak plasmas can exhibit intrinsic toroidal rotation driven by the residual stress. While most studies have attributed the residual stress to the parallel-momentum flux from the turbulent E × B motion, the parallel-momentum flux from the drift-orbit motion (denoted $Π^D_\parallel$) and the E × B-momentum flux from the E × B motion (denoted $Π_{E×B}$) are often neglected. Here, we use the global total-f gyrokinetic code XGC to study the residual stress in the core and the edge of a DIII-D H-mode plasma. Numerical results show that both $Π^D_\parallel$ and $Π_{E×B}$ make up a significant portion of the residual stress. In particular, $Π^D_\parallel$ in the core is higher than the collisional neoclassical level in the presence of turbulence, while in the edge it represents an outflux of countercurrent momentum even without turbulence. Using a recently developed “orbit-flux” formulation, we show that the higher-than-neoclassical-level $Π^D_\parallel$ in the core is driven by turbulence, while the outflux of countercurrent momentum from the edge is mainly due to collisional ion orbit loss. In conclusion, these results suggest that $Π^D_\parallel$ and $Π_{E×B}$ can be important for the study of intrinsic toroidal rotation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Safety Analysis Report for Packagings – Model 9979 Type AF Shipping Package

This Safety Analysis Report for Packaging (SARP) documents the performance of the Department of Energy (DOE) Model USA/9979/AF-96 shipping package and compliance with the regulatory safety requirements of Title 10 and Title 49 of the Code of Federal Regulations. The SARP is prepared in accordance with U.S. Department of Energy (DOE) Order 460.1D. The Model 9979 Type A Fissile shipping package replaces the Department of Transportation (DOT) UN1A2 Specification Packaging. Formerly, the UN1A2 was authorized under Title 49 of the Code of Federal Regulations, Section 173.417 (2004) for transport of Type A quantities of fissile material over public highways. In accordance with the Federal Register, final rule making issued on January 26, 2004, the Nuclear Regulatory Commission, in a joint effort with the Department of Transportation, initiated a planned phase out of multiple Specification Packagings from 49 CFR 173, the UN1A2 being one. The rulemaking mandated, by law, removal of these specification packagings from service no later than September 30, 2008. The 9979 is a Performance Packaging and replaces the UN1A2 Specification Packaging. The Packaging Technology and Transportation Engineering (PT&TE) organization of the Savannah River National Laboratory (SRNL) is the Design Authority and Design Agency for the 9979.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

Technical Interface Coordination with SHINE (FY2020 Summary)

In FY20, Los Alamos National Laboratory (LANL) performed technical interface coordination with SHINE including coordination of SHINE CFD and thermal hydraulics conference calls, national laboratory conference calls, and program management activities such as reporting and tracking requirements for current FY20 activities. SHINE did not request any general technology consultation support for emergent issues. LANL has participated in the bi-weekly conference calls with SHINE and other national laboratories regarding LANL’s current activities. Activities discussed include the Supo 3-D CFD model geometry validation, releasing the CFD-MCNP coupled codes for the generic solution vessel simulations, and the ANL bubble experiment #2 design, preparation, and collaboration.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2023

P2M Simulation Exercise on Past Fuel Melting Irradiation Experiments

This paper presents the results of the Power To Melt and Maneuverability (P2M) Simulation Exercise on past fuel melting irradiation experiments, organized within the Organisation for Economic Co-operation and Development/Nuclear Energy Agency Framework for IrraDiation ExperimentS (FIDES) framework by the Core Group (CEA, EDF, and SCK·CEN) and open to all FIDES members. The exercise consisted in simulating two past power ramps where fuel melting was detected: (1) the xM3 staircase power transient [ramp terminal level (RTL) 70 kW·m -1 , average burnup 27 GWd·tU -1 ], carried out in 2005 in the R2 reactor at Studsvik (Sweden), where the rodlet maintained its integrity, and (2) the HBC4 fast power transient (RTL 66 kW·m -1 , average burnup 48 GWd·tU -1 ), carried out in 1987 in the BR2 reactor at SCK·CEN (Belgium), where the cladding failed during the experiment. The exercise was joined by 13 organizations from 9 countries using 11 different fuel performance codes. In this paper, the main results of the Simulation Exercise are presented and compared to available postirradiation examinations (PIE) or on-line measurements during the power ramps (fuel and clad diameters, rod elongation, pellet-clad gap, and fission gas release). Since the focus of the Simulation Exercise is on fuel melting assessment, determination of the boundary between melted/nonmelted fuel and the consequent definition of a melting radius from PIE are first discussed. During the HBC4 ramp, fuel melting was predicted by most of the codes despite differences in the melting models. Higher discrepancies were observed for the xM3 rod that can be attributed partly to power uncertainty and partly to the limited capability of the models to describe partial melting of the fuel during this ramp. Finally, possible code developments to improve simulation results are presented.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Runaway electron deconfinement in SPARC and DIII-D by a passive 3D coil

The operation of a 3D coil—passively driven by the current quench (CQ) loop voltage—for the deconfinement of runaway electrons (REs) is modeled for disruption scenarios in the SPARC and DIII-D tokamaks. Nonlinear magnetohydrodynamic (MHD) modeling is carried out with the NIMROD code including time-dependent magnetic field boundary conditions to simulate the effect of the coil. Further modeling in some cases uses the ASCOT5 code to calculate advection and diffusion coefficients for REs based on the NIMROD-calculated fields, and the DREAM code to compute the runaway evolution in the presence of these transport coefficients. Compared with similar modeling in Tinguely et al (2021 Nucl. Fusion 61 124003), considerably more conservative assumptions are made with the ASCOT5 results, zeroing low levels of transport, particularly in regions in which closed flux surfaces have reformed. Of three coil geometries considered in SPARC, only the n = 1 coil is found to have sufficient resonant components to suppress the runaway current growth. Without the new conservative transport assumptions, full suppression of the RE current is maintained when the thermal quench MHD is included in the simulation or when the RE current is limited to 250kA, but when transport in closed flux regions is fully suppressed, these scenarios allow RE beams on the order of 1–2 MA to appear. Additional modeling is performed to consider the effects of the close ideal wall. In DIII-D, the CQ is modeled for both limited and diverted equilibrium shapes. In the limited shape, the onset of stochasticity is found to be insensitive to the coil current amplitude and governed largely by the evolution of the safety-factor profile. In both devices, prediction of the q-profile evolution is seen to be critical to predicting the later time effects of the coil.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Validation of EDGDE2D-EIRENE predicted 2D distributions of electron temperature and density against divertor Thomson scattering measurements in the low-field side divertor leg in DIII-D

Comparisons of profiles of the electron temperature (T e ), density (n e ) and pressure (p e ) measured with Divertor Thomson Scattering in DIII-D low-confinement mode discharges to predictions from the edge fluid code EDGE2D-EIRENE show that the models implemented in EDGE2D-EIRENE predict the measurements within their collective uncertainties if the T e at the separatrix (T e,sep ) is 10 eV, or higher. The simulations do not predict, however, the peaked T e and ne profiles measured adjacent to the target plate when T e,sep , is below 10 eV, i.e., for the plasma downstream from the region of ionization of deuterium atoms. Inclusion of cross-field drifts and a five-fold reduction of radial transport cannot reconcile the discrepancy between the measurements and predictions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Thermal Neutron Scattering Law (TNSL) Implementation and Testing in FUDGE

Thermal neutron scattering law processing capabilities have been recently implemented in the LLNL code FUDGE (For Updating Data and Generating Evaluations). FUDGE is now capable of producing processed TNSL data for use in Monte Carlo or deterministic transport. To test this new capability, LLNL scientists performed an extensive intercomparison between the Mercury, Ardra, COG and MCNP transport codes. This intercomparison helped probe differences between how TNSL data are handled by the transport codes and by two different processing codes (LLNL’s FUDGE and LANL’s NJOY). This report summarizes recent improvements in FUDGE TNSL capabilities as well as results from the code intercomparison.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Deinococcus radiodurans UWO298 Dependence on Background Radiation for Optimal Growth

Ionizing radiation is a major environmental variable for cells on Earth, and so organisms have adapted to either prevent or to repair damages caused by it, primarily from the appearance and accumulation of reactive oxygen species (ROS). In this study, we measured the differential gene expression in Deinococcus radiodurans UWO298 cultures deprived of background ionizing radiation (IR) while growing 605 m underground at the Waste Isolation Pilot Plant (WIPP), reducing the dose rate from 72.1 to 0.9 nGy h –1 from control to treatment, respectively. This reduction in IR dose rate delayed the entry into the exponential phase of the IR-shielded cultures, resulting in a lower biomass accumulation for the duration of the experiment. The RNASeq-based transcriptome analysis showed the differential expression of 0.2 and 2.7% of the D. radiodurans genome after 24 and 34 h of growth in liquid culture, respectively. Gene expression regulation after 34 h was characterized by the downregulation of genes involved in folding newly synthesized and denatured/misfolded proteins, in the assimilation of nitrogen for amino acid synthesis and in the control of copper transport and homeostasis to prevent oxidative stress. We also observed the upregulation of genes coding for proteins with transport and cell wall assembly roles. These results show that D. radiodurans is sensitive to the absence of background levels of ionizing radiation and suggest that its transcriptional response is insufficient to maintain optimal growth.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Code Verification and Solution Verification framework in pin-resolved neutron transport code MPACT

Program verification in scientific computing encompasses the application of formal and mathematical techniques to a scientific computing code for its credibility, accuracy, and validity. Code Verification identifies bugs and performance issues in the software development stage. Solution Verification assesses the applicability of the code and the accuracy of the solution to problems of interest. Both activities utilize application cases and quantify the error against prescribed acceptance criteria. However, simply executing more application cases does not guarantee stronger or more comprehensive credibility. Here, we establish a verification framework that involves Code Verification and Solution Verification, both of which work together such that the overarching goal of “converge to the correct answer for the intended application” can be reasonably inferred. The application of such a verification framework is demonstrated using the pin-resolved neutron transport code MPACT, where standard unit tests and regression tests are covered, and where the Method of Exact Solutions and the Method of Manufactured Solutions are successfully used. Additionally, the applicability of Method of Manufactured Solutions is extended to the OECD/NEA C5G7 benchmark problems of practical material and geometric configurations. Solution Verification activities are demonstrated on a practical hierarchy of application models of increasing complexity ranging from 2D pin cell problems to 3D assembly problems. The convergence behavior and rate of convergence with respect to each individual variable are studied and provided. This framework can be adapted broadly to other fields involving scientific computing codes.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Progress in DIII-D towards validating divertor power exhaust predictions

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

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Radiation-Induced Error Mitigation by Read-Retry Technique for MLC 3-D NAND Flash Memory

Here, we have evaluated the Read-Retry (RR) functionality of the 3-D NAND chip of multilevel-cell (MLC) configuration after total ionization dose (TID) exposure. The RR function is typically offered in the high-density state-of-the-art NAND memory chips to recover data once the default memory read method fails to correct data with error correction codes (ECCs). In this work, we have applied the RR method on the irradiated 3-D NAND chip that was exposed with a Co-60 gamma-ray source for TID up to 50 krad (Si). Based on our experimental evaluation results, we have proposed an algorithm to efficiently implement the RR method to extend the radiation tolerance of the NAND memory chip. Our experimental evaluation shows that the RR method coupled with ECC can ensure data integrity of MLC 3-D NAND for TID up to 50 krad (Si).

3-D NAND↗
DOE OSTI2023

ELM sputter erosion modeling of a tungsten coated small angle slot divertor in DIII-D

We modeled plasma edge localized mode (ELM) sputter erosion for a Small Angle Slot divertor with a tungsten coated region (SAS-VW), designed for experiments in the DIII-D tokamak, and proposed for use in future advanced tokamaks. The simulations use a free-streaming, 1000 eV, C +6 and D +1 ELM impingement model, with SOLPS-ITER, ITMC-DYN, and REDEP/WBC code packages for background plasma, material response, and erosion/redeposition respectively. The results show ELM'ing plasma gross and net tungsten erosion fluxes of the mixed-material C/W surface peaking at the slot entrance region, and an order of magnitude higher than for non-ELMs. The per-pulse erosion, however, remains low, of order 0.5 nm, due to expected moderate ELM frequencies and duration in DIII-D. Here, the ELMs result in a ~25x higher peak sputtered W current leaving the divertor slot region, towards the core plasma, compared to the ELM-free plasma case. The time-integrated escape current, however, may not significantly affect core plasma high-Z contamination concerns, for a 1% ELM duty factor, but may be an issue for higher frequency ELMs. In general, the modeling results appear favorable for effective testing of the SAS-VW divertor in DIII-D, and extrapolation to innovative divertor designs in future ITER-like and DEMO fusion devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Assessing the impact of alpha particles on thermal confinement in JET D-T plasmas through global GENE-Tango simulations

The capability of the global, electromagnetic gyrokinetic GENE code interfaced with the transport Tango solver is exploited to address the impact of fusion alpha particles (in their dual role of fast particles and heating source) on plasma profiles and performance at JET in the discharges with the highest quasi-stationary peak fusion power during the DTE2 experimental campaigns. Employing radially global nonlinear electromagnetic GENE-Tango simulations, we compare results with/without alpha particles and alpha heating. Our findings reveal that alpha particles have a negligible impact on turbulent transport, with GENE-Tango converging to similar plasma profiles regardless of their inclusion as a kinetic species in GENE. On the other hand, alpha heating is found to contribute to the peaking of the electron temperature profiles, leading to a 1 keV drop on the on-axis electron temperature when alpha heating is neglected in Tango. The minimal impact of alpha particles on turbulent transport in this JET discharge–despite this being the shot with the highest fusion output–is attributed to the low content of fusion alpha in this discharge. To assess the potential impact of alpha particles on turbulent transport in regimes with higher alpha particle density, as expected in ITER and fusion reactors, we artificially increased the alpha particle concentration to levels expected for ITER. By performing global nonlinear GENE standalone simulations, we found that increasing the alpha particle density beyond five times the nominal value lead to significant overall turbulence destabilization. These results demonstrate that an increased alpha particle concentration can significantly impact transport properties under simulated JET experimental conditions. However, these findings cannot be directly extrapolated to ITER due to the substantial differences in parameters such as plasma size, magnetic field, plasma current, and thermal pressure.

energetic particles↗
DOE OSTI2022

Modeling of a Generic Pebble Bed High-temperature Gas-cooled Reactor (PB-HTGR) with SAM

This report presents the modeling of the core of a generic pebble-bed reactor (PBR) at the system level using the System Analysis Module (SAM) code. This work is an extension of a previous work by the authors (Ooi et al. (2021)) that used the so-called 2-D ring model approach to model the PBMR-400. With the new approach, the pebble bed of the reactor is modeled with multiple PBCoreChannel components with spherical heat structures which allows the code to calculate thermal fluid parameters with built-in closure relations. The new core-channel approach is an improvement to the 2-D ring model approach as it does not introduce geometric distortions to the model and thus reduces the uncertainties of the predictions. In addition to thermal fluid simulations, point kinetics (PKE) are included to the model. Simulations are performed under a steady-state normal operation condition and a load-following transient scenario. This particular transient scenario is chosen as it tests both the thermal fluid and neutronics aspects of the model. The predicted results from both the steady-state and transient scenarios are compared with the results by Stew- art et al. (2021) who performed similar simulations with a Griffin-Pronghorn coupled tool. Despite the differences between the codes, with SAM being a system-analysis code and Pronghorn being a porous-medium code, both sets of results compare favorably. The over- all profiles and trends of the predicted temperatures and reactivities from the SAM and Griffin/Pronghorn simulations are similar, with some differences in their predicted values. The first part of the report covers the significance of a relatively fast-running approach that is capable of modeling the pebble bed reactor at the system-level while simultaneously capturing the radial thermal behavior of the core. Then, the modeling approach used in this work is discussed in details. Lastly, the results and comparisons with the Griffin/Pronghorn simulation are presented.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

LANL ASC FY23 L2 Milestone Description and Completion Criteria [Slides]

In early 2021, the ASC Program decided that Ristra should provide a production code (Moya) for LANL’s mission in the Low Energy-Density Physics application space. At the same time, Ristra should design an environment for rapid development of new codes that can answer new questions on new hardware. We started FY22 by expanding upon our R&D efforts of high-order DG methods for hydrodynamics; Technical & staffing challenges caused us to abandon this for a traditional FV SGH approach, but leveraging modern techniques. Significant changes in FleCSI (v1.4→2.x) have us writing code anew, in a co-design effort with the FleCSI and FleCSI specialization projects.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Computational exploration of DIII-D QH-mode with extended MHD

This work will enable improved understanding of the scientific basis for DIII-D Quiescent H-mode (QH-mode) by leveraging recent computational developments of the NIMROD code. Experimentally, it is established that the QH-mode tokamak operation regime is free of edge-localized modes (ELMs) with good energy confinement. ELMs produce an impulsive heat-flux on the divertor that is unacceptable under burning-plasma conditions. Beyond DIII-D in future high-power-density tokamaks, accommodation of a high heat flux at the divertor represents a key challenge to overcome. Proposed solutions involve operation at high divertor density to achieve divertor detachment, operation with a liquid lithium wall and/or introduction of impurities to enhance radiative losses. Additionally, the first wall material is likely to be a high-Z material such as Tungsten in future devices as materials such as graphite tiles restrict the lifetime of the first wall due to a high erosion rate. While the atomic species composition of the plasma is known to impact the edge-pedestal structure and stability, the specifics of the underlying physics basis for this interaction is not firm. Studying this interaction is the goal of this proposal without focusing specifically on the fine details of the first wall and divertor solution. It is the goal is to lay the ground work to establish the compatibility of the divertor solutions and wall material with pedestal solutions that eliminate ELMs by leveraging present DIII-D experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗