Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “simulation analytics”

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,063 records · Page 59

Dark Energy Survey Year 3 results: Cosmology with moments of weak lensing mass maps

We present a cosmological analysis using the second and third moments of the weak lensing mass (convergence) maps from the first three years of data (Y3) data of the Dark Energy Survey. The survey spans an effective area of 4139 square degrees and uses the images of over 100 million galaxies to reconstruct the convergence field. The second moment of the convergence as a function of smoothing scale contains information similar to standard shear 2-point statistics. The third moment, or the skewness, contains additional non-Gaussian information. The data is analyzed in the context of the Λ CDM model, varying five cosmological parameters and 19 nuisance parameters modeling astrophysical and measurement systematics. Our modeling of the observables is completely analytical, and has been tested with simulations in our previous methodology study. We obtain a 1.7% measurement of the amplitude of fluctuations parameter S 8 ≡ σ 8 ( Ω m / 0.3 ) 0.5 = 0.784 ± 0.013 . The measurements are shown to be internally consistent across redshift bins, angular scales, and between second and third moments. In particular, the measured third moment is consistent with the expectation of gravitational clustering under the Λ CDM model. The addition of the third moment improves the constraints on S 8 and Ω m by ∼ 15 % and ∼ 25 % compared to an analysis that only uses second moments. We compare our results with Planck constraints from the cosmic microwave background, finding a 2.2 – 2.8 σ tension in the full parameter space, depending on the combination of moments considered. The third moment, independently, is in 2.8 σ tension with Planck, and thus provides a cross-check on the analyses of 2-point correlations.

79 ASTRONOMY AND ASTROPHYSICS↗

Kinetic relaxation and nucleation of Bose stars in self-interacting wave dark matter

We revisit kinetic relaxation and soliton/boson star nucleation in fuzzy scalar dark matter featuring short-ranged self-interactions H int = − λ | ψ | 4 / 2 m 2 , alongside gravitational self-interactions. We map out the full curve of nucleation timescale for both repulsive ( λ < 0 ) and attractive ( λ > 0 ) short-ranged self-interaction strength and in doing so reveal two new points. Firstly, besides the two usual terms, ∝ G 2 and ∝ λ 2 , in the total relaxation rate Γ relax , there is an additional cross term ∝ G λ arising due to interference between gravitational and short-ranged self-interaction scattering amplitudes. This yields a critical repulsive interaction strength λ cr ≃ − 2 π G m 2 / v 0 2 , at which the relaxation rate is smallest and serves as the transition point between typical net attractive self-interaction ( λ ≳ λ cr ) and net repulsive self-interaction ( − λ ≳ − λ cr ). Secondly, while in the net attractive regime, nucleation timescale is similar to inverse relaxation timescale τ nuc ∼ Γ relax − 1 ; in the net repulsive regime, nucleation occurs at a delayed time τ nuc ∼ ( λ / λ cr ) Γ relax − 1 . We confirm our analytical understanding by performing 3D field simulations with varying average mass density ρ ¯ , box size L and grid size N . Published by the American Physical Society 2024

Astronomy & Astrophysics↗

Dynamic Power Network State Estimation with Asynchronous Measurements

The operation of distribution networks is becoming increasingly volatile, due to fast variations of renewables and, hence, net-loading conditions. To perform a reliable state estimation under these conditions, this paper considers the case where measurements from meters, phasor measurement units, and distributed energy resources are collected and processed in real time to produce estimates of the state at a fast time scale. Streams of measurements collected in real time and at heterogenous rates render the underlying processing asynchronous, and poses severe strains on workhorse state estimation algorithms. In this work, a real-time state estimation algorithm is proposed, where data are processed on the fly. Starting from a regularized least-squares model, and leveraging appropriate linear models, the proposed scheme boils down to a linear dynamical system where the state is updated based on the previous estimate and on the measurement gathered from a few available sensors. The estimation error is shown to be always bounded under mild condition. Numerical simulations are provided to corroborate the analytical findings.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Improving Progressive Retrieval for HPC Scientific Data using Deep Neural Network

