On the compactness of quasi-conforming element spaces and the convergence of quasi-conforming element method

In this paper, the compactness of quasi-conforming element spaces and the—convergence of quasi-conforming element method are discussed. The well-known Rellich compactness theorem is generalized to the sequences of quasi-conforming element spaces with certain properties, and the generalized Poincare inequality. The generalized Friedrichs inequality and the generalized inequality of Poincare-Friedrichs are proved true for them. The error estimates are also given. It is shown that the quasi-conforming element method is convergent if the quasi-conforming element spaces have the approximability and the strong continuity, and satisfy the rank condition of element and pass the test IPT. As practical examples, 6-parameter, 9-paramenter, 12-paramenter, 15-parameter, 18-parameter and 21-paramenter quasi-conforming elements are shown to be convergent, and their L²²()-errors are O(h), O(h), O(h ² ), O(h ² ), O(h ), and O(h ⁴ ) respectively.     

Viscoelastic properties of semidilute poly(methyl methacrylate) solutions

Viscoelastic properties were examined for semidilute solutions of poly(methyl methacrylate) (PMMA) and polystyrene (PS) in chlorinated biphenyl. The number of entanglement per molecule, N, was evaluated from the plateau modulus, G _N . Two time constants, _s and ₁, respectively, characterizing the glass-to-rubber transition and terminal flow regions, were evaluated from the complex modulus and the relaxation modulus. A time constant _k supposedly characterizing the shrink of an extended chain, was evaluated from the relaxation modulus at finite strains. The ratios ₁/ _s and _k / _s were determined solely by N for each polymer species. The ratio ₁/ _s was proportional to N ^4.5 and N ^3.5 for PMMA and PS solutions, respectively. The ratio _k / _s was approximately proportional to N ^2.0 in accord with the prediction of the tube model theory, for either of the polymers. However, the values for PMMA were about four times as large as those for PS. The result is contrary to the expectation from the tube model theory that the viscoelasticity of a polymeric system, with given molecular weight and concentration, is determined if two material constants _s and G _N are known.     

Effect of suspended particles on the stability of plane poiseuille flow

The behavior of the neutral stability curves is investigated for various values of the particle relaxation time and mass concentration 0 100 and 0 f 0.1. It is shown that as increases from zero the flow is at first destabilized and then at >6 becomes stable, while at >40 the stabilizing effect of the dispersed phase grows weaker. It is found that there is a certain interval 10< <40 on which the flow is most stable.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 46–53, January–February, 1986.     

The Preston tube in adiabatic compressible flow

A new procedure for the reduction of Preston tube data is introduced, based on the van Driest transformation. It appears to give results agreeing with the better calibration experiments, although a significant assumption in its derivation is violated.List of Symbols M _s Mach number sensed by Preston tube - M Friction Mach number (=u/wall sound speed) - R Gas constant - T _w Wall temperature - d Diameter of Preston tube - h Height of effective centre of Preston tube - p Preston tube pressure difference reading - p _i Equivalent incompressible Preston tube reading - p _w Wall pressure - r Recovery factor (=0.896) - u Friction velocity (=[_w/wall density]^1/2) - Empirical constant allowing for departure from Crocco temperature-velocity correlation (=0.975) - Specific heat ratio - Fluid kinematic viscosity - _w Wall shear-stress     

Testing viscometric predictions of the internal viscosity model,using dilute viscous theta solutions

A stress-symmetrized internal viscosity (I.V.) model for flexible polymer chains, proposed by Bazua and Williams, is scrutinized for its theoretical predictions of complex viscosity ^* () = – i and non-Newtonian viscosity (), where is frequency and is shear stress. Parameters varied are the number of submolecules,N (i.e., molecular weightM = NM _s ); the hydrodynamic interaction,h ^*; and/f, where andf are the I.V. and friction coefficients of the submolecule. Detailed examination is made of the eigenvalues _p (N, h ^*) and how they can be estimated by various approximations, and property predictions are made for these approximations.Comparisons are made with data from our preceding companion paper, representing intrinsic properties [], [], [] in very viscous theta solutions, so that theoretical foundations of the model are fulfilled. It is found that [ ()] data can be predicted well, but that [ ()] data cannot be matched at high. The latter deficiency is attributed in part to unrealistic predictions of coil deformation in shear.     

