SEARCH · Engineering Papers

Results for “phase velocity”

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 451 records · Page 25

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Integration of the Biot–Gassmann Fluid Substitution Method and Machine Learning-Based Velocity–Stress Relationship for Estimating In Situ Stresses

Recent advancements have shown that in situ stresses can be reliably estimated through an integrated machine/deep learning (ML/DL)-based framework, which relies on models trained and validated using true triaxial ultrasonic velocity (TUV) experimental data that involve measurements of ultrasonic velocity in saturated rocks under varying stress configurations. However, when the goal is to interpret lower frequency measurements, it may be more appropriate to run experiments on dry rocks and then obtain Biot–Gassmann-derived equivalent saturated velocities (low-frequency approximation) and employ these quantities for training ML/DL models to predict in situ stress. Whether the dispersion effect of frequency on the velocity–stress relationship substantially impacts in situ stress prediction is an important and unresolved question. This work presents an enhancement of ML/DL-based workflow by training and implementing ML/DL models using equivalent saturated acoustic velocities (low-frequency) obtained by applying Biot–Gassmann fluid substitution on the ultrasonic velocities of dry cores. The models were trained on TUV data sets derived from three subsurface cores extracted from the geothermal well 16B(78)-32 at the Utah FORGE site. Each core was subjected to 75 unique stress configurations for velocity measurement in the dry state. The ML/DL trained on the TUV data set with equivalent saturated velocities demonstrated promising performance to predict in situ stress in subsurface geological rocks using velocity–stress relationships with R 2 of 0.86, 0.971, and 0.975 and root mean squared error (RMSE) of 2.59, 1.92, and 1.80 for validation/testing phases of vertical, minimum horizontal, and maximum horizontal stress models, respectively. Additionally, interpretation and explanation by Shapley additive explanations (SHAP) analysis further improved scientific validation and model reliability for estimating in situ stresses.

colloids↗

DOE OSTI2020

Using Ambient Noise Seismic Interferometry and Local and Teleseismic Earthquakes to Determine Crustal Thickness and Moho Structure of the Northwestern Gulf of Mexico Margin

Abstract The northwestern part of the Gulf of Mexico has undergone two episodes of continental rifting and collision and produced structural artifacts that are now buried under many kilometers of sediments, complicating investigations of the region. The deep sedimentary package precludes outcrops and points to a need for the application of seismic techniques, but low rates of seismicity in the region and sparse seismic monitoring limit the utility of traditional seismic methods. We therefore use diverse data to perform two‐dimensional seismic tomography across the dry land portion of the margin. Data are gleaned from teleseismic and regional earthquakes, postcritical SsPmp arrivals, and direct P wave energy identified with seismic interferometry that were recorded by a broadband, three‐component array and partially overlapping short‐period, vertical‐component array. The Pn and postcritical SsPmp phase help constrain the Moho discontinuity, which a previous receiver function study suggested was absent beneath the seaward portion of this transect. A high‐velocity body is observed in the crust at the same location as the Houston Magnetic Anomaly, possibly marking rocks from the Alleghenian continental assembly. The crust thins from NW to SE, indicating that extension occurred mostly to the south of the Ouachita orogeny. Our model indicates that the margin's sediment package reaches a maximum thickness of ~15 km at the coast and becomes unresolvably thin near the Llano Uplift.

58 GEOSCIENCES↗

DOE OSTI2021

Nonlinear kinetic simulation study of the ion–ion streaming instability in single- and multi-ion species plasmas

The nonlinear evolution of the ion–ion streaming instability (IISI) is studied using numerical techniques novel to this problem that afford direct insight into the evolution of the particle distributions of each species. Here, during the linear phase of the instability, we demonstrate quantitative agreement with linear kinetic theory. Subsequently, the electrostatic field generated by the IISI causes ring-like velocity distributions of ions to form that are both heated and slowed to varying degrees relative to their initial flows. Due to variation in the trapping conditions for ion species of differing charge-to-mass ratio, when flows of multiple species interact, the nonlinear evolution of each species can be starkly different: we show a case where a lighter ion species is completely stopped by a heavier ion species via the IISI alone (i.e., without collisions) and, for the first time, demonstrate how the IISI can introduce a relative flow between ion species that initially have the same flow velocities, thereby separating them.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Converting existing optical detectors into fast x-ray detectors

