Second order phase equilibrium for classical bodies

Summary In this note the Author, recalling a previous work[15], gives a new formulation of second order phase equilibria for classical bodies such as those defined by Truesdell and Toupin in[8].The Author arrives at three equivalent systems of partial differential equations (generalized Ehrenfest equations), the conditions for whose integration are shown to be always satisfied.Finally, as particular cases, the equations ruling the phase equilibria for classical fluids and for n component-classical fluid mixtures are given.

---

Sommario In questa nota l'Autore, rifacendosi ad un lavoro precedente[15], presenta una nuova formulazione degli equilibri di fase del secondo ordine per corpi classici come quelli definiti da Truesdell e Toupin in[8].L'Autore perviene a tre sistemi equivalenti di equazioni alle derivate parziali (equazioni di Ehrenfest generalizzate) dei quali viene dimostrata la integrabilità.Infine, come casi particolari, si ottengono le equazioni che governano l'equilibrio di fase per fluidi classici e per miscele fluide classiche ad n componenti.

---

---

This work was supported by the Gruppo Nazionale per la Fisica Matematica of C.N.R.     

Deformed particle image pattern matching in Particle Image Velocimetry

The deformation of particle image patterns due to velocity gradients causes errors of velocity measurements and false velocity detections in PIV (Particle Image Velocimetry). A novel technique to overcome those limitations inherent in the conventional PIV by correcting the particle image pattern according to the local velocity gradients in two dimensional flows, i.e. u/x, u/y, v/x and v/y, is proposed and successfully applied to a water flow downstream of a backward facing step.     

Theoretical derivation of tangential velocity profiles in a flat vortex chamber-influence of turbulence and wall friction

Summary As part of a study on the hydrodynamics of a cyclone separator, a theoretical investigation of the flow pattern in a flat box cyclone (vortex chamber) has been carried out. Expressions have been derived for the tangential velocity profile as influenced by internal friction (eddy viscosity) and wall friction. The most important parameter controlling the tangential velocity profile is = –u ₀ R/(v+ ), where u ₀ is the radial velocity at the outer radius R of the cyclone, the kinematic liquid viscosity and is the kinematic eddy viscosity. For values of greater than about 10 the tangential velocity profile is nearly hyperbolic, for smaller than 1 the tangential velocity even decreases towards the centre. It is shown how and also the wall friction coefficient may be obtained from experimental velocity profiles with the aid of suitable graphs. Because of the close relation between eddy viscosity and eddy diffusion, measurements of velocity profiles in flat box cyclones will also provide information on the eddy motion of particles in a cyclone, a motion reducing its separation efficiency.List of symbols A cross-sectional area of cyclone inlet - h height of cyclone - p static pressure in cyclone - p static pressure difference in cyclone between two points on different radius - r radius in cyclone - r ₁ radius of cyclone outlet - R radius of cyclone circumference - u radial velocity in cyclone - u ₀ radial velocity at circumference of flat box cyclone - v tangential velocity - v ₀ tangential velocity at circumference of flat box cyclone - w axial velocity - z axial co-ordinate in cyclone - friction coefficient in flat box cyclone (for definition see § 5) - ₁ value of friction coefficient for ₁<< ₂ - ₂ value of friction coefficient for ₂<<1 - = - ₁ value of for ₁<< ₂ - ₂ value of for ₂<<1 - thickness of laminar boundary layer - =/h - turbulent kinematic viscosity - ratio of z to h - k ratio of height of cyclone to radius R of cyclone - parameter describing velocity profile in cyclone =–u ₀ R/(+) - kinematic viscosity of fluid - density of fluid - ratio of r to R - ₁ value of at outlet of cyclone - ₂ value of at inner radius of cyclone inlet - _w shear stress at cyclone wall - angular momentum in cyclone/angular momentum in cyclone inlet - ₁ value of at = ₁ - ₂ value of at = ₂     

Exact Solutions of the Hydrodynamic Equations Derived from Partially Invariant Solutions

