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.     

The mirror relations and nonlinear viscoelasticity of polymer melts

An attempt is made to incorporate into a quasilinear viscoelastic constitutive equation of the Boltzmann superposition type the two mirror relations of Gleissle, as well as his relation between the steady-state first normal-stress difference and the shear viscosity curve. It is shown that the three relations can hold separately within this constitutive model, but not simultaneously, because they require a different nonlinear strain measure, namelyS ₁₂ () = – a ( – 1) (a = 0 for 1,a = 1 for 1) for the mirroring of the viscosities,S ₁₂ () = – a (–k ²/) (a = 0 for k, a = 1 for k) for the mirroring of the first normal-stress coefficients, and for the third relation. Here denotes the shear strain and erf the error function. Experimental data on melts of a low-density polyethylene, a high-density polyethylene and a polypropylene show that the mirror relations are passable approximations, but that the third relation meets reality surprisingly close if the right value ofk is used.     

On the Ellipticity of the Middle Surface of a Shell and its Application to the Asymptotic Analysis of ‘Membrane’ Shells

We consider a surface S = (), where ² is a bounded, connected, open set with a smooth boundary and : ³ is a smooth map; let () denote the components of the two-dimensional linearized strain tensor of S and let 0 with length 0 > 0. We assume the the norm ,|| ()||0, in the space V0() = { H¹() × H¹() × L²(); = 0 on 0 } is equivalent to the usual product norm on this space. We then establish that this assumption implies that the surface S is uniformly elliptic and that we necessarily have 0 = .     

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.     

Nonlinear rheological behavior of a concentrated spherical silica suspension

Linear and nonlinear viscoelastic properties were examined for a 50 wt% suspension of spherical silica particles (with radius of 40 nm) in a viscous medium, 2.27/1 (wt/wt) ethylene glycol/glycerol mixture. The effective volume fraction of the particles evaluated from zero-shear viscosities of the suspension and medium was 0.53. At a quiescent state the particles had a liquid-like, isotropic spatial distribution in the medium. Dynamic moduli G^* obtained for small oscillatory strain (in the linear viscoelastic regime) exhibited a relaxation process that reflected the equilibrium Brownian motion of those particles. In the stress relaxation experiments, the linear relaxation modulus G(t) was obtained for small step strain (0.2) while the nonlinear relaxation modulus G(t, ) characterizing strong stress damping behavior was obtained for large (>0.2). G(t, ) obeyed the time-strain separability at long time scales, and the damping function h() (–G(t, )/G(t)) was determined. Steady flow measurements revealed shear-thinning of the steady state viscosity () for small shear rates (< ^–1; = linear viscoelastic relaxation time) and shear-thickening for larger (>^–1). Corresponding changes were observed also for the viscosity growth and decay functions on start up and cessation of flow, + (t, ) and _– (t, ). In the shear-thinning regime, the and dependence of ₊(t,) and _–(t,) as well as the dependence of () were well described by a BKZ-type constitutive equation using the G(t) and h() data. On the other hand, this equation completely failed in describing the behavior in the shear-thickening regime. These applicabilities of the BKZ equation were utilized to discuss the shearthinning and shear-thickening mechanisms in relation to shear effects on the structure (spatial distribution) and motion of the suspended particles.Dedicated to the memory of Prof. Dale S. Parson     

The value of the critical exponent for reaction-diffusion equations in cones

Let D R ^N be a cone with vertex at the origin i.e., D = (0, )x where S ^N–1 and x D if and only if x = (r, ) with r=¦x¦, . We consider the initial boundary value problem: u _t = u+u ^p in D×(0, T), u=0 on Dx(0, T) with u(x, 0)=u ₀(x) 0. Let ₁ denote the smallest Dirichlet eigenvalue for the Laplace-Beltrami operator on and let ₊ denote the positive root of (+N–2) = ₁. Let p ^* = 1 + 2/(N + ₊). If 1 < p < p ^*, no positive global solution exists. If p>p ^*, positive global solutions do exist. Extensions are given to the same problem for u _t=+¦x¦ u ^p .This research was supported in part by the Air Force Office of Scientific Research under Grant # AFOSR 88-0031 and in part by NSF Grant DMS-8 822 788. The United States Government is authorized to reproduce and distribute reprints for governmental purposes not withstanding any copyright notation therein.     

