Search NASA⌕ Search

SEARCH · Search NASA

Results for “Numerical performance”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Nonlinear power flow control for networked AC/DC microgrids

A method for designing feedforward and feedback controllers for integration of stochastic sources and loads into a nonlinear networked AC/DC microgrid system is provided. A reduced order model for general networked AC/DC microgrid systems is suitable for HSSPFC control design. A simple feedforward steady state solution is utilized for the feedforward controls block. Feedback control laws are provided for the energy storage systems. A HSSPFC controller design is implemented that incorporates energy storage systems that provides static and dynamic stability conditions for both the DC random stochastic input side and the AC random stochastic load side. Transient performance was investigated for the feedforward/feedback control case. Numerical simulations were performed and provided power and energy storage profile requirements for the networked AC/DC microgrid system overall performance. The HSSPFC design can be implemented in the Matlab/Simulink environment that is compatible with real time simulation/controllers.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

Modeling and Simulation of Air-Source CO2 Heat Pump Water Heater

Carbon dioxide (CO2) has been widely used as working fluid for the vapor-compression refrigeration systems in large marine device. Due to the potential energy efficiency and the favorable environmental properties of CO2 as a working fluid, CO2 heat pump water heater (HPWH) systems are regarded a promising technology for centralized domestic hot water (DHW) heating in residential and commercial buildings. However, there is still at the early stage of appropriately optimizing and improving the energy performance of CO2 HPWH. This requires CO2 HPWH simulation tools capable of capturing the accurate impact of the emerging compressor, throttle device, and heat exchanger technology on CO2 heat transfer and energy efficiency. In this study, high efficiency components (compressors, pumps, fans, heat exchangers) were identified and applied to the state-of-art CO2 HPWH designs and analyzed their performance by using numerical simulation. This was done by simulating the performance of CO2 HPWH using ACMODEL design model combined with the component models developed at Oak Ridge National Laboratory (ORNL) for orifice tube, map-based compressor, and tube-in-tube gas cooler. ACMODEL is an equipment design model for CO2-based air conditioners and heat pumps developed by Purdue University to account for the details of each component. The simulated CO2 HPWH performance was then compared with the heat pump water heater using conventional refrigerants.

Gao, Zhiming↗
DOE OSTI2024

Single-shot Quantum Signal Processing Interferometry

Quantum systems of infinite dimension, such as bosonic oscillators, provide vast resources for quantum sensing. Yet, a general theory on how to manipulate such bosonic modes for sensing beyond parameter estimation is unknown. We present a general algorithmic framework, quantum signal processing interferometry (QSPI), for quantum sensing at the fundamental limits of quantum mechanics by generalizing Ramsey-type interferometry. Our QSPI sensing protocol relies on performing nonlinear polynomial transformations on the oscillator's quadrature operators by generalizing quantum signal processing (QSP) from qubits to hybrid qubit-oscillator systems. We use our QSPI sensing framework to make efficient binary decisions on a displacement channel in the single-shot limit. Theoretical analysis suggests the sensing accuracy, given a single-shot qubit measurement, scales inversely with the sensing time or circuit depth of the algorithm. We further concatenate a series of such binary decisions to perform parameter estimation in a bit-by-bit fashion. Numerical simulations are performed to support these statements. Our QSPI protocol offers a unified framework for quantum sensing using continuous-variable bosonic systems beyond parameter estimation and establishes a promising avenue toward efficient and scalable quantum control and quantum sensing schemes beyond the NISQ era.

Physics↗
DOE OSTI2023

A hybrid analytical and numerical model for cross-over and performance decay in a unit cell vanadium redox flow battery

