A finite volume technique to simulate the flow of a viscoelastic fluid

Hydrodynamically developing flow of Oldroyd B fluid in the planar die entrance region has been investigated numerically using SIMPLER algorithm in a non-uniform staggered grid system. It has been shown that for constant values of the Reynolds number, the entrance length increases as the Weissenberg number increases. For small Reynolds number flows the center line velocity distribution exhibit overshoot near the inlet, which seems to be related to the occurrence of numerical breakdown at small values of the limiting Weissenberg number than those for large Reynolds number flows. The distributions of the first normal stress difference display clearly the development of the flow characteristics from extensional flow to shear flow.List of symbols D rate of strain tensor - L slit halfheight - P pressure, indeterminate part of the Cauchy stress tensor - R the Reynolds number - t time - U average velocity in the slit - u velocity vector - u,v velocity components - W the Weissenberg number based on the difference between stress relaxation time and retardation time - W ₁ the Weissenberg number based on stress relaxation time - x,y rectangular Cartesian coordinates - ratio of retardation time to stress relaxation time - zero-shear-rate viscosity, ₁ + ₂ - ₁ non-Newtonian contribution to - ₂ Newtonian contribution to - ₁ stress relaxation time - ₂ retardation time - density - (, , ) xx, yy and xy components of ₁, respectively - determinate part of the Cauchy stress tensor - ₁ non-Newtonian contribution to - ₂ Newtonian contribution to      

Shear viscosity of liquid4He and3He-4He mixtures,especially near the superfluid transition

High-resolution measurements of are reported for liquid⁴He and³He-⁴He mixtures at saturated vapor pressures between 1.2 and 4.2 K with particular emphasis on the superfluid transition. Here is the mass density, the shear viscosity, and in the superfluid phase both and are the contributions from the normal component of the fluid ( _n and _n ). The experiments were performed with a torsional oscillator operating at 151 Hz. The mole fraction X of³He in the mixtures ranged from 0.03 to 0.65. New data for the total density and data for _n by various authors led to the calculation of . For⁴He, the results for are compared with published ones, both in the normal and superfluid phases, and also with predictions in the normal phase both over a broad range and close to T. The behavior of and of in mixtures if presented. The sloped/dT near T and its change at the superfluid transition are found to decrease with increasing³He concentration. Measurements at one temperature of versus pressure indicate a decreasing dependence of on molar volume asX(³He) increases. Comparison of at T, the minimum of _n in the superfluid phase and the temperature of this minimum is made with previous measurements. Thermal conductivity measurements in the mixtures, carried out simultaneously with those of , revealed no difference in the recorded superfluid transition, contrary to earlier work. In the appendices, we present data from new measurements of the total density for the same mixtures used in viscosity experiments. Furthermore, we discuss the data for _n determined for⁴He and for³He-⁴He mixtures, and which are used in the analysis of the data.     

The Viscosity Surfaces of Propane in the Form of Multilayer Feed Forward Neural Networks

The present work focuses on the development of a viscosity equation =(,T) for propane through a multilayer feedforward neural network (MLFN) technique. Having been successfully applied to a variety of fluids so far, the proposed technique can be regarded as a general approach to viscosity modeling. The MLFN viscosity equation has been based on the available experimental data for propane: validation on the 969 primary data shows an average absolute deviation (AAD) of 0.29% in the temperature, pressure, and density range of applicability, i.e., 90 to 630 K, 0 to 60 MPa, and 0 to 730 kgm^–3. This result is very promising, especially when compared with experimental data uncertainty. The minimum amount of required data for setting up the MLFN has been investigated, to explore the minimum cost of the model. Comparisons with other viscosity models are presented regarding amount of input data, claimed accuracy, and range of applicability, with the aim of providing a guideline when viscosity has to be calculated for engineering purposes. A high accuracy equation of state for the conversion of variables from experimental P,T to operative ,T has to be provided. To overcome this requirement, two viscosity explicit equations in the form =(P,T) are also developed, for the liquid and for the vapor phases. The respective AADs are 0.58 and 0.22%, comparable with those of the former =(,T) equation. Finally, the trend of the experimental viscosity second virial coefficient is reproduced and compared with that obtained from the MLFN.     

