Search NASA⌕ Search

SEARCH · Search NASA

Results for “viscous effects”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Viscous effects on plasmoid formation from nonlinear resistive tearing growth in a Harris sheet

Abstract In this work, the evolution of a highly unstable m = 1 resistive tearing mode, leading to plasmoid formation in a Harris sheet, is studied in the framework of full MHD model using the Non-Ideal Magnetohydrodynamics with Rotation, Open Discussion simulation. Following the initial nonlinear growth of the primary m = 1 island, the X-point develops into a secondary elongated current sheet that eventually breaks into plasmoids. Two distinctive viscous regimes are found for the plasmoid formation and saturation. In the low viscosity regime (i.e. P r ≲ 1), the plasmoid width increases sharply with viscosity, whereas in the viscosity dominant regime (i.e. P r ≳ 1), the plasmoid size gradually decreases with viscosity. Such a finding quantifies the role of viscosity in modulating the plasmoid formation process through its effects on the plasma flow and the reconnection itself.

Physics↗
DOE OSTI2022

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↗
DOE OSTI2021

Molecular Viscosity and Diffusivity Effects in Transitional and Shock-Driven Mixing Flows

Here, this paper investigates the importance of molecular viscosity and diffusivity for the prediction of transitional and shock-driven mixing flows featuring high and low Reynolds and Mach number regions. Two representative problems are computed with implicit large-eddy simulations using the inviscid Euler equations (EE) and viscous Navier-Stokes equations (NSE): the Taylor-Green vortex at Reynolds number Re = 3000 and initial Mach number Ma = 0.28, and an air- SF 6 -air gas curtain subjected to two shock waves at Ma = 1.2. The primary focus is on differences between NSE and EE predictions due to viscous effects. The outcome of the paper illustrates the advantages of utilizing NSE. In contrast to the EE, where the effective viscosity decreases upon grid refinement, NSE predictions can be assessed for simulations of flows with transition to turbulence at prescribed constant Re. The NSE can achieve better agreement between solutions and reference data, and the results converge upon grid refinement. On the other hand, the EE predictions do not converge with grid refinement, and can only exhibit similarities with the NSE results at coarse grid resolutions. We also investigate the effect of viscous effects on the dynamics of the coherent and turbulent fields, as well as on the mechanisms contributing to the production and diffusion of vorticity. The results show that nominally inviscid calculations can exhibit significantly varying flow dynamics driven by changing effective resolution-dependent Reynolds number, and highlight the role of viscous processes affecting the vorticity field. These tendencies become more pronounced upon grid refinement. The discussion of the results concludes with the assessment of the computational cost of inviscid and viscous computations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Two-particle transverse momentum correlations in p p and p -Pb collisions at energies available at the CERN Large Hadron Collider

Two-particle transverse momentum differential correlators, recently measured in Pb-Pb collisions at energies available at the CERN Large Hadron Collider (LHC), provide an additional tool to gain insights into particle production mechanisms and infer transport properties, such as the ratio of shear viscosity to entropy density, of the medium created in Pb-Pb collisions. The longitudinal long-range correlations and the large azimuthal anisotropy measured at low transverse momenta in small collision systems, namely pp and p-Pb, at LHC energies resemble manifestations of collective behavior. This suggests that locally equilibrated matter may be produced in these small collision systems, similar to what is observed in Pb-Pb collisions. In this work, the same two-particle transverse momentum differential correlators are exploited in pp and p-Pb collisions at $\sqrt{s}$ = 7 TeV and $\sqrt{s_{NN}}$ = 5.02 TeV, respectively, to seek evidence for viscous effects. Specifically, the strength and shape of the correlators are studied as a function of the produced particle multiplicity to identify evidence for longitudinal broadening that might reveal the presence of viscous effects in these smaller systems. The measured correlators and their evolution from pp and p-Pb to Pb-Pb collisions are additionally compared to predictions from Monte Carlo event generators, and the potential presence of viscous effects is discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Temperature Effects on Droplet Oscillation Decay with Application to Fuel Property Measurement

