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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Massively parallel GPU enabled third-order cluster perturbation excitation energies for cost-effective large scale excitation energy calculations

We present here a massively parallel implementation of the recently developed CPS(D-3) excitation energy model that is based on cluster perturbation theory. The new algorithm extends the one developed in Baudin et al. [J. Chem. Phys., 150, 134110 (2019)] to leverage multiple nodes and utilize graphical processing units for the acceleration of heavy tensor contractions. Furthermore, we show that the extended algorithm scales efficiently with increasing amounts of computational resources and that the developed code enables CPS(D-3) excitation energy calculations on large molecular systems with a low time-to-solution. More specifically, calculations on systems with over 100 atoms and 1000 basis functions are possible in a few hours of wall clock time. This establishes CPS(D-3) excitation energies as a computationally efficient alternative to those obtained from the coupled-cluster singles and doubles model.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2026

Minimizing the Electromechanical Stresses in Poloidal Field Coils by Optimizing their Numbers and Locations using FREDA Framework

Poloidal field (PF) and central solenoid (CS) coils play a crucial role in sustaining the equilibrium and preserving the shape of highly confined tokamak plasmas. Ensuring that PF coil current and mechanical stress stay within superconducting and structural limitations is an important check in the design assessment. Minimizing the PF coil currents and mechanical stresses influences reliability, cost, and performance. A free-boundary MHD equilibrium code—FreeGS is employed within the fusion reactor design and assessment (FREDA) whole facility modeling (WFM) framework to construct the plasma equilibrium based on the configuration and currents in the PF coils. Here, we present the capability of the FreeGS code to minimize the currents, forces, and electromagnetic stresses on the PF coils by optimizing their number, sizes, structures, and locations while maintaining an MHD stable plasma configuration with a large confinement factor. The workflow is initialized with a configuration of plasma parameters and coils’ locations from the 0-D tokamak build systems code in the FREDA framework. Then, FreeGS is called to calculate the initial equilibrium at the minimum total current in PF coils. Thereafter, FreeGS’s internal optimizer minimizes the currents and hoop and central forces on the PF coils while maintaining the reference equilibrium. Finally, the input configuration is updated with the optimized parameters for equilibria over the ramp-up phase of a burning-plasma operation. FREDA’s whole facility optimization capability, which includes all magnetic field coil systems, blanket, vacuum vessel (VV), first wall, divertor, etc., is under development and out of the scope for this study.

Hassan, Ehab [ORNL] (ORCID:0000000181060301)↗
DOE OSTI2021

Bottomonium suppression in an open quantum system using the quantum trajectories method

We solve the Lindblad equation describing the Brownian motion of a Coulombic heavy quark-antiquark pair in a strongly coupled quark-gluon plasma using the highly efficient Monte Carlo wave-function method. The Lindblad equation has been derived in the framework of pNRQCD and fully accounts for the quantum and non-Abelian nature of the system. The hydrodynamics of the plasma is realistically implemented through a 3+1D dissipative hydrodynamics code. We compute the bottomonium nuclear modification factor and compare with the most recent LHC data. The computation does not rely on any free parameter, as it depends on two transport coefficients that have been evaluated independently in lattice QCD. Our final results, which include late-time feed down of excited states, agree well with the available data from LHC 5.02 TeV PbPb collisions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Nonlinear models for coupling the effects of stimulated Raman scattering to inertial confinement fusion codes

