Search NASA⌕ Search

SEARCH · Search NASA

Results for “Critical 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.

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

At least 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Thermodynamics and collisionality in firehose-susceptible high- β plasmas

We study the evolution of collisionless plasmas that, due to their macroscopic evolution, are susceptible to the firehose instability, using both analytic theory and hybrid-kinetic particle-in-cell simulations. We establish that, depending on the relative magnitude of the plasma β, the characteristic time scale of macroscopic evolution and the ion-Larmor frequency, the saturation of the firehose instability in high-β plasmas can result in three qualitatively distinct thermodynamic (and electromagnetic) states. By contrast with the previously identified ‘ultra-high-beta’ and ‘Alfvén-inhibiting’ states, the newly identified ‘Alfvén-enabling’ state, which is realised when the macroscopic evolution time τ exceeds the ion-Larmor frequency by a β-dependent critical parameter, can support linear Alfvén waves and Alfvénic turbulence because the magnetic tension associated with the plasma’s macroscopic magnetic field is never completely negated by anisotropic pressure forces. We characterise these states in detail, including their saturated magnetic-energy spectra. The effective collision operator associated with the firehose fluctuations is also described; we find it to be well approximated in the Alfvén-enabling state by a simple quasi-linear pitch-angle scattering operator. The box-averaged collision frequency is ν eff ∼ β/τ, in agreement with previous results, but certain subpopulations of particles scatter at a much larger (or smaller) rate depending on their velocity in the direction parallel to the magnetic field. Our findings are essential for understanding low-collisionality astrophysical plasmas including the solar wind, the intracluster medium of galaxy clusters and black hole accretion flows. We show that all three of these plasmas are in the Alfvén-enabling regime of firehose saturation and discuss the implications of this result.

astrophysical plasmas↗
DOE OSTI2022

Role of the third dimension in searching for Majorana fermions in α-RuCl 3 via phonons

Understanding phonons in α-RuCl 3 is critical to analyze the controversy around the observation of the half-integer thermal quantum Hall effect. While many studies have focused on the magnetic excitations in α-RuCl 3 , its vibrational excitation spectrum has remained relatively unexplored. We investigate the phonon structure of α-RuCl 3 via inelastic neutron-scattering experiments and density-functional-theory calculations. Our results show excellent agreement between experiment and first-principles calculations. After validating our theoretical model, we extrapolate the low-energy phonon properties. We find that the phonons in α-RuCl 3 that either propagate or vibrate in the out-of-plane direction have significantly reduced velocities and therefore have the potential to dominate the observability of the elusive half-integer plateaus in the thermal Hall conductance. In addition, we use low-energy interlayer phonons to resolve the low-temperature stacking structure of our large crystal of α-RuCl 3 , which we find to be consistent with that of the R$\bar{3}$ space group, in agreement with neutron diffraction.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

High Temperature Erosion Modeling in Particle Based CSP Systems

