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 991 records · Page 55
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

Upgrade and benchmarking of a 4D treatment planning system for scanned ion beam therapy

Purpose: Upgrade and benchmarking of a research 4D treatment planning system (4DTPS) suitable for realistic patient treatment planning and treatment simulations taking into account specific requirements for scanned ion beam therapy, i.e., modeling of dose heterogeneities due to interplay effects and range changes caused by patient motion and dynamic beam delivery. Methods: The 4DTPS integrates data interfaces to 4D computed tomography (4DCT), deformable image registration and clinically used motion monitoring devices. The authors implemented a novel data model for 4D image segmentation using Boolean mask volume datasets and developed an algorithm propagating a manually contoured reference contour dataset to all 4DCT phases. They further included detailed treatment simulation and dose reconstruction functionality, based on the irregular patient motion and the temporal structure of the beam delivery. The treatment simulation functionality was validated against experimental data from irradiation of moving radiographic films in air, 3D moving ionization chambers in a water phantom, and moving cells in a biological phantom with a scanned carbon ion beam. The performance of the program was compared to results obtained with predecessor programs. Results: The measured optical density distributions of the radiographic films were reproduced by the simulations to (-2 {+-} 12)%. Compared to earlier versions of the 4DTPS, the mean agreement improved by 2%, standard deviations were reduced by 7%. The simulated dose to the moving ionization chambers in water showed an agreement with the measured dose of (-1 {+-} 4)% for the typical beam configuration. The mean deviation of the simulated from the measured biologically effective dose determined via cell survival was (617 {+-} 538) mGy relative biological effectiveness corresponding to (10 {+-} 9)%. Conclusions: The authors developed a research 4DTPS suitable for realistic treatment planning on patient data and capable of simulating dose delivery to a moving patient geometry for scanned ion beams. The accuracy and reliability of treatment simulations improved considerably with respect to earlier versions of the 4DTPS.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

A Bayesian modeling framework for crash severity effects of active traffic management systems

Transportation agencies utilize Active traffic management (ATM) systems to dynamically manage recurrent and non-recurrent congestion based on real-time conditions. While these systems have been shown to have some safety benefits, their impact on injury severity outcomes is currently uncertain. In this paper, we used full Bayesian mixed logit models to quantify the impact that ATM deployment had on crash severities. The estimation results revealed lower severities with ATM deployment. Marginal effects for ATM deployments that featured hard shoulder running (HSR) revealed lower likelihoods for severe and moderate injury crashes of 15.9 % and for minor injury crashes of 10.1 %. The likelihood of severe and moderate injury crashes and minor injury crashes reduced by 12.4 % and 8.33 % with ATM without HSR. The models were observed to be temporally transferable and had forecast error of 0.301 and 0.304 for the two models, revealing better performance with validation data. These results have implications for improving freeway crash risk at critical locations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Dust Polarization Maps from TIGRESS: E/B Power Asymmetry and TE Correlation

We present the first large set of all-sky synthetic dust polarization maps derived directly from a self-consistent magnetohydrodynamics simulation using the TIGRESS framework. Turbulence in this simulation is predominantly driven by supernova explosions, with rates that are self-consistently regulated by feedback loops. The simulation covers both the outer scale and inertial range of turbulence with uniformly high resolution. The shearing-box utilized in the simulation, in concert with resolved supernova-driven turbulence, enables the capturing of generation, growth, and saturation of both turbulent and mean magnetic fields. We construct polarization maps at 353 GHz, as seen by observers inside a model of the multiphase, turbulent, magnetized interstellar medium (ISM). To fully sample the simulated ISM state, we use 350 snapshots spanning over $\sim 350\,\mathrm{Myr}$ (more than six feedback loops) and nine representative observers. The synthetic skies show a prevalent E/B power asymmetry (${EE}\gt {BB}$) and positive TE correlation in broad agreement with observations by the Planck satellite. However, the ranges of ${EE}/{BB}\sim 1.4\mbox{-}1.7$ and ${TE}/{({TT}\cdot {EE})}^{1/2}\sim 0.2\mbox{-}0.3$ are generally lower than those measured by Planck. We find large fluctuations of E/B asymmetry and TE correlation depending on the observer's position and temporal fluctuations of ISM properties due to bursts of star formation. The synthetic maps are made publicly available to provide novel models of the microwave sky.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