Laser plasma instabilities (LPI) reduce driver-target coupling, alter implosion symmetry, and therefore can fundamentally limit fusion performance in inertial confinement fusion (ICF). Developing a predictive modeling capability for LPI effects can critically advance the success of the field. We perform vector particle-in-cell simulations of multi-speckled laser beams undergoing stimulated Raman scattering (SRS) at various densities and intensities relevant to mainly indirectly driven and a subset of parameter space for directly driven ICF systems, focusing on the regimes with intensities above the onset of electron trapping. Based on the wavenumber of the SRS daughter electron plasma wave, we identify several regions with underpinning SRS saturation physics: the electron-trapping dominated region with intermediate kλD values, the strong Landau damping region at larger kλD values, and the region where the Langmuir decay instability arises at lower kλD values. We develop a nonlinear SRS reflectivity model that features the base trapping-dominated scaling of (kλD)−4 and its modifications. Electron trapping deforms the initialized electron distribution functions, and we have developed a new δf-Gaussian-mixture algorithm for an accurate characterization of the trapped hot electron population. With this SRS hot electron description, we construct a nonlinear energy deposition model and a hot electron source model—based on a modified Manley–Rowe relation—suitable for including SRS effects as a sub-grid module in a high-fidelity ICF design code.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Radiation burnthrough measurements to infer opacity at conditions close to the solar radiative zone–convective zone boundary

Recent measurements at the Sandia National Laboratory of the x-ray transmission of iron plasma have inferred opacities much higher than predicted by theory, which casts doubt on modeling of iron x-ray radiative opacity at conditions close to the solar convective zone-radiative zone boundary. An increased radiative opacity of the solar mixture, in particular iron, is a possible explanation for the disagreement in the position of the solar convection zone-radiative zone boundary as measured by helioseismology and predicted by modeling using the most recent photosphere analysis of the elemental composition. In this report we present data from radiation burnthrough experiments, which do not support a large increase in the opacity of iron at conditions close to the base of the solar convection zone and provide a constraint on the possible values of both the mean opacity and the opacity in the x-ray range of the Sandia experiments. The data agree with opacity values from current state-of-the-art opacity modeling using the CASSANDRA opacity code.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

McMillan electron lens in a system with space charge

Space charge (SC) forces of a circulating beam in a ring have both linear (defocusing) and nonlinear components, due to a nonuniform beam distribution. The linear component of SC forces produces a betatron tune shift, which is the largest for a zero-amplitude particle, while the nonlinear component produces an amplitude-dependent betatron tune spread. These SC effects are responsible for several undesirable phenomena in accelerators: emittance growth, particle losses, beam halo, etc. In this paper, we investigate the possibility to mitigate the distributed SC forces by a thin McMillan lens, providing an axially-symmetric kick, which is qualitatively opposite to the accumulated effect of beam's own SC. Experimentally, the proposed concept can be tested in Fermilab's IOTA ring. A thin McMillan lens can be implemented by a short (70) insertion of an electron beam with a specifically chosen density distribution in transverse directions. In this article, to test if McMillan lenses can reduce the tune spread induced by SC, we make several simulations with a 6-D particle tracking code, Synergia. We choose such beam and lattice parameters that the SC tune spread is roughly 0.5 and the emittance growth due to the half-integer resonance is clearly observed without the SC compensation. Then, we focus on reducing the emittance growth by adjusting the bare betatron tunes using the ring quadrupoles, and reducing the tune spread by the McMillan lenses. The results of reducing a large tune spread (≈ 0.5), reported here, are not perfect, but substantial. There is still room for further investigation. The simulations performed so far in this work indicate that McMillan lenses can cope with an SC tune spread of ≈ 0.1 per lens.

43 PARTICLE ACCELERATORS↗
DOE OSTI2020

Anomalous Absorption by the Two-Plasmon Decay Instability

Radiation-hydrodynamic simulations of directly driven fusion experiments at the Omega Laser Facility predict absorption accurately when targets are driven at low overlapped laser intensity. Discrepancies appear at increased intensity, however, with higher-than-expected laser absorption on target. Strong correlations with signatures of the two-plasmon decay (TPD) instability—including half-harmonic and hard-x-ray emission—indicate that TPD is responsible for this anomalous absorption. Scattered light data suggest that up to ≈ 30 % of the laser power reaching quarter-critical density can be absorbed locally when the TPD threshold is exceeded. A scaling of absorption versus TPD threshold parameter was empirically determined and validated using the laser–plasma simulation environment code.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

TOUGHREACT V4.12