Here, the very short burn time and small size of burning plasmas created at advanced laser-fusion facilities will require high-spatial-resolution imaging diagnostics with fast time resolution. These instruments will need to function in an environment of extremely large neutron fluxes that will cause conventional diagnostics to fail because of radiation damage and induced background levels. One solution to this challenge is to perform an ultrafast conversion of the x-ray signals into the optical regime before the neutrons are able to reach the detector and then to relay image the signal out of the chamber and into a shielded bunker, protected from the effects of these neutrons. With this goal in mind, the OMEGA laser was used to demonstrate high-temporal-resolution x-ray imaging by using an x-ray snout to image an imploding backlighter capsule onto a semiconductor. The semiconductor was simultaneously probed with the existing velocity interferometry system for any surface reflector (VISAR) diagnostic, which uses an optical streak camera and provided a one-dimensional image of the phase in the semiconductor as a function of time. The phase induced in the semiconductor was linearly proportional to the x-ray emission from the backlighter capsule. This approach would then allow a sacrificial semiconductor to be attached at the end of an optical train with the VISAR and optical streak camera placed in a shielded bunker to operate in a high neutron environment and obtain time-dependent one-dimensional x-ray images or time-dependent x-ray spectra from a burning plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Edge-driven transition between extended quantum anomalous Hall crystal and fractional Chern insulator in rhombohedral graphene multilayers

Fractional Chern insulators (FCI) with fractionally quantized Hall conductance at fractional fillings and an extended quantum anomalous Hall (EQAH) crystal with an integer quantized Hall conductance over an extended region of doping were recently observed in pentalayer graphene. One particularly puzzling observation is the transition between the EQAH and FCI regimes, driven either by temperature or electrical current. Here, in this study, we propose a scenario to understand these transitions based on the topologically protected gapless edge modes that are present in both the FCI and EQAH phases and should be most relevant at temperature scales below the energy gap. Our consideration is based on the simple assumption that the edge velocity in FCI is smaller than that in EQAHE and thus contributes to a higher entropy. We further argue that domains with opposite fractionally quantized Hall conductance are ubiquitous in the devices due to disorder, which gives rise to a network of edge modes. The velocity of the edge modes between domains is further reduced due to edge reconstruction. The edge velocity can also be reduced by current when the occupation of the edge mode approaches the gap edge. The edge entropy therefore drives the transition from EQAH to FCI either by temperature or current at a nonzero temperature.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2020

Steam oxidation of ytterbium disilicate environmental barrier coatings with and without a silicon bond coat

The current generation of multilayer Si/Yb 2 Si 2 O 7 environmental barrier coatings (EBCs) are temperature limited by the melting point of Si, 1414°C. To investigate higher temperature EBCs, the cyclic steam oxidation of EBCs comprised of a single layer of ytterbium disilicate (YbDS) was compared to multilayered Si/YbDS EBCs, both deposited on SiC substrates using atmospheric plasma spray. In this work, after 500 1-h cycles at 1300°C in 90 vol%H2O-10 vol%air with a gas velocity of 1.5 cm/s, both multilayer Si/YbDS and single layer YbDS grew thinner silica scales than bare SiC, with the single layer YbDS forming the thinnest scale. Both coatings remained fully adherent and showed no signs of delamination. Silica scales formed on the single layer coating were significantly more homogeneous and possessed a markedly lower degree of cracking compared to the multilayered EBC. The single layer EBC also was exposed at 1425°C in steam with a gas velocity of 14 cm/s in an alumina reaction tube. The EBC reduced specimen mass loss compared to bare SiC but formed an extensive 2nd phase aluminosilicate reaction product. A similar reaction product was observed to form on some regions of the bare SiC specimen and appeared to partially inhibit silica volatilization. The 1425°C steam exposures were repeated with a SiC reaction tube and no 2nd phase reaction product was observed to form on the single layer EBC or bare SiC.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Fluid-Kinetic Coupling: Advanced Discretizations for Simulations on Emerging Heterogeneous Architectures (LDRD FY20-0643)

