DOE Office of Scientific and Technical Information (OSTI.GOV)

Matyas, Josef; Cooley, Scott K.; Sundaram, S. K.

Slags of low viscosity readily penetrate the refractory lining in slagging gasifiers, causing rapid and severe corrosion called spalling. In addition, a low-viscosity slag that flows down the gasifier wall forms a relatively thin layer of slag on the refractory surface, allowing the corrosive gases in the gasifier to participate in the chemical reactions between the refractory and the slag. In contrast, a slag viscosity of <25 Pa•s at 1400°C is necessary to minimize the possibility of plugging the slag tap. There is a need to predict and optimize slag viscosity so slagging gasifiers can operate continuously at temperatures rangingmore » from 1300 to 1650°C. The approach adopted in this work was to statistically design and prepare simulated slags, measure the viscosity as a function of temperature, and develop a model to predict slag viscosity based on slag composition and temperature. Statistical design software was used to select compositions from a candidate set of all possible vertices that will optimally represent the composition space for 10 main components. A total of 21 slag compositions were generated, including 5 actual coal slag compositions. The Arrhenius equation was applied to measured viscosity versus temperature data of tested slags, and the Arrhenius coefficients (A and B in ln(vis) = A + B/T) were expressed as linear functions of the slag composition. The viscosity model was validated using 1) data splitting approach, and 2) viscosity/temperature data of selected slag compositions from the literature that were formulated and melted at Pacific Northwest National Laboratory. The capability of the model to predict the viscosity of coal slags was compared with the model developed by Browning et al. because this model can predict the viscosity of slags from coal ash better than the most commonly used empirical models found in the literature.« less

Structure–property reduced order model for viscosity prediction in single-component CO 2 -binding organic liquids

DOE Office of Scientific and Technical Information (OSTI.GOV)

Cantu, David C.; Malhotra, Deepika; Koech, Phillip K.

2016-01-01

CO2 capture from power generation with aqueous solvents remains energy intensive due to the high water content of the current technology, or the high viscosity of non-aqueous alternatives. Quantitative reduced models, connecting molecular structure to bulk properties, are key for developing structure-property relationships that enable molecular design. In this work, we describe such a model that quantitatively predicts viscosities of CO2 binding organic liquids (CO2BOLs) based solely on molecular structure and the amount of bound CO2. The functional form of the model correlates the viscosity with the CO2 loading and an electrostatic term describing the charge distribution between the CO2-bearingmore » functional group and the proton-receiving amine. Molecular simulations identify the proton shuttle between these groups within the same molecule to be the critical indicator of low viscosity. The model, developed to allow for quick screening of solvent libraries, paves the way towards the rational design of low viscosity non-aqueous solvent systems for post-combustion CO2 capture. Following these theoretical recommendations, synthetic efforts of promising candidates and viscosity measurement provide experimental validation and verification.« less

Nonequilibrium viscosity of glass

NASA Astrophysics Data System (ADS)

Mauro, John C.; Allan, Douglas C.; Potuzak, Marcel

2009-09-01

Since glass is a nonequilibrium material, its properties depend on both composition and thermal history. While most prior studies have focused on equilibrium liquid viscosity, an accurate description of nonequilibrium viscosity is essential for understanding the low temperature dynamics of glass. Departure from equilibrium occurs as a glass-forming system is cooled through the glass transition range. The glass transition involves a continuous breakdown of ergodicity as the system gradually becomes trapped in a subset of the available configurational phase space. At very low temperatures a glass is perfectly nonergodic (or “isostructural”), and the viscosity is described well by an Arrhenius form. However, the behavior of viscosity during the glass transition range itself is not yet understood. In this paper, we address the problem of glass viscosity using the enthalpy landscape model of Mauro and Loucks [Phys. Rev. B 76, 174202 (2007)] for selenium, an elemental glass former. To study a wide range of thermal histories, we compute nonequilibrium viscosity with cooling rates from 10-12 to 1012K/s . Based on these detailed landscape calculations, we propose a simplified phenomenological model capturing the essential physics of glass viscosity. The phenomenological model incorporates an ergodicity parameter that accounts for the continuous breakdown of ergodicity at the glass transition. We show a direct relationship between the nonequilibrium viscosity parameters and the fragility of the supercooled liquid. The nonequilibrium viscosity model is validated against experimental measurements of Corning EAGLE XG™ glass. The measurements are performed using a specially designed beam-bending apparatus capable of accurate nonequilibrium viscosity measurements up to 1016Pas . Using a common set of parameters, the phenomenological model provides an accurate description of EAGLE XG™ viscosity over the full range of measured temperatures and fictive temperatures.

Inference of viscosity jump at 670 km depth and lower mantle viscosity structure from GIA observations

NASA Astrophysics Data System (ADS)

Nakada, Masao; Okuno, Jun'ichi; Irie, Yoshiya

2018-03-01

A viscosity model with an exponential profile described by temperature (T) and pressure (P) distributions and constant activation energy (E_{{{um}}}^{{*}} for the upper mantle and E_{{{lm}}}^* for the lower mantle) and volume (V_{{{um}}}^{{*}} and V_{{{lm}}}^*) is employed in inferring the viscosity structure of the Earth's mantle from observations of glacial isostatic adjustment (GIA). We first construct standard viscosity models with an average upper-mantle viscosity ({\\bar{η }_{{{um}}}}) of 2 × 1020 Pa s, a typical value for the oceanic upper-mantle viscosity, satisfying the observationally derived three GIA-related observables, GIA-induced rate of change of the degree-two zonal harmonic of the geopotential, {\\dot{J}_2}, and differential relative sea level (RSL) changes for the Last Glacial Maximum sea levels at Barbados and Bonaparte Gulf in Australia and for RSL changes at 6 kyr BP for Karumba and Halifax Bay in Australia. Standard viscosity models inferred from three GIA-related observables are characterized by a viscosity of ˜1023 Pa s in the deep mantle for an assumed viscosity at 670 km depth, ηlm(670), of (1 - 50) × 1021 Pa s. Postglacial RSL changes at Southport, Bermuda and Everglades in the intermediate region of the North American ice sheet, largely dependent on its gross melting history, have a crucial potential for inference of a viscosity jump at 670 km depth. The analyses of these RSL changes based on the viscosity models with {\\bar{η }_{{{um}}}} ≥ 2 × 1020 Pa s and lower-mantle viscosity structures for the standard models yield permissible {\\bar{η }_{{{um}}}} and ηlm (670) values, although there is a trade-off between the viscosity and ice history models. Our preferred {\\bar{η }_{{{um}}}} and ηlm (670) values are ˜(7 - 9) × 1020 and ˜1022 Pa s, respectively, and the {\\bar{η }_{{{um}}}} is higher than that for the typical value of oceanic upper mantle, which may reflect a moderate laterally heterogeneous upper-mantle viscosity. The mantle viscosity structure adopted in this study depends on temperature distribution and activation energy and volume, and it is difficult to discuss the impact of each quantity on the inferred lower-mantle viscosity model. We conclude that models of smooth depth variation in the lower-mantle viscosity following η ( z ) ∝ {{ exp}}[ {( {E_{{{lm}}}^* + P( z )V_{{{lm}}}^*} )/{{R}}T( z )} ] with constant E_{{{lm}}}^* and V_{{{lm}}}^* are consistent with the GIA observations.

Thermal conductivity and viscosity measurements of ethylene glycol-based Al2O3 nanofluids

NASA Astrophysics Data System (ADS)

Pastoriza-Gallego, María José; Lugo, Luis; Legido, José Luis; Piñeiro, Manuel M.

2011-12-01

The dispersion and stability of nanofluids obtained by dispersing Al2O3 nanoparticles in ethylene glycol have been analyzed at several concentrations up to 25% in mass fraction. The thermal conductivity and viscosity were experimentally determined at temperatures ranging from 283.15 K to 323.15 K using an apparatus based on the hot-wire method and a rotational viscometer, respectively. It has been found that both thermal conductivity and viscosity increase with the concentration of nanoparticles, whereas when the temperature increases the viscosity diminishes and the thermal conductivity rises. Measured enhancements on thermal conductivity (up to 19%) compare well with literature values when available. New viscosity experimental data yield values more than twice larger than the base fluid. The influence of particle size on viscosity has been also studied, finding large differences that must be taken into account for any practical application. These experimental results were compared with some theoretical models, as those of Maxwell-Hamilton and Crosser for thermal conductivity and Krieger and Dougherty for viscosity.

Compositional dependence of lower crustal viscosity

NASA Astrophysics Data System (ADS)

Shinevar, William J.; Behn, Mark D.; Hirth, Greg

2015-10-01

We calculate the viscosity structure of the lower continental crust as a function of its bulk composition using multiphase mixing theory. We use the Gibbs free-energy minimization routine Perple_X to calculate mineral assemblages for different crustal compositions under pressure and temperature conditions appropriate for the lower continental crust. The effective aggregate viscosities are then calculated using a rheologic mixing model and flow laws for the major crust-forming minerals. We investigate the viscosity of two lower crustal compositions: (i) basaltic (53 wt % SiO2) and (ii) andesitic (64 wt % SiO2). The andesitic model predicts aggregate viscosities similar to feldspar and approximately 1 order of magnitude greater than that of wet quartz. The viscosity range calculated for the andesitic crustal composition (particularly when hydrous phases are stable) is most similar to independent estimates of lower crust viscosity in actively deforming regions based on postglacial isostatic rebound, postseismic relaxation, and paleolake shoreline deflection.

Viscosity of alumina nanoparticles dispersed in car engine coolant

DOE Office of Scientific and Technical Information (OSTI.GOV)

Kole, Madhusree; Dey, T.K.

The present paper, describes our experimental results on the viscosity of the nanofluid prepared by dispersing alumina nanoparticles (<50 nm) in commercial car coolant. The nanofluid prepared with calculated amount of oleic acid (surfactant) was tested to be stable for more than 80 days. The viscosity of the nanofluids is measured both as a function of alumina volume fraction and temperature between 10 and 50 C. While the pure base fluid display Newtonian behavior over the measured temperature, it transforms to a non-Newtonian fluid with addition of a small amount of alumina nanoparticles. Our results show that viscosity of themore » nanofluid increases with increasing nanoparticle concentration and decreases with increase in temperature. Most of the frequently used classical models severely under predict the measured viscosity. Volume fraction dependence of the nanofluid viscosity, however, is predicted fairly well on the basis of a recently reported theoretical model for nanofluids that takes into account the effect of Brownian motion of nanoparticles in the nanofluid. The temperature dependence of the viscosity of engine coolant based alumina nanofluids obeys the empirical correlation of the type: log ({mu}{sub nf}) = A exp(BT), proposed earlier by Namburu et al. (author)« less

Reynolds Stress Distributions and the Measurement and Calculation of Eddy Viscosity in Gravity Currents

NASA Astrophysics Data System (ADS)

Kelly, R. W.; Chalk, C.; Dorrell, R. M.; Peakall, J.; Burns, A. D.; Keevil, G. M.; Thomas, R. E.; Williams, G.

2016-12-01

In the natural environment, gravity currents transport large volumes of sediment great distances and are often considered one of the most important mechanisms for sediment transport in ocean basins. Deposits from many individual submarine gravity currents, turbidites, ultimately form submarine fan systems. These are the largest sedimentary systems on the planet and contain valuable hydrocarbon reserves. Moreover, the impact of these currents on submarine technologies and seafloor infrastructure can be devastating and therefore they are of significant interest to a wide range of industries. Here we present experimental, numerical and theoretical models of time-averaged turbulent shear stresses, i.e. Reynolds stresses. Reynolds stresses can be conceptually parameterised by an eddy viscosity parameter that relates chaotic fluid motion to diffusive type processes. As such, it is a useful parameter for indicating the extent of internal mixing and is used extensively in both numerical and analytical modelling of both open-channel and gravity driven flows. However, a lack of knowledge of the turbulent structure of gravity currents limits many hydro- and morphodynamic models. High resolution 3-dimensional experimental velocity data, gathered using acoustic Doppler profiling velocimetry, enabled direct calculation of stresses and eddy viscosity. Comparison of experimental data to CFD and analytical models allowed the testing of eddy viscosity-based turbulent mixing models. The calculated eddy viscosity profile is parabolic in nature in both the upper and lower shear layers. However, an apparent breakdown in the Boussinesq hypothesis (used to calculate the eddy viscosity and upon which many numerical models are based) is observed in the region of the current around the velocity maximum. With the use of accompanying density data it is suggested that the effect of stratification on eddy viscosity is significant and alternative formulations may be required.

On the resolution of shallow mantle viscosity structure using post-earthquake relaxation data: Application to the 1999 Hector Mine, California, earthquake

USGS Publications Warehouse

Pollitz, Fred F.; Thatcher, Wayne R.

2010-01-01

Most models of lower crust/mantle viscosity inferred from postearthquake relaxation assume one or two uniform-viscosity layers. A few existing models possess apparently significant radially variable viscosity structure in the shallow mantle (e.g., the upper 200 km), but the resolution of such variations is not clear. We use a geophysical inverse procedure to address the resolving power of inferred shallow mantle viscosity structure using postearthquake relaxation data. We apply this methodology to 9 years of GPS-constrained crustal motions after the 16 October 1999 M = 7.1 Hector Mine earthquake. After application of a differencing method to isolate the postearthquake signal from the “background” crustal velocity field, we find that surface velocities diminish from ∼20 mm/yr in the first few months to ≲2 mm/yr after 2 years. Viscoelastic relaxation of the mantle, with a time-dependent effective viscosity prescribed by a Burgers body, provides a good explanation for the postseismic crustal deformation, capturing both the spatial and temporal pattern. In the context of the Burgers body model (which involves a transient viscosity and steady state viscosity), a resolution analysis based on the singular value decomposition reveals that at most, two constraints on depth-dependent steady state mantle viscosity are provided by the present data set. Uppermost mantle viscosity (depth ≲ 60 km) is moderately resolved, but deeper viscosity structure is poorly resolved. The simplest model that explains the data better than that of uniform steady state mantle viscosity involves a linear gradient in logarithmic viscosity with depth, with a small increase from the Moho to 220 km depth. However, the viscosity increase is not statistically significant. This suggests that the depth-dependent steady state viscosity is not resolvably different from uniformity in the uppermost mantle.

Searching for 3D Viscosity Models of Glacial Isostatic Adjustment in Support of the Global ICE-6G_C Ice History Model

NASA Astrophysics Data System (ADS)

LI, T., II; Wu, P.; Steffen, H.; Wang, H.

2017-12-01

The global ice history model ICE-6G_C was constructed based on the laterally homogeneous earth model VM5a. The combined model of glacial isostatic adjustment (GIA) called ICE-6G_C (VM5a) fits global observations of GIA simultaneously well. However, seismic and geological observations clearly show that the Earth's mantle is laterally heterogeneous. Our aim therefore is to search for the best laterally heterogeneous viscosity models with ICE-6G_C ice history that is able to fit the global relative sea-level (RSL) data, the peak uplift rates (from GNSS) and peak g-dot rates (from the GRACE satellite mission) in Laurentia and Fennoscandia simultaneously. The Coupled Laplace-Finite Element Method is used to compute gravitationally self-consistent sea levels with time dependent coastlines and rotational feedback in addition to changes in deformation, gravity and the state of stress. As a start, the VM5a Earth model is used as the radial background viscosity structure but other radial background viscosity models will also be investigated. Lateral mantle viscosity structure is obtained by the superposition of the radial background viscosity and the lateral viscosity perturbations logarithmically. The latter is inferred from a seismic tomography model using a scaling relationship that takes into account the effects of anharmonicity, anelasticity and non-thermal effects. We will show that several laterally heterogeneous mantle viscosity models can fit the global sea level, GPS and GRACE data better than laterally homogeneous models, provided that the scaling relationship for mantle heterogeneity under northern Europe is allowed to be different from that under Laurentia. In addition, the effects of laterally heterogeneous lithosphere, as inferred from seismic tomography, and the lateral changes in sub-lithospheric properties will also be presented.

A 3D Finite Element Model with Improved Spatial Resolution to Investigate the Effect of Varying Viscosity on Antarctica

NASA Astrophysics Data System (ADS)

Blank, B.; van der Wal, W.; Pappa, F.; Ebbing, J.

2017-12-01

B. Blank1, H. Hu1, W. van der Wal1, F Pappa2, J. Ebbing21Delft University of Technology 2Christian-Albrechts-University of KielSince the beginning of the 2000's time-variable gravity data from GRACE has proved to be an effective method for mapping ice mass loss in Antarctica. However, Glacial Isostatic Adjustment (GIA) models are required to correct for GIA induced mass changes. While most GIA models have adopted an Earth model that only varies radially in parameters, it has long been clear that the Earth structure also varies with longitude and latitude. For this study a new global 3D GIA model has been developed within the finite element software package ABAQUS, which can be modified to operate on a spatial resolution down to 50 km locally. The model is being benchmarked against normal model models for surface loading. It will be used to investigate the effects of a 3D varying lithosphere and upper asthenosphere in Antarctica. Viscosity which will be computed from temperature estimates with laboratory based flow laws. A new 3D temperature map of the Antarctic lithosphere has been developed within ESA's GOCE+ project based on seismic data as well as on GOCE and GRACE inferred gravity gradients. Output from the GIA model with this new temperature estimates will be compared to that of 1D viscosity profiles and other recent 3D viscosity models based on seismic data. From these side to side comparisons we want to investigate the influence of the viscosity map on uplift rates and horizontal movement. Finally the results can be compared to GPS measurement to investigate the validity of all models.

Glass viscosity calculation based on a global statistical modelling approach

DOE Office of Scientific and Technical Information (OSTI.GOV)

Fluegel, Alex

2007-02-01

