Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “calculator”

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 541 records · Page 30

COMPASS-U Global Heat Balance Calculations

COMPASS-U is a medium size, high magnetic field experimental tokamak (R = 0.9 m, B t = 5 T, and I p = 2 MA), built at the Institute of Plasma Physics, Czech Academy of Sciences (IPP-CAS). This global heat balance calculation for COMPASS-U was done at Princeton Plasma Physics Laboratory (PPPL). Based on our previous experience of building global thermal model for National Spherical Torus Experiment-Upgrade (NSTX-U) at PPPL, this 2-D global thermal model geometry represents a typical cross section of COMPASS-U machine. The model includes thermal radiation, conduction, and convection among components and also between the machine and outer environment. Helium gas heating/cooling was modeled with fluid element and surface element. This model was used to calculate component temperatures and heat distribution during heat up, cryogenic cool down, normal operation, and fault operation scenarios. Definition of 14 main thermal scenarios was provided by IPP-CAS. First nine are normal operation scenarios. Scenarios #10–#14 are fault scenarios. Results of thermal scenario #5 will be given and discussed in this article, including peak temperatures, temperature ratcheting, energy distribution, and required cooling power.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

New calculations for photosynthesis measurement systems: what’s the impact for physiologists and modelers?

Measurement of carbon dioxide (CO 2 ) and water vapor (H 2 O) exchanged by leaves has been tremendously important for advancing understanding of plant physiology. Typically, photosynthesis measurement systems estimate the concentration of CO 2 and H 2 O in air that has passed over a leaf inside a leaf chamber. Together with additional information such as flow rate, leaf area, and measurement of environmental variables, the calculations of von Caemmerer & Farquhar (1981), hereafter vCF1981, are used to estimate CO 2 and H 2 O fluxes between the leaf and atmosphere. Furthermore, the combination of the vCF1981 calculations and photosynthesis measurement systems has enabled significant advances in the understanding and model representation of leaf-level photosynthesis and respiration.

54 ENVIRONMENTAL SCIENCES↗

Development of a multiphase equation of state for gallium with experiments and ab initio free-energy calculations

Here, we present a five-phase equation of state (EOS) for elemental gallium (Ga) that is developed using both experimental data and new theoretical predictions. Four experimentally observed solid phases (Ga-I, Ga-II, Ga-III, and Ga-IV) and one liquid phase are included. To improve our understanding of the thermal behavior of Ga and its phase boundaries under compression, we have performed ab initio density functional theory (DFT) free-energy calculations for the Ga-III and liquid phases, which enables us to determine the melt temperature to be 2214 ± 100 K at 110 GPa, extending significantly beyond the existing experimental melt data, which are limited to 25 GPa. In order to best describe the electron-thermal contribution, which dominates the liquid free energy at high temperatures, we have carried out averaged-atom-in-jellium DFT calculations to cover the entire temperature and density range of the EOS. The resulting multiphase Ga EOS is able to accurately reproduce a diverse variety of data, including known phase boundaries, the principal Hugoniot, low-pressure liquid isobars, and diamond-anvil-cell isotherm measurements at high pressures. It agrees more closely with key experimentally measured properties than other Ga EOS models targeted for high-pressure applications.

