Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “VISCOSITY”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Accretion of Jupiter-mass planets in the limit of vanishing viscosity

In the core-accretion model, the nominal runaway gas-accretion phase brings most planets to multiple Jupiter masses. However, known giant planets are predominantly Jupiter mass bodies. Obtaining longer timescales for gas accretion may require using realistic equations of states, or accounting for the dynamics of the circumplanetary disk (CPD) in the low-viscosity regime, or both. Here we explore the second way by using global, three-dimensional isothermal hydrodynamical simulations with eight levels of nested grids around the planet. In our simulations, the vertical inflow from the circumstellar disk (CSD) to the CPD determines the shape of the CPD and its accretion rate. Even without a prescribed viscosity, Jupiter's mass-doubling time is ∼10{sup 4} yr, assuming the planet at 5.2 AU and a Minimum Mass Solar Nebula. However, we show that this high accretion rate is due to resolution-dependent numerical viscosity. Furthermore, we consider the scenario of a layered CSD, viscous only in its surface layer, and an inviscid CPD. We identify two planet-accretion mechanisms that are independent of the viscosity in the CPD: (1) the polar inflow—defined as a part of the vertical inflow with a centrifugal radius smaller than two Jupiter radii and (2) the torque exerted by the star on the CPD. In the limit of zero effective viscosity, these two mechanisms would produce an accretion rate 40 times smaller than in the simulation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Projecting the likely importance of weak-interaction-driven bulk viscosity in neutron star mergers

ABSTRACT In this work, we estimate how much bulk viscosity driven by Urca processes is likely to affect the gravitational wave signal of a neutron star coalescence. In the late inspiral, we show that bulk viscosity affects the binding energy at fourth post-Newtonian order. Even though this effect is enhanced by the square of the gravitational compactness, the coefficient of bulk viscosity is likely too small to lead to observable effects in the waveform during the late inspiral, when only considering the orbital motion itself. In the post-merger, however, the characteristic time-scales and spatial scales are different, potentially leading to the opposite conclusion. We post-process data from a state-of-the-art equal-mass binary neutron star merger simulation to estimate the effects of bulk viscosity (which was not included in the simulation itself). In that scenario, we find that bulk viscosity can reach high values in regions of the merger. We compute several estimates of how much it might directly affect the global dynamics of the considered merger scenario, and find that it could become significant. Even larger effects could arise in different merger scenarios or in simulations that include non-linear effects. This assessment is reinforced by a quantitative comparison with relativistic heavy-ion collisions where such effects have been explored extensively.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

Influence of resonant charge exchange on the viscosity of partially ionized plasma in a magnetic field

The influence of resonant charge exchange for ion-atom interaction on the viscosity of partially ionized plasma embedded in the magnetic field is investigated. The general system of equations used to derive the viscosity coefficients for an arbitrary plasma component in the 21-moment approximation of Grad’s method is presented. The expressions for the coefficients of total and partial viscosities of a multicomponent partially ionized plasma in the magnetic field are obtained. As an example, the coefficients of the parallel and transverse viscosities for the ionic and neutral components of the partially ionized hydrogen plasma are calculated. It is shown that the account for resonant charge exchange can lead to a substantial change of the parallel and transverse viscosity of the plasma components in the region of low degrees of ionization on the order of 0.1.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

The Viscosity and Atomic Structure of Volatile-Bearing Melilititic Melts at High Pressure and Temperature and the Transport of Deep Carbon

Understanding the viscosity of mantle-derived magmas is needed to model their migration mechanisms and ascent rate from the source rock to the surface. High pressure–temperature experimental data are now available on the viscosity of synthetic melts, pure carbonatitic to carbonate–silicate compositions, anhydrous basalts, dacites and rhyolites. However, the viscosity of volatile-bearing melilititic melts, among the most plausible carriers of deep carbon, has not been investigated. In this study, we experimentally determined the viscosity of synthetic liquids with ~31 and ~39 wt% SiO 2 , 1.60 and 1.42 wt% CO 2 and 5.7 and 1 wt% H 2 O, respectively, at pressures from 1 to 4.7 GPa and temperatures between 1265 and 1755 °C, using the falling-sphere technique combined with in situ X-ray radiography. Our results show viscosities between 0.1044 and 2.1221 Pa·s, with a clear dependence on temperature and SiO 2 content. The atomic structure of both melt compositions was also determined at high pressure and temperature, using in situ multi-angle energy-dispersive X-ray diffraction supported by ex situ microFTIR and microRaman spectroscopic measurements. Our results yield evidence that the T–T and T–O (T = Si,Al) interatomic distances of ultrabasic melts are higher than those for basaltic melts known from similar recent studies. Based on our experimental data, melilititic melts are expected to migrate at a rate ~from 2 to 57 km·yr -1 in the present-day or the Archaean mantle, respectively.