1. Introduction. Wear and erosion damage of materials from solid particle and surface interactions is a major issue in various industries. Although more common in mining and oil and gas production, erosion is becoming a critical issue in renewable technologies as well such as particle based concentrated solar thermal power (CSP) systems. In particle based CSP systems, solid particles are used to absorb solar energy and as thermal storage. However, these particles may cause significant amount of wear to system components while through the system at high temperatures. This damage can be costly, and therefore, requires a greater understanding of solid particle erosion in CSP systems. Throughout the years, models and tools have been developed to predict and control erosion in industries such as oil and gas production. However, these erosion models and erosion prediction tools have been mainly developed based on erosion data for much higher velocities and lower temperatures, i.e. the operating conditions for which the existing models have been built are not pertinent to those expected in CSP system. It is known that erosion depends on many factors such as material properties, erodent particle properties, and particle impact speed and angle. However, the effect of temperature on erosion is not vastly investigated. In this work, a temperature-based correlation is introduced that will be used to modify the existing erosion models to predict erosion rates at velocities and temperatures relevant to Gen3 CSP systems. The existing models are first validated against erosion experiments run at low temperatures and low velocities. Subsequently, a temperature correction term is developed that can extend the existing models to high temperatures, based on the available experimental data showing the effect of temperature on impact erosion. 2. Erosion Models. Continuing sequence of impacts from solid particles on surfaces would result in loss of material due to mechanical interaction between solid surface and particles. Erosion equations have been developed to predict erosion under different conditions and erosion mechanism, including cutting and deformation erosion. These models are mainly either mechanistic, empirical, and semi-mechanistic models, when the latter combine the theory of the erosion mechanism and particle motion with the available experimental data. One of the first empirical correlation was introduced by American Petroleum Institute (API) Recommended Practice (RP) 14E [1]. This correlation was very conservative, and several improved empirical and semi-mechanistic equations were developed in the years following. More recently, the models introduced by Erosion/Corrosion Research Center (E/CRC) at the University of Tulsa and Arabnejad. et al. [2] are frequently used in the literature and industry, as they account for many parameters affecting erosion including particle impact speed and angle, material density and hardness, and particle size and shape. However, both these models were developed based on data at relatively high velocities and also do not take into account any mechanical changes in material associated with thermal cycling at high temperatures. 3. Erosion Prediction at High Temperatures. In this work, the accuracy of Arabnejad et al. model and E/CRC model is investigated to predict impact erosion at conditions relevant to Gen3 CSP systems. Two significant deviations expected in Gen3 systems compared to operating conditions used to build these models are particle impact velocities and system temperature. As a first step, the performance of these models was validated at low velocities and temperatures. Impact erosion experiments were conducted on SS316 coupons using HSP 40/70 ceramic particles. Assuming particle rate of 1 kg/s/m, 0.0254 m of particle curtain thickness, and particle velocity of 1.5 m/s in the system, an overall erosion of 0.59 mm/year was calculated. The experimental results were subsequently compared to those from computational simulations and erosion of 0.47 mm/year was obtained using the Arabnejad et al. model. The results indicated that the models work well at low velocities and low temperature conditions. To account for temperature effect, a mathematical correlation was developed using data published by DUCOM [3] for Inconel 600 eroded by alumina particles at high velocities. The correlation was applied to both Arabnejad et. al model and the E/CRC model. The prediction results from these modified models were within 20% of this experimental data. Predictions of of erosion by the modified E/CRC model at three different temperatures are made. Similar to the calculation at low temperature, a typical CSP system with 1 kg/s/m of particle flow rate per unit length of the particle curtain and a curtain thickness of 0.0254 m is assumed. The annual thickness loss calculations were run for different particle impact velocities, assuming a uniform particle impact area equal to the cross-sectional area of the curtain (particle-particle interactions and dispersion of particles are not considered in obtaining the results). It is observed that erosion increases exponentially as temperature increases. Furthermore, it is also noted that, erosion changes non-linearly with impact velocity. Currently, experiments are also being conducted to measure erosion of SS316 at 800 ºC with HSP 40/70 particles at low impact velocities. We expect to use the results from high temperature testing to further improve the temperature correlation function. Similar models are also being developed for abrasion erosion resulting from particle sliding along the surfaces as well as attrition from particle to particle and particle to surface interactions. References. [1] Institute, A. P. (1991). API Recommended Practice for Design and Installation of Offshore Production Platform Piping System, API RP 14E. [2] Arabnejad, H., Mansouri, A., Shirazi, S. A., and McLaury, B. S. (2015a). Development of mechanistic erosion equation for solid particles. Wear, 332–333, 1044–1050. http://doi.org/10.1016/j.wear.2015.01.031. [3] https://ducom.com/high-temperature-erosion-evaluating-sample-wear/

14 SOLAR ENERGY↗
DOE OSTI2024

Exploring ship track spreading rates with a physics-informed Langevin particle parameterization

