Viscometric properties of viscous theta solutions of monodisperse polystyrene

Theta solvents for polystyrene are prepared from high-viscosity blends of styrene and low-molecular-weight polystyrene, and then used to make dilute solutions with monodisperse polystyrene solutes of high-M = 2.3, 6.0, 9.0, 18.0 · 10⁵. A Weissenberg rheogoniometer is used to measure the non-Newtonian viscosity as a function of shear stress, for low values, and also the complex viscosity components and as functions of frequency. A capillary viscometer is used for high- measurements of(). Viscometric properties, at room temperature, are analyzed as functions of high-molecular-weight solute concentrationc with parameters of constant or to obtain [()], [ ()], and [ ()]. Such a collection of data has apparently not previously been available for polymers in theta solvents (in which Gaussian chain statistics prevail). Also unique is the achievement of high stress ( = 2 10⁴ Pa) at low shear rate, by virtue of high solvent viscosity which is not characteristic of other known theta solvents.     

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.     

New calculations on the Faraday effect in wave guides

Summary The propagation of electromagnetic waves is investigated theoretically for a round wave guide, containing a gyroelectrie-gyromagnetic medium with gyroaxis parallel to the guide in the form of a cylindrical shell of thickness, adjacent to the wall of the guide. An equation is set up, permitting to compute the change in the propagation constant due to the presence of the shell, including terms proportional to ². Assuming only the presence of gyromagnetism, the change ₁ of first order in for TE-waves is determined and is found to be the same fpr right- and left-circular polarization. The second order difference ₂ ⁺ — ₂ ^- for the two senses of polarization, however, appears to have a non-vanishing value which, just like ₁ can be expressed in terms of the radius of the guide, the frequency, the dielectric constant and the elements of the gyromagnetic permeability tensor which characterize the medium of the shell.     

Experimental study of mutual interference between two spheres placed on a plane boundary

This paper describes an experimental study of the mutual interference between two spheres placed on a plane boundary. The experiment was carried out in an N. P. L. type wind-tunnel having a working section of 500×500×2000 mm³ in size at a Reynolds number of 4.74×10⁴. The surface-pressure distributions of two spheres were measured for the various relative positions of two spheres and the drag, side-force, and lift coefficients were determined from surface-pressure distributions. The separation of the flow and the formation of vortices were observed by the method of visualization. The distributions of velocities, and turbulent intensities of the flow past two spheres were measured. The experimental results for two spheres were compared with those of a single sphere.List of symbols C _D drag coefficient - C _L lift coefficient - C _p surface-pressure coefficient of sphere=(P-P )(qU ² ) - C _s coefficients of side force - D diameter of sphere [mm] - P static pressure [Pa] - P static pressure in free stream [Pa] - Re Reynolds number= DU/v - S spacing between the centers of two adjoining spheres in plane view [mm] - U time-mean velocity in X-direction [m/s] - [m/s] free stream velocity [m/s] - u, v, w X, Y and Z-components of velocity fluctuation [m/s] - X, Y, Z coordinate axes with origin at the bottom center of test sphere, X, Y, Z axis being taken in the streamwise, lateral and vertical directions respectively [mm] (Fig. 1) - latitude angle [°] - longitude angle [°] - angle between the line connected with the centers of two spheres and wind direction [°] (Fig. 2) - kinematic viscosity of air [m²/s] - density of air [N/m³] This paper was presented at the 10th Symposium on Turbulence, University of Missouri-Rolla, Sept. 22–24, 1986     

Flow in porous media II: The governing equations for immiscible,two-phase flow

