The rotation of a suspended axisymmetric ellipsoid in a magnetic field

Under the influence of a uniform and parallel magnetic field, a ferromagnetic fiber suspended in a Newtonian fluid rotates to align with the field direction. This study examines the field-induced rotation process for an individual non-Brownian axisymmetric ellipsoid suspended in a stagnant Newtonian fluid. Theoretical predictions are derived by a perturbation analysis for the limiting case where the strength of the applied magnetic field far exceeds the saturation magnetization of the ellipsoid. Numerical calculations are performed for the more general problem of an ellipsoid with known isotropic, non-hysteretic magnetic properties, using nickel and a stainless steel as examples. The analysis encompasses materials with field-induced, nonlinear magnetic properties, distinguishing these results from the simpler cases where the particle magnetization is either independent of, or linearly dependent on, the strength of the applied external field. In this study, predictions indicate that when the ellipsoid is magnetically saturated, the particle rotation is governed by the magnetoviscous time constant, _MV = _s/₀ M _s ² . It is found that the rotation rate depends strongly on the aspect ratio,a/b, of the ellipsoid, but only weakly on the dimensionless magnetization,M _s/H ₀. A geometric parameter for an ellipsoid, defined in eq. (2.5) - a, b major, minor semi-axes of an axisymmetric ellipsoid - D demagnetization tensor for an ellipsoid - D ^M magnetometric demagnetization tensor, the volume-average ofD ^P (r) - D ^P (r) position dependent demagnetization tensor, implicitly defined in eq. (2.12) - D _xx,D _yy,D _zz demagnetization factors, the diagonal elements ofD. Values for ellipsoids are defined in eq. (2.15) - F ^(m) magnetic force exerted on a body in a magnetic field - H _i ;H _i magnetic field inside a ferromagnetic body; magnitude ofH _i - H ₀;H ₀ magnetic field applied by external sources; magnitude ofH ₀ - h _i ;h _ix,h _iy Cartesian components of dimensionless internal magnetic field,h _i =H _i /H ₀ - I moment of inertia tensor - k geometric parameter for hydrodynamic resistance of a body rotating in a Newtonian fluid given in eq. (2.3) - L ^(h);L _z ^(h) hydrodynamic torque exerted on a rotating body; thez-component of the hydrodynamic torque - L ^(m);L _z ^(m) magnetic torque exerted on a magnetic body in a magnetic field, eq. (2.10); thez-component of the magnetic torque - M;M the magnetization, or dipole moment density, of a magnetic material; the magnitude ofM - M _s the saturation value ofM, approached by all ferromagnetic materials asH _i becomes large (figure 3) - m _s the dimensionless saturation magnetization,M _s/H ₀ - r position vector of a point within a ferromagnetic body - s dummy integration variable in eq. (2.5) - t time - U magnetoquasistatic potential energy of a magnetic body in a magnetic field, given in eq. (2.8) - u curve-fitting variable in eq. (4.1);u = logH _i - V volume of a magnetic particle; for an axisymmetric ellipsoid,V = (4/3) ab ² - x, y, z rectangular coordinate axes fixed in the ellipsoid (figure 1) - angle of inclination of the major axis of the ellipsoid with respect toH ₀ - _s viscosity of the Newtonian suspending medium - µ ₀ the magnetic permeability of free space,µ ₀ =4 · 10^–7H/m - _MV the magnetoviscous time constant, a characteristic time for a process involving a competition of viscous and magnetic stresses - the magnetic susceptibility of a magnetic material, = M/H _i - ; _z angular velocity of a rotating body; angular velocity about thez-axis of an ellipsoid, _z=–d/dt     

The shear modulus in the fiber direction of unidirectional fiber composites based on the concept of mesophase

A new and accurate relationship for the shear modulus in the fiber direction of unidirectional fiber composites is derived, based on the Kerner model concepts but not using its approximate relationships. Furthermore, this model is extended by taking into account the mesophase between the fiber and the matrix accomodating smoothly the mechanical properties between neighbouring phases. The introduction of the mesophase results in an improvement of the theoretical predictions, which now approach close to experimental values for the moduli of different fibrous composites.     

Three-dimensional studies on bicomponent extrusion

