Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “temporal computing”

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 973 records · Page 54

Data Quality Monitoring for the Hadron Calorimeters Using Transfer Learning for Anomaly Detection

The proliferation of sensors brings an immense volume of spatio-temporal (ST) data in many domains, including monitoring, diagnostics, and prognostics applications. Data curation is a time-consuming process for a large volume of data, making it challenging and expensive to deploy data analytics platforms in new environments. Transfer learning (TL) mechanisms promise to mitigate data sparsity and model complexity by utilizing pre-trained models for a new task. Despite the triumph of TL in fields like computer vision and natural language processing, efforts on complex ST models for anomaly detection (AD) applications are limited. In this study, we present the potential of TL within the context of high-dimensional ST AD with a hybrid autoencoder architecture, incorporating convolutional, graph, and recurrent neural networks. Motivated by the need for improved model accuracy and robustness, particularly in scenarios with limited training data on systems with thousands of sensors, this research investigates the transferability of models trained on different sections of the Hadron Calorimeter of the Compact Muon Solenoid experiment at CERN. The key contributions of the study include exploring TL’s potential and limitations within the context of encoder and decoder networks, revealing insights into model initialization and training configurations that enhance performance while substantially reducing trainable parameters and mitigating data contamination effects.

47 OTHER INSTRUMENTATION↗

Distortion of Digital Image Correlation (DIC) Displacements and Strains from Heat Waves

“Heat waves” is a colloquial term used to describe convective currents in air formed when different objects in an area are at different temperatures. In the context of Digital Image Correlation (DIC) and other optical-based image processing techniques, imaging an object of interest through heat waves can significantly distort the apparent location and shape of the object. We present that there are many potential heat sources in DIC experiments, including but not limited to lights, cameras, hot ovens, and sunlight, yet error caused by heat waves is often overlooked. This paper first briefly presents three practical situations in which heat waves contributed significant error to DIC measurements to motivate the investigation of heat waves in more detail. Then the theoretical background of how light is refracted through heat waves is presented, and the effects of heat waves on displacements and strains computed from DIC are characterized in detail. Finally, different filtering methods are investigated to reduce the displacement and strain errors caused by imaging through heat waves. The overarching conclusions from this work are that errors caused by heat waves are significantly higher than typical noise floors for DIC measurements, and that the errors are difficult to filter because the temporal and spatial frequencies of the errors are in the same range as those of typical signals of interest. In conclusion, eliminating or mitigating the effects of heat sources in a DIC experiment is the best solution to minimizing errors caused by heat waves.

42 ENGINEERING↗

A low-rank control variate for multilevel Monte Carlo simulation of high-dimensional uncertain systems

Multilevel Monte Carlo (MLMC) is a recently proposed variation of Monte Carlo (MC) simulation that achieves variance reduction by simulating the governing equations on a series of spatial (or temporal) grids with increasing resolution. Instead of directly employing the fine grid solutions, MLMC estimates the expectation of the quantity of interest from the coarsest grid solutions as well as differences between each two consecutive grid solutions. When the differences corresponding to finer grids become smaller, hence less variable, fewer MC realizations of finer grid solutions are needed to compute the difference expectations, thus leading to a reduction in the overall work. Here this paper presents an extension of MLMC, referred to as multilevel control variates (MLCV), where a low-rank approximation to the solution on each grid, obtained primarily based on coarser grid solutions, is used as a control variate for estimating the expectations involved in MLMC. Cost estimates as well as numerical examples are presented to demonstrate the advantage of this new MLCV approach over the standard MLMC when the solution of interest admits a low-rank approximation and the cost of simulating finer grids grows fast.

97 MATHEMATICS AND COMPUTING↗

The influence of cloud cover on the reliability of satellite-based solar resource data