The Stokes flow of two immiscible fluids through a rigid porous medium is analyzed using the method of volume averaging. The volume-averaged momentum equations, in terms of averaged quantities and spatial deviations, are identical in form to that obtained for single phase flow; however, the solution of the closure problem gives rise to additional terms not found in the traditional treatment of two-phase flow. Qualitative arguments suggest that the nontraditional terms may be important when / is of order one, and order of magnitude analysis indicates that they may be significant in terms of the motion of a fluid at very low volume fractions. The theory contains features that could give rise to hysteresis effects, but in the present form it is restricted to static contact line phenomena.Roman Letters (, = , , and ) A interfacial area of the- interface contained within the macroscopic system, m² - A _e area of entrances and exits for the -phase contained within the macroscopic system, m² - A interfacial area of the- interface contained within the averaging volume, m² - A ^* interfacial area of the- interface contained within a unit cell, m² - A _e ^* area of entrances and exits for the-phase contained within a unit cell, m² - g gravity vector, m²/s - H mean curvature of the- interface, m^–1 - H area average of the mean curvature, m^–1 - H – H , deviation of the mean curvature, m^–1 - I unit tensor - K Darcy's law permeability tensor, m² - K permeability tensor for the-phase, m² - K viscous drag tensor for the-phase equation of motion - K viscous drag tensor for the-phase equation of motion - L characteristic length scale for volume averaged quantities, m - characteristic length scale for the-phase, m - n unit normal vector pointing from the-phase toward the-phase (n = –n ) - p _c p – P , capillary pressure, N/m² - p pressure in the-phase, N/m² - p intrinsic phase average pressure for the-phase, N/m² - p – p , spatial deviation of the pressure in the-phase, N/m² - r ₀ radius of the averaging volume, m - t time, s - v velocity vector for the-phase, m/s - v phase average velocity vector for the-phase, m/s - v intrinsic phase average velocity vector for the-phase, m/s - v – v , spatial deviation of the velocity vector for the-phase, m/s - V averaging volume, m³ - V volume of the-phase contained within the averaging volume, m³ Greek Letters V /V, volume fraction of the-phase - mass density of the-phase, kg/m³ - viscosity of the-phase, Nt/m² - surface tension of the- interface, N/m - viscous stress tensor for the-phase, N/m² - / kinematic viscosity, m²/s     

On the boundary conditions at the macroscopic level

We study the problem of the boundary conditions specified at the boundary of a porous domain in order to solve the macroscopic transfer equations obtained by means of the volume-averaging method. The analysis is limited to the case of conductive transport but the method can be extended to other cases. A numerical study enables us to illustrate the theoretical results in the case of a model porous medium. Roman Letters _sf interfacial area of the s-f interface contained within the macroscopic system m² - A _sf interfacial area of the s-f interface contained within the averaging volume m² - C _p mass fraction weighted heat capacity, kcal/kg/K - d _s , d _f microscopic characteristic length m - g vector that maps to _s, m - h vector that maps to _f , m - K _eff effective thermal conductivity tensor, kcal/m s K - l REV characteristic length, m - L macroscopic characteristic length, m - n _fs outwardly directed unit normal vector for the f-phase at the f-s interface - n _e outwardly directed unit normal vector at the dividing surface - T ^* macroscopic temperature field obtained by solving the macroscopic equation (3), K - V averaging volume, m³ - V _s , V _f volume of the considered phase within the averaging volume, m³ - volume of the macroscopic system, m³ - _s , _f volume of the considered phase within the volume of the macroscopic system, m³ - dividing surface, m² Greek Letters _s , _f volume fraction - ratio of thermal conductivities - _s , _f thermal conductivities, kcal/m s K - spatial average density, kg/m³ - microscopic temperature, K - ^* microscopic temperature corresponding to T ^* , K - spatial deviation temperature K - error on the temperature due to the macroscopic boundary conditions, K - spatial average - ^s , ^f intrinsic phase average     

Einige Beziehungen zwischen Transportkoeffizienten und thermodynamischen Zustandsgrö\en

Zusammenfassung Eine früher mitgeteilte Beziehung [1] zwischen der ViskositÄt und dem isenthalpen Joule-Thomson-Koeffizienten _h wird für kleine Drücke theoretisch begründet und an sieben Stoffen nachgeprüft. Die WärmeleitfÄhigkeit wird als Funktion von c_v für drei Stoffe dargestellt.

---

Some relations between transport coefficients and thermodynamical properties

A formerly given relation [1] between viscosity and isenthalpic Joule-Thomson-coefficient _h is proved theoretically for small pressures and checked with seven substances. The heat conductivity is presented as a function ofc_v for three substances.

---