A geospatial model for the analysis of time-dependent land subsidence induced by reservoir depletion

Land subsidence due to fluid depletion is an outcome of physical processes operating across a wide range of time and length scales. Although geomechanical models are crucial to simulate reservoir compaction and predict its long-term fate, their use across large regions often bears prohibitive computational costs. To overcome this obstacle, this paper proposes a simplified modelling framework consisting of (i) a near-field numerical solver simulating the coupling between fluid flow and rock deformation with reference to a simplified one-dimensional geometry and (ii) a far-field geospatial algorithm mapping ground settlements across a region through the superposition of poroelastic computations at multiple wells. The model computes the delay between depletion history and reservoir compaction in proximity of a producing well by assuming basal depletion of a fluid-saturated deformable disk, while the Geertsma solution of nucleus of strain is used to extrapolate the impact of such time-varying reservoir compaction around the well. This approach has been used to back-analyze the spatio-temporal progression of subsidence at the Groningen gas field. The results are compared against measurements collected over 50 year of production at 25 benchmark locations scattered over an area of 900 km 2 . It is shown that coupled simulations based on average values of rock compressibility and permeability lead to nonlinear trends of subsidence evolution in good agreement with field measurements, while uncoupled analyses overpredict settlements by more than 70%. Lastly, synthetic forecasts based on different rates of depletion were provided. The results suggest that slower depletion rates lead to lower subsidence at a given time, and that residual subsidence may continue to develop for several decades after interruption of production activities.

58 GEOSCIENCES↗
DOE OSTI2013

GRB060218 AS A TIDAL DISRUPTION OF A WHITE DWARF BY AN INTERMEDIATE-MASS BLACK HOLE

The highly unusual pair of a gamma-ray burst (GRB) GRB060218 and an associated supernova, SN2006aj, has puzzled theorists for years. A supernova shock breakout and a jet from a newborn stellar mass compact object have been proposed to explain this pair's multiwavelength signature. Alternatively, we propose that the source is naturally explained by another channel: the tidal disruption of a white dwarf (WD) by an intermediate-mass black hole (IMBH). This tidal disruption is accompanied by a tidal pinching, which leads to the ignition of a WD and a supernova. Some debris falls back onto the IMBH, forms a disk, which quickly amplifies the magnetic field, and launches a jet. We successfully fit soft X-ray spectra with the Comptonized blackbody emission from a jet photosphere. The optical/UV emission is consistent with self-absorbed synchrotron emission from the expanding jet front. The temporal dependence of the accretion rate M-dot (t) in a tidal disruption provides a good fit to the soft X-ray light curve. The IMBH mass is found to be about 10{sup 4} M{sub Sun} in three independent estimates: (1) fitting the tidal disruption M-dot (t) to the soft X-ray light curve, (2) computing the jet base radius in a jet photospheric emission model, and (3) inferring the mass of the central black hole based on the host dwarf galaxy's stellar mass. The position of the supernova is consistent with the center of the host galaxy, while the low supernova ejecta mass is consistent with that of a WD. The high expected rate of tidal disruptions in dwarf galaxies is consistent with one source observed by the Swift satellite over several years at a distance of 150 Mpc measured for GRB060218. Encounters with WDs provide much fuel for the growth of IMBHs.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Bayesian Framework for Multi-Timescale State Estimation in Low-Observable Distribution Systems