TOUGHREACT V4.12 is a 3-D reactive transport code that handles to non-isothermal, multiphase/multicomponent fluid flow, heat transport, aqueous and gaseous species advection-diffusion, and equilibrium/kinetic water-gas-rock-biological reactions. TOUGHREACT V4.12 is based on TOUGHREACT V3.3, with several new features including: an option to compute activity coefficients using the Pitzer ion interaction model; implementation of ECO2N V2.0 for temperatures up to 300 C; options to simulate forward and reverse osmosis, including 2D Poiseuille flow; addition of EOS8 and revisions for EOS5 and EWASG with H2(g); addition of gas species decay, various new output formats and time step control options; many bug fixes.

Spycher, Nicolas↗
DOE OSTI2020

Development of a Conditioned System-Level Groundwater Model to Evaluate Long-Term Groundwater Impacts within a Performance Assessment - 20115

A preliminary performance assessment (PA) of single-shell tank Waste Management Area (WMA) A-AX located at the U.S. Department of Energy (DOE)'s Hanford Site in southeastern Washington is being conducted to satisfy the requirements of DOE Order 435.1 [1] as it relates to closure of the single-shell radioactive waste tanks in the WMA A-AX tank farms. A PA assesses the fate, transport, and impacts of radionuclides within a low-level radioactive waste disposal facility in its assumed closure configuration and the subsequent potential doses to humans over a 1,000-year compliance period and a 10,000-year performance evaluation period. The WMA A-AX preliminary PA evaluation is structured around the complementary use of process-level and system-level models to calculate the facility performance against established DOE Order 435.1 [1] performance objectives. Process-level models are those that represent a detailed phenomenological representation of processes of concern in the PA. Process models typically only represent one or a few of the components of the PA, such as groundwater flow and transport, and must be integrated with other modeling elements to perform PA calculations. System-level models are those that are abstracted from the process models, retaining the essential features of the process model, while allowing integration of all aspects of the PA in a single modeling framework. System-level models are often characterized by coarser numerical discretization, lower dimensionality, or other similar simplifications compared to the process-level model. Traditionally, a three-dimensional (3-D) process-level model is utilized primarily to evaluate the long-term impact on groundwater and the potential doses to individuals who consume contaminated groundwater. System-level models are also utilized to evaluate the groundwater pathway in PAs. These models typically have reduced dimensionality (1-D) and are conditioned utilizing flow fields (Darcy fluxes) and moisture content distributions that are abstracted from the process level models. The abstraction approach assures that the flow field in both models is consistent for a specific set of input parameters for flow, differing only in the discretization and dimensionality of the two models. The preliminary WMA A-AX PA incorporates a detailed representation of the geological system and hydraulic properties within the 3-D model STOMP{sup C} numerical code so that the effects of relevant features and processes on water flow and radionuclide transport in the subsurface can be evaluated. The complementary system-level model is developed utilizing the GoldSim{sup C} code to implement a simplified, 1-D equivalent model to represent the groundwater pathway. Contaminant transport through the vadose zone and unconfined aquifer for Tc-99 and I-129 were evaluated in each of the models. Adjustments to the saturated portion of the GoldSim{sup C} 1-D model were required to mimic the dispersion effect captured with the 3-D STOMP{sup C} model. Once this conditioning was conducted, highly similar results for the transport of Tc-99 and I-129 were achieved at the point of calculation, located 100 m downgradient from the WMA A-AX fenceline. These radionuclides represent elements that are regarded as primary dose drivers in the PA groundwater pathway analysis. The high degree of conformance between the two models suggests that the use of the equivalent 1-D system model is suitable for evaluating the full suite of radionuclides that will be released from the tank sources within WMA A-AX over the 1,000-year compliance period and over the 10,000-year evaluation period. (authors)

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

Fierro Version 2.x