Observation of oscillation decay in droplets has been shown to be an effective approach in determining physical properties such as viscosity and surface tension for emerging biofuels using μl quantities. Herein this work extends the approach to higher temperatures relevant to fuel injection conditions for internal combustion engines. Experiments are conducted that use high-resolution strobed imaging of moving, heated, μm-sized, fuel droplets to capture shape oscillation decay through image processing and analysis. Two fuels are investigated, iso-butanol and a primary reference fuel (PRF 84), which is a mixture of iso-octane and n-heptane. A piezoelectric droplet generator is used to generate a continuous train of single droplets, which are given an initial perturbation and subsequently undergo damped oscillations, captured using strobed imaging. Surface tension and viscosity calculated at four different temperatures ranging from 30°C-55°C using the frequency and decay time associated with the fundamental mode are found to be within ~ 10% of reference values obtained from the literature. Complementary numerical simulations are performed that utilize a volume-of-fluid approach to track transient droplet oscillation phenomena along with heat and mass transfer in a 2D axisymmetric domain. Simulations, where droplet size and temperature can be independently varied, capture an expected decrease in oscillation frequency and increase in decay time, with increase in fuel temperature. The simulations are further used to study the relative contribution to deviations in surface tension and viscosity predictions due to heat and mass transfer from the droplet, as well as viscous effects violating the inviscid flow assumption in the droplet oscillation theory. Mass loss effects are found to be negligible. Temperature change due to heat transfer is found to have the next highest sensitivity, particularly for the more volatile PRF 84, which has a lower viscosity and associated Ohnesorge number. For isobutanol, which has a higher viscosity and Ohnesorge number, viscous effects contribute the most to deviation in fuel property predictions.

09 BIOMASS FUELS↗
DOE OSTI2020

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↗
DOE OSTI2023

Molecular-gas-dynamics simulations of turbulent Couette flow over a mean-free-path-scale permeable substrate

Here we report flow statistics and visualizations from molecular-gas-dynamics simulations using the direct simulation Monte Carlo (DSMC) method for turbulent Couette flow in a minimal domain where the lower wall is replaced by an idealized permeable fibrous substrate representative of thermal-protection-system materials for which the Knudsen number is O(10 -1 ). Comparisons are made with smooth-wall DSMC simulations and smooth-wall direct numerical simulations (DNS) of the Navier-Stokes equations for the same conditions. Roughness, permeability, and noncontinuum effects are assessed. In the range of Reynolds numbers considered herein, the scalings of the skin friction on the permeable substrate and of the mean flow within the substrate suggest that they are dominated by viscous effects. While the regenerative cycle characteristic of smooth-wall turbulence remains intact for all cases considered, we observe that the near-wall velocity fluctuations are modulated by the permeable substrate with a wavelength equal to the pore spacing. Additionally, the flow within the substrate shows significant rarefaction effects, resulting in an apparent permeability that is 13% larger than the intrinsic permeability. In contrast, the smooth-wall DSMC and DNS simulations exhibit remarkably good agreement for the statistics examined, despite the Knudsen number based on the viscous length scale being as large as O(10 -1 ). This latter result is at variance with classical estimates for the breakdown of the continuum assumption and calls for further investigations into the interaction of noncontinuum effects and turbulence.

42 ENGINEERING↗
DOE OSTI2021

Distinctive rheological and temporal viscoelastic behaviour of alkali-activated fly ash/slag pastes: A comparative study with cement paste

This paper presents the results of the investigation on the fundamental differences in the viscoelasticity between alkali-activated fly ash/slag materials (AAMs) and cement pastes. The effects of the precursor and the activator on the rheological behaviour of AAM pastes were studied. Given the specific precursor, the activator viscosity significantly affected the AAM paste viscosity. The high viscous activator increased the plastic viscosity of AAM pastes by 4–8 times higher than that of cement paste and drastically decreased the yield stress due to viscous effects and weak colloidal interactions. Temporal changes in viscoelasticity showed that the negligible colloidal interactions between particles in AAM paste made the system non-percolated until the initial setting. This resulted in very different viscoelastic behaviour compared to the one from the well-percolated network in cement paste. Considering all results obtained in this study, the paper describes the short-term evolution of the AAM paste from the fresh condition to the initial setting.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Far-from-equilibrium search for the QCD critical point

Initial conditions for relativistic heavy-ion collisions may be far from equilibrium (i.e., there are large initial contributions from the shear stress tensor and bulk pressure) but it is expected that on very short time scales the dynamics converge to a universal attractor that defines hydrodynamic behavior. Thus far, studies of this nature have only considered an idealized situation at LHC energies (high temperatures T and vanishing baryon chemical potential µ B = 0) but, in this work, we investigate for the first time how far-from-equilibrium effects may influence experimentally driven searches for the Quantum Chromodynamic critical point at RHIC. We find that the path to the critical point is heavily influenced by far from equilibrium initial conditions where viscous effects lead to dramatically different {T, µ B } trajectories through the QCD phase diagram. We compare hydrodynamic equations of motion with shear and bulk coupled together at finite µB for both DNMR and phenomenological Israel-Stewart equations of motion and discuss their influence on potential attractors at finite µ B and their corresponding {T, µ B } trajectories.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Effects of electron viscosity on resonant layer responses to non-axisymmetric magnetic perturbations

