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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Analytical Modeling of Metal Foam Composite Phase Change Materials (PCM) in Thermal Energy Storage Using Asymptotic Analysis

The use of phase change materials (PCMs) for thermal energy storage can release or absorb a significant amount of latent heat during the freezing or melting process, offering a higher energy storage density. One of the main drawbacks of PCMs is their low thermal conductivity, resulting in poor thermal performance. Recent research has attempted to enhance heat transfer and increase the thermal conductivity of PCMs, including the use of metal foams. However, modeling the metal foam composite PCM using conventional methods is computationally expensive. This paper proposes an asymptotic solution for a Stefan-like problem subject to a convective boundary for outward solidification in a hollow cylinder, capable of predicting the freeze-melt cycle of the metal foam composite PCM. Specifically, three temporal regimes and four spatial layers are considered in the asymptotic analysis for each phase change process. The thermal conductivity is calculated by a theoretical three-dimensional tetrakaidecahedron model, while other thermophysical properties are obtained using the method of volume averaging. The results are verified with numerical data and validated against experimental data in the literature. The presented analytical modeling framework could have the potential to be applied to other types of composite PCMs with considerably lower computational costs compared with conventional methods.

analytical model↗
DOE OSTI2025

Efficiently Computable Limits on EPR Pair Generation in Quantum Broadcast Channels

We investigate the generation of EPR pairs between three observers in a general causally structured setting, where communication occurs via a noisy quantum broadcast channel. The most general quantum codes for this setup take the form of tripartite quantum channels. Since the receivers are constrained by causal ordering, additional temporal relationships naturally emerge between the parties. These causal constraints enforce intrinsic no-signalling conditions on any tripartite operation, ensuring that it constitutes a physically realizable quantum code for a quantum broadcast channel. We analyze these constraints and, more broadly, characterize the most general quantum codes for communication over such channels. We examine the capabilities of codes that are fully no-signalling among the three parties, positive partial transpose (PPT)-preserving, or both, and derive simple semidefinite programs to compute the achievable entanglement fidelity. We then establish a hierarchy of semidefinite programming converse bounds -- both weak and strong -- for the capacity of quantum broadcast channels for EPR pair generation, in both one-shot and asymptotic regimes. Notably, in the special case of a point-to-point channel, our strong converse bound recovers and strengthens existing results. Finally, we demonstrate how the PPT-preserving codes we develop can be leveraged to construct PPT-preserving entanglement combing schemes, and vice versa.

FOS: Physical sciences↗
DOE OSTI2015

Tungsten dust impact on ITER-like plasma edge

The impact of tungsten dust originating from divertor plates on the performance of edge plasma in ITER-like discharge is evaluated using computer modeling with the coupled dust-plasma transport code DUSTT-UEDGE. Different dust injection parameters, including dust size and mass injection rates, are surveyed. It is found that tungsten dust injection with rates as low as a few mg/s can lead to dangerously high tungsten impurity concentrations in the plasma core. Dust injections with rates of a few tens of mg/s are shown to have a significant effect on edge plasma parameters and dynamics in ITER scale tokamaks. The large impact of certain phenomena, such as dust shielding by an ablation cloud and the thermal force on tungsten ions, on dust/impurity transport in edge plasma and consequently on core tungsten contamination level is demonstrated. It is also found that high-Z impurities provided by dust can induce macroscopic self-sustained plasma oscillations in plasma edge leading to large temporal variations of edge plasma parameters and heat load to divertor target plates.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Vortex breakdown in closed containers with polygonal cross sections

The vortex breakdown bubble in the confined flow generated by a rotating lid in closed containers with polygonal cross sections was analysed both experimentally and numerically for the height/radius aspect ratio equal to 2. The stagnation point locations of the breakdown bubble emergence and the corresponding Reynolds number were determined experimentally and in addition computed numerically by STAR-CCM+ CFD software for square, pentagonal, hexagonal, and octagonal cross section configurations. The flow pattern and the velocity were observed and measured by combining the seeding particle visualization and the temporal accuracy of laser Doppler anemometry. The vortex breakdown size and position on the container axis were determined for Reynolds numbers, ranging from 1450 to 2400. The obtained results were compared with the flow structure in the closed container of cubical and cylindrical configurations. It is shown that the measured evolution of steady vortex breakdown is in close agreement with the numerical results.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Quantum Electromagnetic Transients Program