Shear viscosity measurements of liquid 4He and 3He-4He mixtures near the tricritical point

A torsional oscillator cell is described, by means of which simultaneous precision measurements of () and of the molar volume can be made in liquid ⁴He-⁴He mixtures over the temperature range between 0.5 and 3 K. Here is the mass density, the shear viscosity and in the superfluid phase they become the contributions _n and _n of the normal component. The results of for ⁴He near the superfluid transition are compared with the predictions by Schloms, Pankert and Dohm, and by Ferrell. Measurements of () are reported for mixtures with 0.64X0.74, where X is the ³He mole fraction. Those for X = 0.67 and 0.70 are compared with data by Lai and Kitchens. The viscosity experiments show no evidence of a weak singularity at the tricritical point.     

Viscosity of steam in the critical region

The shear viscosity of fluids exhibits an anomalous enhancement in the close vicinity of the critical point. A detailed experimental study of the viscosity of steam in the critical region has been reported by Rivkin and collaborators. A reanalysis of the experimental data indicates that the behavior of the viscosity of steam near the critical point is similar to that observed for other fluids near the critical point. An interpolating equation for the viscosity of water and steam is presented that incorporates the critical viscosity enhancement.Nomenclature a critical region equation of state parameter - a _k coefficients in equation for ₀ - a _ij coefficients in equation for ¯ - b critical region equation of state parameter - c _p specific heat at constant pressure - c _v specific heat at constant volume - k critical region equation of state parameter - k _B Boltzmann constant - P pressure - P _r 22.115 MPa - P ^* P/P _r - P _c critical pressure - P _i coefficients in critical region equation of state - R~P (P-P _c )/P _c - q parameter in equation for critical viscosity enhancement - r parametric variable in critical region equation of state - T temperature in K (IPTS-48) - T _r 647.27 K - T ^* T/T _r - T _c critical temperature - T (T–T _c )/T _c - V volume - critical exponent of specific heat - critical exponent of coexistence curve - critical exponent of compressibility - critical exponent of chemical potential at T=T _c - dynamic viscosity - ₀ lim ₀ - ¯ normal viscosity - critical viscosity enhancement - ¯ thermal conductivity - normal thermal conductivity - critical thermal conductivity enhancement - parametric variable in critical region equation of state - correlation length - ₀ correlation length amplitude above T _c at = _c - critical exponent of correlation length - density - _r 317.763 kg/m³ - ^* / _r - _c critical density - (– _c )/ _c - _p estimated error of pressure - _T estimated error of temperature - estimated error of viscosity - exponent of critical viscosity enhancement - _t (/P) _T symmetrized compressibility - _T ^* _T P _r / _r ² - _{t t} P _c / _c ²     

p-Version hierarchical three dimensional curved shell element for heat conduction

This paper presents a finite element formulation for a three dimensional nine node p-version hierarchical curved shell element for heat conduction where the element temperature approximation can be of arbitrary order p , p , and p in the , and directions. This is accomplished by first, constructing one dimensional hierarchical approximation functions and the corresponding nodal variable operators for each of the three directions , and using Lagrange interpolating polynomials and then taking their products (sometimes also called tensor products). The element approximation functions as well as the nodal variables are hierarchical and therefore the element matrices and the equivalent heat vectors are hierarchical also i.e. the element properties corresponding to polynomial orders p , p , and p are a subset of those corresponding to (p +1), (p +1), and (p +1). The element formulation ensures C ⁰ continuity. The curved shell geometry is constructed in the usual way by taking the coordinates of the nodes lying on the middle surface of the element (=0) and the nodal thickness vectors. The element properties i.e. element matrices and the equivalent heat vectors are derived using weak formulation (or quadratic functional) of the three dimensional F ourier heat conduction equation and the hierarchical element temperature approximation. The element formulation is equally effective for very thin as well as extremely thick shells. Numerical examples are presented to demonstrate the accuracy, efficiency, modeling convenience, faster rate of convergence and over all superiority of the present formulation. The h-approximation results are presented for comparison purposes.     

The viscosity and some related properties of liquid3He at the melting curve between 1 and 100 mK