A global statistical glass viscosity model was developed for predicting the complete viscosity curve, based on more than 2200 composition-property data of silicate glasses from the scientific literature, including soda-lime-silica container and float glasses, TV panel glasses, borosilicate fiber wool and E type glasses, low expansion borosilicate glasses, glasses for nuclear waste vitrification, lead crystal glasses, binary alkali silicates, and various further compositions from over half a century. It is shown that within a measurement series from a specific laboratory the reported viscosity values are often over-estimated at higher temperatures due to alkali and boron oxide evaporation during the measurementmore » and glass preparation, including data by Lakatos et al. (1972) and the recently published High temperature glass melt property database for process modeling by Seward et al. (2005). Similarly, in the glass transition range many experimental data of borosilicate glasses are reported too high due to phase separation effects. The developed global model corrects those errors. The model standard error was 9-17°C, with R^2 = 0.985-0.989. The prediction 95% confidence interval for glass in mass production largely depends on the glass composition of interest, the composition uncertainty, and the viscosity level. New insights in the mixed-alkali effect are provided.« less

Length-scales of Slab-induced Asthenospheric Deformation from Geodynamic Modeling, Mantle Deformation Fabric, and Synthetic Shear Wave Splitting

NASA Astrophysics Data System (ADS)

Jadamec, M. A.; MacDougall, J.; Fischer, K. M.

2017-12-01

The viscosity structure of the Earth's interior is critically important, because it places a first order constraint on plate motion and mantle flow rates. Geodynamic models using a composite viscosity based on experimentally derived flow laws for olivine aggregates show that lateral viscosity variations emerge in the upper mantle due to the subduction dynamics. However, the length-scale of this transition is still not well understood. Two-dimensional numerical models of subduction are presented that investigate the effect of initial slab dip, maximum yield stress (slab strength), and viscosity formulation (Newtonian versus composite) on the emergent lateral viscosity variations in the upper-mantle and magnitude of slab-driven mantle flow velocity. Significant viscosity reductions occur in regions of large flow velocity gradients due to the weakening effect of the dislocation creep deformation mechanism. The dynamic reductions in asthenospheric viscosity (less than 1018 Pa s) occur within approximately 500 km from driving force of the slab, with peak flow velocities occurring in models with a lower yield stress (weaker slab) and higher stress exponent. This leads to a sharper definition of the rheological base of the lithosphere and implies lateral variability in tractions along the base of the lithosphere. As the dislocation creep mechanism also leads to mantle deformation fabric, we then examine the spatial variation in the LPO development in the asthenosphere and calculate synthetic shear wave splitting. The models show that olivine LPO fabric in the asthenosphere generally increases in alignment strength with increased proximity to the slab, but can be transient and spatially variable on small length scales. The vertical flow fields surrounding the slab tip can produce shear-wave splitting variations with back-azimuth that deviate from the predictions of uniform trench-normal anisotropy, a result that bears on the interpretation of complexity in shear-wave splitting observed in real subduction zones.

The Influence of Wall Conductivity of Film Condensation with Integral Fin Tubes

DTIC Science & Technology

1993-09-23

tube based on Nusselt theory ) dynamic viscosity, kg/(m*s) Mf dynamic viscosity of the condensate film, kg/(m*s) Aw dynamic viscosity of the cooling...improve the simple model of Nusselt to predict the heat transfer 2 coefficient for condensation on horizontal tubes. Nusselt’s theory was based on a plain...be developed and utilized. 1. Norisontal Smooth Tubes Nusselt [Ref. 16] developed the foundation for the study of filmwise condensation on horizontal

Thermal conductivity and viscosity measurements of ethylene glycol-based Al2O3 nanofluids

PubMed Central

2011-01-01

The dispersion and stability of nanofluids obtained by dispersing Al2O3 nanoparticles in ethylene glycol have been analyzed at several concentrations up to 25% in mass fraction. The thermal conductivity and viscosity were experimentally determined at temperatures ranging from 283.15 K to 323.15 K using an apparatus based on the hot-wire method and a rotational viscometer, respectively. It has been found that both thermal conductivity and viscosity increase with the concentration of nanoparticles, whereas when the temperature increases the viscosity diminishes and the thermal conductivity rises. Measured enhancements on thermal conductivity (up to 19%) compare well with literature values when available. New viscosity experimental data yield values more than twice larger than the base fluid. The influence of particle size on viscosity has been also studied, finding large differences that must be taken into account for any practical application. These experimental results were compared with some theoretical models, as those of Maxwell-Hamilton and Crosser for thermal conductivity and Krieger and Dougherty for viscosity. PMID:21711737

Shear Rheology of Suspensions of Porous Zeolite Particles in Concentrated Polymer Solutions

NASA Astrophysics Data System (ADS)

Olanrewaju, Kayode O.; Breedveld, Victor

2008-07-01

We present experimental data on the shear rheology of Ultem (polyetherimide)/NMP(l-methyl-2-pyrrolidinone) solutions with and without suspended surface-modified porous/nonporous zeolite (ZSM-5) particles. We found that the porous zeolite suspensions have relative viscosities that significantly exceed the Krieger-Dougherty predictions for hard sphere suspensions. The major origin of this discrepancy is the selective absorption of NMP solvent into the zeolite pores, which raises both the polymer concentration and the particle volume fraction, thus enhancing both the viscosity of the continuous phase Ultem/NMP polymer solution and the particle contribution to the suspension viscosity. Other factors, such as zeolite non-sphericity and specific interactions with Ultem polymer, contribute to the suspension viscosity to a lesser extent. We propose a predictive model for the viscosity of porous zeolite suspensions by incorporating an absorption parameter, α, into the Krieger-Dougherty model. We also propose independent approaches to determine α. The first one is indirect and based on zeolite density/porosity data, assuming that all pores will be filled with solvent. The other method is based on our experimental data, by comparing the viscosity data of porous versus non-porous zeolite suspensions. The different approaches are compared.

Microfluidics for simultaneous quantification of platelet adhesion and blood viscosity

PubMed Central

Yeom, Eunseop; Park, Jun Hong; Kang, Yang Jun; Lee, Sang Joon

2016-01-01

Platelet functions, including adhesion, activation, and aggregation have an influence on thrombosis and the progression of atherosclerosis. In the present study, a new microfluidic-based method is proposed to estimate platelet adhesion and blood viscosity simultaneously. Blood sample flows into an H-shaped microfluidic device with a peristaltic pump. Since platelet aggregation may be initiated by the compression of rotors inside the peristaltic pump, platelet aggregates may adhere to the H-shaped channel. Through correlation mapping, which visualizes decorrelation of the streaming blood flow, the area of adhered platelets (APlatelet) can be estimated without labeling platelets. The platelet function is estimated by determining the representative index IA·T based on APlatelet and contact time. Blood viscosity is measured by monitoring the flow conditions in the one side channel of the H-shaped device. Based on the relation between interfacial width (W) and pressure ratio of sample flows to the reference, blood sample viscosity (μ) can be estimated by measuring W. Biophysical parameters (IA·T, μ) are compared for normal and diabetic rats using an ex vivo extracorporeal model. This microfluidic-based method can be used for evaluating variations in the platelet adhesion and blood viscosity of animal models with cardiovascular diseases under ex vivo conditions. PMID:27118101

Experience in using a numerical scheme with artificial viscosity at solving the Riemann problem for a multi-fluid model of multiphase flow

NASA Astrophysics Data System (ADS)

Bulovich, S. V.; Smirnov, E. M.

2018-05-01

The paper covers application of the artificial viscosity technique to numerical simulation of unsteady one-dimensional multiphase compressible flows on the base of the multi-fluid approach. The system of the governing equations is written under assumption of the pressure equilibrium between the "fluids" (phases). No interfacial exchange is taken into account. A model for evaluation of the artificial viscosity coefficient that (i) assumes identity of this coefficient for all interpenetrating phases and (ii) uses the multiphase-mixture Wood equation for evaluation of a scale speed of sound has been suggested. Performance of the artificial viscosity technique has been evaluated via numerical solution of a model problem of pressure discontinuity breakdown in a three-fluid medium. It has been shown that a relatively simple numerical scheme, explicit and first-order, combined with the suggested artificial viscosity model, predicts a physically correct behavior of the moving shock and expansion waves, and a subsequent refinement of the computational grid results in a monotonic approaching to an asymptotic time-dependent solution, without non-physical oscillations.

Integrated Solvent Design for CO 2 Capture and Viscosity Tuning

DOE PAGES

Cantu, David C.; Malhotra, Deepika; Koech, Phillip K.; ...

2017-08-18

We present novel design strategies for reduced viscosity single-component, water-lean CO 2 capture organic solvent systems. Through molecular simulation, we identify the main molecular-level descriptor that influences bulk solvent viscosity. Upon loading, a zwitterionic structure forms with a small activation energy of ca 16 kJ/mol and a small stabilization of ca 6 kJ/mol. Viscosity increases exponentially with CO 2 loading due to hydrogen-bonding between neighboring Zwitterions. We find that molecular structures that promote internal hydrogen bonding (within the same molecule) and suppress interactions with neighboring molecules have low viscosities. In addition, tuning the acid/base properties leads to a shift ofmore » the equilibrium toward a non-charged (acid) form that further reduces the viscosity. Here, based on the above structural criteria, a reduced order model is also presented that allows for the quick screening of large compound libraries and down selection of promising candidates for synthesis and testing.« less

Dynamic contact angle of water-based titanium oxide nanofluid

PubMed Central

2013-01-01

This paper presents an investigation into spreading dynamics and dynamic contact angle of TiO2-deionized water nanofluids. Two mechanisms of energy dissipation, (1) contact line friction and (2) wedge film viscosity, govern the dynamics of contact line motion. The primary stage of spreading has the contact line friction as the dominant dissipative mechanism. At the secondary stage of spreading, the wedge film viscosity is the dominant dissipative mechanism. A theoretical model based on combination of molecular kinetic theory and hydrodynamic theory which incorporates non-Newtonian viscosity of solutions is used. The model agreement with experimental data is reasonable. Complex interparticle interactions, local pinning of the contact line, and variations in solid–liquid interfacial tension are attributed to errors. PMID:23759071

Experimental and Theoretical Investigations on Viscosity of Fe-Ni-C Liquids at High Pressures

NASA Astrophysics Data System (ADS)

Chen, B.; Lai, X.; Wang, J.; Zhu, F.; Liu, J.; Kono, Y.

2016-12-01

Understanding and modeling of Earth's core processes such as geodynamo and heat flow via convection in liquid outer cores hinges on the viscosity of candidate liquid iron alloys under core conditions. Viscosity estimates from various methods of the metallic liquid of the outer core, however, span up to 12 orders of magnitude. Due to experimental challenges, viscosity measurements of iron liquids alloyed with lighter elements are scarce and conducted at conditions far below those expected for the outer core. In this study, we adopt a synergistic approach by integrating experiments at experimentally-achievable conditions with computations up to core conditions. We performed viscosity measurements based on the modified Stokes' floating sphere viscometry method for the Fe-Ni-C liquids at high pressures in a Paris-Edinburgh press at Sector 16 of the Advanced Photon Source, Argonne National Laboratory. Our results show that the addition of 3-5 wt.% carbon to iron-nickel liquids has negligible effect on its viscosity at pressures lower than 5 GPa. The viscosity of the Fe-Ni-C liquids, however, becomes notably higher and increases by a factor of 3 at 5-8 GPa. Similarly, our first-principles molecular dynamics calculations up to Earth's core pressures show a viscosity change in Fe-Ni-C liquids at 5 GPa. The significant change in the viscosity is likely due to a liquid structural transition of the Fe-Ni-C liquids as revealed by our X-ray diffraction measurements and first-principles molecular dynamics calculations. The observed correlation between structure and physical properties of liquids permit stringent benchmark test of the computational liquid models and contribute to a more comprehensive understanding of liquid properties under high pressures. The interplay between experiments and first-principles based modeling is shown to be a practical and effective methodology for studying liquid properties under outer core conditions that are difficult to reach with the current static high-pressure capabilities. The new viscosity data from experiments and computations would provide new insights into the internal dynamics of the outer core.

Shear Viscosity of Aluminium under Shock Compression

NASA Astrophysics Data System (ADS)

Liu, Fu-Sheng; Yang, Mei-Xia; Liu, Qi-Wen; Chen, Jun-Xiang; Jing, Fu-Qian

2005-03-01

Based on the Newtonian viscous fluid model and the analytic perturbation theory of Miller and Ahrens for the oscillatory damping of a sinusoidal shock front, a flyer-impact technique is developed to investigate the effective viscosity of shocked aluminium. The shear viscosity coefficient is determined to be about 5000 poises at 42 GPa with strain rate of 1.27×106 s-1, which is a reasonable estimation compared with the results of other measurement methods.

In vitro digestion of short-dough biscuits enriched in proteins and/or fibres, using a multi-compartmental and dynamic system (1): viscosity measurement and prediction.

PubMed

Villemejane, C; Wahl, R; Aymard, P; Denis, S; Michon, C

2015-09-01

The effects of biscuit composition on the viscosity generated during digestion were investigated. A control biscuit, one with proteins, one with fibres, and one with both proteins and fibres were digested under the same conditions, using the TNO intestinal model (TIM-1). The TIM-1 is a multi-compartmental and dynamic in vitro system, simulating digestion in the upper tract (stomach and small intestine) of healthy adult humans. Digesta were collected at different times, in the different compartments of the TIM-1 (stomach, duodenum, jejunum and ileum) and viscosity was measured with a dynamic rheometer. Results showed a marked effect of biscuit composition on chyme viscosity. Highest viscosity was obtained with biscuits containing viscous soluble fibres, followed by those enriched in both proteins and fibres, then by protein-enriched and control biscuits. The viscosity was maintained throughout the gut up to the ileal compartment. A prediction of the evolution of the chyme viscosity in each compartment of the TIM-1 was built, based on model curves describing the evolution of the viscosity as a function of biscuit concentration, and on dilution factors measured by spectrophotometry on a blank digestion. Copyright © 2015 Elsevier Ltd. All rights reserved.

Effect of Qingnao tablet on blood viscosity of rat model of blood stasis induced by epinephrine

NASA Astrophysics Data System (ADS)

Xie, Guoqi; Hao, Shaojun; Ma, Zhenzhen; Liu, Xiaobin; Li, Jun; Li, Wenjun; Zhang, Zhengchen

2018-04-01

To establish a rat model of blood stasis with adrenaline (Adr) subcutaneous injection and ice bath stimulation. The effects of different doses on the blood viscosity of blood stasis model rats were observed. The rats were randomly divided into 6 groups: blank control group (no model), model group, positive control group, high, middle and low dose group. The whole blood viscosity and plasma viscosity were detected by blood viscosity instrument. Compared with the blank group, model group, high shear, low shear whole blood viscosity and plasma viscosity were significantly increased, TT PT significantly shortened, APTT was significantly prolonged, FIB increased significantly, indicating that the model was successful. Compared with the model group, can significantly reduce the Naoluotong group (cut, low cut). Qingnaopian high dose group (low cut), middle dose group (cut, low shear blood viscosity) (P<0.01), Can significantly reduce Naoluotong qingnaopian group, high dose group (P<0.01), plasma viscosity decreased qingnaopian plasma viscosity in low dose group (P<0.05). Conclusion: qingnaopian could improve the blood rheology of blood stasis mice abnormal index, decrease the blood viscosity, blood stasis has certain hemostatic effect.

High-Resolution Lithosphere Viscosity and Dynamics Revealed by Magnetotelluric Imaging

NASA Astrophysics Data System (ADS)

Liu, L.; Hasterok, D. P.

2016-12-01

An accurate viscosity structure is critical to truthfully modeling continental lithosphere dynamics, especially at spatial scales of <200 km where active tectonic deformation and volcanism occur. However, the effective viscosity structure of the lithosphere remains a key challenge in geodynamics due to the intimate involvement of viscosity with time and its dependence on many factors including strain rate, plastic failure, composition, and grain size. Current efforts on inferring the detailed lithosphere viscosity structure are sparse and large uncertainties and discrepancies still exist. Here we report an attempt to infer the effective lithospheric viscosity from a high-resolution magnetotelluric (MT) survey across the western United States. The high sensitivity of MT fields to the presence of electrically conductive fluids makes it a promising proxy for determining mechanical strength variations throughout the lithosphere. We demonstrate how a viscosity structure, approximated from electrical resistivity, results in a geodynamic model that successfully predicts short-wavelength surface topography, lithospheric deformation, and mantle upwelling beneath recent volcanism. The results indicate that lithosphere viscosity structure rather than the buoyancy structure is the dominant controlling factor for short-wavelength topography and intra-plate deformation in tectonically active regions. We further show that this viscosity is consistent with and more effective than that derived from laboratory-based rheology. We therefore propose that MT imaging provides a practical observational constraint for quantifying the dynamic evolution of the continental lithosphere.

Modelling of wicking and moisture interactions of flax and viscose fibres.

PubMed

Stuart, T; McCall, R D; Sharma, H S S; Lyons, G

2015-06-05

Methods for assessing the wicking properties of individual fibre bundles have been developed from models based on the original Washburn equation (WE) and the modified Washburn equation (MWE), which also accounts for swelling. Both models gave indication of differences in wicking properties of flax and the viscose fibres, though MWE gave additional information that could be interpreted in terms of the physical model. Wicking of the viscose fibres is mainly via inter-fibre capillaries while that of flax is a combination of inter-fibre capillaries and lumen present in some elementary fibres. The degree of swelling and associated rotation of flax fibre in a vapour pressure range of 1-6torr were monitored using an environmental scanning electron microscope (ESEM). Viscose fibre exhibited swelling under the same conditions but did not rotate. The two techniques highlighted different mechanisms of wicking which can be used for monitoring moisture uptake/swelling of treated fibres for fabrication of composites. Copyright © 2015 Elsevier Ltd. All rights reserved.

