Search NASA⌕ Search

SEARCH · Search NASA

Results for “power-law fluid”

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.

69 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Superdiffusion of quantized vortices uncovering scaling laws in quantum turbulence

Significance Quantum turbulence (QT) can appear in the presence of a chaotic tangle of quantized vortices in various quantum-fluid systems, including atomic Bose–Einstein condensates, superfluid helium, superfluid neutron stars, etc. Insights into the generic scaling behaviors of tangled vortices are crucial in developing an advanced statistical model of QT. By tracking tracer particles trapped on vortices in superfluid 4 He, we report the observation of an apparent superdiffusion of the vortices in QT. Our analysis shows that this superdiffusion is not due to Lévy flights, i.e., long-distance hops that are responsible for superdiffusion of random walkers. Instead, a power-law scaling of the vortex–velocity correlation is identified as the cause. This finding may motivate extensive future research on hidden scaling laws in QT.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Shivalik Plasma Device-I, a glow discharge device to study the collective dynamics of dusty plasma

Using a uniquely configured glow discharge-based Shivalik Plasma Device-I, we demonstrate a variety of collective phenomena in dusty plasma away from the glow discharge region. The cylindrical glass device produces plasma using parallel disc-shaped electrodes with a smaller anode size than the cathode. The dust microparticles are initially sprinkled over the grounded cathode. These particles acquire a significant negative charge upon plasma formation, resulting in their levitation due to the balance between the Coulomb force and gravity. The new device supports the levitation of a big-sized (10 × 8 × 5 cm 3 ) three-dimensional dust cloud over the glass surface. It contrasts the dusty plasma formations in-between electrodes reported earlier. As the discharge voltage varies from high to low, the dust cloud travels from over the glass surface to between the electrodes. A complex interplay of dust void over the cathode, a sharp density gradient, and gravity lead to self-excitation of collective dust phenomena. It includes dust density waves (phase velocity, v ph ~4 cm/s), dust cloud oscillation (frequency, f = 5 Hz), sheared flow (flow velocity, v f ~ 1 cm/s), and multiple-sized dust vortices. These dust vortices provided an excellent platform for studying turbulent mixing phenomena. The power spectrum analysis agreed with two-dimensional Kolmogorov power-law scaling. This is an ideal dusty plasma apparatus where we can create or move the dust cloud to a location of choice from the glass surface to in-between the electrodes and excite one among many collective dust dynamics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Adiabatic behavior of an elliptical vortex in a time-dependent external strain flow

The adiabatic behavior of two-dimensional (2D) fluid vortices subject to slowly varying, time-dependent external strain flows is studied theoretically, and experimentally using magnetized non-neutral electron plasmas. Here the E × B drift dynamics of the electrons perpendicular to the magnetic field are analogous to the motion of fluid vorticity under the 2D Euler equations describing ideal fluids. A low-dimensional elliptical vortex patch model is used along with a WKB approximation to derive a formula for the breaking of an adiabatic invariant due to the changing external flow. The invariant is interpreted as the amplitude of a perturbative oscillation about a stable fixed point corresponding to a steady elliptical vortex. Smooth, hyperbolic tangent and piecewise linear ramp functions are considered for the external strain time dependence. Standard exponential breaking is observed in the former case, whereas the latter exhibits a power-law breaking curve with periodic modulation. It is found that a driving term in the equations of motion contributes the majority of the breaking whereas the frequency variation plays a weaker, but significant role. The experimental data agree closely with the theoretical model. The most significant deviation is due to inviscid damping behavior associated with the smooth edges of the experimental vortices, which tends to reduce the amplitude of the oscillation. Finally, the results are compared and contrasted with other related experimental, numerical, and theoretical work.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Spatial Intermittency of Particle Distribution in Relativistic Plasma Turbulence