Plasma physics simulations are vital for a host of Sandia mission concerns, for fundamental science, and for clean energy in the form of fusion power. Sandia's most mature plasma physics simulation capabilities come in the form of particle-in-cell (PIC) models and magnetohydrodynamics (MHD) models. MHD models for a plasma work well in denser plasma regimes when there is enough material that the plasma approximates a fluid. PIC models, on the other hand, work well in lower-density regimes, in which there is not too much to simulate; error in PIC scales as the square root of the number of particles, making high-accuracy simulations expensive. Real-world applications, however, almost always involve a transition region between the high-density regimes where MHD is appropriate, and the low-density regimes for PIC. In such a transition region, a direct discretization of Vlasov is appropriate. Such discretizations come with their own computational costs, however; the phase-space mesh for Vlasov can involve up to six dimensions (seven if time is included), and to apply appropriate homogeneous boundary conditions in velocity space requires meshing a substantial padding region to ensure that the distribution remains sufficiently close to zero at the velocity boundaries. Moreover, for collisional plasmas, the right-hand side of the Vlasov equation is a collision operator, which is non-local in velocity space, and which may dominate the cost of the Vlasov solver. The present LDRD project endeavors to develop modern, foundational tools for the development of continuum-kinetic Vlasov solvers, using the discontinuous Petrov-Galerkin (DPG) methodology, for discretization of Vlasov, and machine-learning (ML) models to enable efficient evaluation of collision operators. DPG affords several key advantages. First, it has a built-in, robust error indicator, allowing us to adapt the mesh in a very natural way, enabling a coarse velocity-space mesh near the homogeneous boundaries, and a fine mesh where the solution has fine features. Second, it is an inherently high-order, high-intensity method, requiring extra local computations to determine so-called optimal test functions, which makes it particularly suited to modern hardware in which floating-point throughput is increasing at a faster rate than memory bandwidth. Finally, DPG is a residual-minimizing method, which enables high-accuracy computation: in typical cases, the method delivers something very close to the $L^2$ projection of the exact solution. Meanwhile, the ML-based collision model we adopt affords a cost structure that scales as the square root of a standard direct evaluation. Moreover, we design our model to conserve mass, momentum, and energy by construction, and our approach to training is highly flexible, in that it can incorporate not only synthetic data from direct-simulation Monte Carlo (DSMC) codes, but also experimental data. We have developed two DPG formulations for Vlasov-Poisson: a time-marching, backward-Euler discretization and a space-time discretization. We have conducted a number of numerical experiments to verify the approach in a 1D1V setting. In this report, we detail these formulations and experiments. We also summarize some new theoretical results developed as part of this project (published as papers previously): some new analysis of DPG for the convection-reaction problem (of which the Vlasov equation is an instance), a new exponential integrator for DPG, and some numerical exploration of various DPG-based time-marching approaches to the heat equation. As part of this work, we have contributed extensively to the Camellia open-source library; we also describe the new capabilities and their usage. We have also developed a well-documented methodology for single-species collision operators, which we applied to argon and demonstrated with numerical experiments. We summarize those results here, as well as describing at a high level a design extending the methodology to multi-species operators. We have released a new open-source library, MLC, under a BSD license; we include a summary of its capabilities as well.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

On a Critical Acceleration Scale of Dark Matter in ΛCDM and Dynamical Dark Energy

