The rheology of polymeric liquid crystals

The molecular theory of Doi has been used as a framework to characterize the rheological behavior of polymeric liquid crystals at the low deformation rates for which it was derived, and an appropriate extension for high deformation rates is presented. The essential physics behind the Doi formulation has, however, been retained in its entirety. The resulting four-parameter equation enables prediction of the shearing behavior at low and high deformation rates, of the stress in extensional flows, of the isotropic-anisotropic phase transition and of the molecular orientation. Extensional data over nearly three decades of elongation rate (10^–2–10¹) and shearing data over six decades of shear rate (10^–2–10⁴) have been correlated using this analysis. Experimental data are presented for both homogeneous and inhomogeneous shearing stress fields. For the latter, a 20-fold range of capillary tube diameters has been employed and no effects of system geometry or the inhomogeneity of the flow-field are observed. Such an independence of the rheological properties from these effects does not occur for low molecular weight liquid crystals and this is, perhaps, the first time this has been reported for polymeric lyotropic liquid crystals; the physical basis for this major difference is discussed briefly. A Semi-empirical constant in eq. (18), N/m² - c rod concentration, rods/m³ - c ^* critical rod concentration at which the isotropic phase becomes unstable, rods/m³ - C interaction potential in the Doi theory defined in eq. (3) - d rod diameter, m - D semi-empirical constant in eq. (19), s^–1 - D _r lumped rotational diffusivity defined in eq. (4), s^–1 - rotational diffusivity of rods in a concentrated (liquid crystalline) system, s^–1 - D _ro rotational diffusivity of a dilute solution of rods, s^–1 - f distribution function defining rod orientation - F tensorial term in the Doi theory defined in eq. (7) (or eq. (19)), s^–1 - G tensorial term in the Doi theory defined in eq. (8) - K _B Boltzmann constant, 1.38 × 10^–23 J/K-molecule - L rod length, m - S scalar order parameter - S tensor order parameter defined in eq. (5) - t time, s - T absolute temperature, K - u unit vector describing the orientation of an individual rod - rate of change ofu due to macroscopic flow, s^–1 - v fluid velocity vector, m/s - v velocity gradient tensor defined in eq. (9), s^–1 - V mean field (aligning) potential defined in eq. (2) - x coordinate direction, m - Kronecker delta (= 0 if = 1 if = ) - _r ratio of viscosity of suspension to that of the solvent at the same shear stress - _s solvent viscosity, Pa · s - ^* viscosity at the critical concentrationc ^*, Pa · s - v ₁, v₂ numerical factors in eqs. (3) and (4), respectively - deviatoric stress tensor, N/m² - volume fraction of rods - ⁰ constant in eq. (16) - ^* volume fraction of rods at the critical concentrationc ^* - average over the distribution functionf(u, t) (= d ²u f(u, t)) - gradient operator - d ²u integral over the surface of the sphere (|u| = 1)     

On natural convection from a vertical plate with a prescribed surface heat flux in porous media

This paper presents a theoretical and numerical investigation of the natural convection boundary-layer along a vertical surface, which is embedded in a porous medium, when the surface heat flux varies as (1 +x ²)), where is a constant andx is the distance along the surface. It is shown that for > -1/2 the solution develops from a similarity solution which is valid for small values ofx to one which is valid for large values ofx. However, when -1/2 no similarity solutions exist for large values ofx and it is found that there are two cases to consider, namely < -1/2 and = -1/2. The wall temperature and the velocity at large distances along the plate are determined for a range of values of .Notation g Gravitational acceleration - k Thermal conductivity of the saturated porous medium - K Permeability of the porous medium - l Typical streamwise length - q _w Uniform heat flux on the wall - Ra Rayleigh number, =gK(q _w /k)l/(v) - T Temperature - Too Temperature far from the plate - u, v Components of seepage velocity in the x and y directions - x, y Cartesian coordinates - Thermal diffusivity of the fluid saturated porous medium - The coefficient of thermal expansion - An undetermined constant - Porosity of the porous medium - Similarity variable, =y(1+x ) ^/3/x ^1/3 - A preassigned constant - Kinematic viscosity - Nondimensional temperature, =(T – T )Ra^1/3 k/qw - Similarity variable, = =y(log_e x)^1/3/x ^2/3 - Similarity variable, =y/x ^2/3 - Stream function     