Satellite-based solar resource data are often developed and validated by using binary cloudiness categories: clear sky or overcast cloudy sky. To investigate the reliability of solar resource data in partially cloudy conditions, we estimate cloud fraction using two distinct algorithms: a physical retrieval model using surface observed global horizontal irradiance (GHI) and direct normal irradiance (DNI) and a temporal average of cloud mask data estimated by the observed DNI. Our analysis reveals a significant presence of scattered clouds, broken clouds, and mismatches between satellite- and surface-based cloud data at 17 surface sites across the contiguous United States, though confidently clear and cloudy conditions collectively account for more than 70 % of the data. Solar radiation is computed using the National Solar Radiation Database (NSRDB) algorithm and validated using surface observations. Here, our findings suggest that, in the presence of scattered clouds, NSRDB data for clear-sky conditions can be subject to significant overestimation. In cloudy-sky conditions classified by satellite data, DNI computed by the Fast All-sky Radiation Model for Solar applications with DNI (FARMS-DNI) can be underestimated when limited clouds are detected by surface observations. The bias observed in several cloudiness categories indicates that the NSRDB is exceptionally accurate in confidently clear conditions. However, clear-sky conditions with scattered clouds and mismatched cloud data contribute significantly to the overall uncertainties in the NSRDB. Therefore, future improvements in solar resource data should involve development and implementation of satellite-derived cloud fraction and should consider a novel radiative transfer model accounting for amplified cloud reflection. The evaluation within cloudiness categories also provides a physical rationale for the superior performance of FARMS-DNI compared to the Direct Insolation Simulation Code (DISC) in both cloudy-sky and all-sky conditions.

14 SOLAR ENERGY↗

Anticipating decoherence in quantum systems

Large-scale quantum technologies require coherence across distant nodes, necessitating indistinguishable quantum states. However, environmental disorder, including dephasing, spectral diffusion, and spin-bath interactions, undermines coherence. Using statistical methods, we uncover correlations in decoherence channels induced by slowly varying environments. Spectral diffusion serves as a representative demonstration case that can be extended to other remote, disordered systems such as spins in nitrogen-vacancy centers and quantum-dot spin qubits, as well as flux noise in superconducting qubits. In this work, we employ replica-theory-inspired trajectory analysis to reveal predictable temporal structures in decoherence dynamics, and validate these through an anticipatory systems framework with internal prediction of unseen spectral dynamics in multiple quantum systems, showing that this framework could, if implemented, reduce spectral shift by average factors of approximately 2 to 19, depending on emitter stability, thereby enabling enhanced coherence and multi-node synchronization for scalable quantum communication, computation, imaging, and sensing.

Maan, Pranshu [Purdue University]↗

Disruption of electrostatic contacts in the HNH nuclease from a thermophilic Cas9 rewires allosteric motions and enhances high-temperature DNA cleavage

Allosteric signaling within multidomain proteins is a driver of communication between spatially distant functional sites. Understanding the mechanism of allosteric coupling in large multidomain proteins is the most promising route to achieving spatial and temporal control of the system. The recent explosion of CRISPR-Cas9 applications in molecular biology and medicine has created a need to understand how the atomic level protein dynamics of Cas9, which are the driving force of its allosteric crosstalk, influence its biophysical characteristics. In this study, we used a synergistic approach of nuclear magnetic resonance (NMR) and computation to pinpoint an allosteric hotspot in the HNH domain of the thermostable GeoCas9. We show that mutation of K597 to alanine disrupts a salt-bridge network, which in turn alters the structure, the timescale of allosteric motions, and the thermostability of the GeoHNH domain. This homologous lysine-to-alanine mutation in the extensively studied mesophilic S. pyogenes Cas9 similarly alters the dynamics of the SpHNH domain. We have previously demonstrated that the alteration of allostery via mutations is a source for the specificity enhancement of SpCas9 (eSpCas9). Hence, this may also be true in GeoCas9.

Chemistry↗

Loss tangent fluctuations due to two-level systems in superconducting microwave resonators