This work devises a quantum electromagnetic transients program (QEMTP) which is the first attempt to tackle the computational challenges in solving EMTP through quantum computing. The main contributions lie in: (1) A quantum-enabled EMTP formulation with Dommel's model encoded in quantum states; 2) An Harrow-Hassidim-Lloyd (HHL)-based QEMTP to solve the discrete-time nodal equations; 3) An iteration-based HHL revision for mitigating temporal errors to achieve high accuracy of QEMTP with limited quantum resources. Case studies verify the correctness and efficacy of QEMTP in simulating both fast and slow dynamics.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2015

SU-E-J-113: Effects of Deformable Registration On First-Order Texture Maps Calculated From Thoracic Lung CT Scans

Purpose: To determine the stability of eight first-order texture features following the deformable registration of serial computed tomography (CT) scans. Methods: CT scans at two different time points from 10 patients deemed to have no lung abnormalities by a radiologist were collected. Following lung segmentation using an in-house program, texture maps were calculated from 32×32-pixel regions of interest centered at every pixel in the lungs. The texture feature value of the ROI was assigned to the center pixel of the ROI in the corresponding location of the texture map. Pixels in the square ROI not contained within the segmented lung were not included in the calculation. To quantify the agreement between ROI texture features in corresponding pixels of the baseline and follow-up texture maps, the Fraunhofer MEVIS EMPIRE10 deformable registration algorithm was used to register the baseline and follow-up scans. Bland-Altman analysis was used to compare registered scan pairs by computing normalized bias (nBias), defined as the feature value change normalized to the mean feature value, and normalized range of agreement (nRoA), defined as the range spanned by the 95% limits of agreement normalized to the mean feature value. Results: Each patient’s scans contained between 6.8–15.4 million ROIs. All of the first-order features investigated were found to have an nBias value less than 0.04% and an nRoA less than 19%, indicating that the variability introduced by deformable registration was low. Conclusion: The eight first-order features investigated were found to be registration stable. Changes in CT texture maps could allow for temporal-spatial evaluation of the evolution of lung abnormalities relating to a variety of diseases on a patient-by-patient basis. SGA and HA receives royalties and licensing fees through the University of Chicago for computer-aided diagnosis technology. Research reported in this publication was supported by the National Institute Of General Medical Sciences of the National Institutes of Health under Award Number R25GM109439.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2018

A hydrological emulator for global applications – HE v1.0.0

Abstract. While global hydrological models (GHMs) are very useful in exploring water resources and interactions between the Earth and human systems, their use often requires numerous model inputs, complex model calibration, and high computation costs. To overcome these challenges, we construct an efficient open-source and ready-to-use hydrological emulator (HE) that can mimic complex GHMs at a range of spatial scales (e.g., basin, region, globe). More specifically, we construct both a lumped and a distributed scheme of the HE based on the monthly abcd model to explore the tradeoff between computational cost and model fidelity. Model predictability and computational efficiency are evaluated in simulating global runoff from 1971 to 2010 with both the lumped and distributed schemes. The results are compared against the runoff product from the widely used Variable Infiltration Capacity (VIC) model. Our evaluation indicates that the lumped and distributed schemes present comparable results regarding annual total quantity, spatial pattern, and temporal variation of the major water fluxes (e.g., total runoff, evapotranspiration) across the global 235 basins (e.g., correlation coefficient r between the annual total runoff from either of these two schemes and the VIC is >0.96), except for several cold (e.g., Arctic, interior Tibet), dry (e.g., North Africa) and mountainous (e.g., Argentina) regions. Compared against the monthly total runoff product from the VIC (aggregated from daily runoff), the global mean Kling–Gupta efficiencies are 0.75 and 0.79 for the lumped and distributed schemes, respectively, with the distributed scheme better capturing spatial heterogeneity. Notably, the computation efficiency of the lumped scheme is 2 orders of magnitude higher than the distributed one and 7 orders more efficient than the VIC model. A case study of uncertainty analysis for the world's 16 basins with top annual streamflow is conducted using 100 000 model simulations, and it demonstrates the lumped scheme's extraordinary advantage in computational efficiency. Our results suggest that the revised lumped abcd model can serve as an efficient and reasonable HE for complex GHMs and is suitable for broad practical use, and the distributed scheme is also an efficient alternative if spatial heterogeneity is of more interest.

