Ocean Ultrasonic Shear and Compression Viscosities

[EN]A comprehensive description of ocean molecular flow and deformation is provided with the help of hydrodynamic and ultrasonic principles. Hydrodynamic computation of true or natural viscosities shows that ocean shear viscosity (?G), compression viscosity (?K), and extensional viscosity (?E) are interrelated. There are no experimental methods available for the in situ measurement of these viscosities. Sound absorption coefficients (? obs) allow to know the ultrasonic shear (?UG), compression (?UK), and longitudinal (?L) viscosities, which decrease with increasing frequency and increase with increasing temperature, the flow activation energies having nearly equivalent values; pressure (depth) increase/decrease them at low/high frequencies. The viscosities ?* UG, ?* UK, ?* L are approached at about 1000 KHz. They decrease with temperature and pressure, and increase with salinity. The ?*UG becomes equal to the true shear viscosity ? G at the viscosity ratio ? = ?UK / ?UG = 0.

Pendant drop formation of shear-thinning and yield stress fluids

A volume-of-fluid numerical method is used to predict the dynamics of shear-thinning liquid drop formation in air from a circular orifice. The validity of the numerical calculation is confirmed for a Newtonian liquid by comparison with experimental measurements. For particular values of Weber number and Froude number, predictions show a more rapid pinch-off, and a reduced number of secondary droplets, with increasing shear-thinning. Also a minimum in the limiting drop length occurs for the smallest Weber number as the zero-shear viscosity is varied. At the highest viscosity, the drop length is reduced due to shear-thinning, whereas at lower viscosities there is little effect of shear-thinning. The evolution of predicted drop shape, drop thickness and length, and the configuration at pinch-off are discussed for shear-thinning drops. The evolution of a drop of Bingham yield stress liquid is also considered as a limiting case. In contrast to the shear-thinning cases, it exhibits a plug flow prior to necking, an almost step-change approach to pinch-off of a torpedo shaped drop following the onset of necking, and a much smaller neck length; no secondary drops are formed. The results demonstrate the potential of the numerical model as a design tool in tailoring the fluid rheology for controlling drop formation behaviour. (c) 2006 Elsevier Inc. All rights reserved.

Drop formation and breakup of low viscosity elastic fluids: Effects of molecular weight and concentration

The dynamics of drop formation and pinch-off have been investigated for a series of low viscosity elastic fluids possessing similar shear viscosities, but differing substantially in elastic properties. On initial approach to the pinch region, the viscoelastic fluids all exhibit the same global necking behavior that is observed for a Newtonian fluid of equivalent shear viscosity. For these low viscosity dilute polymer solutions, inertial and capillary forces form the dominant balance in this potential flow regime, with the viscous force being negligible. The approach to the pinch point, which corresponds to the point of rupture for a Newtonian fluid, is extremely rapid in such solutions, with the sudden increase in curvature producing very large extension rates at this location. In this region the polymer molecules are significantly extended, causing a localized increase in the elastic stresses, which grow to balance the capillary pressure. This prevents the necked fluid from breaking off, as would occur in the equivalent Newtonian fluid. Alternatively, a cylindrical filament forms in which elastic stresses and capillary pressure balance, and the radius decreases exponentially with time. A (0+1)-dimensional finitely extensible nonlinear elastic dumbbell theory incorporating inertial, capillary, and elastic stresses is able to capture the basic features of the experimental observations. Before the critical "pinch time" t(p), an inertial-capillary balance leads to the expected 2/3-power scaling of the minimum radius with time: R-min similar to(t(p)-t)(2/3). However, the diverging deformation rate results in large molecular deformations and rapid crossover to an elastocapillary balance for times t>t(p). In this region, the filament radius decreases exponentially with time R-min similar to exp[(t(p)-t)/lambda(1)], where lambda(1) is the characteristic time constant of the polymer molecules. Measurements of the relaxation times of polyethylene oxide solutions of varying concentrations and molecular weights obtained from high speed imaging of the rate of change of filament radius are significantly higher than the relaxation times estimated from Rouse-Zimm theory, even though the solutions are within the dilute concentration region as determined using intrinsic viscosity measurements. The effective relaxation times exhibit the expected scaling with molecular weight but with an additional dependence on the concentration of the polymer in solution. This is consistent with the expectation that the polymer molecules are in fact highly extended during the approach to the pinch region (i.e., prior to the elastocapillary filament thinning regime) and subsequently as the filament is formed they are further extended by filament stretching at a constant rate until full extension of the polymer coil is achieved. In this highly extended state, intermolecular interactions become significant, producing relaxation times far above theoretical predictions for dilute polymer solutions under equilibrium conditions. (C) 2006 American Institute of Physics