Abstract. The rate at which aerosols spread from a point source injection, such as from a ship or other stationary pollution source, is critical for accurately representing subgrid plume spreading in a climate model. Such climate model results will guide future decisions regarding the feasibility and application of large-scale intentional marine cloud brightening (MCB). Prior modeling studies have shown that the rate at which ship plumes spread may be strongly dependent on meteorological conditions, such as precipitating versus non-precipitating boundary layers and shear. In this study, we apply a Lagrangian particle model (PM-ABL v1.0), governed by a Langevin stochastic differential equation, to create a simplified framework for predicting the rate of spreading from a ship-injected aerosol plume in sheared, precipitating, and non-precipitating boundary layers. The velocity and position of each stochastic particle is predicted with the acceleration of each particle being driven by the turbulent kinetic energy, dissipation rate, momentum variance, and mean wind. These inputs to the stochastic particle velocity equation are derived from high-fidelity large-eddy simulations (LES) equipped with a prognostic aerosol–cloud microphysics scheme (UW-SAM) to simulate an aerosol injection from a ship into a cloud-topped marine boundary layer. The resulting spreading rate from the reduced-order stochastic model is then compared to the spreading rate in the LES. The stochastic particle velocity representation is shown to reasonably reproduce spreading rates in sheared, precipitating, and non-precipitating cases using domain-averaged turbulent statistics from the LES.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Feasibility study for test rig assessments of fish passage conditions in a Kaplan turbine

The assessment of fish passage conditions in hydroelectric turbines consists of identifying and quantifying physical magnitudes leading to increased risks of injury of fish passing through turbines in operation. Such assessments are usually carried out either with the use of computer-based methods during design or with field testing of live fish and sensors passing through prototypes. A method in between consists of test rig experimentation, which is critical for testing fish-focused design concepts and offers the opportunity for implementing the most effective design measures for improved fish survivability. However, fish-related assessments in test rigs are not sufficiently documented for industrial applications. This work presents the main findings of an experimental campaign to quantify fish-related hydraulic magnitudes in a physical model of a Kaplan turbine in a commercial test rig. Two operating conditions were tested by releasing miniaturized autonomous sensor devices (termed Sensor Fish Mini) at the turbine intake flow, passing them through the runner in motion and recovering them at the draft tube exit. During passage, time series of acceleration, absolute pressure and rotational velocity were recorded. The recordings were then interpreted to determine the magnitude and likely location of hydraulic stressors hazardous to fish. The statistical tests on the reported measurements indicated that low pressure, collisions on the runner and rotations in the draft tube were not different between the two tested operating points. On the other hand, pressure drop and collision rates on the distributor differed considerably as a function of net head. The outcomes of this investigation showed that test rig evaluations of fish-related properties with Sensor Fish Mini can contribute to an evidence-based development of turbine geometries designed for providing safer passage conditions. Future work will investigate the scaling of test rig measurements to hydraulically equivalent magnitudes in the prototype and their biological consequences.

13 HYDRO ENERGY↗
DOE OSTI2021

Dynamic Phase Alignment in Inertial Alfven Turbulence

In weakly-collisional plasma environments with sufficiently low electron beta, Alfvenic turbulence transforms into inertial Alfvenic turbulence at scales below the electron skin-depth. We argue that, in inertial Alfvenic turbulence, both energy and generalized kinetic helicity exhibit direct cascades. We demonstrate that the two cascades are compatible due to the existence of a strong scale-dependence of the phase alignment angle between velocity and magnetic field fluctuations, with the phase alignment angle scaling as proportional to the perpendicular wavenumber to the power of -1. The kinetic and magnetic energy spectra scale as the perpendicular wavenumber to the power of -5/3 and -11/3, respectively. As a result of the dual direct cascade, the generalized-helicity spectrum scales as the perpendicular wavenumber to the power of -5/3, implying progressive balancing of the turbulence as the cascade proceeds to smaller scales in the sub-electron-skin-depth range. Turbulent eddies exhibit a phase-space anisotropy, with the parallel wavenumber scaling as the perpendicular wavenumber to the power of 5/3, consistent with critically-balanced inertial Alfven fluctuations. Our results may be applicable to a variety of geophysical, space, and astrophysical environments, including the Earth's magnetosheath and ionosphere, solar corona, non-relativistic pair plasmas, as well as to strongly rotating non-ionized fluids.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Time-resolved site-selective imaging of predissociation and charge transfer dynamics: the CH 3 I B-band