Symmetries and Global Solvability of the Isothermal Gas Dynamics Equations

We study the Cauchy problem associated with the system of two conservation laws arising in isothermal gas dynamics, in which the pressure and the density are related by the -law equation p() with =1. Our results complete those obtained earlier for >1. We prove the global existence and compactness of entropy solutions generated by the vanishing viscosity method. The proof relies on compensated compactness arguments and symmetry group analysis. Interestingly, we make use here of the fact that the isothermal gas dynamics system is invariant modulo a linear scaling of the density. This property enables us to reduce our problem to that with a small initial density.One symmetry group associated with the linear hyperbolic equations describing all entropies of the Euler equations gives rise to a fundamental solution with initial data imposed on the line =1. This is in contrast to the common approach (when >1) which prescribes initial data on the vacuum line =0. The entropies we construct here are weak entropies, i.e., they vanish when the density vanishes.Another feature of our proof lies in the reduction theorem, which makes use of the family of weak entropies to show that a Young measure must reduce to a Dirac mass. This step is based on new convergence results for regularized products of measures and functions of bounded variation.Acknowledgement P.G.L. and V.S. were supported by a grant from INTAS (01-868). The support and hospitality of the Isaac Newton Institute for Mathematical Sciences, University of Cambridge, where part of this research was performed during the Semester Program Nonlinear Hyperbolic Waves in Phase Dynamics and Astrophysics (January to July 2003) is also gratefully acknowledged. P.G.L. was also supported by the Centre National de la Recherche Scientifique (CNRS).     

Construction of exact numerical solutions of the stationary traveling wave type for viscous thin films

An effective numerical procedure, based on the Galerkin method, for finding solutions of the stationary traveling wave type in the complete formulation is proposed for the case of viscous liquid films. Examples of a viscous film flowing freely down a vertical surface have been calculated. The calculations have been made for various values of the dimensionless surface tension , including =0. The method makes it possible to predict a number of bifurcations that occur as decreases. The existence of numerous families of stationary traveling waves when 1 was demonstrated in [6]. The present study shows that as 1 all but one of these families of wave solutions disappear. The shape of the periodic and solitary waves and the pressure distribution in the film are found for various . When =0 and the wave number is fairly small, the periodic solution has a singularity, as predicted in [14]: at the crest of the wave a corner point appears; the angle between the tangents at this point =140–150. The method proposed can be used to calculate other wavy film flows.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 94–100, May–June, 1990.     

An example of a quasiconvex function that is not polyconvex in two dimensions

We study the different notions of convexity for the function f () = ||² (||² – 2 det ) where ^2×2, introduced by Dacorogna & Marcellini. We show that f is convex, polyconvex, quasiconvex, rank-one convex, if and only if ¦¦ 2/3 2, 1, 1+ (for some >0), 2/3, respectively.     

On the damping function of shear relaxation modulus for entangled polymers

Published data of the damping function of the shear relaxation modulus, h(), are reviewed. This is the ratio of the relaxation modulus measured at a finite magnitude of shear, , to that at the limit of = 0. Majority of the data are in accord with the universal function of the Doi-Edwards tube model theory, in which the damping or the decrease of h() is attributed to the contraction along the tube of extended polymer chains. The weaker damping seems to be attributed to 1) comb-branching such as in LDPE; 2) lack of entanglement in too short chains; 3) bimodal molecular weight distribution. However, a star-branching does not cause a deviation from the tube model theory and a broadness of molecular weight distribution is not a major origin of a weaker damping. A star-branched polystyrene with 15 arms exhibits no strain dependence: h() = 1. For highly entangled systems with more than 50 entanglement points per molecule, the strain dependence is stronger than that of the Doi-Edwards theory. This could be due to a slip or an instability of deformation in the material.     

Shear softening and thixotropic properties of wheat flour doughs in dynamic testing at high shear strain

Shear softening and thixotropic properties of wheat flour doughs are demonstrated in dynamic testing with a constant stress rheometer. This behaviour appears beyond the strictly linear domain (strain amplitude ₀ 0.2%),G,G and |^*| decreasing with ₀, the strain response to a sine stress wave yet retaining a sinusoidal shape. It is also shown thatG recovers progressively in function of rest time. In this domain, as well as in the strictly linear domain, the Cox-Merz rule did not apply but() and | ^*())| may be superimposed by using a shift factor, its value decreasing in the former domain when ₀ increases. Beyond a strain amplitude of about 10–20%, the strain response is progressively distorted and the shear softening effects become irreversible following rest.     