Abstract Universal acceleration a 0 emerges in various empirical laws, yet its fundamental nature remains unclear. Using Illustris and Virgo N -body simulations, we focus on the velocity and acceleration fluctuations in collisionless dark matter involving long-range gravity. For comparison, in the kinetic theory of gases, molecules undergo random elastic collisions involving short-range interactions, where only velocity fluctuations are relevant. Hierarchical structure formation proceeds through the merging of smaller halos to form larger halos, which facilitates a continuous energy cascade from small to large halos at a constant rate ε u ≈ −10 −7 m 2 s −3 . Velocity fluctuations involve a critical velocity u c ∝ (1 + z ) −3/4 . Acceleration fluctuations involve a critical acceleration a c ∝ (1 + z ) 3/4 . Two critical quantities are related by the rate of energy cascade ε u ≈ − a c u c /[2(3 π ) 2 ], where factor 3 π is from the angle of incidence during merging. With critical velocity u c on the order of 300 km s −1 at z = 0, the critical acceleration is determined to be a c 0 ≡ a c ( z = 0) ≈ 10 −10 m s −2 , suggesting a c might explain the universal acceleration a 0 ≈ 10 −10 m s −2 in the empirical Tully–Fisher relation or modified Newtonian dynamics. The redshift evolution a c ∝ (1 + z ) 3/4 is in good agreement with Magneticum and EAGLE simulations and in reasonable agreement with limited observations. This suggests a larger a 0 at a higher redshift such that galaxies of fixed mass rotate faster at a higher redshift. Note that for dark energy (DE) density ρ DE 0 ≈ a c 0 2 / G = 1 0 − 10 J m −3 , we postulate an entropic origin of the DE from acceleration fluctuations of dark matter, analogous to the gas pressure from velocity fluctuations. This leads to a dynamical DE coupled to the structure evolution involving a relatively constant DE density followed by a slow weakening phase, suggesting possible deviations from the standard ΛCDM paradigm.

N-body simulations↗

DOE OSTI2021

A differential emissivity imaging technique for measuring hydrometeor mass and type

Abstract. The Differential Emissivity Imaging Disdrometer (DEID) is a new evaporation-based optical and thermal instrument designed to measure the mass, size, density and type of individual hydrometeors as well as their bulk properties. Hydrometeor spatial dimensions are measured on a heated metal plate using an infrared camera by exploiting the much higher thermal emissivity of water compared with metal. As a melted hydrometeor evaporates, its mass can be directly related to the loss of heat from the hotplate assuming energy conservation across the hydrometeor. The heat loss required to evaporate a hydrometeor is found to be independent of environmental conditions including ambient wind velocity, moisture level and temperature. The difference in heat loss for snow vs. rain for a given mass offers a method for discriminating precipitation phase. The DEID measures hydrometeors at sampling frequencies of up to 1 Hz with masses and effective diameters greater than 1 µg and 200 µm, respectively, determined by the size of the hotplate and the thermal camera specifications. Measurable snow water equivalent (SWE) precipitation rates range from 0.001 to 200 mm h−1, as validated against a standard weighing bucket. Preliminary field experiment measurements of snow and rain from the winters of 2019 and 2020 provided continuous automated measurements of precipitation rate, snow density and visibility. Measured hydrometeor size distributions agree well with canonical results described in the literature.

47 OTHER INSTRUMENTATION↗

DOE OSTI2022

Validation of seismic velocity models in southern California with full-waveform simulations

SUMMARY Crustal seismic velocity models provide essential information for many applications including earthquake source properties, simulations of ground motion and related derivative products. We present a systematic workflow for assessing the accuracy of velocity models with full-waveform simulations. The framework is applied to four regional seismic velocity models for southern California: CVM-H15.11, CVM-S4.26, CVM-S4.26.M01 that includes a shallow geotechnical layer, and the model of Berg et al. For each model, we perform 3-D viscoelastic wave propagation simulations for 48 virtual seismic noise sources (down to 2 s) and 44 moderate-magnitude earthquakes (down to 2 s generally and 0.5 s for some cases) assuming a minimum shear wave velocity of 200 m s–1. The synthetic waveforms are compared with observations associated with both earthquake records and noise cross-correlation data sets. We measure, at multiple period bands for well-isolated seismic phases, traveltime delays and normalized zero-lag cross-correlation coefficients between the synthetic and observed data. The obtained measurements are summarized using the mean absolute derivation of time delay and the mean correlation coefficient. These two metrics provide reliable statistical representations of model quality with consistent results in all data sets. In addition to assessing the overall (average) performance of different models in the entire study area, we examine spatial variations of the models’ quality. All examined models show good phase and waveform agreements for surface waves at periods longer than 5 s, and discrepancies at shorter periods reflecting small-scale heterogeneities and near-surface structures. The model performing best overall is CVM-S4.26.M01. The largest misfits for both body and surface waves are in basin structures and around large fault zones. Inaccuracies generated in these areas may affect tomography and model simulation results at other regions. The seismic velocity models for southern California can be improved by adding better resolved structural representations of the shallow crust and volumes around the main faults.