Abstract Relativistic magnetically dominated turbulence is an efficient engine for particle acceleration in a collisionless plasma. Ultrarelativistic particles accelerated by interactions with turbulent fluctuations form nonthermal power-law distribution functions in the momentum (or energy) space, f ( γ ) d γ ∝ γ − α d γ , where γ is the Lorenz factor. We argue that in addition to exhibiting non-Gaussian distributions over energies, particles energized by relativistic turbulence also become highly intermittent in space. Based on particle-in-cell numerical simulations and phenomenological modeling, we propose that the bulk plasma density has lognormal statistics, while the density of the accelerated particles, n , has a power-law distribution function, P ( n ) dn ∝ n − β dn . We argue that the scaling exponents are related as β ≈ α + 1, which is broadly consistent with numerical simulations. Non-space-filling, intermittent distributions of plasma density and energy fluctuations may have implications for plasma heating and for radiation produced by relativistic turbulence.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Phase separation and ripening in a viscoelastic gel

The process of phase separation in elastic solids and viscous fluids is of fundamental importance to the stability and function of soft materials. We explore the dynamics of phase separation and domain growth in a viscoelastic material such as a polymer gel. Using analytical theory and Monte Carlo simulations, we report a domain growth regime in which the domain size increases algebraically with a ripening exponent α that depends on the viscoelastic properties of the material. For a prototypical Maxwell material, we obtain α = 1 , which is markedly different from the well-known Ostwald ripening process with α = 1 / 3 . We generalize our theory to systems with arbitrary power-law relaxation behavior and discuss our findings in the context of the long-term stability of materials as well as recent experimental results on phase separation in cross-linked networks and cytoskeleton.

Science & Technology - Other Topics↗
DOE OSTI2021

Yielding Behavior of Bottlebrush and Linear Block Copolymers

Block copolymers can exhibit a pronounced yield stress, but the impact of molecular architecture, chemistry, and self-assembly on macroscopic rheology remains poorly understood. Here, we study the linear-viscoelastic and yield-stress fluid behavior of two architectures—bottlebrush copolymers (with statistical or blocky sequences) and linear diblocks—that self-assemble into body-centered cubic (BCC) spheres and hexagonally close-packed cylinders (HEX). The dynamic properties of these polymers were probed by oscillatory frequency and amplitude sweeps at temperatures well below the order–disorder transition (TODT) to furnish insights into the yielding transition. All BCC-forming polymers have a similar signature of yielding: smaller yield strains (γ yBCC ≈ 0.053 < γ y,HEX ≈ 0.18), sharper solid–liquid transitions, and better reversibility than HEX. Statistical bottlebrushes show the most frequency-independent structural modulus (G 0 ) and no signs of defect relaxation. A simple power-law relationship captures the dependence of the normalized structural modulus (G 0 /RT) on the inter-micelle distance (d) across different architectures and morphologies [G 0 /(RT) = 1.31 × 10 4 (nm 2.6 mol/m 3 ) d –2.6 ]. These studies establish quantitative structure–property relationships that are relevant in contemporary applications, for example, extrusion-based 3D printing.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Transition to turbulence in viscoelastic channel flow of dilute polymer solutions

The transition to turbulence in a plane Poiseuille flow of dilute polymer solutions is studied by direct numerical simulations of a finitely extensible nonlinear elastic fluid with the Peterlin closure. The range of Reynolds number ($Re$)$2000 \le Re \le 5000$is studied but with the same level of elasticity in viscoelastic flows. The evolution of a finite-amplitude perturbation and its effects on the transition dynamics are investigated. A viscoelastic flow begins transition at an earlier time than its Newtonian counterparts, but the transition time appears to be insensitive to polymer concentration in the dilute and semi-dilute regimes studied. Increasing polymer concentration, however, decreases the maximum attainable energy growth during the transition process. The critical or minimum perturbation amplitude required to trigger transition is computed. Interestingly, both Newtonian and viscoelastic flows follow almost the same power-law scaling of$Re^\gamma$with the critical exponent$\gamma \approx -1.25$, which is in close agreement with previous studies. However, a shift downward is observed for viscoelastic flow, suggesting that smaller perturbation amplitudes are required for the transition. A mechanism of the early transition is investigated by the evolution of wall-normal and spanwise velocity fluctuations and flow structure. The early growth of these fluctuations and the formation of quasi-streamwise vortices around low-speed streaks are promoted by polymers, hence causing an early transition. These vortical structures are found to support the critical exponent$\gamma \approx -1.25$. Once the transition process is completed, polymers play a role in dampening the wall-normal and spanwise velocity fluctuations and vortices to attain a drag-reduced state in viscoelastic turbulent flows.

