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 937 records · Page 52

A Systematic Approach to Light Load Calculation for DG Un-Intentional Islanding

With the influx of distributed generation (DG) penetration, power utilities have adopted active anti-islanding protection methods to prevent an unintentional islanding condition. One commonly used protection scheme is implemented through the use of a transfer trip from the upline connected substation to the DG facility and is required when the capacity of connected DGs exceed a certain threshold of local minimum load. The threshold established by standard IEEE 1547.2-2009 states that a 3:1 minimum load to generation ratio is acceptable to ensure the DG will not sustain an unintentional island. However, the term “minimum load” is not precisely defined in any standard, and in practice, a light load value is chosen by engineers through manual methods that can be inconsistent, not representative and change drastically based on various cases. To solve this problem, this paper proposes a new light load calculation method creating consistency and better accuracy; the method has been implemented in. NET C# application and is used within Dominion Energy.

Chen, Le↗

Simple Building Calculator

Whole building energy modeling is a powerful tool for analyzing energy use in buildings. For large buildings benefits of modeling services can easily be quantified due to the short payback period of the implemented measures. For small buildings, the upfront costs of developing robust energy models often deter building owners from investing in energy modeling services. To address this shortcoming, a ”Simple Building Calculator” (SBC) was developed. SBC uses pre-simulated results of a range of common measures over a wide range of efficiency inputs. It combines whole building simulation results with statistical modeling techniques to predict energy impact of measures very quickly. In this paper we present the modeling methodology used to develop the data supporting SBC.

Building Energy Simulation, commercial building, E↗

Electricity Bill Calculator Updates

A summary of updates to the electricity bill calculator in the most recent version of the System Advisor Model including changing the true-up month and new compensation options for metering and billing. Also included a preview of features in the upcoming release, such as billing demand for energy charges.

behind-the-meter↗

Euler-Rodrigues Parameters: A Quantum Circuit to Calculate Rigid-Body Rotations

The use of vectorial parameterization to create geometrical representations in computational models has a large number of applications. One particular application is the calculation of the 3D rotational motion of rigid bodies, that could be used for the spatial location estimation from objects. Provided the algebraic nature of this problem, it could benefit from Quantum Computing, in particular several vectors could be superposed to be transformed with a single operation, providing a quantum processing advantage. In this article, we propose an implementation of a Quantum Computing algorithm to compute Euler-Rodrigues Parameters to model rigid body rotations to transform arbitrary functions, rotating multiple vectors in superposition. We developed this algorithm using Qiskit, taking into account the limitations imposed by the current Noisy Intermediate Scale Quantum (NISQ) devices, such as the reduced number of qubits available and the limited coherence time.

Pelaez, Emilio↗

Calculating the QED correction to the hadronic vacuum polarisation on the lattice

Isospin-breaking corrections to the hadron vacuum polarization component of the anomalous magnetic moment of the muon are needed to ensure the theoretical precision of $g_\mu -2$ is below the experimental precision. We describe the status of our work calculating, using lattice QCD, the QED correction to the light and strange connected hadronic vacuum polarization in a Dashen scheme. We report results using physical $N_f=2+1+1$ HISQ ensembles at three lattice spacings and three heavier-than-light valence quark masses.

Ray, Gaurav↗

Cryogenics Python Coding Calculations for Nitrogen & Helium Loss

Helium and Nitrogen are two significant chemical elements that are used to advance cryogenic scientific breakthroughs. Helium is a colorless, odorless, and tasteless gas that becomes liquid at (−452 °F). Nitrogen is a gas at its standard conditions; It becomes solid once temperatures reach below (-346 °F), and it boils and becomes a gas once temperatures reach above (-320 °F). During this research project I used a web-based interactive computing platform called Jupyter Lab to modify various python programs that calculate, and display helium and nitrogen inventory lost from three cryo plants

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Hybrid power system control and operating strategy based on power system state vector calculation

Controlling a hybrid power system includes calculating a power system state vector based on energy demand and a stored data array including a matrix defined by a power system hardware configuration. The control further includes producing a power request based on the power system state vector, and varying a flow of energy amongst energy devices using drive linkages in the hybrid power system based on the power request. Related apparatus, control logic and controller structure is disclosed.

Guo, Fang↗

Status of the OECD/NEA Watts Bar unit 1 benchmark and calculations of local reactor core power data of the Zero Power physics tests

