Reservoir transport equations by compositional approach

The objective of this paper is to present an overview of the fundamental equations governing transport phenomena in compressible reservoirs. A general mathematical model is presented for important thermo-mechanical processes operative in a reservoir. Such a formulation includes equations governing multiphase fluid (gas-water-hydrocarbon) flow, energy transport, and reservoir skeleton deformation. The model allows phase changes due to gas solubility. Furthermore, Terzaghi's concept of effective stress and stress-strain relations are incorporated into the general model. The functional relations among various model parameters which cause the nonlinearity of the system of equations are explained within the context of reservoir engineering principles. Simplified equations and appropriate boundary conditions have also been presented for various cases. It has been demonstrated that various well-known equations such as Jacob, Terzaghi, Buckley-Leverett, Richards, solute transport, black-oil, and Biot equations are simplifications of the compositional model.Notation List B reservoir thickness - B formation volume factor of phase - Cⁱ mass of component i dissolved per total volume of solution - C ⁱ mass fraction of component i in phase - C heat capacity of phase at constant volume - C_p heat capacity of phase at constant pressure - D ⁱ hydrodynamic dispersion coefficient of component i in phase - D_MTf thermal liquid diffusivity for fluid f - F = F(x, y, z, t) defines the boundary surface - f_p fractional flow of phase - g gravitational acceleration - H_p enthalpy per unit mass of phase - J_p volumetric flux of phase - k_rf relative permeability to fluid f - k₀ absolute permeability of the medium - M_p ⁱ mass of component i in phase - n porosity - N rate of accretion - P_f pressure in fluid f - p_ca capillary pressure between phases and =p-p - Rⁱ rate of mass transfer of component i from phase to phase - Rⁱ _source source rate of component i within phase - S saturation of phase - s gas solubility - T temperature - t time - U displacement vector - u velocity in the x-direction - v velocity in the y-direction - V volume of phase - V_s velocity of soil solids - W_i body force in coordinate direction i - x horizontal coordinate - z vertical coordinate Greek Letters _p volumetric coefficient of compressibility - _T volumetric coefficient of thermal expansion - _ij Kronecker delta - volumetric strain - _m thermal conductivity of the whole matrix - internal energy per unit mass of phase - gf suction head - density of phase - _ij tensor of total stresses - _ij tensor of effective stresses - volumetric content of phase - f viscosity of fluid f     

The dynamic and steady shear properties of synovial fluid and of the components making up synovial fluid

Steady-shear and dynamic properties of a pooled sample of cattle synovial fluid have been measured using techniques developed for low viscosity fluids. The rheological properties of synovial fluid were found to exhibit typical viscoelastic behaviour and can be described by the Carreau type A rheological model. Typical model parameters for the fluid are given; these may be useful for the analysis of the complex flow problems of joint lubrication.The two major constituents, hyaluronic acid and proteins, have been successfully separated from the pooled sample of synovial fluid. The rheological properties of the hyaluronic acid and the recombined hyaluronic acid-protein solutions of both equal and half the concentration of the constituents found in the original synovial fluid have been measured. These properties, when compared to those of the original synovial fluid, show an undeniable contribution of proteins to the flow behaviour of synovial fluid in joints. The effect of protein was found to be more prominent in hyaluronic acid of half the normal concentration found in synovial fluid, thus providing a possible explanation for the differences in flow behaviour observed between synovial fluid from certain diseased joints compared to normal joint fluid.Nomenclature A Ratio of angular amplitude of torsion head to oscillation input signal - G Storage modulus - G Loss modulus - I Moment of inertia of upper platen — torsion head assembly - K Restoring constant of torsion bar - N ₁ First normal-stress difference - R Platen radius - S ⁽ⁱ⁾ Geometric factor in the dynamic property analysis - t ₁ Characteristic time parameter of the Carreau model - X, Y Carreau model parameters - Z () Reimann Zeta function of - Carreau model parameter - Shear rate - Apparent steady-shear viscosity - ^* Complex dynamic viscosity - Dynamic viscosity - Imaginary part of the complex dynamic viscosity - ₀ Zero-shear viscosity - ₀ Cone angle - Carreau model characteristic time - Density of fluid - Shear stress - Phase difference between torsion head and oscillation input signals - ₀ Zero-shear rate first normal-stress coefficient - Oscillatory frequency     