36 MATERIALS SCIENCE↗
DOE OSTI2013

Entropy viscosity method applied to Euler equations

The entropy viscosity method [4] has been successfully applied to hyperbolic systems of equations such as Burgers equation and Euler equations. The method consists in adding dissipative terms to the governing equations, where a viscosity coefficient modulates the amount of dissipation. The entropy viscosity method has been applied to the 1-D Euler equations with variable area using a continuous finite element discretization in the MOOSE framework and our results show that it has the ability to efficiently smooth out oscillations and accurately resolve shocks. Two equations of state are considered: Ideal Gas and Stiffened Gas Equations Of State. Results are provided for a second-order time implicit schemes (BDF2). Some typical Riemann problems are run with the entropy viscosity method to demonstrate some of its features. Then, a 1-D convergent-divergent nozzle is considered with open boundary conditions. The correct steady-state is reached for the liquid and gas phases with a time implicit scheme. The entropy viscosity method correctly behaves in every problem run. For each test problem, results are shown for both equations of state considered here. (authors)

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

Artificial viscosity: back to the basics

Summary In this paper, we take a different perspective on the derivation of artificial viscosity. Heretofore, the development of artificial viscosity has been based on the paper published in Journal of Applied Physics in 1950 authored by John von Neumann and Robert Richtmyer [1]. Earlier, in 1948, Richtmyer published a report at Los Alamos Scientific Laboratory documenting the original concept [2]. This report was the true origin of shock capturing methods and contains several key ideas that are conceptually different than the 1950 journal article. Unfortunately, this report (LA‐671) was classified until 1993. This has resulted in two issues: the misattribution of the invention of artificial viscosity as primarily being the work of von Neumann and the loss of the structurally different ideas in the original report. We seek to right the record of history here and use the ideas contained in Richtmyer's report to good effect in deriving a new shock viscosity. The focus of previous development has been the Hugoniot curve describing the locus of states connected by a single shock wave. Here we follow a path more focused upon the Rayleigh line, which is strongly guided by Richtmyer's line of development of the original artificial viscosity formulation. We provide an implementation of the method resulting from this perspective and computational results for simple shock problems. Copyright © 2014 John Wiley & Sons, Ltd.

Mattsson, Ann E.↗
DOE OSTI2013

Volume-translated cubic EoS and PC-SAFT density models and a free volume-based viscosity model for hydrocarbons at extreme temperature and pressure conditions

This research focuses on providing the petroleum reservoir engineering community with robust models of hydrocarbon density and viscosity at the extreme temperature and pressure conditions (up to 533 K and 276 MPa, respectively) characteristic of ultra-deep reservoirs, such as those associated with the deepwater wells in the Gulf of Mexico. Our strategy is to base the volume-translated (VT) Peng–Robinson (PR) and Soave–Redlich–Kwong (SRK) cubic equations of state (EoSs) and perturbed-chain, statistical associating fluid theory (PC-SAFT) on an extensive data base of high temperature (278–533 K), high pressure (6.9–276 MPa) density rather than fitting the models to low pressure saturated liquid density data. This high-temperature, high-pressure (HTHP) data base consists of literature data for hydrocarbons ranging from methane to C{sub 40}. The three new models developed in this work, HTHP VT-PR EoS, HTHP VT-SRK EoS, and hybrid PC-SAFT, yield mean absolute percent deviation values (MAPD) for HTHP hydrocarbon density of ~2.0%, ~1.5%, and <1.0%, respectively. An effort was also made to provide accurate hydrocarbon viscosity models based on literature data. Viscosity values are estimated with the frictional theory (f-theory) and free volume (FV) theory of viscosity. The best results were obtained when the PC-SAFT equation was used to obtain both the attractive and repulsive pressure inputs to f-theory, and the density input to FV theory. Both viscosity models provide accurate results at pressures to 100 MPa but experimental and model results can deviate by more than 25% at pressures above 200 MPa.

02 PETROLEUM↗
DOE OSTI2020

Viscosity-temperature relation based on the evolution of medium-range structures of silica