Liu, Yaling↗
DOE OSTI2024

Multitiered computational methodology for extracting three-dimensional rotational diffusion coefficients from x-ray photon correlation spectroscopy data without structural information

X-ray photon correlation spectroscopy (XPCS) is a powerful technique for analyzing particle systems by investigating their dynamics in suspensions across a broad range of temporal and spatial scales. This is done by illuminating samples with coherent x-ray beams and calculating the correlation function of the obtained x-ray scattering images. XPCS is uniquely suited for studying Brownian dynamics, consisting of translational and rotational diffusion. While traditional XPCS image analysis techniques can extract translational diffusion components, they are unable to estimate rotational diffusion coefficients. Here, we introduce a methodology that combines the angular-temporal cross-correlation analysis and a algorithmic framework called Multi-Tiered Estimation for Correlation Spectroscopy in 3D for estimating three-dimensional rotational diffusion coefficients from XPCS images of three-dimensional particle systems. We demonstrate our methodology for extracting rotational diffusion coefficients from XPCS data by applying it to simulated noisy x-ray images of systems of crossing nanotubes and proteins that evolve under translational and rotational Brownian motion for different diffusion rates. Furthermore, our results show that our approach determines rotational diffusion coefficients within a few percent error.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Cross-scale modeling of storm surge, tide, and inundation in Mid-Atlantic Bight and New York City during Hurricane Sandy, 2012

Driven by high-resolution NAM (North America Mesoscale) and ECMWF (European Centre for Medium-Range Weather Forecasts) atmospheric fields, a 3-D unstructured grid SCHISM (Semi-implicit Cross-scale Hydroscience Integrated System Model) was applied to simulate total water level in the Mid-Atlantic Bight during Hurricane Sandy in 2012. The simulated storm surge, tide, significant wave height, and wave peak period results compared favorably with NOAA observations along US East Coast, Long Island Sound, and New York Harbor. The maximum total water level at The Battery in the New York Harbor was accurately simulated with an absolute error of less than 0.08 m (out of 2.9 m peak surge) and a timing difference less than 10 min, which includes the effects on the order of 0.1–0.3 m (5–10%) from the coastal setup induced by the surface waves. The scenario comparisons of (1) “NAM” versus “ECMWF”, (2) “2-D” versus “3-D”, and (3) “with” versus “without” wind wave model were examined. The 3-D barotropic model forced by ECMWF including the effects of wind waves performs the best, attributed to wave-induced radiation stress and reduction of bottom stress when the 3-D version is used. Simultaneously, the ELCIRC-Sub model, using a 5-m topography/bathymetry sub-grid and a regular 200-m resolution finite volume computational model grid, was developed to simulate the street-level inundation in New York City. The momentum and mass fluxes calculated by the coarser base grid model were effectively coupled with the sub-grid so that running a full-blown high-resolution base model is not required. The ELCIRC-Sub model uses MPI (Message Passing Interface) parallel computing to enlarge the coverage of the land surface and an efficient non-linear solver to improve the accuracy of the wetting-and-drying processes. The temporal comparisons of modeled water level with NOAA's tidal stations and USGS' rapid-deployed gauges showed overall performance with an average error of 0.1 m. Particularly, ELCIRC-Sub captured the profile of the highest peak surge (3.9 m), which rose 2 m within 3 hour at Kings Point, NY. The spatial comparisons of the modeled peak water level at 80 surveyed locations showed an average error less than 0.13 m. The modeled maximum inundation extent also matched well with an 80% hit rate against FEMA's Hurricane Sandy maximum flooding extent.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