The paper aims to provide an update on the status of the Organization for Economic Cooperation and Development (OECD) / the Nuclear Energy Agency (NEA) Tennessee Valley Authority (TVA) Watts Bar 1 (WB1) Multi-Physics Multi-Cycle depletion benchmark and the results associated with Exercise 1. The benchmark relies on the set of benchmark progression problems developed by the Department of Energy (DOE) Consortium for the Advanced Simulation of Light Water Reactors (CASL) for the Virtual Environment for Reactor (VERA), which are based on real plant design and operational data. The OECD/NEA TVA WB1 benchmark is designed for validation of both traditional and novel high-fidelity multi-physics codes to analyze Pressurized Water Reactors (PWR) depletion cycles. The activities are conducted under the Expert Group on Reactor Systems Multi-Physics (EGMUP) at NEA/OECD. In this work, we analyze Exercise 1 of the benchmark: stand-alone Three-Dimensional (3D) neutronics at Start-up Zero Power Physics Test (ZPPT) at Hot Zero Power Conditions (HZP) and the resulting power maps. The code used to model the exercise is the continuous energy Monte Carlo code Serpent 2.1.31 along with the nuclear data library ENDF/B-VII.1. Serpent results are compared with the high-fidelity deterministic code MPACT, which is part of VERA. The paper presents three selected results from power calculations required output. The results compared are the normalized axially integrated radial core power maps, the normalized axial averaged core power shapes, and the normalized core hottest and coldest assemblies radial power maps. The Root Mean Square Deviation (RMSD) between Serpent and MPACT is 0.51 % for the axially integrated radial core power map, 1.41 % for the axial averaged core power shape, 2.42 % for the hottest assembly axial power shape, 0.88 % for the coldest assembly axial integrated hottest assembly radial pin power map, and 0.17 % for the axially integrated coldest assembly radial pin power map. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Extension of the high-resolution thermal-hydraulics code ESCOT to hexagonal core geometries for multi-physics calculations

The extension of the capabilities of the pin-level nuclear reactor core thermal-hydraulics (T/H) code ESCOT to analyze hexagonal fueled cores and its performance are presented. ESCOT is an accurate yet fast core thermal-hydraulics solution aiming at high-fidelity and high-resolution multi-physics core analysis in the framework of massively parallel computing platforms. Its algorithm solution is based on the four-equation drift-flux model for two-phase calculations, these are numerically solved by applying the Finite Volume Method (FVM) and the Semi-Implicit Method for Pressure-Linked Equation (SIMPLE)-like algorithm in a staggered grid system. Constitutive models such as turbulent mixing, pressure drop, and vapor generation are employed to simulate key phenomena in subchannel-scale analysis. ESCOT is parallelized by a double (radial and axial) domain decomposition that enables its highly parallelized execution. The coupling of the code with the neutronics whole core solver for hexagonal geometries nTRACER is described. The newly implemented ESCOT features are validated by comparing single assembly and full core steady state nTRACER-ESCOT solutions with nTRACER standalone internal one-dimensional T/H solver results. The validation problems are based on the VVER 440 and VVER 1000 cores. ESCOT results show differences within an acceptable range with respect to the simple 1D nTRACER built-in solver. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DFT Calculation of Thermal Expansion within Debye-Grüneisen Framework Made Simple and Accurate

The density functional theory (DFT) based Debye-Grüneisen approach is a computational method used to calculate the thermodynamic properties of materials by relating the volume dependence of the Debye frequency to the Grüneisen parameter. In this study, we propose a new model to describe the volume dependence of the Grüneisen parameter. This model offers a simple yet accurate approach to predict the thermal expansion of materials within the framework of the Debye-Grüneisen theory. Specifically, for pure elements or species used as a reference system, we calibrated the parameters to reproduce the experimental coefficient of thermal expansion, heat capacity, and entropy. For compounds with arbitrary compositions, the parameters are determined by taking the arithmetic or geometric average of those obtained from the pure elements or species based on their composition. We then compared the model prediction with experiments. This allowed for a comprehensive characterization of the thermal expansion behavior of various materials.

Wang, Yi↗

A COUPLED DETERMINISTIC TRANSPORT CALCULATION IN COMSOL USING PROPER ORTHOGONAL DECOMPOSITION

A coupled physics problem consisting of radiation transport and heat transfer was modeled using COMSOL Multiphysics to simulate a criticality accident. Reduced order models (ROM) were created to accelerate the radiation transport calculation through Proper Orthogonal Decomposition (POD). Additionally, on-the-fly neutron cross sections were generated from a non-linear function through a feed forward neural network. The Livelink for MATLAB module allowed for data transfer between an otherwise independent solver and in-house modules. A data driven approach to POD was developed to act as a surrogate model to couple with heat transfer, with the POD module written in MATLAB coupled to an isolated heat transfer solve in COMSOL. A coupled transient solution resulted in a good approximation of the flux and temperature, with the maximum error being 3E-2 and 2.5E-3 respectively.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Generating An Advanced Cross-section Library For HTGR Pebble Bed Depletion Calculations Using Reduced-Order Model Generation Techniques

