SEARCH · Engineering Papers

Results for “Timescale”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 397 records · Page 22

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Thermal Gradient Effects on Local Hotspot Ignition in 1,3,5,7‐Tetranitro‐1,3,5,7‐tetrazocane (HMX)

Understanding hotspot ignition, growth, and criticality, as well as the timescales of each, is crucial for parameterizing mesoscale and continuum‐level models that rely on a statistical understanding of hotspots. However, these models often consider hotspots to have uniform temperatures or for hotspots of a given temperature and size to always behave the same. Therefore, using molecular dynamics simulations, we assess the influence of thermal distribution effects on hotspot local ignition and time to ignition in 1,3,5,7‐tetranitro‐1,3,5,7‐tetrazocane (HMX) nanoscale hotspots. Finally, by assessing hotspots with a gradient driven from an initial two‐temperature core–shell setup, we show that small increases in the core temperature under a constant average temperature can lead to order‐of‐magnitude effects on reaction and local ignition timescales.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Imaging Bias‐Driven Domain Wall Motion With Scanning Oscillator Piezoresponse Force Microscopy

Understanding ferroelectric domain wall dynamics at the nanoscale across a broad range of timescales requires measuring domain wall position under different applied electric fields. The success of piezoresponse force microscopy (PFM) as a tool to apply local electric fields at different positions and imaging their changing position, together with the information obtained from associated switching spectroscopies has fueled numerous studies of the dynamics of ferroelectric domains to determine the impact of intrinsic parameters such as crystalline order, defects and pinning centers, as well as boundary conditions such as environment. However, the investigation of sub-coercive reversible domain wall vibrational modes requires the development of new tools that enable visualizing domain wall motion under varying applied fields with high temporal and spatial resolution while also accounting for spurious electrostatic effects. Here, scanning oscillator piezoresponse force microscopy extends the investigation of domain wall dynamics to new regimes, providing direct visualization of domain wall position as a function of an external electric field that varies in time and location. In conclusion, this enables studying the energetics of field-driven ferroelectric domain wall motion, which is shown to obey a thermally activated flow regime in the millisecond timescale.

Piezoresponse Force Microscopy↗

DOE OSTI2020

Identification of major moisture sources across the Mediterranean Basin

We employ a Lagrangian based moisture back trajectory method on an ensemble of four reanalysis datasets to provide a comprehensive understanding of moisture sources over the Mediterranean land region (30° N–49.5° N and 9.75° W–61.5° E) at seasonal timescales for 1980–2013 period. Using a source region between 10° S–71.35° N along the latitude and 80° W–84.88° E along the longitude that is subdivided into ten complimentary sub-regions, our analyses is able to backtrack up to > 90% of seasonal precipitation at each grid point within the target region. Our results indicate a significant role of moisture advected from the North Atlantic and Mediterranean Sea, and locally recycled moisture over the target region in shaping the spatial organization of seasonal precipitation. However, a clear east–west contrast is witnessed in determining the relative importance of each of these major moisture sources where the North Atlantic dictates the moisture supply over the western Mediterranean while moisture from Mediterranean Sea and local recycling play a key role over the eastern Mediterranean. Our analyses also demonstrate a major footprint of the North Atlantic Oscillation (NAO) on precipitation variability over the Mediterranean land as dynamic and thermodynamic anomalies during the negative phase of NAO match with those during wet years and vice versa. The findings reported here are generally consistent across the four reanalysis datasets. Overall, this study establishes the relative roles of adjacent and far-off oceanic and terrestrial evaporative sources over the Mediterranean land and should help in understanding the drivers of precipitation variability and change at varying timescales.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

A regional (land–ocean) comparison of the seasonal to decadal variability of the Northern Hemisphere jet stream 1871–2011