The paper proposes a heuristic approach to constructing exact solutions of the hydrodynamic equations based on the specificity of these equations. A number of systems of hydrodynamic equations possess the following structure: they contain a reduced system of n equations and an additional equation for an extra function w. In this case, the reduced system, in which w = 0, admits a Lie group G. Taking a certain partially invariant solution of the reduced system with respect to this group as a seed:rdquo; solution, we can find a solution of the entire system, in which the functional dependence of the invariant part of the seed solution on the invariants of the group G has the previous form. Implementation of the algorithm proposed is exemplified by constructing new exact solutions of the equations of rotationally symmetric motion of an ideal incompressible liquid and the equations of concentrational convection in a plane boundary layer and thermal convection in a rotating layer of a viscous liquid.     

The Thermodynamic Significance of the Local Volume Averaged Temperature

Since the temperature is not an additive function, the traditional thermodynamic point of view suggests that the volume integral of the temperature has no precise physical meaning. This observation conflicts with the customary analysis of non-isothermal catalytic reactors, heat pipes, driers, geothermal processes, etc., in which the volume averaged temperature plays a crucial role. In this paper we identify the thermodynamic significance of the volume averaged temperature in terms of a simple two-phase heat transfer process. Given the internal energy as a function of the point temperature and the density

On the Propagation of Long‐Wave Perturbations in a Two‐Layer Free‐Boundary Rotational Fluid

A mathematical model for the propagation of longwave perturbations in a freeboundary shear flow of an ideal stratified twolayer fluid is considered. The characteristic equation defining the velocity of perturbation propagation in the fluid is obtained and studied. The necessary hyperbolicity conditions for the equations of motion are formulated for flows with a monotonic velocity profile over depth, and the characteristic form of the system is calculated. It is shown that the problem of deriving the sufficient hyperbolicity conditions is equivalent to solving a system of singular integral equations. The limiting cases of weak and strong stratification are studied. For these models, the necessary and sufficient hyperbolicity conditions are formulated, and the equations of motion are reduced to the Riemann integral invariants conserved along the characteristics.     

Some characteristic features of diffusion of a magnetic field into a moving conductor

The study of the diffusion of a magnetic field into a moving conductor is of interest in connection with the production of ultra-high-strength magnetic fields by rapid compression of conducting shells [1,2]. In [3,4] it is shown that when a magnetic field in a plane slit is compressed at constant velocity, the entire flux enters the conductor. In the present paper we formulate a general result concerning the conservation of the sum current in the cavity and conductor for arbitrary motion of the latter. We also consider a special case of conductor motion when the flux in the cavity remains constant despite the finite conductivity of the material bounding the magnetic field.Notation ₁, * flux which has diffused into the conductor - ₂ flux in the cavity - ₀ sum flux - r radius - r* cavity boundary - thickness of the skin layer - (r) delta function of r - t time - q intensity of the fluid sink - v velocity - flux which has diffused to a depth larger than r - x self-similar variable - dimensionless fraction of the flux which has diffused to a depth larger than r - * fraction of the flux which has diffused into the conductor - a conductivity - c electrodynamic constant - R_m magnetic Reynolds number - dimensionless parameter     

Optical co-ordinates in general relativity 2 — the problem of distance

Summary Let denote the congruence of null geodesics associated with a given optical observer inV ₄. We prove that determines a unique collection of vector fieldsM₍₎ ( =1, 2, 3) and ₍₀₎ overV ₄, satisfying a weak version of Killing's conditions.This allows a natural interpretation of these fields as the infinitesimal generators of spatial rotations and temporal translation relative to the given observer. We prove also that the definition of the fieldsM₍₎ and ₍₀₎ is mathematically equivalent to the choice of a distinguished affine parameter f along the curves of, playing the role of a retarded distance from the observer.The relation between f and other possible definitions of distance is discussed.

---