Compliance of the Stokes-Einstein model and breakdown of the Stokes-Einstein-Debye model for a urea-based supramolecular polymer of high viscosity.

PubMed

Świergiel, Jolanta; Bouteiller, Laurent; Jadżyn, Jan

2014-11-14

Impedance spectroscopy was used for the study of the static and dynamic behavior of the electrical conductivity of a hydrogen-bonded supramolecular polymer of high viscosity. The experimental data are discussed in the frame of the Stokes-Einstein and Stokes-Einstein-Debye models. It was found that the translational movement of the ions is due to normal Brownian diffusion, which was revealed by a fulfillment of Ohm's law by the electric current and a strictly exponential decay of the current after removing the electric stimulus. The dependence of the dc conductivity on the viscosity of the medium fulfills the Stokes-Einstein model quite well. An extension of the model, by including in it the conductivity relaxation time, is proposed in this paper. A breakdown of the Stokes-Einstein-Debye model is revealed by the relations of the dipolar relaxation time to the viscosity and to the dc ionic conductivity. The importance of the C=O···H-N hydrogen bonds in that breakdown is discussed.

Viscosity Measurements of SiO2-"FeO"-MgO System in Equilibrium with Metallic Fe

NASA Astrophysics Data System (ADS)

Chen, Mao; Raghunath, Sreekanth; Zhao, Baojun

2014-01-01

The present study delivers the measurements of viscosities in the SiO2-"FeO"-MgO system in equilibrium with metallic Fe. The rotational spindle technique was used for the measurements at the temperature range of 1523 K to 1773 K (1250 °C to 1500 °C). Molybdenum crucibles and spindles were employed in all measurements. The viscosity measurements were carried out at 31 to 47 mol pct SiO2 and up to 18.8 mol pct MgO. Analysis of the quenched sample by Electron probe X-ray microanalysis after the viscosity measurement enables the composition and microstructure of the slag to be directly linked with the viscosity. The replacement of "FeO" by MgO was found to increase viscosity and activation energy of the SiO2-"FeO"-MgO slags. The modified Quasi-chemical Viscosity Model was further optimized in this system based on the current viscosity measurements.

Modeling and control of plasma rotation and βn for NSTX-U using Neoclassical Toroidal Viscosity and Neutral Beam Injection

NASA Astrophysics Data System (ADS)

Goumiri, Imene; Rowley, Clarence; Sabbagh, Steven; Gates, David; Gerhardt, Stefan; Boyer, Mark

2015-11-01

A model-based system is presented allowing control of the plasma rotation profile in a magnetically confined toroidal fusion device to maintain plasma stability for long pulse operation. The analysis, using NSTX data and NSTX-U TRANSP simulations, is aimed at controlling plasma rotation using momentum from six injected neutral beams and neoclassical toroidal viscosity generated by three-dimensional applied magnetic fields as actuators. Based on the momentum diffusion and torque balance model obtained, a feedback controller is designed and predictive simulations using TRANSP will be presented. Robustness of the model and the rotation controller will be discussed.

A Hermite-based lattice Boltzmann model with artificial viscosity for compressible viscous flows

NASA Astrophysics Data System (ADS)

Qiu, Ruofan; Chen, Rongqian; Zhu, Chenxiang; You, Yancheng

2018-05-01

A lattice Boltzmann model on Hermite basis for compressible viscous flows is presented in this paper. The model is developed in the framework of double-distribution-function approach, which has adjustable specific-heat ratio and Prandtl number. It contains a density distribution function for the flow field and a total energy distribution function for the temperature field. The equilibrium distribution function is determined by Hermite expansion, and the D3Q27 and D3Q39 three-dimensional (3D) discrete velocity models are used, in which the discrete velocity model can be replaced easily. Moreover, an artificial viscosity is introduced to enhance the model for capturing shock waves. The model is tested through several cases of compressible flows, including 3D supersonic viscous flows with boundary layer. The effect of artificial viscosity is estimated. Besides, D3Q27 and D3Q39 models are further compared in the present platform.

Chitosan solutions as injectable systems for dermal filler applications: Rheological characterization and biological evidence.

PubMed

Halimi, C; Montembault, A; Guerry, A; Delair, T; Viguier, E; Fulchiron, R; David, L

2015-01-01

A new generation of dermal filler for wrinkle filler based on chitosan was compared to current hyaluronic acid-based dermal fillers by using a new rheological performance criterion based on viscosity during injection related to Newtonian viscosity. In addition an in vivo evaluation was performed for preclinical evidence of chitosan use as dermal filler. In this way, biocompatibility and dermis reconstruction was evaluated on a pig model.

Viscosity and diffusivity in melts: from unary to multicomponent systems

NASA Astrophysics Data System (ADS)

Chen, Weimin; Zhang, Lijun; Du, Yong; Huang, Baiyun

2014-05-01

Viscosity and diffusivity, two important transport coefficients, are systematically investigated from unary melt to binary to multicomponent melts in the present work. By coupling with Kaptay's viscosity equation of pure liquid metals and effective radii of diffusion species, the Sutherland equation is modified by taking the size effect into account, and further derived into an Arrhenius formula for the convenient usage. Its reliability for predicting self-diffusivity and impurity diffusivity in unary liquids is then validated by comparing the calculated self-diffusivities and impurity diffusivities in liquid Al- and Fe-based alloys with the experimental and the assessed data. Moreover, the Kozlov model was chosen among various viscosity models as the most reliable one to reproduce the experimental viscosities in binary and multicomponent melts. Based on the reliable viscosities calculated from the Kozlov model, the modified Sutherland equation is utilized to predict the tracer diffusivities in binary and multicomponent melts, and validated in Al-Cu, Al-Ni and Al-Ce-Ni melts. Comprehensive comparisons between the calculated results and the literature data indicate that the experimental tracer diffusivities and the theoretical ones can be well reproduced by the present calculations. In addition, the vacancy-wind factor in binary liquid Al-Ni alloys with the increasing temperature is also discussed. What's more, the calculated inter-diffusivities in liquid Al-Cu, Al-Ni and Al-Ag-Cu alloys are also in excellent agreement with the measured and theoretical data. Comparisons between the simulated concentration profiles and the measured ones in Al-Cu, Al-Ce-Ni and Al-Ag-Cu melts are further used to validate the present calculation method.

Computing the Viscosity of Supercooled Liquids: Markov Network Model

PubMed Central

Li, Ju; Kushima, Akihiro; Eapen, Jacob; Lin, Xi; Qian, Xiaofeng; Mauro, John C.; Diep, Phong; Yip, Sidney

2011-01-01

The microscopic origin of glass transition, when liquid viscosity changes continuously by more than ten orders of magnitude, is challenging to explain from first principles. Here we describe the detailed derivation and implementation of a Markovian Network model to calculate the shear viscosity of deeply supercooled liquids based on numerical sampling of an atomistic energy landscape, which sheds some light on this transition. Shear stress relaxation is calculated from a master-equation description in which the system follows a transition-state pathway trajectory of hopping among local energy minima separated by activation barriers, which is in turn sampled by a metadynamics-based algorithm. Quantitative connection is established between the temperature variation of the calculated viscosity and the underlying potential energy and inherent stress landscape, showing a different landscape topography or “terrain” is needed for low-temperature viscosity (of order 107 Pa·s) from that associated with high-temperature viscosity (10−5 Pa·s). Within this range our results clearly indicate the crossover from an essentially Arrhenius scaling behavior at high temperatures to a low-temperature behavior that is clearly super-Arrhenius (fragile) for a Kob-Andersen model of binary liquid. Experimentally the manifestation of this crossover in atomic dynamics continues to raise questions concerning its fundamental origin. In this context this work explicitly demonstrates that a temperature-dependent “terrain” characterizing different parts of the same potential energy surface is sufficient to explain the signature behavior of vitrification, at the same time the notion of a temperature-dependent effective activation barrier is quantified. PMID:21464988

On the penetration of a hot diapir through a strongly temperature-dependent viscosity medium

NASA Technical Reports Server (NTRS)

Daly, S. F.; Raefsky, A.

1985-01-01

The ascent of a hot spherical body through a fluid with a strongly temperature-dependent viscosity has been studied using an axisymmetric finite element method. Numerical solutions range over Peclet numbers of 0.1 - 1000 from constant viscosity up to viscosity variations of 100,000. Both rigid and stress-free boundary conditions were applied at the surface of the sphere. The dependence of drag on viscosity variation was shown to have no dependence on the stress boundary condition except for a Stokes flow scaling factor. A Nusselt number parameterization based on the stress-free constant viscosity functional dependence on the Peclet number scaled by a parameter depending on the viscosity structure fits both stress-free and rigid boundary condition data above viscosity variations of 100. The temperature scale height was determined as a function of sphere radius. For the simple physical model studied in this paper pre-heating is required to reduce the ambient viscosity of the country rock to less than 10 to the 22nd sq cm/s in order for a 10 km diapir to penetrate a distance of several radii.

Effect of matrix elasticity on the continuous foaming of food models.

PubMed

Narchi, I; Vial, Ch; Djelveh, G

2008-12-01

The aim is to understand the effect of matrix elasticity on continuous foaming using food models based on glucose syrup. This was modified by adding polyacrylamide (PAA) with 2% whey protein isolate (WPI) or Tween 80 as foaming agents. Foaming was conducted in a stirred column. Rotation speed N and gas-to-liquid flow ratio (G/L) were varied. Overrun, average bubble size d (32), texture and stability were measured using densimetry, image analysis, and rheometry, respectively. Experimental results showed that 0.01% PAA did not modify the viscosity of 2% WPI models, but conferred low elastic behavior. PAA (0.05%) doubled matrix viscosity and drastically increased elasticity. The increase of elasticity became slower for further PAA addition. Foaming experiments demonstrated that theoretical overrun could not be achieved for inelastic WPI models in two cases: for high viscosity and low N, as dispersion effectiveness was reduced; for high G/L and N because of enhanced coalescence. Matrix elasticity was shown to increase overrun at constant viscosity for high G/L by enhancing interface stabilization. However, in elastic models, gas dispersion was more difficult and d (32) was higher than in inelastic fluids of similar viscosity. Finally, when the limiting step was dispersion, foaming was shown to be negatively affected by matrix elasticity.

GIA models with composite rheology and 3D viscosity: effect on GRACE mass balance in Antarctica

NASA Astrophysics Data System (ADS)

van der Wal, Wouter; Whitehouse, Pippa; Schrama, Ernst

2014-05-01

Most Glacial Isostatic Adjustment (GIA) models that have been used to correct GRACE data for the influence of GIA assume a radial stratification of viscosity in the Earth's mantle (1D viscosity). Seismic data in Antarctica indicate that there are large viscosity variations in the horizontal direction (3D viscosity). The purpose of this research is to determine the effect of 3D viscosity on GIA model output, and hence mass balance estimates in Antarctica. We use a GIA model with 3D viscosity and composite rheology in combination with ice loading histories ICE-5G and W12a. From comparisons with uplift and sea-level data in Fennoscandia and North America three preferred viscosity models are selected. For two of the 3D viscosity models the maximum gravity rate due to ICE-5G forcing is located over the Ronne-Filchner ice shelf. This is in contrast with the results obtained using a 1D model, in which the maximum gravity rate due to ICE-5G forcing is always located over the Ross ice shelf. This demonstrates that not all 3D viscosity models can be approximated with a 1D viscosity model. Using CSR release 5 GRACE data from February 2003 to June 2013 mass balance estimates for the three preferred viscosity models are -131 to -171 Gt/year for the ICE-5G model, and -48 to -57 Gt/year for the W12a model. The range due to Earth model uncertainty is larger than the error bar for GRACE (10 Gt/year), but smaller than the range resulting from the difference in ice loading histories.

Secular variations in zonal harmonics of Earth's geopotential and their implications for mantle viscosity and Antarctic melting history due to the last deglaciation

NASA Astrophysics Data System (ADS)

Nakada, Masao; Okuno, Jun'ichi

2017-06-01

Secular variations in zonal harmonics of Earth's geopotential based on the satellite laser ranging observations, {\\dot{J}_n}, contain important information about the Earth's deformation due to the glacial isostatic adjustment (GIA) and recent melting of glaciers and the Greenland and Antarctic ice sheets. Here, we examine the GIA-induced {\\dot{J}_n}, \\dot{J}_n^{GIA} (2 ≤ n ≤ 6), derived from the available geopotential zonal secular rate and recent melting taken from the IPCC 2013 Report (AR5) to explore the possibility of additional information on the depth-dependent lower-mantle viscosity and GIA ice model inferred from the analyses of the \\dot{J}_2^{GIA} and relative sea level changes. The sensitivities of the \\dot{J}_n^{GIA} to lower-mantle viscosity and GIA ice model with a global averaged eustatic sea level (ESL) of ∼130 m indicate that the secular rates for n = 3 and 4 are mainly caused by the viscous response of the lower mantle to the melting of the Antarctic ice sheet regardless of GIA ice models adopted in this study. Also, the analyses of the \\dot{J}_n^{GIA} based on the available geopotential zonal secular rates indicate that permissible lower-mantle viscosity structure satisfying even zonal secular rates of n = 2, 4 and 6 is obtained for the GIA ice model with an Antarctic ESL component of ∼20 or ∼30 m, but there is no viscosity solution satisfying \\dot{J}_3^{GIA} and \\dot{J}_5^{GIA} values. Moreover, the inference model for the lower-mantle viscosity and GIA ice model from each odd zonal secular rate is distinctly different from that satisfying GIA-induced even zonal secular rate. The discrepancy between the inference models for the even and odd zonal secular rates may partly be attributed to uncertainties of the geopotential zonal secular rates for n > 2 and particularly those for odd zonal secular rates due to weakness in the orbital geometry. If this problem is overcome at least for the secular rates of n < 5, then the analyses of the \\dot{J}_n^{GIA} would make it possible to put more convincing constraints on the lower-mantle viscosity structure and GIA ice model, particularly for the controversial Antarctic melting history in GIA community.

Modelling droplet collision outcomes for different substances and viscosities

NASA Astrophysics Data System (ADS)

Sommerfeld, Martin; Kuschel, Matthias

2016-12-01

The main objective of the present study is the derivation of models describing the outcome of binary droplet collisions for a wide range of dynamic viscosities in the well-known collision maps (i.e. normalised lateral droplet displacement at collision, called impact parameter, versus collision Weber number). Previous studies by Kuschel and Sommerfeld (Exp Fluids 54:1440, 2013) for different solution droplets having a range of solids contents and hence dynamic viscosities (here between 1 and 60 mPa s) revealed that the locations of the triple point (i.e. coincidence of bouncing, stretching separation and coalescence) and the critical Weber number (i.e. condition for the transition from coalescence to separation for head-on collisions) show a clear dependence on dynamic viscosity. In order to extend these findings also to pure liquids and to provide a broader data basis for modelling the viscosity effect, additional binary collision experiments were conducted for different alcohols (viscosity range 1.2-15.9 mPa s) and the FVA1 reference oil at different temperatures (viscosity range 3.0-28.2 mPa s). The droplet size for the series of alcohols was around 365 and 385 µm for the FVA1 reference oil, in each case with fixed diameter ratio at Δ= 1. The relative velocity between the droplets was varied in the range 0.5-3.5 m/s, yielding maximum Weber numbers of around 180. Individual binary droplet collisions with defined conditions were generated by two droplet chains each produced by vibrating orifice droplet generators. For recording droplet motion and the binary collision process with good spatial and temporal resolution high-speed shadow imaging was employed. The results for varied relative velocity and impact angle were assembled in impact parameter-Weber number maps. With increasing dynamic viscosity a characteristic displacement of the regimes for the different collision scenarios was also observed for pure liquids similar to that observed for solutions. This displacement could be described on a physical basis using the similarity number and structure parameter K which was obtained through flow process evaluation and optimal proportioning of momentum and energy by Naue and Bärwolff (Transportprozesse in Fluiden. Deutscher Verlag für Grundstoffindustrie GmbH, Leipzig 1992). Two correlations including the structure parameter K could be derived which describe the location of the triple point and the critical We number. All fluids considered, pure liquids and solutions, are very well fitted by these physically based correlations. The boundary model of Jiang et al. (J Fluid Mech 234:171-190, 1992) for distinguishing between coalescence and stretching separation could be adapted to go through the triple point by the two involved model parameters C a and C b, which were correlated with the relaxation velocity u_{{relax}} = {σ/μ}. Based on the predicted critical Weber number, denoting the onset of reflexive separation, the model of Ashgriz and Poo (J Fluid Mech 221:183-204, 1990) was adapted accordingly. The proper performance of the new generalised models was validated based on the present and previous measurements for a wide range of dynamic viscosities (i.e. 1-60 mPa s) and liquid properties. Although the model for the lower boundary of bouncing (Estrade et al. in J Heat Fluid Flow 20:486-491, 1999) could be adapted through the shape factor, it was found not suitable for the entire range of Weber numbers and viscosities.

Can a grain size-dependent viscosity help yielding realistic seismic velocities of LLSVPs?

NASA Astrophysics Data System (ADS)

Schierjott, J.; Cheng, K. W.; Rozel, A.; Tackley, P. J.

2017-12-01