Rheology of particle-loaded semi-dilute polymer solutions

We present measurements of the rheology of suspensions of rigid spheres in a semi-dilute polymer solution from experiments of steady and oscillatory shear. For a given value of the shear rate gamma, addition of particles enhances the viscosity and the first normal stress difference but decreases the magnitude of the second normal stress difference. The viscosity eta exhibits a power law variation in gamma for a range of gamma that grows with phi. The first normal stress N-1 is positive and its value grows with phi; it exhibits a clear power law variation for the entire range of gamma that was studied. The second normal stress difference N-2 is negative for the pure polymer solution and much smaller in magnitude than N-1; on addition of particles, its magnitude further decreases, and it appears to change sign at large phi. The behavior of N-1 and N-2 is at odds with the findings of recent studies on particle-loaded dilute polymer solutions and polymer melts. The small-amplitude oscillatory shear experiments show the linear viscoelastic properties, G(') and G('), increasing with phi at a given value of the angular frequency omega. The dynamic viscosity of the suspension differs substantially from its steady shear viscosity, and the difference increases as gamma, omega -> 0.

Re-entrant Phase Behavior of a Concentrated Anionic Surfactant System with Strongly Binding Counterions

The phase behavior of the anionic surfactant sodium dodecyl sulfate (SDS) in the presence of the strongly binding counterion p-toluidine hydrochloride (PTHC) has been examined using small-angle X-ray diffraction and polarizing microscopy. A hexagonal-to-lamellar transition on varying the PTHC to SDS molar ratio (alpha) occurs through a nematic phase of rodlike micelles (N-C) -> isotropic (I) -> nematic of disklike micelles (N-D) at a fixed surfactant concentration (phi). The lamellar phase is found to coexist with an isotropic phase (l') over a large region of the phase diagram. Deuterium nuclear magnetic resonance investigations of the phase behavior at phi = 0.4 confirm the transition from N-C to N-D on varying alpha. The viscoelastic and flow behaviors of the different phases were examined. A decrease in the steady shear viscosity across the different phases with increasing alpha suggests a decrease in the aspect ratio of the micellar aggregates. From the transient shear stress response of the N-C and N-D nematic phases in step shear experiments, they were characterized to be tumbling and now aligning, respectively. Our studies reveal that by tuning the morphology of the surfactant micelles strongly binding counterions modify the phase behavior and rheological properties of concentrated surfactant solutions.

The harmonic torsional oscillations of a thin disk submerged in a fluid with a surfactant surface layer

A formulation in terms of a Fredholm integral equation of the first kind is given for the axisymmetric problem of a disk oscillating harmonically in a viscous fluid whose surface is contaminated with a surfactant film. The equation of the first kind is converted to a pair of coupled integral equations of the second kind, which are solved numerically. The resistive torque on the disk is evaluated and surface velocity profiles are computed for varying values of the ratio of the coefficient of surface shear viscosity to the coefficient of viscosity of the substrate fluid, and the depth of the disk below the surface.

Oral medication delivery in impaired swallowing: Thickening liquid medications for safe swallowing alters dissolution characteristics

Acetaminophen (paracetamol) is available in a wide range of oral formulations designed to meet the needs of the population across the age-spectrum, but for people with impaired swallowing, i.e. dysphagia, both solid and liquid medications can be difficult to swallow without modification. The effect of a commercial polysaccharide thickener, designed to be added to fluids to promote safe swallowing by dysphagic patients, on rheology and acetaminophen dissolution was tested using crushed immediate-release tablets in water, effervescent tablets in water, elixir and suspension. The inclusion of the thickener, comprised of xanthan gum and maltodextrin, had a considerable impact on dissolution; acetaminophen release from modified medications reached 12-50% in 30 minutes, which did not reflect the pharmacopeia specification for immediate release preparations. Flow curves reflect the high zero-shear viscosity and the apparent yield stress of the thickened products. The weak gel nature, in combination with high G’ values compared to G” (viscoelasticity) and high apparent yield stress, impact drug release. The restriction on drug release from these formulations is not influenced by the theoretical state of the drug (dissolved or dispersed), and the approach typically used in clinical practice (mixing crushed tablets into pre-prepared thickened fluid) cannot be improved by altering the order of incorporation or mixing method.

Correlations and the Ring-Kinetic Equation in Dense Sheared Granular Flows

A formal way of deriving fluctuation-correlation relations in dense sheared granular media, starting with the Enskog approximation for the collision integral in the Chapman-Enskog theory, is discussed. The correlation correction to the viscosity is obtained using the ring-kinetic equation, in terms of the correlations in the hydrodynamic modes of the linearised Enskog equation. It is shown that the Green-Kubo formula for the shear viscosity emerges from the two-body correlation function obtained from the ring-kinetic equation.

Effect of screening of intermicellar interactions on the linear and nonlinear rheology of a viscoelastic gel

We report our studies of the linear and nonlinear rheology of aqueous solutions of the surfactant cetyl trimethylammonium tosylate (CTAT) with varying amounts of sodium chloride (NaCl). The CTAT concentration is fixed at 42 mM, and the salt concentration is varied between 0 and 120 mM. On increasing the salt (NaCl) concentration, we see three distinct regimes in the zero-shear viscosity and the high-frequency plateau modulus data. In regime 1, the zero-shear viscosity shows a weak increase with salt concentration due to enhanced micellar growth. The decrease in the zero-shear viscosities with salt concentration in regimes II and III can be explained in terms of intermicellar branching. The most intriguing feature of our data, however, is the anomalous behavior of the high-frequency plateau modulus in regime II (0.12 less than or equal to [NaCl]/[CTAT] less than or equal to 1.42). In this regime, the plateau modulus increases with an increase in NaCl concentration. This is highly interesting, since the correlation length of concentration fluctuations and hence the plateau modulus G(0) are not expected to change appreciably in the semidilute regime. We propose to explain the changes in regime II in terms of a possible unbinding of the organic counterions (tosylate) from the CTA(+) surfaces on the addition of NaCl. In the nonlinear flow curves of the samples with high salt content, significant deviations from the predictions of the Giesekus model for entangled micelles are observed.

Second order transport from anomalies

We study parity odd transport at second order in derivative expansion for a non-conformal charged fluid. We see that there are 27 parity odd transport coefficients, of which 12 are non-vanishing in equilibrium. We use the equilibrium partition function method to express 7 of these in terms of the anomaly, shear viscosity, charge diffusivity and thermodynamic functions. The remaining 5 are constrained by 3 relations which also involve the anomaly. We derive Kubo formulae for 2 of the transport coefficients and show these agree with that derived from the equilibrium partition function.

Holographic charge transport in non commutative gauge theories

In this paper, based on the holographic techniques, we explore the hydrodynamics of charge diffusion phenomena in non commutative N = 4 SYM plasma at strong coupling. In our analysis, we compute the R charge diffusion rates both along commutative as well as the non commutative coordinates of the brane. It turns out that unlike the case for the shear viscosity, the DC conductivity along the non commutative direction of the brane differs significantly from that of its cousin corresponding to the commutative direction of the brane. Such a discrepancy however smoothly goes away in the limit of the vanishing non commutativity.

Thermalization of Green functions and quasinormal modes

We develop a new method to study the thermalization of time dependent retarded Green function in conformal field theories holographically dual to thin shell AdS Vaidya space times. The method relies on using the information of all time derivatives of the Green function at the shell and then evolving it for later times. The time derivatives of the Green function at the shell is given in terms of a recursion formula. Using this method we obtain analytic results for short time thermalization of the Green function. We show that the late time behaviour of the Green function is determined by the first quasinormal mode. We then implement the method numerically. As applications of this method we study the thermalization of the retarded time dependent Green function corresponding to a minimally coupled scalar in the AdS 3 and AdS 5 thin Vaidya shells. We see that as expected the late time behaviour is determined by the first quasinormal mode. We apply the method to study the late time behaviour of the shear vector mode in AdS 5 Vaidya shell. At small momentum the corresponding time dependent Green function is expected to relax to equilibrium by the shear hydrodynamic mode. Using this we obtain the universal ratio of the shear viscosity to entropy density from a time dependent process.

Microfluidic Rheometry

The development and growth of microfluidics has stimulated interest in the behaviour of complex liquids in micro-scale geometries and provided a rich platform for rheometric investigations of non-Newtonian phenomena at small scales. Microfluidic techniques present the rheologist with new opportunities for material property measurement and this review discusses the use of microfluidic devices to measure bulk rheology in both shear and extensional flows. Capillary, stagnation and contraction flows are presented in this context and developments, limitations and future perspectives are examined. (C) 2008 Elsevier Ltd. All rights reserved.

The Equation of State and Petrogenesis of Komatiite

(1) Equation of State of Komatiite

The equation of state (EOS) of a molten komatiite (27 wt% MgO) was detennined in the 5 to 36 GPa pressure range via shock wave compression from 1550°C and 0 bar. Shock wave velocity, U_S, and particle velocity, U_P, in km/s follow the linear relationship U_S = 3.13(±0.03) + 1.47(±0.03) U_P. Based on a calculated density at 1550°C, 0 bar of 2.745±0.005 glee, this U_S-U_P relationship gives the isentropic bulk modulus K_S = 27.0 ± 0.6 GPa, and its first and second isentropic pressure derivatives, K'_S = 4.9 ± 0.1 and K"_S = -0.109 ± 0.003 GPa^-1.

The calculated liquidus compression curve agrees within error with the static compression results of Agee and Walker [1988a] to 6 GPa. We detennine that olivine (FO₉₄) will be neutrally buoyant in komatiitic melt of the composition we studied near 8.2 GPa. Clinopyroxene would also be neutrally buoyant near this pressure. Liquidus garnet-majorite may be less dense than this komatiitic liquid in the 20-24 GPa interval, however pyropic-garnet and perovskite phases are denser than this komatiitic liquid in their respective liquidus pressure intervals to 36 GPa. Liquidus perovskite may be neutrally buoyant near 70 GPa.

At 40 GPa, the density of shock-compressed molten komatiite would be approximately equal to the calculated density of an equivalent mixture of dense solid oxide components. This observation supports the model of Rigden et al. [1989] for compressibilities of liquid oxide components. Using their theoretical EOS for liquid forsterite and fayalite, we calculate the densities of a spectrum of melts from basaltic through peridotitic that are related to the experimentally studied komatiitic liquid by addition or subtraction of olivine. At low pressure, olivine fractionation lowers the density of basic magmas, but above 14 GPa this trend is reversed. All of these basic to ultrabasic liquids are predicted to have similar densities at 14 GPa, and this density is approximately equal to the bulk (PREM) mantle. This suggests that melts derived from a peridotitic mantle may be inhibited from ascending from depths greater than 400 km.

The EOS of ultrabasic magmas was used to model adiabatic melting in a peridotitic mantle. If komatiites are formed by >15% partial melting of a peridotitic mantle, then komatiites generated by adiabatic melting come from source regions in the lower transition zone (≈500-670 km) or the lower mantle (>670 km). The great depth of incipient melting implied by this model, and the melt density constraint mentioned above, suggest that komatiitic volcanism may be gravitationally hindered. Although komatiitic magmas are thought to separate from their coexisting crystals at a temperature =200°C greater than that for modern MORBs, their ultimate sources are predicted to be diapirs that, if adiabatically decompressed from initially solid mantle, were more than 700°C hotter than the sources of MORBs and derived from great depth.

We considered the evolution of an initially molten mantle, i.e., a magma ocean. Our model considers the thermal structure of the magma ocean, density constraints on crystal segregation, and approximate phase relationships for a nominally chondritic mantle. Crystallization will begin at the core-mantle boundary. Perovskite buoyancy at > 70 GPa may lead to a compositionally stratified lower mantle with iron-enriched mangesiowiistite content increasing with depth. The upper mantle may be depleted in perovskite components. Olivine neutral buoyancy may lead to the formation of a dunite septum in the upper mantle, partitioning the ocean into upper and lower reservoirs, but this septum must be permeable.

(2) Viscosity Measurement with Shock Waves

We have examined in detail the analytical method for measuring shear viscosity from the decay of perturbations on a corrugated shock front The relevance of initial conditions, finite shock amplitude, bulk viscosity, and the sensitivity of the measurements to the shock boundary conditions are discussed. The validity of the viscous perturbation approach is examined by numerically solving the second-order Navier-Stokes equations. These numerical experiments indicate that shock instabilities may occur even when the Kontorovich-D'yakov stability criteria are satisfied. The experimental results for water at 15 GPa are discussed, and it is suggested that the large effective viscosity determined by this method may reflect the existence of ice VII on the Rayleigh path of the Hugoniot This interpretation reconciles the experimental results with estimates and measurements obtained by other means, and is consistent with the relationship of the Hugoniot with the phase diagram for water. Sound waves are generated at 4.8 MHz at in the water experiments at 15 GPa. The existence of anelastic absorption modes near this frequency would also lead to large effective viscosity estimates.

(3) Equation of State of Molybdenum at 1400°C

Shock compression data to 96 GPa for pure molybdenum, initially heated to 1400°C, are presented. Finite strain analysis of the data gives a bulk modulus at 1400°C, K'_S. of 244±2 GPa and its pressure derivative, K'_OS of 4. A fit of shock velocity to particle velocity gives the coefficients of U_S = C_O+S U_P to be C_O = 4.77±0.06 km/s and S = 1.43±0.05. From the zero pressure sound speed, C_O, a bulk modulus of 232±6 GPa is calculated that is consistent with extrapolation of ultrasonic elasticity measurements. The temperature derivative of the bulk modulus at zero pressure, θK_OSθT|_P, is approximately -0.012 GPa/K. A thermodynamic model is used to show that the thermodynamic Grüneisen parameter is proportional to the density and independent of temperature. The Mie-Grüneisen equation of state adequately describes the high temperature behavior of molybdenum under the present range of shock loading conditions.