The present work is concerned with the mathematical modelling and numerical simulation of three-dimensional (3-D) bicomponent extrusion. The objective is to provide an understanding of the flow phenomena involved and to investigate their impact on the free surface shape and interface configuration of the extruded article. A finite element algorithm for the 3-D numerical simulation of bicomponent stratified free surface flows is described. The presence of multiple free surfaces (layer interface and external free surfaces) requires special free surface update schemes. The pressure and viscous stress discontinuity due to viscosity mismatch at the interface between the two stratified components is handled with both a double node (u–v–w–P ₁ –P ₂ –h ₁ –h ₂) formulation and a penalty function (u–v–w–P–h ₁ –h ₂) formulation.The experimentally observed tendency of the less viscous layer to encapsulate the more viscous layer in stratified bicomponent flows of side-by-side configuration is established with the aid of a fully 3-D analysis in agreement with experimental evidence. The direction and degree of encapsulation depend directly on the viscosity ratio of the two melts. For shear thinning melts exhibiting a viscosity crossover point, it is demonstrated that interface curvature reversal may occur if the shearing level is such that the crossover point is exceeded. Extrudate bending and distortion of the bicomponent system because of the viscosity mismatch is shown. For flows in a sheath-core configuration it is shown that the viscosity ratio may have a severe effect on the swelling ratio of the bicomponent system.Modelling of the die section showed that the boundary condition imposed at the fluid/fluid/wall contact point is critical to the accuracy of the overall solution.     

Accurate evaluation of the transverse shear modulus in unidirectional fibrous composites

A new relation for the prediction of the transverse shear modulus in unidirectional fiber composites has been derived. The theoretical results of this relationship are in better agreement with the experiments than those of other relations, existing in the literature. The discrepancies, which are observed among the theoretical predictions and the experimental values, are explained by the consideration of the boundary layers existing between the matrix and the fibers of the composite. A new model, which includes the intermediate phase between the matrix and the fiber, called the mesophase, is considered in order to take into account the above-mentioned layers.     

A boundary element simulation of the unsteady motion of a sphere in a cylindrical tube containing a viscoelastic fluid

A boundary element method is used to simulate the unsteady motion of a sphere falling under gravity along the centreline of a cylindrical tube containing a viscoelastic fluid. The fluid is modelled by the upper-convected Maxwell constitutive equation. Results show that the viscoelasticity of the liquid leads to a damped oscillation in sphere velocity about its terminal value. The maximum sphere velocity, which occurs in the first overshoot, is approximately proportional to the square root of the Weissenberg number when the ratio of the sphere radius to the tube radius is sufficiently small. Particular attention is also paid to the wall effects. It is shown that a closer wall reduces the oscillatory amplitude of the sphere velocity but increases its frequency. The results suggest that the falling-ball technique, which is now widely used for viscosity measurement, might also be used for the determination of a relaxation time for a viscoelastic fluid.     

On-line measurements of the rheological properties of fermentation broth

The paper is concerned with measurements of rheological properties of fermentation broth. An on-line rheometer, Rheohelix-1, based on the application of a helical screw impeller rotating in a draught tube has been constructed. The instrument was used for measurements of the rheological parameters of fermentation broth of Aspergillus niger in a submerged fermentation process. The results of rheological and some standard measurements have been compared and proved the applicability of the instrument.Presented at the Golden Jubilee Conference of the British Society of Rheology and Third European Rheology Conference, Edinburgh, 3–7 September, 1990.     

Generalized Cole-Cole behavior and its rheological relevance

Starting from an analysis of the rheological behavior of the complex modulus predicted by the Cole-Cole formalism, a generalized Cole-Cole ansatz is suggested in order to overcome the related difficulties. The corresponding rheological constitutive equation with fractional derivatives belonging to the generalized Cole-Cole respondance is stated and the characteristic material functions of the linear viscoelasticity theory (like the dynamic modulus and compliance, the relaxation and ratardation functions, the spectra, etc.) are derived. Model predictions of these functions will be compared with experimental results from dynamical measurements and creep data on different polymer systems which show cooperative phenomena (polymeric glasses and gelling systems). One can see that the modified ansatz fits the data very well, in spite of its relative simplicity.     

An anomalous phenomenon occurring in the flow of viscoelastic fluids out of ducts

When viscoelastic fluids flow out from horizontal ducts many cracks and protruding ridges are formed on the jet surface. As far as the authors know, this phenomenon has not yet been reported. The occurrence of the anomalous phenomenon is not affected by inlet flow condition or duct shape.The phenomenon may be divided into three regimes, namely, a stable state, a breakage state, and a multiplication state. In the breakage state one ridge divides into two parts after its growth, and in the multiplication state a new ridge is suddenly formed at a crack point. With increasing shear rate, the flow pattern of the jet changes from the stable state to the breakage state, and then to the multiplication state.Furthermore, a recoil effect is observed. In this effect many air bubbles rush into the duct from the exit when the flow is quickly stopped.     