Seismic studies show two antipodal regions of low shear velocity at the core-mantle boundary (CMB), one beneath the Pacific and one beneath Africa. These regions, called Large Low Shear Velocity Provinces (LLSVPs), are thought to be thermally and chemically distinct and thus have a different density and viscosity. Whereas there is some general consensus about the density of the LLSVPs the viscosity is still a very debated topic. So far, in numerical studies the viscosity is treated as either depth- and/or temperature- dependent but the potential grain size- dependence of the viscosity is neglected most of the time. In this study we use a self-consistent convection model which includes a grain size- dependent rheology based on the approach by Rozel et al. (2011) and Rozel (2012). Further, we consider a primordial layer and a time-dependent basalt production at the surface to dynamically form the present-day chemical heterogeneities, similar to earlier studies, e.g by Nakagawa & Tackley (2014). With this model we perform a parameter study which includes different densities and viscosities of the imposed primordial layer. We detect possible thermochemical piles based on different criterions, compute their average effective viscosity, density, rheology and grain size and investigate which detecting criterion yields the most realistic results. Our preliminary results show that a higher density and/or viscosity of the piles is needed to keep them at the core-mantle boundary (CMB). Relatively to the ambient mantle grain size is high in the piles but due to the temperature at the CMB the viscosity is not remarkably different than the one of ordinary plumes. We observe that grain size is lower if the density of the LLSVP is lower than the one of our MORB material. In that case the average temperature of the LLSVP is also reduced. Interestingly, changing the reference viscosity is responsible for a change in the average viscosity of the LLSVP but not for a different average grain size. Finally, we compare the numerical results with seismological observations by computing 1D seismic velocity profiles (p-wave, shear-wave and bulk velocities) inside and outside our detected piles using thermodynamic data calculated from Perple_X .

Correlation of shear and dielectric ion viscosity of dental resins - Influence of composition, temperature and filler content.

PubMed

Steinhaus, Johannes; Hausnerova, Berenika; Haenel, Thomas; Selig, Daniela; Duvenbeck, Fabian; Moeginger, Bernhard

2016-07-01

Shear viscosity and ion viscosity of uncured visible light-curing (VLC) resins and resin based composites (RBC) are correlated with respect to the resin composition, temperature and filler content to check where Dielectric Analysis (DEA) investigations of VLC RBC generate similar results as viscosity measurements. Mixtures of bisphenol A glycidyl methacrylate (Bis-GMA) and triethylene glycol dimethacrylate (TEGDMA) as well as the pure resins were investigated and compared with two commercial VLC dental resins and RBCs (VOCO, Arabesk Top and Grandio). Shear viscosity data was obtained using a Haake Mars III, Thermo Scientific. Ion viscosity measurements performed by a dielectric cure analyzer (DEA 231/1 Epsilon with Mini IDEX-Sensor, Netzsch-Gerätebau). Shear viscosity depends reciprocally on the mobility of molecules, whereas the ion viscosity also depends on the ion concentration as it is affected by both ion concentration and mixture viscosity. Except of pure TEGDMA, shear and ion viscosities depend on the resin composition qualitatively in a similar manner. Furthermore, shear and ion viscosities of the commercial VLC dental resins and composites exhibited the same temperature dependency regardless of filler content. Application of typical rheological models (Kitano and Quemada) revealed that ion viscosity measurements can be described with respect to filler contents of up to 30vol.%. Rheological behavior of a VLC RBC can be characterized by DEA under the condition that the ion concentration is kept constant. Both methods address the same physical phenomenon - motion of molecules. The proposed relations allows for calculating the viscosity of any Bis-GMA-TEGDMA mixture on the base of the viscosities of the pure components. This study demonstrated the applicability of DEA investigations of VLC RBCs with respect to quality assurance purposes. Copyright © 2016 The Academy of Dental Materials. Published by Elsevier Ltd. All rights reserved.

Viscosity of a concentrated suspension of rigid monosized particles

NASA Astrophysics Data System (ADS)

Brouwers, H. J. H.

2010-05-01

This paper addresses the relative viscosity of concentrated suspensions loaded with unimodal hard particles. So far, exact equations have only been put forward in the dilute limit, e.g., by Einstein [A. Einstein, Ann. Phys. 19, 289 (1906) (in German); Ann. Phys. 34, 591 (1911) (in German)] for spheres. For larger concentrations, a number of phenomenological models for the relative viscosity was presented, which depend on particle concentration only. Here, an original and exact closed form expression is derived based on geometrical considerations that predicts the viscosity of a concentrated suspension of monosized particles. This master curve for the suspension viscosity is governed by the relative viscosity-concentration gradient in the dilute limit (for spheres the Einstein limit) and by random close packing of the unimodal particles in the concentrated limit. The analytical expression of the relative viscosity is thoroughly compared with experiments and simulations reported in the literature, concerning both dilute and concentrated suspensions of spheres, and good agreement is found.

Correlation between the critical viscosity and ash fusion temperatures of coal gasifier ashes

DOE PAGES

Hsieh, Peter Y.; Kwong, Kyei-Sing; Bennett, James

2015-09-27

Coal gasification yields synthesis gas, an important intermediate in chemical manufacturing. It is also vital to the production of liquid fuels through the Fischer-Tropsch process and electricity in Integrated Gasification Combined Cycle power generation. Minerals naturally present in coal become molten in entrained-flow slagging gasifiers. Molten coal ash slag penetrates and dissolves refractory bricks, leading to costly plant shutdowns. The extent of coal ash slag penetration and refractory brick dissolution depends on the slag viscosity, the gasification temperature, and the composition of slag and bricks. We measured the viscosity of several synthetic coal ash slags with a high-temperature rotary viscometermore » and their ash fusion temperatures through optical image analysis. All measurements were made in a carbon monoxide-carbon dioxide reducing atmosphere that approximates coal gasification conditions. Empirical correlation models based on ash fusion temperatures were used to calculate critical viscosity temperatures based on the coal ash compositions. These values were then compared with those obtained from thermodynamic phase-transition models. An understanding of slag viscosity as a function of ash composition is important to reducing refractory wear in slagging coal gasifiers, which would help to reduce the cost and environmental impact of coal for chemical and electricity production.« less

Correlation between the critical viscosity and ash fusion temperatures of coal gasifier ashes

DOE Office of Scientific and Technical Information (OSTI.GOV)

Hsieh, Peter Y.; Kwong, Kyei-Sing; Bennett, James

Coal gasification yields synthesis gas, an important intermediate in chemical manufacturing. It is also vital to the production of liquid fuels through the Fischer-Tropsch process and electricity in Integrated Gasification Combined Cycle power generation. Minerals naturally present in coal become molten in entrained-flow slagging gasifiers. Molten coal ash slag penetrates and dissolves refractory bricks, leading to costly plant shutdowns. The extent of coal ash slag penetration and refractory brick dissolution depends on the slag viscosity, the gasification temperature, and the composition of slag and bricks. We measured the viscosity of several synthetic coal ash slags with a high-temperature rotary viscometermore » and their ash fusion temperatures through optical image analysis. All measurements were made in a carbon monoxide-carbon dioxide reducing atmosphere that approximates coal gasification conditions. Empirical correlation models based on ash fusion temperatures were used to calculate critical viscosity temperatures based on the coal ash compositions. These values were then compared with those obtained from thermodynamic phase-transition models. An understanding of slag viscosity as a function of ash composition is important to reducing refractory wear in slagging coal gasifiers, which would help to reduce the cost and environmental impact of coal for chemical and electricity production.« less

Development of Viscosity Model for Petroleum Industry Applications

NASA Astrophysics Data System (ADS)

Motahhari, Hamed reza

Heavy oil and bitumen are challenging to produce and process due to their very high viscosity, but their viscosity can be reduced either by heating or dilution with a solvent. Given the key role of viscosity, an accurate viscosity model suitable for use with reservoir and process simulators is essential. While there are several viscosity models for natural gases and conventional oils, a compositional model applicable to heavy petroleum and diluents is lacking. The objective of this thesis is to develop a general compositional viscosity model that is applicable to natural gas mixtures, conventional crudes oils, heavy petroleum fluids, and their mixtures with solvents and other crudes. The recently developed Expanded Fluid (EF) viscosity correlation was selected as a suitable compositional viscosity model for petroleum applications. The correlation relates the viscosity of the fluid to its density over a broad range of pressures and temperatures. The other inputs are pressure and the dilute gas viscosity. Each fluid is characterized for the correlation by a set of fluid-specific parameters which are tuned to fit data. First, the applicability of the EF correlation was extended to asymmetric mixtures and liquid mixtures containing dissolved gas components. A new set of mass-fraction based mixing rules was developed to calculate the fluid-specific parameters for mixtures. The EF correlation with the new set of mixing rules predicted the viscosity of over 100 mixtures of hydrocarbon compounds and carbon dioxide with overall average absolute relative deviations (AARD) of less than 10% either with measured densities or densities estimated by Advanced Peng-Robinson equation of state (APR EoS). To improve the viscosity predictions with APR EoS-estimated densities, general correlations were developed for non-zero viscosity binary interaction parameters. The EF correlation was extended to non-hydrocarbon compounds typically encountered in natural gas industry. It was demonstrated that the framework of the correlation is valid for these compounds, except for compounds with strong hydrogen bonding such as water. A temperature dependency was introduced into the correlation for strongly hydrogen bonding compounds. The EF correlation fit the viscosity data of pure non-hydrocarbon compounds with AARDs below 6% and predicted the viscosity of sour and sweet natural gases and aqueous solutions of organic alcohols with overall AARDs less than 9%. An internally consistent estimation method was also developed to calculate the fluid-specific parameters for hydrocarbons when no experimental viscosity data are available. The method correlates the fluid-specific parameters to the molecular weight and specific gravity. The method was evaluated against viscosity data of over 250 pure hydrocarbon compounds and petroleum distillations cuts. The EF correlation predictions were found to be within the same order of magnitude of the measurements with an overall AARD of 31%. A methodology was then proposed to apply the EF viscosity correlation to crude oils characterized as mixtures of the defined components and pseudo-components. The above estimation methods are used to calculate the fluid-specific parameters for pseudo-components. Guidelines are provided for tuning of the correlation to available viscosity data, calculating the dilute gas viscosities, and improving the densities calculated with the Peng-Robinson EoS. The viscosities of over 10 dead and live crude oils and bitumen were predicted within a factor of 3 of the measured values using the measured density of the oils as the input. It was shown that single parameter tuning of the model improved the viscosity prediction to within 30% of the measured values. Finally, the performance of the EF correlation was evaluated for diluted heavy oils and bitumens. The required density and viscosity data were collected for over 20 diluted dead and live bitumen mixtures using an in-house capillary viscometer also equipped with an in-line density-meter at temperatures and pressures up to 175 °C and 10 MPa. The predictions of the correlation were found within the same order of magnitude of the measured values with overall AARDs less than 20%. It was shown that the predictions of the correlation with generalized non-zero interaction parameters for the solvent-oil pairs were improved to overall AARDs less than 10%.

Torque Transient of Magnetically Drive Flow for Viscosity Measurement

NASA Technical Reports Server (NTRS)

Ban, Heng; Li, Chao; Su, Ching-Hua; Lin, Bochuan; Scripa, Rosalia N.; Lehoczky, Sandor L.

2004-01-01

Viscosity is a good indicator of structural changes for complex liquids, such as semiconductor melts with chain or ring structures. This paper discusses the theoretical and experimental results of the transient torque technique for non-intrusive viscosity measurement. Such a technique is essential for the high temperature viscosity measurement of high pressure and toxic semiconductor melts. In this paper, our previous work on oscillating cup technique was expanded to the transient process of a magnetically driven melt flow in a damped oscillation system. Based on the analytical solution for the fluid flow and cup oscillation, a semi-empirical model was established to extract the fluid viscosity. The analytical and experimental results indicated that such a technique has the advantage of short measurement time and straight forward data analysis procedures

Lattice Boltzmann simulations of immiscible displacement process with large viscosity ratios

NASA Astrophysics Data System (ADS)

Rao, Parthib; Schaefer, Laura

2017-11-01

Immiscible displacement is a key physical mechanism involved in enhanced oil recovery and carbon sequestration processes. This multiphase flow phenomenon involves a complex interplay of viscous, capillary, inertial and wettability effects. The lattice Boltzmann (LB) method is an accurate and efficient technique for modeling and simulating multiphase/multicomponent flows especially in complex flow configurations and media. In this presentation we present numerical simulation results of displacement process in thin long channels. The results are based on a new psuedo-potential multicomponent LB model with multiple relaxation time collision (MRT) model and explicit forcing scheme. We demonstrate that the proposed model is capable of accurately simulating the displacement process involving fluids with a wider range of viscosity ratios (>100) and which also leads to viscosity-independent interfacial tension and reduction of some important numerical artifacts.

The Importance of Lower Mantle Structure to Plate Stresses and Plate Motions

NASA Astrophysics Data System (ADS)

Holt, W. E.; Wang, X.; Ghosh, A.

2016-12-01

Plate motions and plate stresses are widely assumed as the surface expression of mantle convection. The generation of plate tectonics from mantle convection has been studied for many years. Lithospheric thickening (or ridge push) and slab pull forces are commonly accepted as the major driving forces for the plate motions. However, the importance of the lower mantle to plate stresses and plate motions remains less clear. Here, we use the joint modeling of lithosphere and mantle dynamics approach of Wang et al. (2015) to compute the tractions originating from deeper mantle convection and follow the method of Ghosh et al. (2013) to calculate gravitational potential energy per unit area (GPE) based on Crust 1.0 (Laske et al., 2013). Absolute values of deviatoric stresses are determined by the body force distributions (GPE gradients and traction magnitudes applied at the base of the lithosphere). We use the same relative viscosity model that Ghosh et al. (2013) used, and we solve for one single adjustable scaling factor that multiplies the entire relative viscosity field to provide absolute values of viscosity throughout the lithosphere. This distribution of absolute values of lithosphere viscosities defines the magnitudes of surface motions. In this procedure, the dynamic model first satisfies the internal constraint of no-net-rotation of motions. The model viscosity field is then scaled by the single factor until we achieve a root mean square (RMS) minimum between computed surface motions and the kinematic no-net-rotation (NNR) model of Kreemer et al. (2006). We compute plate stresses and plate motions from recently published global tomography models (over 70 based on Wang et al., 2015). We find that RMS misfits are significantly reduced when details of lower mantle structure from the latest tomography models are added to models that contain only upper and mid-mantle density distributions. One of the key reasons is that active upwelling from the Large Low Shear Velocity Provinces (LLSVPs) in the lower mantle in Pacific (Frost and Rost, 2014) provides important components of mantle flow affecting plate stresses and motions. We demonstrate in this paper how lower mantle density heterogeneity has a marked influence on plate stresses and plate motions.

Titin Based Viscosity in Ventricular Physiology: An Integrative Investigation of PEVK-Actin Interactions

PubMed Central

Chung, Charles S; Methawasin, Methajit; Nelson, O Lynne; Radke, Michael H; Hidalgo, Carlos G; Gotthardt, Michael; Granzier, Henk L

2011-01-01

Viscosity is proposed to modulate diastolic function, but only limited understanding of the source(s) of viscosity exists. In-vitro experiments have shown that the proline-glutamic acid-valine-lysine (PEVK) rich element of titin interacts with actin, causing a viscous force in the sarcomere. It is unknown whether this mechanism contributes to viscosity in-vivo. We tested the hypothesis that PEVK-actin interaction causes cardiac viscosity and is important in-vivo via an integrative physiological study on a unique PEVK-knockout (KO) model. Both skinned cardiomyocytes and papillary muscle fibers were isolated from wildtype (WT) and PEVK KO mice and passive viscosity was examined using stretch-hold-release and sinusoidal analysis. Viscosity was reduced by ~60% in KO myocytes and ~50% in muscle fibers at room temperature. The PEVK-actin interaction was not modulated by temperature or diastolic calcium, but was increased by lattice compression. Stretch-hold and sinusoidal frequency protocols on intact isolated mouse hearts showed a smaller, 30–40% reduction in viscosity, possibly due to actomyosin interactions, and showed that microtubules did not contribute to viscosity. Transmitral Doppler echocardiography similarly revealed a 40% decrease in LV chamber viscosity in the PEVK KO in-vivo. This integrative study is the first to quantify the influence of a specific molecular (PEVK-actin) viscosity in-vivo and shows that PEVK-actin interactions are an important physiological source of viscosity. PMID:21708170

One-dimensional wave bottom boundary layer model comparison: specific eddy viscosity and turbulence closure models

USGS Publications Warehouse

Puleo, J.A.; Mouraenko, O.; Hanes, D.M.

2004-01-01

Six one-dimensional-vertical wave bottom boundary layer models are analyzed based on different methods for estimating the turbulent eddy viscosity: Laminar, linear, parabolic, k—one equation turbulence closure, k−ε—two equation turbulence closure, and k−ω—two equation turbulence closure. Resultant velocity profiles, bed shear stresses, and turbulent kinetic energy are compared to laboratory data of oscillatory flow over smooth and rough beds. Bed shear stress estimates for the smooth bed case were most closely predicted by the k−ω model. Normalized errors between model predictions and measurements of velocity profiles over the entire computational domain collected at 15° intervals for one-half a wave cycle show that overall the linear model was most accurate. The least accurate were the laminar and k−ε models. Normalized errors between model predictions and turbulence kinetic energy profiles showed that the k−ω model was most accurate. Based on these findings, when the smallest overall velocity profile prediction error is required, the processing requirements and error analysis suggest that the linear eddy viscosity model is adequate. However, if accurate estimates of bed shear stress and TKE are required then, of the models tested, the k−ω model should be used.

A New Ductility Exhaustion Model for High Temperature Low Cycle Fatigue Life Prediction of Turbine Disk Alloys

NASA Astrophysics Data System (ADS)

Zhu, Shun-Peng; Huang, Hong-Zhong; Li, Haiqing; Sun, Rui; Zuo, Ming J.

2011-06-01

