Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Operator approximation”

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 1,225 records · Page 68

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↗

Dual-module humidity pump for efficient air dehumidification: Demonstration and performance limitations

Condensation dehumidification in conventional air conditioning technologies is energy-intensive, accounting for up to 50% of building cooling energy used in some climates. Selective vacuum membrane dehumidification (VMD) is one of the leading alternative dehumidification technologies due to its potential for significant energy savings, and the “dual-module humidity pump” is one of the most promising VMD concepts. Here, this work is the first to provide experimental proof-of-concept for the dual-module humidity pump system and provides the first thermodynamic modeling framework that accounts for realistic steady-state operating limitations, both of which are lacking in the current literature. Additionally, this work is the first to provide a system design solution that overcomes practical challenges associated with air accumulation in the vacuum channels. The experimental results in this work show that the current prototype can remove up to 45% of the water vapor in the air stream, and the vapor pressure difference in the vapor rejection module needs to be approximately 2–4 times greater than that of the dehumidification module in order to maintain balanced mass transfer. The thermodynamic model applied to typical air conditioning conditions shows that the ideal dehumidification (latent) COPs can reach up to 40, but practical COPs are limited to approximately 10. Furthermore, the model shows that the overall energy efficiency increases as the membrane air selectivity increases, though this improvement gradually starts to diminish when the membrane selectivity is increased above 10,000.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Reduced-order modeling on a near-term quantum computer

Quantum computing is an advancing area of research in which computer hardware and algorithms are developed to take advantage of quantum mechanical phenomena. In recent studies, quantum algorithms have shown promise in solving linear systems of equations as well as systems of linear ordinary differential equations (ODEs) and partial differential equations (PDEs). Reducedorder modeling (ROM) algorithms for studying fluid dynamics have shown success in identifying linear operators that can describe flowfields, where dynamic mode decomposition (DMD) is a particularly useful method in which a linear operator is identified from data. In this work, DMD is reformulated as an optimization problem to propagate the state of the linearized dynamical system on a quantum computer. This reformulation was chosen as a means of facilitating implementation on a near-term quantum computer. Quadratic unconstrained binary optimization (QUBO), a technique for optimizing quadratic polynomials in binary variables, allows for quantum annealing algorithms to be applied. A quantum circuit model (quantum approximation optimization algorithm, QAOA) is utilized to obtain predictions of the state trajectories. Results are shown for the quantum-ROM predictions for flow over a 2D cylinder at Re = 220 and flow over a NACA0009 airfoil at Re = 500 and α = 15°. The quantum-ROM predictions are found to depend on the number of bits utilized for a fixed point representation and the truncation level of the DMD model. Comparisons with DMD predictions from a classical computer algorithm are made, as well as an analysis of the computational complexity and prospects for future, more fault-tolerant quantum computers.

97 MATHEMATICS AND COMPUTING↗

Photoinduced Charge Carrier Dynamics and Electron Injection Efficiencies in Au Nanoparticle-Sensitized TiO 2 Determined with Picosecond Time-Resolved X-ray Photoelectron Spectroscopy

Progress in the development of plasmon-enabled light-harvesting technologies requires a better understanding of their fundamental operating principles and current limitations. In this work, we employ picosecond time-resolved X-ray photoemission spectroscopy to investigate photoinduced electron transfer in a plasmonic model system composed of 20 nm sized gold nanoparticles (NPs) attached to a nanoporous film of TiO 2 . The measurement provides direct, quantitative access to transient local charge distributions from the perspectives of the electron donor (AuNP) and the electron acceptor (TiO 2 ). On average, approximately two electrons are injected per NP, corresponding to an electron injection yield per absorbed photon of 0.1%. Back electron transfer from the perspective of the electron donor is dominated by a fast recombination channel proceeding on a time scale of 60 ± 10 ps and a minor contribution that is completed after ~1 ns. The findings provide a detailed picture of photoinduced charge carrier generation in this NP–semiconductor junction, with important implications for understanding achievable overall photon-to-charge conversion efficiencies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Kohn-Sham Solver (KSSOLV) v2.0