Rubber-like elasticity by boundary element method: Finite deformation of a circular elastic slice

We briefly review the phenomenological theory of rubber-like elasticity and report a microstructural model that leads us to eventually adopt a particular constitutive equation, which includes the Neo-Hookean and the Mooney materials. A numerical implementation of the Boundary Element method for solving a general two-dimensional or axisymmetric finite deformation problem is described and tested with some simple deformations. The resulting program is used to analyse the finite deformation of a circular elastic slice perfectly bonded to two parallel rigid end plates; the bottom plate is stationary and the top plate is given a constant displacement. The problem has a free surface which must be found as part of the solution. The results indicate that the Boundary Element method can be an efficient tool for stress-strain analyses with rubber-like materials.     

Transient rheology in a new type of couette apparatus. Application to blood

By using our new air-bearing viscometer different types of transient flow can be studied; in the present work this viscometer was specifically applied to non-Newtonian blood suspensions. To observe the influence of both the aggregation and the deformation of red blood cells (RBCs), different concentrations of fibrinogen and dextran were used: the suspended RBCs were simply washed or rigidified with diamide. From the data three rheological parameters were determined, two of which are relative to the behavior at low and at high shear gradients, respectively. Their values were related to the theory of Taylor, and the internal viscosity of RBCs was estimated to be around 3 cp.     

Der elektroviskose Effekt als Folge elektrostatischer Kraft

Zusammenfassung Die Viskosität vieler Flüssigkeiten ändert sich als Folge elektrischer Felder. Die größten elektroviskosen Effekte findet man unter äußerem Feld bei Suspensionen aus dielektrischen Flüssigkeiten und feinen Festkörperteilchen.Die Flüssigkeiten zeigen ausgeprägtes Bingham-Verhalten. Solange die auf die Wand ausgeübte Scherspannung kleiner als die von der Feldstärke abhängige Haftspannung ist, verhalten sich die Materialien wie Festkörper. Man beobachtet, daß Elektroviskosität proportional zur Feldstärke im Quadrat und umgekehrt proportional zur Schergeschwindigkeit ist, und daß sich die Teilchen der Flüssigkeiten in Richtung der Feldstärke orientieren und faserige Strukturen bilden. In dieser Arbeit erklären wir diese Eigenschaften durch ein einfaches Modell. In diesem Modell resultiert die Scherkraftzunahme aus der elektrostatischen Kraft, die sich als Gradient der im Feld abgespeicherten Energie ergibt.The electroviscosity effect following electrostatic force

---

Apparent viscosity of many fluids changes under an electric field. The largest effects are found with dispersions of fine particles in dielectric liquids subjected to an externally applied field. The large effect is apparently due to induced fibration when the particles align themselves as fibers in field direction. The dispersions exhibit Bingham-behavior, i.e., as long as the shear stress is smaller than the field-dependent yield stress, the material behaves as a solid, a behavior which can be made plausible as the electrodes are mechanically linked together by the fibers (Winslow 1949). In this paper arguments are put forward that the induced fibration is the consequence of electrostatic forces acting on the particles and pulling them into regions of higher field strength. The resulting configuration is an equilibrium configuration and it is postulated that even under shear load, the electric field tends to maintain this configuration such that only the particles close to the electrodes are sheared off, which then form a suspension near the electrode surfaces. A phenomenological theory based on the assumption that the fibers are columns of rectangular particles gives quantitative results for the electroviscosity and the yield stress which agree with experimental observations.

     

Measurement of rheological and ultrasonic properties of food and synthetic polymer melts

A slit die viscometer has been used in conjunction with a co-rotating twin screw extruder to study the rheological behaviour of maize grits, potato powder and low density polyethylene, as a function of feed rate, screw speed and temperature. The shear viscosity of both maize and potato decreased with increasing feed rate. Increasing the temperature or screw speed at any given feed rate also reduced the viscosity. The ultrasonic velocity through the material has also been shown to be sensitive to the extruder operating conditions. Overall, the ultrasonic velocity decreased as screw speed and temperature increased and feed rate decreased.     

The unraveling of a polymer chain in a strong extensional flow

By assuming that in a strong extensional flow a polymer molecule in dilute solution is quickly driven into a folded or kinked state in which drag and entropic elastic forces dominate over the Brownian force, we derive kink dynamics equations that describe the unraveling of the molecule in the extensional flow. Solving these equations numerically, we find that although the ends of the chain move, on average, affinely in the flow field until the chain is unfolded to about 1/3 of its fully extended length, large viscous stresses are produced because the solvent must flow around nonextending strands of polymer that lie between neighboring kinks. These predictions are compared with available experimental data and with other theoretical models.     