For code development, Advanced Reactor Technologies - Gas Cooled Reactors Program (ART-GCR) rely on a collaboration with the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, but the cross sections generation and the methodology definition is part of this program area goals. Based on previous studies in FY23, the size of microscopic cross section libraries increases rapidly with the number of tabulations, requiring significant amount of memory and drastically slowing down the Griffin calculations when evaluating cross sections via the multivariate linear interpolation approach. Rising to these challenges, this work investigates constructing Reduced-order Models (ROMs) for the multi-group microscopic cross sections to accelerate the cross section evaluation in Griffin. A database of multigroup cross sections is first collected considering all possible parameters that a designer could change for optimization. Down-selection of the ROM techniques afterward shows Deep Neural Network (DNN) as the best candidate when jointly consider memory efficiency, predictive accuracy, computational cost, scalability, flexibility and ease of implementation of the algorithms in comparison to the multidimensional interpolation. This work develops a specific interface that enables the cross section predictions using pre-trained DNN models into Griffin leveraging the existing ROM capabilities. DNNs have been trained for all isotopes for use in Griffin. Preliminary Griffin testing shows that DNNs exhibit exceptional predictive accuracy and the use of DNNs provides orders of magnitude improvement in memory efficiency compared to conventional interpolation techniques. With such ROM techniques, it holds great promise to further increase the fidelity of the Pebble Bed Reactor (PBR) simulation by increasing the number of tabulations/state variables during cross section evaluation, while maintaining the computational cost affordable in Griffin.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

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 paper presents 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↗

Excited-state uncertainties in lattice-QCD calculations of hadron masses and scattering phase shifts

Lattice QCD has historically produced energy results interpretable as either estimates relying on implicit assumptions about asymptotic behavior or one-sided upper bounds. New Lanczos methods providing two-sided bounds with less-restrictive assumptions are introduced and quantified in a high-statistics calculation with unphysical quark masses. Two-sided bounds without spectral assumptions provide sub-percent constraints on the nucleon mass. Other bounds, which assume all states in a given energy window are resolved, provide meaningful two-sided constraints on nucleon-nucleon scattering phase shifts.

Detmold, William [MIT, Cambridge, CTP]↗

New Developments of MOOSE/FENIX Capabilities for Fusion Neutronics Calculations

This poster summarizes internship work to support new developments of MOOSE/FENIX Capabilities for fusion neutronics calculations. A new feature was added to allow coupling OpenMC models with thermomechanics models and was verified. Two coupled multiphysics models are demonstrated for a tokamak model and monoblock divertor model. Results show heating results, temperature distributions, tritium production and transport.

99 - GENERAL AND MISCELLANEOUS↗

New Developments of FENIX Capabilities for Fusion Neutronics Calculations

This presentation summarizes internship work to support new developments of FENIX Capabilities for fusion neutronics calculations. Two coupled multiphysics models are demonstrated for a tokamak model and monoblock divertor model. Results show heating results, temperature distributions, tritium production and transport. A new feature was added to allow coupling OpenMC models with thermomechanics models and was verified.

99 - GENERAL AND MISCELLANEOUS↗

High-pressure phase diagram of beryllium from ab initio free-energy calculations

In this report we use first-principles molecular dynamics simulations coupled with the thermodynamic integration method to study the hexagonal close-packed (hcp) to body-centered cubic (bcc) transition and melting of beryllium up to a pressure of 1600 GPa. We derive the melting line by equating solid and liquid Gibbs free energies and represent it by a Simon-Glatzel fit T m = 1564 K [1 + P/(15.6032 GPa)] 0.383 , which is in good agreement with previous two-phase simulations <6000 K. We also derive the hcp-bcc solid-solid phase boundary and show that the quasiharmonic approximation underestimates the stability of the hcp structure, predicting lower transition pressures between hcp and bcc phases. Our results are consistent with the stability regime predicted by the phonon quasiparticle method. We also predict that the hcp-bcc-liquid triple point is located at 164.7 GPa and 4314 K. In addition, we compute the shock Hugoniot curve and show that it is in good agreement with experiments, intersecting our derived melting curve at ~235 GPa and 4900 K. Finally, we make predictions for future ramp compression experiments. Starting with an isentropic compression of the liquid, we predict the path to intersect the melting line at low pressure and temperature, then to continue along the melting line over a large temperature interval of 7000 K as the sample remains in the mixed solid-liquid state before it enters the solid phase.

36 MATERIALS SCIENCE↗

Comparison of the canonical transformation and energy functional formalisms for ab initio calculations of self-localized polarons

In materials with strong electron-phonon (e-ph) interactions, charge carriers can distort the surrounding lattice and become trapped, forming self-localized (small) polarons. We recently developed an ab initio approach based on canonical transformations to efficiently compute the formation and energetics of small polarons A different approach based on a Landau-Pekar energy functional has been proposed in the recent literature. Here in this work, we analyze and compare these two methods in detail. We show that the small polaron energy is identical in the two formalisms when using the same polaron wave function. We also show that our canonical transformation formalism can predict polaron band structures and can properly treat zero- and finite-temperature lattice vibration effects, although at present using a fixed polaron wave function. Conversely, the energy functional approach can compute the polaron wave function, but as we show here, it neglects lattice vibrations and cannot address polaron self-localization and thermal band narrowing. Taken together, this work relates two different methods developed recently to study polarons from first-principles, highlighting their merits and shortcomings and discussing them both in a unified formalism.

36 MATERIALS SCIENCE↗