Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “G 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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Detecting Magnetic Fields in Exoplanets with Spectropolarimetry of the Helium Line at 1083 nm

The magnetic fields of the solar system’s planets provide valuable insights into their interiors and can have dramatic consequences for the evolution of their atmospheres and interaction with the solar wind. However, we have little direct knowledge of magnetic fields in exoplanets. Here we present a method for detecting magnetic fields in the atmospheres of close-in exoplanets based on spectropolarimetric transit observations at the wavelength of the helium line at 1083 nm. This methodology has been successfully applied for exploring magnetic fields in solar coronal filaments. Strong absorption signatures (transit depths on the order of a few percent) in the 1083 nm line have recently been observed for several close-in exoplanets. We show that in the conditions in these escaping atmospheres, metastable helium atoms should be optically pumped by the starlight and, for field strengths more than a few × 10 -4 G, should align with the magnetic field. This results in linearly polarized absorption at 1083 nm that traces the field direction (the Hanle effect), which we explore by both analytic computation and the HAZEL numerical code. The linear polarization $\sqrt{{Q}^{2}+{U}^{2}}/I$ ranges from ~10 -3 in optimistic cases down to a few × 10 -5 for particularly unfavorable cases, with very weak dependence on field strength. The line-of-sight component of the field results in a slight circular polarization (the Zeeman effect), also reaching $V/I\sim \mathrm{few}\times {10}^{-5}({B}_{\parallel }/10\,{\rm{G}})$. We discuss the detectability of these signals with current (SPIRou) and future (Extremely Large Telescope) high-resolution infrared spectropolarimeters, and we briefly comment on possible sources of astrophysical contamination.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Reactor physics benchmark of Westinghouse PWR core design suite for high burnup/high enrichment fuel - Part 2: 2D and 3D core simulations

Westinghouse has performed a comprehensive set of code-to-code benchmark comparisons to corroborate application of its state-of-the-art NEXUS-based code suite to high burnup and high enrichment (HB/HE) fueled PWR core design and reload analysis. The benchmark consisted of a progression in complexity from modeling single 2D pin cells, through 2D lattices and finally to whole core 2D and then 3D scenarios. The accompanying paper (Part 1) presented the results for a variety of pin cell and lattices representative of high energy cores and established the ability of the MPACT code to accurately match the predictions of the Monte-Carlo code SERPENT2 such that the former may be used as a surrogate reference code in spatially complex whole core modeling applications where using Monte-Carlo methods is somewhat impractical. This paper focuses on comparisons at the whole core level using a representative high power density 3-loop Westinghouse NSSS core operating for 24 month cycles employing high enriched fuel with use of combined burnable absorbers to hold down excess reactivity. The evolution of the core reactivity and power distribution characteristics throughout the cycle are calculated using the Westinghouse NEXUS-based code suite and compared to the results using MPACT. The results of this whole core benchmark confirm the conclusions that the Westinghouse code suite accurately models HB/HE fuel designs, with very good prediction of core reactivity and power distribution throughout the cycle. The industry can be confident that the NEXUS code system can be deployed to support PWR reload safety and operational analyses for cores featuring HB/HE fuel and 24-month cycle operation. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

All Hands On Deck: Raising the Bar on Whole Building Performance-Based Code Compliance and Above-Code Programs

In 2019 a large-scale national effort was launched with support from U.S. Department of Energy, Pacific Northwest National Laboratory and Northwest Energy Efficiency Alliance (NEEA). The project aimed to develop a long-term vision and a roadmap for achieving a practical application of whole building performance-based code and above code program compliance in commercial buildings. The work was aided by engaging over 70 stakeholders across the country representing jurisdictions, above-code programs such as LEED and ENERGY STAR Multifamily Program, members of the ASHRAE Standard 90.1 and other Standard committees, IBPSA and software developers. The performance path requires whole building energy modeling and is gaining momentum. Many see it as the future of the commercial energy codes and the main pathway for achieving zero energy buildings. The paper describes the findings from a stakeholder survey funded by NEEA on the current state of performance-based path, its market penetration, trends, implementation practices and adoption challenges. It discusses the patterns revealed by experience of the states with decades-long history of performance-based compliance such as Florida and California and opportunities for taking advantage of the synergies between code and above-code programs. In addition, the paper provides an overview of the tools and resources developed as part of the effort including the compliance form for ASHRAE Standard 90.1 Section 11 and Appendix G, modeling submittal review manual, and modeler and reviewer qualification requirements. The tools and resources would help remove market barriers, deliver immediate benefits to jurisdictions and rating authorities, and achieve energy efficiency through whole-building performance path.

TSPR, Performance Based Complinace, Washington Sta↗
DOE OSTI2020