New trends in the analysis of masonry structures

The modern theory of masonry structures has been set up on the hypothesis of no-tension behaviour, with the aim of offering a reference model, independent of materials and building techniques employed. This hypothesis gives rise to inequalities which have to be satisfied by the stress tensor components and, as a dual aspect, to the kinematic behaviour characteristics of media which can be classified as lying between solids and fluids: the structure of the masonry material consists of particles reacting elastically only when in contact. An examination of the plane-stress problem leads us to define, within the prescribed domain under admissible loads, three different subdomains with null, regular, or non-regular principal stress tensors, respectively. As the boundaries of such subdomains are not known a priori, the problem can be classified as a free boundary value problem. The analysis concerns mainly the subdomains where the stress tensor is non-regular; and a non-regularity condition det =0 is added to the equilibrium equations. This condition makes the stress problem isostatic and leads to a violation of Saint-Venant's compliance conditions on strains. Hence there is a need to introduce a strain tensor, not related to the stress tensor, which can be decomposed into an extensional component and a shearing component; we prove that such strains, of the class ^c, are similar to those of the theory of plastic flow. From the point of view of computational analysis the anelastic strains are considered as given distortions; they are computed by means of the Haar-Kármán principle, modified for computational purposes by an idea of Prager and Hodge.

---

Sommario La moderna teoria delle strutture murarie, fondata sulla rigorosa non reagenza a trazione del materiale, ha lo scopo di fornire un modello di riferimento indipendente sia dalle caratteristiche del materiale sia dalle techniche costruttive impiegate. L'ipotesi di non reagenza a trazione si traduce in disuguaglianze che le componenti del tensore di stress devono verificare; dualmente il comportamento caratteristico cinematico può esser classificato di confine, come del resto la stessa statica, tra solidi e fluidi: la struttura ipotizzata del materiale muratura consiste di particelle che reagiscono solo se sono in contatto. L'esame del problema piano porta a definire all'interno del dominio di definizione tre differenti tipi di sub-regioni in cui lo stress è nullo, canonico, o singolare. Poiché le frontiere di queste sub-regioni non sono note a priori il problema può anche essere classificato di frontiera libera. L'analisi concerne fondamentalmente la sub-regione in cui il tensore è non regolare, perché deve verificare anche la condizione det =0. Ciò rende isostatico il problema e conduce anche alla violazione della condizione di integrabilità delle deformazioni. Questo passaggio può essere superato introducendo un tensore di deformazioni a tensioni nulle che si può decomporre in una componente estensionale ed in una componente di scorrimento; si dimostra che queste deformazioni sono equivalenti a quelle che intervengono nella Teoria del flusso plastico. Dal punto di vista computazionale le deformazioni anelastiche sono considerate come distorsioni impresse determinate attraverso il principio di Haar-Kármán modificato, per le techniche computazionali, su idee di Prager e Hodge.

     

Application of ‘Operational Quadrature Methods’ in Time Domain Boundary Element Methods

The usual time domain Boundary Element Method (BEM) contains fundamentalsolutions which are convoluted with time-dependent boundary data andintegrated over the boundary surface. Here, a new approach for theevaluation of the convolution integrals, the so-called OperationalQuadrature Methods developed by Lubich, is presented. In thisformulation, the convolution integral is numerically approximated by aquadrature formula whose weights are determined using the Laplacetransform of the fundamental solution and a linear multisep method. Tostudy the behaviour of the method, the numerical convolution of afundamental solution with a unit step function is compared with theanalytical result. Then, a time domain Boundary Element formulationapplying the Operational Quadrature Methods is derived. For thisformulation only the fundamental solutions in Laplace domain arenecessary. The properties of the new formulation are studied with anumerical example.     

Fourier Transform in Bounded Domains

The functional analysis, the concept of distributionsu in the sense of Schwartz [7] andtheir extension given by Gelfand and Shilov [5]to ultradistributions u ,enables us to find by the means of the Fourier transform a secondlanguage to characterize physical behaviour. Almost any expressionwith physical meaning can be transformed, even if it isformulated in domains with complicated boundaries and evenif it is not integrable.Numerical procedures in the transformed space can bedeveloped in analogy to those well known in engineeringmechanics like the methods of Finite or BoundaryElements (FEM or BEM). Basis of the approaches presentedhere is the analytical representation of characteristicdistribution of a domain and the theorem of Parseval whichstates the invariance of energy in respect to thetransformation. In addition, the concept of thecharacteristic distribution leads to a very simplederivation of the Green-Gauss formulas fundamental for theBoundary or Finite Elements (e.g. [6]).     

