Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Calculation”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

An Inverse Heat Conduction Algorithm Used to Calculate the Temperatures on the Inner and Outer Cylindrical Surfaces of an HMX-based PBX Explosive Annulus

In this work, a new Inverse Heat Conduction (IHC) algorithm is applied to estimate the surface temperatures at twelve locations on the inner and outer cylindrical boundaries of an HMX-based Plastic Bonded Explosive (PBX) annulus. This IHC algorithm was developed in references using a set of Direct Heat Conduction (DHC) solutions and a temperature correction method. The DHC solutions were calculated using a Galerkin based finite element (FE) method. This HMX based PBX annulus was used in the Large Scale Annular Cookoff (LSAC) experiment, Shot 5. The reason Shot 5 was chosen as a prototype mathematical model for this study is that the temperature was measured at eighteen locations in the midplane of the HMX-based PBX annulus. In addition, this annulus underwent an experimental thermal ignition and a deflagration that caused a thermal explosion and the disassembly of the experiment. The objective of this study is to describe how the application of the temperature correction algorithm produced the convergence of the DHC solutions to the measured temperatures at twelve internal locations in the midplane of the HMX-based PBX annulus.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Modeling graphene sheet growth and dynamical matrix calculations using molecular dynamics

Molecular dynamics (MD) has been an incredibly useful tool to model physical processes that were synthesized experimentally but not fully understood. MD, through the use of semi-empirical inter-atomic potentials, has allowed understanding of different physical processes in materials science. Yet as well as providing useful insights into materials science, molecular dynamics has a wider range of usability. In this report, I will be detailing how MD can be used to study graphene formation from a carbon liquid which requires high temperatures and pressures. Beyond this, I will describe the usefulness of MD for understanding the physics for phonon transport quantum sensors. To do this, MD was employed to determine the dynamical matrix by treating atoms as coupled oscillators. An accurate understanding of the dynamical matrix of a system is required to calculate the non-equilibrium Green’s function used to describe the phonon transport within phonon wave-guides. I found that, across multiple pressures and temperatures, randomly placed carbon atoms will show evidence of pent-first formation with semi-empirical models. Density functional theory (DFT), on the other hand, was too computationally expensive to use for full scale MD simulations, but we have the possibility of training a machine learned interatomic potential to approximate DFT for carbon in the environments being studied for pent-first graphene sheet formation.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Calculating the Effects of Solids Input and Removal as a Temperature Control in the Advanced Scale Up Reactor Experiment (ASURE) Facility at NETL Using Aspen

The Advanced Scale Up Reactor Experiment (ASURE) facility at NETL is being designed to be a fuel-flexible multi-purpose reactor that can be used for pyrolysis/gasification or evaluation of other high pressure “circulating fluidized bed” (CFB) chemical processes. The initial system design calculations for pyrolysis/gasification are presented in this work showing the expected performance of the ASURE facility when used as a biomass conversion reactor. Several other areas of application include gasification of any carbonaceous fuel including biomass, coal, plastics, and other waste materials. The reactor can therefore be used to produce SYNGAS of various compositions and hydrogen as well as other high value chemicals resulting from a typical tuned gasification process. This paper discusses an ASPEN model of the facility, focusing on the riser of the CFB reactor and the solids recirculation loop. The ASPEN model divides the riser into two sections. A bottom section which receives ash, char and sand which have been recirculated from a return loop. In this section an inert fluidization gas, (N2 or CO2), is introduced which acts as the primary mover of the solids through the system. The bottom section is equipped with a restricted air feed so that the recirculated char can be partially oxidized. This oxidation process along with the inventory of recirculating sand are used to effectively control the temperature in the following two chemical conversion sections of the reactor which are the pyrolysis zone followed by a tar cracking zone. Fresh fuel is added to the pyrolysis zone and undergoes drying and devolatilization. The products ash, char, volatile matter, and water vapor exit the pyrolysis zone and enter the reaction block for tar cracking. Steam and CO2 gasification reactions will be incorporated into the tar cracking zone, however at the design operating temperature, conversion from these reactions is expected to be essentially zero. The unit when completed in 2026 will test mixtures of biomass, plastics, and waste coal. This presentation discusses the basic ASPEN engineering design model for this project and provides preliminary sensitivity studies to determine how the various parts of the reactor will perform.

ASSURE↗
DOE OSTI2025