Seasonal to decadal variations in Northern Hemisphere jet stream latitude and speed over land (Eurasia, North America) and oceanic (North Atlantic, North Pacific) regions are presented for the period 1871–2011 from the Twentieth Century Reanalysis dataset. Significant regional differences are seen on seasonal to decadal timescales. Seasonally, the jet latitude range is lower over the oceans compared to land, reduced from 20° over Eurasia to 10° over the North Atlantic where the ocean meridional heat transport is greatest. The mean jet latitude range is at a minimum in winter (DJF), particularly along the western boundary of the North Pacific and North Atlantic, where the land-sea contrast and SST gradients are strongest. The 141-year trends in jet latitude and speed show differences on a regional basis. The North Atlantic has significant increasing jet latitude trends in all seasons, up to 3° in winter. Eurasia has significant increasing trends in winter and summer, however, no increase is seen across the North Pacific or North America. Jet speed shows significant increases evident in winter (up to 4.7 ms –1 ), spring and autumn over the North Atlantic, Eurasia and North America however, over the North Pacific no increase is observed. Long term trends are generally overlaid by multidecadal variability, particularly evident in the North Pacific, where 20-year variability in jet latitude and jet speed are seen, associated with the Pacific Decadal Oscillation which explains 50% of the winter variance in jet latitude since 1940. The results highlight that northern hemisphere jet variability and trends differ on a regional basis (North Atlantic, North Pacific, Eurasia and North America) on seasonal to decadal timescales, suggesting that different mechanisms are influencing the jet latitude and speed. This is important from a climate modelling perspective and for climate predictions in the near and longer term.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Observation of Fundamental Mechanisms in Compression-Induced Phase Transformations Using Ultrafast X-ray Diffraction

As theoretically hypothesized for several decades in group IV transition metals, we have discovered a dynamically stabilized body-centered cubic (bcc) intermediate state in Zr under uniaxial loading at sub-nanosecond timescales. Under ultrafast shock wave compression, rather than the transformation from α-Zr to the more disordered hex-3 equilibrium ω-Zr phase, in its place we find the formation of a previously unobserved nonequilibrium bcc metastable intermediate. We probe the compression-induced phase transition pathway in zirconium using time-resolved sub-picosecond x-ray diffraction analysis at the Linac Coherent Light Source. We also present molecular dynamics simulations using a potential derived from first-principles methods which independently predict this intermediate phase under ultrafast shock conditions. In contrast with experiments on longer timescale (> 10 ns) where the phase diagram alone is an adequate predictor of the crystalline structure of a material, our recent study highlights the importance of metastability and time dependence in the kinetics of phase transformations.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Convolutional neural networks for intra-hour solar forecasting based on sky image sequences

Accurate and timely solar forecasts play an increasingly critical role in power systems. Compared to longer forecasting timescales, very short-term solar forecasting has lagged behind in both research and practice. In this paper, we propose deep convolutional neural networks (CNNs) to provide operational intra-hour (10-minute-ahead to 60-minute-ahead) solar forecasts. We develop two CNN structures inspired by a widely-used CNN architecture. The CNNs are tailored to our solar forecasting regression tasks and rely solely on sky image sequences. Case studies based on six years of data (over 150,000 data points) demonstrate that the best CNN model has forecast skill scores of 20%-39% over the naive persistence of cloudiness benchmark, even at these very short timescales. The CNNs also have consistently superior performance when compared to shallow machine learning models with meteorological predictors, where the improvement averages around 7%. The sensitivity analyses show that the sky image length, resolution, and weather conditions have impacts on the deep learning model accuracy. In our intra-hour problem with specific setups, two sky images with a 10-minute 128 x 128 resolution yield the most accurate forecasts. Current limitations, future work, and deployment challenges and solutions are also discussed.

14 SOLAR ENERGY↗

DOE OSTI2022

Multiscale simulation of integrated energy system and electricity market interactions

Accelerating the deep decarbonization of the world's electric grids requires the coordination of complex energy systems and infrastructures across timescales from seconds to decades. Here, we present a new multiscale simulation framework that integrates process- and grid-centric modeling paradigms to better design, operate, and control integrated energy systems (IESs), which combine multiple technologies, in wholesale energy markets. Traditionally, IESs are analyzed with a process-centric paradigm such as levelized cost of electricity (LCOE) or annualized net revenue, ignoring important interactions with electricity markets. This framework explicitly models the complex interactions between an IES's bidding, scheduling, and control decisions and the energy market's clearing and settlement processes, while incorporating operational uncertainties. Through two case studies, we show the importance of understanding and quantifying complex resource-grid interactions. In case study 1, we demonstrate that optimized bidding from one resource shifts the profit distribution for all energy systems in the market. This result suggests new and more flexible IES technologies can disrupt the economics of all market participants, possibly leading to accelerated retirements of less flexible resources. Interestingly, the optimized bidding has little impact on grid-level aggregate statistics, such as total generation costs and renewable penetration rate. While aggregate modeling strategies may remain valid under some IES adoption scenarios for analysis focused on regional outcomes, direct comparisons of IES technologies at specific locations without considering these interactions may lead to misleading or incorrect conclusions. In case study 2, we consider the design and flexible operation of IESs that hybridize conventional generators with energy storage. Through a sensitivity analysis, we find that as the size of the storage system increases, the total number of start-ups for coal- and natural gas-based IESs reduced by 25% and 33.6%, and the total thermal generator ramping (i.e., mileage) reduced by 86.5% and 62.5%, respectively. This shows the primary benefit of storage may not be reduced operational costs (which do not change significantly) but fewer start-ups and less ramping, which may greatly simplify the design, operation, and control of carbon capture systems. The new modeling and optimization capabilities from this work enable the coupling of rigorous, dynamic process models with grid-level production cost models to quantitatively identify the nuanced interdependencies across these vast timescales that must be addressed to realize clean, safe, and secure energy production. Moreover, the proposed general multiscale simulation framework is applicable to all IES technologies and can be easily extended to consider other energy carriers (e.g., hydrogen, ammonia) and energy infrastructures (e.g., natural gas pipelines).

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2020