The viscosity of glass is the most important technological property to glass manufactures and various applications. Practically, finding an accurate equation to express the glass viscosity behavior in the entire glass transition temperature range is tremendously challenge because it spans more than ten orders of magnitude. After a brief review of existing empirical viscosity equations, this work focuses on the correlating silica viscosity behavior with the evolution of glass medium-range structure, based on the recently proposed nanoflake model. From this new model, a new equation is constructed, which correctly describes the Arrhenius-type behavior of silica viscosity η above the melting temperature Tm, and non-Arrhenius-type behavior from Tm to a critical temperature Tc. At temperature lower than Tc, the equation predicts Arrhenius-type behavior again for the viscous flow with increased activation energy. The new equation agrees with experimental data in the entire temperature spanned from extremely high to extremely low. Furthermore, the application of the new equation is shown to extend to sodium silicate glasses as well.

36 MATERIALS SCIENCE↗
DOE OSTI2015

A study of overflow simulations using MPAS-Ocean: Vertical grids, resolution, and viscosity

MPAS-Ocean is used to simulate an idealized, density-driven overflow using the dynamics of overflow mixing and entrainment (DOME) setup. Additionally, numerical simulations are carried out using three of the vertical coordinate types available in MPAS-Ocean, including z-star with partial bottom cells, z-star with full cells, and sigma coordinates. The results are first benchmarked against other models, including the MITgcm’s z-coordinate model and HIM’s isopycnal coordinate model, which are used to set the base case used for this work. A full parameter study is presented that looks at how sensitive overflow simulations are to vertical grid type, resolution, and viscosity. Horizontal resolutions with 50 km grid cells are under-resolved and produce poor results, regardless of other parameter settings. Vertical grids ranging in thickness from 15 m to 120 m were tested. A horizontal resolution of 10 km and a vertical resolution of 60 m are sufficient to resolve the mesoscale dynamics of the DOME configuration, which mimics real-world overflow parameters. Mixing and final buoyancy are least sensitive to horizontal viscosity, but strongly sensitive to vertical viscosity. This suggests that vertical viscosity could be adjusted in overflow water formation regions to influence mixing and product water characteristics. In conclusion, the study shows that sigma coordinates produce much less mixing than z-type coordinates, resulting in heavier plumes that go further down slope. Sigma coordinates are less sensitive to changes in resolution but as sensitive to vertical viscosity compared to z-coordinates.

58 GEOSCIENCES↗
DOE OSTI2026

Revisiting the validity of eddy viscosity models for predicting airflow over water waves

In this study, we revisit the validity of eddy viscosity models for predicting wave-induced airflow disturbances over ocean surface waves. We first derive a turbulence curvilinear model for the phase-averaged Navier–Stokes equations, extending the work of Cao, Deng & Shen (2020 J. Fluid Mech. 901, A27), by incorporating turbulence stress terms previously neglected in the linearised viscous curvilinear model. To verify our formulation, we perform a priori tests by numerically solving the model using mean wind and turbulence stress profiles from large-eddy simulations (LES) of airflow over waves across various wave ages. Results show that including turbulence stress terms improves wave-induced airflow predictions compared with the previous viscous curvilinear model. We further show that using a standard mixing-length eddy viscosity yields inaccurate predictions at certain wave ages, as it fails to capture wave-induced turbulence, which fundamentally differs from mean shear-driven turbulence. The LES data show that accurate representations of wave-induced stresses require a complex-valued eddy viscosity. The maximum magnitude of this eddy viscosity scales as ∼𝑢 𝜏 ⁢𝜁 𝑖𝑛𝑛𝑒𝑟 , where 𝑢 𝜏 is the friction velocity and 𝜁 𝑖𝑛𝑛𝑒𝑟 is the inner-layer thickness, the height at which the eddy-turnover time matches the wave advection time scale. This scaling aligns with the prediction by Belcher & Hunt (1993 J. Fluid Mech. 251, 109–148). Overall, the findings demonstrate that traditional eddy viscosity models are inadequate for capturing wave-induced turbulence. More sophisticated turbulence models are essential for the accurate prediction of airflow disturbances and form drag in wind–wave interaction models.

16 TIDAL AND WAVE POWER↗
DOE OSTI2018

Effects of surface tension and viscosity on gold and silver sputtered onto liquid substrates

Here, we report on DC magnetron sputtering of gold and silver onto liquid substrates of varying viscosities and surface tensions. We were able to separate the effects of viscosity from surface tension by depositing the metals onto silicone oils with a range of viscosities. The effects of surface tension were studied by depositing the metals onto squalene, poly(ethylene glycol), and glycerol. It was discovered that dispersed nanoparticles were formed on liquids with low surface tension and low viscosity whereas dense films were formed on liquids with low surface tension and high viscosity. Nanoparticles were formed on both the liquid surface and within the bulk liquid for high surface tension liquids. Our findings can be used to tailor the metal and liquid interaction to fabricate particles and films for various applications in optics, electronics, and catalysis.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Numerical minimization of neoclassical poloidal viscosity for supersonic equilibria in tokamak geometry