Mean and turbulent velocity measurements of supersonic mixing layers

The behavior of supersonic mixing layers under three conditions has been examined by schlieren photography and laser Doppler velocimetry. In the schlieren photographs, some large-scale, repetitive patterns were observed within the mixing layer; however, these structures do not appear to dominate the mixing layer character under the present flow conditions. It was found that higher levels of secondary freestream turbulence did not increase the peak turbulence intensity observed within the mixing layer, but slightly increased the growth rate. Higher levels of freestream turbulence also reduced the axial distance required for development of the mean velocity. At higher convective Mach numbers, the mixing layer growth rate was found to be smaller than that of an incompressible mixing layer at the same velocity and freestream density ratio. The increase in convective Mach number also caused a decrease in the turbulence intensity ( _u/U).List of symbols a speed of sound - b total mixing layer thickness between U ₁ – 0.1 U and U ₂ + 0.1 U - f normalized third moment of u-velocity, f u³/(U)³ - g normalized triple product of u² , g u²/(U)³ - h normalized triple product of u ², h u ²/(U)³ - l _u axial distance for similarity in the mean velocity - l _u axial distance for similarity in the turbulence intensity - M Mach number - M _c convective Mach number (for ₁ = ₂), M _c (U ₁ – U ₂)/(a ₁ + a ₂) - P static pressure - r freestream velocity ratio, r U ₂/U ₁ - Re unit Reynolds number, Re U/ - s freestream density ratio, s ₂/₁ - T _t total temperature - u instantaneous streamwise velocity - u deviation of u-velocity, uu – U - U local mean streamwise velocity - U ₁ primary freestream velocity - U ₂ secondary freestream velocity - average of freestream velocities, (U ₁ + U ₂)/2 - U freestream velocity difference, U U ₁ – U ₂ - instantaneous transverse velocity - v deviation of -velocity, – V - V local mean transverse velocity - x streamwise coordinate - y transverse coordinate - y ₀ transverse location of the mixing layer centerline - ensemble average - ratio of specific heats - boundary layer thickness (y-location at 99.5% of free-stream velocity) - similarity coordinate, (y – y ₀)/b - compressible boundary layer momentum thickness - viscosity - density - standard deviation - dimensionless velocity, (U – U ₂)/U - 1 primary stream - 2 secondary stream A version of this paper was presented at the 11th Symposium on Turbulence, October 17–19, 1988, University of Missouri-Rolla     

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.

     

Die Torsionsspannungen in Wellen mit tiefen Längsnuten

Übersicht MitF(x, y) als Spannungsfunktion einer Welle ohne Nut und(, y) als Potentialfunktion des Quelle-Senke-Systems erhält man Spannungsfunktionen(, y) =F(x, y) –(, y) für Wellen mit tiefen Längsnuten. Es wird gezeigt, daß sich damit die Schubspannungen in den Läufern von Schraubenverdichtern ermitteln lassen.

---

Shearing stresses in shafts with deep longitudinal grooves

Summary The stress functions(, y) of shafts with deep longitudinal grooves may be represented by(, y) =F(x, y) –(, y) whereF(x, y) is the stress function of a cylindrical shaft without grooves and(, y) denotes the potential function of the source-sink system. It is shown that the shearing stresses in rotors of screw-compressors may be obtained in this way.

     

Normal forms for random diffeomorphisms

Given a dynamical system (,, ,) and a random diffeomorphism (): ^{d d} with fixed point at x=0. The normal form problem is to construct a smooth near-identity nonlinear random coordinate transformation h() to make the random diffeomorphism ()=h()^–1() h() as simple as possible, preferably linear. The linearization D(, 0)=:A() generates a matrix cocycle for which the multiplicative ergodic theorem holds, providing us with stochastic analogues of eigenvalues (Lyapunov exponents) and eigenspaces. Now the development runs pretty much parallel to the deterministic one, the difference being that the appearance of turns all problems into infinite-dimensional ones. In particular, the range of the homological operator is in general not closed, making the conceptof-normal form necessary. The stochastic versions of resonance and averaging are developed. The case of simple Lyapunov spectrum is treated in detail.     