An exponential time-integrator scheme for steady and unsteady inviscid flows

In this report an exponential time-integrator scheme of second-order accuracy based on the predictor–corrector methodology, denoted PCEXP, is developed to solve multi-dimensional nonlinear partial differential equations pertaining to fluid dynamics. The effective and efficient implementation of PCEXP is realized by means of the Krylov method. The linear stability and truncation error are analyzed through a one-dimensional model equation. The proposed PCEXP scheme is applied to the Euler equations discretized with a discontinuous Galerkin method in both two and three dimensions. The effectiveness and efficiency of the PCEXP scheme are demonstrated for both steady and unsteady inviscid flows. The accuracy and efficiency of the PCEXP scheme are verified and validated through comparisons with the explicit third-order total variation diminishing Runge–Kutta scheme (TVDRK3), the implicit backward Euler (BE) and the implicit second-order backward difference formula (BDF2). For unsteady flows, the PCEXP scheme generates a temporal error much smaller than the BDF2 scheme does, while maintaining the expected acceleration at the same time. Moreover, the PCEXP scheme is also shown to achieve the computational efficiency comparable to the implicit schemes for steady flows.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Material Fracturing and Failure Simulation Datasets

Fracturing is a fundamental physics phenomena with broad relevance across multiple domains, ranging from infrastructure integrity, aerospace durability, reservoir production, and seismic events. We present a diverse dataset of simulated fracture evolution and material failure generated from two numerical solvers: the phase-field method and the combined finite-discrete element method (FDEM). These solvers differ in formulation, physical fidelity, and computational efficiency. The dataset includes five materials: PBX, anisotropic shale, tungsten, aluminum, and steel. For each, phase-field simulations span 400,000 cases: 200,000 under uniaxial tension and 200,000 under biaxial tension. The computationally expensive FDEM simulations include 90,000 split evenly among PBX, shale, and tungsten under uniaxial loading. All simulations begin with randomized initial fracture patterns. Each entry includes temporal data capturing fracture propagation dynamics. This comprehensive dataset is designed to support the development of foundational or surrogate machine learning approaches for predicting material failure. While no such models are introduced here, the dataset lays a robust foundation for advancing future research and innovation in these areas.

36 MATERIALS SCIENCE↗
DOE OSTI2015

A new model for the radio emission from SN 1994I and an associated search for radio transients in M51

We revisit the exquisite archival radio data for the Type Ic supernova SN 1994I and present a revised model for the supernova (SN) radio emission and a pilot study that aims to constrain the rate of C-band radio transients within the face-on host galaxy, M51 (NGC 5194). We find that the temporal and spectral evolution of the SN 1994I radio emission are well fit by a synchrotron self-absorption model and use this to estimate physical parameters. We compute a pre-explosion mass loss rate of M-dot =3.0×10{sup −5} M{sub ⊙} yr{sup −1} for the progenitor, consistent with those observed from galactic Wolf–Rayet stars. Our model makes different assumptions for the dynamical model for the shockwave interaction than the model previously published by Weiler et al., but our M-dot is consistent with theirs to within errors and assumptions. Drawing from a subset of the archival radio observations from the Very Large Array collected for the monitoring of SN 1994I, we conduct a pilot study to search for previously unidentified transients. Data were primarily taken at a frequency of 4.9 GHz and are logarithmic in cadence, enabling sensitivity to transients with variability timescales ranging from days to months. We find no new transient detections in 31 epochs of data, allowing us to place a 2σ upper limit of 17 deg{sup −2} for the source density of radio transients above 0.5 mJy (L ≳ 4×10{sup 25} erg s{sup −1} Hz{sup −1} at the distance of M51). This study highlights the feasibility of utilizing archival high-cadence radio studies of SN host galaxies to place constraints on the radio transient rate as a function of luminosity in the local universe.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Efficient Uncertainty Quantification in Stochastic Economic Dispatch