Developing an accurate and efficient model for cross-over is critical for improving the long-term performance of redox flow batteries (RFBs). Here, this work presents a hybrid analytical and numerical model that combines a two-dimensional analytical solution to the active species, a one-dimensional analytical model for cross-over mechanisms, and a zero-dimensional numerical model for outlet concentrations of reactants. By comparing with experiment over 41 cycles (ca. 144 h), the model reported a mean voltage difference of 0.0089 V, mean time difference of 14 s per cycle (of 3.5 h), maximum relative difference for capacity and energy of 1.34% and 1.63%, respectively. The predicted mean concentrations for V 2+ , V 3+ , and VO$_{2}^{+}$ in membrane are 65%~77% of measured values from the literature. Upon validation, the model reproduced behaviors in electrolyte imbalance similar to those observed in experiments and numerical models, and revealed the control of cross-over, self-discharge, stoichiometry of side reactions, and Coulombic efficiency on electrolyte imbalance. The model also demonstrates excellent computational efficiency for simulating 41 cycles (around 37,000 points) within 3~4 s. The demonstrated efficiency and accuracy for predicting cross-over and its impacts on voltage, capacity&energy decay, membrane concentrations, and electrolyte imbalance makes it a reliable tool for optimizing RFBs’ long-term performance.

25 ENERGY STORAGE↗
DOE OSTI2022

A Variational Method for Analyzing Vortex Flows in Radar-Scanned Tornadic Mesocyclones. Part III: Sensitivities to Vortex Center Location Errors

Abstract When the vortex center location is estimated from a radar-scanned tornadic mesocyclone, the estimated location is not error-free. This raises an important issue concerning the sensitivities of analyzed vortex flow (VF) fields by the VF-Var (formulated in Part I of this paper series and tested in Part II) to vortex center location errors, denoted by Δ x c . Numerical experiments are performed to address this issue with the following findings: The increase of |Δ x c | from zero to a half of vortex core radius causes large analysis error increases in the vortex core but the increased analysis errors decrease rapidly away from the vortex core especially for dual-Doppler analyses. The increased horizontal-velocity errors in the vortex core are mainly in the Δ x c -normal component, because this component varies much more rapidly than the other component along the Δ x c direction in the vortex core. The vertical variations of Δ x c distort the vertical correlation structure of Δ x c -dislocated VF-dependent background error covariance, which can increase the analysis errors in the vortex core. The dual-Doppler analyses have adequate accuracies outside the vortex core even when |Δ x c | increases to a half of vortex core radius, while single-Doppler analyses can also have adequate accuracies outside the vortex core mainly for the single-Doppler-observed velocity component. The sensitivities to Δ x c are largely unaffected by the vortex slanting. The above findings are important and useful for assessing the accuracies of analyzed VFs for real radar-observed tornadic mesocyclones. Significance Statement When the vortex center location is estimated from a radar-scanned tornadic mesocyclone, the estimated location is not error-free. This raises an issue concerning the sensitivity of analyzed vortex flow (VF) by the VF-Var (formulated in Part I of this paper series and tested with simulated radar observations in Part II) to vortex center location error. This issue and its required investigations are very important for the VF-Var to be applied to real radar-observed tornadic mesocyclones, especially in an operational setting with the WSR-88Ds. Numerical experiments are performed to address this issue. The findings from these experiments are important and useful for assessing the accuracies of VF-Var analyzed VF fields for real radar-observed tornadic mesocyclones.

Xu, Qin↗
DOE OSTI2022

Modeling Astrophysical Explosions and Nucleosynthesis (Final Technical Report)

The goal of this project was to further our theoretical understanding of the extreme physics and heavy element nucleosynthesis occurring in stellar explosions, in particular supernovae and neutron star mergers. To study these systems, we further developed and applied numerical simulations (radiation hydrodynamical models and nuclear reaction networks) run on high performance computers. The numerical calculations and companion theoretical analysis helped to better define the sites of astrophysical nucleosynthesis, and make predictions of the associated observables (nucleosynthetic yields and electromagnetic signals) that can be used to anticipate and interpret experimental data. This project enhances the scientific value of the nuclear data obtained from current and future experimental facilities (e.g., rare isotope beams) by improving our understanding of the astrophysical context and specific physical conditions under which we expect such reactions to take place in Nature.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Non-destructive simulation of node defects in additively manufactured lattice structures