As the disparity between compute and I/O on high-performance computing systems has continued to widen, it has become increasingly difficult to perform post-hoc data analytics on full-resolution scientific simulation data due to the high I/O cost. Error-bounded data decomposition and progressive data retrieval framework has recently been developed to address such a challenge by performing data decomposition before storage and reading only part of the decomposed data when necessary. However, the performance of the progressive retrieval framework has been suffering from the over-pessimistic error control theory, such that the achieved maximum error of recomposed data is significantly lower than the required error. Therefore, more data than required is fetched for recomposition, incurring additional I/O overhead. In order to tackle this issue, we propose a DNN-based progressive retrieval framework that can better identify the minimum amount of data to be retrieved. Our contributions are as follows: 1) We provide an in-depth investigation of the recently developed progressive retrieval framework; 2) We propose two designs of prediction models (named D-MGARD and E-MGARD) to estimate the amount of retrieved data size based on error bounds. 3) We evaluate our proposed solutions using scientific datasets generated by real-world simulations from two domains. Evaluation results demonstrate the effectiveness of our solution in accurately predicting the amount of retrieval data size, as well as the advantages of our solution over the traditional approach to reducing the I/O overhead. Based on our evaluation, our solution is shown to read significantly less data (5% - 40% with D-MGARD, 20% - 80% with E-MGARD).

Wang, Jinzhen↗

Robust Stabilization of Inverter-Based Resources Using Virtual Resistance-Based Control

This letter proposes a virtual resistance-based nonlinear control to stabilize and robustify the current layer of inverter-based resources, subject to the grid voltage disturbances, and the grid parameter uncertainties. A class of virtual resistances is proposed and analyzed using concepts from dissipative systems theory. Moreover, specific nonlinear virtual resistance-based controllers are derived, with their corresponding performance analytically bounded. Here, the theoretical and simulation results show that the proposed nonlinear virtual resistance-based controllers significantly reduce the L 2 gain of the closed-loop error system to grid voltage variation and parametric uncertainties. Significant improvements of the transient and steady-state current responses are also demonstrated over linear virtual resistance counterparts.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Online State Estimation for Time-Varying Systems

The paper investigates the problem of estimating the state of a time-varying system with a linear measurement model; in particular, the paper considers the case where the number of measurements available can be smaller than the number of states. In lieu of a batch linear least-squares (LS) approach well-suited for static networks, where a sufficient number of measurements could be collected to obtain a full-rank design matrix the paper proposes an online algorithm to estimate the possibly time-varying state by processing measurements as and when available. The design of the algorithm hinges on a generalized LS cost augmented with a proximal-point-type regularization. With the solution of the regularized LS problem available in closed-form, the online algorithm is written as a linear dynamical system where the state is updated based on the previous estimate and based on the new available measurements. Conditions under which the algorithmic steps are in fact a contractive mapping are shown, and bounds on the estimation error are derived for different noise models. Numerical simulations are provided to corroborate the analytical findings.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Improving Longitudinal Impedance of Inductive Inserts

In this article, we discuss methods for improving the longitudinal impedance of inductive inserts used for space charge compensation in the proton storage ring (PSR) at Los Alamos Neutron Science Center (LANSCE) Facility. The PSR relies on inductive inserts to counteract the capacitive impedance produced by the space charge of the proton beam, which can otherwise lead to beam debunching. Existing inductive inserts, which utilize Toshiba M4C21A ferrite, have been effective but require heating to reduce resistive components that can cause instabilities at higher frequencies. This article explores the potential of alternative ferrite materials, specifically National Magnetics C2050, to enhance the performance of these inserts by reducing the need for heating and minimizing resistive impedance. Through both analytical modeling and CST Particle Studio simulations, it is demonstrated that optimized ferrite materials and geometries can significantly improve space charge compensation while mitigating instability risks, ultimately enhancing the performance of the PSR.

43 PARTICLE ACCELERATORS↗

Enabling Scalable and Extensible Memory-mapped Datastores in Userspace

Exascale workloads are expected to incorporate data-intensive processing in close coordination with traditional physics simulations. These emerging scientific, data-analytics and machine learning applications need to access a wide variety of datastores in flat files and structured databases. Programmer productivity is greatly enhanced by mapping datastores into the application process's virtual memory space to provide a unified “in-memory” interface. Currently, memory mapping is provided by system software primarily designed for generality and reliability. However, scalability at high concurrency is a formidable challenge on exascale systems. Also, there is a need for extensibility to support new datastores potentially requiring HPC data transfer services. In this article, we present UMap , a scalable and extensible userspace service for memory-mapping datastores. Furthermore, through decoupled queue management, concurrency aware adaptation, and dynamic load balancing, UMap enables application performance to scale even at high concurrency. We evaluate UMap in data-intensive applications, including sorting, graph traversal, database operations, and metagenomic analytics. Our results show that UMap as a userspace service outperforms an optimized kernel-based service across a wide range of intra-node concurrency by 1.22-1.9 × . We performed two case studies to demonstrate UMap 's extensibility. First, a new datastore residing in remote memory is incorporated into UMap as an application-specific plugin. Second, we present a persistent memory allocator Metall built atop UMap for unified storage/memory.