Asymptotic Nusselt numbers for Newtonian flow through a parallel-plate channel with recycle

General expressions for evaluating the asymptotic Nusselt number for a Newtonian flow through a parallel-plate channel with recycle at the ends have been derived. Numerical results with the ratio of thicknesses as a parameter for various recycle ratios are obtained. A regression analysis shows that the results can be expressed by log Nur^0.83=0.3589 (log)² -0.2925 (log) + 0.3348 forR 3, 0.1 0.9; logNu=0.5982(log)² +0.3755 × 10^–2 (log) +0.8342 forR 10^–2, 0.1 0.9.

---

Asymptotische Nusselt-Zahlen für die Newtonsche Strömung durch einen Kanal aus parallelen Platten mit Rückführung

Zusammenfassung In dieser Untersuchung wurden allgemeine Ausdrücke hergeleitet um die asymptotische Nusselt-Zahl für eine Newtonsche Strömung durch einen Kanal aus parallelen Platten mit Rückführung an den Enden berechnen zu können. Es wurden numerische Ergebnisse mit den Dicken-Verhältnissen, als Parameter für verschiedene Rückführungs-verhältnisse, erhalten. Eine Regressionsanalyse zeigt, daß die Ergebnisse wie folgt ausgedrückt werden können: log Nur^0,83=0,3589 (log)² -0,2925 (log) + 0,3348 fürR 3, 0,1 0,9; logNu=0,5982(log)² +0,3755 × 10^–2 (log) + 0,8342 fürR 10^–2, 0,1 0,9.

---

Nomenclature A₁ shooting value,d(0)/d - A₂ shooting value,d(1)/d - B channel width - Gz Graetz number, U_bW²/L - h _m logarithmic average convective heat transfer coefficient - h _x average local convective heat transfer coefficient - k thermal conductivity - L channel length - Nu average local Nusselt number, 2 h_xW/k - Nu _m logarithmic average Nusselt number, 2h_mW/k - R recycle ratio, reverse volume flow rate divided by input volume flow rate - T temperature of fluid - T _m bulk temperature, Eq. (8) - T ₀ temperature of feed stream - T _s wall temperature - U velocity distribution - U _b reference velocity,V/BW - V input volume flow rate - v dimensionless velocity distribution, U/U_b - W channel thickness - x longitudinal coordinate - y transversal coordinate - Z₁-z₆ functions defined in Eq. (A1) - thermal diffusivity - least squares error, Eq. (A7) - weight, Eqs. (A8), (A9) - dimensionless coordinate,y/W - dimensionless coordinate,x/GzL - function, Eq. (7)     

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.     

Der Spannungs- und Verschiebungszustand drehsymmetrischer Membrane mit beliebig gekrümmter Meridiankurve

Zusammenfassung Die Einführung von Zylinderkoordinaten (x, r, ) in die Gleichgewichtsbedingungen der Schnittkräfte bzw. in die Beziehungen zwischen Verzerrung und Verschiebungen am differentialen Schalenabschnitt ermöglicht die Berechnung des Spannungs- und Verschiebungszustandes von drehsymmetrischen Membranen mit beliebig gekrümmter Meridiankurve auf die Integration einer einfachen, linearen partiellen Differentialgleichung zweiter Ordnung für eine charakteristische FunktionF bzw. zurückzuführen. Eine geschlossene Lösung und damit eine Darstellung der Schnittkräfte und Verschiebungen durch explizite Formeln ist bei harmonischer Belastung cosn für zwei Funktionsgruppen=x ² und=x ^–3 möglich. Im Sonderfall der drehsymmetrischen und der antimetrischen Belastung mitn=0 undn=1 gelten die Gleichungen der Schnitt- und Verschiebungsgrößen für eine beliebige Meridianfunktion=(). Die Betrachtungen der Randbedingungen offener Schalen bei harmonischer Belastung geben über die infinitesimalen Deformationen einer drehsymmetrischen Membran mit überall negativer Krümmung Aufschluß.     