58 GEOSCIENCES↗

DOE OSTI2021

High-speed characterization of two-phase flow and bubble dynamics in titanium felt porous media for hydrogen production

Bubble dynamics and two-phase flow phenomena are closely related to the performance of proton exchange membrane electrolyzer cells (PEMECs). This paper reports an in-situ study of the oxygen bubble behavior and associated multiphase evolutions in the anode side of PEMECs with titanium (Ti) felt liquid gas diffusion layers (LGDLs) by a high-speed visualization system. The micro oxygen bubble dynamics was captured and analyzed at different locations and virous operating conditions. The results show that the bubble detachment frequency and detachment diameter greatly increase with the operating current density. Additionally, they are significantly impacted by the local pore structure and morphology of Ti felt LGDLs. In the flow channels, there exist only several discrete micro bubbles at a low current density (0.04 A/cm 2 ) and a large flow velocity (133 mm/s). At a current density (0.2 A/cm 2 ) and a flow velocity (67 mm/s), a number of gas slugs are formed in the follow channels, in addition to discrete micro bubbles. At a high current density (1 A/cm 2 ) and a flow velocity (67 mm/s), more bubbles appear in the channel, and the flow field is dominated by slug or annular flows. These investigations can help to better understand the two-phase flow and bubble detachment mechanism, and provide a foundation for electrochemical reaction, multiphase flow studies and optimize the design of gas diffusion layers and flow fields for PEMECs in the future.

08 HYDROGEN↗

DOE OSTI2020

Role of sheared E × B flow in self-organized, improved confinement states in magnetized plasmas

A great scientific success story of magnetic fusion research in the past several decades is the theoretical development and experimental testing of the process of turbulence decorrelation and stabilization by sheared E×B flow which shows that E×B shear effects are ubiquitous in magnetized plasmas. This concept of turbulence decorrelation and stabilization has the universality needed to explain the H-mode edge transport barriers seen in limiter and divertor tokamaks, stellarators, and mirror machines; the broader edge transport barrier seen in VHmode plasmas; and the core transport barriers formed in tokamaks. Similar effects are seen in linear devices. These examples of confinement improvement are of considerable physical interest; it is not often that a system self-organizes to reduce transport when an additional source of free energy is applied to it. The transport decrease associated with E×B velocity shear is also of great practical benefit to fusion research, since it contributed to substantially increased fusion yield in all DT magnetic fusion experiments conducted to date. The fundamental physics involved in transport reduction is the effect of E×B shear on the growth, radial extent, and phase correlation of turbulent eddies in the plasma. The same basic transport reduction process can be operational in various portions of the plasma because there are a number of ways to change the radial electric field Er. An important secondary theme in this area is the synergistic effect of E×B velocity shear and magnetic shear. Although the E×B velocity shear appears to have an effect on broader classes of microturbulence, magnetic shear can mitigate some potentially harmful effects of E×B velocity shear and facilitate turbulence stabilization. Our present understanding in this area is the result of a multi-decade, intertwined effort in theory, modeling and diagnostic development combined with continuing experimental investigations. These experiments have clearly demonstrated that increased E×B shear causes reductions in turbulence and transport. Finally, the experimental findings are generally consistent with the basic theoretical models although considerable work remains to be done before we have a fully predictive theory of transport in magnetized plasmas including E×B shear effects.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

Atomically Resolved Acoustic Dynamics Coupled with Magnetic Order in a Van der Waals Antiferromagnet

