On elastic-plastic structures with associated stress-strain relations allowing for work softening

Summary For the elementary (finite or infinitesimal) constituents of the structure are assumed generalized stress-strain relations which satisfy the condition of normality but which may exhibit work-softening, concavity of yield surfaces, variation of elastic coefficients with stress and/or plastic strain.The usual phenomena of geometric instability are excluded. Sufficient conditions are formulated for overall stability in spite of the presence of unstable elements, and for uniqueness of the incremental boundary-values problem. Conditions are discussed with a view to applications and expressed in terms of positive definiteness of appropriate quadratic forms.Finally, yield surfaces and flow laws for the structure are examined, and among other things their necessary association is shown.

---

Sommario Per i costituenti elementari (finiti o infinitesimi) della struttura si assumono legami incrementali tra sforzi e deformazioni generalizzati che soddisfano alla condizione di normalità ma che per il resto sono generici, cioè tali da presentare eventualmente incrudimento negativo, concavità del campo elastico, variazione dei coefficienti elastici con gli sforzi e/o con le deformazioni plastiche.Esclusi per il sistema i fenomeni usuali di instabilità geometrica si formulano condizioni sufficienti per la stabilità del complesso nonostante la presenza di parti a funzionamento instabile, e per l'unicità del problema incrementale al contorno. Le condizioni sono discusse in vista delle applicazioni ed espresse in termini di definizione positiva di forme quadratiche opportune.Si esaminano infine le superfici di snervamento e le leggi di scorrimento per l'intera struttura e, tra l'altro, se ne dimostra la necessaria associazione.

---

---

First published in Italian in Rendiconti dell'Istituto Lombardo Classe Scienze e Lettere, A 100, 1966.The present investigation has been promoted and financed by the Consiglio Nazionale delle Ricerche (C.N.R.) at the Istituto di Scienza delle Costruzioni of the Facoltà di Architettura, Politecnico of Milano. Part of it was carried out at Brown University (Providence R. I. - U.S.A.). The author acknowledges with tanks the NATO Fellowship assigned by the C.N.R. in 1964, the encouragement and kind hospitality he received from Professor D.C. Drucker, and the interesting discussions he had with Professor Drucker and Dr. A. C. Palmer.     

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     

Die Temperaturabhängigkeit der Oberflächenspannung von reinen Kältemitteln vom Tripelpunkt bis zum kritischen Punkt

Zusammenfassung Die Oberflächenspannung von sechs reinen Substanzen — SF₆, CCl₃F, CCl₂F₂, CClF₃, CBrF₃ und CHClF₂ — wurde mit Hilfe einer modifizierten Kapillarmethode gemessen. Die zur Berechnung der Oberflächenspannung erforderlichen Sättigungsdichten und wurden teils aus vorhandenen Zustandsgleichungen, teils aus ebenfalls gemessenen Brechungsindizes bestimmt. Die Temperaturabhängigkeit der Oberflächenspannung läßt sich durch einen erweiterten Ansatz nach van der Waals =_O (T_c-T)(1+...) darstellen, wobei bei einfachen Stoffen ein eingliedriger, bei polaren und assoziierenden Stoffen ein zweigliedriger Ansatz notwendig und ausreichend ist. Für den kritischen Exponenten der Oberflächenspannung wurde ein von der molekularen Substanz weitgehend unabhängiger Wert von =1.284±0.005 gefunden.

---

Temperature dependence of surface tension of pure refrigerants from triple point up to the critical point

The surface tension of six fluids (SF₆, CCl₃F, CCl₂F₂, CClF₃, CBrF₃, CHClF₂) have been measured by means of a modified capillary rise method. The liquid vapor densities, which are needed to calculate the surface tension, have partly been determined by means of refractive indices simultaneously measured in the same apparatus. The temperature dependence of the surface tension is described by an extended van der Waals power law =_O(T_c-T)(1+...). For simple fluids one term and for polar and associating fluids two terms are necessary and sufficient. The critical exponent is found to be 1.284 ± 0.005 and nearly independent of the molecular structure.

---

Formelzeichen a² Laplace-Koeffizient - a Parameter - B_O, B_on Koeffizient der Koexistenzkurve - g Erdbeschleunigung - H Höhe, kapillare Steighöhe - LL Lorentz-Lorenz-Funktion oder Refraktionskonstante - M molare Masse - M Zahl der Meßwerte - N Zahl der unbekannten Parameter - n Brechungsindex - p Druck - R,r Radius - s Entropie - SD Standardabweichung - T, t Temperatur - u innere Energie Griechische Formelzeichen Exponent des Laplace-Koeffizienten - Exponent der Koexistenzkurve - 2. Exponent der Oberflächenspannung - Wellenlänge des Lichts - Exponent der Oberflächenspannung - _D Dipolmoment - , Dichte der Flüssigkeit bzw. des Dampfes - Oberflächenspannung - reduzierte Temperatur (1-T/T_c) - ² gewichtete Varianz Indizes c kritischer Zustand - D Differenz - m Mittelwert - T Isotherme - t Zustand am Tripelpunkt - S Zustand am Schmelzpunkt - bezogen auf Oberfläche     

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.     