Bezeichnungen B ^* dimensionsloser zweiter Virialkoeffizient - B ₁ ^* Ableitung vonB ^* nachT ^*.B ₁ ^*=T ^* dB ^*/dT ^* - c _v isochore spezifische WärmekapazitÄt - C _p ^o isobare molare WärmekapazitÄt fürp 0 - h spezifische Enthalpie - k Boltzmann-Konstante.k=R/N _A - M molare Masse - N _A Avogadro-Konstante - p Druck - R molare Gaskonstante - R _i spezifische Gaskonstante des Stoffesi - it Celsius-Temperatur - T Kelvin-Temperatur - T ^* dimensionslose Temperatur.T ^*=kT/ - _h isenthalper Joule-Thomson-Koeffizient. _h=(T/p)_h - , Konstanten der Potentialfunktion - ViskositÄt - WärmeleitfÄhigkeit - (^2,2)^* dimensionsloses Sto\integral     

Schwere symmetrische Kreisel kleiner Nutationen

Zusammenfassung Zur Integration der Eulerschen Bewegungsgleichungen schwerer symmetrischer Kreisel werden der Winkel (t) (Abb. 1) durch (t)=₀+(t) ersetzt und in sämtlichen Reihenentwicklungen von abhängiger Funktionen die Potenzen höheren als zweiten Grades vernachlässigt. Dadurch ist es möglich, die Eulerschen Winkel (t), (t) und (t) durch elementare Formeln zu beschreiben und somit sind die wesentlichsten Erscheinungen im Bewegungsablauf der schweren symmetrischen Kreisel einfach zu übersehen.     

Remarks on the bifurcation of solutions of a nonlinear eigenvalue problem

Conclusions We have investigated solutions of equation (3) when ² is an eigenvalue of the linearized operator (13) and when it is not. In Section 4 we have shown that for 0 and ² = _i ² we have exactly two nontrivial solutions which bifurcate to the right of _i ² ; these solutions are shown to exist in an interval ( _i ² , _i ² + ₀). The method of Section 3 may then be used to extend these two solutions to the right of _i ² + ₀ providing that ²= _i ² + ₀ is not an eigenvalue of the linear operator (13) evaluated at = ± ₁. Either a solution can be uniquely extended, or there exists a value of ²where the bifurcation method must be applied again³.While the method used here gives the exact number of solutions bifurcating from _i ² , other problems remain open; for example, it is still not proven that the two bifurcating branches have i zeros, as is the case for Hammerstein operators with oscillation kernels [4]. The conjecture of Odeh and Tadjbakhsh that there are exactly 2(i+1) nontrivial solutions in the interval _i ² < _i ^+1/2 remains un-answered, although it would be proven if one could show that there is no secondary bifurcation as in the cases of Kolodner [7] and Coffman [8].     

Existence theorem for a minimum problem with free discontinuity set

We study the variational problem Where is an open set in ⁿ ,n2gL ^q () L (), 1q<+, O<, <+ andH _n–1 is the (n–1)-dimensional Hausdorff Measure.     

Messung des binären Diffusionskoeffizienten in einem Entmischungssystem mit Hilfe der Photonen-Korrelationsspektroskopie

Zusammenfassung Die dynamische Lichtstreuung in Form der Photonen-Korrelationsspektroskopie wird in einem Entmischungssystem beispielhaft zur Messung des binären Diffusionskoeffizienten eingesetzt. Mit einem Versuchsaufbau wird nahe dem kritischen Entmischungspunkt mit der Homodyn-Methode und weit davon entfernt mit der Heterodyn-Methode gearbeitet. Ein Verfahren ermöglicht die Korrektur der Störeinflüsse des Heterodyn-Anteils bei der Homodyn-Messung. Fürn-Hexan/Nitrobenzol wird der Diffusionskoeffizient für vier unterschiedliche Konzentrationswerte als Funktion der Temperatur ausgemessen. Bei der kritischen Konzentration zeigt er bei gleicher Temperatur den kleinsten Wert und läßt sich nahe dem kritischen Entmischungspunkt durch einen einfachen Potenzansatz mit der reduzierten Temperaturdifferenz =T (– T_c)/T_c ausdrücken. Die statistische Genauigkeit ist besser als 1 %. Die Übereinstimmung mit Literaturwerten ist gut.