Based on ductility exhaustion theory and the generalized energy-based damage parameter, a new viscosity-based life prediction model is introduced to account for the mean strain/stress effects in the low cycle fatigue regime. The loading waveform parameters and cyclic hardening effects are also incorporated within this model. It is assumed that damage accrues by means of viscous flow and ductility consumption is only related to plastic strain and creep strain under high temperature low cycle fatigue conditions. In the developed model, dynamic viscosity is used to describe the flow behavior. This model provides a better prediction of Superalloy GH4133's fatigue behavior when compared to Goswami's ductility model and the generalized damage parameter. Under non-zero mean strain conditions, moreover, the proposed model provides more accurate predictions of Superalloy GH4133's fatigue behavior than that with zero mean strains.

Vertical structure of mean cross-shore currents across a barred surf zone

USGS Publications Warehouse

Haines, John W.; Sallenger, Asbury H.

1994-01-01

Mean cross-shore currents observed across a barred surf zone are compared to model predictions. The model is based on a simplified momentum balance with a turbulent boundary layer at the bed. Turbulent exchange is parameterized by an eddy viscosity formulation, with the eddy viscosity Aυ independent of time and the vertical coordinate. Mean currents result from gradients due to wave breaking and shoaling, and the presence of a mean setup of the free surface. Descriptions of the wave field are provided by the wave transformation model of Thornton and Guza [1983]. The wave transformation model adequately reproduces the observed wave heights across the surf zone. The mean current model successfully reproduces the observed cross-shore flows. Both observations and predictions show predominantly offshore flow with onshore flow restricted to a relatively thin surface layer. Successful application of the mean flow model requires an eddy viscosity which varies horizontally across the surf zone. Attempts are made to parameterize this variation with some success. The data does not discriminate between alternative parameterizations proposed. The overall variability in eddy viscosity suggested by the model fitting should be resolvable by field measurements of the turbulent stresses. Consistent shortcomings of the parameterizations, and the overall modeling effort, suggest avenues for further development and data collection.

Nonlinear ELM simulations based on a nonideal peeling–ballooning model using the BOUT++ code

DOE PAGES

Xu, X. Q.; Dudson, B. D.; Snyder, P. B.; ...

2011-09-23

A minimum set of equations based on the peeling–ballooning (P–B) model with nonideal physics effects (diamagnetic drift, E × B drift, resistivity and anomalous electron viscosity) is found to simulate pedestal collapse when using the BOUT++ simulation code, developed in part from the original fluid edge code BOUT. Linear simulations of P–B modes find good agreement in growth rate and mode structure with ELITE calculations. The influence of the E × B drift, diamagnetic drift, resistivity, anomalous electron viscosity, ion viscosity and parallel thermal diffusivity on P–B modes is being studied; we find that (1) the diamagnetic drift and Emore » × B drift stabilize the P–B mode in a manner consistent with theoretical expectations; (2) resistivity destabilizes the P–B mode, leading to resistive P–B mode; (3) anomalous electron and parallel ion viscosities destabilize the P–B mode, leading to a viscous P–B mode; (4) perpendicular ion viscosity and parallel thermal diffusivity stabilize the P–B mode. With addition of the anomalous electron viscosity under the assumption that the anomalous kinematic electron viscosity is comparable to the anomalous electron perpendicular thermal diffusivity, or the Prandtl number is close to unity, it is found from nonlinear simulations using a realistic high Lundquist number that the pedestal collapse is limited to the edge region and the ELM size is about 5–10% of the pedestal stored energy. Furthermore, this is consistent with many observations of large ELMs. The estimated island size is consistent with the size of fast pedestal pressure collapse. In the stable α-zones of ideal P–B modes, nonlinear simulations of viscous ballooning modes or current-diffusive ballooning mode (CDBM) for ITER H-mode scenarios are presented.« less

Numerical dissipation vs. subgrid-scale modelling for large eddy simulation

NASA Astrophysics Data System (ADS)

Dairay, Thibault; Lamballais, Eric; Laizet, Sylvain; Vassilicos, John Christos

2017-05-01

This study presents an alternative way to perform large eddy simulation based on a targeted numerical dissipation introduced by the discretization of the viscous term. It is shown that this regularisation technique is equivalent to the use of spectral vanishing viscosity. The flexibility of the method ensures high-order accuracy while controlling the level and spectral features of this purely numerical viscosity. A Pao-like spectral closure based on physical arguments is used to scale this numerical viscosity a priori. It is shown that this way of approaching large eddy simulation is more efficient and accurate than the use of the very popular Smagorinsky model in standard as well as in dynamic version. The main strength of being able to correctly calibrate numerical dissipation is the possibility to regularise the solution at the mesh scale. Thanks to this property, it is shown that the solution can be seen as numerically converged. Conversely, the two versions of the Smagorinsky model are found unable to ensure regularisation while showing a strong sensitivity to numerical errors. The originality of the present approach is that it can be viewed as implicit large eddy simulation, in the sense that the numerical error is the source of artificial dissipation, but also as explicit subgrid-scale modelling, because of the equivalence with spectral viscosity prescribed on a physical basis.

Carbohydrate composition, viscosity, solubility, and sensory acceptance of sweetpotato- and maize-based complementary foods

PubMed Central

Amagloh, Francis Kweku; Mutukumira, Anthony N.; Brough, Louise; Weber, Janet L.; Hardacre, Allan; Coad, Jane

2013-01-01

Background Cereal-based complementary foods from non-malted ingredients form a relatively high viscous porridge. Therefore, excessive dilution, usually with water, is required to reduce the viscosity to be appropriate for infant feeding. The dilution invariably leads to energy and nutrient thinning, that is, the reduction of energy and nutrient densities. Carbohydrate is the major constituent of food that significantly influences viscosity when heated in water. Objectives To compare the sweetpotato-based complementary foods (extrusion-cooked ComFa, roller-dried ComFa, and oven-toasted ComFa) and enriched Weanimix (maize-based formulation) regarding their 1) carbohydrate composition, 2) viscosity and water solubility index (WSI), and 3) sensory acceptance evaluated by sub-Sahara African women as model caregivers. Methods The level of simple sugars/carbohydrates was analysed by spectrophotometry, total dietary fibre by enzymatic-gravimetric method, and total carbohydrate and starch levels estimated by calculation. A Rapid Visco™ Analyser was used to measure viscosity. WSI was determined gravimetrically. A consumer sensory evaluation was used to evaluate the product acceptance of the roller-dried ComFa, oven-toasted ComFa, and enriched Weanimix. Results The sweetpotato-based complementary foods were, on average, significantly higher in maltose, sucrose, free glucose and fructose, and total dietary fibre, but they were markedly lower in starch content compared with the levels in the enriched Weanimix. Consequently, the sweetpotato-based complementary foods had relatively low apparent viscosity, and high WSI, than that of enriched Weanimix. The scores of sensory liking given by the caregivers were highest for the roller-dried ComFa, followed by the oven-toasted ComFa, and, finally, the enriched Weanimix. Conclusion The sweetpotato-based formulations have significant advantages as complementary food due to the high level of endogenous sugars and low starch content that reduce the viscosity, increase the solubility, impart desirable sensory characteristics, and potentially avoid excessive energy and nutrient thinning. PMID:23516115

Perfluoroalkyl phosphonic and phosphinic acids as proton conductors for anhydrous proton-exchange membranes.

PubMed

Herath, Mahesha B; Creager, Stephen E; Kitaygorodskiy, Alex; DesMarteau, Darryl D

2010-09-10

A study of proton-transport rates and mechanisms under anhydrous conditions using a series of acid model compounds, analogous to comb-branch perfluorinated ionomers functionalized with phosphonic, phosphinic, sulfonic, and carboxylic acid protogenic groups, is reported. Model compounds are characterized with respect to proton conductivity, viscosity, proton, and anion (conjugate base) self-diffusion coefficients, and Hammett acidity. The highest conductivities, and also the highest viscosities, are observed for the phosphonic and phosphinic acid model compounds. Arrhenius analysis of conductivity and viscosity for these two acids reveals much lower activation energies for ion transport than for viscous flow. Additionally, the proton self-diffusion coefficients are much higher than the conjugate-base self-diffusion coefficients for these two acids. Taken together, these data suggest that anhydrous proton transport in the phosphonic and phosphinic acid model compounds occurs primarily by a structure-diffusion, hopping-based mechanism rather than a vehicle mechanism. Further analysis of ionic conductivity and ion self-diffusion rates by using the Nernst-Einstein equation reveals that the phosphonic and phosphinic acid model compounds are relatively highly dissociated even under anhydrous conditions. In contrast, sulfonic and carboxylic acid-based systems exhibit relatively low degrees of dissociation under anhydrous conditions. These findings suggest that fluoroalkyl phosphonic and phosphinic acids are good candidates for further development as anhydrous, high-temperature proton conductors.

What's Cooler Than Being Cool? Ice-Sheet Models Using a Fluidity-Based FOSLS Approach to Nonlinear-Stokes Flow

NASA Astrophysics Data System (ADS)

Allen, Jeffery M.

This research involves a few First-Order System Least Squares (FOSLS) formulations of a nonlinear-Stokes flow model for ice sheets. In Glen's flow law, a commonly used constitutive equation for ice rheology, the viscosity becomes infinite as the velocity gradients approach zero. This typically occurs near the ice surface or where there is basal sliding. The computational difficulties associated with the infinite viscosity are often overcome by an arbitrary modification of Glen's law that bounds the maximum viscosity. The FOSLS formulations developed in this thesis are designed to overcome this difficulty. The first FOSLS formulation is just the first-order representation of the standard nonlinear, full-Stokes and is known as the viscosity formulation and suffers from the problem above. To overcome the problem of infinite viscosity, two new formulation exploit the fact that the deviatoric stress, the product of viscosity and strain-rate, approaches zero as the viscosity goes to infinity. Using the deviatoric stress as the basis for a first-order system results in the the basic fluidity system. Augmenting the basic fluidity system with a curl-type equation results in the augmented fluidity system, which is more amenable to the iterative solver, Algebraic MultiGrid (AMG). A Nested Iteration (NI) Newton-FOSLS-AMG approach is used to solve the nonlinear-Stokes problems. Several test problems from the ISMIP set of benchmarks is examined to test the effectiveness of the various formulations. These test show that the viscosity based method is more expensive and less accurate. The basic fluidity system shows optimal finite-element convergence. However, there is not yet an efficient iterative solver for this type of system and this is the topic of future research. Alternatively, AMG performs better on the augmented fluidity system when using specific scaling. Unfortunately, this scaling results in reduced finite-element convergence.

Bingham fluid behavior of plagioclase-bearing basaltic magma: Approach from laboratory viscosity measurements

NASA Astrophysics Data System (ADS)

Ishibashi, H.; Sato, H.

2010-12-01

Datasets of one atmosphere high temperature rotational viscometry of the Fuji 1707 basalt (Ishibashi, 2009) were analyzed based on the Bingham fluid model, and both yield stress and Bingham viscosity were determined. Reproducibility of the dataset by the Bingham fluid model was slightly better than that by the power law fluid modes adopted in our previous study although both the fluid models well represent the dataset in practical perspective. The relation between Bingham viscosity and crystallinity was compared with the Krieger-Dougherty equation, and both the maximum packing fraction of crystals and intrinsic viscosity for Bingham viscosity were determined ca. 0.45 and ca. 5.25, respectively, revealing that the maximum packing fraction decreased and intrinsic viscosity increased concomitantly with the increase in shape-anisotropy of crystals. However, the obtained value of the product of the maximum packing fraction and intrinsic viscosity (= ca. 2.36) was similar to that of uniform, isotropic-shaped particles (= 2.5), indicating that the effect of crystal shape-anisotropy on Bingham viscosity might be predicted only by change of the maximum packing fraction. Finite yield stress was detected for crystallinity larger than 0.133; it increased with crystallinity which suggests that critical crystallinity for onset of yield stress is at least lower than 0.133. The upper limit value of the critical crystallinity resembles the value calculated numerically for randomly oriented uniform particles by Saar et al. (2001) (0.10-0.15 for width/length ratio of 0.1-0.2, which is similar to the ratios in the basalt) whereas crystals in the basalt were moderately parallel arranged and their sizes vary significantly. That fact might be explained as follows; effects of parallel arrangement and size variation of crystals on the critical crystallinity are offset by the effect of variation in crystal shape-anisotropy, which suggests that shape-anisotropy distribution of crystals must be a critical factor for the onset of yield stress. Keywords: magma, viscosity, Bingham fluid, yield stress, plagioclase

Electroosmotic flow in capillary channels filled with nonconstant viscosity electrolytes: exact solution of the Navier-Stokes equation.

PubMed

Otevrel, Marek; Klepárník, Karel

2002-10-01

The partial differential equation describing unsteady velocity profile of electroosmotic flow (EOF) in a cylindrical capillary filled with a nonconstant viscosity electrolyte was derived. Analytical solution, based on the general Navier-Stokes equation, was found for constant viscosity electrolytes using the separation of variables (Fourier method). For the case of a nonconstant viscosity electrolyte, the steady-state velocity profile was calculated assuming that the viscosity decreases exponentially in the direction from the wall to the capillary center. Since the respective equations with nonconstant viscosity term are not solvable in general, the method of continuous binding conditions was used to solve this problem. In this method, an arbitrary viscosity profile can be modeled. The theoretical conclusions show that the relaxation times at which an EOF approaches the steady state are too short to have an impact on a separation process in any real systems. A viscous layer at the wall affects EOF significantly, if it is thicker than the Debye length of the electric double layer. The presented description of the EOF dynamics is applicable to any microfluidic systems.

An experimental study on rheological behavior of a nanofluid containing oxide nanoparticle and proposing a new correlation

NASA Astrophysics Data System (ADS)

Saeedi, Amir Hussein; Akbari, Mohammad; Toghraie, Davood

2018-05-01

In this paper, the nanofluid dynamic viscosity composed of CeO2- Ethylene Glycol is examined within 25-50 °C with 5 °C intervals and at six volume fractions (0.05, 0.1, 0.2, 0.4, 0.8 and 1.2%) experimentally. The nanofluid was exposed to ultrasound waves for various durations to study the effect of this parameter on dynamic viscosity of the fluid. We found that at a constant temperature, nanofluid viscosity increases with increases in the volume fraction of the nanoparticles. Also, at a given volume fraction, nanofluid viscosity decreases when temperature is increased. Maximum increase in nanofluid viscosity compared to the base fluid viscosity occurs at 25 °C and volume fraction of 1.2%. It can be inferred that the obtained mathematical relationship is a suitable predicting model for estimating dynamic viscosity of CeO2- Ethylene Glycol (EG) at different volume fractions and temperatures and its results are consistent to laboratory results in the set volume fraction and temperature ranges.

A Non-Arrhenian Viscosity Model for Natural Silicate Melts with Applications to Volcanology

NASA Astrophysics Data System (ADS)

Russell, J. K.; Giordano, D.; Dingwell, D. B.

2005-12-01

Silicate melt viscosity is the most important physical property in volcanic systems. It governs styles and rates of flow, velocity distributions in flowing magma, rates of vesiculation, and, ultimately, sets limits on coherent(vs. fragmented or disrupted) flow. The prediction of melt viscosity over the range of conditions found on terrestrial planets remains a challenge. However, the extraordinary increase in number and quality of published measurements of melt viscosity suggests the possibility of new models. Here we review the attributes of previous models for silicate melt viscosity and, then, present a new predictive model natural silicate melts. The importance of silicate melt viscosity was recognized early [1] and culminated in 2 models for predicting silicate melt viscosity [2,3]. These models used an Arrhenian T-dependence; they were limited by a limited experimental database dominated by high-T measurements. Subsequent models have aimed to: i) extend the compositional range of Arrhenian T-dependent models [4,5]; ii) to develop non-Arrhenian models for limited ranges of composition [6,7,8], iii) to develop new strategies for modelling the composition and T-dependence of viscosity [9,10,11], and, finally, to create chemical models for the non-Arrhenian T-dependence of natural melts [12]. We present a multicomponent model for the compositional and T dependence of silicate melt viscosity based on data spanning a wide range of anhydrous melt compositions. The experimental data include micropenetration and concentric cylinder viscometry measurements covering a viscosity range of 10-1 to 1012 Pa s and a T-range from 700 to 1650°C. These published data provide a high- quality database comprising ~ 800 experimental data on 44 well-characterized melt compositions. Our model uses the Adam-Gibbs equation to capture T-dependence: log η = A + B/[T · log (T/C)] where A, B, and C are adjustable parameters that vary for different melt compositions. We assume that all silicate melts converge to a common, but unknown, high-T limit (e.g., A) and that all compositional dependence is accommodated for by B and C. We adopt a linear compositional dependence for B and C: B = σi=1..n [xi βi] C = σi=1..n [xi γi] where xi's are the mole fractions of oxide components (n=8) and βi and γi are adjustable parameters. The model, therefore, comprises 2 · n+1 adjustable parameters which are optimized for against the experimental database including a common value of A and compositional coefficeints for B and C. The new model reproduces the original database to within experimental uncertainty and can predict the viscosity of silicate melts across the full range of conditions found in Nature. References Cited: [1] Friedman et al., 1963. J Geophys Res 68, 6523-6535. [2] Bottinga Y & Weill D 1972. Am J Sci 272, 438- 475. [3] Shaw HR 1972. Am J Sci 272, 438- 475. [4] Persikov ES 1991. Adv Phys Geochem 9, 1-40. [5] Prusevich AA 1988. Geol Geofiz 29, 67-69. [6] Baker DR 1996. Am Min 81, 126-134. [7] Hess KU & Dingwell DB 1996. Am Min 81, 1297- 1300. [8] Zhang, et al. 2003. Am min 88, 1741- 1752. [9] Russell et al. 2002. Eur J Min 14, 417-428. [10] Russell et al. 2003. Am Min 8, 1390- 1394. [11] Russell JK & Giordano D In Press. Geochim Cosmochim Acta. [12] Giordano D & Dingwell DB 2003. Earth Planet. Sci. Lett. 208, 337-349.

