Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “TIME DEPENDENCE”

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,171 records · Page 65
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Iodate interactions with calcite: implications for natural attenuation

Solid-phase interactions and speciation are important to radioiodine transport in groundwater. At the Hanford Site in Southeastern Washington State, iodate (IO 3 - ) is the main aqueous species in dilute radioiodine groundwater plumes. Like other oxyanions, IO 3 - may be incorporated into and/or adsorbed onto calcite, a common mineral at Hanford, decreasing its mobility in the environment. A series of macroscale batch experiments combined with solid phase characterization were conducted to identify variables impacting time-dependent aqueous IO 3 - removal via calcite precipitation and determine the location of IO 3 - within the calcite crystal structure. Results demonstrated 11.5-97% aqueous IO 3 - removal during initial rapid calcite precipitation. Incorporation was apparently the main removal mechanism, although later slower precipitation and/or adsorption may have also contributed to IO 3 - removal. Using a higher concentration of the calcite forming solutions (i.e., using 1M vs. 0.1M concentrations) resulted in an increase in the amount of precipitated calcite and a greater percentage of IO 3 - removed; however, calcite formed with lower molarity solutions resulted in higher IO 3 - mass (µg/g) removal. Solubility testing of laboratory produced calcites showed only small differences in solubility for calcite with and without IO 3 - incorporated into the structure. Evidence collected from SEM/FIB and TEM/SAED suggested that the IO 3 - incorporated into calcite was present in regions close to surface (implying potential easy release upon calcite dissolution).

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

Machine learning a molecular Hamiltonian for predicting electron dynamics

We develop a computational method to learn a molecular Hamiltonian matrix from matrix-valued time series of the electron density. As we demonstrate for three small molecules, the resulting Hamiltonians can be used for electron density evolution, producing highly accurate results even when propagating 1,000 time steps beyond the training data. As a more rigorous test, we use the learned Hamiltonians to simulate electron dynamics in the presence of an applied electric field, extrapolating to a problem that is beyond the field-free training data. We find that the resulting electron dynamics predicted by our learned Hamiltonian are in close quantitative agreement with the ground truth. Our method relies on combining a reduced-dimensional, linear statistical model of the Hamiltonian with a time-discretization of the quantum Liouville equation within time-dependent Hartree Fock theory. In conclusion, we train the model using a least-squares solver, avoiding numerous, CPU-intensive optimization steps. For both field-free and field-on problems, we quantify training and propagation errors, highlighting areas for future development.

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

Optical Response of Soda-Lime Glass Shocked to 14 GPa

To investigate the potential use of soda-lime glass (SLG) as an optical window in shock wave experiments, laser interferometry (VISAR) was used to determine the window correction and the refractive index in SLG samples subjected to plane shock wave compression. Symmetric plate impact experiments were conducted at peak stresses varying between 4 and 14 GPa and laser interferometry ($λ_0$ = 532$\textit {nm}$ wavelength) measurements were used to monitor the particle velocity histories at the impact surface. The higher stress experiments displayed a constant state prior to stress unloading. In contrast, the lower stress experiments (4.3 and 7.8 GPa) displayed an initial constant state followed by weakly time-dependent response prior to stress unloading. Because the actual particle velocities were independently known from the projectile velocity measurements, the window correction factor and the refractive index changes were determined over the stress range examined. Above 7.5 GPa, our window correction and refractive index results show significant differences with previously published data for the same material. Potential reasons for these differences are indicated.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Autoencoder-Based Anomaly Detection System for Online Data Quality Monitoring of the CMS Electromagnetic Calorimeter

The CMS detector is a general-purpose apparatus that detects high-energy collisions produced at the LHC. Online data quality monitoring of the CMS electromagnetic calorimeter is a vital operational tool that allows detector experts to quickly identify, localize, and diagnose a broad range of detector issues that could affect the quality of physics data. A real-time autoencoder-based anomaly detection system using semi-supervised machine learning is presented enabling the detection of anomalies in the CMS electromagnetic calorimeter data. A novel method is introduced which maximizes the anomaly detection performance by exploiting the time-dependent evolution of anomalies as well as spatial variations in the detector response. The autoencoder-based system is able to efficiently detect anomalies, while maintaining a very low false discovery rate. The performance of the system is validated with anomalies found in 2018 and 2022 LHC collision data. In addition, the first results from deploying the autoencoder-based system in the CMS online data quality monitoring workflow during the beginning of Run 3 of the LHC are presented, showing its ability to detect issues missed by the existing system.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