Leaf mass per area and leaf water content measurements from field survey in association with NEON AOP survey, East River, CO 2018

Leaf area, wet, and dry weights for leaf samples from sites where field sampling was conducted during the 2018 National Ecological Observatory Network (NEON) Airborne Observation Platform (AOP) imaging spectroscopy and lidar surveys in Gunnison County, Colorado. The sampling sites were located across East River, Washington Gulch, Slate River, and Coal Creek watersheds and contained a mixture of meadow, shrub, and tree sampling sites. This data package contains the leaf area and wet and dry weights for leaves from tree and shrub sampling sites, as well as samples for meadow species present within a given sampling area. Associated data packages contain species composition information at the plot level (DOI: 10.15485/1618130), as well as additional biogeochemical, microbial, and geophysical data as they become available. For full documentation, please see associated reference.This package contains zipped folders of all scans that were used for area calculations for both sampling sites (tree/shrub, site_scans.zip) as well as representative samples from a daily sampling area of meadow species that were found in the 1 m plots (meadow_area_scans.zip). It also contains the macro code that was used in ImageJ for processing of the area calculations. Finally, it includes excel files with the area (cm2), wet weight (g), dry weight (g), leaf mass per area (LMA, g/m2), and leaf water content (%) for each set of scans.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Importance of 3-D S{sub N} depletion in non-proliferation using BSOLVE

We present a single Pressurized Water Reactor (PWR) 3-D fuel rod design for depletion analysis using BSOLVE, our newly developed Runge-Kutta-Fehlberg based depletion code. BSOLVE is coupled with the deterministic 3-D S{sub N} particle transport code, PENTRAN, applied here with a 4-neutron energy group comparison to Continuous Energy (C/E) SERPENT2 Monte Carlo results. Differences are expected, as PENTRAN+BSOLVE retains full (multi-group) energy information for reactions, nuclide specific fission contributions, and energy dependent fission yields, using the latest available ENDF-BVIII data, important to retain accurate burned fuel inventories; SERPENT2 collapses burnup reactions to a single energy value. For depletion times up to ∼ 700 days and typical PWR power densities, relative differences between multigroup 3-D S{sub N} with full energy data and Monte Carlo one group burnup for trans-uranium nuclide concentrations and fission products are up to ∼20%. System eigenvalues are consistent, but with differences early and late in the cycle attributed to multigroup vs. C/E Monte Carlo cross sections. This work highlights the importance of low variance transport driven burnup for non-proliferation concerns, since plutonium quality varies significantly along axial lengths, and is more challenging to converge using Monte Carlo; details of depletion steps with spatial/zone dependent plutonium quality are provided. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2020

Metastable states from multinucleon excitations in 202 Tl and 203 Pb

The excited level structures of 202 Tl and 203 Pb, above the 7 + and 29/2 – isomers, respectively, have been studied. An isomer with I π = 20 + and T 1/2 = 215(10) μs has been established in 202 Tl, and the level scheme extended from I = 10 to 20ℏ with the placement of fifteen new transitions. In 203 Pb, the I π = 37/2 + state is established to be metastable, with T 1/2 = 2.5(3) ns. Levels in both nuclei arise from intrinsic excitations, with likely particle-hole character for the higher-lying states in 203 Pb. The 20 + isomer in 202 Tl is most likely associated with a $πh$ $^{–1}_{11/2}$ Ⓧ ν($i$ $^{–2}_{13/2}$, $f$ $^{–1}_{5/2}$) configuration, while the 37/2 + state in 203 Pb results from the excitation of five neutrons. Furthermore, calculations, using both an empirical approach and the oxbash code, have been performed to aid in the description of the excited level structure.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Sensor Distortion Effects in Photon Monte Carlo Simulations

We present a detailed method to simulating sensor distortions using a photon and electron Monte Carlo method. We use three-dimensional electrostatic simulations to parameterize the perturbed electric field profile for nonideal sensor details. We follow the conversion of simulated photons, and the subsequent response of the converted electrons to the electric field pattern. These nonideal sensor details can be implemented efficiently in a Monte Carlo approach. We demonstrate that the nonideal sensor distortions have a variety of observable consequence including the modification of the astrometric pattern; the distortion of the electron diffusion size and shape; and the distortion of flats. We show analytic validation of the diffusion physics, reproduce two kinds of edge distortion, and show qualitative validation of field-free regions, lithography errors, and fringing. We also demonstrate that there are two related effects of doping variation having different observable consequences. We show that field distortions from accumulated electrons lead to intensity-dependent point-spread functions and the sublinear variance in flats. The method is implemented in the Photon Simulator and the code is publicly available.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

The Athena++ Adaptive Mesh Refinement Framework: Design and Magnetohydrodynamic Solvers