Die swell of filled polymer melts

The Barus effect in polypropylene and polystyrene blended with a variety of fillers at various concentrations was investigated using a capillary extrusion rheometer. If the die swell is defined as the square of the ratio of the extrudate diameterd to the die diameterD, it is found to depend on the apparent shear stress _W . Below a certain value of _w the relation =B _B ^A applies. The die swell, _M , of a filled polymer depends on the type, size and volume fraction of the filler. In particular,A increases as the volume fraction increases and is largest for powders, smaller for flakes and smallest for fibres, whereasB shows the opposite trend but to a lesser extent.     

Stability of a radially bounded rotating fluid heated from below

The stability of a bounded rotating cylinder of fluid heated from below is treated mathematically under the assumptions of stationary onset and axisymmetry. Critical Rayleigh numbers are computed by Galerkin's method as a function of the Taylor number and cylinder aspect ratio for Taylor numbers,10⁶. The constraining effect of the side walls is shown to decrease with either increasing or increasing radius/height ratios. For>10⁶, most cylinders, excluding extremely tall ones, will appear infinite in horizontal extent as far as stability characteristics are concerned. The form of the motion at onset is discussed in relation to previous work.     

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     

Turbulent boundary layer over a smooth wall with widely separated transverse square cavities

Laser-Doppler velocimetry (LDV) measurements and flow visualizations are used to study a turbulent boundary layer over a smooth wall with transverse square cavities at two values of the momentum thickness Reynolds number (R =400 and 1300). The cavities are spaced 20 element widths apart in the streamwise direction. Flow visualizations reveal a significant communication between the cavities and the overlying shear layer, with frequent inflows and ejections of fluid to and from cavities. There is evidence to suggest that quasi-streamwise near-wall vortices are responsible for the ejections of fluid out of the cavities. The wall shear stress, which is measured accurately, increases sharply immediately downstream of the cavity. This increase is followed by a sudden decrease and a slower return to the smooth wall value. Integration of the wall shear stress in the streamwise direction indicates that there is an increase in drag of 3.4% at bothR .Nomenclature C _f skin friction coefficient - C _fsw friction coefficient for a continuous smooth wall - k height of the cavity - k ⁺ ku / - R Reynolds number based on momentum thickness (U ₁ /v) - Rx Reynolds number based on streamwise distance (U ₁ x/) - s streamwise distance between two cavities - t time - t ⁺ tu ₂ / - U ₁ freestream velocity - mean velocity inx direction - u,v,w rms turbulent intensities inx,y andz directions - u local friction velocity - u _sw friction velocity for a continuous smooth wall - w width of the cavity - x streamwise co-ordinate measured from the downstream edge of the cavity - y co-ordinate normal to the wall - z spanwise co-ordinate - y ⁺ yu / - boundary layer thickness - ₀ boundary layer thickness near the upstream edge of the cavity - _i thickness of internal layer - kinematic viscosity of water - ⁺ zu / - momentum thickness     

An analytical solution for the velocity and shear rate distribution of non-ideal Bingham fluids in concentric cylinder viscometers

An analytical solution is presented for the calculation of the flow field in a concentric cylinder viscometer of non-ideal Bingham-fluids, described by the Worrall-Tuliani rheological model. The obtained shear rate distribution is a function of the a priori unknown rheological parameters. It is shown that by applying an iterative procedure experimental data can be processed in order to obtain the proper shear rate correction and the four rheological parameters of the Worrall-Tuliani model as well as the yield surface radius. A comparison with Krieger's correction method is made. Rheometrical data for dense cohesive sediment suspensions have been reviewed in the light of this new method. For these suspensions velocity profiles over the gap are computed and the shear layer thicknesses were found to be comparable to visual observations. It can be concluded that at low rotation speeds the actually sheared layer is too narrow to fullfill the gap width requirement for granular suspensions and slip appears to be unavoidable, even when the material is sheared within itself. The only way to obtain meaningfull measurements in a concentric cylinder viscometer at low shear rates seems to be by increasing the radii of the viscometer. Some dimensioning criteria are presented.Notation A, B Integration constants - C Dimensionless rotation speed = µ/_y - c = 2µ - d = ₀ ²–2c_y - f() = (–₀)²+2c(–_y) - r Radius - r _b Bob radius - r _c Cup radius - r _y Yield radius - r ₀ Stationary surface radius - r Rotating Stationary radius - Y ₀ Shear rate parameter = /µ Greek letters Shear rate - = (r _y /r _b )²– 1 - µ Bingham viscosity - µ₀ Initial differential viscosity - µ µ₀-µ - Rotation speed - Angular velocity - Shear stress - _b Bob shear stress - _B Bingham stress - _y (True) yield stress - ₀ Stress parameter = _B +µ Y ₀ - _B - _y      