The predissociation dynamics of the 6s (B 2 E) Rydberg state of gas-phase CH 3 I were investigated by time-resolved Coulomb-explosion imaging using extreme ultraviolet (XUV) free-electron laser pulses. Inner-shell ionization at the iodine 4d edge was utilized to provide a site-specific probe of the ensuing dynamics. The combination of a velocity-map imaging (VMI) spectrometer coupled with the pixel imaging mass spectrometry (PImMS) camera permitted three-dimensional ionic fragment momenta to be recorded simultaneously for a wide range of iodine charge states. In accord with previous studies, initial excitation at 201.2 nm results in internal conversion and subsequent dissociation on the lower-lying A-state surface on a picosecond time scale. Examination of the time-dependent yield of low kinetic energy iodine fragments yields mechanistic insights into the predissociation and subsequent charge transfer following multiple ionization of the iodine products. The effect of charge transfer was observed through differing delay-dependencies of the various iodine charge states, from which critical internuclear distances for charge transfer could be inferred and compared to a classical over-the-barrier model. Time-dependent photofragment angular anisotropy parameters were extracted from the central slice of the Newton sphere, without Abel inversion, and highlight the effect of rotation of the parent molecule before dissociation, as observed in previous works. Our results demonstrate the ability to perform three-dimensional ion imaging at high event rates and showcase the potential benefits of this approach, particularly in relation to further time-resolved studies at free-electron laser facilities.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Turbulence and Control of Wind Farms

The dynamics of the turbulent atmospheric boundary layer play a fundamental role in wind farm energy production, governing the velocity field that enters the farm as well as the turbulent mixing that regenerates energy for extraction at downstream rows. Understanding the dynamic interactions among turbines, wind farms, and the atmospheric boundary layer can therefore be beneficial in improving the efficiency of wind farm control approaches. Anticipated increases in the sizes of new wind farms to meet renewable energy targets will increase the importance of exploiting this understanding to advance wind farm control capabilities. This review discusses approaches for modeling and estimation of the wind farm flow field that have exploited such knowledge in closed-loop control, to varying degrees. We focus on power tracking as an example application that will be of critical importance as wind farms transition into their anticipated role as major suppliers of electricity. The discussion highlights the benefits of including the dynamics of the flow field in control and points to critical shortcomings of the current approaches.

Automation & Control Systems↗
DOE OSTI2024

Analytical model for the motion and interaction of two-dimensional active nematic defects

Here, we develop an approximate, analytical model for the velocity of defects in active nematics by combining recent results for the velocity of topological defects in nematic liquid crystals with the flow field generated from individual defects in active nematics. Importantly, our model takes into account the long-range interactions between defects that result from the flows they produce as well as the orientational coupling between defects inherent in nematics. Our work complements previous studies of active nematic defect motion by introducing a linear approximation that allows us to treat defect interactions as two-body interactions and incorporates the hydrodynamic screening length as a tuning parameter. We show that the model can analytically predict bound states between two +1/2 winding number defects, effective attraction between two –1/2 defects, and the scaling of a critical unbinding length between ±1/2 defects with activity. The model also gives predictions for the trajectories of defects, such as the scattering of +1/2 defects by –1/2 defects at a critical impact parameter that depends on activity. In the presence of circular confinement, the model predicts a braiding motion for three +1/2 defects that was recently seen in experiments, as well as stable and ergodic trajectories for four or more defects.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Development and performance of high-temperature superconducting CORC ® cables with CFD REBCO tapes