Application of the method of integral diffusion to investigation of turbulent transfer

A number of authors have critically examined semiempirical mixing length theories [1]. A defect of these theories is connected with the fact that the magnitude of the mixing length, which is assumed to be small in constructing the theory, turns out in experiments to be comparable with the characteristic dimensions of the flow region. Thus, the concept of volume convection [2–4] or integral diffusion [5], which is understood to be a transfer mechanism in which the friction stress is not expressed in terms of the velocity gradient, is introduced along with the concept of gradient diffusion. In addition, there are a number of experimental papers [6] in which it is shown that the turbulent friction stress cannot be equal to zero at the place in the flow where the derivative of the velocity is equal to zero. Mixing length theory does not describe this effect.It is possible to generalize mixing length theory [7–9] in a way which eliminates these defects. Flow of an incompressible fluid is considered.     

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     

Some problems concerning nonstationary currents in dense plasma systems confined by walls

Nonstationary currents are examined in a dense magnetized plasma with 1, in which energy release and heat loss by thermal conduction and radiation are possible. Solutions are found in two limiting cases: ¦f¦ ¦ div (T)¦ and ¦f¦ ¦ div(T)¦ (f is the radiation intensity, is the coefficient of heat conduction, and T is the temperature). In the first case a solution was obtained of some problems of the cooling and heating of a plasma illustrated in part by the evolution in time of the temperature profile in the boundary layer. In the second case an isomorphic solution was found for an arbitrary dependence of the coefficient of heat conduction on the temperature, pressure, and magnetic field.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 1, pp. 3–8, January–February, 1972.The author is grateful to G. I. Budker for formulating the problem.     

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ß.     

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     

Second order gravitational field theories with hilbert gauge

Summary We consider, in the field-theoretical approach, a class of gravitational theories deducible by a variational principle in the unrenormalized pseudo-Euclidean space-time. At first order in the coupling constant f we require the theories to coincide with the Einstein one. Moreover we assume the Hilbert gauge which assure the exclusion of the vector component of the gravitational potential . To get the higher order consistency we substitute the most general energy-momentum tensor for the particle tensorT ^(p) in the field equations. Requiring the latter to be deducible by a variational principle varying the potentials , we get a Lagrangian which, varying the particle coordinates, gives the equations of motion. So we get a class of theories depending on 5 arbitrary parameters. To have observable quantities we have to renormalize. So we realize that, to satisfy the equivalence principle, we have to put one of the arbitrary parameters equal to zero. With this choice the class of theories coincides at second order with general relativity.

---

Sommario Si vuole ottenere una classe di teorie gravitazionali deducibili da un principio variazionale, nell'ambito della teoria dei campi e nello spazio-tempo pseudoeuclideo non-rinormalizzato. Si richiede che tali teorie coincidano, al primo ordine nella costante di accoppiamento f, con la teoria di Einstein. Si assume inoltre la gauge di Hilbert al fine di escludere la presenza della componente vettoriale del potenziale . Per ottenere la consistenza al secondo ordine delle equazioni di campo, si sostituisce, in queste ultime, al tensore della particellaT ^(p) il più generale tensore energia-quantità-di-moto . Imponendo alle equazioni di campo di essere deducibili mediante un principio variazionale ove si varino i potenziali , si ottiene una lagrangiana che, ove si varino le coordinate della particella di prova, dà le equazioni di moto. In tal modo si ottiene una classe di teorie dipendenti da 5 parametri arbitrari. Per un confronto con i dati sperimentali è necessario rinormalizzare, onde esprimere quantità osservabili. Si dimostra così che per soddisfare il principio di equivalenza al secondo ordine è necessario porre uno dei 5 parametri uguale a zero e che, con tale scelta, l'intera classe di teorie coincide, al secondo ordine, con la relatività generale.

---

---

Research sponsored by the CNR, Gruppi di ricerca Matematica     

Effect of Porous Coatings on Stability of Hypersonic Boundary Layers