The viscosity of liquid ³ He has been measured along the melting curve from 1 to 100 mK by means of a vibrating wire viscometer. In the normal Fermiliquid region we find 1/T² = 1.17–3.10T, where is in P and T in K. At the transition temperature T _A = 2.6 mK a rapid decrease occurs in _n , the viscosity of the normal component. Within 0.3 mK below T _A , _n decreases to about 25% of _A, but then becomes essentially constant. In the B phase _n first decreases to 20% of _A and then seems to increase below 1.4 mK. Data on _n , the density of the normal component, are also presented in the A and B phases. The results show that viscous flow is accompanied by a flow of zero dissipation, thus proving superfluidity in the A and B phases. The viscosity data at magnetic fields up to 0.9T have been related to theoretical calculations of the energy gap of superfluid ³ He near T _A . The splitting of the A transition and the suppression of the B phase in an external field were also measured.     

Non-Newtonian properties of emulsions in solutions of surface-active agents

The adsorption of surface-active agents (surfactants) on channels changes the effective viscosity of an emulsion and gives it non-Newtonian properties.Notation a drop radius - C and C° concentration of surfactant near a drop and averaged over the medium - D diffusion coefficient of the surfactant - I second-rank unit tensor - n_i components of the unit normal vector to the drop surface - r radius vector directed from the center of the drop - s surface area of the drop occupied by one molecule of surfactant - t_j and T characteristic times - and sorption and desorption constants - gG and ° true and equilibrium surface concentrations of surfactant - , , and gh₁ effective viscosity and the viscosities of the disperse phase and dispersion medium - volume concentration of the disperse phase - surface tension of the drop - T_j, and characteristic times Indices + and * quantities near and inside the drop - t tangential components of vectors and tensors. The operators div_s and grad_s have the same meaning as the ordinary divergence and gradient operators, but with fixed values of the radius vector ra Translated from Inzhenerno-Fizicheskii Zhurnal, Vol. 56, No. 5, pp. 787–793, May, 1989.     

Complex microstructural changes in as-cast eutectoid Zn-Al alloy

Phase transformations and microstructural changes of an as-cast eutectoid Zn-Al alloy (ZnAl₂₂Cu₂) were investigated during isothermal holding. The typical dendritic structures consisted of _s phase as a core with the edge of decomposed _s phase and decomposed _s in the interdendritic regions. A series of complex phase transformations was observed. Both decompositions of _s and _s were determined at an early stage of ageing and a four-phase transformation, _f+T+, was observed at the boundaries of _f phase and the phase, instead of clearly observed at the boundaries of phase, in a solution-treated Zn-Al alloy during prolonged ageing.     

Effects of non-steady state nucleation in the kinetics of crystallization of the amorphous alloy Fe80B20

The existence of non-steady state nucleation in the isothermal crystallization of the amorphous alloy Fe₈₀B₂₀ is shown. The incubation time ₀ of isothermal volume crystallization of the alloy is investigated in a wide range as a function of temperature. From these data an equation for the temperature dependence of the viscosity = (T) is derived: = 6.62 exp (2526.97T) × exp [836.52/(T – 530)].[/p]     

Critical Behavior of 2 − Ions in Argon Gas

We measured the drift mobility of ₂ ^– ions in argon gas close to the critical point for (1.005 < T/T _c < 1.04) above T _c 150.7 K in the density range (0.025 <N/N _c< 1.733) around the bulk critical density N _c = 8.08 atoms . nm^–3. The density-normalized zero-field mobility ₀ N of the ions shows a deep minimum as a function of the gas density N as T T _c ⁺ This anomalous reduction of ₀ N occurs at a density N _m 0.76N _c. We believe that this behavior is due to the strong electrostriction exerted by the ion on the highly compressible gas. By introducing suitable contributions to the effective ion radius R due to the large gas compressibility and taking into account short-range local density and viscosity augmentation due to electrostriction, the hydrodynamic Stokes formula ₀ = e/6R, where is the gas viscosity, is in good agreement with the experimental data.     

Separable functions in load separation for the ηpl and ηpl CMOD plastic factors estimation