The resonant field penetration to magnetic islands is the central MHD mechanism of non-axisymmetric plasma responses in a tokamak such as disruptive locking or favorable ELM stabilization. The resonant field penetration can be induced by any non-ideal processes as manifested in the delicate balance under the generalized Ohm's law. Here, we show that the viscous effects by electrons are not ignorable in the field penetration unlike previous presumption, even if the electron viscosity is as small as the square root of its mass compared to the ions. It is clear that its effects become only bigger if the electron viscosity becomes anomalously large. The work strictly follows the three-field model in the linear regime targeting the prediction of the onset of the field penetration and successfully extending it with electron viscosity and identifying new regimes. The results also indicate that the error field thresholds become more strongly dependent on plasma density than ones predicted in the linear regimes without the electron viscosity, which is consistent with experimental observations and thus a significant implication.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

The Viscoplastic Behavior of Natural Hydrate‐Bearing Sandy‐Silts Under Uniaxial Strain Compression ( K 0 Loading)

Abstract The in‐situ stress state and geomechanical properties of hydrate‐bearing sediments impact hydrate formation and gas production strategies. We explore the uniaxial strain compression and stress evolution of natural hydrate‐bearing sandy‐silts from Green Canyon Block 955 in the deep‐water Gulf of Mexico. We performed constant rate of strain uniaxial strain experiments, interrupted by periods where we held the axial stress constant, to explore the vertical deformation and the evolution of the ratio of lateral to axial effective stress ( K 0 ) with time. The hydrate‐bearing sandy‐silt is stiffer and has a larger K 0 than the equivalent hydrate‐free sediment upon loading. During stress holds, the void ratio decreases sigmoidally with the log of time, and K 0 converges to isotropic conditions. We interpret that during loading, the hydrate bears the load and deforms. With time, the hydrate redistributes the load and K 0 increases. We used a viscoelastic model to describe the behavior. The model accurately captures deformation and K 0 trends but does not reproduce all the complex interactions of the hydrate with the porous skeleton. We anticipate that viscous effects within hydrate sediments will impact reservoir compression and stresses during production (hours to days), result in isotropic stress state over geological timescales, and explain the creeping movement in submarine landslides.

58 GEOSCIENCES↗
DOE OSTI2022

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↗
DOE OSTI2024

An extension of Thwaites’ method for turbulent boundary layers

Thwaites developed an approximate method for determining the evolution of laminar boundary layers. The approximation follows from an assumption that the growth of a laminar boundary layer in the presence of pressure gradients could be parameterized solely as a function of the Holstein–Bohlen flow parameter, thus reducing the von Kármán momentum integral to a first-order ordinary differential equation. This method is useful for the analysis of laminar flows, and in computational potential flow solvers to account for the viscous effects. In this work, an approximate method for determining the momentum thickness of a two-dimensional, turbulent boundary layer is proposed following Thwaites’ work. It is shown that the method provides good estimates of the momentum thickness for multiple boundary layers, including both favourable and adverse pressure gradient effects, up to the point of separation. In the limit of high Reynolds numbers, it is possible to derive a criterion for the onset of separation from the proposed model, which is shown to be in agreement with prior empirical observations. The sensitivity of the separation location with respect to upstream perturbations is also analysed through this model for the NASA/Boeing speed bump and the transonic Bachalo–Johnson bump.

17 WIND ENERGY↗
DOE OSTI2023

Dynamics of multilayer Rayleigh–Taylor instability at moderately high Atwood numbers

This paper investigates the multilayer Rayleigh–Taylor instability (RTI) using statistically stationary experiments conducted in a gas tunnel. Employing diagnostics such as particle image velocimetry (PIV) and planar laser induced fluorescence (PLIF), we make simultaneous velocity–density measurements to study how dynamics and mixing are linked in this variable density flow. Experiments are conducted in a newly built, blow-down three-layer gas tunnel facility. Mixing between three gas streams is studied, where the top and bottom streams are comprised of air, and the middle stream is an air–helium mixture. Shear is minimized between these streams by matching their inlet velocities. The four experimental conditions investigated here consist of two different density ratios (Atwood numbers 0.3 and 0.6), each investigated at two instability development times (or equivalently, two streamwise locations), and all experiments are with the same middle stream thickness of 3 cm. The growth of the middle layer is measured using laser-based planar Mie scattering visualization. The mixing width is found to grow linearly with time at late times. Various quantitative measures of molecular mixing indicate a very high degree of molecular mixing at late times in the multilayer RTI flow. The vertical turbulent mass flux $a_y$ is calculated. In addition to mostly negative values of $a_y$ , typical of buoyancy-dominated flows due to negative correlation between velocity and density fluctuations, positive regions are also observed in profiles of $a_y$ due to entrainment and erosion at the lower edge of the mixing region. Global energy budgets are calculated for the multilayer RTI flow at late times and it is found that the majority of potential energy released has been dissipated due to viscous effects, and a large value of mixing efficiency ( $\sim$ 60 %) is observed.