The influence of ultrasoundabsorbing coatings on stability of hypersonic boundary layers is considered. Two types of coatings were used in experiments: feltmetal with a random porous microstructure and a sheet perforated by blind cylindrical microchannels. The experiments were performed in a wind tunnel at a Mach number M = 5.95 on sharp cones with a 7° apex halfangle. Evolution of natural disturbances and artificially induced wave packets in the boundary layer was studied with the help of hotwire anemometry. Spatial characteristics of artificial disturbances were obtained. It is demonstrated that such coatings exert a stabilizing effect on secondmode disturbances.     

On an inversion theorem for conical regions in elasticity theory

Summary Two-dimensional stress singularities in wedges have already drawn attention since a long time. An inverse square-root stress singularity (in a 360° wedge) plays an important role in fracture mechanics.Recently some similar three-dimensional singularities in conical regions have been investigated, from which one may be also important in fracture mechanics.Spherical coordinates are r, , . The conical region occupied by the elastic homogeneous body (and possible anisotropic) has its vertex at r=0. The mantle of the cone is described by an arbitrary function f(, )=0. The displacement components be u. For special values of (eigenvalues) there exist states of displacements (eigenstates) % MathType!MTEF!2!1!+-% feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqef0uAJj3BZ9Mz0bYu% H52CGmvzYLMzaerbd9wDYLwzYbItLDharqqr1ngBPrgifHhDYfgasa% acOqpw0xe9v8qqaqFD0xXdHaVhbbf9v8qqaqFr0xc9pk0xbba9q8Wq% Ffea0-yr0RYxir-Jbba9q8aq0-yq-He9q8qqQ8frFve9Fve9Ff0dme% GabaqaaiGacaGaamqadaabaeaafiaakabbaaa6daaahjxzL5gapeqa% aiaadwhadaWgaaWcbaGaeqOVdGhabeaakiabg2da9iaadkhadaahaa% WcbeqaaiabeU7aSbaakiaadAgadaWgaaWcbaGaeqOVdGhabeaakiaa% cIcacqaH7oaBcaGGSaGaeqiUdeNaaiilaiabfA6agjaacMcaaaa!582B!\[u_\xi = r^\lambda f_\xi (\lambda ,\theta ,\Phi )\],which may satisfy rather arbitrary homogeneous boundary conditions along the generators.The paper brings a theorem which expresses that if is an eigenvalue, then also-1- is an eigenvalue. Though the theorem is related to a known theorem in Potential Theory (Kelvin's theorem), the proof has to be given along quite another line.

---

Zusammenfassung Zwei-dimensionale Spannungssingularitäten in keilförmigen Gebieten sind schon längere Zeit untersucht worden und neuerdings auch ähnliche drei-dimensionale Singularitäten in konischen Gebieten.Kugelkoordinaten sind r, , . Das konische Gebiet hat seine Spitze in r=0. Der Mantel des Kegels lässt sich beschreiben mittels einer willkürlichen Funktion f(, )=0. Die Verschiebungskomponenten seien u. Für spezielle Werte von (Eigenwerte) bestehen Verschiebunszustände % MathType!MTEF!2!1!+-% feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqef0uAJj3BZ9Mz0bYu% H52CGmvzYLMzaerbd9wDYLwzYbItLDharqqr1ngBPrgifHhDYfgasa% acOqpw0xe9v8qqaqFD0xXdHaVhbbf9v8qqaqFr0xc9pk0xbba9q8Wq% Ffea0-yr0RYxir-Jbba9q8aq0-yq-He9q8qqQ8frFve9Fve9Ff0dme% GabaqaaiGacaGaamqadaabaeaafiaakabbaaa6daaahjxzL5gapeqa% aiaadwhadaWgaaWcbaGaeqOVdGhabeaakiabg2da9iaadkhadaahaa% WcbeqaaiabeU7aSbaakiaadAgadaWgaaWcbaGaeqOVdGhabeaakiaa% cIcacqaH7oaBcaGGSaGaeqiUdeNaaiilaiabfA6agjaacMcaaaa!582B!\[u_\xi = r^\lambda f_\xi (\lambda ,\theta ,\Phi )\],welche homogene Randwerte der Beschreibenden des Kegels entlang genügen.Das Bericht bringt ein Theorem, welches aussagt, das und =–1– beide Eigenwerte sind.

     

Some self-similar solutions of the Leibenzon equation

Self-similar one-dimensional solutions of the Leibenzon equation c²_t= _zz ^k (z 0, k 2) are considered. Approximate solutions are constructed for the two cases in which the initial value = ₁ = const > 0 and on the boundary either a constant value = ₂ < ₁ is maintained or the flow (directed outwards) is given. In the first problem the dependence of the boundary flow on the governing parameters is determined. A characteristic property of the types of motion in question is the existence near the boundary of a region, expanding with time, in which the flow is almost independent of the coordinate.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 145–150, September–October, 1991.     

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.     

Heat Transfer on the Surface on a Spherically Blunted Cone Exposed to a Supersonic Spatial Flow with Injection of a Coolant Gas

The heattransfer processes in a supersonic spatial flow around a spherically blunted cone with allowance for heat overflow along the longitudinal and circumferential coordinates and injection of a coolant gas are studied numerically. The prospects of using highly heatconducting materials and injection of a coolant gas for reduction of the maximum temperatures at the body surface are demonstrated. The solutions of the direct and inverse problems in one, two, and threedimensional formulations for different shell materials are compared. The error of the thinwall method in determining the heat flux on the heatloaded boundary of the body is estimated.     

Exact Solutions of the Hydrodynamic Equations Derived from Partially Invariant Solutions

The paper proposes a heuristic approach to constructing exact solutions of the hydrodynamic equations based on the specificity of these equations. A number of systems of hydrodynamic equations possess the following structure: they contain a reduced system of n equations and an additional equation for an extra function w. In this case, the reduced system, in which w = 0, admits a Lie group G. Taking a certain partially invariant solution of the reduced system with respect to this group as a seed:rdquo; solution, we can find a solution of the entire system, in which the functional dependence of the invariant part of the seed solution on the invariants of the group G has the previous form. Implementation of the algorithm proposed is exemplified by constructing new exact solutions of the equations of rotationally symmetric motion of an ideal incompressible liquid and the equations of concentrational convection in a plane boundary layer and thermal convection in a rotating layer of a viscous liquid.     

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).     

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.     