Leveraging operator learning to accelerate convergence of the preconditioned conjugate gradient method

We propose a new deflation strategy to accelerate the convergence of the preconditioned conjugate gradient (PCG) method for solving parametric large-scale linear systems of equations. Unlike traditional deflation techniques that rely on eigenvector approximations or recycled Krylov subspaces, we generate the deflation subspaces using operator learning, specifically the Deep Operator Network (DeepONet). To this aim, we introduce two complementary approaches for assembling the deflation operators. The first approach approximates near-null space vectors of the discrete PDE operator using the basis functions learned by the DeepONet. The second approach directly leverages solutions predicted by the DeepONet. To further enhance convergence, we also propose several strategies for prescribing the sparsity pattern of the deflation operator. Here, a comprehensive set of numerical experiments encompassing steady-state, time-dependent, scalar, and vector-valued problems posed on both structured and unstructured geometries is presented and demonstrates the effectiveness of the proposed DeepONet-based deflated PCG method, as well as its generalization across a wide range of model parameters and problem resolutions.

Deflation↗
DOE OSTI2025

2032 and 2036 risk enhancement from NEOs in the Taurid stream: Is there a significant coherent component to impact risk?

Impact risk is normally quantified by summing the product of the probability of an event and some measure of its consequences over the set of all possible events. The probability factor is considered to be more objective and is based on the size frequency distribution of Near Earth objects (NEOs) and an implicit assumption of randomness, which can be described as “stochastic risk”. Impact frequency does change with time, however, and there have been episodes in the deep geological past when the flux has been much higher. The hypothesis of “coherent catastrophism” suggests large variations on shorter timescales. It postulates the existence of a “Taurid resonant swarm” (TRS) of debris associated with Comet Encke that is stabilized by Jupiter and in a 7:2 resonance with it. The hypothetical cluster orbits in the broad Taurid stream, which crosses Earth's orbit twice a year at its nodes. Whereas the most extreme and fanciful versions of coherent catastrophism in recent geologic history have been comprehensively refuted, the possibility of a significant component of coherent (time-dependent) risk, associated with non-random correlated events, remains. This paper incorporates recently published data from observational campaigns associated with the Taurid stream. There is no evidence for objects in the Taurid stream that are above the global catastrophe threshold, but the possibility of a few large objects and a significant population of objects the size of the Chelyabinsk or Tunguska bodies has not been eliminated. Eyewitness accounts and comparison of airburst models to ground truth suggest that the Tunguska object was a Beta Taurid. The Tunguska event is probably an outlier when compared to the size frequency distributions under stochastic assumptions, but if there is a significant coherent component it may be representative of events that take place far more frequently, and risk assessments may have underestimated the contribution from airbursts. If so, the Earth will experience increased risk from objects this size, peaking in November 2032 and June 2036. We recommend targeted survey campaigns during these hypothetical TRS node crossings to quantify the population, search for potentially hazardous objects, and identify imminent impactors.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Facile healing of cracks in organic–inorganic halide perovskite thin films

Organic-inorganic halide perovskite (OIHP) thin films at the heart of the burgeoning thin-film perovskite solar cells (PSCs) technology possess poor mechanical properties, which is likely to limit the long-term reliability of PSCs as they are poised for commercialization. In an effort to address this issue, in this work we demonstrate that through-thickness cracks induced by bending-tension in thin films of the two prototypical OlHPs, methylammonium lead triiodide (MAPbl(3)) and formamidinium lead triiodide (FAPbl(3)), can be healed easily. This is through the application of either a moderate compressive stress (bending-compression) at room temperature or a simple heat-treatment at modest temperatures. The crack-healing process is found to be time-dependent, which indicates that facile mass-transport in OlHPs plays a key role in this phenomenon. An explanation for this phenomenon is provided, one based on the fundamentals of brittle fracture. This discovery has broad implications for the prevention and/or restoration of the overall performance, environmental stability, and mechanical reliability of PSCs, and other devices.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Effect of chemical disorder on the electronic stopping of solid solution alloys