FIERRO is a parallel C++ code designed to simulate fluid mechanics, heat transfer, and solid mechanics in two- and three-dimensional space. FIERRO is written to run on homogeneous (CPU) and heterogeneous (CPU+GPU) high performance computing machines. Fierro can aid a) modeling and design efforts that have historically relied on commercial implicit and explicit finite element codes, b) numerical methods research, c) manufacturing research, and d) computer science research. The code contains diverse numerical methods to solve the governing physics equations for both quasi-static and dynamic problems. Mathematical optimization solvers are coupled to the numerical methods to research topology and shape optimization that has application to additive manufacturing, and to create novel numerical approaches. Phase-field methods with micromechanical solvers are provided to simulate microstructure formation and evolution in manufacturing processes. The micromechanical solvers can also help research efforts create continuum-scale constitutive models for solids, as a function of the microstructure, in situ in a calculation or in a stand-alone manner. No physical data exists within the code.

Morgan, Nathaniel↗
DOE OSTI2023

Fierro

FIERRO is a parallel C++ code designed to simulate fluid mechanics, heat transfer, and solid mechanics in two- and three dimensional space. FIERRO is written to run on homogeneous (CPU) and heterogeneous (CPU+GPU) high performance computing machines. Fierro can aid a) modeling and design efforts that have historically relied on commercial implicit and explicit finite element codes, b) numerical methods research, c) manufacturing research, and d) computer science research. The code contains diverse numerical methods to solve the governing physics equations for both quasi-static and dynamic problems. Mathematical optimization solvers are coupled to the numerical methods to research topology and shape optimization that has application to additive manufacturing, and to create novel numerical approaches. Phase-field methods with micromechanical solvers are provided to simulate microstructure formation and evolution in manufacturing processes. The micromechanical solvers can also help research efforts create continuum-scale constitutive models for solids, as a function of the microstructure, in situ in a calculation or in a stand-alone manner. No physical data exists within the code.

Morgan, Nathaniel↗
DOE OSTI2026

Linear simulations of a wide pedestal quiescent H-mode plasma with the extended-MHD code NIMROD

We present linear simulations of a Wide Pedestal QH (WPQH)-mode DIII-D plasma with the NIMROD code to address the role of ion multispecies collisionality on WPQH-modes. We use a full extended-MHD model that captures Alfvénic and magnetoacoustic waves, including two-fluid, ion gyroviscosity and cross heat flux. We show that two-fluid and ion gyroviscous effects destabilize electron-directed peeling–ballooning modes in the pedestal. The linear growth rates of the modes decrease when carbon (C) is added to a deuterium plasma because the multispecies collisional effects increase the plasma resistivity. When replacing C with tungsten (W), the multispecies collisional effects increase the resistivity further, as well as the stabilizing effect on the instabilities. The resistivity impact on mode stability is reversed when both two-fluid and ion gyroviscous effects are ignored, which is consistent with previous works, proving that multispecies collisional effects, together with two-fluid and ion gyroviscosity, should be considered to address the stability of plasmas with significant impurity content.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Clifford Hierarchy Stabilizer Codes: Transversal Non-Clifford Gates and Magic States

A fundamental problem in fault-tolerant quantum computation is the tradeoff between universality and dimensionality, exemplified by the the Bravyi-König bound for $n$-dimensional topological stabilizer codes. In this work, we extend topological Pauli stabilizer codes to a broad class of $n$-dimensional Clifford hierarchy stabilizer codes. These codes correspond to the $(n+1)$D Dijkgraaf-Witten gauge theories with non-Abelian topological order. We construct transversal non-Clifford gates through automorphism symmetries represented by cup products. In 2D, we obtain the first transversal non-Clifford logical gates including T and CS for Clifford stabilizer codes, using the automorphism of the twisted $\mathbb{Z}_2^3$ gauge theory (equivalent to $\mathbb{D}_4$ topological order). We also combine it with the just-in-time decoder to fault-tolerantly prepare the logical T magic state in $O(d)$ rounds via code switching. In 3D, we construct a transversal logical $\sqrt{\text{T}}$ gate in a non-Clifford stabilizer code at the third level of the Clifford hierarchy, located on a tetrahedron corresponding to a twisted $\mathbb{Z}_2^4$ gauge theory. Furthermore, our constructions surpass the Bravyi-König bound by achieving the logical gates in the $(n+1)$-th level of Clifford hierarchy in $n$ spatial dimension.