Wu, Christine J. [Lawrence Livermore National Labo↗

Calculated fission-fragment mass yields and average total kinetic energies of heavy and superheavy nuclei

Fission-fragment mass and total-kinetic-energy (TKE) distributions following fission of even-even nuclides in the region \(74 \le Z \le 126\) and \(92 \le N \le 230\) , comprising 896 nuclides have been calculated using the Brownian shape-motion method. The emphasis is the region of superheavy nuclei. To show compatibility with earlier results the calculations are extended to include earlier studied regions. An island of asymmetric fission is obtained in the superheavy region, \(106\le Z\le 114\) and \(162\le N\le 176\) , where the heavy fragment is found to be close to \(^{208}\) Pb and the light fragment adjusts accordingly. Most experimentally observed \(\alpha \) -decay chains of superheavy nuclei with \(Z > 113 \) terminate by spontaneous fission in our predicted region of asymmetric fission. In these cases, the pronounced large asymmetry is accompanied by a low TKE value compatible with measurements.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Faster spectral density calculation using energy moments

Abstract Accurate predictions of inclusive scattering cross sections in the linear response regime require efficient and controllable methods to calculate the spectral density in a strongly-correlated many-body system. In this work we reformulate the recently proposed Gaussian Integral Transform technique in terms of Fourier moments of the system Hamiltonian which can be computed efficiently on a quantum computer. One of the main advantages of this framework is that it allows for an important reduction of the computational cost by exploiting previous knowledge about the energy moments of the spectral density. For a simple model of medium mass nucleus like $$^{40}$$ 40 Ca and target energy resolution of 1 MeV we find an expected speed-up of $$\approx 125$$ ≈ 125 times for the calculation of the giant dipole response and of $$\approx 50$$ ≈ 50 times for the simulation of quasi-elastic electron scattering at typical momentum transfers.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

eBPF Community ID Calculation

The code titled “eBPF Community ID Calculation” allows easy correlation of monitoring logs with existing network forensics logs on high-performance computing (HPC) systems. This code implements the Community ID Flow Hashing standard to calculate community ID Flow identifiers. These identifiers will be used to correlate host forensics logs with network forensics logs. The correlation provides valuable insight and the ability to easily pivot between disparate data sources. The computer code removes the libpcap dependency from the code base.

Maerz, StefanA↗

Guideline 36 Savings Calculator v0.1

This software package contains a simple calculator that can be used to estimate savings from implementing a retrofit of HVAC controls to ASHRAE Guideline 36 in a multi-zone office building. Implemented in Excel, the calculator has a user-interface front-end to configure parameters of the case study, such as location, building size, and existing control strategies, and a back-end containing results from a large number of pre-run annual energy simulations. The software package also contains the simulation models and associated pre- and post-processing scripts to generate these results.

Blum, David↗

Risk-Based Dispatch Bound Calculator

This calculator outputs time-varying lower and upper bounds (matrices), called dispatch bounds, for the controllable portions of the power injections at each node in a distribution network to guarantee that the voltage magnitudes stay within prescribed bounds with at least a given probability for any set of power injections that are within their dispatch bounds. The calculator is meant to be a tool for a central controller of a distribution network to set dispatch bounds for each node that allows each node to independently decide their actual power injection within their dispatch bound. The user can choose whether to assume the forecast errors of the uncontrollable power injections are gaussian or are worst-case (distributionally robust). The user can choose whether the dispatch bounds include the controllable reactive power or just the active power under a constant power factor assumption.

Comden, Joshua↗

Calculating topological properties of artificial graphene in B-field

SAND2025-03258O Calculating topological properties of artificial graphene in B-field is a user-friendly tool designed to analyze artificial graphene systems influenced by magnetic fields. It helps researchers understand the unique properties of these materials by calculating the local Chern marker, a key indicator of their behavior. With just one file, this software runs easily on any device with Matlab, making it accessible for scientists and engineers. It provides valuable insights into the electronic characteristics of artificial graphene, which supports advancements in material science and technology, paving the way for innovative applications in electronics and beyond. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

Spataru, Dan [Sandia National Lab. (SNL-CA), Liver↗

An initial investigation of accuracy required for the identification of small molecules in complex samples using quantum chemical calculated NMR chemical shifts

The majority of primary and secondary metabolites in nature have yet to be identified, representing a major challenge for metabolomics studies that currently require reference libraries from analyses of authentic compounds. Using currently available analytical methods, complete chemical characterization of metabolomes is infeasible for both technical and economic reasons. For example, unambiguous identification of metabolites is limited by the availability of authentic chemical standards, which, for the majority of molecules, do not exist. Computationally predicted or calculated data are a viable solution to expand the currently limited metabolite reference libraries, if such methods are shown to be sufficiently accurate. For example, determining nuclear magnetic resonance (NMR) spectroscopy spectra in silico has shown promise in the identification and delineation of metabolite structures. Many researchers have been taking advantage of density functional theory (DFT), a computationally inexpensive yet reputable method for the prediction of carbon and proton NMR spectra of metabolites. However, such methods are expected to have some error in predicted 13 >C and 1 H NMR spectra with respect to experimentally measured values. This leads us to the question–what accuracy is required in predicted 13 C and 1 H NMR chemical shifts for confident metabolite identification? Using the set of 11,716 small molecules found in the Human Metabolome Database (HMDB), we simulated both experimental and theoretical NMR chemical shift databases. We investigated the level of accuracy required for identification of metabolites in simulated pure and impure samples by matching predicted chemical shifts to experimental data. We found 90% or more of molecules in simulated pure samples can be successfully identified when errors of 1 H and 13 C chemical shifts in water are below 0.6 and 7.1 ppm, respectively, and below 0.5 and 4.6 ppm in chloroform solvation, respectively. In simulated complex mixtures, as the complexity of the mixture increased, greater accuracy of the calculated chemical shifts was required, as expected. However, if the number of molecules in the mixture is known, e.g., when NMR is combined with MS and sample complexity is low, the likelihood of confident molecular identification increased by 90%.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Decay Dose Shielding Analysis with Hybrid Unstructured Mesh/Constructive Solid Geometry Monte Carlo Calculation and ADVANTG Acceleration

The Oak Ridge National Laboratory (ORNL) Second Target Station (STS) neutron production facility is an accelerator driven pulsed neutron source that is currently being actively developed at ORNL. The neutrons are produced by proton-induced spallation reactions. A proton beam of 700 kW power is delivered to a spallation target in short, less than 1 µs long pulses, with 15 Hz repetition rate. The spallation target of ORNL STS is a rotating water-cooled tungsten target with tantalum cladding housed in a stainless-steel shroud. It is divided into 21 segments. These segments become highly activated due to spallation reactions or nuclei transmutation by the emitted neutrons. The radioactive nuclides continue to decay after ceasing operation. The decay dose rates generated from the target segments once they are removed from their operational location within the core vessel must be accurately quantified to determine the shielding configurations of remote handling tools and transport casks and to aid in planning maintenance events. To determine the shielding configurations needed for an activated target segment after ceasing operation, both the hybrid unstructured mesh (UM)/constructive solid geometry (CSG) approach that was previously utilized for STS analyses [1] and the ADVANTG code [2] were used. Even though the ADVANTG code does not include UM capability, the utilization of its advanced variance reduction technique was crucial to accelerate the extremely difficult final photon transport calculation in this analysis. This paper also describes the procedures taken to mitigate the convergence issues that occur when ADVANTG uses a source definition that does not match the source of the final Monte Carlo (MC) calculation. These convergence issues often occur because of inconsistencies between source and transport biasing parameters.

Ibrahim, Ahmad↗

Thermal Neutron Scattering Law Evaluation for Zirconium Carbide and Critical Mass Calculations

Zirconium carbide (ZrC) is a candidate material for use in advanced high temperature reactors, including space nuclear thermal propulsion applications. Thermal neutron scattering laws (TSLs) are generated for carbon bound in ZrC, C(ZrC), and zirconium bound in ZrC, Zr(ZrC), using ab initio lattice dynamics methods. These evaluations are to be submitted for inclusion in ENDF/B-VIII.1 and use the incoherent approximation for inelastic scattering as well as the new mixed elastic scattering treatment. The application of disordered alloy theory is introduced to appropriately capture the isotopic composition of Zr and C in the elastic scattering cross section. Localized higher energy vibrations in the C(ZrC) phonon density of states that are separated from lower energy modes result in quantized oscillations in the inelastic contributions to the TSL with a significant likelihood of large energy down-scattering and up-scattering interactions, where the latter increases in probability with temperature. The quanta of energy transfer during neutron thermalization is substantially greater than classically expected within the thermal neutron energy range. MC21 critical mass calculations of ZrC mixtures with high-enriched uranium demonstrate an impact of the TSLs when compared to free-gas treatment for 235 U concentrations less than 0.2 g/cm 3 . Additional MC21 critical mass calculations with homogenous mixtures of ZrC and reactor-grade graphite also demonstrate sensitivity to the ZrC TSL for thermal spectrum driven fission systems.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Stability analysis of Monte Carlo burnup calculation based on PWR fuel rods

The existing Monte Carlo burnup algorithms may give some unstable results when using multiple burnable cells during burnup calculations. A very small time step is required to maintain the stability of the system, which pose some limits for the realistic applications. In this paper, we firstly repeat the instability issue by using a Monte Carlo code based on a typical PWR fuel rod model. Then, a theoretical stability analysis was carried out by using the perturbation theory in order to get a better understanding about the property of system stability. We determine a perturbation matrix which relates some main influence factors, such as the time step and the height of fuel rod, it can help people estimate the stability of burnup calculation system and the results are quantitatively consistent with the Monte Carlo code. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

TRIPOLI-4 neutronics calculations for IAEA-CRP benchmark of CEFR start-up tests using new libraries JEFF-3.3 and ENDF/B-VIII

Sodium-cooled Fast Reactors (SFR) are one of the promising Generation IV fast reactors. The China Experimental Fast Reactor (CEFR) is a 65 MWth pool-type SFR with a high neutron leakage core using high enriched uranium oxide fuel. The CEFR start-up tests in 2010 consist of series of neutronics experiments. Essential experimental data are available from the 2018 IAEA-CEFR-CRP benchmark document and useful for the validation of neutron transport codes and nuclear data libraries. The TRIPOLI-4 Monte Carlo transport code is a general-purpose neutronics code using continuous-energy nuclear data libraries. It has a rich validation database covering different computational and experimental benchmark data sets to assure the accuracy and credibility of numerical studies. Previous TRIPOLI-4 SFR core physics calculations used mainly MOX fuels. SFR control rod worth studies with TRIPOLI-4 were mainly on big-size cores. To test the modeling capability of TRIPOLI-4 on fuel loading patterns and on variable control rods positions for SFRs and to check recent nuclear data libraries, it is interesting to investigate the CEFR start-up tests in this work, including core states from subcritical to supercritical, control rods and rod group worth, point kinetics parameters, and radial foil activation measurements. TRIPOLI-4 calculation results using new data libraries JEFF-3.3, ENDF/B-VIII, and those of ENDF/B-VII.1 for different CEFR core states are reported here. By means of different options of TRIPOLI-4, control-rod SAs reactivity worth, point kinetics parameters, and radial {sup 237}Np(n,f) fission rate distribution were successfully evaluated. Both ENDF/B-VIII and JEFF-3.3 nuclear data libraries provided reliable results with a difference of 170 +/- 11 pcm in k{sub eff}.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Analytic Insights into the Neutronic Characteristics of Neutron Moderators from MCNP Calculations

In the course of re-processing the ENDF/B VIII.0 thermal scattering law data, formulas and data checking insights were developed from MCNP calculations into the neutronic characteristics of moderators. Infinite medium test calculations were carried out in MCNP as part of the verification effort for the ACE formatted thermal scattering files. Edit results from MCNP were combined with derived results to formulate several equations which give analytic insights into the neutronic behavior of moderators.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Water chemistry data and calculations for Fracture Sustainability Tests

Measured cation, anion, and silica are presented as well as calculations evaluating mass balance performed for 4 fracture sustainability tests under EGS conditions: Stripa granite, Metamudstone, Rhyolitic Ashflow Tuff, and Silicified Rhyolitic Tuff. Calculations of charge balance are included for water samples in the 4 tests.

15 GEOTHERMAL ENERGY↗

Proton pulse charge calculation algorithm in Beam Power Limiting System at Spallation Neutron Source

A proton pulse charge calculation algorithm in the Beam Power Limiting System (BPLS) at the Spallation Neutron Source (SNS) was developed and implemented in an FPGA. The algorithm calculates one-minute running average of the pulse charges and issues a fault to the Personal Protection System (PPS) and the Machine Protection System (MPS) when a limit is reached. A bit-accurate model of the algorithm was first developed and tested in Matlab® and then implemented and simulated in VHDL using Vivado® design environment. Finally, the algorithm was verified on a µTCA-based hardware platform.

Bobrek, Miljko [ORNL] (ORCID:0000000332763451)↗

interflow: A Python package to organize, calculate, and visualize sectoral interdependency flow data

Many economic sectors rely on an uninterrupted “upstream” supply of a resource to conduct their primary functions, leaving them vulnerable to adverse effects should that resource flow be interrupted or compromised (OECD, 2017; U.S. EPA, 2010). Well-known examples of these relationships include water demand by the energy sector (e.g., thermoelectric cooling for nuclear generation) (Grubert & Sanders, 2018; Webber, 2017) and energy demand by the water sector (e.g., electricity required to treat or move water in the public water sector) (Congressional Research Service, 2017) though many others exist. Being able to calculate and document these interdependencies and evaluate where the greatest cross-sectoral intensities and flows exist can reveal opportunities to enhance the overall network. Despite the implications and potential impacts, however, these interconnections and flows have been historically complex to analyze and understand. The interflow package provides a flexible tool to organize, calculate, and visualize (using Sankey diagrams and other visualizations) sectoral interdependency flows for multiple subsectors and resources. This tool can help decision-makers, researchers, and other audiences more easily pull meaning from these interdependencies to reveal multi-faceted opportunities and risks. interflow can help investigate questions such as (1) which sectors have high cross-resource dependencies, (2) how does demand for a resource in various sectors compare across regions, and (3) where the sectoral and regional opportunities are for enhanced efficiency, security, and resiliency.

97 MATHEMATICS AND COMPUTING↗