On properties of solutions of functional differential equations continuously differentiable on (0, +∞)

We establish the saddle-point property of the system of functional differential equations (t) = Ax(t) + Bx((t)) + C ((t)) + f (x(t), x((t))), (0) = 0.Translated from Neliniini Kolyvannya, Vol. 7, No. 3, pp. 302–310, July–September, 2004.     

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.     

A general model for the effective viscosity of pseudoplastic and dilatant fluids

Summary A new and very general expression is proposed for correlation of data for the effective viscosity of pseudoplastic and dilatant fluids as a function of the shear stress. Most of the models which have been proposed previously are shown to be special cases of this expression. A straightforward procedure is outlined for evaluation of the arbitrary constants.

---

Zusammenfassung Eine neue und sehr allgemeine Formel wird für die Korrelation der Werte der effektiven Viskosität von strukturviskosen und dilatanten Flüssigkeiten in Abhängigkeit von der Schubspannung vorgeschlagen. Die meisten schon früher vorgeschlagenen Methoden werden hier als Spezialfälle dieser Gleichung gezeigt. Ein einfaches Verfahren für die Auswertung der willkürlichen Konstanten wird beschrieben.

---

Nomenclature b arbitrary constant inSisko model (eq. [5]) - n arbitrary exponent in eq. [1] - x independent variable - y(x) dependent variable - y ₀(x) limiting behavior of dependent variable asx 0 - y(x) limiting behavior of dependent variable asx - z original dependent variable - arbitrary constant inSisko model (eq. [5]) andBird-Sisko model (eq. [6]) - arbitrary exponent in eqs. [2] and [8] - effective viscosity = shear stress/rate of shear - _A effective viscosity at = _A - _B empirical constant in eqs. [2] and [8] - ₀ limiting value of effective viscosity as 0 - ₀() limiting behavior of effective viscosity as 0 - limiting value of effective viscosity as - () limiting behavior of effective viscosity as - rate of shear - arbitrary constant inBird-Sisko model (eq.[6]) - shear stress - _A arbitrary constant in eqs. [2] and [8] - ₀ shear stress at inBingham model - _1/2 shear stress at = ( ₀ + )/2 With 8 figures     

Correlations between porous medium flow data and turbulent tube flow results for aqueous hydroxypropylguar solutions (Part 2)

Turbulent tube flow and the flow through a porous medium of aqueous hydroxypropylguar (HPG) solutions in concentrations from 100 wppm to 5000 wppm is investigated. Taking the rheological flow curves into account reveals that the effectiveness in turbulent tube flow and the efficiency for the flow through a porous medium both start at the same onset wall shear stress of 1.3 Pa. The similarity of the curves = ( _w ) and = ( _w ), respectively, leads to a simple linear relation / =k, where the constantk or proportionality depends uponc. This offers the possibility to deduce (for turbulent tube flow) from (for flow through a porous medium). In conjunction with rheological data, will reveal whether, and if yes to what extent, drag reduction will take place (even at high concentrations).The relation of our treatment to the model-based Deborah number concept is shown and a scale-up formula for the onset in turbulent tube flow is deduced as well.     

Bifurcation Phenomena from Standing Pulse Solutions in Some Reaction-Diffusion Systems

Bifurcation phenomena from standing pulse solutions of the problem is considered. (>0) is a sufficiently small parameter and is a positive one. It is shown that there exist two types of destabilization of standing pulse solutions when decreases. One is the appearance of travelling pulse solutions via the static bifurcation and the other is that of in-phase breathers via the Hopf bifurcation. Furthermore which type of destabilization occurs first with decreasing is discussed for the piecewise linear nonlinearities f and g.     