The electronic stopping power of nickel-based equiatomic solid solutions alloys NiCr, NiFe and NiCo for protons and alpha projectiles is investigated in detail using real-time time-dependent density functional theory over a wide range of velocities. Recently developed numerical electronic structure methods are used to probe fundamental aspects of electron-ion coupling non-perturbatively and in a fully atomistic context, capturing the effect of the atomic scale disorder. The effects of particular electronic band structures and density of states reflect in the low velocity limit behavior. We compare our results for the alloys with those of a pure nickel target to understand how alloying affects the electronic stopping. We discover that NiCo and NiFe have similar stopping behavior as Ni while NiCr has an asymptotic stopping power that is more than a factor of two larger than its counterparts for velocities below 0.1 a.u.. Overall, we show that the low-velocity limit of electronic stopping power can be manipulated by controlling the broadening of the d-band through the chemical disorder. In this regime, the Bragg’s additive rule for the stopping of composite materials also fails for NiCr.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

High throughput automated characterization of enamel microstructure using synchrotron tomography and optical flow imaging

The remarkable damage-tolerance of enamel has been attributed to its hierarchical microstructure and the organized bands of decussated rods. A thorough characterization of the microscale rod evolution within the enamel is needed to elucidate this complex structure. While prior efforts in this area have made use of single particle tracking to track a single rod evolution to various degrees of success, such a process can be both computationally and labor intensive, limited to the evolution path of a single rod, and is therefore prone to error from potentially tracking outliers. Particle image velocimetry (PIV) is a well-established algorithm to derive field information from image sequences for processes that are time-dependent, such as fluid flows and structural deformation. In this work, we demonstrate the use of PIV in extracting the full-field microstructural distribution of rods within the enamel. Further, enamel samples from a wild African lion were analyzed using high-energy synchrotron X-ray micro-tomography. Results from the PIV analysis provide sufficient full-field information to reconstruct the growth of individual rods that can potentially enable rapid analysis of complex microstructures from high resolution synchrotron datasets. Such information can serve as a template for designing damage-tolerant bioinspired structures for advanced manufacturing.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Prediction of residual stresses in additively manufactured parts using lumped capacitance and classical lamination theory

Several industries are interested in Laser Powder Bed Fusion (L-PBF) Additively Manufactured (AM) metal parts because their designs can be made arbitrarily complex while retaining bulk-type material properties. However, the residual stresses (RS) and distortions caused by the heat gradients inherent to L-PBF processes are detrimental to the structural integrity of the parts and must be taken into consideration during the part design cycle. Predicting the state of stresses in as-built 3D printed parts is a difficult problem that is typically approached with the use of transient thermomechanical Finite Element Models (FEMs). However, the nonlinearities associated with AM processes are difficult to capture in these FEMs without increasing the computational cost of the simulation, limiting their ability to be incorporated into practical design cycles. This work presents a novel analytical framework that combines lumped capacitance nonlinear heat transfer with time dependent classical lamination theory to efficiently and accurately predict RS in as-built L-PBF parts without the need of FEMs. The simulation was compared to Neutron Diffraction (ND) residual strain measurements taken at Oak Ridge National Laboratories (ORNL) as well as Synchrotron X-ray Diffraction (XRD) strain data published by the National Institute of Standards and Technology (NIST). The simulation predictions and the experimental data showed excellent agreement for the in-plane strain directions, and general agreement for the out of plane strain component, highlighting an area where further development can be implemented.

42 ENGINEERING↗
DOE OSTI2021

Improved MPACT energy deposition and explicit heat generation coupling with CTF

The CASL neutronics code MPACT assumes local energy deposition with equilibrium delayed energy for steady-state and transient calculations. These approximations limit the accurate representation of the heat generation in space and its variations over time, both of which are essential for power distribution and thermal–hydraulic coupling analyses. This paper presents an explicit energy deposition model considering the neutron and gamma heating of all regions and the time-dependent delayed energy. The MPACT/CTF interfaces are updated to transfer the explicit heat sources. The new model agrees well with MCNP for problems without thermal–hydraulic feedback. The MPACT/CTF coupled calculation for a hot-full-power quarter-core case using the new model exhibits a reduction of peak pin power by 2.3% and fuel centerline temperature by 17 K. The new model also shows a keff difference up to -100 pcm in depletion and a 30 K peak fuel temperature reduction in a reactivity-initiated-accident problem.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2023