Stochastic economic dispatch models usually address uncertainties in forecasts of renewable generation output by considering a finite number of realizations drawn from a stochastic process model, typically via Monte Carlo sampling. Accurate evaluations of expectations or higher order moments for quantities of interest, e.g., generating cost, can require a prohibitively large number of samples. We propose an alternative to Monte Carlo sampling based on polynomial chaos expansions. These representations enable efficient and accurate propagation of uncertainties in model parameters, using sparse quadrature methods. Furthermore, we use Karhunen-Loève expansions for efficient representation of uncertain renewable energy generation that follows geographical and temporal correlations derived from historical data at each wind farm. Considering expected production cost, we show that the proposed approach can yield several orders of magnitude reduction in computational cost for solving stochastic economic dispatch relative to Monte Carlo sampling, for a given target error threshold.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

Dynamic sparse x-ray nanotomography reveals ionomer hydration mechanism in polymer electrolyte fuel-cell catalyst

Tomographic imaging of time-evolving samples is a challenging yet important task for various research fields. At the nanoscale, current approaches face limitations of measurement speed or resolution due to lengthy acquisitions. We developed a dynamic nanotomography technique based on sparse dynamic imaging and 4D tomography modeling. We demonstrated the technique, using ptychographic x-ray computed tomography as its imaging modality, on resolving the in situ hydration process of polymer electrolyte fuel cell (PEFC) catalyst. The technique provides a 40-time increase in temporal resolution compared to conventional approaches, yielding 28 nm half-period spatial and 12 min temporal resolution. The results allow a quantitative characterization of the water intake process inside PEFC catalysts with nanoscale resolution, which is crucial for understanding their electrochemical mechanisms and optimizing their performance. Our technique enables high-speed operando nanotomography studies and paves the way for wider application of dynamic tomography at the nanoscale.

Science & Technology - Other Topics↗
DOE OSTI2025

PaleoSTeHM v1.0: a modern, scalable spatiotemporal hierarchical modeling framework for paleo-environmental data

Abstract. Geological records of past environmental change provide crucial insights into long-term climate variability, trends, non-stationarity, and nonlinear feedback mechanisms. However, reconstructing spatiotemporal fields from these records is statistically challenging due to their sparse, indirect, and noisy nature. Here, we present PaleoSTeHM, a scalable and modern framework for spatiotemporal hierarchical modeling of paleo-environmental data. This framework enables the implementation of flexible statistical models that rigorously quantify spatial and temporal variability from geological data while clearly distinguishing measurement and inferential uncertainty from process variability. We illustrate its application by reconstructing temporal and spatiotemporal paleo-sea-level changes across multiple locations. Using various modeling and analysis choices, PaleoSTeHM demonstrates the impact of different methods on inference results and computational efficiency. Our results highlight the critical role of model selection in addressing specific paleo-environmental questions, showcasing the PaleoSTeHM framework's potential to enhance the robustness and transparency of paleo-environmental reconstructions.

58 GEOSCIENCES↗
DOE OSTI2016

Fiber Optic Distributed Sensors for High-resolution Temperature Field Mapping