---

Measuring of the binary diffusion coefficient in a separation system with photon-correlation-spectroscopy

The dynamic light scattering in form of photon-correlation spectroscopy is examplary used in a separation system for measuring the binary diffusion coefficient. In a test setup the homodyntechnique is used near the critical separation point and in distance the heterodyn-technique is used. A special method allows the correction of the disturbing influences of the heterodyn-part using the homodyn-measuring. Forn-hexane/nitrobenzene the diffusion coefficient is measured for four different concentration values as a function of temperature. At the critical concentration with constant temperature the coefficient shows the minimum value and it is expressed near the critical separation point with an elementary exponential equation with the reduced temperature difference =T (– T_c)/T_c. The statistical precision is better than 1%. The conformity with the literature is well.

---

Abbreviations

Formelzeichen A Konstante - B Konstante - b Konstante - b ₁ Konstante - b₂ Konstante - C Konstante - c Konzentration - c _c kritische Konzentration - D Konstante - D ₁₂ binärer Diffusionskoeffizient - E ₀ elektrisches Feld des Laserlichts - E _s elektrisches Feld des Streulichts - G () Korrelationsfunktion - I ₀ Intensität des Referenzlichts - _s mittlere Intensität des Streulichts - k ₀ Wellenvektor des Laserlichts - k _s Wellenvektor des Streulichts - n Brechungsindex der zu untersuchenden Flüssigkeit - p Druck - q Streuvektor - R Ortsvektor - r Ortsvektor - T Temperatur, Zeit - t Zeit - x Molenbruch Griechische Buchstaben reduzierte Temperaturdifferenz; Dielektrizität - Frequenz des Laserlichts - statistische Schwankungen der Dielektrizitätskonstante - ₀ Wellenlänge des Laserlichts - Streuwinkel - Zeit - _c Zeitkonstante - kritischer Exponent     

On Perturbative Expansions to the Stochastic Flow Problem

When analyzing stochastic steady flow, the hydraulic conductivity naturally appears logarithmically. Often the log conductivity is represented as the sum of an average plus a stochastic fluctuation. To make the problem tractable, the log conductivity fluctuation, f, about the mean log conductivity, lnK _G, is assumed to have finite variance, _f ². Historically, perturbation schemes have involved the assumption that _f ²<1. Here it is shown that _f may not be the most judicious choice of perturbation parameters for steady flow. Instead, we posit that the variance of the gradient of the conductivity fluctuation, _f ², is more appropriate hoice. By solving the problem withthis parameter and studying the solution, this conjecture can be refined and an even more appropriate perturbation parameter, , defined. Since the processes f and f can often be considered independent, further assumptions on f are necessary. In particular, when the two point correlation function for the conductivity is assumed to be exponential or Gaussian, it is possible to estimate the magnitude of f in terms of _f and various length scales. The ratio of the integral scale in the main direction of flow ( _x ) to the total domain length (L^*), _x ²=_x/L^*, plays an important role in the convergence of the perturbation scheme. For _x smaller than a critical value _c, _x < _c, the scheme's perturbation parameter is =_f/_x for one- dimensional flow, and =_f/_x ² for two-dimensional flow with mean flow in the x direction. For _x > _c, the parameter =_f/_x ³ may be thought as the perturbation parameter for two-dimensional flow. The shape of the log conductivity fluctuation two point correlation function, and boundary conditions influence the convergence of the perturbation scheme.     

A noninvasive method for the measurement of flow-induced surface displacement of a compliant surface