Abstract Increasing the normal zone propagation velocity (NZPV) in superconducting cables based on rare earth barium copper oxide (REBCO) tapes is expected to cause a paradigm shift by enabling the development of faster and more reliable quench detection systems to mitigate the development of destructive hot spots. Furthermore, a higher NZPV is beneficial in terms of fault current limitation capabilities by accelerating the homogenization of the quench in superconducting power devices. One way envisioned to increase the NZPV of cables based on REBCO tapes is the current flow diverter (CFD) concept. A Conductor on Round Core (CORC ® ) cable made with CFD REBCO tapes, called CFD CORC ® cable, and a CORC ® cable made with regular REBCO tapes, were fabricated and tested. The critical current of the CFD and regular CORC ® cables were obtained at temperatures ranging from 67 to 77 K in self-field. Measurements showed that the NZPV was increased by a factor of 4–4.5 in the case of the CFD CORC ® cable. Furthermore, the results suggest that the NZPV depends only on the applied current, similar to what has been observed previously on single REBCO tapes. These results demonstrate the successful integration of CFD REBCO tapes in CORC ® cables, without compromising their superconducting properties or their enhanced NZPV. The possibility of enhancing the NZPV of REBCO-based cables could facilitate quench detection, which remains an important issue in HTS magnet applications. This advancement also holds promise for the fault current limitation functionality of CORC ® cables.

Ben Saâd, Haïfa (ORCID:0000000216205364)↗
DOE OSTI2020

Nonlinear dynamics and dissipation of a curvilinear vortex driven by a strong time-dependent Meissner current

Here, we report numerical simulations of large-amplitude oscillations of a trapped vortex line under a strong ac magnetic field $\textit{H(t) = H}$ sin$\textit{ωt}$ parallel to the surface. The power dissipated by an oscillating vortex segment driven by the surface ac Meissner currents was calculated by taking into account the nonlinear vortex line tension, vortex mass, and a nonlinear Larkin-Ovchinnikov (LO) viscous drag coefficient $\textit{η(v)}$. We show that the LO decrease of $\textit{η(v)}$ with the vortex velocity $\textit{v}$ can radically change the field dependence of the surface resistance R i ( H ) caused by trapped vortices. At low frequencies, $R_i(H)$ exhibits a conventional increase with $\textit{H}$, but as $\textit{ω}$ increases, the surface resistance becomes a nonmonotonic function of $\textit{H}$ which decreases with $\textit{H}$ at higher fields. The effects of frequency, pin spacing, and the mean-free path $l_i$ on the field dependence of $R_i(H)$ were calculated. It is shown that, as the surface gets dirtier and $l_i$ decreases, the anomalous drop of $R_i(H)$ with $\textit{H}$ shifts to lower fields which can be much smaller than the lower critical magnetic field. Our numerical simulations also show that the LO decrease of $\textit{η(v)}$ with $\textit{v}$ can cause a vortex bending instability at high field amplitudes and frequencies, giving rise to the formation of dynamic kinks along the vortex. Measurements of $R_i(H)$ caused by sparse vortices trapped perpendicular to the surface can offer opportunities to investigate an extreme nonlinear dynamics of vortices driven by strong current densities up to the depairing limit at low temperatures. The behavior of $R_i(H)$, which can be tuned by varying the rf frequency or concentration of nonmagnetic impurities, is not masked by strong heating effects characteristic of dc or pulse transport measurements.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Experimental Comparison of Elastomeric Materials for Hydraulic Seal Durability under Reciprocating Conditions

Wave Energy Converters (WECs) rely on hydraulic Power Take-Off (PTO) systems to transform the kinetic energy of ocean waves into electricity, where hydraulic seals play a critical role in preventing fluid leakage and ensuring watertight operation under harsh marine conditions. However, seal degradation due to wear, fatigue, and corrosion remains a key challenge limiting system reliability and increasing maintenance costs. In this study, the tribological performance of elastomeric materials for hydraulic sealing applications that are based on commercial polyurethane and custom-compounded Ethylene Propylene Diene Monomer (EPDM) is investigated. Using a CSM tribometer, reciprocating wear tests were conducted in accordance with ASTM G133 under dry sliding conditions, with normal loads ranging from 3 to 10 N, and sliding speeds of 10 to 30 mm/s. Coefficient of friction (COF), specific wear rate, and worn surface morphology are evaluated to understand how load, velocity, and material composition affect seal durability. Improved wear stability and high Shore A hardness (93–94) were demonstrated by the commercial thermoplastic polyurethane (Duralast 4203 and 4758). EPDM-based compounds showed different filler-dependent behaviors, formulations filled with carbon black (E3) had the lowest wear rate (1.45 × 10?4 mm³ N?¹ m?¹) and increased load-bearing capacity, while samples modified with silica and plasticizer (E4, E5) had better conformability and damping but less hardness. Compounded EPDMs showed wider deformation zones because of their lower modulus and filler-matrix interfacial fluctuations, while commercial grades exhibited shallower wear scars (= 380 µm), according to optical profilometry. The Archard model-based lifetime estimation showed that, under 10 N and 30 mm/s, operational endurance ranged from 3.1 to 6.2 years. Because of its low wear rate and steady interfacial behavior, E3 had the longest anticipated life. For medium-load hydraulic sealing applications, these results demonstrate the possibility of customized EPDM formulations as affordable substitutes for commercial polyurethanes. They also offer a quantitative foundation for material selection and lifespan prediction.