To support the smart grid paradigm, there has been a significant increase in sensor deployments and metering infrastructure in distribution systems. However, the measurements provided by these sensors and metering devices are typically sampled at different rates and could suffer from losses during the aggregation process. It is crucial to effectively reconcile the time-series measurements for a reliable state estimation. While weighted least squares has been the traditional approach for state estimation, sparsity-based approaches like matrix completion have become popular due to their superior performance in low-observability conditions. This paper proposes a Bayesian framework for both multi-timescale data aggregation and matrix completion based state estimation. Specifically, the multiscale time-series data aggregated from heterogenous sources are reconciled using a multitask Gaussian process that exploits the spatio-temporal correlations. Here, the resulting consistent timeseries alongwith the confidence bound on the imputations are fed into a Bayesian matrix completion method augmented with linearized power-flow constraints to accurately estimate the states in low-observability conditions. Results on three phase unbalanced IEEE 37 and IEEE 123 bus test systems reveal the superior performance of the proposed Bayesian framework. The computational complexity for the proposed Bayesian framework is also quantified.

42 ENGINEERING↗
DOE OSTI2023

Interpretation of Fracture Initiation Points by In-Well Low-Frequency Distributed Acoustic Sensing in Horizontal Wells

Summary Low-frequency distributed acoustic sensing (LF-DAS) exploits the optical phase shift of Rayleigh backscatter in fiber-optic cables to obtain distributed measurements of changes in strain and temperature. Fiber-optic cables are often installed for multistage hydraulic fracture diagnostics in horizontal wells. LF-DAS in an untreated well provides far-field strain measurements, while offset wells are hydraulically fractured. Such a configuration is called crosswell LF-DAS sensing. Crosswell LF-DAS measurements have proved useful in diagnosing fracture hits, fracture azimuth, planarity, cluster efficiency, fracture propagation rates, and the dynamic distance to the fracture front. In contrast, in-well LF-DAS is conducted on the actively fractured well. Due to cool fracture fluid being injected at high injection rates, the strain component of the LF-DAS response is largely obscured by temperature changes. In permanent fiber-optic cable installations, distributed temperature sensing (DTS) is often conducted simultaneously with LF-DAS. An opportunity exists to decouple the temperature and strain components of LF-DAS sensors to observe strain changes on in-well LF-DAS. The LF-DAS response is modeled as linearly dependent on strain and temperature changes. Theoretical LF-DAS temperature and strain sensitivity coefficients are derived based on the changes to the index of refraction and length of the fiber. Using the DTS measurements, temperature changes are computed, smoothed, filtered, and compared to the LF-DAS response. Crossplots of the in-well LF-DAS measurements and temperature changes from DTS measurements far from the actively fractured region are used to validate the theoretical sensitivity coefficients. Uncertainty in the temperature component of the LF-DAS response is quantified. The difficulty in corresponding the different spatial and temporal resolutions of the DTS and LF-DAS measurements is overcome by comparing the responses over a moving temporal and spatial window. If the LF-DAS response at the center of the window agrees with the DTS response within uncertainty, the measurement is filtered out. After filtering, the remaining nonzero in-well LF-DAS measurements are due to changes in strain. The data are then visualized in waterfall plots. The results indicate that the theoretical and observed strain and temperature coefficients agree within 10%. After the temperature component of the in-well LF-DAS response is extracted, the remaining nonzero measurements are located primarily within the actively treated region. Locations with peaks in the strain response are interpreted to indicate fracture initiation points. These fracture initiation points are compared with in-well DTS and high-frequency DAS noise measurements across multiple stages to better understand fracture initiation along the horizontal well.

Engineering↗
DOE OSTI2020

Scaling of pulsed nanosecond capillary plasmas at different specific energy deposition