A noninvasive optical method is described which allows the measurement of the vertical component of the instantaneous displacement of a surface at one or more points. The method has been used to study the motion of a passive compliant layer responding to the random forcing of a fully developed turbulent boundary layer. However, in principle, the measurement technique described here can be used equally well with any surface capable of scattering light and to which optical access can be gained. The technique relies on the use of electro-optic position-sensitive detectors; this type of transducer produces changes in current which are linearly proportional to the displacement of a spot of light imaged onto the active area of the detector. The system can resolve displacements as small as 2 m for a point 1.8 mm in diameter; the final output signal of the system is found to be linear for displacements up to 200 m, and the overall frequency response is from DC to greater than 1 kHz. As an example of the use of the system, results detailing measurements obtained at both one and two points simultaneously are presented.List of symbols C _t elastic transverse wave speed = (G/)^1/2 - d ⁺ spot diameter normalized by viscous length scale - G frequency average of G() - G() shear storage modulus - G() shear loss modulus - l. viscous length scale = v/u _* - N total number of sampled data values - r separation vector for 2-point measurements = (, ) - rms root-mean-square value - R momentum thickness Reynolds number = U _t8/v - t time - u (y) mean streamwise component of velocity in boundary layer - u _* friction velocity = (t _w/)^1/2 - U free-stream velocity - x, y, z longitudinal, normal and spanwise directions - y _o undisturbed surface position - vertical component of compliant surface displacement - ₉₉ boundary layer thickness for which u(y) = 0.99 U _t8 - _l viscous sublayer thickness 5 l _* - frequency average of G()/ - boundary layer momentum thicknes = - fluid dynamic viscosity - v fluid kinematic viscosity = / - , longitudinal, spanwise components of separation vector r - fluid density - time delay - _w wall shear stress     

Peristaltic pumping of a second-order fluid in a planar channel

The peristaltic motion of a non-Newtonian fluid represented by the constitutive equation for a second-order fluid was studied for the case of a planar channel with harmonically undulating extensible walls. A perturbation series for the parameter ( half-width of channel/wave length) obtained explicit terms of 0(²), 0(²Re²) and 0(₁Re²) respectively representing curvature, inertia and the non-Newtonian character of the fluid. Numerical computations were performed and compared to the perturbation analysis in order to determine the range of validity of the terms.Presented at the second conference Recent Developments in Structured Continua, May 23–25, 1990, in Sherbrooke, Québec, Canada     

Chaos and integrability in a nonlinear wave equation

We consider the parametrized family of equations _tt ,u- _xx u-au+u ₂ ² u=O,x(0,L), with Dirichlet boundary conditions. This equation has finite-dimensional invariant manifolds of solutions. Studying the reduced equation to a four-dimensional manifold, we prove the existence of transversal homoclinic orbits to periodic solutions and of invariant sets with chaotic dynamics, provided that =2, 3, 4,.... For =1 we prove the existence of infinitely many first integrals pairwise in involution.     

Production of temperature fluctuations in grid turbulence: Wiskind's experiment revisited

Measurements have been made in nearly-isotropic grid turbulence on which is superimposed a linearly-varying transverse temperature distribution. The mean-square temperature fluctuations, , increase indefinitely with streamwise distance, in accordance with theoretical predictions, and consistent with an excess of production over dissipation some 50% greater than values recorded in previous experiments. This high level of production has the effect of reducing the ratio,r, of the time scales of the fluctuating velocity and temperature fields. The results have been used to estimate the coefficient,C, in Monin's return-to-isotropy model for the slow part of the pressure terms in the temperature-flux equations. An empirical expression by Shih and Lumley is consistent with the results of earlier experiments in whichr 1.5, C 3.0, but not with the present data where r 0.5, C 1.6. Monin's model is improved when it incorporates both time scales.List of symbols C coefficient in Monin model, Eq. (5) - M grid mesh length - m exponent in power law for temperature variance, x ^m - n turbulence-energy decay exponent,q ² x ^-n - p production rate of - p pressure - q ² - R microscale Reynolds number - r time-scale ratiot/t - T mean temperature - U mean velocity - mean-square velocity fluctuations (turbulent energy components) - turbulent temperature flux - x, y, z spatial coordinates - temperature gradient dT/dy - thermal diffusivity - dissipation rate ofq ²/2 - dissipation rate of - Taylor microscale (₂=5q₂/) - temperature microscale - _v temperature-flux correlation coefficient, /v - dimensionless distance from the grid,x/M     

Inertial Manifolds on Squeezed Domains