Mechanics↗
DOE OSTI2021

Structure and Shape of Surface-Mediated Assembly of Surfactants

Achieving controls on the self-assembly of cationic and block copolymer micelles has significant implications for advancing novel material systems for energy, environmental, and biological applications and for tuning fluid flow behavior in subsurface geologic environments. While it is known that the aggregation of micelles is influenced by the composition of the surrounding fluid, the influence of a solid interface has been less explored. In this study, we probe the organization of cetyltrimethylammonium bromide (CTAB) micelles in the absence and presence of the Pluronic P123 block copolymer and quartz substrate using transmission and grazing-incidence small-angle X-ray scattering measurements and classical molecular dynamics (MD) simulations. In the absence of the quartz interface, CTAB with and without P123 molecules assemble as ellipsoid core–shell micelles. Densification of the core on adding P123 is noted from the decrease in the core radius from 15.9 to 14.6 Å. Here, the presence of a quartz substrate causes the micelles to elongate, which is noted by the emergence of a power-law slope in the low q region (<0.02 Å –1 ). Moreover, in the presence of both P123 and the quartz substrate, the micelle shape changes from ellipsoid core–shell to cylindrical core–shell, and a significantly higher number of aggregates are formed. The higher number of aggregates, and faster aggregation kinetics are linked to the organization of the solvent structures as noted from the MD simulations. These findings demonstrate the effects of adding a block copolymer and including a quartz substrate on the self-assembly of cationic (CTAB) micelles, which can be used to inform the design of functional material systems that harness surfactant chemistries.

42 ENGINEERING↗
DOE OSTI2021

Pulsating Poiseuille flow of a cement slurry

In this paper, we investigate the pulsatile flow of a cement slurry in a pipe. The constitutive relation for the viscous stress tensor is based on the power-law model, where the shear viscosity not only depends on the shear rate but also the volume fraction of the cement particles. To solve for the volume fraction field, a convection-diffusion equation is used. The dimensionless form of the governing equations and the boundary conditions are solved numerically at different pulsatile cycles. Further, a parametric study is conducted to study the effect of different dimensionless numbers at different cycles on the velocity and the volume fraction profiles.

36 MATERIALS SCIENCE↗
DOE OSTI2024

On the subject of large-scale pool fires and turbulent boundary layer interactions

The role to which a realistic inflow turbulent boundary layer (TBL) influences transient and mean large-scale pool fire quantities of interest (QoIs) is numerically investigated. High-fidelity, low-Mach large-eddy simulations that activate low-dissipation, unstructured numerics are conducted using an unsteady flamelet combustion modeling approach with mutiphysics coupling to soot and participating media radiation transport. Three inlet profile configurations are exercised for a large-scale, high-aspect rectangular pool that is oriented perpendicular to the flow direction: a time-varying, TBL inflow profile obtained from a periodic precursor simulation, the time-mean of the transient TBL, and a steady power-law inflow profile that replicates the mean TBL crosswind velocity of 10.0 m/s at a vertical height of 10 m. Results include both qualitative transient flame evolution and quantitative flame shape with ground-level temperature and convective/radiative heat flux profiles. While transient fire events, which are driven by burst-sweep TBL coupling, such as blow-off and reattachment are vastly different in the TBL case (contributing to increased root mean square QoI fluctuation prediction and disparate flame lengths), mean surface QoI magnitudes are similar. Quadrant analysis demonstrates that the TBL configuration modifies burst-sweep phenomena at windward pool locations, while leeward recovery is found. In conclusion, positive fluctuations of convective heat flux correlate with fast moving fluid away from the pool surface due to intermittent combustion events.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1985