Sommario Sia la congruenza di geodetiche nulle associata ad un osservatore ottico assegnato nello spazio-tempoV ₄. Dimostriamo che determina un'unica collezione di campi vettorialiM₍₎ ( =1, 2, 3) e ₍₀₎ inV ₄ che soddisfano una versione in forma debole delle equazioni di Killing. Ciò suggerisce una naturale interpretazione di questi campi come generatori infinitesimi di rotazioni spaziali e traslazioni temporali relative all'osservatore assegnato. Dimostriamo anche che la definizione dei campiM₍₎, ₍₀₎ è matematicamente equivalente alla scelta di un parametro affine privilegiato f lungo le curve di, che gioca il ruolo di distanza ritardata dall'osservatore. Successivamente si esaminano i legami tra f ed altre possibili definizioni di distanza in grande.

---

---

Work performed in the sphere of activity of: Gruppo Nazionale per la Fisica Matematica del CNR.     

Critical heat fluxes in boiling of liquid nitrogen underheated to saturation point in forced-flow conditions

This paper gives the results of experimental determinations of the critical heat fluxes in the boiling of Liquid nitrogen in forced-flow conditions in the mass velocity range 2 · 10³-40 · 10³ kg/m² · sec, pressure range 29 · 10⁴–245 · 10⁴ N/m², and at underheatings corresponding to the onset of normal boiling crises.Notation q₀ critical heat flux - r heat of vaporization - i enthalpy of flow corresponding to saturation point - i enthalpy of flow corresponding to liquid temperature - surface tension - density of liquid - density of saturated vapor - C _f friction factor - W_g mass velocity - Fr_* Froude number - g acceleration due to gravity     

Supersonic air flow past a sphere with account for vibrational relaxation

A numerical solution is obtained for the problem of air flow past a sphere under conditions when nonequilibrium excitation of the vibrational degrees of freedom of the molecular components takes place in the shock layer. The problem is solved using the method of [1]. In calculating the relaxation rates account was taken of two processes: 1) transition of the molecular translational energy into vibrational energy during collision; 2) exchange of vibrational energy between the air components. Expressions for the relaxation rates were computed in [2]. The solution indicates that in the state far from equilibrium a relaxation layer is formed near the sphere surface. A comparison is made of the calculated values of the shock standoff with the experimental data of [3].Notation uV_max, vV_max velocity components normal and tangential to the sphere surface - V_max maximal velocity - P V _max ² pressure - density - TT temperature - e_viRT vibrational energy of the i-th component per mole (i=–O₂, N₂) - =rb–1 shock wave shape - a _f the frozen speed of sound - HRT/m gas total enthalpy     

The dynamic performance of the Weissenberg rheogoniometer III. Large amplitude oscillatory shearing — harmonic analysis

The harmonic content of the nonlinear dynamic behaviour of 1% polyacrylamide in 50% glycerol/water was studied using a standard Model R 18 Weissenberg Rheogoniometer. The Fourier analysis of the Oscillation Input and Torsion Head motions was performed using a Digital Transfer Function Analyser.In the absence of fluid inertia effects and when the amplitude of the (fundamental) Oscillation Input motion _I is much greater than the amplitudes of the Fourier components of the Torsion Head motion _Tn empirical nonlinear dynamic rheological propertiesG _n (, ⁰),G _n (, ⁰) and/or _n (, ⁰), _n (, ⁰) may be evaluated without a-priori-knowledge of a rheological constitutive equation. A detailed derivation of the basic equations involved is presented.Cone and plate data for the third harmonic storage modulus (dynamic rigidity)G ₃ (, ⁰), loss modulusG ₃ (, ⁰) and loss angle ₃ (, ⁰) are presented for the frequency range 3.14 × 10^–2 1.25 × 10² rad/s at two strain amplitudes, _CP ⁰ = 2.27 and 4.03. Composite cone and plate and parallel plates data for both the third and fifth harmonic dynamic viscosities ₃ (, ⁰), _S (, ⁰) and dynamic rigiditiesG ₃ (, ⁰),G ₅ (, ⁰) are presented for strain amplitudes in the ranges 1.10 _CP ⁰ 4.03 and 1.80 _PP ⁰ 36 for a single frequency, = 3.14 × 10^–1 rad/s. Good agreement was obtained between the results from both geometries and the absence of significant fluid inertia effects was confirmed by the superposition of the data for different gap widths.     