Ermittlung des dehnrheologischen Verhaltens thermoplastischer Kautschuke

Zusammenfassung Mit Hilfe eines einfachen Dehnrheometers wurden Dehnfließversuche mit thermoplastischen Kautschukblends (Ibuflex-SEH und Santoprene) durchgeführt. Die auf diese Weise erhaltenen Dehnviskositätsdaten wurden mit den von kautschukmodifizierten Thermoplasten (PVC, PE und PS) verglichen. Nach der Anlaufphase zeigt Ibuflex-SEH ein stationäres Fließverhalten, während die Dehnviskosität von Santoprene als Funktion der Zeit monoton ansteigt, ohne einen stationären Zustand zu erreichen. Die Ursachen des Unterschieds im Fließverhalten der untersuchten Polymermaterialien wurden diskutiert.

---

Extensional flow experiments were carried out on a simple extensional rheometer with the thermoplastic elastomer blends: Ibuflex-SEH and Santoprene. The extensional viscosity data obtained were comparable with those of elastomer modified PVC, PE and PS. After the transient stage Ibuflex-SEH shows a steady-state flow behaviour, whereas the extensional viscosity of Santoprene as a function of time increases monotonically without reaching a stationary state. The causes for these differences were discussed.

Erweiterte Fassung eines Vortrags anläßlich der Jahrestagung der Deutschen Rheologischen Gesellschaft vom 13.–15. Mai 1985 in Berlin     

Experimental investigations of the stability limit of the helical flow of pseudoplastic liquids

The stability of the laminar helical flow of pseudoplastic liquids has been investigated with an indirect method consisting in the measurement of the rate of mass transfer at the surface of the inner rotating cylinder. The experiments have been carried out for different values of the geometric parameter = R ₁/R ₂ (the radius ratio) in the range of small values of the Reynolds number,Re < 200. Water solutions of CMC and MC have been used as pseudoplastic liquids obeying the power law model. The results have been correlated with the Taylor and Reynolds numbers defined with the aid of the mean viscosity value. The stability limit of the Couette flow is described by a functional dependence of the modified critical Taylor number (including geometric factor) on the flow indexn. This dependence, general for pseudoplastic liquids obeying the power law model, is close to the previous theoretical predictions and displays destabilizing influence of pseudoplasticity on the rotational motion. Beyond the initial range of the Reynolds numbers values (Re>20), the stability of the helical flow is not affected considerably by the pseudoplastic properties of liquids. In the range of the monotonic stabilization of the helical flow the stability limit is described by a general dependence of the modified Taylor number on the Reynolds number. The dependence is general for pseudoplastic as well as Newtonian liquids.Nomenclature C _i concentration of reaction ions, kmol/m³ - d = R ₂ –R ₁ gap width, m - F _M () Meksyn's geometric factor (Eq. (1)) - F ₀ Faraday constant, C/kmol - i _l density of limit current, A/m³ - k _c mass transfer coefficient, m/s - n flow index - R ₁,R ₂ inner, outer radius of the gap, m - Re = V _m ·2d·/µ _m Reynolds number - Ta _c = _c ·d^3/2·R ₁ ^1/2 ·/µ _m Taylor number - Z _i number of electrons involved in electrochemical reaction - = R ₁/R ₂ radius ratio - µ apparent viscosity (local), Ns/m² - µ _m mean apparent viscosity value (Eq. (3)), Ns/m² - µ _i apparent viscosity value at a surface of the inner cylinder, Ns/m² - density, kg/m³ - _c angular velocity of the inner cylinder (critical value), 1/s     

Distribution of stresses in a pre-cracked bone specimen

Of concern in the paper is an analysis of the stress-distribution in a specimen of bone having an exterior crack. In conformity to experimental findings, the osseous medium is treated as a composite transversely isotropic material. By considering a numerical example, values of the longitudinal stress at different locations of the bone specimen are obtained by computing the analytically derived expressions, using experimental data for the various physical constants of osseous tissues.     

On the role of stress-induced migration on time-dependent terminal velocities of falling spheres

There is experimental evidence to suggest that even under steady-state conditions the velocity of solid spheres or bubbles moving through viscoelastic fluids can become time dependent. One of the possible explanations offered for interpreting this phenomenon has been the generation of a polymer depleted layer in the line of passage of the particles, which disappears due to the counterbalancing effect of molecular diffusion in the long range. We have done some careful experiments and measured these concentrations to show that no such depletion layers are formed. Alternative explanations of the phenomenon have been examined and the importance of the possible effects of microstructures generated through temporary associations has been emphasised.     

