Search NASA⌕ Search

SEARCH · Search NASA

Results for “Very viscous flows”

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 235 records · Page 13

A method for predicting failure statistics for steady state elevated temperature structural components

This paper presents the initial development of a high temperature life prediction method that accounts for the variability in the material properties of Grade 91 steel. The method accounts for material variability by fitting a variable 3-parameter Weibull distribution to experimental rupture data and accounts for the variability of creep deformation on the steady-state stresses via a Monte Carlo approach. To ensure reasonable computational times, the model represents the material as an extremely viscous Stokes fluid with a non-Newtonian viscosity, therefore solving the stress relaxation problem with a steady, static, instead of transient, analysis. Furthermore, the complete statistical analysis combines this model for creep deformation with a probabilistic model for creep rupture to evaluate the probability of premature failure for a set of sample problems, comparing the predicted failure statistics to the design life predicted by the ASME Boiler and Pressure Vessel Code rules.

42 ENGINEERING↗

Interplay of the forces governing steroid hormone micropollutant adsorption in vertically-aligned carbon nanotube membrane nanopores

Vertically-aligned carbon nanotube (VaCNT) membranes allow water to conduct rapidly at low pressures and open up the possibility for water purification and desalination, although the ultralow viscous stress in hydrophobic and low-tortuosity nanopores prevents surface interactions with contaminants. In this experimental investigation, steroid hormone micropollutant adsorption by VaCNT membranes is quantified and explained via the interplay of the hydrodynamic drag and friction forces acting on the hormone, and the adhesive and repulsive forces between the hormone and the inner carbon nanotube wall. It is concluded that a drag force above 2.2 × 10 —3 pN overcomes the friction force resulting in insignificant adsorption, whereas lowering the drag force from 2.2 × 10 —3 to 4.3 × 10 —4 pN increases the adsorbed mass of hormones from zero to 0.4 ng cm —2 . At a low drag force of 1.6 × 10 —3 pN, the adsorbed mass of four hormones is correlated with the hormone–wall adhesive (van der Waals) force. These findings explain micropollutant adsorption in nanopores via the forces acting on the micropollutant along and perpendicular to the flow, which can be exploited for selectivity.

36 MATERIALS SCIENCE↗

High-Temperature Heat Transfer Fluid Circulator for Concentrating Solar Power Systems

Concentrating Solar Power (CSP) plants with Thermal Energy Storage (TES), require compressors capable of circulating supercritical carbon dioxide (sCO2) heat transfer fluids (HTFs) at temperatures from 500 to 700+°C and pressures from 8 to 25 MPa (80 to 250 Bar or 1160 to 3626 psia) in order to achieve thermal-to-electric efficiencies >50%. A review of papers relating to SCO2 compressors reveals that in almost every case the conditions being considered either allow for cooling of the SCO2 working fluid or are considering operating conditions close to the critical point at a temperature and pressure of approximately 31°C and 1070 psi). While the technology reported in these papers is relevant to most turbomachinery, the current CSP gas pathway circulator operating pressures are well above the CO2 critical point and do not allow for cooling of the main HTF flow. The purpose of this paper was to describe the key system issues affecting the design of a HTF circulator that could be scaled from a small size research and proof of concept circulator compressor to sizes supporting power plants producint upwards of 100 to 300 Mega Watts of electrical power. Due to the high temperatures and pressures the research scale HTF circulator was designed to be hermetically sealed to prevent leakage of the super critical carbon dioxide to the environment. In order to achieve hermetic sealing the drive motor and centrifugal impeller were fully integrated in a single shaft system configuration. From a preliminary assessment of compression power as a function of operating speed and impeller diameter it was determined that a smaller, higher speed impeller would result in a more efficient circulator. However, due to the high pressure and corresponding density of the sCO2 coupled with high operating speeds, large windage losses (i.e., viscous drag on the high speed rotating shaft) would occur which could result in a substantial reduction in efficiency of the system. Therefore tradoff studies were complerted to identify the operating speeds and system component sizes to achieve high efficiency. Based on the trade studies, the final design was a 55 mm diamter compressor impeller operating at speeds to 85,000 rpm. The unit will be built and tested in the future. Key Words: Super Critical Carbon Dioxide, Heat Transfer Fluid, Turbomachinery, Compressor, Circulator, Critical Point, Concentrating Solar Power, Solar Energy