97 MATHEMATICS AND COMPUTING↗

Enabling Grid-Forming Control with Fault Ride-Through in Unbalanced Distribution Networks

Distribution networks are often unbalanced, causing oscillatory responses in inverter control designed for balanced conditions. Here, to address this problem, this paper proposes a novel time-domain transformation appropriate for inverter control and enables the decomposition of three-phase unbalanced signals into constant positive and negative components. Relations useful for calculating unbalanced active and reactive power are derived from first principle, providing insight into vector products of unbalanced three-phase signals. Furthermore, a grid-forming control effective under unbalanced conditions is developed, which delivers superior performance while meeting UNIFI1 specifications for grid-forming control under unbalanced conditions. specifications applicable to category 4 inverter-based resource, like setting and regulating frequency/voltage, providing voltage support, sharing active power, injecting negative sequence current, and riding through faults. A current limiter is proposed for safe fault ride-through and integrates with the grid-forming control featuring frequency/voltage droop controllers and current and voltage control loops. The transformation of interconnected inverters is formulated and stability of the proposed control analyzed to support robust parameter selections. The effectiveness of the proposed transformation and grid-forming control is demonstrated through analytical results and real-time simulation of a IEEE 123 distribution network on the Real-Time Digital Simulator. Comparison with existing methods shows that the proposed strategy satisfies the UNIFI specifications with a much better performance.

24 - POWER TRANSMISSION AND DISTRIBUTION↗

Connecting Material Properties and Redox Flow Cell Cycling Performance through Zero-Dimensional Models

Improvements in redox flow battery (RFB) performance and durability can be achieved through the development of new active materials, electrolytes, and membranes. While a rich design space exists for emerging materials, complex tradeoffs challenge the articulation of unambiguous target criteria, as the relationships between component selection and cycling performance are multifaceted. Here, we derive zero-dimensional, analytical expressions for mass balances and cell voltages under galvanostatic cycling, enabling direct connections between material/electrolyte properties, cell operating conditions, and resulting performance metrics (e.g., energy efficiency, capacity fade). To demonstrate the utility of this modeling framework, we highlight several considerations for RFB design, including upper bound estimation, active species decay, and membrane/separator conductivity-selectivity tradeoffs. Furthermore, we also discuss modalities for extending this framework to incorporate kinetic losses, distributed ohmic losses, and multiple spatial domains. Importantly, because the mass balances are solved analytically, hundreds of cycles can be simulated in seconds, potentially facilitating detailed parametric sweeps, system optimization, and parameter estimation from cycling experiments. More broadly, this approach provides a means for assessing the impact of cell components that simultaneously influence multiple performance-defining processes, aiding in the elucidation of key descriptors and the identification of favorable materials combinations for specific applications.

25 ENERGY STORAGE↗

Spectral signatures of transverse optical modes in semiconductor nanowires

In order to build nanophotonic devices, it is important to understand and ultimately control the optical mode structure within potential components such as nanoscale waveguides. However, experimental characterization of such modes in the optical regime is difficult due to the nanoscale dimensions of such components and the perturbations that would be induced by a near-field probe. Here, we demonstrate a probe-free, far-field method to characterize the optical modes within GaN nanowires (NWs) based on a novel off-axis scanning confocal microscope system. Using this microscope, we observe spectral signatures resulting from lateral leakage of waveguide modes when they exceed their respective cutoff limits. We identify these modes within hyperspectral images using an analytical model coupled with finite element simulations. The model can also be used to predict the spectral signatures for given geometrical and dielectric parameters, which enabled us to deduce the transverse dimension of the NW from hyperspectral images with an accuracy of ∼ <#comment/> 30 nm .

Wang, Maoji↗

An investigation into the potential for wind turbines to cause barotrauma in bats

