Search NASA⌕ Search

SEARCH · Search NASA

Results for “front propagation”

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

Quantification of the Crack Evolution Process by Extracting Relevant Signal Components from Wave Propagation and Diffusive Transport Front Measurements

Wave propagation and diffusive transport phenomena in a geological rock sample undergoing crack evolution process are expected to interact with the mechanical discontinuities in the medium. The measurements of the signals associated with these phenomena can be used to assess and monitor the crack-driven micromechanical alterations in the rock. Different wave/diffusion phenomena, such as sonic propagation, pressure diffusion, and acoustic emission (AE), are sensitive to different elements of the mechanical discontinuities generated during the evolution of the crack clusters from initiation to coalescence. Sonic propagation, AE, and pressure diffusion monitoring have the potential to map the crack evolution because the transmitter-receiver arrays can be designed, arranged and tuned to (1) achieve maximum recovery of the scattered waveforms and travel times, (2) capture the later arrivals and multiple reflections, and (3) illuminate large rock volume. However, the structural/topological complexities of the mechanical discontinuities, complex distribution of the stress fields, complex mechanical alterations in media, and fluid redistribution in the crack system pose serious challenges for the detection and modeling of the crack evolution process (from here on, we will use the term ‘crack evolution process’ to mean that the crack evolution occurred under shallow crustal conditions). For purposes of accurately accounting such complexities and heterogeneities in the absence of reliable physical laws, simulation methods, and signal processing techniques, my early-career research proposal will develop and apply novel data-driven machine learning methods to: (1) extract signal components relevant to the various phases of crack evolution and (2) generate a 2D visual map of the crack evolution process.

58 GEOSCIENCES↗
DOE OSTI2022

The evolution of curvature in planar, photoionization-driven heat fronts

Photoionized plasmas are common in astrophysics and cosmology, especially in space near compact objects, and there are effects from photoionization in high-energy-density plasmas due to the large radiation fields present. Photoionized plasmas are an active area of laboratory research and there are currently experiments to study photoionization-supported heat fronts. These photoionization fronts differ from the physics of diffusive radiation waves, commonly called Marshak waves, that are also an active area of research. This work uses a geometric argument to describe the expected evolution of the photoionization front curvature, in a planar geometry. It then compares this curvature to that of a Marshak wave as a method of diagnosing a heat front experiment. It is found that while the curvature of a planar Marshak wave increases in time, it decreases for a photoionization front. A comparison of radiation energy and electron heat fluxes through the container for the heat front propagating medium demonstrates that the geometric argument for the photoionization front curvature is sufficient. Furthermore, this comparison also demonstrates that wall losses are not significant in a photoionization front because the post-front region is very optically thin. A discussion of the implication this work has on material choice in the targets for an experiment follows.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Initiator‐Driven Tailoring of Reactivity and Inhibition in Frontally Polymerized Dicyclopentadiene‐Co‐Dihydrofuran Thermosets

Incorporating 2,3‐dihydrofuran (DHF) into dicyclopentadiene (DCPD) enables frontally polymerizable, deconstructible thermosets via acid hydrolysis. At 5–20 mol% DHF and low initiator concentration (~100 ppm), DHF inhibits front propagation, lowers the glass transition temperature ( T g ) by up to 20°C, and quenches the front, especially at 20% DHF. This paper demonstrates that increasing the initiator concentration to 250–1000 ppm overcomes inhibition, enhances front speeds by ~50%, and suppresses quenching at the expense of pot‐life. Above 5 mol% DHF, an emergence of a second exotherm is observed, attributed to increased rate dependence on the inhibition of the initiator by DHF. Void formation from volatile DHF is suppressed by applying 110 kPa nitrogen, though at the cost of reduced front speed. Dynamic mechanical analysis shows that higher initiator concentration produced more highly cross‐linked polymers, with a ~20% increase in storage modulus at 25°C and a 20°C increase in T g .

deconstructible thermosets↗
DOE OSTI2025

Urban Impacts on Convective Squall Lines over Chicago in the Warm Season—Part I: Observations of Multi-Scale Convective Evolution