Super Critical Carbon Dioxide, Heat Transfer Fluid↗

Multiphase Processing of the Water-Soluble and Insoluble Phases of Biomass Burning Organic Aerosol

Biomass burning is one of the most significant sources of organic aerosol in the atmosphere. Biomass burning organic aerosol (BBOA) has been observed to undergo liquid– liquid phase separation (LLPS) to give core–shell morphology with the hydrophobic phase encapsulating the hydrophilic phase, potentially impacting the evolution of light-absorbing components, i.e., brown carbon (BrC), through multiphase processes. Here, we demonstrate how multiphase processing differs between the watersoluble (i.e., hydrophilic) and insoluble (i.e., hydrophobic) phases of BBOA in terms of reactive uptake of ozone in a coated-wall flow tube. Effects of relative humidity (RH) and ultraviolet (UV) irradiation were investigated. Experimental timeseries were used to inform simulations using multilayer kinetic modeling. Among non-irradiated thin films, the uptake coefficient was greatest for the water-soluble phase at 75% RH (3 × 10 –5 , corresponding to a diffusion coefficient of BrC, D BrC , of 3 × 10 –9 cm 2 s –1 ) and least for the same phase at 0% RH (1 × 10 –5 , corresponding to D BrC of 1 × 10 –10 cm 2 s –1 ). The uptake coefficient for the water-insoluble phase fell between these two (about 1.5 × 10 –5 ), regardless of RH, and the corresponding D BrC increased only slightly (8 × 10 –10 cm 2 s –1 at 0% RH to 9 × 10 –10 cm 2 s –1 at 75% RH). The uptake coefficients of both phases at 0% RH decreased significantly after UV irradiation, consistent with a transition from viscous liquid to solid and supported by qualitative microscopy observations. Modeling multiphase ozone oxidation of primary BrC components in the atmosphere demonstrated, first, that LLPS may extend the lifetime of water-soluble BBOA encapsulated by water-insoluble species by a factor of 1.5 at moderate to high RH and, also, that UV irradiation may extend the lifetime of both phases by more than a factor of 2.5.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Comparison of momentum transport models for numerical relativity

The main problems of nonvacuum numerical relativity, compact binary mergers and stellar collapse, involve hydromagnetic instabilities and turbulent flows, so that kinetic energy at small scales leads to mean effects at large scale that drive the secular evolution. Notable among these effects is momentum transport. We investigate two models of this transport effect, a relativistic Navier-Stokes system and a turbulent mean stress model, that are similar to all of the prescriptions that have been attempted to date for treating subgrid effects on binary neutron star mergers and their aftermath. Our investigation involves both stability analysis and numerical experimentation on star and disk systems. We also begin the investigation of the effects of particle and heat transport on postmerger simulations. Here, we find that correct handling of turbulent heating is crucial for avoiding unphysical instabilities. Given such appropriate handling, the evolution of a differentially rotating star and the accretion rate of a disk are reassuringly insensitive to the choice of prescription. However, disk outflows can be sensitive to the choice of method, even for the same effective viscous strength. We also consider the effects of eddy diffusion in the evolution of an accretion disk and show that it can interestingly affect the composition of outflows.

79 ASTRONOMY AND ASTROPHYSICS↗

Oil Recovery Prediction for Polymer Flood Field Test of Heavy Oil on Alaska North Slope Via Machine Assisted Reservoir Simulation

The first ever polymer flood field pilot to enhance the recovery of heavy oils on the Alaska North Slope is ongoing. This study constructs and calibrates a reservoir simulation model to predict the oil recovery performance of the pilot through machine-assisted reservoir simulation techniques. To replicate the early water breakthrough observed during waterflooding, transmissibility contrasts are introduced into the simulation model, forcing viscous fingering effects. In the ensuing polymer flood, these transmissibility contrasts are reduced to replicate the restoration of injection conformance during polymer flooding, as indicated by a significant decrease in water cut. Later, transmissibility contrasts are reinstated to replicate a water surge event observed in one of the producing wells during polymer flooding. This event may represent decreased injection conformance from fracture overextension; its anticipated occurrence in the other production well is included in the final forecast. The definition of polymer retention in the simulator incorporates the tailing effect reported in laboratory studies; this tailing effect is useful to the simultaneous history match of producing water cut and produced polymer concentration. The top 24 best-matched simulation models produced at each stage of the history matching process are used to forecast oil recovery. The final forecast clearly demonstrates that polymer flooding significantly increases the heavy oil production for this field pilot compared to waterflooding alone. This exercise displays that a simulation model is only valid for prediction if flow behavior in the reservoir remains consistent with that observed during the history matched period. Critically, this means that a simulation model calibrated for waterflooding may not fully capture the benefits of an enhanced oil recovery process such as polymer flooding. Therefore, caution is recommended in using basic waterflood simulation models to scope potential enhanced oil recovery projects.