Extensive experimental evidence has shown that the presence of poloidal flow in tokamaks can dramatically improve transport properties. However, theory indicates that poloidal flows are damped by poloidal viscosity, thus necessitating external drivers, such as neutral beam injection or radio frequency heating. In this work, ideal magnetohydrodynamic equilibria are calculated via the FORTRAN code FLOW [Guazzotto et al., Phys. Plasmas 11, 604 (2004)] and a postprocessor is used to estimate the neoclassical poloidal viscosity. The equilibrium inputs, which correspond to intuitive physical quantities, are then numerically optimized to reduce a viscosity figure of merit. We present supersonic equilibria in tokamak geometry with minimized neoclassical poloidal viscosities for various velocity free function inputs, plasma aspect ratios, and collisionality regimes. Benchmarks are made against an analytic theory as well as a classical expression of poloidal viscosity. Numerical confirmation of the analytic theory is obtained in the high aspect ratio and high collisionality limit. Good agreement is also seen near the plasma core and edge, with discrepancies arising in the intermediate region. Outside of these limits, rotation input function profiles are found that provide ∼order of magnitude improvements over the analytic theory, with additional progress being made toward predictions for tokamak-relevant equilibria.

Physics↗
DOE OSTI2025

Taking control of compressible modes: bulk viscosity and the turbulent dynamo

Many polyatomic astrophysical plasmas are compressible and out of chemical and thermal equilibrium, introducing a bulk viscosity into the plasma via the internal degrees of freedom of the molecular composition, directly impacting the decay of compressible modes, $\mathrm{{\boldsymbol {\mathit {v}}}}_{\parallel }(\boldsymbol {k})$. This is especially important for small-scale, turbulent dynamo processes in the interstellar medium (ISM), which are known to be sensitive to the effects of compression. To control the viscous properties of $\mathrm{{\boldsymbol {\mathit {v}}}}_{\parallel }(\boldsymbol {k})$, we perform trans-sonic, visco-resistive dynamo simulations with additional bulk viscosity $\nu _{\text{bulk}}$, deriving a new $\nu _{\text{bulk}}$ Reynolds number $\text{Re}_{\text{bulk}}$, and viscous Prandtl number $\text{P}\nu \equiv \text{Re}_{\text{bulk}}/ \text{Re}_{\text{shear}}$, where $\text{Re}_{\text{shear}}$ is the shear viscosity Reynolds number. We derive a framework for decomposing $E_{\rm mag}$ growth rates into incompressible and compressible terms via orthogonal tensor decompositions of $\boldsymbol {\nabla }\otimes \mathrm{{\boldsymbol {\mathit {v}}}}$, where $\mathrm{{\boldsymbol {\mathit {v}}}}$ is the fluid velocity. We find that $\mathrm{{\boldsymbol {\mathit {v}}}}_{\parallel }(\boldsymbol {k})$ play a dual role, growing and decaying $E_{\rm mag}$, and that field-line stretching is the main driver of growth, even in compressible dynamos. In the absence of $\nu _{\text{bulk}}$ ($\text{P}\nu \rightarrow \infty$), $\mathrm{{\boldsymbol {\mathit {v}}}}_{\parallel }(\boldsymbol {k})$ pile up on small-scales, creating a spectral bottleneck, which disappears for $\text{P}\nu \approx 1$. As $\text{P}\nu$ decreases, $\mathrm{{\boldsymbol {\mathit {v}}}}_{\parallel }(\boldsymbol {k})$ are dissipated at increasingly larger scales, in turn suppressing incompressible modes through a coupling between high-k modes. We emphasize the importance of further understanding the role of $\nu _{\text{bulk}}$ in compressible astrophysical plasmas, which we estimate could be as strong as the shear viscosity in the cold ISM, and highlight that compressible direct numerical simulations without bulk viscosity have unresolved compressible mode dissipation scales.

MHD↗
DOE OSTI2025

Shear viscosity from perturbative quantum chromodynamics to the hadron resonance gas at finite baryon, strangeness, and electric charge densities

Through model-to-data comparisons from heavy-ion collisions, it has been shown that the quark gluon plasma has an extremely small shear viscosity at vanishing densities. At large baryon densities, significantly less is known about the nature of the shear viscosity from quantum chromodynamics (QCD). Within heavy-ion collisions, there are three conserved charges: baryon number (B), strangeness (S), and electric charge (Q). Here we calculate the shear viscosity in two limits using perturbative QCD (pQCD) and an excluded-volume hadron resonance gas at finite BSQ densities. We then develop a framework that interpolates between these two limits such that shear viscosity is possible to calculate across a wide range of finite BSQ densities. We find that the pQCD and hadron resonance gas calculations have different BSQ density dependencies such that a rather nontrivial shear viscosity appears at finite densities.