Coupling between collective modes in the deformed 98 Zr nucleus: Insights from consistent HFB+QRPA calculations with the Gogny interaction

The Zirconium isotopes exhibit structural properties that present multiple challenges to nuclear theory. Investigations of the coupling present within isoscalar modes and within isovector modes are scarce but important for advancing our understanding of the microscopic picture of nuclei. To explore some of these underlying coupling features, and to test the predictive power of a state-of the-art nuclear structure approach, we provide a detailed analysis of the properties of 90,96,98Zr. This region includes a benchmarking case and offers insights into nuclear deformation phenomena. To investigate the coupling between collective modes in deformed nuclei, we focused our analysis on the ground and excited-state properties of these isotopes, employing a consistent approach with the axially-symmetric deformed Hartree-Fock-Bogoliubov (HFB) and the Quasiparticle Random Phase Approximation (QRPA) framework, both using the Gogny D1M force. This approach effectively describes both low-lying and giant-resonance states. We devoted special attention to the deformed 98Zr nucleus, where we confirm the existence of coupling between monopole and quadrupole excitations through the K π = 0 + QRPA components and demonstrate an analogous dipole-octupole coupling through the K π = 0 − and K π = 1 − components. Intrinsic transition densities and associ ated radial projections illustrate the coupling. Our work complements and extends earlier studies carried out using density-functional-based methods and notably, we included the complete Coulomb interaction also in the pairing fields, i.e. we treat terms exactly that are approximated in typical calculations that use the Gogny D1 and D2 interaction families.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

FSEN Reaction Rate Calculations in MCNP [Slides]

FSEN provides an opportunity to quickly calculate sensitivities to orthogonal measurements during integral experiments. Toy problem has been deployed to expand on verification of FSEN for reaction rate ratios. While there is decent agreement, further investigation must be done on multiplication’s impact on sensitivity vector.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

WHOLESCALE - Water & Hole Observations Leverage Effective Stress Calculations And Lessen Expenses (Final Technical Report 2020 - 2024)

The WHOLESCALE acronym stands for Water & Hole Observations Leverage Effective Stress Calculations and Lessen Expenses. The goal of the WHOLESCALE project is to simulate the spatial distribution and temporal evolution of stress in the geothermal system at San Emidio in Nevada, United States. To reach this goal, the WHOLESCALE team has developed a methodology to incorporate and interpret data from four methods of measurement into a multi-physics model that couples thermal, hydrological, and mechanical (T H-M) processes. The WHOLESCALE team has applied this methodology at the San Emidio geothermal field, located ~100 km north of Reno, Nevada in the northwestern Basin and Range province. The WHOLESCALE team includes 30 individuals working at two universities, two national laboratories, and one industry partner. Two master-degree students and five post-doctoral researchers have gained professional experience and earned partial financial support via the WHOLESCALE project. The WHOLESCALE team has taken advantage of the perturbations created by changes in pumping operations during planned shutdowns in 2016, 2021, and 2022 to infer temporal changes in the state of stress in the geothermal system at San Emidio, Nevada, U.S. The WHOLESCALE results support the working hypothesis that increasing pore-fluid pressure reduces the effective normal stress acting across fault zones. During normal operations, pumping in deep production wells decreases fluid pressures and thus increases the effective normal stresses on faults, reducing microseismicity. During planned shutdowns, the cessation of production increases pore-fluid pressure and reduces effective normal stress. The WHOLESCALE products generated during the 4-year period between 2020 and 2024 include: three articles published in the open-access, peer-reviewed scientific literature, two master’s theses, 20 presentations or papers at scientific conferences, and 17 data sets available on public repositories. The WHOLESCALE project has been completed in two phases that included three performance periods separated by two Go/No-go Stage Gate Reviews. Tasks were classified by data type (i.e., Geologic Structure, Borehole, Geodesy, Hydrology, Seismology, and Modeling). The first phase of the project started July 31, 2020 and included ongoing project coordination (Task 1), a project kickoff (Task 2), analysis of existing data (Task 3), development of the initial stress model & deployment design (Task 4), and Go/No-go Decision Point #1 (Task 5). Phase II began with implementing the 2022 deployment (Task 6), followed by Go/No-go Decision Point #2 (Task 7) The remainder of Phase II consisted of analyzing data collected during deployment (Task 8), calibration of the stress model on all observations (Task 9), and the Final Review (August 23, 2024) & Reporting (Task 10).