Representation of the rheological properties of polymer melts in terms of complex fluidity

In dynamic rheological experiments melt behavior is usually expressed in terms of complex viscosity ^* () or complex modulusG ^* (). In contrast, we attempted to use the complex fluidity ^* () = 1/µ ^* () to represent this behavior. The main interest is to simplify the complex-plane diagram and to simplify the determination of fundamental parameters such as the Newtonian viscosity or the parameter of relaxation-time distribution when a Cole-Cole type distribution can be applied. ^* () complex shear viscosity - () real part of the complex viscosity - () imaginary part of the complex viscosity - G ^* () complex shear modulus - G() storage modulus in shear - G() loss modulus in shear - J ^* () complex shear compliance - J() storage compliance in shear - J() loss compliance in shear - shear strain - rate of strain - angular frequency (rad/s) - shear stress - loss angle - ^* () complex shear fluidity - () real part of the complex fluidity - () imaginary part of the complex fluidity - ₀ zero-viscosity - ₀ average relaxation time - h parameter of relaxation-time distribution     

On shock layer method for the hypersonic flow problems

Chernyi’s series method[1] is not proper for the case that(γ-l)/(γ+l)<<2/(γ+1)×M²sin²β (γ=c_p/c_v-adiabatic index number, M-Much number, β-shock incidence). In this paper, we only suppose that in the neighbour of the shock, there exists a shock layer in which the density of the gas is very big, but we do not remove the case that (γ-1)/(γ+1)<<2/(γ+1)M²sin²β.     

Hydrodynamic instability of the contact zone between accelerated gases

An experimental apparatus for investigating Rayleigh-Taylor instability in the transition layer between two gases at accelerations g 10⁵g₀ (g₀ is the acceleration of gravity) is described. The constantly acting acceleration is communicated to the contact zone by the compression wave formed ahead of a flame front. The linear stage of development is investigated together with the effect of the thickness of the contact zone. It is shown that on the interval 0.3 < <- ( is the wavelength of the disturbance at the edge of the contact zone) the rate of growth of the perturbation amplitude 0.5₀, where ₀ is the amplitude growth rate for media separated by an interface with a discontinuous change of density.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 15–21, November–December, 1991.     

An “ideal equivalent gas method” for the study of shock waves in supersonic real gas flows

Summary A method developed by the author for the systematic study of the thermodynamic and dynamic properties of the gas behind a shock wave is reported.The method is applicable to supersonic flow regimes for which the excitation, dissociation and ionization effects invalidate the usually adopted hypothesis of ideal gas.An Ideal Equivalent Gas, having the ratio of the specific heats _s dependent on Mach number and altitude of flight is postulated.On the basis of the mass, momentum and energy conservation equations, valid through the shock wave, the relations defining the thermodynamic and dynamic state of the gas behind the shock wave are derived. These relations establish an extension of the classic relations valid for the ideal gas and reduce to them identically for _s=.The dependence of the ratio of specific heats _s of the Ideal Equivalent Gas on Mach number and altitude has been established, over a wide range, on the basis of the real gas solutions derived by Huber.

---

Sommario Nella presente nota viene esposto un metodo sviluppato dall'autore per lo studio sistematico dello stato termodinamico e dinamico del gas a valle di un'onda d'urto in regime supersonico, allorchè cioè gli effetti dell'eccitazione dei gradi di libertà vibrazionali delle molecole e della loro dissociazione e successiva ionizzazione invalidano l'ipotesi di gas ideale generalmente adottata.Viene definito un gas ideale equivalente avente rapporto dei calori specifici _s funzione del numero di Mach e della quota di volo ed in base alle equazioni di conservazione della massa, della quantità di moto e dell'energia, valide attraverso all'onda d'urto, vengono derivate delle relazioni definenti lo stato termodinamico e dinamico del gas a valle dell'onda d'urto. Tali relazioni costituiscono una estensione delle classiche relazioni dell'urto valide per il gas ideale alle quali si riducono per _s=.La dipendenza del rapporto dei calori specifici _s del gas ideale equivalente, dal numero di Mach e dalla quota è stata stabilita sulla base delle soluzioni ottenute da Huber per il gas reale.