A Dynamic Eddy Viscosity Model for the Shallow Water Equations Solved by Spectral Element and Discontinuous Galerkin Methods

NASA Astrophysics Data System (ADS)

Marras, Simone; Suckale, Jenny; Giraldo, Francis X.; Constantinescu, Emil

2016-04-01

We present the solution of the viscous shallow water equations where viscosity is built as a residual-based subgrid scale model originally designed for large eddy simulation of compressible [1] and stratified flows [2]. The necessity of viscosity for a shallow water model not only finds motivation from mathematical analysis [3], but is supported by physical reasoning as can be seen by an analysis of the energetics of the solution. We simulated the flow of an idealized wave as it hits a set of obstacles. The kinetic energy spectrum of this flow shows that, although the inviscid Galerkin solutions -by spectral elements and discontinuous Galerkin [4]- preserve numerical stability in spite of the spurious oscillations in the proximity of the wave fronts, the slope of the energy cascade deviates from the theoretically expected values. We show that only a sufficiently small amount of dynamically adaptive viscosity removes the unwanted high-frequency modes while preserving the overall sharpness of the solution. In addition, it yields a physically plausible energy decay. This work is motivated by a larger interest in the application of a shallow water model to the solution of tsunami triggered coastal flows. In particular, coastal flows in regions around the world where coastal parks made of mitigation hills of different sizes and configurations are considered as a means to deviate the power of the incoming wave. References [1] M. Nazarov and J. Hoffman (2013) "Residual-based artificial viscosity for simulation of turbulent compressible flow using adaptive finite element methods" Int. J. Numer. Methods Fluids, 71:339-357 [2] S. Marras, M. Nazarov, F. X. Giraldo (2015) "Stabilized high-order Galerkin methods based on a parameter-free dynamic SGS model for LES" J. Comput. Phys. 301:77-101 [3] J. F. Gerbeau and B. Perthame (2001) "Derivation of the viscous Saint-Venant system for laminar shallow water; numerical validation" Discrete Contin. Dyn. Syst. Ser. B, 1:89?102 [4] F. X. Giraldo and M. Restelli (2010) "High-order semi-implicit time-integrators for a triangular discontinuous Galerkin oceanic shallow water model. Int. J. Numer. Methods Fluids, 63:1077-1102

Dynamics of Compressible Convection and Thermochemical Mantle Convection

NASA Astrophysics Data System (ADS)

Liu, Xi

The Earth's long-wavelength geoid anomalies have long been used to constrain the dynamics and viscosity structure of the mantle in an isochemical, whole-mantle convection model. However, there is strong evidence that the seismically observed large low shear velocity provinces (LLSVPs) in the lowermost mantle are chemically distinct and denser than the ambient mantle. In this thesis, I investigated how chemically distinct and dense piles influence the geoid. I formulated dynamically self-consistent 3D spherical convection models with realistic mantle viscosity structure which reproduce Earth's dominantly spherical harmonic degree-2 convection. The models revealed a compensation effect of the chemically dense LLSVPs. Next, I formulated instantaneous flow models based on seismic tomography to compute the geoid and constrain mantle viscosity assuming thermochemical convection with the compensation effect. Thermochemical models reconcile the geoid observations. The viscosity structure inverted for thermochemical models is nearly identical to that of whole-mantle models, and both prefer weak transition zone. Our results have implications for mineral physics, seismic tomographic studies, and mantle convection modelling. Another part of this thesis describes analyses of the influence of mantle compressibility on thermal convection in an isoviscous and compressible fluid with infinite Prandtl number. A new formulation of the propagator matrix method is implemented to compute the critical Rayleigh number and the corresponding eigenfunctions for compressible convection. Heat flux and thermal boundary layer properties are quantified in numerical models and scaling laws are developed.

Structure-Property Relationships of Architectural Coatings by Neutron Methods

NASA Astrophysics Data System (ADS)

Nakatani, Alan

2015-03-01

Architectural coatings formulations are multi-component mixtures containing latex polymer binder, pigment, rheology modifiers, surfactants, and colorants. In order to achieve the desired flow properties for these formulations, measures of the underlying structure of the components as a function of shear rate and the impact of formulation variables on the structure is necessary. We have conducted detailed measurements to understand the evolution under shear of local microstructure and larger scale mesostructure in model architectural coatings formulations by small angle neutron scattering (SANS) and ultra small angle neutron scattering (USANS), respectively. The SANS results show an adsorbed layer of rheology modifier molecules exist on the surface of the latex particles. However, the additional hydrodynamic volume occupied by the adsorbed surface layer is insufficient to account for the observed viscosity by standard hard sphere suspension models (Krieger-Dougherty). The USANS results show the presence of latex aggregates, which are fractal in nature. These fractal aggregates are the primary structures responsible for coatings formulation viscosity. Based on these results, a new model for the viscosity of coatings formulations has been developed, which is capable of reproducing the observed viscosity behavior.

Phase-field modeling of mixing/demixing of regular binary mixtures with a composition-dependent viscosity

NASA Astrophysics Data System (ADS)

Lamorgese, A.; Mauri, R.

2017-04-01

We simulate the mixing (demixing) process of a quiescent binary liquid mixture with a composition-dependent viscosity which is instantaneously brought from the two-phase (one-phase) to the one-phase (two-phase) region of its phase diagram. Our theoretical approach follows a standard diffuse-interface model of partially miscible regular binary mixtures wherein convection and diffusion are coupled via a nonequilibrium capillary force, expressing the tendency of the phase-separating system to minimize its free energy. Based on 2D simulation results, we discuss the influence of viscosity ratio on basic statistics of the mixing (segregation) process triggered by a rapid heating (quench), assuming that the ratio of capillary to viscous forces (a.k.a. the fluidity coefficient) is large. We show that, for a phase-separating system, at a fixed value of the fluidity coefficient (with the continuous phase viscosity taken as a reference), the separation depth and the characteristic length of single-phase microdomains decrease monotonically for increasing values of the viscosity of the dispersed phase. This variation, however, is quite small, in agreement with experimental results. On the other hand, as one might expect, at a fixed viscosity of the dispersed phase both of the above statistics increase monotonically as the viscosity of the continuous phase decreases. Finally, we show that for a mixing system the attainment of a single-phase equilibrium state by coalescence and diffusion is retarded by an increase in the viscosity ratio at a fixed fluidity for the dispersed phase. In fact, for large enough values of the viscosity ratio, a thin film of the continuous phase becomes apparent when two drops of the minority phase approach each other, which further retards coalescence.

Modeling the viscosity of polydisperse suspensions: Improvements in prediction of limiting behavior

NASA Astrophysics Data System (ADS)

Mwasame, Paul M.; Wagner, Norman J.; Beris, Antony N.

2016-06-01

The present study develops a fully consistent extension of the approach pioneered by Farris ["Prediction of the viscosity of multimodal suspensions from unimodal viscosity data," Trans. Soc. Rheol. 12, 281-301 (1968)] to describe the viscosity of polydisperse suspensions significantly improving upon our previous model [P. M. Mwasame, N. J. Wagner, and A. N. Beris, "Modeling the effects of polydispersity on the viscosity of noncolloidal hard sphere suspensions," J. Rheol. 60, 225-240 (2016)]. The new model captures the Farris limit of large size differences between consecutive particle size classes in a suspension. Moreover, the new model includes a further generalization that enables its application to real, complex suspensions that deviate from ideal non-colloidal suspension behavior. The capability of the new model to predict the viscosity of complex suspensions is illustrated by comparison against experimental data.

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

NASA Astrophysics Data System (ADS)

Reckinger, Shanon M.; Petersen, Mark R.; Reckinger, Scott J.

2015-12-01

MPAS-Ocean is used to simulate an idealized, density-driven overflow using the dynamics of overflow mixing and entrainment (DOME) setup. 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. Lastly, 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.

Characterization of oscillator circuits for monitoring the density-viscosity of liquids by means of piezoelectric MEMS microresonators

NASA Astrophysics Data System (ADS)

Toledo, J.; Ruiz-Díez, V.; Pfusterschmied, G.; Schmid, U.; Sánchez-Rojas, J. L.

2017-06-01

Real-time monitoring of the physical properties of liquids, such as lubricants, is a very important issue for the automotive industry. For example, contamination of lubricating oil by diesel soot has a significant impact on engine wear. Resonant microstructures are regarded as a precise and compact solution for tracking the viscosity and density of lubricant oils. In this work, we report a piezoelectric resonator, designed to resonate with the 4th order out-of-plane modal vibration, 15-mode, and the interface circuit and calibration process for the monitoring of oil dilution with diesel fuel. In order to determine the resonance parameters of interest, i.e. resonant frequency and quality factor, an interface circuit was implemented and included within a closed-loop scheme. Two types of oscillator circuits were tested, a Phase-Locked Loop based on instrumentation, and a more compact version based on discrete electronics, showing similar resolution. Another objective of this work is the assessment of a calibration method for piezoelectric MEMS resonators in simultaneous density and viscosity sensing. An advanced calibration model, based on a Taylor series of the hydrodynamic function, was established as a suitable method for determining the density and viscosity with the lowest calibration error. Our results demonstrate the performance of the resonator in different oil samples with viscosities up to 90 mPa•s. At the highest value, the quality factor measured at 25°C was around 22. The best resolution obtained was 2.4•10-6 g/ml for the density and 2.7•10-3 mPa•s for the viscosity, in pure lubricant oil SAE 0W30 at 90°C. Furthermore, the estimated density and viscosity values with the MEMS resonator were compared to those obtained with a commercial density-viscosity meter, reaching a mean calibration error in the best scenario of around 0.08% for the density and 3.8% for the viscosity.

Viscosity, density, and thermal conductivity of aluminum oxide and zinc oxide nanolubricants

PubMed Central

Kedzierski, M.A.; Brignoli, R.; Quine, K.T.; Brown, J.S.

2017-01-01

This paper presents liquid kinematic viscosity, density, and thermal conductivity measurements of eleven different synthetic polyolester-based nanoparticle nanolubricants (dispersions) at atmospheric pressure over the temperature range 288 K to 318 K. Aluminum oxide (Al2O3) and zinc oxide (ZnO) nanoparticles with nominal diameters of 127 nm and 135 nm, respectively, were investigated. A good dispersion of the spherical and non-spherical nanoparticles in the lubricant was maintained with a surfactant. Viscosity, density, and thermal conductivity measurements were made for the neat lubricant along with eleven nanolubricants with differing nanoparticle and surfactant mass fractions. Existing models were used to predict kinematic viscosity (±20%), thermal conductivity (±1%), and specific volume (±6%) of the nanolubricant as a function of temperature, nanoparticle mass fraction, surfactant mass fraction, and nanoparticle diameter. The liquid viscosity, density and thermal conductivity were shown to increase with respect to increasing nanoparticle mass fraction. PMID:28736463

Computational tool for the early screening of monoclonal antibodies for their viscosities

PubMed Central

Agrawal, Neeraj J; Helk, Bernhard; Kumar, Sandeep; Mody, Neil; Sathish, Hasige A.; Samra, Hardeep S.; Buck, Patrick M; Li, Li; Trout, Bernhardt L

2016-01-01

Highly concentrated antibody solutions often exhibit high viscosities, which present a number of challenges for antibody-drug development, manufacturing and administration. The antibody sequence is a key determinant for high viscosity of highly concentrated solutions; therefore, a sequence- or structure-based tool that can identify highly viscous antibodies from their sequence would be effective in ensuring that only antibodies with low viscosity progress to the development phase. Here, we present a spatial charge map (SCM) tool that can accurately identify highly viscous antibodies from their sequence alone (using homology modeling to determine the 3-dimensional structures). The SCM tool has been extensively validated at 3 different organizations, and has proved successful in correctly identifying highly viscous antibodies. As a quantitative tool, SCM is amenable to high-throughput automated analysis, and can be effectively implemented during the antibody screening or engineering phase for the selection of low-viscosity antibodies. PMID:26399600

The viscosity of magmatic silicate liquids: A model for calculation

NASA Technical Reports Server (NTRS)

Bottinga, Y.; Weill, D. F.

1971-01-01

A simple model has been designed to allow reasonably accurate calculations of viscosity as a function of temperature and composition. The problem of predicting viscosities of anhydrous silicate liquids has been investigated since such viscosity numbers are applicable to many extrusive melts and to nearly dry magmatic liquids in general. The fluidizing action of water dissolved in silicate melts is well recognized and it is now possible to predict the effect of water content on viscosity in a semiquantitative way. Water was not incorporated directly into the model. Viscosities of anhydrous compositions were calculated, and, where necessary, the effect of added water and estimated. The model can be easily modified to incorporate the effect of water whenever sufficient additional data are accumulated.

A new traffic model with a lane-changing viscosity term

NASA Astrophysics Data System (ADS)

Ko, Hung-Tang; Liu, Xiao-He; Guo, Ming-Min; Wu, Zheng

2015-09-01

In this paper, a new continuum traffic flow model is proposed, with a lane-changing source term in the continuity equation and a lane-changing viscosity term in the acceleration equation. Based on previous literature, the source term addresses the impact of speed difference and density difference between adjacent lanes, which provides better precision for free lane-changing simulation; the viscosity term turns lane-changing behavior to a “force” that may influence speed distribution. Using a flux-splitting scheme for the model discretization, two cases are investigated numerically. The case under a homogeneous initial condition shows that the numerical results by our model agree well with the analytical ones; the case with a small initial disturbance shows that our model can simulate the evolution of perturbation, including propagation, dissipation, cluster effect and stop-and-go phenomenon. Project supported by the National Natural Science Foundation of China (Grant Nos. 11002035 and 11372147) and Hui-Chun Chin and Tsung-Dao Lee Chinese Undergraduate Research Endowment (Grant No. CURE 14024).

Internal friction in enzyme reactions.

PubMed

Rauscher, Anna; Derényi, Imre; Gráf, László; Málnási-Csizmadia, András

2013-01-01

The empirical concept of internal friction was introduced 20 years ago. This review summarizes the results of experimental and theoretical studies that help to uncover the nature of internal friction. After the history of the concept, we describe the experimental challenges in measuring and interpreting internal friction based on the viscosity dependence of enzyme reactions. We also present speculations about the structural background of this viscosity dependence. Finally, some models about the relationship between the energy landscape and internal friction are outlined. Alternative concepts regarding the viscosity dependence of enzyme reactions are also discussed. Copyright © 2012 International Union of Biochemistry and Molecular Biology, Inc.

Viscosity Meaurement Technique for Metal Fuels

DOE Office of Scientific and Technical Information (OSTI.GOV)

Ban, Heng; Kennedy, Rory

2015-02-09

Metallic fuels have exceptional transient behavior, excellent thermal conductivity, and a more straightforward reprocessing path, which does not separate out pure plutonium from the process stream. Fabrication of fuel containing minor actinides and rare earth (RE) elements for irradiation tests, for instance, U-20Pu-3Am-2Np-1.0RE-15Zr samples at the Idaho National Laboratory, is generally done by melt casting in an inert atmosphere. For the design of a casting system and further scale up development, computational modeling of the casting process is needed to provide information on melt flow and solidification for process optimization. Therefore, there is a need for melt viscosity data, themore » most important melt property that controls the melt flow. The goal of the project was to develop a measurement technique that uses fully sealed melt sample with no Americium vapor loss to determine the viscosity of metallic melts and at temperatures relevant to the casting process. The specific objectives of the project were to: develop mathematical models to establish the principle of the measurement method, design and build a viscosity measurement prototype system based on the established principle, and calibrate the system and quantify the uncertainty range. The result of the project indicates that the oscillation cup technique is applicable for melt viscosity measurement. Detailed mathematical models of innovative sample ampoule designs were developed to not only determine melt viscosity, but also melt density under certain designs. Measurement uncertainties were analyzed and quantified. The result of this project can be used as the initial step toward the eventual goal of establishing a viscosity measurement system for radioactive melts.« less

An experimental approach to manufacturing technology of historical glass (XIII-XV centuries). Comparison with current glassmaking technology.

NASA Astrophysics Data System (ADS)

Tarrago, Mariona; Gimeno, Domingo; Bazzocchi, Flavia; Garcia-Valles, Maite; Martinez, Salvador

2015-04-01

One of the major and less explored issues in the characterization of historical glasses is the determination of their viscosity as a function of temperature in order to constrain technological aspects of glass production. Until now, assumptions on temperatures have been based on mathematical models based on chemical compositions. Hence, the topic of this work is to explore the technology of stained glass production related to the workability and melting process of the glass by experimental laboratory measurements. This work presents the analysis of viscosity of glasses from different historical sites and chemical compositions: four from Santes Creus (Tarragona, XIII century), two of classical medieval stained glass window from Santa Maria de Pedralbes (Barcelona, mid XIV century), and three of evolved late-medieval type from Santa Maria del Mar (Barcelona first half of XV century), and one sample of soda-lime industrial glass by means of Hot-Stage Microscopy and glass transformation temperature Tg by dilatometry. These data are then compared to the predictions on theoretical viscosity obtained from mathematical models based on chemical composition. These samples are classified according to their major modifier in: Na-rich (12-17% of Na2O, between 65-77% of SiO2 and less than 3 % of K2O); Ca-rich (29% of CaO, 54% of SiO2, 4% of K2O, and 4% of Na2O); K-Ca-rich (17 to 21% of K2O, more than 14% of CaO, 49-55% of SiO2and less than 2% of Na2O); Na-Ca-rich (12-14% of Na2O, 9-15% of CaO, 57-71% of SiO2 and < 6% of K2O). Glass transition temperature (Tg) is correlated to chemical composition: 464-492 °C for Na-rich, 645 °C for Ca-rich, 582-586 °C for K-Ca-rich and 497-542 °C for Na-Ca-rich glasses. Experimental viscosity-temperature curves are traced using Tg and fixed viscosity points measured by Hot-Stage microscopy (according to German standard 51730) in order to provide more accurate insight into the phases of glass production process (melting, working, conditioning and annealing ranges). These results are also compared to mathematical models of glass viscosity based on chemical composition. The annealing range (viscosity between 1013.5 and 1012 Pa-s) is reached at temperatures between 484-633°C (strain point) and 509-664°C (upper limit). The working point (viscosity of 103 Pa-s) has temperature values in the range between 958 and 1097°C.