Keith, Cody Douglas↗

Effects of dissipative baryon current in heavy-ion collisions at RHIC-BES energies

The CLVisc (3+1)D viscous hydrodynamic model is extended to include the equation of net baryon conservation and the Israel-Stewart-like equations for dissipative baryon current. Using the NEOSB equation of state, we simulate the dynamical evolution and collectivity of the quark-gluon plasma with finite chemical potential, assuming smooth energy density and net baryon density distributions at the initial proper time. Numerical results are shown for the impact of net-baryon dissipation on particle yields and p T spectra in heavy-ion collisions at beam energy scan energies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Fluctuation dynamo in a weakly collisional plasma

The turbulent amplification of cosmic magnetic fields depends upon the material properties of the host plasma. In many hot, dilute astrophysical systems, such as the intracluster medium (ICM) of galaxy clusters, the rarity of particle–particle collisions allows departures from local thermodynamic equilibrium. These departures – pressure anisotropies – exert anisotropic viscous stresses on the plasma motions that inhibit their ability to stretch magnetic-field lines. We present an extensive numerical study of the fluctuation dynamo in a weakly collisional plasma using magnetohydrodynamic (MHD) equations endowed with a field-parallel viscous (Braginskii) stress. When the stress is limited to values consistent with a pressure anisotropy regulated by firehose and mirror instabilities, the Braginskii-MHD dynamo largely resembles its MHD counterpart, particularly when the magnetic field is dynamically weak. If instead the parallel viscous stress is left unabated – a situation relevant to recent kinetic simulations of the fluctuation dynamo and, we argue, to the early stages of the dynamo in a magnetized ICM – the dynamo changes its character, amplifying the magnetic field while exhibiting many characteristics reminiscent of the saturated state of the large-Prandtl-number ( ${Pm}\gtrsim {1}$ ) MHD dynamo. We construct an analytic model for the Braginskii-MHD dynamo in this regime, which successfully matches simulated dynamo growth rates and magnetic-energy spectra. A prediction of this model, confirmed by our numerical simulations, is that a Braginskii-MHD plasma without pressure-anisotropy limiters will not support a dynamo if the ratio of perpendicular and parallel viscosities is too small. This ratio reflects the relative allowed rates of field-line stretching and mixing, the latter of which promotes resistive dissipation of the magnetic field. Finally, in all cases that do exhibit a viable dynamo, the generated magnetic field is organized into folds that persist into the saturated state and bias the chaotic flow to acquire a scale-dependent spectral anisotropy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Vortex rings from high energy central 𝑝 + 𝐴 collisions

Relativistic p + A collisions may produce droplets of quark gluon plasma (QGP) that quickly develop a toroidal vortex structure similar to that of an expanding smoke ring. We present viscous relativistic hydrodynamic calculations of ultra-central p + A collisions and develop an experimental observable to probe the structure, correlating the polarization and momentum of hyperons emitted from the collision. This effect is robust against changes in the definition of vorticity used to calculate the polarization. Experiments at RHIC and LHC may test the existence and strength of the vortex toroids, bringing new evidence to bear on the question of collectivity in the smallest QGP droplets.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Heterogeneous oxidation of amorphous organic aerosol surrogates by O 3 , NO 3 , and OH at typical tropospheric temperatures