The objective of this paper is to develop the load separation method for evaluating the _pl and _pl ^CMOD plastic factors used in the J estimation approach based on load versus displacement records. Appropriate forms for the geometry and deformation functions have been suggested from the EPRI Handbook solutions to produce the separable form for the load. The obtained functions are applied to evaluate the _pl and _pl ^CMOD plastic factors for center cracked tension specimen. The present load separation method gave results which are somewhat different from the estimated values of _pl given in the literature. For shallow cracks, the _pl and ^CMOD _pl plastic factors show considerable variation with crack size and the strain hardening exponent. For a deeply cracked CCT specimen, the ^CMOD _pl factor tends to the _pl factor and equals approximately unity. Abbreviations: CCT – center cracked specimen; CMOD – Crack Mouth Opening Displacement; EPRI – Electric Power Research Institute; FEM – Finite Element Method; LLD – Load Line Displacement.     

Ultrasound attenuation peaks due to order parameter collective modes in impure superconductors with strong electron-hole asymmetry

First we develop the general theory of order parameter collective modes and their sound attenuation in dirty superconductors with even-parity pairing. For the BCS state the theory yields a peak in the attenuation of longitudinal ultrasound whose height is proportional to _o[J(s)]² and whose positiont _p 1 -T _p /T _c and width are proportional to _o/ = 2sql(s=v _s /v _F ). Here, is the particle-hole asymmetry parameter,J(s) a given function, _o the sound frequency, and =1/2 _N ^–1 the normal state scattering rate. For T _c /T _o 0.1,s=2/3, and _o/ ql « 1, the peak agrees qualitatively with the peaks found earlier for the polar and axial states. However, at low temperatures the attenuation goes like T⁴ instead ofT ³. We conclude that the observed longitudinal ultrasound peaks in the heavy-fermion superconductors UBe₁₃ and UPt₃ can be explained equally well by any anisotropic state having nodes in the gap, in particular, by a conventional state whose gap has the full symmetry of the Fermi surface. The nature of the state is reflected mainly in the dependence of the peak on the direction of sound propagation.     

Effect of the elastic factor on the hydrodynamic stability of a structurally viscous medium