Phenomenology↗
DOE OSTI2017

Phenomenological consequences of enhanced bulk viscosity near the QCD critical point

In the proximity of the QCD critical point the bulk viscosity of quark-gluon matter is expected to be proportional to nearly the third power of the critical correlation length, and become significantly enhanced. Here, this work is the first attempt to study the phenomenological consequences of enhanced bulk viscosity near the QCD critical point. For this purpose, we implement the expected critical behavior of the bulk viscosity within a non-boost-invariant, longitudinally expanding 1 + 1 dimensional causal relativistic hydrodynamical evolution at nonzero baryon density. We demonstrate that the critically enhanced bulk viscosity induces a substantial nonequilibrium pressure, effectively softening the equation of state, and leads to sizable effects in the flow velocity and single-particle distributions at the freeze-out. In conclusion, the observable effects that may arise due to the enhanced bulk viscosity in the vicinity of the QCD critical point can be used as complementary information to facilitate searches for the QCD critical point.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Bulk viscosity from Urca processes: npeμ matter in the neutrino-trapped regime

In this work, we extend our previous study of the bulk viscosity of hot and dense npe matter induced by the Urca process in the neutrino trapped regime to npeμ matter by adding the muonic Urca processes as well as the purely leptonic electroweak processes involving electron-muon transition. Here, the nuclear matter is modeled in a relativistic density functional approach with two different parametrizations which predict neutrino dominated matter (DDME2 model) and anti-neutrino dominated matter (NL3 model) at temperatures for which neutrinos/anti-neutrinos are trapped. In the case of neutrino-dominated matter, the main equilibration mechanism is lepton capture, whereas in the case of antineutrino-dominated matter this is due to neutron decay. We find that the equilibration rates of Urca processes are higher than that of the pure leptonic processes, which implies that the Urca bulk viscosity can be computed with the leptonic reactions assumed to be frozen. We find that the bulk viscosity decreases with temperature as ζ~T –2 at moderate temperatures. At high temperatures this scaling breaks down by sharp drops of the bulk viscosity close to the temperature where the proton fraction is density-independent and the matter becomes scale-invariant. This occurs also when the matter undergoes a transition from the antineutrino-dominated regime to the neutrino-dominated regime where the bulk viscosity attains a local maximum. We also estimate the bulk viscous dissipation timescales and find that these are in the range ≳ 1 s for temperatures above the neutrino trapping temperature. These timescales would be relevant only for long-lived objects formed in binary neutron star mergers and hot proto-neutron stars formed in core-collapse supernovas.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2018

Effects of particle size and briquetting of soda-lime-silicate glass batch on viscosity during batch-to-melt conversion

To assess the effects of grain size and briquetting on glass-batch melting, we measured the viscosity (?b) of an industrial glass batch during its conversion to molten simulated architectural glass within the temperature interval of 1050 °C to 1300 °C and at a shear rate of 0.02 s-1. On the basis of the mass fraction of undissolved sand (xs), which was determined using X-ray diffraction analysis on quenched batch samples, and the gas-phase fraction (p), which was obtained from volumetric measurements, we established the relationship ln(?b/?m) = a0 + bsxs + app, where ?m is the transition melt viscosity, a0 is the melt homogeneity coefficient, bs is the coefficient for undissolved sand, and ap is the coefficient for bubbles. The obtained values, bs = 8.756 and ap = 0.658, indicate that undissolved silica sand strongly influences batch viscosity. Consequently, compared with coarse-grained batches, fine raw materials and briquetting resulted in lower batch viscosity during conversion. Owing to their lower viscosities and higher heat conductivities, batches from finer raw materials or batches shaped into briquettes are expected to enhance the melting rate of batch blankets in industrial glass-melting furnaces.

Doi, Yoji↗
DOE OSTI2022

Verification study of xRAGE’s multi-ion viscosity model

This report on the adequacy of the plasma viscosity model implemented in LANL’s ASC program code, xRAGE is organized as follows. Section I provides a brief introduction on the kinetic theory which underpins plasma viscosity, and Section II describes how plasma kinetic theory can be used to calculate plasma viscosity coefficients accurate within the hydrodynamic limit. Section III describes the existing model for calculating the plasma viscosity coefficient in xRAGE. Next, Section IV shows tests of the xRAGE viscosity model against established hydrodynamic theory in the case of binary and trinary plasma mixtures.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