Solution of the equations for a regular beam emitted by a curvilinear surface in the nonstationary case

An analytical solution is given of the equations for a regular beam (§1) emitted by an arbitrary surface in the nonstationary case and the - and T-limited states for nonzero initial velocity (§ §2-4). It is assumed that the emitter is the coordinate surface x¹=0 in the orthogonal system x¹ (i=1, 2, 3), and the current density J, the electric field , and the magnetic field H are given functions J(t, x²,x³), (t, x², x³), and H (x¹, x², x³). The solution is given in the form of series in terms of (X¹) with coefficients that are functions of t, x², and x³. These coefficients are determined from recurrence relations ( =1/3, 1/2, 1, depending on the emission conditions). Plane, cylindrical, and spherical diodes are considered in § 5 in the case in which the high-frequency component of the current density J is not small in comparison with its constant components.     

Boundary layer effect on thermal NO decomposition behind incident shock waves

The thermal decomposition of nitric oxide (diluted in Argon) has been measured behind incident shock waves by means of IR diode laser absorption spectroscopy. In two independent runs the diode laser was tuned to the=0 =1² _3/2 R(18.5)-rotational vibrational transition and the=1 =2² _3/2 R(20.5)-rotational vibrational transition of nitric oxide, respectively. These two transitions originating from the vibrational ground state (=0) and the first excited vibrational state (=1) were selected in order to probe the homogeneity along the absorption path. The measured NO decomposition could satisfactorily be described by a chemical reaction mechanism after taking into account boundary layer corrections according to the theory of Mirels. The study forms a further proof of Mirels' theory including his prediction of the laminar-turbulent transition. It also shows, that the inhomogeneities from the boundary layer do not affect the IR linear absorption markedly.This article was processed using Springer-Verlag T_EX Shock Waves macro package 1.0 and the AMS fonts, developed by the American Mathematical Society.     

Backward Uniqueness for Parabolic Equations