Broadening of looplike solar coronal transients

The broadening of the tops of nine outwardly expanding looplike coronal mass-ejection transients observed during 1973-1974 and 1980 has been examined. Five of the nine transients exhibit increased breadth with time (height); although the rates vary widely, no event shows an increase which can be fitted by a power-law form with a greater-than-linear increase with height. The breadth of four events decreases with height. When all nine events are considered together, the transient breadth with height is fitted with an expression h varies as R exp 0.48 with large scatter. In addition, measurements of the total observed loop length of the nine transients show that this length L varies as R exp 1.30. Both of these results are contrary to the predictions of one model of loop expulsion as a result of an assumed azimuthal magnetic field gradient. Finally, when combined with the facts that looplike mass ejections are usually associated with the occurrence of eruptive prominences, and that such mass ejections are generally accelerated through the corona, it is found that the observed variation of transient breadth with height contradicts some characteristics of self-similar fluid flow which have been ascribed to mass-ejection transients.

Macqueen, R. M.↗
DOE OSTI2024

On the Flow of a Cement Suspension: The Effects of Nano-Silica and Fly Ash Particles

Additives such as nano-silica and fly ash are widely used in cement and concrete materials to improve the rheology of fresh cement and concrete and the performance of hardened materials and increase the sustainability of the cement and concrete industry by reducing the usage of Portland cement. Therefore, it is important to study the effect of these additives on the rheological behavior of fresh cement. In this paper, we study the pulsating Poiseuille flow of fresh cement in a horizontal pipe by considering two different additives and when they are combined (nano-silica, fly ash, combined nano-silica, and fly ash). To model the fresh cement suspension, we used a modified form of the power-law model to demonstrate the dependency of the cement viscosity on the shear rate and volume fraction of cement and the additive particles. The convection–diffusion equation was used to solve for the volume fraction. After solving the equations in the dimensionless forms, we conducted a parametric study to analyze the effects of nano-silica, fly ash, and combined nano-silica and fly ash additives on the velocity and volume fraction profiles of the cement suspension. According to the parametric study presented here, larger nano-silica content results in lower centerline velocity of the cement suspension and larger non-uniformity of the volume fraction. Compared to nano-silica, fly ash exhibits an opposite effect on the velocity. Larger fly ash content results in higher centerline velocity, while the effect of the fly ash on the volume fraction is not obvious. For cement suspension containing combined nano-silica and fly ash additives, nano-silica plays a dominant role in the flow behavior of the suspension. The findings of the study can help the design and operation of the pulsating flow of fresh cement mortars and concrete in the 3D printing industry.

36 MATERIALS SCIENCE↗
DOE OSTI2020

A buoyancy–shear–drag-based turbulence model for Rayleigh–Taylor, reshocked Richtmyer–Meshkov, and Kelvin–Helmholtz mixing

A phenomenological turbulence model for Rayleigh–Taylor, reshocked Richtmyer–Meshkov, and Kelvin–Helmholtz instability-induced mixing is developed using a general buoyancy–shear–drag model. Analytic solutions to the simplified model equations corresponding to each instability are derived separately, which are then used to calibrate the model coefficients to predict specific values of the mixing layer growth parameters and exponents. The buoyancy-shear–drag equations for the bubble and spike mixing layer widths are then solved numerically, and a turbulent diffusivity (or viscosity) is constructed dimensionally as the product of the mixing layer width, h, and its time-derivative, dh/dt. Surrogate turbulent fields are then constructed by multiplying a presumed, approximate self-similar spatial profile by appropriate functions of h and dh/dt. Using several simplifying approximations, the turbulent diffusion equations satisfied by the mean mass fraction and mean shear velocity are solved analytically in a reference frame moving with the mean advection velocity, and it is shown that the mean fields evolve in space and time as expected. The explicit modeling and solution of turbulent transport equations (e.g., for the turbulent kinetic energy and its dissipation rate or a lengthscale) are not required. By using separate equations for the bubble and spike evolution, the model naturally captures the asymmetry induced by an increasing density contrast between the heavy and light fluids. The model also includes molecular dissipation and diffusion, and therefore, can describe transition to fully-developed turbulence if the initial turbulent diffusivity/viscosity is much smaller than the molecular diffusivity/viscosity. This model is applied to constant-acceleration Rayleigh–Taylor, impulsively reshocked Richtmyer–Meshkov, and variable-density Kelvin–Helmholtz mixing layers to demonstrate its utility. It is shown that the numerical solutions of the model calibrated using specific values of the instability growth parameters and exponents: (1) produces mixing layer widths in agreement with the expected self-similar growth power laws; (2) gives spatiotemporal profiles of turbulent fields that are expected and consistent with previous results; (3) predicts the expected power-law growths and decays of the spatially-integrated turbulent fields, and; (4) gives spatiotemporal profiles of the mean fields that are expected and consistent with previous results. This relatively simple model can be implemented numerically, and is expected to be useful for fundamental as well as astrophysical and inertial confinement fusion applications.