A non-dualistic unified field theory of gravitation,electromagnetism and spin

The wisdom of classicalunified field theories in the conceptual framework of Weyl, Eddington, Einstein and Schrödinger has often been doubted and in particular there does not appear to be any empirical reason why the Einstein-Maxwell (E-M) theory needs to be geometrized. The crux of the matter is, however not whether the E-M theory is aesthetically satisfactory but whether it answers all the modern questions within the classical context. In particular, the E-M theory does not provide a classical platform from which the Dirac equation can be derived in the way Schrödinger's equation is derived from classical mechanics via the energy equation and the Correspondence Principle. The present paper presents a non-dualistic unified field theory (UFT) in the said conceptual framework as propounded by M. A. Tonnelat. By allowing the metric formds ²=g dx ^v x ^v and the non-degenerate two-formF=(1/2> _l) dx ^vdx ^vto enter symmetrically into the theory we obtain a UFT which contains Einstein's General Relativity and the Born-Infeld electrodynamics as special cases. Above all, it is shown that the Dirac equation describing the electron in an external gravito-electromagnetic field can be derived from the non-dualistic Einstein equation by a simple factorization if the Correspondence Principle is assumed.     

Approximate solutions for drag coefficient of bubbles moving in shear-thinning elastic fluids

The drag coefficient for bubbles with mobile or immobile interface rising in shear-thinning elastic fluids described by an Ellis or a Carreau model is discussed. Approximate solutions based on linearization of the equations of motion are presented for the highly elastic region of flow. These solutions are in reasonably good agreement with the theoretical predictions based on variational principles and with published experimental data. C _D Drag coefficient - E ^* Differential operator [E ^{* 2} = ²/² + (sin/ ²)/(1/sin /)] - El Ellis number - F _D Drag force - K Consistency index in the power-law model for non-Newtonian fluid - n Flow behaviour index in the Carreau and power-law models - P Dimensionless pressure [=(p – p ₀)/₀ (U /R)] - p Pressure - R Bubble radius - Re ₀ Reynolds number [= 2R U /₀] - Re Reynolds number defined for the power-law fluid [= (2R) ⁿ U ^2–n /K] - r Spherical coordinate - t Time - U Terminal velocity of a bubble - u Velocity - Wi Weissenberg number - Ellis model parameter - Rate of deformation - Apparent viscosity - ₀ Zero shear rate viscosity - Infinite shear rate viscosity - Spherical coordinate - Parameter in the Carreau model - ^* Dimensionless time [=/(U /R)] - Dimensionless length [=r/R] - Second invariant of rate of deformation tensors - ^* Dimensionless second invariant of rate of deformation tensors [=/(U /R)²] - Second invariant of stress tensors - ^* Dimensionless second invariant of second invariant of stress tensor [= / ₀ ² (U /R)²] - Fluid density - Shear stress - ^* Dimensionless shear stress [=/ ₀ (U /R)] - _1/2 Ellis model parameter - ₁ ^2/* Dimensionless Ellis model parameter [= _1/2/ ₀(U /R)] - Stream function - ^* Dimensionless stream function [=/U R ²]