Typical tropospheric temperatures render possible phase states of amorphous organic aerosol (OA) particles of solid, semisolid, and liquid. This will affect the multiphase oxidation kinetics involving the organic condensed-phase and gaseous oxidants and radicals. To quantify this effect, we determined the reactive uptake coefficients (γ) of O 3 , NO 3 , and OH by substrate films composed of single and binary OA surrogate species under dry conditions for temperatures from 213 to 313 K. A temperature-controlled coated-wall flow reactor coupled to a chemical ionization mass spectrometer was applied to determine γ with consideration of gas diffusion transport limitation and gas flow entrance effects, which can impact heterogeneous reaction kinetics. The phase state of the organic substrates was probed via the poke-flow technique, allowing the estimation of the substrates' glass transition temperatures. γ values for O 3 and OH uptake to a canola oil substrate, NO 3 uptake to a levoglucosan and a levoglucosan / xylitol substrate, and OH uptake to a glucose and glucose / 1,2,6-hexanetriol substrate have been determined as a function of temperature. We observed the greatest changes in γ with temperature for substrates that experienced the largest changes in viscosity as a result of a solid-to-liquid phase transition. Organic substrates that maintain a semisolid or solid phase state and as such a relatively higher viscosity do not display large variations in heterogeneous reactivity. From 213 to 293 K, γ values of O 3 with canola oil, of NO 3 with a levoglucosan / xylitol mixture, and of OH with a glucose / 1,2,6-hexanetriol mixture and canola oil, increase by about a factor of 34, 3, 2, and 5, respectively, due to a solid-to-liquid phase transition of the substrate. These results demonstrate that the surface and bulk lifetime of the OA surrogate species can significantly increase due to the slowed heterogeneous kinetics when OA species are solid or highly viscous in the middle and upper troposphere. This experimental study will further our understanding of the chemical evolution of OA particles with subsequent important consequences for source apportionment, air quality, and climate.

54 ENVIRONMENTAL SCIENCES↗

Elliptic flow of charged particles at midrapidity relative to the spectator plane in Pb–Pb and Xe–Xe collisions

Measurements of the elliptic flow coefficient relative to the collision plane defined by the spectator neutrons v 2 {$Ψ_{SP}$} in collisions of Pb ions at center-of-mass energy per nucleon–nucleon pair $\sqrt{s_{NN}}=$ 2.76 TeV and Xe ions at $\sqrt{s_{NN}}=$ 5.44 TeV are reported. The results are presented for charged particles produced at midrapidity as a function of centrality and transverse momentum for the 5–70% and 0.2–6 GeV/c ranges, respectively. The ratio between v 2 {$Ψ_{SP}$} and the elliptic flow coefficient relative to the participant plane v 2 {4}, estimated using four-particle correlations, deviates by up to 20% from unity depending on centrality. This observation differs strongly from the magnitude of the corresponding eccentricity ratios predicted by the T R ENTo and the elliptic power models of initial state fluctuations that are tuned to describe the participant plane anisotropies. The differences can be interpreted as a decorrelation of the neutron spectator plane and the reaction plane because of fragmentation of the remnants from the colliding nuclei, which points to an incompleteness of current models describing the initial state fluctuations. A significant transverse momentum dependence of the ratio v 2 {$Ψ_{SP}$} /v 2 {4} is observed in all but the most central collisions, which may help to understand whether momentum anisotropies at low and intermediate transverse momentum have a common origin in initial state fluctuations. The ratios of v 2 {$Ψ_{SP}$} and v 2 {4} to the corresponding initial state eccentricities for Xe–Xe and Pb–Pb collisions at similar initial entropy density show a difference of (7.0 ± 0.9)% with an additional variation of +1.8% when including RHIC data in the T R ENTo parameter extraction. These observations provide new experimental constraints for viscous effects in the hydrodynamic modeling of the expanding quark–gluon plasma produced in heavy-ion collisions at the LHC.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Running the gamut of high energy nuclear collisions

We present calculations of bulk properties and multiparticle correlations in a large variety of collision systems within a hybrid formalism consisting of IP-Glasma initial conditions, Music viscous relativistic hydrodynamics, and UrQMD microscopic hadronic transport. In particular, we study heavy ion collisions at the Large Hadron Collider (LHC), including Pb + Pb, Xe + Xe, and O + O collisions, and Au + Au, U + U, Ru + Ru, Zr + Zr, and O + O collisions at the Relativistic Heavy Ion Collider (RHIC). We further study asymmetric systems, including p + Au, d + Au, 3 He + Au, and p + Pb collisions at various energies as well as p+p collisions at 0.5 and 13 TeV. Here, we describe experimental observables in all heavy ion systems well with one fixed set of parameters, validating the energy and system dependence of the framework. As a result, many observables in the smaller systems are also well described, although they test the limits of the model.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Experimental Investigation of Polymer Induced Fouling of Heater Tubes in the First-Ever Polymer Flood Pilot on Alaska North Slope