KSSOLV is a MATLAB toolbox for solving Kohn-Sham density functional theory based electronic structure eigenvalue problems. It uses an object oriented features of MATLAB to represent atom, molecules, wavefunctions and Hamiltonians and their operations. It is designed to make it easier for users to prototype and test new algorithms for solving the Kohn-Sham problem. KSSOLV2.0 contains significant improvement over the original KSSOLV described in a paper published in ACM Transaction on Mathematical Software (attached). In addition to performing ground state calculation for small molecules, it can also perform geometry optimization for both molecules and solids. It uses standard pseudopotentials and implements local density approximation, generalized gradient approximation and hybrid functionals. Future releases will also include time-dependent DFT and post DFT calculations such as the GW quasi-particle energy calculation and Bethe-Salpeter equation solver for optical absorption.

Yang, Chao↗

Phenomenological R-Matrix parameterization of direct, doorway, and compound nuclear reactions [Slides]

The presented R-matrix formalism could improve evaluations of resolved resonance nuclear cross section data by parameterizing direct reaction processes via an orthonormal channel mixing matrix and doorway state processes. Together, they can parameterize direct-semidirect capture. A formal similarity among the R-, K-, T- matrix formalisms was leveraged to parameterize direct and doorway processes. The advantage of this R-matrix parameterization is that it is easily extensible due to its application of projection operators directly upon the (H-E) on p. 6.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Identifying Controlling Variables for Mercury Vapors in Alpha-4 at Y-12: Two Year Data Collection Update

Multiple sensor packages were deployed at Alpha-4 by SRNL, in collaboration with United Cleanup Oak Ridge LLC (UCOR), to monitor mercury vapor concentrations and meteorological parameters. These sensors collected data, inside and outside of the legacy-use facility, for approximately two years. Though data gaps still exist, particularly in colder months, several controlling variables were identified that govern mercury vapor concentrations within Alpha-4. Temperature, barometric pressure gradients, humidity, and wind speed have been identified as controlling variables. A strong positive correlation was seen between mercury vapor concentrations and temperature which generally followed diurnal fluctuations. Temperatures below approximately 10 degrees Celsius did not show any spikes above the PEL, indicating more work can be performed at any time during the winter months – more data should be collected to confirm consistency in this finding. Additionally, spikes in mercury vapor concentrations above that of the permissible exposure limit (PEL; 100 µg/m 3 ) occurred primarily in late afternoon or evening/overnight hours (between 3 PM and 6 AM), which suggests D&D operations might be best scheduled during morning or daytime hours prior to the late afternoon. However, a limited number of spikes did occur outside of the identified window, although this may be attributed to disturbances in air flow and mercury vapor release from work activities performed inside of the Alpha-4 building. The analysis conducted allows for a strong predictive capability for estimating mercury vapor concentrations based upon accurate meteorological parameters. Still, additional data collection, particularly in the winter months, could help to strengthen the predictive power and validate the identified data trends. Further, increased temporal resolution could also help to better characterize the incipient stages of the increases and decreases in the mercury vapor concentration. Continued monitoring support by SRNL at Y-12 is underway at Alpha-4 to further close remaining data gaps and support deactivation and decommissioning work. Within a collaborative effort with UCOR, the SRNL team is collecting mercury vapor data to study the efficacy of a novel mercury suppressant, FerroBlack® which was recently deployed at Alpha-4. In addition, modifications to the current monitoring setup to increase measurement resolution is also being investigated.

54 ENVIRONMENTAL SCIENCES↗

A Realistic Full-Scale 3D Modeling of Turning Using Coupled Smoothed Particle Hydrodynamics and Finite Element Method for Predicting Cutting Forces