42 ENGINEERING↗
DOE OSTI2023

Symmetry and scaling in one-dimensional compressible two-phase flow

Investigations of shock compression of heterogeneous materials often focus on the shock front width and overall profile. The number of experiments required to fully characterize the dynamic response of a material often belie the structure–property relationships governing these aspects of a shock wave. Recent observations measured a pronounced shock-front width on the order of 10 s of ns in particulate composites. We focus on particulate composites with disparate densities and investigate whether the mechanical interactions between the phases are adequate to describe this emergent behavior. The analysis proceeds with a general Mie–Grüneisen equation of state for the matrix material, a general drag force law with general power-law scaling for the particle-matrix coupling of the phases, and a volume fraction-dependent viscosity. Lie group analysis is applied to one-dimensional hydrodynamic flow equations for the self-consistent interaction of particles embedded in a matrix material. The particle phase is characterized by a particle size and volume fraction. The Lie group analysis results in self-similar solutions reflecting the symmetries of the flow. The symmetries lead to well-defined scaling laws, which may be used to characterize the propagation of shock waves in particle composites. An example of the derived scaling laws for shock attenuation and rise time is shown for experimental data on shock-driven tungsten-loaded polymers. A key result of the Lie analysis is that there is a relationship between the exponents characterizing the form of the drag force and the exponent characterizing the shock velocity and its attenuation in a particulate composite. Comparison to recent experiments results in a single exponent that corresponds to a conventional drag force.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Turbulence, Waves, and Taylor’s Hypothesis for Heliosheath Observations

Magnetic field fluctuations measured in the heliosheath by the Voyager spacecraft are often characterized as compressible, as indicated by a strong fluctuating component parallel to the mean magnetic field. However, the interpretation of the turbulence data faces the caveat that the standard Taylor’s hypothesis is invalid because the solar wind flow velocity in the heliosheath becomes subsonic and slower than the fast magnetosonic speed, given the contributions from hot pickup ions (PUIs) in the heliosheath. We attempt to overcome this caveat by introducing a 4D frequency-wavenumber spectral modeling of turbulence, which is essentially a decomposition of different wave modes following their respective dispersion relations. Isotropic Alfvén and fast mode turbulence are considered to represent the heliosheath fluctuations. We also include two dispersive fast wave modes derived from a three-fluid theory. We find that (1) magnetic fluctuations in the inner heliosheath are less compressible than previously thought, an isotropic turbulence spectral model with about 25% in compressible fluctuation power is consistent with the observed magnetic compressibility in the heliosheath; (2) the hot PUI component and the relatively cold solar wind ions induce two dispersive fast magnetosonic wave branches in the perpendicular propagation limit, PUI fast wave may account for the spectral bump near the proton gyrofrequency in the observable spectrum; (3) it is possible that the turbulence wavenumber spectrum is not Kolmogorov-like although the observed frequency spectrum has a −5/3 power-law index, depending on the partitioning of power among the various wave modes, and this partitioning may change with wavenumber.

79 ASTRONOMY AND ASTROPHYSICS↗