Mineral fouling in heat exchangers has been extensively investigated by researchers in recent times. The oil and gas industry has a long history of fouling issues in production systems as a result of produced fluids treatment. Due to decline in production rates in oilfields new technologies are being developed and field tested in pilots. Polymer flooding is one such technology that involves addition of polymers to injection fluids to enhance oil production. A polymer flood pilot has been set up in the Schrader Bluff viscous oil reservoir at Milne Point field on the Alaska North Slope (ANS). The results from the pilot are encouraging, however a major concern of the operator is the influence of polymer on the production system after breakthrough, especially the fouling in heat exchangers. This study investigates the propensity of polymer fouling on the heater tubes as a function of different variables, with the ultimate goal of determining safe and efficient operating conditions. This work applies a multi-experimental approach to study the severity of polymer-induced fouling in both dynamic and static states of produced fluids as well as studying the stability of polymer solutions at different temperatures. A unique experimental setup was designed and developed in-house to simulate the fouling process on the heating tube. The influence of heating tube skin temperature, tube material, and polymer concentration on fouling tendency was investigated. Each test was run five times with the same tube, and in each run, the freshly prepared synthetic brine and polymer solution was heated from 77°F to 122°F to mimic field-operating conditions. The heating time and fouling amount were recorded for each run. Dynamic Scale Loop (DSL) tests were conducted to study fouling due to polymer at different temperatures (165°F to 350°F) in a dynamic state of fluid flow where the fluids mimic the residence time of fluids in the heat exchanger on the field pilot. Cloud point measurement has also been conducted to find the critical temperature at which the polymer in solution becomes unstable and precipitates out. The morphology and composition of the deposit samples were analyzed by environmental scanning electron microscopy (ESEM) and X-ray diffraction (XRD), respectively. It was found that the presence of polymer in produced fluids would aggravate the fouling issues on both carbon steel and stainless-steel surfaces at all tested skin temperatures. Only higher skin temperatures of 250°F and 350°F could cause polymer-induced fouling issues on the copper tube surface, and the fouling tendency increased with polymer concentration. At the lower skin temperatures of 165°F, no polymer-induced fouling was identified on the copper tube. A critical temperature that is related to the cloud point of the polymer solution was believed to exist, below which polymer-induced fouling would not occur, and only mineral scale was deposited but above which the polymer would aggravate the fouling issue. The cloud point of the tested polymer solution was determined to be between 220°F and 230°F. In the DSL tests it was found that at higher skin temperatures of 250°F and 350°F tube blocking was observed in the DSL tests whereas the tests at 165°F and 200°F did not show any tube blocking in the same time period. These experiments also manifested the influence of cloud point of the solution as deposit rate increased significantly in both carbon steel and stainless-steel tubes when the skin temperature was higher than the solution cloud point. The results of this study have provided guidance to the operator for the field-operations.

Dhaliwal, Anshul↗

Jet cone radius dependence of R AA and v 2 at PbPb 5.02 TeV from JEWEL+T R ENTo+v-USPhydro

We combine, for the first time, event-by-event T R ENTo initial conditions with the relativistic viscous hydrodynamic model v-USPhydro and the Monte Carlo event generator JEWEL to make predictions for the nuclear modification factor R AA and jet azimuthal anisotropies v n { 2 } in $\sqrt{s_{NN}}$ = 5.02 TeV PbPb collisions for multiple centralities and values of the jet cone radius R. The R-dependence of R AA and v 2 { 2 } strongly depends on the presence of recoiling scattering centers. We find a small jet v 3 { 2 } in mid-central collisions and consistent results in wide jet p T regions and centralities with ATLAS data.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

New scaling law for turbulent boundary layers with high surface mass transfer