Nonlinear simulations of peeling-ballooning modes with anomalous electron viscosity and their role in edge localized mode crashes

DOE PAGES

Xu, X. Q.; Dudson, B.; Snyder, P. B.; ...

2010-10-22

A minimum set of equations based on the peeling-ballooning (P-B) model with nonideal physics effects (diamagnetic drift, E×B drift, resistivity, and anomalous electron viscosity) is found to simulate pedestal collapse when using the new BOUT++ simulation code, developed in part from the original fluid edge code BOUT. Nonlinear simulations of P-B modes demonstrate that the P-B modes trigger magnetic reconnection, which leads to the pedestal collapse. With the addition of a model of the anomalous electron viscosity under the assumption that the electron viscosity is comparable to the anomalous electron thermal diffusivity, it is found from simulations using a realisticmore » high-Lundquist number that the pedestal collapse is limited to the edge region and the edge localized mode (ELM) size is about 5–10% of the pedestal stored energy. Furthermore, this is consistent with many observations of large ELMs.« less

The effects of rotational flow, viscosity, thickness, and shape on transonic flutter dip phenomena

NASA Technical Reports Server (NTRS)

Reddy, T. S. R.; Srivastava, Rakesh; Kaza, Krishna Rao V.

1988-01-01

The transonic flutter dip phenomena on thin airfoils, which are employed for propfan blades, is investigated using an integrated Euler/Navier-Stokes code and a two degrees of freedom typical section structural model. As a part of the code validation, the flutter characteristics of the NACA 64A010 airfoil are also investigated. In addition, the effects of artificial dissipation models, rotational flow, initial conditions, mean angle of attack, viscosity, airfoil thickness and shape on flutter are investigated. The results obtained with a Euler code for the NACA 64A010 airfoil are in reasonable agreement with published results obtained by using transonic small disturbance and Euler codes. The two artificial dissipation models, one based on the local pressure gradient scaled by a common factor and the other based on the local pressure gradient scaled by a spectral radius, predicted the same flutter speeds except in the recovery region for the case studied. The effects of rotational flow, initial conditions, mean angle of attack, and viscosity for the Reynold's number studied seem to be negligible or small on the minima of the flutter dip.

Modeling and control of plasma rotation for NSTX using neoclassical toroidal viscosity and neutral beam injection

NASA Astrophysics Data System (ADS)

Goumiri, I. R.; Rowley, C. W.; Sabbagh, S. A.; Gates, D. A.; Gerhardt, S. P.; Boyer, M. D.; Andre, R.; Kolemen, E.; Taira, K.

2016-03-01

A model-based feedback system is presented to control plasma rotation in a magnetically confined toroidal fusion device, to maintain plasma stability for long-pulse operation. This research uses experimental measurements from the National Spherical Torus Experiment (NSTX) and is aimed at controlling plasma rotation using two different types of actuation: momentum from injected neutral beams and neoclassical toroidal viscosity generated by three-dimensional applied magnetic fields. Based on the data-driven model obtained, a feedback controller is designed, and predictive simulations using the TRANSP plasma transport code show that the controller is able to attain desired plasma rotation profiles given practical constraints on the actuators and the available measurements of rotation.

Modeling and control of plasma rotation for NSTX using neoclassical toroidal viscosity and neutral beam injection

DOE Office of Scientific and Technical Information (OSTI.GOV)

Goumiri, I. R.; Rowley, C. W.; Sabbagh, S. A.

2016-02-19

A model-based feedback system is presented to control plasma rotation in a magnetically confined toroidal fusion device, to maintain plasma stability for long-pulse operation. This research uses experimental measurements from the National Spherical Torus Experiment (NSTX) and is aimed at controlling plasma rotation using two different types of actuation: momentum from injected neutral beams and neoclassical toroidal viscosity generated by three-dimensional applied magnetic fields. Based on the data-driven model obtained, a feedback controller is designed, and predictive simulations using the TRANSP plasma transport code show that the controller is able to attain desired plasma rotation profiles given practical constraints onmore » the actuators and the available measurements of rotation.« less

Non-steady peristaltic propulsion with exponential variable viscosity: a study of transport through the digestive system.

PubMed

Tripathi, Dharmendra; Pandey, S K; Siddiqui, Abdul; Bég, O Anwar

2014-01-01

A theoretical study is presented for transient peristaltic flow of an incompressible fluid with variable viscosity in a finite length cylindrical tube as a simulation of transport in physiological vessels and biomimetic peristaltic pumps. The current axisymmetric analysis is qualitatively similar to two-dimensional analysis but exhibits quantitative variations. The current analysis is motivated towards further elucidating the physiological migration of gastric suspensions (food bolus) in the human digestive system. It also applies to variable viscosity industrial fluid (waste) peristaltic pumping systems. First, an axisymmetric model is analysed in the limit of large wavelength ([Formula: see text]) and low Reynolds number ([Formula: see text]) for axial velocity, radial velocity, pressure, hydromechanical efficiency and stream function in terms of radial vibration of the wall ([Formula: see text]), amplitude of the wave ([Formula: see text]), averaged flow rate ([Formula: see text]) and variable viscosity ([Formula: see text]). Subsequently, the peristaltic flow of a fluid with an exponential viscosity model is examined, which is based on the analytical solutions for pressure, wall shear stress, hydromechanical efficiency and streamline patterns in the finite length tube. The results are found to correlate well with earlier studies using a constant viscosity formulation. This study reveals some important features in the flow characteristics including the observation that pressure as well as both number and size of lower trapped bolus increases. Furthermore, the study indicates that hydromechanical efficiency reduces with increasing magnitude of viscosity parameter.

Experimental and Molecular Modeling Evaluation of the Physicochemical Properties of Proline-Based Deep Eutectic Solvents.

PubMed

van den Bruinhorst, Adriaan; Spyriouni, Theodora; Hill, Jörg-Rüdiger; Kroon, Maaike C

2018-01-11

The liquid range and applicability of deep eutectic solvents (DESs) are determined by their physicochemical properties. In this work, the physicochemical properties of glycolic acid:proline and malic acid:proline were evaluated experimentally and with MD simulations at five different ratios. Both DESs exhibited esterification upon preparation, which affected the viscosity in particular. In order to minimize oligomer formation and water release, three different experimental preparation methods were explored, but none could prevent esterification. The experimental and calculated densities of the DESs were found to be in good agreement. The measured and modeled glass transition temperature showed similar trends with composition, as did the experimental viscosity and the calculated diffusivities. The MD simulations provided additional insight at the atomistic level, showing that at acid-rich compositions, the acid-acid hydrogen bonding (HB) interactions prevail. Malic acid-based DESs show stronger acid-acid HB interactions than glycolic acid-based ones, possibly explaining its extreme viscosity. Upon the addition of proline, the interspecies interactions become predominant, confirming the formation of the widely assumed HB network between the DESs constituents in the liquid phase.

Effect of a weak layer at the base of an oceanic plate on subduction dynamics

NASA Astrophysics Data System (ADS)

Carluccio, Roberta; Kaus, Boris

2017-04-01

The plate tectonics model relies on the concept of a relatively rigid lithospheric lid moving over a weaker asthenosphere. In this frame, the lithosphere asthenosphere boundary (LAB) is a first-order discontinuity that accommodates differential motions between tectonic plates and the underlying mantle. Recent seismic studies have revealed the existence of a low velocity and high electrical conductivity layer at the base of subducting tectonic plates. This thin layer has been interpreted as being weak and slightly buoyant and was suggested to affect the dynamics of subducting plates. However, geodynamically, the role of a weak layer at the base of the lithosphere remains poorly studied, especially at subduction zones. Therefore, we here use numerical models to investigate the first-order effects of a weak buoyant layer at the LAB on subduction dynamics. We employ both 2-D and 3-D models in which the slab and mantle are either linear viscous or have a more realistic temperature-dependent visco-elastic-plastic rheology. Results show that a weak layer affects the dynamics of the plates, foremost by increasing the subduction speed. The impact of this effect depends on the thickness of the layer and the viscosity contrast between the mantle and the weak layer. For moderate viscosity contrasts (<100) and a layer thickness of 30 km, it increases the plate velocity but not the overall shape of the slab. However, for larger viscosity contrasts (>1000), it can also change the morphology of the subduction itself, perhaps because this changes the overall effective viscosity contrast between the slab the and the mantle. For thinner layers, the overall effect is reduced. Yet, if seismological observations are correct that suggests that this layer is 10 km thick and partially molten, such that the viscosity is 1000 times lower than that of the mantle, our models suggest that this effect should be measurable. Some of our models also show a pile-up of weak material in the bending zone of the subducting plate, consistent with recent seismological observations.

Geothermal Heat Flux and Upper Mantle Viscosity across West Antarctica: Insights from the UKANET and POLENET Seismic Networks

NASA Astrophysics Data System (ADS)

O'Donnell, J. P.; Dunham, C.; Stuart, G. W.; Brisbourne, A.; Nield, G. A.; Whitehouse, P. L.; Hooper, A. J.; Nyblade, A.; Wiens, D.; Aster, R. C.; Anandakrishnan, S.; Huerta, A. D.; Wilson, T. J.; Winberry, J. P.

2017-12-01

Quantifying the geothermal heat flux at the base of ice sheets is necessary to understand their dynamics and evolution. The heat flux is a composite function of concentration of upper crustal radiogenic elements and flow of heat from the mantle into the crust. Radiogenic element concentration varies with tectonothermal age, while heat flow across the crust-mantle boundary depends on crustal and lithospheric thicknesses. Meanwhile, accurately monitoring current ice mass loss via satellite gravimetry or altimetry hinges on knowing the upper mantle viscosity structure needed to account for the superimposed glacial isostatic adjustment (GIA) signal in the satellite data. In early 2016 the UK Antarctic Network (UKANET) of 10 broadband seismometers was deployed for two years across the southern Antarctic Peninsula and Ellsworth Land. Using UKANET data in conjunction with seismic records from our partner US Polar Earth Observing Network (POLENET) and the Antarctic Seismographic Argentinian Italian Network (ASAIN), we have developed a 3D shear wave velocity model of the West Antarctic crust and uppermost mantle based on Rayleigh and Love wave phase velocity dispersion curves extracted from ambient noise cross-correlograms. We combine seismic receiver functions with the shear wave model to help constrain the depth to the crust-mantle boundary across West Antarctica and delineate tectonic domains. The shear wave model is subsequently converted to temperature using a database of densities and elastic properties of minerals common in crustal and mantle rocks, while the various tectonic domains are assigned upper crustal radiogenic element concentrations based on their inferred tectonothermal ages. We combine this information to map the basal geothermal heat flux variation across West Antarctica. Mantle viscosity depends on factors including temperature, grain size, the hydrogen content of olivine and the presence of melt. Using published mantle xenolith and magnetotelluric data to constrain grain size and hydrogen content, respectively, we use the temperature model to estimate the regional upper mantle viscosity structure. The viscosity information will be incorporated in a 3D GIA model that will better constrain estimates of current ice loss from the West Antarctic Ice Sheet.

Sensitivity of viscosity Arrhenius parameters to polarity of liquids

NASA Astrophysics Data System (ADS)

Kacem, R. B. H.; Alzamel, N. O.; Ouerfelli, N.

2017-09-01

Several empirical and semi-empirical equations have been proposed in the literature to estimate the liquid viscosity upon temperature. In this context, this paper aims to study the effect of polarity of liquids on the modeling of the viscosity-temperature dependence, considering particularly the Arrhenius type equations. To achieve this purpose, the solvents are classified into three groups: nonpolar, borderline polar and polar solvents. Based on adequate statistical tests, we found that there is strong evidence that the polarity of solvents affects significantly the distribution of the Arrhenius-type equation parameters and consequently the modeling of the viscosity-temperature dependence. Thus, specific estimated values of parameters for each group of liquids are proposed in this paper. In addition, the comparison of the accuracy of approximation with and without classification of liquids, using the Wilcoxon signed-rank test, shows a significant discrepancy of the borderline polar solvents. For that, we suggested in this paper new specific coefficient values of the simplified Arrhenius-type equation for better estimation accuracy. This result is important given that the accuracy in the estimation of the viscosity-temperature dependence may affect considerably the design and the optimization of several industrial processes.

Predicting the viscosity of solids using steady-state creep behavior of the fibrous composites semi-theoretically

NASA Astrophysics Data System (ADS)

Monfared, Vahid

A semi-analytical formulation is presented for obtaining the viscosity of solids (such as metals) using the steady state creep model of the short-fiber composites. For achieving this aim, fluid mechanics theory is used for determining the viscosity. Sometimes, obtaining the viscosity is experimentally difficult and intricate. So, the present model may be beneficial to obtain the viscosity of metals.

Temperature-Dependent Conformations of Model Viscosity Index Improvers

DOE Office of Scientific and Technical Information (OSTI.GOV)

Ramasamy, Uma Shantini; Cosimbescu, Lelia; Martini, Ashlie

2015-05-01

Lubricants are comprised of base oils and additives where additives are chemicals that are deliberately added to the oil to enhance properties and inhibit degradation of the base oils. Viscosity index (VI) improvers are an important class of additives that reduce the decline of fluid viscosity with temperature [1], enabling optimum lubricant performance over a wider range of operating temperatures. These additives are typically high molecular weight polymers, such as, but not limited to, polyisobutylenes, olefin copolymer, and polyalkylmethacrylates, that are added in concentrations of 2-5% (w/w). Appropriate polymers, when dissolved in base oil, expand from a coiled to anmore » uncoiled state with increasing temperature [2]. The ability of VI additives to increase their molar volume and improve the temperature-viscosity dependence of lubricants suggests there is a strong relationship between molecular structure and additive functionality [3]. In this work, we aim to quantify the changes in polymer size with temperature for four polyisobutylene (PIB) based molecular structures at the nano-scale using molecular simulation tools. As expected, the results show that the polymers adopt more conformations at higher temperatures, and there is a clear indication that the expandability of a polymer is strongly influenced by molecular structure.« less

Structural and dynamical properties of liquid Al-Au alloys

NASA Astrophysics Data System (ADS)

Peng, H. L.; Voigtmann, Th.; Kolland, G.; Kobatake, H.; Brillo, J.

2015-11-01

We investigate temperature- and composition-dependent structural and dynamical properties of Al-Au melts. Experiments are performed to obtain accurate density and viscosity data. The system shows a strong negative excess volume, similar to other Al-based binary alloys. We develop a molecular-dynamics (MD) model of the melt based on the embedded-atom method (EAM), gauged against the available experimental liquid-state data. A rescaling of previous EAM potentials for solid-state Au and Al improves the quantitative agreement with experimental data in the melt. In the MD simulation, the admixture of Au to Al can be interpreted as causing a local compression of the less dense Al system, driven by less soft Au-Au interactions. This local compression provides a microscopic mechanism explaining the strong negative excess volume of the melt. We further discuss the concentration dependence of self- and interdiffusion and viscosity in the MD model. Al atoms are more mobile than Au, and their increased mobility is linked to a lower viscosity of the melt.

Thermal evolution of flattening slab and formation of wet plume: Insight into intraplate volcanism in East Asia

NASA Astrophysics Data System (ADS)

He, L.

2016-12-01

Geophysical observations imply the intraplate volcanism in East Asia is related to dehydration of slab stagnating in the transition zone. To better understand the dynamics of such process, a thermochemical mantle convection model is constructed to simulate numerically the thermal evolution of slab and the transportation of water in the process of subduction. Equation of water transfer is explicitly included, and water effects on density and viscosity are considered. Modeling results indicate that behavior of water transport relates closely to the transient thermal state and viscosities both of the slab and the surrounding mantle. Generally, initiation of wet plume is mainly influenced by the viscosity of the wet layer in the uppermost slab, whereas the horizontal distance of water transport and its ascending rate is affected strongly by the viscosity of the big mantle wedge. Whether water can be carried successfully by slab into the mantle transition zone and trigger wet plume at the surface of flattening slab depends on the viscosity contrast between wet layer and surrounding mantle. The complex fluid flow superposed by corner flow and free thermal convection controls the water transport pattern in the upper mantle. Modeling results together with previous modeling infer three stages of water circulation in the big mantle wedge: 1) water is brought into the mantle transition zone by downward subducting slab when water layer viscosity is much higher than the wedge viscosity, otherwise water is released at shallow depth near wedge tip; 2) wet plume generates from surface of warm flattening slab if containing water, which arrives at the lithospheric base and induces melting; and 3) water spreads all over the big mantle wedge, mantle convection within the big mantle wedge becomes more active, leading to upwelling of asthenosphere and erosion of the overriding continental lithosphere. Wet plume from the flattening Pacific Plate can explain the intraplate Cenozoic volcanoes in East Asia.

Deep particle bed dryout model based on flooding

DOE Office of Scientific and Technical Information (OSTI.GOV)

Kuan, P.

1987-01-01