On surface waves in generalized thermoelasticity

Knowles' representation theorem for harmonically time-dependent free surface waves on a homogeneous, isotropic elastic half-space is extended to include harmonically time-dependent free processes for thermoelastic surface waves in generalized thermoelasticity of Lord and Shulman and of Green and Lindsay.r , , r , , .This work was done when author was unemployed.     

Free convection heat transfer from an isothermal plate with arbitrary inclination

This paper presents a new formulation for the laminar free convection from an arbitrarily inclined isothermal plate to fluids of any Prandtl number between 0.001 and infinity. A novel inclination parameter is proposed such that all cases of the horizontal, inclined and vertical plates can be described by a single set of transformed equations. Moreover, the self-similar equations for the limiting cases of the horizontal and vertical plates are recovered from the transformed equations by setting=0 and=1, respectively. Heated upward-facing plates with positive and negative inclination angles are investigated. A very accurate correlation equation of the local Nusselt number is developed for arbitrary inclination angle and for 0.001 Pr .

---

Wärmeübertragung bei freier Konvektion an einer isothermen Platte mit beliebiger Neigung

Zusammenfasssung Diese Untersuchung stellt eine neue Formulierung der laminaren freien Konvektion von Flüssigkeiten mit einer Prandtl-Zahl zwischen 0,001 und unendlich an einer beliebig schräggestellten isothermen Platte dar. Ein neuer Neigungsparameter wird eingeführt, so daß alle Fälle der horizontalen, geneigten oder vertikalen Platte von einem einzigen Satz transformierter Gleichungen beschrieben werden können. Die unabhängigen Gleichungen für die beiden Fälle der horizontalen and vertikalen Platte wurden für=0 und=1 aus den transformierten Gleichungen wieder abgeleitet. Es wurden erwärmte aufwärtsgerichtete Platten mit positiven und negativen Neigungswinkeln untersucht. Eine sehr genaue Gleichung wurde für die lokale Nusselt-Zahl bei beliebigen Neigungswinkeln und für 0,001 Pr entwickelt.

---