Computational modelling is an effective technique for understanding the complex physics of machining. Large deformations, material separation, and high computational requirements are the key challenges faced while simulating machining. This work introduces a full-scale three-dimensional model of turning operations using a combined approach based on the Smoothed Particle Hydrodynamics (SPH) and Finite Element (FE) methods. By exploiting the advantages of each method, this approach leads to high-fidelity coupled SPH-FE machining models. Cutting forces and chip morphology are the primary results of interest. The machining models are validated with the results of turning experiments. Two-dimensional machining model underpredicts the cutting force and feed force by approximately 49% and 70%, respectively. Moreover, passive force cannot be predicted using the two-dimensional model. On the other hand, with the three-dimensional models developed in this manuscript, the difference between the total simulated force and experimentally measured force is ∼17%. The chip morphologies correlate with experiments in terms of the direction of the chip movement and the “long” continuous chips observed while turning Al 6061. This work expands the realm of machining simulations from two-dimensional orthogonal machining or sectional three-dimensional model to a full-scale realistic simulation. The encouraging simulation results show the potential to study more complex phenomena, such as machining stability and tool path modulation.

42 ENGINEERING↗

A digital twin solution for floating offshore wind turbines validated using a full-scale prototype

Abstract. In this work, we implement, verify, and validate a physics-based digital twin solution applied to a floating offshore wind turbine. The digital twin is validated using measurement data from the full-scale TetraSpar prototype. We focus on the estimation of the aerodynamic loads, wind speed, and section loads along the tower, with the aim of estimating the fatigue lifetime of the tower. Our digital twin solution integrates (1) a Kalman filter to estimate the structural states based on a linear model of the structure and measurements from the turbine, (2) an aerodynamic estimator, and (3) a physics-based virtual sensing procedure to obtain the loads along the tower. The digital twin relies on a set of measurements that are expected to be available on any existing wind turbine (power, pitch, rotor speed, and tower acceleration) and motion sensors that are likely to be standard measurements for a floating platform (inclinometers and GPS sensors). We explore two different pathways to obtain physics-based models: a suite of dedicated Python tools implemented as part of this work and the OpenFAST linearization feature. In our final version of the digital twin, we use components from both approaches. We perform different numerical experiments to verify the individual models of the digital twin. In this simulation realm, we obtain estimated damage equivalent loads of the tower fore–aft bending moment with an accuracy of approximately 5 % to 10 %. When comparing the digital twin estimations with the measurements from the TetraSpar prototype, the errors increased to 10 %–15 % on average. Overall, the accuracy of the results is promising and demonstrates the possibility of using digital twin solutions to estimate fatigue loads on floating offshore wind turbines. A natural continuation of this work would be to implement the monitoring and diagnostics aspect of the digital twin to inform operation and maintenance decisions. The digital twin solution is provided with examples as part of an open-source repository.

17 WIND ENERGY↗

Carnot Analysis of Heat Pump Drying: Ideal Efficiency and Dry Time

Drying processes are important in appliances and in industry, and clothes drying accounts for approximately 3% of residential primary energy consumption in the US. Globally, heat pump tumble dryers (HPD) are increasing their market share against the ubiquitous electric resistance tumble dryer (ERD). In this work, efficiency and dry time limits are defined for ideal HPDs, for both closed air cycle (unvented) and open air cycle (vented). These limits are compared with the limits for ERD. The traditional Carnot efficiency limit for an ideal heat pump, operating between a hot and a cold thermal reservoir, does not apply directly to clothes dryers. One reason dryers require a novel analysis is the presence of additional degrees of freedom, since the hot and cold temperatures are floating, unfixed by ambient conditions. Furthermore, dryers can operate in a closed or open air cycle, and each requires a different analysis. In the closed (unvented) case, the hot and cold temperatures are coupled to each other; while in the open (vented) case, the hot and cold temperatures are both independent free variables. This paper provides an analysis of the fundamental efficiency limits of ERDs and HPDs, which can inform the design and performance limits of evaporative drying technology.

Gluesenkamp, Kyle↗

Analysis of SPND Material Candidates for Real-Time Neutron Spectrum Monitoring Poster