Mechanics↗
DOE OSTI2025

Assessment of integral models for non-Boussinesq lazy plumes using numerical simulations

Integral modelling of turbulent buoyant plumes is crucial for rapid predictions of plume characteristics. While the governing equations are typically derived using self-similarity and a Boussinesq approximation, these assumptions may not hold for plumes originating from finite-area sources with large density ratios. Here, this work evaluates the accuracy of integral-scale models for non-Boussinesq lazy plumes using high-fidelity numerical simulations of turbulent helium plumes. We analyse the plume kinematics by computing vertical fluxes, plume radius and radial profiles, establishing some disparities between common practice and physical accuracy. We identify how the definition of the plume radius changes the perception of the plume structure when the flow is not self-similar and derive a relationship between the flux-based and threshold-based definitions without requiring self-similarity. We then examine the plume dynamics by evaluating the source terms from the governing plume equations. Our results support neglecting diffusive and viscous effects but emphasise the importance of the mean pressure gradient, even in the self-similar regime. Two coefficients need to be modelled: the well-known entrainment coefficient and the lesser-known momentum correction coefficient, which is a correction required for the momentum equation to account for self-similar and slender approximations. The momentum correction coefficient is found to be approximately constant and slightly greater than the assumed value of 1. The standard entrainment coefficient models perform well up to a local Richardson number three times the asymptotic value but overpredict entrainment for larger Richardson numbers. We propose a correction using the known finite limit of entrainment at infinite Richardson number.

Meehan, Michael Alexander [Sandia National Laborat↗
DOE OSTI2023

Non-ideal MHD growth of current interchange tearing modes at plasma edge and response to externally-imposed flow

Current interchange tearing modes (CITMs) have been proposed to explain an intermittent eruption of electric current toward a divertor at an edge region of a tokamak device. Extended magnetohydrodynamic (MHD) simulations combined with the Scrape-Off-Layer diffusion model, which has been developed by the authors [Miura et al., Phys. Plasmas 24, 092111 (2017)], are carried out in order to study the growth of a CITM under influences of two-fluid and gyro-viscous effects, or non-ideal MHD effects. Numerical simulations both with and without an externally imposed azimuthal flow show that the growth of a CITM is observed for a relatively small flow velocity, and that a CITM can be a candidate mechanism for the intermittent eruption of the current in a tokamak, whether the growth is under influences of non-ideal MHD effects or not. Here, the growth can be suppressed when an externally imposed azimuthal flow is sufficiently large to cause a finite radial displacement because such a displacement prevents a transition from an interchange to a tearing mode that is the key process of a CITM. Furthermore, it is also found that a stripe pattern similar to so-called streamers is formed in the course of the CITM growth. The basic nature of this streamer-like structure is also presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Validation of predictive performance models for supersonic gas-jet nozzles at the Laboratory for Laser Energetics

We present results characterizing the neutral-density distributions produced by the supersonic nozzles used in experiments on the OMEGA-60 and OMEGA-EP laser systems at the University of Rochester’s Laboratory for Laser Energetics (LLE). Axisymmetric Fluent® simulations using LLE nozzle specifications capture the viscous effects, gas expansion, and shock waves that complicate flow predictions for offsets above the nozzle exit. Here, these simulations show good agreement with neutral-density measurements obtained using a four-wave shearing interferometer. An analytical form is given for the plateau length. Fits to simulation data for boundary layer thickness, mean plateau density, and density ramps are given as functions of nozzle offset and nozzle backing pressure for a number of nozzles and gases.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

Rippled shock propagation in a laser-driven target at multimegabar pressures

The evolution of non-uniform shocks produced by modulated laser irradiation or surface perturbations is relevant to studies of inertial confinement fusion and material properties at high-energy-density conditions. We present results from an experiment conducted at the OMEGA EP laser facility, where a 300 GPa shock was driven into a fused silica sample with pre-fabricated single-mode surface modulations. Using time-resolved optical velocimetry, we captured the continuous evolution of rippled shock motion, enabling a comprehensive mapping of the spatial amplitude history from formation to phase reversal in a single experiment. Initially, the ablation-driven shock inherits a fraction of the surface modulation amplitude from the sample, which subsequently grows before decaying, ultimately leading to the flattening of the rippled shock and a phase reversal. We find that two-dimensional inviscid hydrodynamic simulation of the experiment is able to qualitatively capture many aspects of the rippled shock evolution but over-predicts the initial amplitude growth. This experimental platform, capable of accommodating varying ripple wavelengths, lays the groundwork for a potential viscometry method at extreme pressures, where viscous effects manifest as differences in shock flattening times between rippled shocks of two distinct wavelengths propagating through the sample.

36 MATERIALS SCIENCE↗