Magnetoelastic coupling in van der Waals (vdW) magnetic materials enables a unique interplay between the spin and lattice degrees of freedom. Characterizing the elastic responses with atomic and femtosecond resolution across the magnetic transition is essential for guiding the design of magnetically tunable actuators and strain-mediated spintronic devices. Here, ultrafast X-ray diffraction employed at a free-electron laser reveals that the atomic displacements, wave vectors, and dispersion relations of acoustic phonon modes in a vdW antiferromagnet FePS3 are coupled with the magnetic order, by tracking both in-plane and out-of-plane Bragg peaks upon optical excitation across the Néel temperature (TN). One transverse mode shows that a quasi-out-of-plane atomic displacement undergoes a significant directional change across TN. Its quasi-in-plane wave vector is derived by comparing the measured sound velocity and the first-principles calculations. The other transverse mode is an interlayer shear acoustic mode whose amplitude is strongly enhanced in the antiferromagnetic phase, exhibiting eight times stronger amplitude than the longitudinal acoustic mode below TN. The atomically resolved characterization of acoustic phonon dynamics that couple with magnetic ordering opens opportunities for harnessing unique magnetoelastic coupling in vdW magnets on ultrafast timescales.

coherent acoustic phonons↗

DOE OSTI2021

An Eulerian Vlasov-Fokker–Planck algorithm for spherical implosion simulations of inertial confinement fusion capsules

Here, we present a numerical algorithm that enables a phase-space adaptive Eulerian Vlasov–Fokker–Planck (VFP) simulation of inertial confinement fusion (ICF) capsule implosions. The approach relies on extending a recent mass, momentum, and energy conserving phase-space moving-mesh adaptivity strategy to spherical geometry. In configuration space, we employ a mesh motion partial differential equation (MMPDE) strategy while, in velocity space, the mesh is expanded/contracted and shifted with the plasma’s evolving temperature and drift velocity. The mesh motion is dealt with by transforming the underlying VFP equations into a computational (logical) coordinate, with the resulting inertial terms carefully discretized to ensure conservation. To deal with the spatial and temporally varying dynamics in a spherically imploding system, we have developed a novel nonlinear stabilization strategy for MMPDE in the configuration space. The strategy relies on a nonlinear optimization procedure that optimizes between mesh quality and the volumetric rate change of the mesh to ensure both accuracy and stability of the solution. Implosions of ICF capsules are driven by several boundary conditions: (1) an elastic moving wall boundary; (2) a time-dependent Maxwellian Dirichlet boundary; and (3) a pressure-driven Lagrangian boundary. Of these, the pressure-driven Lagrangian boundary driver is new to our knowledge. The implementation of our strategy is verified through a set of test problems, including the Guderley and Van-Dyke implosion problems — the first-ever reported using a Vlasov–Fokker–Planck model.

1D2V↗

DOE OSTI2022

Electric-field control of the nucleation and motion of isolated three-fold polar vertices

Recently various topological polar structures have been discovered in oxide thin films. Despite the increasing evidence of their switchability under electrical and/or mechanical fields, the dynamic property of isolated ones, which is usually required for applications such as data storage, is still absent. Here, we show the controlled nucleation and motion of isolated three-fold vertices under an applied electric field. At the PbTiO 3 /SrRuO 3 interface, a two-unit-cell thick SrTiO 3 layer provides electrical boundary conditions for the formation of three-fold vertices. Utilizing the SrTiO 3 layer and in situ electrical testing system, we find that isolated three-fold vertices can move in a controllable and reversible manner with a velocity up to ~629 nm s -1 . Microstructural evolution of the nucleation and propagation of isolated three-fold vertices is further revealed by phase-field simulations. This work demonstrates the ability to electrically manipulate isolated three-fold vertices, shedding light on the dynamic property of isolated topological polar structures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2021

The effect of entrance flow development on vortex formation and wall shear stress in a curved artery model