The future of the nuclear industry will be shaped by advanced reactor designs requiring innovative instrumentation for safe operation. A key area of innovation is the ability to accurately and reliably monitor the neutron energy spectrum during operations, which would also enhance the experimental capabilities of existing research reactors. To address this need, our research focuses on the potential use of Self-Powered Neutron Detectors (SPNDs) for real-time neutron spectrum monitoring. We investigated various candidate materials for SPNDs by evaluating their energy-dependent reaction rates. Promising materials were further modeled using MCNP to determine the free electron production per neutron across the neutron energy spectrum, approximating their energy-dependent signal sensitivity. Our findings identified Samarium (Sm) and Gadolinium (Gd) as top candidates for thermal neutrons, Hafnium (Hf) and Thulium (Tm) for epithermal neutrons, and Tantalum (Ta) for fast neutrons. These materials will undergo future experimental validation to confirm the modeled responses and assess their feasibility for real-time neutron spectrum measurements.

46 - INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AN↗

Resonance compensation at the CERN PS booster aided by Bayesian optimization and BOBYQA

The CERN Proton Synchrotron Booster (PSB) operation involves the crossing of multiple resonance lines in the tune diagram. Loss maps from dynamic tune scans are a helpful way to visualize and quantify the strength of such resonances. Sextupole and octupole correctors can be used in order to partially or fully compensate multiple resonance lines, i.e., third and fourth order lines. The following work explores the application of advanced optimization algorithms such as Bayesian Optimization and Bound Optimization By Quadratic Approximation (BOBYQA) in order to compensate these resonance lines with available correctors.

43 PARTICLE ACCELERATORS↗

Multifidelity Active Learning for Failure Estimation of TRISO Nuclear Fuel

The Tristructural isotropic (TRISO)-coated particle fuel is a robust nuclear fuel proposed to be used for multiple modern nuclear technologies. Therefore, characterizing its safety is vital for the reliable operation of nuclear technologies. However, the TRISO fuel failure probabilities are small and the computational model is time consuming to evaluate them using traditional Monte Carlo-type approaches. In the paper, we present a multifidelity active learning approach to efficiently estimate small failure probabilities given an expensive computational model. Active learning suggests the next best training set for optimal subsequent predictive performance and multifidelity modeling uses cheaper low-fidelity models to approximate the high-fidelity model output. After presenting the multifidelity active learning approach, we apply it to efficiently predict TRISO failure probability and make comparisons to the reference results.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Self-interaction-corrected Kohn–Sham effective potentials using the density-consistent effective potential method

Density functional theory (DFT) and beyond-DFT methods are often used in combination with photoelectron spectroscopy to obtain physical insights into the electronic structure of molecules and solids. The Kohn–Sham eigenvalues are not electron removal energies except for the highest occupied orbital. The eigenvalues of the highest occupied molecular orbitals often underestimate the electron removal or ionization energies due to the self-interaction (SI) errors in approximate density functionals. Here, we adapt and implement the density-consistent effective potential method of Kohut, Ryabinkin, and Staroverov (2014) to obtain SI-corrected local effective potentials from the SI-corrected Fermi–Löwdin orbitals and density in the Fermi–Löwdin orbital self-interaction correction scheme. The implementation is used to obtain the density of states (photoelectron spectra) and HOMO–LUMO gaps for a set of molecules and polyacenes. Good agreement with experimental values is obtained compared to a range of SI uncorrected density functional approximations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A Model for Evaluating the Configuration and Dispatch of PV Plus Battery Power Plants

An open-source model was developed to optimize energy storage operation for photovoltaic- (PV-) plus-battery systems with AC-coupled and DC-coupled configurations. It includes the ability to use forecast energy prices to optimize battery charge and discharge on a rolling time horizon. The model allows for exploration of different configurations, including capital costs, inverter performance, dispatch flexibility, and capturing otherwise clipped energy for the DC-coupled system. The model can run 20 full years of hourly data in approximately two seconds, allowing comparison of a large number of configurations. We applied the model in a test case demonstrating reduced inverter clipping for DC-coupled systems and yielded slightly higher overall value than AC-coupled systems, with an approximately 2 percent increase in internal rate of return or benefit/cost ratio. Our results show that at current estimated prices for lithium-ion battery systems, large-scale PV-plus-battery plants are economically viable under the right conditions, with the configuration playing a role in system flexibility and performance. This model provides the ability for project developers, industry professionals, and researchers to use readily available software to quickly evaluate and design these systems.