Nano-second, capillary discharges (nCDs) are unique plasma sources in their ability to sustain high specific energy deposition ω dep approaching 10 eV/molecule in molecular gases. This high energy loading on short timescales produces both high plasma densities and high densities of molecular exited states. These high densities of electrons and excited states interact with each other during the early afterglow through electron collision quenching and associative ionization. In this paper we discuss results from a two-dimensional computational investigation of a nCD sustained in air at a pressure of 28.5 mbar and with a voltage amplitude 20 kV. Discharges were investigated for two circuit configurations—a floating low voltage electrode and with the low voltage electrode connected to ground through a ballast resistor. The first configuration produced a single ionization wave from the high to low voltage electrode. The second produced converging ionization waves beginning at both electrodes. With a decrease of the tube radius, the velocity of the ionization fronts decreased while the shape of the ionization wave changed from the electron density being distributed smoothly in the radial direction, to being hollow shaped where there is a higher electron density near the tube wall. For sufficiently small tubes, the near-wall maxima merge to have the higher density on the axis of the capillary tube. In the early afterglow, the temporal and radial behavior of the N 2 (C 3 Π u ) density is a sensitive function of ω dep due to electron collision quenching. Finally, these trends indicate that starting from ω dep $\geqslant$ 0.3 eV/molecule, it is necessary to take into account interactions of electrons with electronically excited species during the discharge and early afterglow.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

High-order shock-fitted detonation propagation in high explosives

In this paper, a highly accurate numerical shock and material interface fitting scheme composed of fifth-order spatial and third- or fifth-order temporal discretizations is applied to the two-dimensional reactive Euler equations in both slab and axisymmetric geometries. High rates of convergence are not typically possible with shock-capturing methods as the Taylor series analysis breaks down in the vicinity of discontinuities. Furthermore, for typical high explosive (HE) simulations, the effects of material interfaces at the charge boundary can also cause significant computational errors. Fitting a computational boundary to both the shock front and material interface (i.e. streamline) alleviates the computational errors associated with captured shocks and thus opens up the possibility of high rates of convergence for multi-dimensional shock and detonation flows. Several verification tests, including a Sedov blast wave, a Zel'dovich–von Neumann–Döring (ZND) detonation wave, and Taylor–Maccoll supersonic flow over a cone, are utilized to demonstrate high rates of convergence to nontrivial shock and reaction flows. Comparisons to previously published shock-capturing multi-dimensional detonations in a polytropic fluid with a constant adiabatic exponent (PF-CAE) are made, demonstrating significantly lower computational error for the present shock and material interface fitting method. For an error on the order of 10 m / s , which is similar to that observed in experiments, shock-fitting offers a computational savings on the order of 1000. In addition, the behavior of the detonation phase speed is examined for several slab widths to evaluate the detonation performance of PBX 9501 while utilizing the Wescott–Stewart–Davis (WSD) model, which is commonly used in HE modeling. It is found that the thickness effect curve resulting from this equation of state and reaction model using published values is dramatically more steep than observed in recent experiments. Utilizing the present fitting strategy, in conjunction with a nonlinear optimizer, a new set of reaction rate parameters improves the correlation of the model to experimental results. In conclusion, this new model is tested against two dimensional slabs as a validation test.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2016

Analysis of the accelerated crucible rotation technique applied to the gradient freeze growth of cadmium zinc telluride

Here we employ finite-element modeling to assess the effects of the accelerated crucible rotation technique (ACRT) on cadmium zinc telluride (CZT) crystals grown from a gradient freeze system. Via consideration of tellurium segregation and transport, we show, for the first time, that steady growth from a tellurium-rich melt produces persistent undercooling in front of the growth interface, likely leading to morphological instability. The application of ACRT rearranges melt flows and tellurium transport but, in contrast to conventional wisdom, does not altogether eliminate undercooling of the melt. Rather, a much more complicated picture arises, where spatio-temporal realignment of undercooled melt may act to locally suppress instability. A better understanding of these mechanisms and quantification of their overall effects will allow for future growth optimization.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2015

Numerical modeling of the SNS H{sup −} ion source