The high rates of bat mortality caused by operating wind turbines is a concern for wind energy and wildlife stakeholders. One theory that explains the mortality is that bats are not only killed by impact trauma, but also by barotrauma that results from exposure to the pressure variations caused by rotating turbine blades. To date, no published research has calculated the pressure changes that bats may be exposed to when flying near wind turbines and then used these data to estimate the likelihood that turbines cause barotrauma in bats. To address this shortcoming, we performed computational fluid dynamics simulations of a wind turbine and analytical calculations of blade-tip vortices to estimate the characteristics of the sudden pressure changes bats may experience when flying near a utility-scale wind turbine. Because there are no data available that characterize the pressure changes that cause barotrauma in bats, we compared our results to changes in pressure levels that cause barotrauma and mortality in other mammals of similar size. This comparison shows that the magnitude of the low-pressures bats experience when flying near wind turbines is approximately 8 times smaller than the pressure that causes mortality in rats, the smallest mammal for which data are available. The magnitude of the high-pressures that bats may experience are approximately 80 times smaller than the exposure level that causes 50% mortality in mice, which have a body mass similar to several bat species that are killed by wind turbines. Further, our results show that for a bat to experience the largest possible magnitude of low- and high-pressures, they must take very specific and improbable flight paths that skim the surface of the blades. Even a small change in the flight path results in the bat being hit by the blade or experiencing a much smaller pressure change. Accordingly, if bats have a physiological response to rapid low- and high-pressure exposure that is similar to other mammals, we conclude that it is unlikely that barotrauma is responsible for a significant number of turbine-related bat fatalities, and that impact trauma is the likely cause of the majority of wind-turbine-related bat fatalities.

17 WIND ENERGY↗

Numerical Modelling of the Optical Stochastic Cooling Experiment at IOTA

A proof-of-principle optical-stochastic cooling (OSC) experiment is currently in its commissioning phase at the Fermilab’s IOTA ring. In support of this experiment, we recently implemented an OSC element in the ELEGANT tracking program. The model, based on a semi-analytic description of OSC [*], supports the simulation of a large number of macroparticles (10⁴-10⁶) over many turns (10⁶). This paper showcases the simulation capabilities to investigate the beam dynamics in the presence of cooling (or self-interacting radiation field in general) and quantify the impact of various sources of error (e.g. transverse and phase jitter), guide data analysis.

43 PARTICLE ACCELERATORS↗

Microstructure-Based Modeling of Inner Oxygen Pressure in Solid Oxide Electrolysis Cells

One major degradation mechanism during long-term operation of solid oxide electrolysis cells (SOECs) is delamination of oxygen electrodes (OEs). The driving force for the electrode delamination could be the generated high inner oxygen pressure near the electrode-electrolyte interface during operation. However, the effects of transport properties and electrode thickness on the inner oxygen partial pressure are not well understood. Here, a microstructure-based electrochemical model which includes the conduction of electrons and oxygen ions coupled with Butler-Volmer-type chemical reactions at triple-phase-boundaries (TPBs), is employed to investigate the oxygen pressure in lanthanum strontium manganate (LSM)-based SOECs. The model is applied to both two-dimensional (2D) prototype microstructures and three-dimensional (3D) realistic microstructures, and the oxygen pressure is analyzed as a function of transport properties and electrode thickness under both potentiostatic and galvanostatic operations. The simulation results suggest strategies to suppress electrode delamination. The simulation results are compared to an analytical solution, and the discrepancies are attributed to the Butler-Volmer-type kinetics included in the microstructure-based model.

Xue, Fei↗

Modeling The Nucleosynthetic Imprint of Stellar Merger Phenomena - Final Technical Report

This DOE CAREER project developed a comprehensive, multi-physics framework for modeling stellar mergers and their observable consequences. The work successfully integrated analytical models, three-dimensional hydrodynamic simulations, stellar evolution calculations, nucleosynthesis, and radiation transport into a unified pipeline. Key scientific advances include demonstrating stellar mergers as a unifying explanation for systems such as Betelgeuse and R Coronae Borealis stars, quantifying merger-driven chemical signatures, and establishing the connection between merger physics and circumstellar environments. A major outcome of the project is the development of SuperLite, an open-source Monte Carlo radiation transport code that enables the generation of synthetic spectra for astrophysical transients. This work bridges dynamical, thermal, and radiative timescales and provides a powerful, DOE-relevant computational framework for interpreting observations of supernovae and related transient phenomena, while also contributing to workforce development through the training of postdoctoral researchers and graduate students.

Chatzopoulos, Emmanouil [Louisiana State Universit↗

Online Control of Process Variance Using Feedback

For many stochastic control problem formulations, design objectives typically include constraining the covariance of the outputs of the system. Recent engineering challenges have demonstrated an interest in controlling system output covariances to a setpoint. However, open-loop control of system output covariance can be exceptionally challenging when the system has significant plant uncertainties. This paper explores closed-loop control of the system’s output covariance using a windowed variance function. The stochastic nature of the windowed variance feedback means that root locus and other classical control design techniques cannot be immediately used. Expected value analysis is used to generate transfer functions for the average response of such systems, and the results are used to design controllers. Stochastic simulations are presented to verify the analytical results and make other observations concerning closed-loop variance control.

Bieniek, Mitchell↗