It is shown that a function u satisfying |_t+u|M(|u|+|u|), |u(x, t)|Me^M|x|2 in (ⁿ \ (B_R) × [0, T] and u(x, 0) = 0 for xⁿ \ B_R must vanish identically in ⁿ \ B_R×[0, T].     

Turbulent heat and mass transfer for Pr ≫ 1 and law of turbulent diffusion damping at a solid wall

The question of determining the law of damping for the turbulent diffusion coefficient at a smooth wall according to data on mass and heat transfer for Pr 1 is discussed. It is proved that the hypothesis that this law is determined by the first member of the Taylor series expansion of , namely, / = yⁿ ₊ is valid in the Pr range from 10³ to 10⁵ only under the assumption that the subsequent terms in the expansion have smaller coefficients. A statistical analysis of electrochemical and other experiments devoted to this problem shows that apparently n = 3, but singularities in the experimental results do not permit making a final conclusion. Requirements on a conclusive experiment are formulated on the basis of the analysis made.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 172–175, March–April, 1977.     

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     

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.     

On forced vibration of anisotropic cylinders

Summary The dynamic response of a circular cylinder with thick walls of transverse curvilinear isotropy subjected to a uniformly distributed pressure varying periodically with time is analyzed by means of the Laplace transformation, and the exact solution is obtained in closed form. The previously obtained solutions for forced vibrations with isotropy, and free vibrations with transverse curvilinear isotropy are included as special cases of the general results reported here.Nomenclature t time - r, , z cylindrical coordinates - _ii components of normal strain - _ii components of normal stress - u radial displacement - c _ij elastic constant - mass density - c ² c ₁₁/ - ² c ₂₂/c ₁₁ - a, b inner, outer radius of the cylinder - , A, B constants - forced angular frequency - function defined by (9) - p, real, complex variables - constant defined by (14) - real number - , Lamé elastic constants - J (x) Bessel function of first kind - Y (x) Bessel function of second kind - I (x) modified Bessel function of first kind - K (x) modified Bessel function of second kind     

On the Jacobian conjecture for global asymptotic stability

An old conjecture says that, for the two-dimensional system of ordinary differential equationsx=f(x), wheref: ^{2 2},f C ¹, andf(0)=0 the originx=0 should beglobally asymptotically stable (i.e., a stable equilibrium and all trajectoriesx(t) converge to it ast +) whenever the following conditions on the Jacobian matrixJ(x) off hold: trJ(x) < 0, detJ(x) > 0, x ² It is known that if such anf is globallyone-to-one as a mapping of the plane into itself, then the origin is a globally asymptotically stable equilibrium point for the systemx =f(x). In this paper we outline a new strategy to tackle the injectivity off, based on anauxiliary boundary value problem. The strategy is shown to be successful if the norm of the matrixJ(x) ^T J(x)t/det J(x) is bounded or, at least, grows slowly (for instance, linearly) as ¦x¦ t.     

Flow in porous media I: A theoretical derivation of Darcy's law

Stokes flow through a rigid porous medium is analyzed in terms of the method of volume averaging. The traditional averaging procedure leads to an equation of motion and a continuity equation expressed in terms of the volume-averaged pressure and velocity. The equation of motion contains integrals involving spatial deviations of the pressure and velocity, the Brinkman correction, and other lower-order terms. The analysis clearly indicates why the Brinkman correction should not be used to accommodate ano slip condition at an interface between a porous medium and a bounding solid surface.The presence of spatial deviations of the pressure and velocity in the volume-averaged equations of motion gives rise to aclosure problem, and representations for the spatial deviations are derived that lead to Darcy's law. The theoretical development is not restricted to either homogeneous or spatially periodic porous media; however, the problem ofabrupt changes in the structure of a porous medium is not considered.Roman Letters 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² - Ae area of entrances and exits for the -phase contained within a unit cell, m² - B second order tensor used to represent the velocity deviation (see Equation (3.30)) - b vector used to represent the pressure deviation (see Equation (3.31)), m^–1 - d distance between two points at which the pressure is measured, m - g gravity vector, m/s² - K Darcy's law permeability tensor, m² - L characteristic length scale for volume averaged quantities, m - characteristic length scale for the -phase (see Figure 2), m - characteristic length scale for the -phase (see Figure 2), m - n unit normal vector pointing from the -phase toward the -phase (n =–n ) - n _e unit normal vector for the entrances and exits of the -phase contained within a unit cell - 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 and radius of a capillary tube, m - 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² - arbitrary function used in the representation of the velocity deviation (see Equations (3.11) and (B1)), m/s - arbitrary function used in the representation of the pressure deviation (see Equations (3.12) and (B2)), s^–1     

Impacts of local dispersion and first-order decay on solute transport in randomly heterogeneous porous media

Stochastic subsurface transport theories either disregard local dispersion or take it to be constant. We offer an alternative Eulerian-Lagrangian formalism to account for both local dispersion and first-order mass removal (due to radioactive decay or biodegradation). It rests on a decomposition of the velocityv into a field-scale componentv , which is defined on the scale of measurement support, and a zero mean sub-field-scale componentv _s , which fluctuates randomly on scales smaller than. Without loss of generality, we work formally with unconditional statistics ofv _s and conditional statistics ofv . We then require that, within this (or other selected) working framework,v _s andv be mutually uncorrelated. This holds whenever the correlation scale ofv is large in comparison to that ofv _s . The formalism leads to an integro-differential equation for the conditional mean total concentration c which includes two dispersion terms, one field-scale and one sub-field-scale. It also leads to explicit expressions for conditional second moments of concentration cc. We solve the former, and evaluate the latter, for mildly fluctuatingv by means of an analytical-numerical method developed earlier by Zhang and Neuman. We present results in two-dimensional flow fields of unconditional (prior) mean uniformv . These show that the relative effect of local dispersion on first and second moments of concentration dies out locally as the corresponding dispersion tensor tends to zero. The effect also diminishes with time and source size. Our results thus do not support claims in the literature that local dispersion must always be accounted for, no matter how small it is. First-order decay reduces dispersion. This effect increases with time. However, these concentration moments c and cc of total concentrationc, which are associated with the scale below, cannot be used to estimate the field-scale concentrationc directly. To do so, a spatial average over the field measurement scale is needed. Nevertheless, our numerical results show that differences between the ensemble moments ofc and those ofc are negligible, especially for nonpoint sources, because the ensemble moments ofc are already smooth enough.     

The wedge subjected to tractions: a paradox resolved

The classical two-dimensional solution provided by Lévy for the stress distribution in an elastic wedge, loaded by a uniform pressure on one face, becomes infinite when the opening angle 2 of the wedge satisfies the equation tan 2 = 2. Such pathological behavior prompted the investigation in this paper of the stresses and displacements that are induced by tractions of O(r ^–) as r0. The key point is to choose an Airy stress function which generates stresses capable of accommodating unrestricted loading. Fortunately conditions can be derived which pre-determine the form of the necessary Airy stress function. The results show that inhomogeneous boundary conditions can induce stresses of O(r ^–), O(r ^– ln r), or O(r ^– ln² r) as r0, depending on which conditions are satisfied. The stress function used by Williams is sufficient only if the induced stress and displacement behavior is of the power type. The wedge loaded by uniform antisymmetric shear tractions is shown in this paper to exhibit stresses of O(ln r) as r0 for the half-plane or crack geometry. At the critical opening angle 2, uniform antisymmetric normal and symmetric shear tractions also induce the above type of stress singularity. No anticipating such stresses, Lévy used an insufficiently general Airy stress function that led to the observed pathological behavior at 2.     

Flow in porous media III: Deformable media

Stokes flow in a deformable medium is considered in terms of an isotropic, linearly elastic solid matrix. The analysis is restricted to steady forms of the momentum equations and small deformation of the solid phase. Darcy's law can be used to determine the motion of the fluid phase; however, the determination of the Darcy's law permeability tensor represents part of the closure problem in which the position of the fluid-solid interface must be determined.Roman Letters A interfacial area of the- interface contained within the macroscopic system, m² - A interfacial area of the- interface contained within the averaging volume, 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 a unit cell, m² - A _e ^* area of entrances and exits for the-phase contained within a unit cell, m² - E Young's modulus for the-phase, N/m² - e i unit base vectors (i = 1, 2, 3) - g gravity vector, m²/s - H height of elastic, porous bed, m - k unit base vector (=e ₃) - characteristic length scale for the-phase, m - L characteristic length scale for volume-averaged quantities, m - n unit normal vector pointing from the-phase toward the-phase (n = -n ) - p pressure in the-phase, N/m² - P p – g·r, N/m² - r ₀ radius of the averaging volume, m - r position vector, m - t time, s - T total stress tensor in the-phase, N/m² - T ⁰ hydrostatic stress tensor for the-phase, N/m² - u displacement vector for the-phase, m - V averaging volume, m₃ - V volume of the-phase contained within the averaging volume, m³ - v velocity vector for the-phase, m/s Greek Letters V /V, volume fraction of the-phase - mass density of the-phase, kg/m³ - shear coefficient of viscosity for the-phase, Nt/m² - first Lamé coefficient for the-phase, N/m² - second Lamé coefficient for the-phase, N/m² - bulk coefficient of viscosity for the-phase, Nt/m² - T –T ⁰ , a deviatoric stress tensor for the-phase, N/m²     

Extended Korn's inequalities and the associated best possible constants

We are concerned with the coerciveness of the strain energy E(u) (in linear elasticity) associated with a displacement vector u on the Sobolev space H¹ () or its subspaces, a domain in ⁿ representing an isotropic elastic body—certain specific cases are called Korn's inequalities. Sufficient (and necessary) conditions on the Lamé moduli for E(·) to be coercive (or uniformly positive) on such spaces are given, and the associated best possible constants are obtained for some cases.     

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.