High-fidelity simulations of the run-in process for a pebble-bed reactor

Pebble-bed reactors (PBRs) rely on a continual feed of fuel pebbles being cycled through the core. As a result, they require a “run-in” period in order to reach an equilibrium state. The run-in period for a PBR is a complex, time-dependent problem that requires the injection of new fuel, different types of fuel, and power increases. This complexity in the run-in makes it important to capture the physical processes in order to generate an accurate representation. The present work details the creation of a high-fidelity Monte Carlo methodology for analyzing the run-in and subsequent approach to equilibrium for PBRs. The methodology entails a Python module wrapped around Serpent so as to perform neutronics calculations, move pebbles, refuel the core, and discharge pebbles, thereby modeling the explicit behavior of the PBR run-in. Further, three run-in simulations (a constant temperature profile, a linear temperature profile, and a constant temperature profile using control rods) were examined in order to identify the key physical phenomena present in the run-in process. Utilizing kugelpy, we found the inclusion of a temperature profile to be important for accurately capturing a discharge burnup (around 141 MWd/kg), a consistent k-eff (around 1.005), and an average pebble power (around 2.5 kW/pebble) that all fall within acceptable limits.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

New capabilities for rapid depletion analysis of pebble-bed reactors in SCALE

The SCALE Leap-In for Cores at Equilibrium (SLICE) method leverages capabilities available in the SCALE modeling and simulation suite to facilitate rapid estimation of equilibrium core inventories for flowing-pebble reactor systems in support of a variety of fuel cycle applications. New capabilities that compliment the SLICE method have been developed for the ORIGAMI interface to ORIGEN in SCALE to facilitate rapid depletion calculations for flowing-pebble systems. New features include a more generalized and flexible means of specifying interpolation dimensions, the ability to evaluate pebbles moving through user-defined “axial zones” in each of the pebble’s multiple passes through the core, and the treatment of differential velocities across radial channels (representing both pebble-to-pebble and wall-to-pebble friction effects). ORIGAMI thus provides an efficient user interface to define a pebble’s path and irradiation histories of its multiple passes through the core to calculate the pebble’s time-dependent inventories, which can be useful in various applications domains such as safeguards, criticality safety, and disposal analyses. The consistency of new ORIGAMI approach was verified against the solution obtained using the SLICE method with ORIGEN-ARP, showing excellent agreement.