Displacement type unified equation for bending, buckling and vibration of conical shells and its applications

By using Donnell's simplication and starting from the displacement type equations of conical shells, and introducing a displacement functionU(s,,) (In the limit case, it will be reduced to cylindrical shell displacement function introduced by V. S. Vlasov) and a generalized loadq,(s,,),the equations of conical shells are changed into an eighth—order solvable partial differential equation about the displacement functionU(s,,). As a special case, the general bending problem of conical shells on Winkler foundation has been studied. Detailed numerical results and boundary coefficients for edge unit loads are obtained.The project supported by the National Natural Science Foundation of China.     

Quantitative analysis of flow visualizations in an unsteady channel flow

Quantitative results concerning the modulation of the ejection and bursting frequency in an unsteady channel flow obtained by flow visualizations are presented and compared with probe measurements. The frequency of the imposed velocity oscillations f covers a large range going from the quasi steady limit to the time mean bursting frequency in the corresponding steady flow. The imposed amplitudes of the velocity oscillations are 13% and 20% of the centerline velocity. The bursting process is identified by the intermittent lift up of the dye injected at the wall. Qualitative analysis of the flow visualizations show that the ejection activity at a given phase of the oscillation cycle is repetitive from one cycle to the other. The modulation amplitude of the ejection frequency f _e is sensitive to the imposed frequency. At low imposed frequency f _e is modulated as the wall shear stress, but the inner scaling does not hold when f ⁺ is high. Here, (+) corresponds to the quantities normalized with the inner variables, i.e. the friction velocity u and the viscosity . The grouping of the ejections into bursts show the coexistence of two categories of events which react differently to the forcing. The groups of ejections (Multiple Ejection Bursts) are governed by the modulation of the wall shear stress in the whole imposed frequency range. The solitary ejections (or the Single Ejection Bursts) have modulation amplitudes and phases which differ significantly from those of in the intermediate and high imposed frequency range. There is a good agreement between the flow visualization data and the probe measurements.     

Effect of surface layers on the constriction resistance of an isothermal spot. Part II: Analytical results for thick layers

Solution of a non-homogeneous Fredholm integral equation of the second kind [1], which forms the basis for the evaluation of the constriction resistance of an isothermal circular spot on a half-space covered with a surface layer of different material, is considered for the case when the ratio, , of layer thickness to spot radius is larger than unity. The kernel of the integral equation is expanded into an infinite series in ascending odd-powers of (1/) and an approximate kernel accurate to (^–(2M+1)) is derived therefrom by terminating the series after an arbitrary but finite number of terms, M. The approximate kernel is rearranged into a degenerate form and the integral equation with this approximate kernel is reduced to a system of M linear equations. An explicit analytical solution is obtained for a four-term approximation of the kernel and the resulting constriction resistance is shown to be accurate to (^–9). Solutions of lower orders of accuracy with respect to (1/) are deduced from the four-term solution. The analytical approximations are compared with very accurate numerical solutions and it is shown that the (^–9)-approximation predicts the constriction resistance exceedingly well for any 1 over a four orders of magnitude variation of layer-to-substrate conductivity ratio for both conducting and insulating layers. It is further shown that, for all practical purposes, an (^–3)-approximation gives results of adequate accuracy for > 2.     

The eddy viscosity of turbulent equilibrium layers

Summary Similarity laws for the mean flow and scaling laws for the turbulent motion are used in an attempt to obtain a general expression for the eddy viscosity of equilibrium layers. It is found that =0.09 w ² /w*, in which w ² is a Reynolds stress representative for the region of overlap between the law of the wall and the velocity-defect law, while w* is the logarithmic slope of the mean velocity profile in that region. The distinction between w and w* is related to the strong inhomogeneity of the mean rate of strain in the inner layer. The results of the theory agree with experimental evidence obtained from transpired equilibrium layers.     

Filament stretching equation in Lagrangian coordinates

The general integral of the very simple equation ^21/n/() was found to describe the cross sectional area of filaments of isothermal power law fluids while in transient stretching where is time and is the initial location of fluid molecules at time = 0 given as the distance from a reference point fixed in space. Any such stretching transient given as a solution of the above equation is physically realizable subject to the restrictions > 0 and/ < 0.