---

     

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.     

Classical limits for the coefficient of variation for the normal distribution

The exact classical limits for the coefficient of variation c for the normal distribution are derived. The hand-calculating approximated classical limits for c having high accuracy are given to meet practical engineering needs. Using Odeh and Owen's computational method and Brent's algorithm, the tables for the r-upper exact classical limits of coefficient of variation for normal distribution are calculated for the different confidence coefficient , the sample size n=1 (1)30,40,60,120, the sample coefficient of variation =0.01(0.01)0.20. It is shown that if n8, 0.20, then the -upper exact classical limits c_u for c are slightly higher than the exact fiducial limits c_u,F for c. if n>8, 0.20, then c_u–c_u,F<5×10^–6.     

Nonlinear viscoelastic properties and change in entanglement structure of linear polymer

Simultaneous measurements of stress relaxation and differential dynamic modulus were made at 268 K over a time scale of 10 to 10⁴⁵ s for nearly monodisperse polybutadiene (M _w =2.2x10⁵, 1,2-structure 70%, M _e =3600) and also one having coarse cross-linking (M _c =29000). Static shear strain ranged from 0.1 to 2.0. In a long-time region (t> _k ), the relaxation modulus G (; t) could be expressed by the product G (0; t) h (y). The observed h() agreed well with the Doi-Edwards theory without use of IA approximation. Both the cured and uncured samples showed initial drop of the differential storage modulus G (), ; t) followed by gradual recovery, but did not attain the value before shearing G (, ; t) for the uncured sample showed smaller values than that for the cured one in the whole measured time scale at the higher strain, confirming the two origins of nonlinear viscoelasticity of well entangled polymer; induced chain anisotropy and induced decrement in entanglement density. G (, ; t) curves for the cured sample agreed well with the BKZ predictions. But the curves for the uncured sample agreed well with the BKZ prediction only at the time scale of t< _k . BKZ prediction showed significant upward deviations at t> _k . Such the differences are discussed in terms of the two origins.Dedicated to Prof. John D. Ferry on the occasion of his 85th birthday.     

Instability of a plane axisymmetric flow with a shear discontinuity in the velocity profile

It is proposed to investigate the stability of a plane axisymmetric flow with an angular velocity profile (r) such that the angular velocity is constant when r < r_O – L and r > r_O + L but varies monotonically from ₁ to ₂ near the point r_O, the thickness of the transition zone being small L r_O, whereas the change in velocity is not small ¦₂ – ₁¦ ₂, ₁. Obviously, as L O short-wave disturbances with respect to the azimuthal coordinate (k=m/r_O 1/r_O) will be unstable with a growth rate-close to the Kelvin—Helmholtz growth rate. In the case L=O (i.e., for a profile with a shear-discontinuity) we find the instability growth rate _O and show that where the thickness of the discontinuity L is finite (but small) the growth rate does not differ from _O up to terms proportional to kL 1 and 1/m 1. Using this example it is possible to investigate the effect of rotation on the flow stability. It is important to note that stabilization (or destabilization) of the flow in question by rotation occurs only for three-dimensional or axisymmetric perturbations.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 111–114, January–February, 1985.     

Asymptotic solution to the problem of hypersonic flow over bodies in the neighborhood of the separation point of a thin shock layer

A study is made of the problem of hypersonic flow of an inviscid perfect gas over a convex body with continuously varying curvature. The solution is sought in the framework of the asymptotic theory of a strongly compressed gas [1–4] in the limit M when the specific heat ratio tends to 1. Under these assumptions, the disturbed flow is situated in a thin shock layer between the body and the shock wave. At the point where the pressure found by the Newton-Buseman formula vanishes there is separation of the flow and formation of a free layer next to the shock wave [1–4]. The singularity of the asymptotic expansions with respect to the parameter ₁ = ( –1)/( + 1) associated with separation of the strongly compressed layer has been investigated previously by various methods [3–9]. Local solutions to the problem valid in the neighborhood of the singularity have been obtained for some simple bodies [3–7]. Other solutions [7, 9] eliminate the singularity but do not give the transition solution entirely. In the present paper, an asymptotic solution describing the transition from the attached to the free layer is constructed for a fairly large class of flows.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 99–105, January–February, 1982.