Modeling of Hidden Structures Using Sparse Chemical Shift Data from NMR Relaxation Dispersion

NMR relaxation dispersion measurements report on conformational changes occurring on the μs-ms timescale. Chemical shift information derived from relaxation dispersion can be used to generate structural models of weakly populated alternative conformational states. Current methods to obtain such models rely on determining the signs of chemical shift changes between the conformational states, which are difficult to obtain in many situations. Here, we use a “sample and select” method to generate relevant structural models of alternative conformations of the C-terminal-associated region of Escherichia coli dihydrofolate reductase (DHFR), using only unsigned chemical shift changes for backbone amides and carbonyls ( 1 H, 15 N, and 13 C'). We find that CS-Rosetta sampling with unsigned chemical shift changes generates a diversity of structures that are sufficient to characterize a minor conformational state of the C-terminal region of DHFR. The excited state differs from the ground state by a change in secondary structure, consistent with previous predictions from chemical shift hypersurfaces and validated by the x-ray structure of a partially humanized mutant of E. coli DHFR (N23PP/G51PEKN). The results demonstrate that the combination of fragment modeling with sparse chemical shift data can determine the structure of an alternative conformation of DHFR sampled on the μs-ms timescale. Such methods will be useful for characterizing alternative states, which can potentially be used for in silico drug screening, as well as contributing to understanding the role of minor states in biology and molecular evolution.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Coarse-grained molecular dynamics simulation of solvent-dependent cellulose nanofiber interactions

Associations between cellulose are important both in biofuel production and in the use of cellulose for biomaterials. Cellulose nanofibers (CNFs) are sustainable, strong, light-weight alternatives to traditional materials in manufacturing, but are challenging to obtain due to irreversible aggregation in solution during preparative fibrillation. Therefore, it is imperative to understand the underlying factors driving aggregation with a view to designing solvents that can effectively compete with interfiber interactions, hence reducing aggregation. Molecular dynamics (MD) simulation at atomic detail can provide useful information on local interactions. However, the length and timescales accessible are too short to fully capture association processes. Here, we provide a method for accessing the longer length and timescales required using coarse-grained (CG) MD simulations with a MARTINI force field to calculate the interaction behavior of CNFs in three selected solvents: NaOH-urea-water, acetone, and neat water. The CG results are consistent with our prior all-atom MD and with previous experimental results. While acetone is found not to be an effective solvent, urea and ionic moieties in NaOH-urea-water not only solvate the fibrils but also improve the confinement of water molecules around them as shown by the solvent residence times and mean-square displacements. Overall, the presence of urea and ions reduces the likelihood of aggregation in multi-CNF systems relative to neat water irrespective of whether the hydrophobic or hydrophilic CNF surfaces are interacting. In conclusion, the CG method shows clear promise for selecting potential high-performance solvents for experimental prioritization in bioenergy and biomaterials research in a relatively fast manner as well as for understanding the aggregation and rheological behavior of CNF-solvent systems.

aggregation↗

DOE OSTI2025

Simulation and generalized Langevin equation study of lipid subdiffusion in biomembrane phases

In biomembranes lipid mobility exhibits deviations from the classical diffusive behavior of Brownian particles, i.e. ‘anomalous’ diffusion. The question arises as to how this anomalous diffusive behavior varies in gel, ripple, and fluid biomembrane phases. Here, in this study, we perform all-atom molecular dynamics (MD) simulations of dimyristoylphosphatidylcholine (DMPC) bilayers in the three different phases and analyze the results using the framework of the Generalized Langevin Equation (GLE). This analysis emphasizes sub-diffusive behavior on the relatively short, picosecond-nanosecond timescales, capturing local molecular constraints and transient caging effects during the crossover of atomic dynamics from vibrational to incipient anharmonic motion. The ripple and gel phases are found to exhibit strong transient caging and prolonged memory effects resulting in distinct sub-diffusive behavior. The role of hydrogen bonding in lipid confinement is also examined, demonstrating its influence on phase-dependent molecular ordering and on short-time diffusional constraints. These findings demonstrate the GLE framework’s utility in characterizing molecular transport and lipid dynamics, with implications for longer timescale membrane dynamics.

Malik, Sheeba [Oak Ridge National Laboratory (ORNL↗

DOE OSTI2023

Mineral detection of neutrinos and dark matter. A whitepaper

Minerals are solid state nuclear track detectors — nuclear recoils in a mineral leave latent damage to the crystal structure. Depending on the mineral and its temperature, the damage features are retained in the material from minutes (in low-melting point materials such as salts at a few hundred °C) to timescales much larger than the 4.5 Gyr-age of the Solar System (in refractory materials at room temperature). The damage features from the O (50) MeV fission fragments left by spontaneous fission of 238U and other heavy unstable isotopes have long been used for fission track dating of geological samples. Laboratory studies have demonstrated the readout of defects caused by nuclear recoils with energies as small as O (1) keV. This whitepaper discusses a wide range of possible applications of minerals as detectors for E R ≳ O (1) keV nuclear recoils: Using natural minerals, one could use the damage features accumulated over O (10) Myr–O (1) Gyr to measure astrophysical neutrino fluxes (from the Sun, supernovae, or cosmic rays interacting with the atmosphere) as well as search for Dark Matter. Using signals accumulated over months to few-years timescales in laboratory-manufactured minerals, one could measure reactor neutrinos or use them as Dark Matter detectors, potentially with directional sensitivity. Research groups in Europe, Asia, and America have started developing microscopy techniques to read out the O (1)–O (100) damage features in crystals left by O (0.1)–O (100) keV nuclear recoils. We report on the status and plans of these programs. Furthermore, the research program towards the realization of such detectors is highly interdisciplinary, combining geoscience, material science, applied and fundamental physics with techniques from quantum information and Artificial Intelligence.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2025

Universal method for the optimization of HDC coating uniformity on non-planar, non-stationary substrates for inertial confinement fusion targets

The thickness uniformity of chemical vapor deposited (CVD) diamond coatings on non-planar, non-stationary substrates depends on both the intrinsic instantaneous coating thickness distribution (ICTD) of the coating conditions used and, if applicable, on the frequency of substrate reorientation. While important for many CVD diamond applications, the relative impact of the ICTD and substrate reorientation on the coating thickness uniformity has not been studied. In this work, we systematically investigate the effect of these factors for microwave-plasma chemical vapor deposition (MPCVD) of diamond (referred to as high density carbon (HDC) in the inertial confinement fusion (ICF) community) coatings on spherical, rolling substrates. This coating technique is used to fabricate capsules for ICF experiments, which require extreme coating uniformity with <0.3 % thickness variation (so-called Mode 1 or M1) to ensure symmetric compression of imploding targets. To extract the otherwise unobservable reorientation timescale (Δt), Monte Carlo simulations were performed using experimental ICTD data as input. This combined approach confirms scaling relationships between the substrate reorientation timescale as well as coating thickness and coating uniformity, as expected from a 3D random walk. Simulations confirm that M1 is Rayleigh-distributed and scales as (Δt) 1/2 , consistent with the randomization of two angles that determine orientation of a sphere. We also demonstrate that, under the conditions studied, Δt is the dominant factor in determining thickness uniformity while the intrinsic ICTD has minimal impact. Finally, experiments show that Δt can be affected by total batch size under constant agitation conditions due to space constraints that limit the capsule reorientation kinetics. In conclusion, this study highlights the utility of a combined experiment-simulation approach as a general methodology for understanding and improving coating uniformity on non-planar, non-stationary substrates.

Capsule↗

DOE OSTI2022

Controlling pollutant emissions in a high-pressure combustor with fuel-diluent blending

This paper investigates the formation of nitric oxide (NO) and carbon monoxide (CO) emissions as a function of nitrogen (N 2 ) or carbon dioxide (CO 2 ) diluent content, added to a premixed reacting jet in crossflow. Here, reaction characteristics of a rich methane-air jet injected into a lean vitiated crossflow were analyzed at an elevated pressure of 5 atm. The jet was pre-heated and enriched with 0%, 15%, 30%, and 50% mass fraction to quantify the effect on pollutant emissions. Simulated results of the full chemistry Star-CCM+CFD model were verified with data taken in the experimental high-pressure combustion facility, which provides pressure, temperature and velocity profiles, as well as line-of-sight chemiluminescence and exit emission measurements. The analysis revealed the significant influence of the diluent to delay axial combustion; this mixing delay increased the axial flame lift-off and controlled any thermal ignition limitations. Increasing the diluent content increased the timescale for the flame to stabilize, which allowed for greater entrainment of crossflow oxidizer into the axial jet stream, and led to decreased pollutant emissions. Hence, crossflow entrainment is a critical driving force at high diluent content. Local nitric oxide (NO) emission formation in the axial stage was predicted numerically, showing the correlation between diluent addition, axial heat release, and the formation of nitric oxide pollution. The high diluent levels leaned out the jet mixture and delayed axial combustion, while reducing prompt and thermal NOx by mitigating flame hotspots and minimizing the timescale that the products remain at high temperatures.

03 NATURAL GAS↗

DOE OSTI2023

Chronology of alkali anorthosite 14304 clast “b” records basin forming impact at ~3.95 Ga

Alkali-suite rocks constitute one of three major suites of lunar crustal rocks. As such, constraining their formation timescales and petrogenesis is important for understanding the earliest magmatic history of the Moon. However, the magmatic history of alkali-suite rocks is partly obscured by superimposed effects of major basin-forming impact events on the lunar nearside. Consequently, unambiguous crystallization ages of samples from this suite of rocks have not been determined. Here, the aim of this study is to better understand the petrogenetic history of the alkali-suite and the potential superimposed effects of impact metamorphism by determining Sm-Nd, Rb-Sr, and 40 Ar/ 39 Ar ages for an alkali anorthosite clast from Apollo 14 lithic breccia 14304 termed clast “b”. The new chronologic measurements of clast “b” yield concordant Sm-Nd, Rb-Sr, and 40 Ar/ 39 Ar ages of 3947 ± 13 Ma, 3975 ± 34 Ma, and 3937 ± 37 Ma respectively, resulting in a weighted mean age of 3949 ± 11 Ma. This age is not interpreted to date an igneous event related to production of the lunar highlands crust and instead the chronology and petrography of clast “b” are most readily explained by an impact event at ~3.95 Ga that caused near-complete isotopic re-equilibration of the Sm-Nd, Rb-Sr, and 40 Ar/ 39 Ar chronometers. The weighted mean age of 3949 ± 11 Ma of clast “b” is several hundred Ma younger than 4.3–4.4 Ga ages typically determined for samples of other crustal rock suites but in very good agreement with independent estimates for the formation of Imbrium basin ejecta and other marginally older impact events which are thought to have been sampled at the Apollo 14 landing site. Thus, although petrologic and geochemical examination suggest that clast “b” is a pristine igneous clast, its age likely records an impact event at the Apollo 14 landing site. Whereas the various determined ages of clast “b” do not reflect the timescales of alkali-suite magmatism, the relatively low initial Sr and Nd isotopic compositions of clast “b” indicate that its protolith evolved with very low 147 Sm/ 144 Nd and 87 Rb/ 86 Sr that are distinct from estimates for urKREEP but similar to that of lunar plagioclase. This implies that the igneous protolith of clast “b” derived from a plagioclase-dominated KREEP-rich source that must have formed after the formation of the urKREEP source at ~4.35 Ga but well before the impact event recorded by clast “b” at ~3.95 Ga.

58 GEOSCIENCES↗

DOE OSTI2025

Evolving physical and mechanical rock properties during exposure at Earth's surface

Rocks at Earth's surface experience progressive damage due to exposure to environmental conditions. Although environmental damage to rocks is recognized as being time-dependent, our quantitative understanding of its effects on physical and mechanical properties remains limited by field and laboratory studies constrained to short timescales. In this study, we aim to rigorously quantify how these properties evolve under longer natural environmental exposure spanning from the Holocene to the late Pleistocene. We conducted a ‘time series’ analysis of key physical and mechanical rock properties using multiple interior samples from 10 cogenetic granitoid boulders exposed for different amounts of time at two sites in Eastern California, USA. We find that over the late Pleistocene (~30 ka), rocks exhibit porosity increases of 40 ± 4 % and seismic velocity decreases of 7 ± 1 % at both sites. Similarly, over the same time period, rock tensile and compressive strengths declined by an average of 11 ± 4 %, and rock permeability increased on average by 104 ± 50 %. Microscopic analysis of thin sections of the same rocks indicates that microfracture density and intensity also progressively increased over the same timescales. These results suggest that time-dependent fracture propagation due to environmental exposure is likely a central contributor to the observed changes in rock properties. This study provides quantitative documentation that pronounced changes in the physical and mechanical properties of rock progressively occur within thousands of years of exposure.

Critical zone↗

DOE OSTI2026

Activation dynamics of a water-soluble human mu-opioid receptor

The mu-opioid receptor (MOR), a class A G protein-coupled receptor mediates opioid analgesia and remains a central target for pain therapeutics. While crystal structures of MOR exist, they provide limited insight into the receptor’s dynamic conformational landscape underlying function. Here, we engineered a thermostable water-soluble MOR variant (wsMOR) that retains native-like ligand-binding and activation dynamics. This variant enables high-yield production and detailed solution-phase structural studies that are challenging with membrane-embedded MOR, providing a valuable tool for studying receptor activation and aqueous-phase drug screening. Using a combined computational and experimental approach, we performed long-timescale all-atom molecular dynamics simulations together with neutron scattering and single-molecule FRET, revealing a structurally stable receptor with a diverse ensemble of conformations at different temporal resolutions. In the ligand-free state, wsMOR displayed high conformational flexibility, which decreased upon agonist binding, particularly in transmembrane helix 6, a hallmark of G protein-coupled receptor activation. Positive allosteric modulation and G protein binding further stabilized active-like states. These findings highlight wsMOR’s conformational plasticity across picosecond to millisecond timescales and provide a foundation for structure-guided development of next-generation opioid ligands with improved efficacy and safety.

E, Agyemang [University of Tennessee Knoxville]↗

DOE OSTI2024

One-sweep moment-based semi-implicit-explicit integration for gray thermal radiation transport

Thermal radiation transport (TRT) is a time dependent, high dimensional partial integro-differential equation. In practical applications such as inertial confinement fusion, TRT is coupled to other physics such as hydrodynamics, plasmas, etc., and the timescales one is interested in capturing are often much slower than the radiation timescale. As a result, TRT is treated implicitly, and due to its stiffness and high dimensionality, is often a dominant computational cost in multiphysics simulations. Here we develop a new approach for implicit-explicit (IMEX) integration of gray TRT in the deterministic SN setting, which requires only one sweep per stage, with the simplest first-order method requiring only one sweep per time step. The partitioning of equations is done via a moment-based high-order low-order formulation of TRT, where the streaming operator and first two moments are used to capture the asymptotic stiff regimes of the streaming limit and diffusion limit. Absorption-reemission is treated explicitly, and although stiff, is sufficiently damped by the implicit solve that we achieve stable accurate time integration without incorporating the coupling of the high order and low order equations implicitly. Due to nonlinear coupling of the high-order and low-order equations through temperature-dependent opacities, to facilitate IMEX partitioning and higher-order methods, we use a semi-implicit integration approach amenable to nonlinear partitions. In conclusion, results are demonstrated on thick Marshak and crooked pipe benchmark problems, demonstrating orders of magnitude improvement in accuracy and wallclock compared with the standard first-order implicit integration typically used.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

Comparison of continuum and cross-core theories of dynamic strain aging

Dynamic strain aging (DSA) is the process of solute atoms segregating around dislocations on the timescale of loading. Continuum theories of DSA derived from elasticity theory have been shown to severely overpredict both the timescale and strengthening of DSA. Recently, cross-core theory was developed to reconcile this gap, invoking a special single-atomic-hop diffusion mechanism across the core of an extended dislocation. In this work, we show that the classical continuum theory expression for the rate of solute segregation is in error. After correcting this error, we show that continuum theory predictions match cross-core theory when the elevated diffusivity near the dislocation core is accounted for. Finally, our findings indicate that continuum theory is still a useful tool for studying dislocation-solute interactions.

36 MATERIALS SCIENCE↗