In this work, we numerically investigate the effect of entrance condition on the spatial and temporal evolution of multiple three-dimensional vortex pairs and the wall shear stress distribution in a curved artery model. We perform this study using a Newtonian blood-analog fluid subjected to a pulsatile flow with two inflow conditions. The first flow condition is fully developed while the second condition is undeveloped (i.e., uniform). We discuss the connection along the axial direction between regions of organized vorticity observed at various cross sections of the model and compare results between the different entrance conditions. We model a human artery with a simple, rigid 180° curved pipe with a circular cross section and constant curvature, neglecting the effects of taper, torsion, and elasticity. Numerical results are computed from a discontinuous high-order spectral element flow solver. The flow rate used in this study is physiological. We observe differences in secondary flow patterns, especially during the deceleration phase of the physiological waveform where multiple vortical structures of both Dean-type and Lyne-type coexist. The results indicate that decreased axial velocities under an undeveloped condition produce smaller secondary flows that ultimately inhibit growth of any interior flow vortices. We highlight the effect of the entrance condition on the formation of these structures and subsequent appearance of abnormal inner wall shear stresses, which suggest there may be a lower prevalence of cardiovascular disease in curved arteries where the flow is rather undeveloped—a potentially physiologically significant result to help understand the influence of blood flow development on disease.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2025

Pressure-temperature equation of state of Al 2 ⁢O 3 up to 14 Mbar and 40 kK

Sapphire (Al 2 ⁢O 3 ), known for its remarkable incompressibility at ambient conditions, plays a pivotal role in both static and dynamic compression research. Accurately characterizing its equation of state (EoS) is essential for these applications. Here, we present a complete Hugoniot of Al 2 ⁢ O 3 as locus of experimentally assessed, high-precision, pressure, density and temperature states up to 14 Mbar and 43 kK. The Hugoniot is established with single shock experiments using magnetically launched hyper velocity flyers on the Z Accelerator at Sandia National Laboratories. We explore principal Hugoniot states at very high shock 𝑇 and 𝑝 in the solid phase, tracking the solid-liquid boundary and culminating at 2.4-fold compression, where data provides a direct constraint on the liquid phase. Corresponding shock release data probe thermodynamic states complementary to the Hugoniot and place additional constraints on tabular EoS models. Our findings indicate a significant deviation from existing tabular EoS models for Al 2 ⁢ O 3 dictating a comprehensive overhaul. We develop two advanced EoSs for Al 2 ⁢ O 3 the SESAME 97412 model, featuring an extensive phase diagram that includes three solid phases and the liquid phase, and the updated LEOS 2200m2 model. EoS development is assisted with Quantum Molecular Dynamics simulations. Our experimental data allows for stringent testing of our EoSs. Both models accurately capture the Hugoniot of Al 2 ⁢O 3 up to the highest pressures and temperatures. Rigorous experimental determination of extreme pressures and temperatures, paired with sophisticated models, advances the frontier of EoS development beyond 1 terapascal.

Kalita, Patricia [Sandia National Laboratories (SN↗

DOE OSTI2025

Charges on a suspended silicon nitride membrane under a high-energy electron beam

Thin silicon nitride (SiN x ) membranes are widely used in gas and liquid phase transmission electron microscopy (TEM) and as phase plates to enhance imaging contrast. SiN x contains trap sites for both positive and negative charges, which can be manipulated by high-energy electron irradiation, external potential biasing, or light exposure. Charge accumulation on the membrane can significantly affect in situ TEM processes, including chemical and electrochemical reactions, nanoparticle dynamics, and catalytic activity, or introduce unwanted phase shifts when used as a phase plate. Here, in this study, charge accumulation on suspended SiN x membranes was investigated using off-axis electron holography combined with model-free charge analysis, supported by custom finite element analysis (FEA) simulations. An average residual positive charge density of approximately 2.8 × 10 −4 C m −2 was measured. Localized and stable regions of both positive and negative charges were identified on the membrane. The global positive and localized positive/negative charges give rise to strong electric fields and electroosmotic slip velocities at the membrane surface, which are sufficient to induce non-Brownian particle behavior and directional fluid flow, offering a physical explanation for previously observed anomalies in particle dynamics, nucleation, and growth during gas and liquid phase TEM experiments. These results provide a benchmark for understanding charge behavior at SiNx interfaces in gas and liquid phase TEM. Furthermore, the FEA simulations establish a framework for future investigations into charge distribution, electrostatic potentials, and electrical double layers at solid–liquid interfaces, particularly in complex geometries and chemically dynamic environments.

Suspended silicon↗