Effects of MHD free convection and mass transfer on the flow past an oscillating infinite coaxial vertical circular cylinder

The effects of MHD free convection and mass transfer are taken into account on the flow past oscillating infinite coaxial vertical circular cylinder. The analytical expressions for velocity, temperature and concentration of the fluid are obtained by using perturbation technique.

---

Einwirkungen von freier MHD-Konvektion und Stoffübertragung auf eine Strömung nach einem schwingenden unendlichen koaxialen vertikalen Zylinder

Zusammenfassung Die Einwirkungen der freien MHD-Konvektion und Stoffübertragung auf eine Strömung nach einem schwingenden, unendlichen, koaxialen, vertikalen Zylinder wurden untersucht. Die analytischen Ausdrücke der Geschwindigkeit, Temperatur und Fluidkonzentration sind durch die Perturbationstechnik erhalten worden.

---

Nomenclature C _p Specific heat at constant temperature - C the species concentration near the circular cylinder - C _w the species concentration of the circular cylinder - C the species concentration of the fluid at infinite - * dimensionless species concentration - D chemical molecular diffusivity - g acceleration due to gravity - Gr Grashof number - Gm modified Grashof number - K thermal conductivity - Pr Prandtl number - r _a ,r _b radius of inner and outer cylinder - a, b dimensionless inner and outer radius - r,r coordinate and dimensionless coordinate normal to the circular cylinder - Sc Schmidt number - t time - t dimensionless time - T temperature of the fluid near the circular cylinder - T _w temperature of the circular cylinder - T temperature of the fluid at infinite - u velocity of the fluid - u dimensionless velocity of the fluid - U ₀ reference velocity - z,z coordinate and dimensionless coordinate along the circular cylinder - coefficient of volume expansion - * coefficient of thermal expansion with concentration - dimensionless temperature - H ₀ magnetic field intensity - coefficient of viscosity - _e permeability (magnetic) - kinematic viscosity - electric conductivity - density - M Hartmann number - dimensionless skin-friction - frequency - dimensionless frequency     

Binary laminar boundary layer on a vertical surface in the presence of combined free and forced convection