41 EE - Solar Energy Technologies Office (EE-4S)↗

Seismic tomography of compressional wave velocity and attenuation structure for Makushin Volcano, Alaska

Located in the central portion of the Alaska-Aleutian subduction zone, Makushin Volcano is among the most active volcanoes in the United States and has been classified as a very high threat by the National Volcano Early Warning System based on its eruptive history and proximity to the City of Unalaska and international air routes. In 2015, five standalone seismic stations and three mini seismic arrays of 15 stations each were deployed on Unalaska Island to supplement the Alaska Volcano Observatory permanent seismic network. This temporary array was operational for one year. Taking advantage of the increased azimuthal coverage and the array's increased earthquake detection capability, we developed body-wave Vp and Vp/Vs seismic images of the velocity structure beneath the volcano. Our results show a complex structure with the upper 4 km of the crust. The shallow high-Vp features possibly delineate remnant magma pathways or conduits. Low-Vp regions are found east of the caldera at approximately 4–8 km depth. This is in agreement with previous seismic tomographic work and InSAR-based geodetic models, which had identified this region as a possible long-term source of magma. A three-dimensional compressional wave attenuation model complemented the velocity models. In general, we observe regions of corresponding low P-wave velocity and high attenuation east of the summit caldera between 5 and 7 km depth, supporting the presence of a melt-rich zone. The absolute Vp values (5.5–6.0 km/s) in this region are consistent with the presence of a magma body with a mafic signature. Here, we suggest that a triggering mechanism that involves mafic injections into a more evolved shallow magma chamber is representative of the current magmatic system below the volcano, similarly to what petrological studies have inferred for past eruptions.

58 GEOSCIENCES↗

Toroidal Field Reversing Switch

A motor-driven, high-current switch is being designed by Stäubli Electrical Connectors in collaboration with Princeton Plasma Physics Laboratory to reverse the DIII-D Toroidal field coil current. The existing field reversing procedure requires approximately 4 h of labor performed between run days to reconfigure the bus work, thus requiring a separate run day to complete experimental scans with both field directions. This typically results in different wall conditions which complicate the results or make them invalid. The proposed switch will complete the reversal process in approximately 4-5 min, between shots, so that wall conditions are not a factor, and easily fits within the 10-12 min DIII-D shot cycle. The switch is a compact design fitting into a 157.5 × 41 cm envelope. The switch uses two copper pins that slide between Stäubli MULTILAM contacts. The switch is rated for total current of 180-kA dc for 10 s with a 12-min repetition rate; the maximum operating coil current is 124-kA dc. The switch diagnostics include: voltage drop, current monitor, thermal monitors, position sensors, and motor torque. Here, the switch is rated for 300 cycles per year, and 3000 cycles overall before contact replacement, with contact inspection yearly. The design details and operation of the switch are covered in this article.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Stochastic tensor contraction for quantum chemistry

Many computational methods in ab initio quantum chemistry are formulated in terms of high-order tensor contractions, whose cost determines the size of system that can be studied. We introduce stochastic tensor contraction to perform such operations with greatly reduced cost, and present its application to the gold-standard quantum chemistry method, coupled cluster theory with up to perturbative triples. For total energy errors more stringent than chemical accuracy, we reduce the computational scaling to that of mean-field theory, while starting to approach the mean-field absolute cost, thereby challenging the existing cost-to-accuracy landscape. Benchmarks against state-of-the-art local correlation approximations further show that we achieve an order-of-magnitude improvement in both total computation time and error, with significantly reduced sensitivity to system dimensionality and electron delocalization. We conclude that stochastic tensor contraction is a powerful computational primitive to accelerate a wide range of quantum chemistry.

Chemical Physics (physics.chem-ph)↗