Durability↗
DOE OSTI2021

Kink-antikink collisions and multi-bounce resonance windows in higher-order field theories

In this work, we study collisions of coherent structures in higher-order field-theoretic models, such as the Φ 8 , Φ 10 and Φ 12 ones. The main distinguishing feature, of the example models considered herein, is that the collision arises due to the long-range interacting algebraic tails of these solitary waves. We extend the approach to suitably initialize the relevant kinks, in the additional presence of finite initial velocity, in order to minimize the dispersive wave radiation potentially created by their slow spatial decay. We find that, when suitably initialized, these models still feature the multi-bounce resonance windows earlier found in models in which the kinks bear exponential tails, such as the Φ 4 and Φ 6 field theories among others. Also present is the self-similar structure of the associated windows with three- and more-bounce windows at the edges of two- and lower-bounce ones. Moreover, phenomenological but highly accurate (and predictive), scaling relations are derived for the dependence of the time between consecutive collisions and, e.g., the difference in kinetic energy between the incoming one and the critical one for one-bounces. Such scalings are traced extensively over two-bounce collision windows throughout the three models, hinting at the possibility of an analytical theory in this direction.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

TESS Hunt for Young and Maturing Exoplanets (THYME). III. A Two-planet System in the 400 Myr Ursa Major Group

Exoplanets can evolve significantly between birth and maturity, as their atmospheres, orbits, and structures are shaped by their environment. Young planets (<1 Gyr) offer an opportunity to probe the critical early stages of this evolution, where planets evolve the fastest. However, most of the known young planets orbit prohibitively faint stars. We present the discovery of two planets transiting HD 63433 (TOI 1726, TIC 130181866), a young Sun-like (M{sub ∗}=0.99±0.03) star. Through kinematics, lithium abundance, and rotation, we confirm that HD 63433 is a member of the Ursa Major moving group (τ = 414 ± 23 Myr). Based on the TESS light curve and updated stellar parameters, we estimate that the planet radii are 2.15 ± 0.10 R {sub ⊕} and 2.67 ± 0.12 R {sub ⊕}, the orbital periods are 7.11 and 20.55 days, and the orbital eccentricities are lower than about 0.2. Using High Accuracy Radial velocity Planet Searcher for the Northern hemisphere velocities, we measure the Rossiter–McLaughlin signal of the inner planet, demonstrating that the orbit is prograde. Since the host star is bright (V = 6.9), both planets are amenable to transmission spectroscopy, radial velocity measurements of their masses, and more precise determination of the stellar obliquity. This system is therefore poised to play an important role in our understanding of planetary system evolution in the first billion years after formation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Robust 3D multi-material hydrodynamics using discontinuous Galerkin methods

A high-order discontinuous Galerkin (DG) method is presented for nonequilibrium multi-material (m ≥ 2) flow with sharp interfaces. Material interfaces are reconstructed using the algebraic THINC approach, resulting in a sharp interface resolution. The system assumes stiff velocity relaxation and pressure nonequilibrium. The presented DG method uses Dubiner's orthogonal basis functions on tetrahedral elements. This results in a unique combination of sharp multimaterial interfaces and high-order accurate solutions in smooth single-material regions. A novel shock indicator based on the interface conservation condition is introduced to mark regions with discontinuities. Slope limiting techniques are applied only in these regions so that nonphysical oscillations are eliminated while maintaining high-order accuracy in smooth regions. A local projection is applied on the limited solution to ensure discrete closure law preservation. The effectiveness of this novel limiting strategy is demonstrated for complex three-dimensional multi-material problems, where robustness of the method is critical. The presented numerical problems demonstrate that more accurate and efficient multi-material solutions can be obtained by the DG method, as compared to second-order finite volume methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Urban wash-off of tire wear particles