Superconducting microwave resonators are critical to quantum computing and sensing technologies. Additionally, they are common proxies for superconducting qubits when determining the effects of performance-limiting loss mechanisms such as from two-level systems (TLSs). The extraction of these loss mechanisms is often performed by measuring the internal quality factor Qi as a function of power or temperature. In this work, we investigate large temporal fluctuations of Qi at low powers over periods of 12–16 h (relative standard deviation σQi/Qi=13%). These fluctuations are ubiquitous across multiple resonators, chips, and cooldowns. We are able to attribute these fluctuations to variations in the TLS loss tangent due to two main indicators. First, measured fluctuations decrease as power and temperature increase. Second, for interleaved measurements, we observe correlations between low- and medium-power Qi fluctuations and an absence of correlations with high-power fluctuations. Agreement with the TLS loss tangent mean is obtained by performing measurements over a time span of a few hours. We hypothesize that, in addition to decoherence, due to coupling to individual near-resonant TLS, superconducting qubits are affected by these observed TLS loss tangent fluctuations.

Vallières, André [Northwestern U.; Fermilab] (ORCI↗

LATERAL DOWNFLOWS IN SUNSPOT PENUMBRAL FILAMENTS AND THEIR TEMPORAL EVOLUTION

We study the temporal evolution of downflows observed at the lateral edges of penumbral filaments in a sunspot located very close to the disk center. Our analysis is based on a sequence of nearly diffraction-limited scans of the Fe i 617.3 nm line taken with the CRisp Imaging Spectro-Polarimeter instrument at the Swedish 1 m Solar Telescope. We compute Dopplergrams from the observed intensity profiles using line bisectors and filter the resulting velocity maps for subsonic oscillations. Lateral downflows appear everywhere in the center-side penumbra as small, weak patches of redshifts next to or along the edges of blueshifted flow channels. These patches have an intermittent life and undergo mergings and fragmentations quite frequently. The lateral downflows move together with the hosting filaments and react to their shape variations, very much resembling the evolution of granular convection in the quiet Sun. There is a good relation between brightness and velocity in the center-side penumbra, with downflows being darker than upflows on average, which is again reminiscent of quiet Sun convection. These results point to the existence of overturning convection in sunspot penumbrae, with elongated cells forming filaments where the flow is upward but very inclined, and weak lateral downward flows. In general, the circular polarization profiles emerging from the lateral downflows do not show sign reversals, although sometimes we detect three-lobed profiles that are suggestive of opposite magnetic polarities in the pixel.

79 ASTRONOMY AND ASTROPHYSICS↗

Transient Simulation of the Multi-SERTTA Experiment with MAMMOTH

This work details the MAMMOTH reactor physics simulations of the Static Environment Rodlet Transient Test Apparatus (SERTTA) conducted at Idaho National Laboratory in FY-2017. TREAT static-environment experiment vehicles are being developed to enable transient testing of Pressurized Water Reactor (PWR) type fuel specimens, including fuel concepts with enhanced accident tolerance (Accident Tolerant Fuels, ATF). The MAMMOTH simulations include point reactor kinetics as well as spatial dynamics for a temperature-limited transient. The strongly coupled multi-physics solutions of the neutron flux and temperature fields are second order accurate both in the spatial and temporal domains. MAMMOTH produces pellet stack powers that are within 1.5% of the Monte Carlo reference solutions. Some discrepancies between the MCNP model used in the design of the flux collars and the Serpent/MAMMOTH models lead to higher power and energy deposition values in Multi-SERTTA unit 1. The TREAT core results compare well with the safety case computed with point reactor kinetics in RELAP5-3D. The reactor period is 44 msec, which corresponds to a reactivity insertion of 2.685% delta k/k$. The peak core power in the spatial dynamics simulation is 431 MW, which the point kinetics model over-predicts by 12%. The pulse width at half the maximum power is 0.177 sec. Subtle transient effects are apparent at the beginning insertion in the experimental samples due to the control rod removal. Additional difference due to transient effects are observed in the sample powers and enthalpy. The time dependence of the power coupling factor (PCF) is calculated for the various fuel stacks of the Multi-SERTTA vehicle. Sample temperatures in excess of 3100 K, the melting point UO$_2$, are computed with the adiabatic heat transfer model. The planned shaped-transient might introduce additional effects that cannot be predicted with PRK models. Future modeling will be focused on the shaped-transient by improving the control rod models in MAMMOTH and adding the BISON thermo-elastic models and thermal-fluids heat transfer.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Closures and Simulation for Thermal Radiation Transport in Stochastic Media with Nonlinear Temperature Dependence

Because of the practical challenge of rendering very complex realistic spatial structures for numerical work, it is common practice to resort to characterizing such media as stochastic mixtures of materials, ideally parametrized with low order statistics such as the mean, variance, and correlation functions of the now random material properties. This enables realizations of the medium to be repeatedly generated and radiation transport computations to, in principle, be performed for a large ensemble of these realizations to obtain a statistically well-characterized radiation field. Statistical post-processing yields desired quantities such as conditional and unconditional mean radiation flux and probability distributions of transmitted radiation. However, such computations prove expensive for all but the simplest stochastic geometries and are most suited for benchmarking approximate models. The most common approximations lead to homogenized media so that transport computations are required only on a single medium realization but by construct provide only limited statistical information on the radiation field. Almost all approximate approaches to this problem attempt to develop equations for low order moments of the radiation intensity (mean, second moment, correlation function) but inevitably encounter a closure problem: the equation for any statistical moment will contain terms depending on unknown higher-order moments. Thus, the challenge shifts to one of developing closure relations that relate the unknown moments to the lower order moments. Under very special conditions, an exact closure can be derived but in general closures are heuristically stated constitutive relations. Also, closure approaches depend on whether the mixing statistics are spatially and/or temporally continuous as in fluctuating turbulent fields, or discontinuous as in randomly mixed solid chunks of material. Thus, unconditional averaging is generally applied in the former case but conditional averaging is more appropriate when the mixing is discontinuous. In this work, the emphasis is on binary statistical mixtures of immiscible fluids as as such the quantities of interest are averages (flux, temperature) conditioned on the material type.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

SU-E-J-60: Evaluation of Temporal Lag in Radiotherapy Gating for Tumor Motion Trajectories

Purpose: Evaluating timing differences between LINAC beam ON/OFF and the estimation of tumor positioning using gating systems is essential for establishing confidence when treating with gating during radiotherapy, and is an annual requirement of TG-142. Temporal discrepancies between the trajectories of external marker surrogates and beam delivery may vary depending upon the type of external marker motion, which is quantified in this work for several trajectories. Methods: A precise robotic 3D motion stage performed several trajectories typically used for gating phantoms, including sinusoidal and Lujan-type motion; a commercial respiratory motion simulator was also employed. The true motions were monitored using variable resistors. The beam ON/OFF was controlled separately by two RPM (Varian) systems, an integrated version delivered by a Varian Truebeam LINAC and version 1.6 delivered by a Varian Trilogy, and measured using a diode. The resistor and diode signals were read by a multichannel digital oscilloscope, and timing differences between beam ON/OFF as measured by the diode and the phantom motion were determined using a peak detection algorithm. Results: Timing differences between beam ON/OFF and 3D stage motion peaks (diode—true motion timing) were computed to be 79.4 & 57.7ms for sinusoidal motion and 109.1 & 63.6ms for Lujan-type motion on the Truebeam LINAC, for beam ON and OFF, respectively. Timing differences for the Trilogy LINAC were 34.4 & 55.2ms for the sinusoidal motion and 29.0 & 26.3ms for the Lujan-type motion, for beam ON and OFF, respectively. With the commercial motion simulator, the timing differences were found to be −9ms and −78ms for beam ON/OFF, respectively, with the Truebeam, and −97.6ms and −60.9ms for beam ON/OFF, respectively, with the Trilogy. Conclusion: Setup-dependent temporal lags were found using this methodology. These discrepancies have the potential to influence quality assurance on gating systems and ultimately the treatment dose during radiotherapy. NIH National Institute of Biomedical Imaging and Bioengineering grant T32 EB002103-21.

60 APPLIED LIFE SCIENCES↗

Data-Driven Whole-Genome Clustering to Detect Geospatial, Temporal, and Functional Trends in SARS-CoV-2 Evolution

Current methods for defining SARS-CoV-2 lineages ignore the vast majority of the SARS-CoV-2 genome. We develop and apply an exhaustive vector comparison method that directly compares all known SARS-CoV-2 genome sequences to produce novel lineage classifications. We utilize data-driven models that (i) accurately capture the complex interactions across the set of all known SARS-CoV-2 genomes, (ii) scale to leadership-class computing systems, and (iii) enable tracking how such strains evolve geospatially over time. We show that during the height of the original Omicron surge, countries across Europe, Asia, and the Americas had a spatially asynchronous distribution of Omicron sub-strains. Moreover, neighboring countries were often dominated by either different clusters of the same variant or different variants altogether throughout the pandemic. Analyses of this kind may suggest a different pattern of epidemiological risk than was understood from conventional data, as well as produce actionable insights and transform our ability to prepare for and respond to current and future biological threats.

Jacobson, Daniel↗

Scalable subsurface inverse modeling of huge data sets with an application to tracer concentration breakthrough data from magnetic resonance imaging

When characterizing subsurface properties is crucial for reliable and cost-effective groundwater supply management and contaminant remediation. With recent advances in sensor technology, large volumes of hydro-geophysical and geochemical data can be obtained to achieve high-resolution images of subsurface properties. However, characterization with such a large amount of information requires prohibitive computational costs associated with “big data” processing and numerous large-scale numerical simulations. To tackle such difficulties, the Principal Component Geostatistical Approach (PCGA) has been proposed as a “Jacobian-free” inversion method that requires much smaller forward simulation runs for each iteration than the number of unknown parameters and measurements needed in the traditional inversion methods. PCGA can be conveniently linked to any multi-physics simulation software with independent parallel executions. In our paper, we extend PCGA to handle a large number of measurements (e.g. 106 or more) by constructing a fast preconditioner whose computational cost scales linearly with the data size. For illustration, we characterize the heterogeneous hydraulic conductivity (K) distribution in a laboratory-scale 3-D sand box using about 6 million transient tracer concentration measurements obtained using magnetic resonance imaging. Since each individual observation has little information on the K distribution, the data was compressed by the zero-th temporal moment of breakthrough curves, which is equivalent to the mean travel time under the experimental setting. Moreover, only about 2,000 forward simulations in total were required to obtain the best estimate with corresponding estimation uncertainty, and the estimated K field captured key patterns of the original packing design, showing the efficiency and effectiveness of the proposed method. This article is protected by copyright. All rights reserved.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Brownian dynamics without Green's functions

We develop a Fluctuating Immersed Boundary (FIB) method for performing Brownian dynamics simulations of confined particle suspensions. Unlike traditional methods which employ analytical Green's functions for Stokes flow in the confined geometry, the FIB method uses a fluctuating finite-volume Stokes solver to generate the action of the response functions “on the fly.” Importantly, we demonstrate that both the deterministic terms necessary to capture the hydrodynamic interactions among the suspended particles, as well as the stochastic terms necessary to generate the hydrodynamically correlated Brownian motion, can be generated by solving the steady Stokes equations numerically only once per time step. This is accomplished by including a stochastic contribution to the stress tensor in the fluid equations consistent with fluctuating hydrodynamics. We develop novel temporal integrators that account for the multiplicative nature of the noise in the equations of Brownian dynamics and the strong dependence of the mobility on the configuration for confined systems. Notably, we propose a random finite difference approach to approximating the stochastic drift proportional to the divergence of the configuration-dependent mobility matrix. Through comparisons with analytical and existing computational results, we numerically demonstrate the ability of the FIB method to accurately capture both the static (equilibrium) and dynamic properties of interacting particles in flow.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Faster network disruption from layered oscillatory dynamics

Nonlinear complex network-coupled systems typically have multiple stable equilibrium states. Following perturbations or due to ambient noise, the system is pushed away from its initial equilibrium, and, depending on the direction and the amplitude of the excursion, it might undergo a transition to another equilibrium. It was recently demonstrated [M. Tyloo, J. Phys. Complex. 3 03LT01 (2022)] that layered complex networks may exhibit amplified fluctuations. Here, I investigate how noise with system-specific correlations impacts the first escape time of nonlinearly coupled oscillators. Interestingly, I show that, not only the strong amplification of the fluctuations is a threat to the good functioning of the network but also the spatial and temporal correlations of the noise along the lowest-lying eigenmodes of the Laplacian matrix. Finally, I analyze first escape times on synthetic networks and compare noise originating from layered dynamics to uncorrelated noise.

97 MATHEMATICS AND COMPUTING↗

Numerical Methodology for Optimization of Compression-Ignited Engines Considering Combustion Noise Control

It is challenging to develop highly efficient and clean engines while meeting user expectations in terms of performance, comfort, and drivability. One of the critical aspects in this regard is combustion noise control. Combustion noise accounts for about 40 percent of the overall engine noise in typical turbocharged diesel engines. The experimental investigation of noise generation is difficult due to its inherent complexity and measurement limitations. Therefore, it is important to develop efficient numerical strategies in order to gain a better understanding of the combustion noise mechanisms. In this work, a novel methodology was developed, combining computational fluid dynamics (CFD) modeling and genetic algorithm (GA) technique to optimize the combustion system hardware design of a high-speed direct injection (HSDI) diesel engine, with respect to various emissions and performance targets including combustion noise. The CFD model was specifically set up to reproduce the unsteady pressure field inside the combustion chamber, thereby allowing an accurate prediction of the acoustic response of the combustion phenomena. The model was validated by simulating several steady operating conditions and comparing the numerical results against experimental data, in both temporal and frequency domains. Thereafter, a GA optimization was performed with the goal of minimizing indicated specific fuel consumption (ISFC) and combustion noise, while restricting pollutant (soot and NOx) emissions to their respective baseline values. Eight design variables were selected pertaining to piston bowl geometry, nozzle inclusion angle, number of injector nozzle holes, and in-cylinder swirl. An objective merit function (MF) based on the emissions, ISFC, and combustion noise was constructed to quantify the strength of the engine designs and was determined using the CFD model as the function evaluator. The in-cylinder noise level was characterized by the total resonance energy of local pressure oscillations. The optimum engine configuration thus obtained showed a significant improvement in terms of efficiency and combustion noise compared to the baseline system, along with both soot and NOx emissions within their respective constraints. This optimum configuration included a deeper and tighter bowl geometry with higher swirl and larger number of nozzle holes. Subsequently, a more detailed acoustics analysis based on proper orthogonal decomposition (POD) technique was carried out to further explore the combustion noise benefits achieved by the GA optimum. In conclusion, this computational study is a first of its kind (to the best of the authors’ knowledge), which demonstrates a comprehensive framework to incorporate combustion noise into a numerical optimization strategy for engine design.

42 ENGINEERING↗

Quantifying the atomistic free-volume morphology of materials with graph theory

Here, we introduce a new computational methodology for the identification and characterization of free volume within/around atomistic configurations. This scheme employs a three-stage workflow, by which spheres are iteratively grown inside of voxels, and ultimately converted to planar graphs, which are then characterized via a graph-based order parameter. Our approach is computationally efficient, physically intuitive, and universally transferable to any material system. Validation of our methodology is performed on several sets of materials problems: (1) classification of unique free volumes in various crystal phases, (2) autonomous detection and classification of complex surface defects during epitaxial growth simulations, (3) characterization of free volume defects in metals/alloys, and (4) quantification of the spatio-temporal behavior of nano-scale free volume morphologies as a function of both temperature and free-volume size. Our method accurately identifies and characterizes unique free volumes over a multitude of systems and length scales, indicating its potential for future use in understanding the relationship between free volume morphology and material properties under both static and dynamic conditions.

36 MATERIALS SCIENCE↗