The design and implementation of a new framework for adaptive mesh refinement calculations are described. It is intended primarily for applications in astrophysical fluid dynamics, but its flexible and modular design enables its use for a wide variety of physics. The framework works with both uniform and nonuniform grids in Cartesian and curvilinear coordinate systems. It adopts a dynamic execution model based on a simple design called a "task list" that improves parallel performance by overlapping communication and computation, simplifies the inclusion of a diverse range of physics, and even enables multiphysics models involving different physics in different regions of the calculation. We describe physics modules implemented in this framework for both nonrelativistic and relativistic magnetohydrodynamics (MHD). These modules adopt mature and robust algorithms originally developed for the Athena MHD code and incorporate new extensions: support for curvilinear coordinates, higher-order time integrators, more realistic physics such as a general equation of state, and diffusion terms that can be integrated with super-time-stepping algorithms. The modules show excellent performance and scaling, with well over 80% parallel efficiency on over half a million threads. The source code has been made publicly available.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

symPACK: A GPU-Capable Fan-Out Sparse Cholesky Solver

Sparse symmetric positive definite systems of equations are ubiquitous in scientific workloads and applications. Parallel sparse Cholesky factorization is the method of choice for solving such linear systems. Therefore, the development of parallel sparse Cholesky codes that can efficiently run on today’s large-scale heterogeneous distributed-memory platforms is of vital importance. Modern supercomputers offer nodes that contain a mix of CPUs and GPUs. To fully utilize the computing power of these nodes, scientific codes must be adapted to offload expensive computations to GPUs. We present symPACK, a GPU-capable parallel sparse Cholesky solver that uses one-sided communication primitives and remote procedure calls provided by the UPC++ library. We also utilize the UPC++ "memory kinds" feature to enable efficient communication of GPU-resident data. We show that on a number of large problems, symPACK outperforms comparable state-of-the-art GPU-capable Cholesky factorization codes by up to 14x on the NERSC Perlmutter supercomputer.

Bellavita, Julian↗
DOE OSTI2020

Code-to-code-to-experiment validation of LES-ALM wind farm simulators

The aim of this work is to present a detailed code-to-code comparison of two Large-Eddy Simulation (LES) solvers. Corresponding experimental measurements are used as a reference to validate the quality of the CFD simulations. The comparison highlights the effects of solver order on the solutions, and it tries to answer the question of whether a high order solver is necessary to capture the main characteristics of a wind farm. Both solvers were used on different grids to study their convergence behavior. While both solvers show a good match with experimental measurements, it appears that the low order solver is more accurate and substantially cheaper in terms of computational cost.

17 WIND ENERGY↗
DOE OSTI2022

Modeling of the TRIGA IPR-R1 research reactor with the Serpent2/RINNOVO Nodal core analysis package

The Serpent2/RINNOVO nodal core analysis code system, a dedicated tool for modeling research reactors, has the capability to accurately predict important core physics parameters involved in the safety of reactor operation, such as various reactivity coefficients, control rod and bank worths including the shutdown margin, power distributions as well as local neutron flux predictions at various core locations of high importance, e.g., at irradiation rigs. In this work, a hexagonal model of the unrodded initial core of the Brazilian IPR-R1 Mark I type TRIGA nuclear reactor has been created using the Serpent2/RINNOVO code system. The choice of employing a hexagonal core geometry representation has mainly been made to facilitate subsequent fuel shuffling operations and core follow calculations of this reactor. Numerical results in terms of the core eigenvalue and the assembly power distribution have then been compared against corresponding full core Serpent2 results to prove feasibility of using the Serpent2/RINNOVO code system for modeling small and highly heterogeneous TRIGA reactors. Overall, RINNOVO predicts the power distribution very accurately but the eigenvalue error still remains quite large. In light of being the very first evaluation of a TRIGA reactor with the RINNOVO nodal core simulator, these results are also considered to be very preliminary. Improved accuracy is expected by incorporating a proper methodology in the current code system for computing discontinuity factors for multi-assembly configurations. Furthermore, improved results are expected by increasing the number of energy groups used by RINNOVO in these core calculations.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Assessment of BISON capabilities for component-level prediction of tritium transport in fusion and fission applications

Here we assessed the hydrogen migration and redistribution model implemented in the BISON fuel performance code for component-level tritium transport applications in fusion and fission example scenarios. We developed BISON models of an ITER heat exchanger, light water reactor fuel cladding, and a fluoride salt cooled high temperature reactor heat exchanger. Reasonable agreement between reported values in reference studies and BISON predictions demonstrated the ability of the BISON models to predict tritium transport behavior through steel materials. Next, sensitivity and uncertainty analysis was used to understand key parameter sensitivities within the BISON model for evaluation of tritium migration and redistribution. . The sensitivity analysis showed that the diffusivity activation energy was the most important parameter and that the heat of transport was the least important parameter for tritium migration in steels. Overall, the capabilities of the BISON code for component level modeling of tritium transport are promising.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Liquid hydrogen storage system for heavy duty trucks: Configuration, performance, cost, and safety