Skutnik, Steven E. [Oak Ridge National Laboratory ↗
DOE OSTI2023

New static solutions of symmetric $\phi^4$ equation

In this paper, we provide new exact solutions of nonlinear Klein–Gordon ( $\phi^4$) equation in 1 + 1-dimension. For simplicity, we focus on the static equation and ignore the time-dependence. The symmetric $\phi^4$ equation has played an important role in several areas of physics. We obtain several novel non-singular solutions of the symmetric $\phi^4$ model in terms of the Jacobi elliptic functions and compare them with the well-known solutions. Finally, we categorize these solutions in terms of the potential parameters.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Insights into Methodologies and Operational Details of Resource Adequacy Assessment: A Case Study with Application to a Broader Flexibility Framework

Assessing and maintaining resource adequacy (RA) is a core pillar of power systems. However, recent changes in the physical makeup of these systems and the conditions under which these systems must operate have yielded a renewed interest in the methods, metrics, and assumptions that underpin RA assessments. In this paper, we systematically explore a wide range of RA modeling dimensions, including: the objective function and level of operational detail in the underlying model formulation; the quantity (look-ahead) and quality (accuracy) of data that is available for making operational decisions within those models; and the physical configuration of solar photovoltaics (PV) with battery storage hybrid resources. We apply a set of probabilistic RA tools and production cost modeling tools to a realistic test system based loosely on a future Electric Reliability Council of Texas power system dominated by solar PV resources. Under the assumptions of our system and models, we find that multi-stage probabilistic assessments may provide a more robust evaluation of RA by capturing a wider range of operational and system interactions, but this comes at a computational cost of 1-2 orders of magnitude longer run time depending on the specific configuration. In addition, the information on thermal generator availability impacts RA performance by an order of magnitude more than solar resource forecasts, which is driven by the comparatively larger magnitude of thermal outages than solar forecast errors within our test system. Lastly, the flexibility provided by hybrid and other resources can help reduce system load-shedding event frequencies and enable the system to be more robust to inaccurate forecast information, and alternative hybrid inverter sizes can impact RA levels by 1-2 orders of magnitude. Our results point to the importance of a broader flexibility framework to describe the interaction between (1) flexibility "supply" from both physical resource capabilities and operational constraints considered in the modeling, and (2) flexibility "demand" from forecast errors, thermal generator outages, and other sources of uncertainty, as well as their RA impacts. Results are likely sensitive to the system buildout explored; future work could consider additional system configurations and conditions.

ENERGY PLANNING, POLICY, AND ECONOMY,POWER TRANSMI↗
DOE OSTI2024

Simulation of gas mixture dynamics in a pipeline network using explicit staggered-grid discretization

Here we develop an explicit staggered finite difference discretization scheme for simulating the transport of highly heterogeneous gas mixtures through pipeline networks. This study is motivated by the proposed blending of hydrogen into natural gas pipelines to reduce end use carbon emissions while using existing pipeline systems throughout their planned lifetimes. Our computational method accommodates an arbitrary number of constituent gases with very different physical properties that may be injected into a network with significant spatiotemporal variation. In this setting, the gas flow physics are highly location- and time- dependent, so that local composition and nodal mixing must be accounted for. The resulting conservation laws are formulated in terms of pressure, partial densities and flows, and volumetric and mass fractions of the constituents. We include non-ideal equations of state that employ linear approximations of gas compressibility factors, so that the pressure dynamics propagate locally according to a variable wave speed that depends on mixture composition and density. We derive compatibility relationships for network edge boundary values that are more complex than for a homogeneous gas. The simulation method is evaluated on initial boundary value problems for a single pipe and a small network, is cross-validated with a lumped element simulation, and used to demonstrate a local monitoring and control policy for maintaining allowable concentration levels.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Steric effects of central dogma processes on the compaction and segregation of bacterial nucleoids

The bacterial cytoplasm is characterized by a distinctive membrane-less organelle, the nucleoid, which harbors the chromosomal DNA. Here, we investigate the steric effects of dynamic processes associated with transcription and translation on the structure of this organelle using coarse-grained molecular dynamics simulations that incorporate out-of-equilibrium reactions. Our model captures the scale of the entire cell and incorporates a reaction-diffusion system for ribosomes and polyribosomes, coupling their nonequilibrium kinetics to DNA excluded volume interactions. Our findings demonstrate that out-of-equilibrium reactions increase the size of the nucleoid and the number of ribosomes in this subcellular region. In addition, we show that nucleoid size is proportional to transcriptional activity. Our model reproduces the time-dependent change in nucleoid size observed in rifampicin treatment experiments, where the pool of polyribosomes is depleted. Furthermore, we find that these active processes are essential for complete sister chromosome separation and correct nucleoid positioning within the cell. Overall, our study reveals the effects of the central dogma processes on the internal organization and localization of bacterial nucleoids.

Chang, Mu-Hung [Univ. of Tennessee, Knoxville, TN ↗
DOE OSTI2021

Fast multiscale contrast independent preconditioners for linear elastic topology optimization problems

The goal of this work is to present a fast and viable approach for the numerical solution of the high-contrast state problems arising in topology optimization. The optimization process is iterative, and the gradients are obtained by an adjoint analysis, which requires the numerical solution of large high-contrast linear elastic problems with features spanning several length scales. The size of the discretized problems forces the utilization of iterative linear solvers with solution time dependent on the quality of the preconditioner. The lack of clear separation between the scales, as well as the high-contrast, imposes severe challenges on the standard preconditioning techniques. Thus, here we propose new methods for the high-contrast elasticity equation with performance independent of the high-contrast and the multi-scale structure of the elasticity problem. The solvers are based on two-levels domain decomposition techniques with a carefully constructed coarse level to deal with the high-contrast and multi-scale nature of the problem. The construction utilizes spectral equivalence between scalar diffusion and each displacement block of the elasticity problems and, in contrast to previous solutions proposed in the literature, is able to select the appropriate dimension of the coarse space automatically. The new methods inherit the advantages of domain decomposition techniques, such as easy parallelization and scalability. Finally, the presented numerical experiments demonstrate the excellent performance of the proposed methods.

97 MATHEMATICS AND COMPUTING↗