Additive Manufacturing (AM), commonly referred to as 3D printing, offers the ability to not only fabricate geometrically complex lattice structures but parts in which lattice topologies in-fill volumes bounded by complex surface geometries. However, current AM processes produce defects on the strut and node elements which make up the lattice structure. This creates an inherent difference between the as-designed and as-fabricated geometries, which negatively affects predictions (via numerical simulation) of the lattice’s mechanical performance. Although experimental and numerical analysis of an AM lattice’s bulk structure, unit cell and struts have been performed, there exists almost no research data on the mechanical response of the individual as-manufactured lattice node elements. Here we propose a methodology that, for the first time, allows non-destructive quantification of the mechanical response of node elements within an as-manufactured lattice structure. A custom-developed tool is used to extract and classify each individual node geometry from micro-computed tomography scans of an AM fabricated lattice. Voxel-based finite element meshes are generated for numerical simulation and the mechanical response distribution is compared to that of the idealised computer-aided design model. The method demonstrates compatibility with Uncertainty Quantification methods that provide opportunities for efficient prediction of a population of nodal responses from sampled data. Overall, the non-destructive and automated nature of the node extraction and response evaluation is promising for its application in qualification and certification of additively manufactured lattice structures.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Modeling heat transport processes in enhanced geothermal systems: Validation study from EGS Collab Experiment 1

Heat recovery from enhanced geothermal systems (EGS) is a complex process involving heat transport in both fracture networks and rock formations. A comprehensive understanding of and the ability to model the underlying heat transport mechanisms is important for the success of EGS but remains challenging in practice due to the generally insufficient characterization of EGS reservoirs. In the present study, we analyze an extensively monitored intermediate-scale EGS field experiment performed in a well-characterized testbed. The high-resolution, high-quality measurements from the field experiment enable the development of a high-fidelity model incorporating a well-constrained fracture network. Based on the field experiment, we investigate the complex heat transport processes in an EGS-relevant environment and validate the capability of a numerical approach in simulating these inherently coupled heat transport processes. A series of numerical simulations were performed to study the effects of different heat transport mechanisms, including thermal convection with fracture flow, thermal conduction in rock formations, and the Joule-Thomson effect. The agreement of thermal responses between field measurements and simulation results indicates that our numerical approach can appropriately model the heat transport processes pertaining to heat recovery from EGS reservoirs.

Wu, Hui↗
DOE OSTI2024

Numerical eigen-spectrum slicing, accurate orthogonal eigen-basis, and mixed-precision eigenvalue refinement using OpenMP data-dependent tasks and accelerator offload

Performing a variety of numerical computations efficiently and, at the same time, in a portable fashion requires both an overarching design followed by a number of implementation strategies. All of these are exemplified below as we present transitioning the PLASMA numerical library from relying on dependence-driven large tasks to achieving utilization of fine grain tasking and offload to hardware accelerators while keeping its core dependence sets: OpenMP source code pragmas and runtime for most system-level functionality and basic low-level numerical kernels provided directly by hardware vendors or open source projects with vendor contributions. We also present new algorithmic methods and their efficient parallel implementations including fine grained tasking for eigen-spectrum slicing and offload for mixed-precision eigenvalue refinement. We provide performance, scaling, and numerical results showing sizable gains over the available solutions from either the open source and vendor-provided packages.

Luszczek, Piotr↗
DOE OSTI2024

A B-spline based gradient-enhanced micropolar implicit material point method for large localized inelastic deformations

The quasi-brittle response of cohesive-frictional materials in numerical simulations is commonly represented by softening plasticity or continuum damage models, either individually or in combination. However, classical models, particularly when coupled with non-associated plasticity, often suffer from ill-posedness and a lack of objectivity in numerical simulations. Moreover, the performance of the finite element method significantly degrades in simulations involving finite strains when mesh distortion reaches excessive levels. This represents a challenge for modeling cohesive-frictional materials, given their tendency to experience strongly localized deformations, such as those occurring during shear band dominated failure. Hence, accurate modeling of the response of cohesive-frictional solids is a demanding task. To address these challenges, we present an extension of the material point method (MPM) for the unified gradient-enhanced micropolar continuum, aiming at the analysis of finite localized inelastic deformations in cohesive-frictional materials. The generalized gradient-enhanced micropolar continuum formulation is employed to tackle challenges related to localization and softening material behavior, while the MPM addresses issues arising from excessive deformations. The method utilizes a B-spline formulation for the rigid background mesh to mitigate the well-known cell crossing errors of the MPM. To demonstrate the performance of the method, 2D and 3D numerical studies on localized failure in sandstone in plane strain compression and triaxial extension tests are presented. A comparison with finite element results confirms the suitability of the formulation. Moreover, an efficient numerical implementation of the formulation is presented, and it is demonstrated that the additional MPM specific overhead is negligible.