A fundamental understanding of the phenomena occurring in the turbulent boundary layer in the presence of surface blowing is limited, and considerable disagreements persist even in describing primary quantities, such as the boundary layer profile. The theories based on the linear boundary layer equations show that the thickness of the sublayer increases in the presence of surface blowing; therefore, the viscous sublayer and law of the wall modify. In this study, direct numerical simulations (DNS) of turbulent boundary layers with uniform surface mass transfers are carried out in order to scale the velocity profile. Emphasis is placed on moderate to high mass transfer rates, which are relevant to the most common hybrid rockets configuration. DNS data are used to establish a functional law of the wall and a law of wake by means of the relation between the wall shear stress and surface mass transfer. Analysis of the mean kinetic energy budget shows that the magnitude of turbulent kinetic energy increases by surface mass transfer, and the production rate extends significantly in the inner layer as the injection rate increases. DNS data of various surface blowing are used to complete the closure of turbulence kinetic energy equation and develop an eddy viscosity model. The predicted turbulent kinetic energy and eddy viscosity agree with DNS data for moderate to high blowing rates.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Computing the viscous effect in early-time drop impact dynamics

The impact of a liquid drop on a solid surface involves many intertwined physical effects, and is influenced by drop velocity, surface tension, ambient pressure and liquid viscosity, among others. Experiments by Kolinski et al. ( Phys. Rev. Lett. , vol. 112, no. 13, 2014 b , p. 134501) show that the liquid–air interface begins to deviate away from the solid surface even before contact. They found that the lift-off of the interface starts at a critical time that scales with the square root of the kinematic viscosity of the liquid. To understand this, we study the approach of a liquid drop towards a solid surface in the presence of an intervening gas layer. We take a numerical approach to solve the Navier–Stokes equations for the liquid, coupled to the compressible lubrication equations for the gas, in two dimensions. With this approach, we recover the experimentally captured early time effect of liquid viscosity on the drop impact, but our results show that lift-off time and liquid kinematic viscosity have a more complex dependence than the square-root scaling relationship. We also predict the effect of interfacial tension at the liquid–gas interface on the drop impact, showing that it mediates the lift-off behaviour.

42 ENGINEERING↗

Co-flow injection for serial crystallography at X-ray free-electron lasers

Serial femtosecond crystallography (SFX) is a powerful technique that exploits X-ray free-electron lasers to determine the structure of macromolecules at room temperature. Despite the impressive exposition of structural details with this novel crystallographic approach, the methods currently available to introduce crystals into the path of the X-ray beam sometimes exhibit serious drawbacks. Samples requiring liquid injection of crystal slurries consume large quantities of crystals (at times up to a gram of protein per data set), may not be compatible with vacuum configurations on beamlines or provide a high background due to additional sheathing liquids present during the injection. Proposed and characterized here is the use of an immiscible inert oil phase to supplement the flow of sample in a hybrid microfluidic 3D-printed co-flow device. Co-flow generation is reported with sample and oil phases flowing in parallel, resulting in stable injection conditions for two different resin materials experimentally. A numerical model is presented that adequately predicts these flow-rate conditions. The co-flow generating devices reduce crystal clogging effects, have the potential to conserve protein crystal samples up to 95% and will allow degradation-free light-induced time-resolved SFX.

3D printing↗

Universal Anomalous Diffusion of Quantized Vortices in Ultraquantum Turbulence

In classical viscous fluids, turbulent eddies are known to be responsible for the rapid spreading of embedded particles. However, in inviscid quantum fluids where the turbulence is induced by a chaotic tangle of quantized vortices, dispersion of the particles can be achieved via a nonclassical mechanism, i.e., their binding to the evolving vortices. However, knowledge on how the vortices diffuse and spread in quantum-fluid turbulence is very limited, especially for the so-called ultraquantum turbulence (UQT) generated by a random tangle of vortices. Here we report a systematic numerical study of the apparent diffusion of vortices in UQT in superfluid helium-4 using the full Biot-Savart simulation. We reveal that the vortices in the superfluid exhibit a universal anomalous diffusion (superdiffusion) at small times, which transits to normal diffusion at large times. This behavior is found to be the result of a generic scaling property of the vortex velocity. Our simulation at finite temperatures also nicely reproduces recent experimental observations. Lastly, the knowledge obtained from this study may form the base for understanding turbulent transport and universal vortex dynamics in various quantum fluids.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