The reliability of computational fluid dynamics (CFD) codes is checked by comparing simulations with experimental data. A typical data set consists chiefly of velocity and temperature readings, both ideally having high spatial and temporal resolution to facilitate rigorous code validation. While high resolution velocity data is readily obtained through optical measurement techniques such as particle image velocimetry, it has proven difficult to obtain temperature data with similar resolution. Traditional sensors such as thermocouples cannot fill this role, but the recent development of distributed sensing based on Rayleigh scattering and swept-wave interferometry offers resolution suitable for CFD code validation work. Thousands of temperature measurements can be generated along a single thin optical fiber at hundreds of Hertz. Sensors function over large temperature ranges and within opaque fluids where optical techniques are unsuitable. But this type of sensor is sensitive to strain and humidity as well as temperature and so accuracy is affected by handling, vibration, and shifts in relative humidity. Such behavior is quite unlike traditional sensors and so unconventional installation and operating procedures are necessary to ensure accurate measurements. This paper demonstrates implementation of a Rayleigh scattering-type distributed temperature sensor in a thermal mixing experiment involving two air jets at 25 and 45 °C. We present criteria to guide selection of optical fiber for the sensor and describe installation setup for a jet mixing experiment. We illustrate sensor baselining, which links readings to an absolute temperature standard, and discuss practical issues such as errors due to flow-induced vibration. This material can aid those interested in temperature measurements having high data density and bandwidth for fluid dynamics experiments and similar applications. We highlight pitfalls specific to these sensors for consideration in experiment design and operation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2017

A historical survey of algorithms and hardware architectures for neural-inspired and neuromorphic computing applications

In this study, biological neural networks continue to inspire new developments in algorithms and microelectronic hardware to solve challenging data processing and classification problems. Here in this research, we survey the history of neural-inspired and neuromorphic computing in order to examine the complex and intertwined trajectories of the mathematical theory and hardware developed in this field. Early research focused on adapting existing hardware to emulate the pattern recognition capabilities of living organisms. Contributions from psychologists, mathematicians, engineers, neuroscientists, and other professions were crucial to maturing the field from narrowly-tailored demonstrations to more generalizable systems capable of addressing difficult problem classes such as object detection and speech recognition. Algorithms that leverage fundamental principles found in neuroscience such as hierarchical structure, temporal integration, and robustness to error have been developed, and some of these approaches are achieving world-leading performance on particular data classification tasks. Additionally, novel microelectronic hardware is being developed to perform logic and to serve as memory in neuromorphic computing systems with optimized system integration and improved energy efficiency. Key to such advancements was the incorporation of new discoveries in neuroscience research, the transition away from strict structural replication and towards the functional replication of neural systems, and the use of mathematical theory frameworks to guide algorithm and hardware developments.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Towards a novel stellar opacity measurement scheme using stability properties of double ablation front structures

A novel method for bringing sample elements to hydrodynamic conditions relevant to the base of the solar convection zone is investigated. The method is designed in the framework of opacity measurements and exploits the temporal and spatial stability of hydrodynamic parameters in counter-propagating Double Ablation Front (DAF) structures. The physics of symmetric DAF structures is first studied in 1D geometries to assess the influence of tracer layers in the target. These results are used to propose an experimental design compatible with the OMEGA [Boehly et al., Opt. Commun. 133(1-6), 495–506 (1997)] laser. Radiative-hydrodynamic simulations conducted using the Chic code [Breil et al., Comput. Fluids 46, 161–167 (2011).] in 2D-axisymmetric geometries suggest that a Fe sample can be brought to an electron temperature of ∼160 eV and electron number density of ∼1.35 × 1023 cm−3. These parameters are reached during a 500 ps window with temporal variations of the order of 10 eV and 1022 cm−3, respectively. This allows for potential time-integrated spectral measurements. During that time, the sample is almost at local thermal equilibrium and 2D spatial gradients in the sample are less than 5% in a 360 μm diameter cylindrical volume, including the potential effects of Hot Electrons (HE) and typical uncertainties related to target fabrication and laser performances. The effects of HEs are assessed using an inline model in Chic. The HEs are found to deposit most of their energy in the cold and dense ablator between the two fronts, leading to a small efficiency loss on the DAF parameters. The calculations also suggest that negligible amounts of unabsorbed HEs are present in the probed volume, thus not affecting the atomic properties of the sample. Potential extensions of the current design to higher sample temperatures within the OMEGA capabilities are briefly discussed.

Colaïtis, A.↗