B-spline↗
DOE OSTI2025

A hybrid numerical and machine learning framework for evaluating the performance of a 780 cm 2 aqueous organic redox flow battery

Aqueous organic redox flow battery (AORFB) is a promising cost-competitive technology for large-scale energy storage. Among existing work, the dihydroxyphenazine (DHP)-based AORFB has demonstrated high energy density and low capacity degradation in 10 cm2 cells during lab tests. However, its commercial-scale performance in more complex environments remains unknown, posing a barrier for commercialization. To address this gap, this work presents a comprehensive performance evaluation of a 780 cm 2 DHP-based AORFB by combining physics-based numerical model, machine learning (ML)-based surrogate models, and ML-derived sensitivity quantification. Specifically, we first select 12 key battery parameters that include 10 physicochemical quantities and 2 operation quantities, then select 6 performance metrics that include energy efficiency (EE), discharging capacity, charging energy, and power losses due to concentration, activation, and ohmic over-potentials. With such selection, 12800 combinations of the 12 parameters are subsequently generated using the Latin Hypercube Sampling method. These combinations, together with 38 pre-defined State of Charge, are then integrated to a validated AORFB model developed in COMSOL to compute the performance metrics. With both input parameters and performance metrics, 60 deep neural network (DNN) surrogate models are then trained to approximate the relationship between the 10 physicochemical quantities and 6 performance metrics at each flow rate and current density. Sensitivity scores are then calculated based on the DNN models. Two additional sensitivity analysis tools, i.e., MARS, and SHAP, are also used to cross-validate the sensitivity scores from the DNN. The results demonstrate that 1) the standard potential ranks the first in controlling EE and charging energy, 2) the membrane conductivity is most critical for power loss and EE, and 3) specific area and reaction rate control activation power loss.

25 ENERGY STORAGE↗
DOE OSTI2023

Optimization of metal-supported solid oxide electrolysis cells with infiltrated catalysts

We report metal-supported solid oxide electrolysis cells (MS-SOECs) are being developed for steam-to-hydrogen electrolysis, especially for utilization of dynamic or intermittent electrical power from renewable sources. Various aspects of the electrocatalyst processing and composition, and metal support structure were explored. Catalyst materials, infiltration temperature and infiltration cycles were optimized for high performance and durability. Numerous catalyst materials were screened for both oxygen and steam electrodes. The oxygen catalyst had moderate impact on both initial cell performance and durability. Reducing Ni content in the steam electrode had little effect on durability, but reduced initial performance. Ex-situ XRD analysis and cell assessment of catalyst infiltration temperature revealed that the optimal range is 750–850 °C. The best cell performance and durability was achieved with LSCF-SDC oxygen electrocatalyst and SDC-Ni (60:40 vol%) steam electrocatalyst infiltrated 11 times at 800 °C and operated at 700 °C. At low steam content, a significant mass transport limitation on the steam side results in limiting current behavior. Thinner and more porous metal supports were implemented, and found to improve steam mass transport at low steam content, relevant for SOECs operating under high H 2 recycle rate or high steam utilization.

08 HYDROGEN↗
DOE OSTI2024

Evidence of scaling advantage for the quantum approximate optimization algorithm on a classically intractable problem