Examination of the damaged Three Mile island Unit 2 (TMI-2) reactor indicates that a deep (approx. 1-m) bed of relatively large (approx. 1-mm) particles was formed in the core. Cooling of such beds is crucial to the arrest of core damage progression. The Lipinski model, based on flows in the bed, has been used to predict the coolability, but uncertainties exist in the turbulent permeability. Models based on flooding at the top of the bed either have a dimensional viscosity term, or no viscosity dependence, thus limiting their applicability. This paper presents a dimensionless correlation based on flooding data thatmore » involves a liquid Reynolds number. The derived dryout model from this correlation is compared with data for deep beds of large particles at atmospheric pressure, and with other models over a wide pressure range. It is concluded that the present model can give quite accurate predictions for the dryout heat flux of particle beds formed during a light water reactor accident and it is easy to use and agrees with the Lipinski n = 5 model, which requires iterative calculations.« less

Constraints on mantle viscosity from convection models with plate motion history

NASA Astrophysics Data System (ADS)

Mao, W.; Zhong, S.

2017-12-01

The Earth's long-wavelength geoid and dynamic topography are mainly controlled by the mantle buoyancy and viscosity structure. Previous dynamical models for the geoid provide constraints on the 1-D mantle viscosity, using mantle buoyancy derived from seismic topography models. However, it is a challenge in these studies on how to convert seismic velocity to density anomalies and mantle buoyancy. Furthermore, these studies provide constraints only on relative viscosity variations but not on absolute magnitude of viscosity. In this study, we formulate time-dependent 3-D spherical mantle convection models with imposed plate motion history and seek constraints on mantle viscosity structure for both its radial relative variations and its absolute magnitude (i.e., Rayleigh number), using the geoid from the convection models. We found that the geoid at intermediate wavelengths of degrees 4-9 is mainly controlled by the subducted slabs in the upper mantle and the upper part of lower mantle that result from subduction from the last 50 Myr or the Cenozoic. To fit the degrees 4-9 geoid, we need viscosity contrast β defined as the ratio of the lower mantle viscosity and the asthenospheric viscosity to be larger than 2000 and Ra to be 1e8 (defined by the Earth's radius). The best fit model leads to 57% variance reduction and 76% correlation between the model and the observations. However, the long-wavelength geoid at degrees 2-3 is controlled by the lower mantle structure which requires much longer time scale to develop, as seen from our modeling. The preferred viscosity structure and Rayleigh number as constrained by the Cenozoic plate motion and the degrees 4-9 geoid no longer provide adequate fit to the geoid in models with the plate motion history for the last 450 Myr. The degrees 4-9 geoid amplitude is smaller for the models with longer plate motion history and a smaller Ra is required to fit the observation. In order to satisfy the relative amplitude between degrees 2-3 and degrees 4-9 geoid, either a gradually increase of viscosity in the upper part of lower mantle or larger thermal expansivity in the lower mantle is needed. We also consider thermo-chemical models to examine the effects of the African and Pacific thermochemical piles (i.e., LLSVPSs) on the geoid and the inferred mantle viscosity and Ra.

Two-Point Microrheology of Phase-Separated Domains in Lipid Bilayers

PubMed Central

Hormel, Tristan T.; Reyer, Matthew A.; Parthasarathy, Raghuveer

2015-01-01

Though the importance of membrane fluidity for cellular function has been well established for decades, methods for measuring lipid bilayer viscosity remain challenging to devise and implement. Recently, approaches based on characterizing the Brownian dynamics of individual tracers such as colloidal particles or lipid domains have provided insights into bilayer viscosity. For fluids in general, however, methods based on single-particle trajectories provide a limited view of hydrodynamic response. The technique of two-point microrheology, in which correlations between the Brownian dynamics of pairs of tracers report on the properties of the intervening medium, characterizes viscosity at length-scales that are larger than that of individual tracers and has less sensitivity to tracer-induced distortions, but has never been applied to lipid membranes. We present, to our knowledge, the first two-point microrheological study of lipid bilayers, examining the correlated motion of domains in phase-separated lipid vesicles and comparing one- and two-point results. We measure two-point correlation functions in excellent agreement with the forms predicted by two-dimensional hydrodynamic models, analysis of which reveals a viscosity intermediate between those of the two lipid phases, indicative of global fluid properties rather than the viscosity of the local neighborhood of the tracer. PMID:26287625

The viscosity of pāhoehoe lava: In situ syn-eruptive measurements from Kilauea, Hawaii

NASA Astrophysics Data System (ADS)

Chevrel, Magdalena Oryaëlle; Harris, Andrew J. L.; James, Mike R.; Calabrò, Laura; Gurioli, Lucia; Pinkerton, Harry

2018-07-01

Viscosity is one of the most important physical properties controlling lava flow dynamics. Usually, viscosity is measured in the laboratory where key parameters can be controlled but can never reproduce the natural environment and original state of the lava in terms of crystal and bubble contents, dissolved volatiles, and oxygen fugacity. The most promising approach for quantifying the rheology of molten lava in its natural state is therefore to carry out direct field measurements by inserting a viscometer into the lava while it is flowing. Such in-situ syn-eruptive viscosity measurements are notoriously difficult to perform due to the lack of appropriate instrumentation and the difficulty of working on or near an active lava flow. In the field, rotational viscometer measurements are of particular value as they have the potential to measure the properties of the flow interior rather than an integration of the viscosity of the viscoelastic crust + flow interior. To our knowledge only one field rotational viscometer is available, but logistical constraints have meant that it has not been used for 20 yr. Here, we describe new viscosity measurements made using the refurbished version of this custom-built rotational viscometer, as performed on active pāhoehoe lobes from the 61G lava flow of Kilauea's Pu'u 'Ō'ō eruption in 2016. We successfully measured a viscosity of ∼380 Pa s at strain-rates between 1.6 and 5 s-1 and at 1144 °C. Additionally, synchronous lava sampling allowed us to provide detailed textural and chemical characterization of quenched samples. Application of current physico-chemical models based on this characterization (16 ± 4 vol.% crystals; 50 ± 6 vol.% vesicles), gave viscosity estimates that were approximately compatible with the measured values, highlighting the sensitivity of model-based viscosity estimates on the effect of deformable bubbles. Our measurements also agree on the range of viscosities in comparison to previous field experiments on Hawaiian lavas. Conversely, direct comparison with sub-liquidus rheological laboratory measurements on natural lavas was unsuccessful because recreating field conditions (in particular volatile and bubble content) is so far inaccessible in the laboratory. Our work shows the value of field rotational viscometry fully-integrated with sample characterization to quantify three-phase lava viscosity. Finally, this work suggests the need for the development of a more versatile instrument capable of recording precise measurements at low torque and low strain rate, and with synchronous temperature measurements.

Tidal effects in differentiated viscoelastic bodies: a numerical approach

NASA Astrophysics Data System (ADS)

Walterová, M.; Běhounková, M.

2017-09-01

The majority of confirmed terrestrial exoplanets orbits close to their host stars and their evolution was likely altered by tidal interaction. Nevertheless, due to their viscoelastic properties on the tidal frequencies, their response cannot be described exactly by standardly employed constant-lag models. We therefore introduce a tidal model based on the numerical evaluation of a continuum mechanics problem describing the deformation of viscoelastic (Maxwell or Andrade) planetary mantles subjected to external force. We apply the method on a model Earth-size planet orbiting a low-mass star and study the effect of the orbital eccentricity, the mantle viscosity and the chosen rheology on the tidal dissipation, the complex Love numbers and the tidal torque. The number of stable spin states (i.e., zero tidal torque) grows with increasing mantle viscosity, similarly to the analytical model of Correia et al. (Astron Astrophys 571:A50, 2014) for homogeneous bodies. This behavior is only slightly influenced by the rheology used. Similarly, the Love numbers do not distinctly depend on the considered rheological model. The increase in viscosity affects the amplitude of their variations. The tidal heating described by the Maxwell rheology attains local minima associated with low spin-orbit resonances, with depth and shape depending on both the eccentricity and the viscosity. For the Andrade rheology, the minima at low resonances are very shallow and the tidal heating for all viscosities resembles a "fluid limit." The tidal heating is the quantity influenced the most by the rheology, having thus possible impact on the internal thermal evolution.

Temperature Dependence of Density, Viscosity and Electrical Conductivity for Hg-Based II-VI Semiconductor Melts

NASA Technical Reports Server (NTRS)

Li, C.; Ban, H.; Lin, B.; Scripa, R. N.; Su, C.-H.; Lehoczky, S. L.

2004-01-01

The relaxation phenomenon of semiconductor melts, or the change of melt structure with time, impacts the crystal growth process and the eventual quality of the crystal. The thermophysical properties of the melt are good indicators of such changes in melt structure. Also, thermophysical properties are essential to the accurate predication of the crystal growth process by computational modeling. Currently, the temperature dependent thermophysical property data for the Hg-based II-VI semiconductor melts are scarce. This paper reports the results on the temperature dependence of melt density, viscosity and electrical conductivity of Hg-based II-VI compounds. The melt density was measured using a pycnometric method, and the viscosity and electrical conductivity were measured by a transient torque method. Results were compared with available published data and showed good agreement. The implication of the structural changes at different temperature ranges was also studied and discussed.

Simultaneous feedback control of plasma rotation and stored energy on NSTX-U using neoclassical toroidal viscosity and neutral beam injection

NASA Astrophysics Data System (ADS)

Goumiri, I. R.; Rowley, C. W.; Sabbagh, S. A.; Gates, D. A.; Boyer, M. D.; Gerhardt, S. P.; Kolemen, E.; Menard, J. E.

2017-05-01

A model-based feedback system is presented enabling the simultaneous control of the stored energy through βn and the toroidal rotation profile of the plasma in National Spherical Torus eXperiment Upgrade device. Actuation is obtained using the momentum from six injected neutral beams and the neoclassical toroidal viscosity generated by applying three-dimensional magnetic fields. Based on a model of the momentum diffusion and torque balance, a feedback controller is designed and tested in closed-loop simulations using TRANSP, a time dependent transport analysis code, in predictive mode. Promising results for the ongoing experimental implementation of controllers are obtained.

Simultaneous feedback control of plasma rotation and stored energy on NSTX-U using neoclassical toroidal viscosity and neutral beam injection

PubMed Central

Goumiri, I. R.; Sabbagh, S. A.; Boyer, M. D.; Gerhardt, S. P.; Kolemen, E.; Menard, J. E.

2017-01-01

A model-based feedback system is presented enabling the simultaneous control of the stored energy through βn and the toroidal rotation profile of the plasma in National Spherical Torus eXperiment Upgrade device. Actuation is obtained using the momentum from six injected neutral beams and the neoclassical toroidal viscosity generated by applying three-dimensional magnetic fields. Based on a model of the momentum diffusion and torque balance, a feedback controller is designed and tested in closed-loop simulations using TRANSP, a time dependent transport analysis code, in predictive mode. Promising results for the ongoing experimental implementation of controllers are obtained. PMID:28435207

Simultaneous feedback control of plasma rotation and stored energy on NSTX-U using neoclassical toroidal viscosity and neutral beam injection

DOE PAGES

Goumiri, I. R.; Rowley, C. W.; Sabbagh, S. A.; ...

2017-02-23

In this study, a model-based feedback system is presented enabling the simultaneous control of the stored energy through β n and the toroidal rotation profile of the plasma in National Spherical Torus eXperiment Upgrade device. Actuation is obtained using the momentum from six injected neutral beams and the neoclassical toroidal viscosity generated by applying three-dimensional magnetic fields. Based on a model of the momentum diffusion and torque balance, a feedback controller is designed and tested in closed-loop simulations using TRANSP, a time dependent transport analysis code, in predictive mode. Promising results for the ongoing experimental implementation of controllers are obtained.

The Impacts of 3-D Earth Structure on GIA-Induced Crustal Deformation and Future Sea-Level Change in the Antarctic

NASA Astrophysics Data System (ADS)

Powell, E. M.; Hay, C.; Latychev, K.; Gomez, N. A.; Mitrovica, J. X.

2016-12-01

Glacial Isostatic Adjustment (GIA) models used to constrain the extent of past ice sheets and viscoelastic Earth structure, or to correct geodetic and geological observables for ice age effects, generally only consider depth-dependent variations in Earth viscosity and lithospheric structure. A et al. [2013] argued that 3-D Earth structure could impact GIA observables in Antarctica, but concluded that the presence of such structure contributes less to GIA uncertainty than do differences in Antarctic deglaciation histories. New seismic and geological evidence, however, indicates the Antarctic is underlain by complex, high amplitude variability in viscoelastic structure, including a low viscosity zone (LVZ) under West Antarctica. Hay et al. [2016] showed that sea-level fingerprints of modern melting calculated using such Earth models differ from those based on elastic or 1-D viscoelastic Earth models within decades of melting. Our investigation is motivated by two questions: (1) How does 3-D Earth structure, especially this LVZ, impact observations of GIA-induced crustal deformation associated with the last deglaciation? (2) How will 3-D Earth structure affect predictions of future sea-level rise in Antarctica? We compute the gravitationally self-consistent sea level, uplift, and gravity changes using the finite volume treatment of Latychev et al. [2005]. We consider four viscoelastic Earth models: a global 1-D model; a regional, West Antarctic-like 1-D model; a 3-D model where the lithospheric thickness varies laterally; and a 3-D model where both viscosity and lithospheric thickness vary laterally. For our Last Glacial Maximum to present investigations we employ ICE6g [Peltier et al., 2015]. For our present-future investigations we consider a melt scenario consistent with GRACE satellite gravity derived solutions [Harig et al., 2015]. Our calculations indicate that predictions of crustal deformations due to both GIA and ongoing melting are strongly influenced by 3-D lithospheric thickness and viscosity structure. Future sea level change due to ongoing melting is primarily influenced by 3-D viscosity structure. We show that 1-D Earth models built using regional inferences of viscosity and lithospheric thickness do not accurately capture the variability introduced by 3-D Earth structure.

The Impacts of 3-D Earth Structure on GIA-Induced Crustal Deformation and Future Sea-Level Change in the Antarctic

NASA Astrophysics Data System (ADS)

Powell, E. M.; Hay, C.; Latychev, K.; Gomez, N. A.; Mitrovica, J. X.

2017-12-01

Glacial Isostatic Adjustment (GIA) models used to constrain the extent of past ice sheets and viscoelastic Earth structure, or to correct geodetic and geological observables for ice age effects, generally only consider depth-dependent variations in Earth viscosity and lithospheric structure. A et al. [2013] argued that 3-D Earth structure could impact GIA observables in Antarctica, but concluded that the presence of such structure contributes less to GIA uncertainty than do differences in Antarctic deglaciation histories. New seismic and geological evidence, however, indicates the Antarctic is underlain by complex, high amplitude variability in viscoelastic structure, including a low viscosity zone (LVZ) under West Antarctica. Hay et al. [2016] showed that sea-level fingerprints of modern melting calculated using such Earth models differ from those based on elastic or 1-D viscoelastic Earth models within decades of melting. Our investigation is motivated by two questions: (1) How does 3-D Earth structure, especially this LVZ, impact observations of GIA-induced crustal deformation associated with the last deglaciation? (2) How will 3-D Earth structure affect predictions of future sea-level rise in Antarctica? We compute the gravitationally self-consistent sea level, uplift, and gravity changes using the finite volume treatment of Latychev et al. [2005]. We consider four viscoelastic Earth models: a global 1-D model; a regional, West Antarctic-like 1-D model; a 3-D model where the lithospheric thickness varies laterally; and a 3-D model where both viscosity and lithospheric thickness vary laterally. For our Last Glacial Maximum to present investigations we employ ICE6g [Peltier et al., 2015]. For our present-future investigations we consider a melt scenario consistent with GRACE satellite gravity derived solutions [Harig et al., 2015]. Our calculations indicate that predictions of crustal deformations due to both GIA and ongoing melting are strongly influenced by 3-D lithospheric thickness and viscosity structure. Future sea level change due to ongoing melting is primarily influenced by 3-D viscosity structure. We show that 1-D Earth models built using regional inferences of viscosity and lithospheric thickness do not accurately capture the variability introduced by 3-D Earth structure.

Modeling of viscoelastic properties of nonpermeable porous rocks saturated with highly viscous fluid at seismic frequencies at the core scale

NASA Astrophysics Data System (ADS)

Wang, Zizhen; Schmitt, Douglas R.; Wang, Ruihe

2017-08-01

A core scale modeling method for viscoelastic properties of rocks saturated with viscous fluid at low frequencies is developed based on the stress-strain method. The elastic moduli dispersion of viscous fluid is described by the Maxwell's spring-dash pot model. Based on this modeling method, we numerically test the effects of frequency, fluid viscosity, porosity, pore size, and pore aspect ratio on the storage moduli and the stress-strain phase lag of saturated rocks. And we also compared the modeling results to the Hashin-Shtrikman bounds and the coherent potential approximation (CPA). The dynamic moduli calculated from the modeling are lower than the predictions of CPA, and both of these fall between the Hashin-Shtrikman bounds. The modeling results indicate that the frequency and the fluid viscosity have similar effects on the dynamic moduli dispersion of fully saturated rocks. We observed the Debye peak in the phase lag variation with the change of frequency and viscosity. The pore structure parameters, such as porosity, pore size, and aspect ratio affect the rock frame stiffness and result in different viscoelastic behaviors of the saturated rocks. The stress-strain phase lags are larger with smaller stiffness contrasts between the rock frame and the pore fluid. The viscoelastic properties of saturated rocks are more sensitive to aspect ratio compared to other pore structure parameters. The results suggest that significant seismic dispersion (at about 50-200 Hz) might be expected for both compressional and shear waves passing through rocks saturated with highly viscous fluids.