15 GEOTHERMAL ENERGY↗
DOE OSTI2025

Expansion and Update of the R&D Feedstock Carbon Intensity Calculator (FD-CIC) for R&D GREET 2025

The R&D Feedstock Carbon Intensity Calculator (R&D FD-CIC) is a transparent and easy-touse tool to estimate feedstock-specific carbon intensity (CI) of biofuel feedstocks, and to examine potential CI variations of different farming practices to grow agricultural feedstocks for biofuel production. The objective of the tool is to assess potential effects of different farming practices with underlined factors to reduce uncertainties in such assessment and help R&D activities. The original work and early versions of R&D FD-CIC were supported by Advanced Research Projects Agency–Energy (ARPA-E) of the Department of Energy, with recent support by the Bioenergy Technologies Office of the Department of Energy. The tool uses life cycle inventory (LCI) data of key farming inputs from the R&D Greenhouse gases, Regulated Emissions, and Energy use in Technologies (R&D GREET) model. The system boundary of R&D FD-CIC covers the cradle-to-farm-gate activities, including upstream emissions related to farming input manufacturing and feedstock production.

09 BIOMASS FUELS↗
DOE OSTI2025

As-Run Neutronic Calculations of HFIR Mini-Fuel Experiments: Guidelines for Analysts

This report provides guidelines for analysts performing radiation transport calculations for mini-fuel experiments irradiated in the High Flux Isotope Reactor (HFIR). Procedures, practical recommendations, and modeling conventions which are essential to ensuring consistency and high fidelity in as-run simulations have been compiled. The content is intended to be a portion of the official training material for new staff in this area, offering both technical references and practical guidance derived from previous experience. This report is a living document and will be periodically updated to incorporate new developments, methodological improvements, and lessons learned from ongoing and future irradiation campaigns.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Comparison Calculations of Subcritical Limits in ANSI/ANS-8.1 [Slides]

The calculations presented in this paper were performed with ENDF/B-VIII.0 and SCALE 6.3.1 to provide users with perspective on what the current k eff values are for the SCLs listed in ANS-8.1. There are multiple SCLs that exceed the reference value of 0.98. This may indicate that these values need to be adjusted by the WG for the next iteration of the standard.

ENDF↗
DOE OSTI2024

Recommendations for an Applicant to Calculate Activity Data for Greenhouse Gases Estimates

In 2009, the U.S. Nuclear Regulatory Commission (NRC) directed the NRC staff to address climate change issues and consider the impacts of the emissions of carbon dioxide (CO 2 ) and other greenhouse gases (GHGs) in its environmental reviews for major licensing actions (NRC 2009b). To implement this direction from the Commission, the staff issued guidance in 2011 and updated guidance in 2014 in Attachment 1 to Interim Staff Guidance COL/ESP-ISG-026 (NRC 2011; NRC 2014). This guidance provides a simpler method than the method described in RG 4.2 Rev. 3, that an applicant can use to meet the plant parameter envelope (PPE) value from the Generic Environmental Impact Statement for Licensing of New Nuclear Reactors (NR GEIS). NRC staff estimated the 97-year lifecycle GHG emissions from a reference 1000 megawatt electrical (MWe) light-water reactor (LWR) for various activities associated with construction, operation (including uranium fuel cycle), and decommissioning of nuclear power plants and presented the results in Appendix H of the NR GEIS. Appendix H of the NR GEIS includes estimates of direct emissions from construction equipment and emergency diesel engines in a nuclear facility and indirect emissions from workforce vehicular traffic, fuel transportation and the uranium fuel cycle. The NR GEIS Section 3.3 extended the estimates in Appendix H for the installation of two 1000 MWe nuclear reactors on the same site. Scaling factors were used to extrapolate the GHG emissions of a reference 1000 MWe reactor to a two-unit nuclear reactor plant (each reactor unit generating 1000 MWe). GHG emission estimates for building, operation, decommissioning and safe storage (SAFSTOR) for a two-unit nuclear reactor plant would be based on the plant’s physical size, and therefore estimates for these source categories were assumed to be twice the value of the reference 1000 MWe reactor. However, GHG emissions from the fuel cycle (including fuel transportation) were scaled upward by a factor of 3, based on plant efficiencies greater than the 80 percent assumption in Appendix H. Table 1 below shows the PPE emissions for two 1000 MWe nuclear reactors as provided in NR GEIS. The total GHG emissions for two 1000 MWe reactors were calculated as 2,534,000 metric tons (MT) of CO 2 equivalent (CO 2 (e)) based on a 97 year GHG life cycle period. The GHG emissions lifetime of 97 years for a reference nuclear reactor includes a 7-year building phase, 40 years of operation, 10 years of active decommissioning, and 40 years of SAFSTOR operations (NRC 2024). Construction equipment and vehicular traffic from workers commute would contribute to the GHG emissions during a 7-year building phase. Uranium fuel cycle, vehicular traffic, fuel and waste transportation, and testing of standby diesel generators would contribute to GHG emissions during the 40-year operations phase. While NRC’s regulations allow up to 60 years of reactor facility decommissioning, Appendix H estimated that most of the GHGs would occur over an estimated 10-year period during which to the licensee would engage in significant demolition and earth-moving activities, as discussed in Supplement 1 to NUREG-0586 (NRC 2002). Vehicular traffic by the workforce during a 40-year SAFSTOR period would additionally contribute GHG emissions. The carbon footprint for a 40-year SAFSTOR period was separately analyzed from the decommissioning activities as provided in Table YYYY-2 of the staff issued guidance in 2011 (NRC 2011).

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Design and calculations for HWR Strongback