Tire wear particles (TWPs) are an important class of microplastics due to their toxicity and abundance. Because most TWPs are generated on impervious road surfaces, urban wash-off is the critical first phase of waterborne transport from their zone of production to stormwater drainage. However, little is known about the driving factors behind their mobilization. In this study, we use a rainfall simulator to investigate how surface roughness, rainfall intensity, and surface slope affect wash-off behaviors of TWPs. We also analyze how the size and shape of mobilized TWPs change over the course of simulated storm events. We found that low surface roughness, high rainfall intensity (most significant factor), and low slope result in the most rapid conveyance of TWP load. On average, large particles (>1000 µm) travelled faster than small particles (<125 µm). Particle shape explained a very small amount of variance in TWP wash-off velocity but was found to be more important under higher surface roughness conditions. In addition to wash-off velocity, we found similar conditions controlled the percent mobilization of TWPs. Low surface roughness and high rainfall intensity resulting in higher TWP wash-off rates is consistent with mineral sediment wash-off behavior. Conversely, low surface slope and large particle size leading to faster conveyance is directly opposed to mineral sediment wash-off. Our findings suggest drag-dominated flow and that sufficient runoff depth is the most important parameter governing TWP wash-off. These findings are important first steps to understanding wash-off behaviors of TWPs and informing future modeling efforts and mitigation strategies.

13 HYDRO ENERGY↗
DOE OSTI2020

Wettability and capillary effects: Dynamics of pinch-off in unconstricted straight capillary tubes

We study the interfacial evolution of immiscible two-phase flow within a capillary tube in the partial wetting regime using direct numerical simulation. We investigate the flow patterns resulting from the displacement of a more viscous fluid by a less viscous one under a wide range of wettability conditions. We find that beyond a wettability dependent critical capillary number, a uniform displacement by a less viscous fluid can transition into a growing finger that eventually breaks up into discrete blobs by a series of pinch-off events for both wetting and non-wetting contact angles. This study validates previous experimental observations of pinch-off for wetting contact angles and extends those to non-wetting contact angles. We find that the blob length increases with the capillary number. We observe that the time between consecutive pinch-off events decreases with the capillary number and is greater for more wetting conditions in the displaced phase. Here, we further show that the blob separation distance as a function of the difference between the inlet velocity and the contact line speed collapses into two monotonically decreasing curves for wetting and non-wetting contact angles. For the phase separation in the form of pinch-off, this work provides a quantitative study of the emerging length and time scales and their dependence on the wettability conditions, capillary effects, and viscous forces.

42 ENGINEERING↗
DOE OSTI2025

Sample glue layer investigation and mitigation for laser induced prompt impulse experiments

Understanding longer timescale material reactions under dynamic stress loading is critical for applications in materials engineering, shock physics, and planetary science. Prompt impulse experiments generate lower pressures since the ablator—the material first removed by the laser—is thicker and farther from the diagnostic plane, capturing aggregate material responses from the initial shock wave, rarefaction waves, and later time effects. This complexity demands thorough material characterization and simulation support. Since traditional sample construction is specific to supported shock experiments, designing prompt impulse experiments requires reconsideration around target design and sample engineering. Here, we present sample preparation techniques, experimental investigations, and theoretical simulations to investigate glue layer impacts, aiming to standardize samples for consistent data at lower laser fluences. We find that glue layers <30 μm have a minimal impact on peak velocity and pulse shape. The peak velocity scales linearly with glue layer thickness until a glue layer of 75 μm. For glue layers >75 μm, the peak velocity no longer scales with thickness; however, the pulse shape continues to degrade as described by simulations.

Lasers↗