Nomenclature C _p specific heat - f reduced stream function - g gravitational acceleration - Gr local Grashof number,g(T _w –T _w ) x³/v² - h local heat transfer coefficient - k thermal conductivity - n constant exponent - Nu local Nusselt number,hx/k - p pressure - Pr Prandtl number, v/ - Ra local Rayleigh number,g(T _w –T )J x³/v - T fluid temperature - T _w wall temperature - T temperature of ambient fluid - u velocity component in x-direction - v velocity component in y-direction - x coordinate parallel to the plate - y coordinate normal to the plate Greek symbols thermal diffusivity - thermal expansion coefficient - (Ra¦sin¦)^1/4/( Ra cos()^1/5 - pseudo-similarity variable, (y/) - dimensionless temperature, (T–T )/(T _w–T ) - ( Ra cos)^1/5+(Rasin)^1/4 - v kinematic viscosity - 1/[1 +(Ra cos)^1/5/( Ra¦sin)^1/4] - density of fluid - Pr/(1+Pr) - _w wall shear stress - angle of plate inclination measured from the horizontal - stream function - dimensionless dynamic pressure     

Injection moulding of thermoplastics: Freezing of variable-viscosity fluids.

The injection moulding of thermoplastics involves, during mould filling, flows of hot polymer melts into mould networks, the walls of which are so cold that frozen layers form on them. An analytical study of such flows is presented here for the case when the Graetz and Nahme numbers are large and the Pearson number is small. Thus the flows are developing and temperature differences due to heat generation by viscous dissipation are sufficiently large to cause significant variations in viscosity (but the difference between the entry temperature of the polymer to a specific part of the mould network and the melting temperature of the polymer is not). Br Brinkman number - Gz Graetz number - h half-height of channel or disc - h ^* half-height of polymer melt region in channel or disc - L length of channel or pipe - m viscosity shear-rate exponent - Na Nahme number - p pressure - P pressure drop - Pe Péclet number - Pn Pearson number - Q volumetric flowrate - r radial coordinate in pipe or disc - R radius of pipe - Re Reynolds number - R _i inner radius of disc - R _o outer radius of disc - R ^* radius of polymer melt region in pipe - T temperature - T _ad adiabatic temperature rise - T _e entry polymer melt temperature - T _m melting temperature of polymer - T _max maximum temperature - T ₀ reference temperature - T _w wall temperature - flow-average temperature rise - u _r radial velocity in pipe or disc - u _x axial velocity in channel - u _y transverse velocity in channel or disc - u _z axial velocity in pipe - w width of channel - x axial coordinate in channel or modified radial coordinate in disc - y transverse coordinate in channel or disc - z axial coordinate in pipe - thermal conductivity of molten polymer - thermal conductivity of frozen polymer - scaled dimensionless axial coordinate in channel or pipe or radial coordinate in disc - ₀ undetermined integration constant - heat capacity of molten polymer - viscosity temperature exponent - dimensionless transverse coordinate in channel or disc - ^* dimensionless half-height of polymer melt region in channel or disc - H ^* scaled dimensionless half-height of polymer melt region in channel or disc or radius of polymer melt region in pipe - dimensionless temperature - ^* dimensionless wall temperature - scaled dimensionless temperature - numerical constant - µ viscosity of molten polymer - µ ₀ consistency of molten polymer - dimensionless pressure gradient - scaled dimensionless pressure gradient - density of molten polymer - dimensionless radial coordinate in pipe or disc - _i dimensionless inner radius of disc - ^* dimensionless radius of polymer melt region in pipe - dimensionless streamfunction - scaled dimensionless streamfunction - dummy variable - streamfunction - similarity variable - similarity variable     

On ‘roller-coaster’ experiments for nonlinear oscillators

We consider the dynamics of roller-coaster type experimental models used as analog devices for nonlinear oscillators. It is shown how to chose the shape of the track in order to achieve a desired oscillator equation, in terms of the are length coordinate or its projection onto the horizontal. Explicit calculations are carried out for the linear oscillator, the so-called escape equation, the two-well Duffing oscillator, and the pendulum.     

Mixed convection flow in narrow vertical ducts

The mixed convection flow in a vertical duct is analysed under the assumption that , the ratio of the duct width to the length over which the wall is heated, is small. It is assumed that a fully developed Poiseuille flow has already been set up in the duct before heat from the wall causes this to be changed by the action of the buoyancy forces, as measured by a buoyancy parameter . An analytical solution is derived for the case when the Reynolds numberRe, based on the duct width, is of 0 (1). This is extended to the case whenRe is 0 (^–1) by numerical integrations of the governing equations for a range of values of representing both aiding and opposing flows. The limiting cases, || 1 andR=Re of 0 (1), andR and both large, with of 0 (R ^1/3) are considered further. Finally, the free convection limit, large with R of 0 (1), is discussed.

---

Mischkonvektion in engen senkrechten Rohren

Zusammenfassung Mischkonvektion in einem senkrechten Rohr wird unter der Voraussetzung untersucht, daß das Verhältnis der Rohrbreite zur Länge, über welche die Wand beheizt wird, klein ist. Es wird angenommen, daß sich bereits eine voll entwickelte Poiseuille-Strömung in dem Rohr eingestellt hat, bevor Antriebskräfte, gemessen mit dem Auftriebsparameter , aufgrund der Wandbeheizung die Strömung verändern. Es wird eine analytische Lösung für den Fall erhalten, daß die mit der Rohrbreite als charakteristische Länge gebildete Reynolds-ZahlRe konstant ist. Dies wird mittels einer numerischen Integration der wichtigsten Gleichungen auf den FallRe =f (^–1) sowohl für Gleich- als auch für Gegenstrom ausgedehnt. Weiterhin werden die beiden Grenzfälle betrachtet, wenn || 1 undR=Re konstant ist, sowieR und beide groß mit proportionalR ^1/3. Schließlich wird der Grenzfall der freien Konvektion, großes mit konstantem R, diskutiert.

---

Nomenclature g acceleration due to gravity - Gr Grashof number - G modified Grashof number - h duct width - l length of the heated section of the duct wall - p pressure - Pr Prandtl number - Q flow rate through the duct - Q ₀ heat transfer on the wally=0 - Q ₁ heat transfer on the wally=1 - Re Reynolds number - R modified Reynolds number - T temperature of the fluid - T ₀ ambient temperature - T applied temperature difference - u, velocity component in thex-direction - v, velocity component in they-direction - x, co-ordinate measuring distance along the duct - y, co-ordinate measuring distance across the duct - buoyancy parameter - ₀ modified buoyancy parameter, ₀=R ^–1/3 - coefficient of thermal expansion - ratio of duct width to heated length, =h/l - (non-dimensional) temperature - _w applied temperature on the wally=0 - kinematic viscosity - density of the fluid - ₀ shear stress on the wally=0 - ₁ shear stress on the wally=1 - stream function     

The virtual origin of riblets in an open channel flow

In this paper, a method using the mean velocity profiles for the buffer layer was developed for the estimation of the virtual origin over a riblets surface in an open channel flow. First, the standardized profiles of the mixing length were estimated from the velocity measurement in the inner layer, and the location of the edge of the viscous layer was obtained. Then, the virtual origins were estimated by the best match between the measured velocity profile and the equations of the velocity profile derived from the mixing length profiles. It was made clear that the virtual origin and the thickness of the viscous layer are the function of the roughness Reynolds number. The drag variation coincided well with other results.Nomenclature f _r skin friction coefficient - f _ro skin friction coefficient in smooth channel at the same flow quantity and the same energy slope - g gravity acceleration - H water depth from virtual origin to water surface - H ⁺ u*H/ - H false water depth from top of riblets to water surface - H ⁺ u*H/ - I _e streamwise energy slope - I _b bed slope - k riblet height - k ⁺ u*k/ - l mixing length - l _s standardized mixing length - Q flow quantity - Re Reynolds number volume flow/unit width/v - s riblet spacing - u mean velocity - u* friction velocity = - u* false friction velocity = - y distance from virtual origin - y distance from top of riblet - y ₀ distance from top of riblet to virtual origin - y _v distance from top of riblet to edge of viscous layer - y ⁺ u*y/ - y ⁺ u*y/ - y ₀ ⁺ u*y ₀/ - u ⁺ u*y/ - shifting coefficient for standardization - thickness of viscous layer=y ₀+y - ⁺ u*/ - ⁺ u*/ - eddy viscosity - ridge angle - v kinematic viscosity - density - shear stress     

Designing multi-link function mechanisms by passing conventional boundaries

The mechanism synthesis concerning pure numerical methods can only be realized below given boundaries — and this because of the high power of the system equations. For some time the author has been using graphical methods to pass those boundaries concerning multilink mechanisms. By means of the Drawing-Following Method developed by the author, simple but highly accurate computer programs are now available. As a fundamental rule, reductions are used by the coincidences of several poles. Additionally, the new methods are based on the first deviation of correlation of point positions and relative angles, i.e. by the use of the momentary velocity ratios. By this the torque and force relationships are also immediately given.This paper is a renewed part of a series of lectures given by the author between 25 and 28 September 1972 at the International Centre for Mechanical Sciences, Udine, Italy. The lectures were entitled Theory of mechanisms, dimensional synthesis of linkages by geometrical and graphical methods (10 lecture hours). In the meantime those lectures have been computer programmed using the so-called Drawing-Following Method developed by the author.     

Nonlinear analysis of the forced response of a beam with three mode interaction

An analysis is presented for the primary resonance of a clamped-hinged beam, which occurs when the frequency of excitation is near one of the natural frequencies,_n . Three mode interaction (₂ 3₁ and _{3 1} + 2₂) is considered and its influence on the response is studied. The case of two mode interaction (₂ 3₁) is also considered to compare it with the case of three mode interaction. The straight beam experiencing mid-plane stretching is governed by a nonlinear partial differential equation. By using Galerkin's method the governing equation is reduced to a system of nonautonomous ordinary differential equations. The method of multiple scales is applied to solve the system. Steady-state responses and their stability are examined. Results of numerical investigations show that there exists no significant difference between both modal interactions' influences on the responses.     

Inhomogeneous ‘longitudinal’ plane waves in a deformed elastic material

By definition, a homogeneous isotropic compressible Hadamard material has the property that an infinitesimal longitudinal homogeneous plane wave may propagate in every direction when the material is maintained in a state of arbitrary finite static homogeneous deformation. Here, as regards the wave, homogeneous means that the direction of propagation of the wave is parallel to the direction of eventual attenuation; and longitudinal means that the wave is linearly polarized in a direction parallel to the direction of propagation. In other words, the displacement is of the form u = ncos k(n · x – ct), where n is a real vector. It is seen that the Hadamard material is the most general one for which a longitudinal inhomogeneous plane wave may also propagate in any direction of a predeformed body. Here, inhomogeneous means that the wave is attenuated, in a direction distinct from the direction of propagation; and longitudinal means that the wave is elliptically polarized in the plane containing these two directions, and that the ellipse of polarization is similar and similarly situated to the ellipse for which the real and imaginary parts of the complex wave vector are conjugate semi-diameters. In other words, the displacement is of the form u = {S exp i(S · x – ct)}, where S is a complex vector (or bivector). Then a Generalized Hadamard material is introduced. It is the most general homogeneous isotropic compressible material which allows the propagation of infinitesimal longitudinal inhomogeneous plane circularly polarized waves for all choices of the isotropic directional bivector. Finally, the most general forms of response functions are found for homogeneously deformed isotropic elastic materials in which longitudinal inhomogeneous plane waves may propagate with a circular polarization in each of the two planes of central circular section of the ⁿ -ellipsoid, where is the left Cauchy-Green strain tensor corresponding to the primary pure homogeneous deformation.     

Alloy separation of a binary mixture in a stressed elastic sphere

We consider the static state of a spherical isotropic binary elastic solid mixture whose boundary is given a uniform radial displacement. The elastic volumetric strain energy is given by the classical quadratic form from linear elasticity theory,% MathType!MTEF!2!1!+-% feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqef0uAJj3BZ9Mz0bYu% H52CGmvzYLMzaerbd9wDYLwzYbItLDharqqr1ngBPrgifHhDYfgasa% acOqpw0xe9v8qqaqFD0xXdHaVhbbf9v8qqaqFr0xc9pk0xbba9q8Wq% Ffea0-yr0RYxir-Jbba9q8aq0-yq-He9q8qqQ8frFve9Fve9Ff0dme% GabaqaaiGacaGaamqadaabaeaafiaakeaacaqGxbGaeyypa0ZaaSaa% aeaacaGGXaaabaGaaiOmaaaacaGG7bGaaCOUdiaacIcacaqGJbGaam% ykaiaadIcacaWG0bGaamOCaiaadwgacaWGPaWaaWbaaSqabeaacaGG% YaaaaOGaey4kaSIaaiOmaiabeY7aTjaacIcacaqGJbGaaiykaiaacY% hacaWGLbGaeyOeI0IaaiiiamaalaaabaGaamymaaqaaiaadodaaaGa% aiikaiaadshacaWGYbGaamyzaiaacMcacaWGXaGaaiiFamaaCaaale% qabaGaaiOmaaaakiaac2hacaGGUaaaaa!63E0!\[{\text{W}} = \frac{1}{2}\{ {\mathbf{\kappa }}({\text{c}})(tre)^2 + 2\mu ({\text{c}})|e - \frac{1}{3}(tre)1|^2 \} .\]Here, e is the infinitesimal strain tensor, c[0, 1] is the volumetric concentration of the mixture, and (·) and (·) are the (positive) bulk and shear material moduli, respectively, which are given functions of the concentration. As a function of c and e, the strain energy function is generally nonconvex. Thus, we consider the nonconvex problem of minimizing the potential energy of the body, among all spatial concentration and displacement fields, subject to a given boundary displacement and a fixed amount of component materials. Assuming spherical symmetry, we find that the two component materials must be separated in the optimal state of minimum potential energy. The harder material forms the central core of the sphere, and the softer material is segregated into a surrounding shell. This behavior is remindful of a general notion in metallurgy that in the casting of materials the harder material tends to migrate toward the center.Partial support of the NSF under grant MSS-9024637 and Alliant Techsystems Inc. is gratefully acknowledged.Professor R. Bartelletti of the Università di Pisa is gratefully acknowledged.