The PIP-II Superconductive Linac is a linear accelerator that provides the first stage of beam acceleration. The velocity increase is achieved through a sequence of different cryomodules, where all the cavities are aligned to form a continuous string. The cryomodule team is currently developing a spare for the HWR cryomodule. This spare will not be an exact replica of the existing one, but rather will be based on the SSR f650 cryomodule design created for PIP-II . Each cavity , magnet , and coupler form a unit, with a total of eight units per string. These components, together with additional equipment (such as piping , the thermal shield , and beamline interconnections ), are supported by the Strongback . The aim of this project is to develop a brand-new Strongback design for HWR cavities, and to perform the related analyses and calculations in order to meet the required specifications.

43 PARTICLE ACCELERATORS↗
DOE OSTI2026

Calculations Document

Figure 1 shows a simplified model schematic of the Tetracage target design. In this figure the green lines represent the Tetracage yarns, which consist of two parallel yarns on the top and two parallel yarns on the bottom with the pairs in a cross hatched pattern. For simplicity, as we are only considering vertical forces, we will perform calculations without consideration for the two pairs of yarns being perpendicular to each other. Also worth noting is the model shows the yarns extending directly to the capsule then closely hugging it – in reality, the yarns behave as you would expect where they hug the capsule only as much as the geometry of their deflection allows, similarly to if you held a tennis ball with a pair of parallel strings anchored at opposing ends. The dashed black lines show the path the yarns take straight across from slot glue point to slot glue point before the target is closed and the capsule deforms them. The variable r corresponds to the capsule outer radius. The variable L corresponds to the hohlraum inner radius. The variable a corresponds to the yarn offset from the center of the diagnostic band. The variable d corresponds to the wire deflection vertically.

42 ENGINEERING↗
DOE OSTI2025

An Educational Guide for 2D Stellar Structure Calculations of Rapidly Rotating Stars using the ESTER code

The Evolution STEllaire en Rotation (ESTER) code is the first 2D stellar structure code to be made open-source and freely available to the astronomy and astrophysics community. An important and novel advancement of this code is that it can reproduce the distorted shape and observable signatures (e.g., gravity darkening) of rapidly rotating stars. ESTER also calculates the steady-state large-scale flows within the star, namely their differential rotation and associated meridional circulation. In this report, we explore and document the physics implemented within version 1.1.0rc2 of the ESTER code, in a way that complements published descriptions. We illustrate this physics by plotting how stellar structure parameters vary through stellar interiors at a range of latitudes and at different angular velocities. We investigate how the thin convective envelopes of intermediate mass stars vary with latitude when rapidly rotating, becoming deeper and thicker near the equator. Simple comparisons of ESTER model predictions (e.g., central temperature and density, luminosity) with the output from the Modules for Experiments in Stellar Astrophysics (MESA) code [Paxton et al., 2010] shows generally good agreement. Additional comparisons provide important benchmarking and verification for ESTER as a comparatively young code. Finally, we provide a guide for installing and running the code on our local university cluster, aimed at helping students to begin work.

79 ASTRONOMY AND ASTROPHYSICS↗