Measurements of slip at the wall during flow of high-density polyethylene through a rectangular conduit

A hot-film probe has been used to measure slip of high-density polyethylene flowing through a conduit with a rectangular cross section. A transition from no slip to a stick-slip condition has been observed and associated with irregular extrudate shape. Appreciable extrudate roughness was initiated at the same flow rate as that at which the relationship between Nusselt number and Péclet number for heat transfer from the probe departed from the behavior expected for a no-slip condition at the conduit wall. A ₁ constant defined by eq. (A3) - C dimensionless group used in eq. (7) - C _p heat capacity - D constant in eq. (13) - f u _s/u - f _lin defined by eq. (A6) - G storage modulus - G loss modulus - k thermal conductivity - L length of hot film in thex-direction - L _eff effective length of large probe found from eq. (A3) - Nu _L Nusselt number, defined for a lengthL by eq. (2) - (Nu _L)₀ value ofNu _L atPe = 0 (eq. (A 1)) - Pe Péclet number,uL/ - Pe ₀ Péclet number in slip flow, eq. (6) - Pe ₁ Péclet number in shear flow, eq. (4) - q _L average heat flux over hot film of lengthL - R _i resistances defined by figure 8 - R _AB correlation coefficient defined by eq. (14) for signalsA andB - T temperature - T _s temperature of probe surface - T ambient temperature - T T _s –T - u average velocity - u _s slip velocity - V _b voltage indicated in figure 8 - W probe dimension (figure 6) - x distance in flow direction (figure 1) - y distance perpendicular to flow direction (figure 1) - thermal diffusivity,k/C _p - wall shear rate - _5% thickness of lubricating layer during probe calibration for introduction of an error no greater than 5%; see Appendix I - ^* complex viscosity - density - time - _c critical shear stress, eq. (13) - _w wall shear stress - frequency characterizing extrudate distortion (figures 12 and 13), or frequency of oscillation during rheometric characterization (figures 18–20) - * quantity obtained from normalized Nusselt number, eq. (A1), or complex viscosity ^* - A actual (small) probe (see Appendix I) - M model (large) probe (see Appendix I)     

Evaluating Herschel-Bulkley fluids with the back extrusion (annular pumping) technique

A mathematical model was developed to describe the behavior of Herschel-Bulkley fluids in a back extrusion (annular pumping) device. A technique was also developed to determine the rheological properties (yield stress, flow behavior index, and consistency coefficient) of these fluids. Mathematical terms were expressed in four dimensionless terms, and graphical aids and tables were prepared to facilitate the handling of the expressions.Nomenclature a radius of the plunger, m - dv/dr shear rate, s^–1 - F force applied to the plunger, N - F _b buoyancy force, N - F _cb force corrected for buoyancy, N - F _T recorded force just before the plunger is stopped, N - F _Te recorded force after the plunger is stopped, N - g acceleration due to gravity, m/s² - H(t) momentary height between plunger and container bottom, m - K a/R, dimensionless - L length of annular region, m - L(t) depth of plunger penetration, m - n flow behavior index, dimensionless - p static pressure, Pa - P _L pressure in excess of hydrostatic pressure at the plunger base, Pa - p ₀ pressure at entrance to annulus, Pa - P pressure drop per unit of length, Pa/m - Q total volumetric flow rate through the annulus, m³/s - r radial coordinate, measured from common axis of cylinder forming annulus, m - R radius of outer cylinder of annulus, m - s reciprocal of n, dimensionless - t time, s - T dimensionless shear stress, defined in Eq. (3) - T ₀ dimensionless yield stress, defined in Eq. (4) - T _w dimensionless shear stress at the plunger wall - _p velocity of plunger, m/s - velocity, m/s - mass density of fluid, kg/m³ - Newtonian viscosity, Pa s - P p ₀ –p _L , Pa - consistency coefficient, Pa sⁿ - value of where shear stress is zero - _–, ₊ limits of the plug flow region (Fig. 1) - r/R - shear stress, Pa - _y yield stress, Pa - _w shear stress at the plunger wall, Pa - dimensionless flow rate defined in Eq. (24) - dimensionless velocity defined by Eq. (5) - _–, ₊ dimensionless velocity outside the plug flow region - _max dimensionless maximum velocity in the plug flow region - _p dimensionless velocity at the plunger wall