Ion source rf antennas that produce H- ions can fail when plasma heating causes ablation of the insulating coating due to small structural defects such as cracks. Reducing antenna failures that reduce the operating capabilities of the Spallation Neutron Source (SNS) accelerator is one of the top priorities of the SNS H- Source Program at ORNL. Numerical modeling of ion sources can provide techniques for optimizing design in order to reduce antenna failures. There are a number of difficulties in developing accurate models of rf inductive plasmas. First, a large range of spatial and temporal scales must be resolved in order to accurately capture the physics of plasma motion, including the Debye length, rf frequencies on the order of tens of MHz, simulation time scales of many hundreds of rf periods, large device sizes on tens of cm, and ion motions that are thousands of times slower than electrons. This results in large simulation domains with many computational cells for solving plasma and electromagnetic equations, short time steps, and long-duration simulations. In order to reduce the computational requirements, one can develop implicit models for both fields and particle motions (e.g. divergence-preserving ADI methods), various electrostatic models, or magnetohydrodynamic models. We have performed simulations using all three of these methods and have found that fluid models have the greatest potential for giving accurate solutions while still being fast enough to perform long timescale simulations in a reasonable amount of time. We have implemented a number of fluid models with electromagnetics using the simulation tool USim and applied them to modeling the SNS H- ion source. We found that a reduced, single-fluid MHD model with an imposed magnetic field due to the rf antenna current and the confining multi-cusp field generated increased bulk plasma velocities of > 200 m/s in the region of the antenna where ablation is often observed in the SNS source. We report here on comparisons of simulated plasma parameters and code performance using more accurate physical models, such as two-temperature extended MHD models, for both a related benchmark system describing a inductively coupled plasma reactor, and for the SNS ion source. We also present results from scaling studies for mesh generation and solvers in the USim simulation code.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2017

Basic Research Needs for Future Nuclear Energy: Report of the Basic Energy Sciences Workshop for Future Nuclear Energy, August 9-11, 2017

Nuclear power today accounts for approximately 20% of the electrical supply for the US. It provides reliable energy that helps to stabilize the national grid. In contrast to current nuclear power facilities, future concepts operate at higher temperatures, are operationally more efficient, more fully utilize the energy stored within the fuel, and reduce the amount of waste produced. Recent advances in chemistry and materials sciences provide an unprecedented opportunity to develop the critical systems—the fuel, coolant, and structural materials—needed for deployment of advanced nuclear reactor designs. To address the challenge of stringent demands on material systems as their properties and performance continuously evolve over the operational lifetime, the frameworks to uncover the underlying processes that cause reduced performance have been developed over the last decade. There is now the prospect of designing materials that report upon their own damage, are self-healing, and leverage what would be considered deleterious chemical processes to increase their performance within advanced reactor designs. Success will provide an energy source needed to secure a stable and safe energy supply for the nation. To identify the key research directions required to discover the chemical processes and material degradation mechanisms that result in loss of performance of the fuel, coolant, and structural materials, nearly 150 theoreticians, computational experts, and experimentalists from industry, national laboratories, academia, and federal agencies participated in a workshop held in August 2017. Based on the discussions held at the workshop, five Priority Research Directions (PRDs) were identified. These PRDs defined the key scientific challenges that must be met as well as the opportunities that can be exploited to achieve a multi-scale spatial and temporal understanding of fundamental processes that govern the properties and performance of the different material systems required for advanced reactors.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

NANT Site - Microwave Radiometer Thermodynamic Retrievals TROPoe / Derived Data

This dataset contains daily files with thermodynamic profiles retrieved with the optimal estimation physical retrieval TROPoe (Turner and Löhnert 2014; Turner and Blumberg 2019; Turner and Löhnert 2021). The profiles are retrieved every 10 minutes from instantaneous brightness temperature measurements at 35 channels observed with a microwave radiometer MP3000A operated by NOAA Physical Sciences Laboratory on Nantucket Island for WFIP3. Additional input data in TROPoe are cloud base height from a collocated ceilometer operated by NOAA GML and temperature, water vapor mixing ratio, and pressure from a sensor attached to the MWR housing. In addition to these temporally resolved input data, TROPoe requires an a priori dataset (prior) that provides mean climatological estimates of thermodynamic profiles and specifies how temperature and humidity covary with height as an input (for details see, e.g., Djalalova et al. 2022). The prior is a key component of the retrieval and provides a constraint on the ill-posed inversion problem. A monthly prior was computed from operational radiosonde launches at Upton, NY.