Kobayashi, Ryohei [Institute for Advanced Study, P↗
DOE OSTI2022

Calculation of Dangerous Values for Radionuclides Considered by the IAEA Code of Conduct

The D-value or dangerous quantity system was designed by the International Commission for Radiological Protection for the determination of source protection categories that can be used to reduce the likelihood of accidents, the consequences of which could result in harm to individuals or costly or expensive cleanup. The process includes multiple scenarios for exposure and two different approaches to the evaluation of detriment. This document provides an example calculation using 137 Cs to walk through the complex process of determining its D-value in the hopes of making the process easily understandable.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2023

Report on FY 2023 Welding R&D at ORNL in support of ASME Alloy 709 Code Case Development

As part of the Alloy 709 (A709) ASME Code Case development effort under the Advanced Reactor Technologies (ART) Program, this work covers the development of the technical basis for weld fabrication and weld qualification of A709. This report summarizes the A709 welding research conducted at Oak Ridge National Laboratory (ORNL) in FY 2023. One of the focuses of the A709 welding research at ORNL is to relax the restriction of the P content in the A709 filler metal weld wire. In this report, we show the success in welding of high phosphorus (P) commercial A709 plates using A709 weld wires with relaxed P content of 80 wppm. A test weld was fabricated with the 80 wppm weld wire on the first commercial heat A709 base metal with 140 wppm P content using gas tungsten arc welding (GTAW). The test weld passed all weld qualification tests without issues. This study concludes that A709 matching filler metal with moderate amount of P of less than 80 wppm could be used to produce code qualified A709 welds with standard welding process without special treatment. It is recommended that the P level in the weld wire to be kept low to avoid hot cracking. Cross-weld specimens were machined from the two A709 production welds.. Cross welding creep rupture Code Case testing on these two production welds was initiated . The preliminary cross-weld creep tests results continue to show little or no creep strength reduction relative to the base metal.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Spectroscopic fingerprints of gapless type-II fracton phases

Fracton phases feature elementary excitations with fractionalized mobility and are exciting interest from multiple areas of theoretical physics. However, the most exotic type-II fracton phases, like the Haah codes, currently have no known experimental diagnostics. Here, we explain how type-II fracton phases with gapless gauge modes, such as the U(1) Haah code, may be identified experimentally. Our analysis makes use of the multipole gauge theory description of type-II fracton phases, which exhibits ultraviolet-infrared (UV-IR) mixing. We show that neutron scattering experiments on gapless type-II fracton phases should generically exhibit exotic pinch points in the structure factor, with distinctive anisotropic contours as a direct consequence of UV-IR mixing. Furthermore, this characteristic pinch point structure provides a clean diagnostic of type-II fracton phases. We also identify distinctive signatures of the (3 + 1)-D U(1) Haah code in the low-temperature specific heat.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Coherent wiggler radiation impedance at the storage ring cooler for the EIC project

The Electron Ion Collider project is presently under design at Brookhaven National Laboratory. One of the options how to achieve an electron-proton high-luminosity of 10 34 cm –2 s –1 range, is the storage ring cooler concept, which is based on employing a significant amount of the damping wigglers. One of the main concerns, in achieving the required beam parameters, is the collective effects, especially the coherent synchrotron radiation impedance produced by the damping wigglers and its effect on the longitudinal beam dynamics. Low energy of the electrons, E o =149.6 MeV, small vacuum chamber aperture, b=15 mm, small bending radius and a big number of poles make the coherent synchrotron radiation simulations for the damping wiggler with D. Zhou’s CSRZ code, pretty challenging. The obtained numerical results have been compare with a theoretical approach of Stupakov and Zhou. The strong narrow-band impedance, due to a presence of the periodic poles and the vacuum chamber, have been identified and classified. To suppress or detune the high-Q resonance peaks, a design of the damping wiggler with a varied period of length is presented and discussed.

43 PARTICLE ACCELERATORS↗