The effect of relaxation phenomena on the hydrodynamic stability of the plane gradient flow of a structurally viscous medium is investigated using linear theory.Notation _ij stress tensor deviator - U_i components of the velocity vector - x_i coordinates - t time - P pressure - =₀L/^*V plasticity parameter - _o limiting shear stress - andc dimensionless wave number and the perturbation frequency - Re=VL/* Reynolds number - density - F_ij deformation rate tensor Translated from Inzhenerno-Fizicheskii Zhurnal, Vol. 35, No. 5, pp. 868–871, November, 1978.     

Coupling of order parameter collective modes to spin density in3He-B

We derive a general expression for the dynamic spin susceptibility of³He-B which is valid for all magnetic fields. The coupling of real and imaginary modes by particle-hole asymmetry is taken into account. Then we calculate the contribution of the mode at frequency =2 – 1/4 ( is the effective Larmor frequency) to the transverse susceptibility. The spectral weight of this mode in magnetic resonance absorption is proportional to (/)^1/2 (–)², where and are particle-hole asymmetry parameters. From the experimental coupling strength of the real squashing mode to sound we estimate (–)²10^–4. The dynamic susceptibility satisfies the sum rules of Leggett. Finally we point out the difficulties in calculating the transverse NMR frequency of³He-B. These difficulties arise from theS _z =0 Cooper pairs and from the coupling ofJ _z =±1 modes forJ=1 andJ=2.     

Shear viscosity of the B phase of superfluid3He

The shear viscosity (T) in the Balian-Werthamer (BW) state of superfluid ³ He is calculated variationally throughout the region 0t 1(t=T/T _c) from the transport equation for Bogoliubov quasiparticles. Coherence factors are treated exactly in the calculation of the collision integral. The numerical result for =_s= _s(T)/_n(T_c) agree very well with experiment in the range 0.8t1.0. Analytic expressions = 0.577 (1–1.0008t) and =1–(23/64) [=(T)/k _B T] are obtained in the low-temperature region and in the vicinity ofT _c, respectively. From the numerical analysis it is shown that the latter equation is valid only in the temperature range 0.9997t1.0.Supported by the Research Institute for Fundamental Physics, Kyoto University.     

Shear viscosity measurements of liquid 3He-4He mixtures above 0.5 K

Simultaneous measurements of () and of the molar volume are reported for liquid mixtures of ³He in ⁴He over the temperature range between 0.5 and 2.5 K. Here is the shear viscosity and is the mass density. In the superfluid phase, the product of the normal components, _n and _n , is measured. The mixtures with ³He molefractions 0.30 < X < 0.80 are studied with emphasis on the region near the superfluid transition T and near the phase-separation curve. Along the latter, they are compared with data by Lai and Kitchens. For X > 0.5, the viscosity singularity near T becomes a faint peak, which however fades into the temperature-dependent background viscosity as X tends to the tricritical concentration X _t. Likewise, no singularity in is apparent when T _t is approached along the phase separation branches _– and ₊. Furthermore, viscosity data are reported for ³He and compared with previous work. Finally, for dilute mixtures with 0.01 X 0.05, the results for are compared with previous data and with predictions.     

Second-order wave force on a large cylinder

Summary A new technique has been developed for estimating the wave loadings on large circular cylinders. The theory, mainly due to Lighthill, has been applied to the case of large cylinders. Comparison with the previously reported experimental results shows favourable agreement for the range of parameters indicated in the graphs.Notations The following symbols are used in this paper a Wave amplitude - b Radius of the cylinder - D Diameter of the cylinder - F _l Linear force - F _d Dynamic force - F _q Quadratic force - F _w Water line force - g Acceleration due to gravity - h Depth of water - H=2a Total wave height - k Wave number,k=2/L - L Wave length - n Outward normal to body surface - n _x Direction cosine between the normaln and the given force - S Body surface - t Time - T Wave period - x, y, z Cartesian coordinates - r, ,z Cylindrical coordinates - W Vertical velocity of fluid - Density of fluid - Wave frequency, =2/T - Total velocity potential - _r , _z Partial derivatives with respect tor andz, respectively - ⁽¹⁾, _l Linear theory velocity potential - ⁽²⁾, _q Second-order velocity potential - Water surface elevation or wave height - _i Linear theory wave height - _q Quadratic theory wave height With 3 Figures     

Hydrodynamic similarity in an oscillating-body viscometer

Hydrodynamic similarity can be used to calibrate simply and accurately an oscillating-body viscometer of arbitrarily complicated geometry. Usually, an explicit hydrodynamic model based on a simple geometry is required to deduce viscosity from the transfer function of an oscillating body such as a vibrating wire or a quartz torsion crystal. However, at low Reynolds numbers the transfer function of any immersed oscillator depends on the fluid's viscosity only through the viscous penetration depth(2/)^1/2. (Here and are the fluid's viscosity and density and/2 is the oscillator's frequency.) This hydrodynamic similarity can be exploited if the oscillator is overdamped and thus is sensitive to viscosity in a broad frequency range. Even an oscillator of poorly known geometry can be characterized over a range of penetration depths by measurements in a fluid of known and over the corresponding range of frequencies. The viscosity of another fluid can then be compared to that of the calibrating fluid with high accuracy by varying the frequency so that the penetration depth falls within the characterized range. In the present work, hydrodynamic similarity was demonstrated with a highly damped viscometer comprised of an oscillating screen immersed in carbon dioxide. The fluid's density was varied between 2 and 295 kg·m^–3 and the fluid's temperature was varied between 25 and 60°C. The corresponding variation of the viscosity was 50%.Paper presented at the Twelfth Symposium on Thermophysical Properties, June 19–24, 1994, Boulder, Colorado, U.S.A.     

On the dispersion of gases under the action of the medium

Some general regularities of dispersion of a gas emerging from a nozzle submerged in a liquid are considered. A condition for establishment of the so-called maximum dispersion state is formulated.Notation ₀ coefficient of surface tension at the liquidgas boundary - contact angle of wetting of the nozzle material surface by the liquid - p_at atmospheric pressure - p air pressure - density of the liquid - g gravitational acceleration - h height of the liquid column - ₁, and _g dynamic viscosity coefficients of the liquid and gas, respectively - R and r radii of the bubble and nozzle, respectively - Q and F dimensionless criteria - , , , , and undetermined coefficients - ratio of the circumference of a circle to its diameter