17 WIND ENERGY↗
DOE OSTI2025

NANT Site - ASSIST / Thermodynamic retrievals TROPoe / Derived Data

This dataset contains daily files with thermodynamic profiles retrieved with the optimal estimation physical retrieval TROPoe (Turner and Löhnert 2014; Turner and Blumberg 2019; Turner and Löhnert 2021). The profiles are retrieved every 10 minutes from instantaneous radiances observed with an Atmospheric Sounder Spectrometer by Infrared Spectral Technology (ASSIST, Rochette et al. 2009) operated by NOAA Physical Sciences Laboratory on Nantucket Island for WFIP3. The spectral bands used in the retrieval are in the wavenumber range from 612 - 905.4 cm-1 and are specified in Turner and Löhnert (2021). Additional input data in TROPoe are cloud base height from a collocated ceilometer operated by NOAA GML and temperature, water vapor mixing ratio, and pressure from a sensor attached to the ASSIST housing. In addition to these temporally resolved input data, TROPoe requires an a priori dataset (prior) that provides mean climatological estimates of thermodynamic profiles and specifies how temperature and humidity covary with height as an input (for details see, e.g., Djalalova et al. 2022). The prior is a key component of the retrieval and provides a constraint on the ill-posed inversion problem. A monthly prior was computed from operational radiosonde launches at Upton, NY.

17 WIND ENERGY↗
DOE OSTI2025

NOAA PSL Microwave Radiometer Block Island / Thermodynamic retrievals TROPoe

This dataset contains daily files with thermodynamic profiles retrieved with the optimal estimation physical retrieval TROPoe (Turner and Löhnert 2014; Turner and Blumberg 2019; Turner and Löhnert 2021). The profiles are retrieved every 10 minutes from instantaneous brightness temperature measurements at 35 channels observed with a microwave radiometer MP3000A operated by NOAA Physical Sciences Laboratory on the Block Island for WFIP3. Additional input data in TROPoe are cloud base height from a collocated ceilometer operated by NOAA GML and temperature, water vapor mixing ratio, and pressure from a sensor attached to the MWR housing. In addition to these temporally resolved input data, TROPoe requires an a priori dataset (prior) that provides mean climatological estimates of thermodynamic profiles and specifies how temperature and humidity covary with height as an input (for details see, e.g., Djalalova et al. 2022). The prior is a key component of the retrieval and provides a constraint on the ill-posed inversion problem. A monthly prior was computed from operational radiosonde launches at Upton, NY.

17 WIND ENERGY↗
DOE OSTI2025

BLOC Site - ASSIST Thermodynamic retrievals TROPoe / Derived Data

This dataset contains daily files with thermodynamic profiles retrieved with the optimal estimation physical retrieval TROPoe (Turner and Löhnert 2014; Turner and Blumberg 2019; Turner and Löhnert 2021). The profiles are retrieved every 10 minutes from instantaneous radiances observed with an Atmospheric Sounder Spectrometer by Infrared Spectral Technology (ASSIST, Rochette et al. 2009) operated by NOAA Physical Sciences Laboratory on Block Island for WFIP3. The spectral bands used in the retrieval are in the wavenumber range from 612 - 905.4 cm-1 and are specified in Turner and Löhnert (2021). Radiances are noise-filtered but not averaged in time to minimize errors due to non-uniform clouds. Additional input data in TROPoe are cloud base height from a collocated ceilometer operated by NOAA GML and temperature, water vapor mixing ratio, and pressure from a sensor attached to the ASSIST housing. In addition to these temporally resolved input data, TROPoe requires an a priori dataset (prior) that provides mean climatological estimates of thermodynamic profiles and specifies how temperature and humidity covary with height as an input (for details see, e.g., Djalalova et al. 2022). The prior is a key component of the retrieval and provides a constraint on the ill-posed inversion problem. A monthly prior was computed from operational radiosonde launches at Upton, NY.

17 WIND ENERGY↗