In this work, we investigate the potential of liquid hydrogen storage (LH 2 ) on-board Class-8 heavy duty trucks to resolve many of the range, weight, volume, refueling time and cost issues associated with 350 or 700-bar compressed H 2 storage in Type-3 or Type-4 composite tanks. We present and discuss conceptual storage system configurations capable of supplying H 2 to fuel cells at 5-bar with or without on-board LH 2 pumps. Structural aspects of storing LH 2 in double walled, vacuum insulated, and low-pressure Type-1 tanks are investigated. Structural materials and insulation methods are discussed for service at cryogenic temperatures and mitigation of heat leak to prevent LH 2 boiloff. Failure modes of the liner and shell are identified and analyzed using the regulatory codes and detailed finite element (FE) methods. The conceptual systems are subjected to a Failure modes and effects analysis (FMEA) and a safety, codes, and standards (SCS) review to rank failures and identify safety gaps. The results indicate that the conceptual systems can reach 19.6% usable gravimetric capacity, 40.9 g-H 2 /L usable volumetric capacity and $174-183/kg-H 2 cost (2016 USD) when manufactured 100,000 systems annually.

08 HYDROGEN↗
DOE OSTI2023

Codes and standards for ceramic composite core materials for High Temperature Reactor applications

Fiber-reinforced ceramic matrix composites are attractive for high-temperature nuclear applications due to excellent thermal and mechanical properties as well as reasonable-to-outstanding radiation resistance. Over the past 20 years, the use of ceramic matrix composite applications expanded to many commercial non-nuclear industries as fabrication and application of the technologies mature. The ASME Boiler Pressure Vessel Code, under Section III Division 5, provides the design and construction rules for High Temperature Reactor components. It published the first rules for ceramic matrix composites to be used for reactor core components. The rules lay out the quality requirements together with the design and materials criteria for the use and application of silicon carbide- and carbon-based matrix material technologies. As with the established graphite rules, the ceramic composite material rules are structured in Subsection HH (from Section III), that addresses the criteria for class SN nonmetallic core components. The code rules rely heavily on the development and publication of standards for composite specification, classification, and testing of mechanical, thermal, and other properties. These test methods are developed in ASTM Committee C28 on Advanced Ceramics, with a current focus on ceramic composite tubes. This article describes the detail of the composites code, the design methodology and similarities to the graphite code, the guidance for the development of specifications for ceramic composites (for nuclear applications) including recent standard developments, and it mentions the next steps to support licensing aspects by validating the code with benchmarking data.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

RASPA3

RASPA3, a molecular simulation code for computing adsorption and diffusion in nanoporous materials and thermodynamic and transport properties of fluids. It implements force field based classical Monte Carlo/molecular dynamics in various ensembles. RASPA3 is rewritten from the ground up in C++23 with speed and code readability in mind. Transition-matrix Monte Carlo is added to compute the density of states and free energies. The Monte Carlo code for rigid molecules is based on quaternions, and the atomic positions needed in the energy evaluation are recreated from the center of mass position and quaternion orientation. The expanded ensemble methodology for fractional molecules, with a scaling parameter λ between 0 and 1, now also keeps track of analytic expressions of dU/dλ, allowing independent verification of the chemical potential using thermodynamic integration. The source code is freely available under the MIT license on GitHub.

Dubbeldam, David↗
DOE OSTI2024

Cross sections for the formation of Rb84m,g, Rb83, and Rb82m in Sr86(d,x) reactions up to deuteron energies of 49 MeV: Competition between α-particle and multinucleon emission processes

Cross sections of Sr86(d,x) reactions leading to the products Rb84m,g, Rb83, and Rb82m were measured by the stacked-sample activation technique up to deuteron energies of 49 MeV. Nuclear model calculations were performed using the codes talys and empire, which combine the statistical, precompound, and direct interaction components. In all cases, the empire results were much higher than the talys calculation. Fairly good agreement was obtained between measured data and the talys calculation after some optimization of the input model parameters. Insight into competition between α-particle and multinucleon emission in the Y88 compound-nucleus system was also gained.

59 ≤ A ≤ 89↗
DOE OSTI2020

Omnibus User Manual

This manual provides instructions for using Omnibus to run the Exnihilo code suite, which contains the Denovo and Shift transport solvers.

97 MATHEMATICS AND COMPUTING↗