In this study, our aim is to diagnose how two quasi-linear convective systems (QLCS) are organized so one can determine the possible role of the city of Chicago, IL, USA, in modifying convective precipitation systems. In this Part I of a two-part study, we employ large-scale analyses, radiosonde soundings, surface observations, and Doppler radar data to diagnose the precursor atmospheric circulations that organize the evolution of two mesoscale convective systems and compare those circulations to radar and precipitation. Several multi-scale processes are found that organize and modify convection over the Chicago metroplex. Two sequential quasi-linear convective systems (QLCS #1 and #2) were organized that propagated over Chicago, IL, USA, during an eight-hour period on 5–6 July 2018. The first squall line (QLCS #1) built from the southwest to the northeast while strengthening as it propagated over the city, and the second (QLCS #2) propagated southeastwards and weakened as it passed over the city in association with a polar cold front. The weak upper-level divergence associated with a diffluent flow poleward of an expansive ridge built over and strengthened a low-level trough and confluence zone, triggering QLCS #1. Convective downdrafts from QLCS #1 produced a cold pool that interacted with multiple confluent low-level jets surrounding and focused on the metroplex urban heat island, thus advecting the convection poleward over the metroplex. The heaviest precipitation occurred just south-southeast of Midway Airport, Chicago. Subsequently, a polar cold front propagated into the metroplex, which triggered QLCS #2. However, the descending air above it under the polar jet and residual cold pool from QLCS #1 rapidly dissipated the cold frontal convection. This represents a case study where very active convection built over the metroplex and was likely modified by it, as evidenced in numerical simulations to be described in Part II.

Kaplan, Michael L.↗
DOE OSTI2022

Robust Carbon Dioxide Plume Imaging Using Joint Tomographic Inversion of Seismic Onset Time and Distributed Pressure and Temperature Measurements (Final Report)

We develop and demonstrate rapid and cost-effective methodologies for spatiotemporal tracking of CO2 plumes during geologic sequestration using joint inversion of seismic data and distributed pressure and temperature measurements. Key elements of our methodology are: (a) a computationally efficient approach to pressure and temperature propagation, (b) analysis of time lapse seismic data using a novel ‘seismic onset time’ approach to detect fluid front propagation, and (c) data assimilation and uncertainty assessment via joint inversion of pressure, temperature and time lapse seismic data, and (d) validating the numerical tomographic inversion using a CO2 injection demonstration projects, specifically data collected from the from the Petra Nova Parish Holdings CCUS project in the West Ranch Field, Texas and the Chester-16 reef CO2 injection site in Northern Michigan which is part of the DOE Midwestern Carbon Sequestration Project. The research team is led by Texas A&M University and includes Battelle as a subcontractor with support from Shell, Anadarko, Chevron and JX Nippon. A carbon dioxide (CO2) water-alternating-gas (WAG) pilot was conducted to gain insights into tertiary oil recovery potential via CO2 flood in the West Ranch Field as part of the Petra Nova project, the world’s largest post-combustion CO2 capture and utilization initiative. With a fluvial formation geology and large contrasts in permeability, this is a challenging and novel application of CO2 enhanced oil recovery (EOR). We build a predictive dynamic model of the subsurface that incorporates the multiphase and compositional data acquired during the pilot operation. The calibrated model is used for the carbon dioxide plume imaging. The study began with an initialization of the pilot sector model extracted from a calibrated full-field model. The pilot model calibration follows a two-step hierarchical workflow. First, we performed a large-scale update of the permeability distribution by integrating available bottomhole pressure and multiphase production data. In the second step, local permeability field is fine-tuned using a streamline-based method to match CO2 breakthrough times at the producers. The predictive capability of the calibrated model was verified through two blind validation tests: (1) the model showed good agreement with saturation logs acquired at two observation wells; and (2) the model reproduced the CO2 recovery as a fraction of the injected CO2. The use of seismic onset times has shown great promise for integrating near-continuous seismic surveys for updating geologic models. In this study, we analyze the impact of seismic survey frequency on the onset time approach aiming to extend the application of onset time to infrequent seismic surveys. In addition, we quantitatively examine the nonlinearity of the onset time method and compare it to the commonly used amplitude inversion method. We carry out a sensitivity analysis of seismic survey frequency based on the complete seismic survey data (over 175 surveys) of steam injection in a heavy oil reservoir (Peace River Unit) in Canada. Our results show that an adequate onset time map can be obtained from the infrequent seismic surveys by interpolation between seismic surveys as long as there is no change in the dominant underlying physics between the successive surveys. The study also shows that nonlinearity of the onset time method can be -smaller than that of the amplitude inversion method by several orders of magnitude. Application to the Brugge benchmark case shows that the onset time method obtains comparable permeability update as the traditional seismic amplitude inversion method with faster computation and improved convergence characteristics. We extend the streamline-based data integration approach to incorporate distributed temperature sensor (DTS) data using the concept of thermal tracer travel time. Then, a hierarchical workflow composed of evolutionary and streamline methods is employed to jointly history match the DTS and pressure data. Finally, CO2 saturation and streamline maps are used to visualize the CO2 plume movement during the sequestration process. The hierarchical workflow is applied to a carbon sequestration project in a carbonate reef reservoir within the Northern Niagaran Pinnacle Reef Trend in Michigan, USA. The monitoring data set consists of distributed temperature sensing (DTS) data acquired at the injection well and a monitoring well, flowing bottom-hole pressure data at the injection well, and time-lapse pressure measurements at several locations along the monitoring well. The history matching results indicate that the CO2 movement is mostly restricted to the intended zones of injection which is consistent with an independent warm-back analysis of the temperature data. In addition to employing simulation models and inverse methods for CO2 plume imaging, we also initialized a data-driven technology for detecting inter-well connectivity based on production and pressure data. Our machine-learning framework is built on the statistical recurrent unit (SRU) model and interprets well-based injection/production data into inter-well connectivity without relying on a geologic model. We test it on synthetic and field-scale CO2 EOR projects utilizing the water-alternating-gas (WAG) process. The validation of the proposed data-driven inter-well connectivity assessment is performed using synthetic data from simulation models where inter-well connectivity can be easily measured using the streamline-based flux allocation. The SRU model is shown to offer excellent prediction performance on the synthetic case. Despite significant measurement noise and frequent well shut-ins imposed in the field-scale case, the SRU model offers good prediction accuracy, the overall relative error of the phase production rates at most producers ranges from 10% to 30%. It is shown that the dominant connections identified by the data-driven method and streamline method are in close agreement. Texas A&M University, the lead organization in the project, was primarily responsible for the development of tomographic approaches for CO2 plume mapping in conjunction with distributed pressure, temperature and seismic onset time data. Battelle, as a subcontractor, was primarily responsible for the development of analytical and empirical methods for analyzing transient injection rate and pressure data from point/line sources such as injection and monitoring wells. An additional area of emphasis for Battelle was the use of machine learning for such tasks as inferring reservoir connectivity information from injection-production data, and identifying variable importance for machine learning-based proxy models developed from full-physics simulations. The two organizations also collaborated on the application of the tomographic inversion methodology for a field data set.

02 PETROLEUM↗
DOE OSTI2020

Imaging perturbed shock propagation in powders

A novel experimental methodology is presented to study the deviatoric response of powders in shock regimes. The powders are confined to a cylindrical wedge volume, and a projectile-driven shock wave with a sinusoidally varying front propagates through the powder. The perturbed shock wave exhibits a damping behavior due to irreversible processes of viscosity and strength (deviatoric) of the powder with propagation through increasing powder thicknesses. The inclined surface of the wedge is polished and coated to establish a diffuse surface suitable for reflecting incident laser light into a high-speed camera imaging at 5 MHz. Images of the contrast loss upon shock wave arrival at the observation surface are post-processed for qualitative and quantitative information. New data of shock damping behavior with parameters of perturbation wavelength and initial shock strength are presented for powders of copper, tantalum, and tungsten carbide as well as their mixtures. We present the first full-field images showing additional spatial disturbances on the perturbed shock front that appear dependent on particle material and morphology.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

Direct Visualization of Current-Stimulated Oxygen Migration in YBa 2 Cu 3 O 7–δ Thin Films

The past years have witnessed major advancements in all-electrical doping control on cuprates. In the vast majority of cases, the tuning of charge carrier density has been achieved via electric field effect by means of either a ferroelectric polarization or using a dielectric or electrolyte gating. Unfortunately, these approaches are constrained to rather thin superconducting layers and require large electric fields in order to ensure sizable carrier modulations. In this work, we focus on the investigation of oxygen doping in an extended region through current-stimulated oxygen migration in YBa 2 Cu 3 O 7–δ superconducting bridges. The underlying methodology is rather simple and avoids sophisticated nanofabrication process steps and complex electronics. Here, a patterned multiterminal transport bridge configuration allows us to electrically assess the directional counterflow of oxygen atoms and vacancies. Importantly, the emerging propagating front of current-dependent doping δ is probed in situ by optical microscopy and scanning electron microscopy. The resulting imaging techniques, together with photoinduced conductivity and Raman scattering investigations, reveal an inhomogeneous oxygen vacancy distribution with a controllable propagation speed permitting us to estimate the oxygen diffusivity. These findings provide direct evidence that the microscopic mechanism at play in electrical doping of cuprates involves diffusion of oxygen atoms with the applied current. The resulting fine control of the oxygen content would permit a systematic study of complex phase diagrams and the design of electrically addressable devices.

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

Impact of Global Warming on U.S. Summertime Mesoscale Convective Systems: A Simple Lagrangian Parcel Model Perspective

Abstract Mesoscale convective systems (MCSs) are the dominant rainfall producer over the United States during the warm season, causing natural disasters and severe weather every year. Global climate models have large uncertainty in projecting precipitation changes in the future climate. Here, a simple Lagrangian parcel model is used to investigate the impact of global warming on MCS initiation and growth. The single-column parcel model projects a mean precipitation decrease over the central United States and an increase to its east, in agreement with the CMIP5 model projection. It also highlights the crucial role of current climate mean-state model bias in exaggerating the change in future mean precipitation projection by 25%. As for convective population, the model captures the decreased occurrence frequency of weak to moderate convection and increased frequency of strong convection due to the increased CAPE and CIN, in agreement with convection-permitting regional simulations. Novel parameterizations of gust-front propagation speed and subsidence strength are developed as guided by cloud-resolving simulations. The multicolumn parcel model employing those parameterizations captures readily the cold pool–induced upscale growth feature. It simulates smaller mesoscale clusters over the central United States under global warming due to gust-front slowdown and subsidence strength enhancement, which are further attributed to land aridity–induced weakening of initial accumulated precipitation and strengthening of updraft speed, respectively. That said, mesoscale clusters could become bigger under more favorable conditions in future climate, including boundary layer moistening, convection lifetime lengthening, and cold pool mechanical-lifting enhancement, which require further investigations to improve mechanistic understanding of future MCS changes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Emergence of liquid following laser melting of gold thin films

X-ray structural science is undergoing a revolution driven by the emergence of X-ray Free-electron Laser (XFEL) facilities. The structures of crystalline solids can now be studied on the picosecond time scale relevant to phonons, atomic vibrations which travel at acoustic velocities. In the work presented here, X-ray diffuse scattering is employed to characterize the time dependence of the liquid phase emerging from femtosecond laser-induced melting of polycrystalline gold thin films using an XFEL. In a previous analysis of Bragg peak profiles, we showed the supersonic disappearance of the solid phase and presented a model of pumped hot electrons carrying energy from the gold surface to scatter at internal grain boundaries. This generates melt fronts propagating relatively slowly into the crystal grains. By conversion of diffuse scattering to a partial X-ray pair distribution function, we demonstrate that it has the characteristic shape obtained by Fourier transformation of the measured F ( Q ). The diffuse signal fraction increases with a characteristic rise-time of 13 ps, roughly independent of the incident pump fluence and consequent final liquid fraction. This suggests the role of further melt-front nucleation processes beyond grain boundaries.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Equilibrium-gated pattern formation: How molecular dissociation thermodynamics drive emergent behavior in dissipative polymeric systems

Emergent patterns in biological systems arise through dissipative processes that balance reaction and transport phenomena, producing highly functional properties from self-regulating mechanisms. Synthetic fabrication, by contrast, often relies on user-controlled, multistep methods that lack the self-organizing capabilities of natural systems. Inspired by nature, we sought chemical systems that integrate strongly coupled reaction and transport phenomena, identifying frontal ring-opening metathesis polymerization (FROMP) as a method capable of creating diverse forms and functions through reactive processing. By employing discrete molecular initiators, FROMP allows precise control of key reaction steps—inhibition, initiation, and propagation. Using an integrated computational and experimental framework, we uncover how near-equilibrium inhibition dynamics, coupled with far-from-equilibrium reaction kinetics, drive pattern formation in frontally polymerized synthetic materials. We propose the concept of equilibrium-gated pattern formation, demonstrating how initiator chemistry can be tuned to achieve programmable macroscale properties. Our study reveals a surprising insight: Emergent behavior in FROMP systems arises from the inhibition-dominated regime of resin composition, expanding prior observations that such behavior is confined to a narrow compositional space near the boundary between front quenching and uniform front propagation. We identify a broader compositional window, far from the quenching regime, where emergent behavior reliably manifests. This expanded design space significantly enhances the operational flexibility of reactive systems and their capacity for self-organization. Furthermore, these insights provide a roadmap for designing bioinspired materials with self-organizing capabilities, unlocking possibilities in synthetic manufacturing.

Architected Polymers↗
DOE OSTI2025

Multi-device analysis of energy loss duration and pellet penetration with implications for shattered pellet injection in ITER

A robust disruption mitigation system (DMS) requires accurate characterization of key disruption timescales, one of the most notable being the thermal quench (TQ). Recent modeling of shattered pellet injection (SPI) into ITER plasmas, using JOREK and INDEX, suggests long TQ durations (6–10 ms) and slow cold front propagation due to the large plasma size. If validated, these predictions would have an impact on the desired pellet parameters and mitigation strategies for the ITER DMS. To resolve these questions, a database of SPI experiments from several small-to-large sized devices (J-TEXT, KSTAR, AUG, DIII-D, and JET) has been compiled under the auspices of the International Tokamak Physics Activity MHD, disruptions, and control topical group. Analysis of the energy loss duration (proxy for the TQ duration) with machine size is presented for both mixed neon/deuterium (Ne/D) SPI and pure deuterium (D) SPI. Several metrics for the energy loss onset (e.g. soft x-ray signal drop, I p dip, and radiation flash) were considered as the conventional metric, electron cyclotron emission, is often cut-off during SPI. Several scalings with different onset metrics showed an increase in energy loss duration with machine size. The energy loss duration was additionally shown to be a function of the ratio between the number of SPI neon atoms injected and the stored energy. Analysis of the pellet shard position relative to the cold front found that in larger devices, pellets are typically found inboard of the q = 2 surface at the energy loss onset. Lastly, the delay between the pellet shards hitting the q = 2 surface and the energy loss onset was additionally found to increase with machine size. This suggests that the pellet shards in large devices will penetrate faster and further than the cooling front.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Spatially Distributed Ramp Reversal Memory in VO 2

Ramp‐reversal memory has recently been discovered in several insulator‐to‐metal transition materials where a non‐volatile resistance change can be set by repeatedly driving the material partway through the transition. This study uses optical microscopy to track the location and internal structure of accumulated memory as a thin film of VO 2 is temperature cycled through multiple training subloops. These measurements reveal that the gain of insulator phase fraction between consecutive subloops occurs primarily through front propagation at the insulator‐metal boundaries. By analyzing transition temperature maps, it is found, surprisingly, that the memory is also stored deep inside both insulating and metallic clusters throughout the entire sample, making the metal‐insulator coexistence landscape more rugged. This non‐volatile memory is reset after heating the sample to higher temperatures, as expected. Diffusion of point defects is proposed to account for the observed memory writing and subsequent erasing over the entire sample surface. By spatially mapping the location and character of non‐volatile memory encoding in VO 2 , this study results enable the targeting of specific local regions in the film where the full insulator‐to‐metal resistivity change can be harnessed in order to maximize the working range of memory elements for conventional and neuromorphic computing applications.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Effects of sub‐mm cylindrical voids on detonation performance in PBX 9501

Abstract Internal features of varying scale and geometry are always present in explosives systems. Below a critical length, dependent on the explosive, these features can operate as a driving force for energy concentration and reaction, known as hotspots. At larger length scales, internal features result in jetting and detonation wave shaping, allowing for bulk work to be done by the explosive as is seen in shape charges. To date, a large volume of work has been performed to simulate hotspot ignition and large‐scale wave shaping. However, little experimental data exists on the effects of features in intermediate length scales, 0.1 to 1 millimeter. It has been observed in many tests that these small‐scale features can influence the high explosive (HE) performance and in some cases cause substantial damage to adjacent systems. This work provides quantitative data on the effects of machined voids moderately above more typical hotspot lengths, 0.3–0.8 mm in diameter, in PBX 9501 pressed to 1.785 g cm −3 ±2.5 mg cm −3 . Streak imaging was used to visualize void collapse, jet velocity, re‐initiation time, and wave shape evolution. Delay in detonation front propagation time was found to be linearly dependent on the void diameter and jet velocity was found to be independent within the tested range and resolution. Cut‐back experiments were used to investigate wave shape distortion and evolution downstream of the void, showing consistent growth and decay shapes across all void sizes. Simulations using CTH, a hydrocode by Sandia National Laboratory, were used to investigate void collapse showing agreement with the trends of experimental results but yielded inaccurate wave shape development delay values. Both simulation and experimental results identified several re‐initiation mechanisms with jetting and subsequent double shocking of localized HE being the dominant mechanisms for the void sizes that were studied.

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

Computer simulation of reaction and transport of core-scale serpentinization of Fe-bearing olivine and evolution of geological hydrogen

We present reactive flow and transport simulations of core-scale serpentinization of forsteritic, Fe-bearing olivine under hydrothermal conditions. The model captures fluid flow, evolving porosity–permeability, and redox-controlled H 2 production in a cylindrical core over 55 days at 245 °C and 37 bar. Key processes include olivine dissolution, precipitation of serpentine and magnetite, and oxidation of Fe2+ to Fe 3+ . Here, results show strong coupling among flow velocity, alteration front propagation, and spatially heterogeneous H 2 generation. Elevated H2 concentrations align with inlet-localized magnetite precipitation, consistent with experiments. Secondary mineral formation reduces porosity and permeability, altering transport pathways. Reactive flow at the fluid–solid interface governs H 2 generation rates and distribution. Dissolved SiO2 promotes serpentine and talc formation by suppressing brucite, while bicarbonate extends reaction duration and moderates surface complexation effects.

Geologic hydrogen↗
DOE OSTI2023

Numerical Study of Coupled Water and Vapor Flow, Heat Transfer, and Solute Transport in Variably-Saturated Deformable Soil During Freeze-Thaw Cycles

As climate change intensifies, soil water flow, heat transfer, and solute transport in the active, unfrozen zones within permafrost and seasonally frozen ground exhibit progressively more complex interactions that are difficult to elucidate with measurements alone. For example, frozen conditions impede water flow and solute transport in soil, while heat and mass transfer are significantly affected by high thermal inertia generated from water-ice phase change during the freeze-thaw cycle. To assist in understanding these subsurface processes, the current study presents a coupled two-dimensional model, which examines heat conduction-convection with water-ice phase change, soil water (liquid water and vapor) and groundwater flow, advective-dispersive solute transport with sorption, and soil deformation (frost heave and thaw settlement) in variably saturated soils subjected to freeze-thaw actions. This coupled multiphysics problem is numerically solved using the finite element method. The model's performance is first verified by comparison to a well-documented freezing test on unsaturated soil in a laboratory environment obtained from the literature. Then based on the proposed model, we quantify the impacts of freeze-thaw cycles on the distribution of temperature, water content, displacement history, and solute concentration in three distinct soil types, including sand, silt and clay textures. The influence of fluctuations in the air temperature, groundwater level, hydraulic conductivity, and solute transport parameters was also comparatively studied. The results show that (a) there is a significant bidirectional exchange between groundwater in the saturated zone and soil water in the vadose zone during freeze-thaw periods, and its magnitude increases with the combined influence of higher hydraulic conductivity and higher capillarity; (b) the rapid dewatering ahead of the freezing front causes local volume shrinkage within the non-frozen region when the freezing front propagates downward during the freezing stage and this volume shrinkage reduces the impact of frost heave due to ice formation. This gradually recovers when the thawed water replenishes the water loss zone during the thawing stage; and (c) the profiles of soil moisture, temperature, displacement, and solute concentration during freeze-thaw cycles are sensitive to the changes in amplitude and freeze-thaw period of the sinusoidal varying air temperature near the ground surface, hydraulic conductivity of soil texture, and the initial groundwater levels. Our modeling framework and simulation results highlight the need to account for coupled thermal-hydraulic-mechanical-chemical behaviors to better understand soil water and groundwater dynamics during freeze-thaw cycles and further help explain the observed changes in water cycles and landscape evolution in cold regions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Emergence of Diverse Failure Patterns in Weathering‐Induced Landslides: Insights From Particle Finite Element Simulations

Weathering is a fundamental driver of landslide evolution over geological timescales. Despite its ubiquity and importance, quantifying how weathering drives the progressive destabilization of rock slopes remains challenging. In this work, we develop a unified computational framework based on the particle finite element method to investigate the evolution of weathering‐induced landslides, from long‐term weathering to short‐term slope failure and runout dynamics. The framework integrates key processes, including weathering front propagation, time‐dependent strength degradation, rupture surface development, and post‐failure runout dynamics. Through numerical simulation experiments, we elucidate how interactions among weathering characteristics (type, intensity, and rate law), bedrock strength, fracture distribution, and slope geometry govern the failure modes and kinematics of weathering‐induced landslides. Simulations show that matrix‐dominated weathering leads to shallow translational failures, whereas fracture‐dominated weathering produces deep‐seated rotational and compound landslides. Pre‐existing fractures and slope morphology also strongly influence the movement of destabilized landmasses, affecting the failure pattern (e.g., kinematic mode and rupture surface geometry) and post‐failure behavior (e.g., runout velocity). We further demonstrate that the failure time and volume of weathered slopes are governed by the competition between gravitational driving forces and cohesive resisting forces during progressive destabilization. These findings provide new insights into the fundamental mechanisms that drive the emergence of diverse failure patterns of weathering‐induced landslides with important implications for landslide hazard assessment.

Wang, Liang [Eidgenoessische Technische Hochschule↗
DOE OSTI2023

Asynchronous domain dynamics and equilibration in layered oxide battery cathode

To improve lithium-ion battery technology, it is essential to probe and comprehend the microscopic dynamic processes that occur in a real-world composite electrode under operating conditions. The primary and secondary particles are the structural building blocks of battery cathode electrodes. Their dynamic inconsistency has profound but not well-understood impacts. In this research, we combine operando coherent multi-crystal diffraction and optical microscopy to examine the chemical dynamics in local domains of layered oxide cathode. Our results not only pinpoint the asynchronicity of the lithium (de)intercalation at the sub-particle level, but also reveal sophisticated diffusion kinetics and reaction patterns, involving various localized processes, e.g., chemical onset, reaction front propagation, domains equilibration, particle deformation and motion. These observations shed new lights onto the activation and degradation mechanisms of state-of-the-art battery cathode materials.

25 ENERGY STORAGE↗