The quantum approximate optimization algorithm (QAOA) is a leading candidate algorithm for solving optimization problems on quantum computers. However, the potential of QAOA to tackle classically intractable problems remains unclear. Here, we perform an extensive numerical investigation of QAOA on the low autocorrelation binary sequences (LABS) problem, which is classically intractable even for moderately sized instances. We perform noiseless simulations with up to 40 qubits and observe that the runtime of QAOA with fixed parameters scales better than branch-and-bound solvers, which are the state-of-the-art exact solvers for LABS. The combination of QAOA with quantum minimum finding gives the best empirical scaling of any algorithm for the LABS problem. We demonstrate experimental progress in executing QAOA for the LABS problem using an algorithm-specific error detection scheme on Quantinuum trapped-ion processors. Our results provide evidence for the utility of QAOA as an algorithmic component that enables quantum speedups.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Overview on Shielding Analyses for the VENUS Instrument at SNS

VENUS, a world‐class versatile neutron imaging instrument, is under construction and is expected to be completed and ready to start commissioning in 2023. The range of cold to epithermal neutrons at SNS will give users of VENUS access to novel imaging methods, as well as to significantly improved existing methods. The instrument is being built on beam line 10 at Spallation Neutron Source (SNS) First Target Station (FTS) facing a decoupled poisoned hydrogen moderator. Instrument design, which includes optics, Front-End components and instrument enclosure components and shape started over 10 years ago and during this time had significant changes. All changes were supported by neutronics analyses to provide adequate shielding for both Front-End and instrument enclosure.Final optics design will provide the field of view (FOV) at the detector position (where the image is formed) to be as high as 0.20 by 0.20 m. Both the beam, which contains a large fraction of high-energy neutrons, and the desire for a large footprint on the detector are challenges for shielding design and budget, because the driving cost for the instrument is the beam line and enclosure shielding. For cost reason VENUS baseline design is optimized with respect to both the instrument cave footprint and its wall thickness, and Front-End shielding is tailored along the beam. Analyses are performed with the Monte Carlo particle transport code MCNPX version 2.7.0 to make choices on materials, thicknesses and configurations of the shielding. Numerous calculations are performed to meet space and coast constrain and to satisfy instrument physics needs and to comply with radiation protection requirements.

Gallmeier, Franz X.↗
DOE OSTI2024

Resolution requirements for high-speed droplet impact along a rigid wall

The mechanical and thermal loads that are produced on high-speed projectiles from a droplet are key in the projectile development.To accurately predict the pressure produced along the wall surface during impact, it is important to correctly model the air/water interface of the droplet. Novel modeling approaches such as the Phase-Field enable control of the (numerical) interface thickness. However, the interplay between the grid resolution (including potential adaptive mesh refinement) and the interface thickness is unclear. In this work, the resolution requirements for numerical simulations of a two-dimensional cylindrical water droplet impacting a rigid wall at Mach four are computationally investigated. Numerical simulations are performed using a second-order accurate scheme with adaptive mesh refinement and a Phase-Field interface model. The choice of numerical interface thickness and resolution parameters affects the maximum pressure and location of the maximum pressure along the wall. The overall maximum pressure does not converge even at high resolutions. The number of levels and thickness (number of cells across the interface) play a role in obtaining a well behaved solution.

Burrell, E.↗
DOE OSTI2022

A two-fluid single-column model of turbulent shallow convection. Part II: Single-column model formulation and numerics

The two-fluid single-column model of Thuburn et al. (Quart. J. R. Meteorol. Soc., 2019, 145, 1535–1550) is extended to include moisture and horizontal wind shear. Turbulent kinetic energy is introduced as a prognostic variable, dependence on a diagnosed boundary-layer height is removed, and subfilter fluxes are approximated using a two-fluid version of a Mellor–Yamada scheme. Three mechanisms for entrainment and detrainment processes are introduced, which represent entrainment of unstable air at the surface, forced detrainment of air at the top of the boundary/cloud layers, and turbulent mixing that relaxes the convective fluid to a reference profile. A semi-implicit Eulerian discretization replaces the semi-implicit semi-Lagrangian implementation of Thuburn et al. (Quart. J. R. Meteorol. Soc., 2019, 145, 1535–1550) to improve numerical stability and conservation. The equations for the implicit time step are solved using a quasi-Newton method, which is shown to perform well in numerical tests for conservation and convergence. The two-fluid single-column model presented in this article will be applied to simulations of shallow cumulus convection in Part III.

54 ENVIRONMENTAL SCIENCES↗