Let be an arbitrary smooth bounded domain in and > 0 be arbitrary. Squeeze by the factor in the y-direction to obtain the squeezed domain = {(x,y)(x,y)}. In this paper we study the family of reaction-diffusion equations

Laser velocimetry measurements in a horizontal gas-solid pipe flow

This paper presents laser measurements of particle velocities in a horizontal turbulent two-phase pipe flow. A phase Doppler particle analyzer, (PDPA), was used to obtain particle size, velocity, and rms values of velocity fluctuations. The particulate phase consisted of glass spheres 50 m in diameter with the volume fraction of the suspension in the range _p=10^-4 to _p=10^-3. The results show that the turbulence increases with particle loading.List of symbols a particle diameter - C _va velocity diameter cross-correlation - d pipe diameter - Fr ² Froude number - g gravitational constant - p(a) Probability density of the particle diameter - Re pipe Reynolds number based on the friction velocity - T characteristic time scale of the energy containing eddies - T _L integral scale of the turbulence sampled along the particle path - u, U, u characteristic fluid velocities: fluctuating, mean and friction - v characteristic velocity of the paricle fluctuations - f expected value of any random variable f - f¦g expected value of f given a value of the random variable g - _p particle volume fraction - _p particle response time - absolute fluid viscosity - v kinematic fluid viscosity - _p, _f densities, particle and fluid - _a ² particle diameter variance - _va ² velocity variance due to the particle diameter variance - _vT ² total particle velocity variance - _vt ² particle velocity variance due to the response to the turbulent field     

Image processing of hydrogen bubble flow visualization for determination of turbulence statistics and bursting characteristics

A technique is described which employs automated image processing of hydrogen-bubble flow visualization pictures to establish local, instantaneous velocity profile information. Hydrogen bubble flow visualization sequences are recorded using a high-speed video system and then digitized, stored, and evaluated by a VAX 11/780 computer. Employing special smoothing and gradient detection algorithms, individual bubble-lines are computer identified, which allows local velocity profiles to be constructed using time-of-flight techniques. It is demonstrated how this techniques may be used to 1) determine local velocity behavior as a function of position and time, 2) evaluate time-averaged turbulence properties, and 3) correlate probe-type turbulent burst detection techniques with the corresponding visualization data.List of symbols Re Reynolds number based on momentum thickness, u / - t ⁺ nondimensional time tu ² / - T VITA variance averaging time period - u shear velocity = - u local instantaneous streamwise velocity,x-direction - u local fluctuating streamwise velocity,x-direction - u ⁺ nondimensional streamwise velocity, /u - local normal velocity,y-direction - w local spanwise velocity,z-direction - x ⁺ nondimensional coordinate in streamwise direction xu /v - y ⁺ nondimensional coordinate normal to wall, yu /v Greek momentum thickness, - kinematic viscosity - _w wall shear stress This paper was presented at the Ninth Symposium on Turbulence, University of Missouri-Rolla, October 1–3, 1984     

Asymptotic analysis of equilibrium states for rotating turbulent flows

The equilibrium states of homogeneous turbulence simultaneously subjected to a mean velocity gradient and a rotation are examined by using asymptotic analysis. The present work is concerned with the asymptotic behavior of quantities such as the turbulent kinetic energy and its dissipation rate associated with the fixed point (/kS)=0, whereS is the shear rate. The classical form of the model transport equation for (Hanjalic and Launder, 1972) is used. The present analysis shows that, asymptotically, the turbulent kinetic energy (a) undergoes a power-law decay with time for (P/)<1, (b) is independent of time for (P/)=1, (c) undergoes a power-law growth with time for 1<(P/)<(C ₂–1), and (d) is represented by an exponential law versus time for (P/)=(C ₂–1)/(C ₁–1) and (/kS)>0 whereP is the production rate. For the commonly used second-order models the equilibrium solutions forP/,II, andIII (whereII andIII are respectively the second and third invariants of the anisotropy tensor) depend on the rotation number when (P/kS)=(/kS)=0. The variation of (P/kS) andII versusR given by the second-order model of Yakhot and Orzag are compared with results of Rapid Distortion Theory corrected for decay (Townsend, 1970).