Heat and mass transfer at a vertical surface is examined in the case of combined free and forced convection. The boundary layer equations, transformed to ordinary differential equations, contain a parameter that determines the effect of free convection on the forced motion. Criteria are offered for differentiating the free-convection, forced-convection, and combined regimes.Notation x, y coordinates - u, v velocity components - g acceleration of gravity - T temperature - kinematic viscosity - coefficient of thermal expansion - a thermal diffusivity - ₁ partial vapor density - D diffusion coefficient - W₂ mass velocity of air - independent variable - _w shear stress at wall - thermal conductivity - r latent heat of phase transition - , dimensionless temperature and partial vapor density - m^* the complex (m ₁–m _1w )/(1–m(_1w ) - c_p specific heat at constant pressure - G Grashof number - R Reynolds number - P Prandtl number - S Schmidt number     

Analytical Solution for the Effect of Viscous Dissipation on Mixed Convection in Saturated Porous Media

A new analytical solution is introduced for the effect of viscous dissipation on mixed convection flow and heat transfer about an isothermal vertical wall embedded in Darcy and non-Darcy porous media with uniform free stream velocity. The effect of viscous dissipation on mixed convection in both regimes has been analyzed for both the aiding and opposing flows using Gebhart number, Ge _x =gx/c _p. The governing parameters are Re, Ra, Pe and Ge _x . The case of Re=0 corresponds to Darcy mixed convection region and Re/Pe is identified as the mixed convection governing parameter, Ra=0 leading to pure forced convection. A good agreement was found between the numerical and analytical solutions. It was found from the Nusselt number results that viscous dissipation lowers the heat transfer rate in both Darcy and Forchheimer flow regimes for aiding as well as opposing flows.     

Asymptotic Modeling of Nonlinear Wave Processes in Shock‐Loaded Elastoplastic Materials

Nonlinear wave processes in shockloaded elastoplastic materials are modeled. A comparison of the results obtained with experimental data and numerical solutions of exact systems of dynamic equations shows that the model equations proposed qualitatively describe the stressdistribution evolution in both the elasticflow and plasticflow regions and can be used to solve one and twodimensional problems of pulsed deformation and fracture of elastoplastic media.     

Axisymmetric spreading of a heavy liquid over a plane

A study is made of the steady flow over a horizontal plane of a heavy inviscid incompressible liquid which flows through the side surface of a circular cylinder which rises above the plane to height h and has a base radius ofa. The motion of the liquid is assumed to be symmetric with respect to the axis of the cylinder; the pressure p is constant (equal to the atmospheric pressure) on the free surface of the liquid. Fora/h = 1, this problem can be regarded as a problem of perturbation of the flow from a flat source by a free surface. Investigation showed that this perturbation problem is essentially nonlinear, and a solution of it in the complete region occupied by the liquid can be obtained only in variables of the boundary layer type. The problem admits linearization under the additional assumption that the parameter = Q²/(8²ga³) is small; here, Q is the constant volume flow rate of the liquid per unit height of the cylinder, and g is the acceleration of free fall. For the case 1, 1 the problem is solved by the method of integral transformations. A noteworthy feature of the solution is the slow damping of the perturbations of the velocity with the depth (inversely proportional to the square of the distance from the free surface), in contrast to the similar problem of the wave motions of a heavy liquid, for which the velocity perturbations are damped exponentially.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 3–7, March–April, 1984.     

The contribution of internal degrees of freedom to the non-Newtonian rheology of model polymer fluids

We report non-equilibrium molecular dynamics simulations of rigid and non-rigid dumbbell fluids to determine the contribution of internal degrees of freedom to strain-rate-dependent shear viscosity. The model adopted for non-rigid molecules is a modification of the finitely extensible nonlinear elastic (FENE) dumbbell commonly used in kinetic theories of polymer solutions. We consider model polymer melts — that is, fluids composed of rigid dumbbells and of FENE dumbbells. We report the steady-state stress tensor and the transient stress response to an applied Couerte strain field for several strain rates. We find that the rheological properties of the rigid and FENE dumbbells are qualitatively and quantitatively similar. (The only exception to this is the zero strain rate shear viscosity.) Except at high strain rates, the average conformation of the FENE dumbbells in a Couette strain field is found to be very similar to that of FENE dumbbells in the absence of strain. The theological properties of the two dumbbell fluids are compared to those of a corresponding fluid of spheres which is shown to be the most non-Newtonian of the three fluids considered.Symbol Definition b dimensionless time constant relating vibration to other forms of motion - F force on center of mass of dumbbell - F _i force on bead i of dumbbell - F force between center of masses of dumbbells and - F _ij force between beads i and j - h vector connecting bead to center of mass of dumbbell - H dimensionless spring constant for dumbbells, in units of / ² - I moment of inertia of dumbbell - J general current induced by applied field - k _B Boltzmann's constant - L angular momentum - m mass of bead, (= m/2) - M mass of dumbbell, g - N number of dumbbells in simulation cell - P translational momentum of center of mass of dumbbell - P pressure tensor - P _xy xy component of pressure tensor - Q separation of beads in dumbbell - Q _eq equilibrium extension of FENE dumbbell and fixed extension of rigid dumbbell - Q ₀ maximum extension of dumbbell - r _ij vector connecting beads i and j - r position vector of center of mass dumbbell - R vector connecting centers of mass of two dumbbells - t time - t ^* dimensionless time, in units of m/ - T ^* dimensionless temperature, in units of /k - u potential energy - u velocity vector of flow field - u _x x component of velocity vector - V volume of simulation cell - X general applied field - strain rate, s^–1 - ^* dimensionless shear rate, in units of /m ² - general transport property - Lennard-Jones potential well depth - friction factor for Gaussian thermostat - shear viscosity, g/cms - ^* dimensionless shear viscosity, in units of m/ ² - ^* dimensionless number density, in units of ^–3 - Lennard-Jones separation of minimum energy - relaxation time of a fluid - angular velocity of dumbbell - orientation angle of dumbbell      

Effect of Free‐Stream Turbulence on the Flow Structure near a Wedge and the Windward Side of an Airfoil

The flow structure behind wire grids is studied for flows with a low subsonic velocity, and the effect of grids on the boundarylayer flow structure is considered. It is shown that the meanvelocity inhomogeneity induced by the grid does not disappear until a distance of 925 calibers downstream of the grid is reached. Liquidcrystal thermography combined with hotwire measurements made it possible to find the source of steady largescale streamwise vortex structures in the boundary layer on a wedge and on an airfoil and to determine the parameters of these structures.     

Flow of second-order fluids over an enclosed rotating disc

Summary This paper is devoted to a study of the flow of a second-order fluid (flowing with a small mass rate of symmetrical radial outflow m, taken negative for a net radial inflow) over a finite rotating disc enclosed within a coaxial cylinderical casing. The effects of the second-order terms are observed to depend upon two dimensionless parameters ₁ and ₂. Maximum values ₁ and ₂ of the dimensionless radial distances at which there is no recirculation, for the cases of net radial outflow (m>0) and net radial inflow (m<0) respectively, decrease with an increase in the second-order effects [represented by T(=₁+₂)]. The velocities at ₁ and ₂ as well as at some other fixed radii have been calculated for different T and the associated phenomena of no-recirculation/recirculation discussed. The change in flow phenomena due to a reversal of the direction of net radial flow has also been studied. The moment on the rotating disc increases with T.Nomenclature , , z coordinates in a cylindrical polar system - z ₀ distance between rotor and stator (gap length) - =/z ₀, dimensionless radial distance - =z/z ₀, dimensionless axial distance - _s = _s/z₀, dimensionless disc radius - V =(u, v, w), velocity vector - dimensionless velocity components - uniform angular velocity of the rotor - , p fluid density and pressure - P =p/(₂ z ₀₂ ² , dimensionless pressure - ₁, ₂, ₃ kinematic coefficients of Newtonian viscosity, elastico-viscosity and cross-viscosity respectively - ₁, _{2 2}/z ₀ ² , resp. ₃/z ₀ ² , dimensionless parameters representing the ratio of second-order and inertial effects - m = , mass rate of symmetrical radial outflow - l a number associated with induced circulatory flow - Rm =m/(z ₀₁), Reynolds number of radial outflow - R _l =l/(z ₀₁), Reynolds number of induced circulatory flow - Rz =z ₀ ² /₁, Reynolds number based on the gap - ₁, ₂ maximum radii at which there is no recirculation for the cases Rm>0 and Rm<0 respectively - ₁(T), ₂(T) ₁ and ₂ for different T - U _1(T) ⁽⁺⁾ = dimensionless radial velocity, Rm>0 - V _1(T) ⁽⁺⁾ = , dimensionless transverse velocity, Rm>0 - U _2(T) ^(–) = , dimensionless radial velocity, Rm=–Rn<0, m=–n - V _2(T) ^(–) = , dimensionless transverse velocity, Rm<0 - C _m moment coefficient     

Saint-Venant's principle in linear two-dimensional elasticity for non-striplike domains

Let be a simply connected domain in the x ₁-x ₂ plane which lies within the strip 0<x ₂, is a simple closed piecewise smooth curve. Let l= [(x ₁, x ₂): (x ₁, x ₂) and x ₁>0], _l = [(x ₁ x ₂): (x ₁ ,x ₂) and x ₁>1>0].Suppose that a two-dimensional homogeneous isotropic elastic body occupies , that a self-equilibrated stress loading is applied to - l, and that l is stress-free. Knowles [2] and Flavin [6] showed that the elastic energy in _l decays exponentially with respect to l with an exponential decay constant of the form k/b, where k is a universal constant. It is shown here that a decay constant of the form c/ may be obtained where c is a universal constant and is a characteristic dimension of , which is more appropriate than b for general non-striplike domains. In addition, an appropriate decay theorem is obtained for coil-like domains.