Deformation to breaking of deep water gravity waves

The results of laboratory observations of the deformation of deep water gravity waves leading to wave breaking are reported. The specially developed visualization technique which was used is described. A preliminary analysis of the results has led to similar conclusions than recently developed theories. As a main fact, the observed wave breaking appears as the result of, first, a modulational instability which causes the local wave steepness to approach a maximum and, second, a rapidly growing instability leading directly to the breaking.List of symbols L total wave length - H total wave height - crest elevation above still water level - trough depression below still water level - wave steepness =H/L - crest steepness =/L - trough steepness =/L - F ₁ forward horizontal length from zero-upcross point (A) to wave crest - F ₂ backward horizontal length from wave crest to zero-downcross point (B) - crest front steepness =/F ₁ - crest rear steepness =/F ₂ - vertical asymmetry factor=F ₂/F ₁ (describing the wave asymmetry with respect to a vertical axis through the wave crest) - µ horizontal asymmetry factor=/H (describing the wave asymmetry with respect to a horizontal axis: SWL) - T ₀ wavemaker period - L ₀ theoretical wave length of a small amplitude sinusoïdal wave generated at T _inf0 ^sup–1 frequency - ₀ average wave height     

Critical phenomena and waves in gas-liquid systems

The results of the hydraulic studies of gas-liquid media, wave processes in two-phase media and critical phenomena are described. Some methodological foundations to describe these media and methods to obtain the basic similarity criteria for the hydraulics and gas-dynamics of bubble suspensions are discussed. A detailed consideration is given for the phase transition processes on interfaces and the interface stability. A relation has been revealed between the wave and critical phenomena in two-phase systems.Nomenclature a thermal diffusivity - Ar Archimedes number - B gas constant - C heat capacity - C _p heat capacity at constant pressure - C _v heat capacity at constant volume - c ₀ acoustic velocity in the mixture - c _l acoustic velocity in the liquid - C _f flow resistance coefficient - G mass rate of flow - g gravitational acceleration - L latent heat of evaporation - l initial perturbation width - M Mach number - Nu Nusselt number - P pressure - Pr Prandtl number - R bubble radius - (3P ₀/R ₀ ² _f )^–1 bubble resonance frequency square - T temperature - U medium motion velocity - W heavy phase velocity - W light phase velocity - We Weber number - heat release coefficient - dispersion coefficient - void fraction - adiabatic index - film thickness - dimensionless film thickness - kinematic viscosity coefficient - dynamical viscosity coefficient - dissipation coefficient in the mixture - dispersion parameter - _f liquid phase density - light phase density - heat conductivity - surface tension - frequency, ₀ ² =3P ₀/ _f R ₀ ²     

Transport in ordered and disordered porous media V: Geometrical results for two-dimensional systems

In this paper we continue the geometrical studies of computer generated two-phase systems that were presented in Part IV. In order to reduce the computational time associated with the previous three-dimensional studies, the calculations presented in this work are restricted to two dimensions. This allows us to explore more thoroughly the influence of the size of the averaging volume and to learn something about the use of anon-representative region in the determination of averaged quantities.

Nomenclature

Roman Letters A interfacial area of the– interface associated with the local closure problem, m² - a _i i=1, 2, gaussian probability distribution used to locate the position of particles - l unit tensor - characteristic length for the-phase particles, m - ⁰ reference characteristic length for the-phase particles, m - characteristic length for the-phase, m - _i i=1,2,3 lattice vectors, m - m convolution product weighting function - m _V special convolution product weighting function associated with a unit cell - n _i i=1, 2 integers used to locate the position of particles - n unit normal vector pointing from the-phase toward the-phase - r _p position vector locating the centroid of a particle, m - r gaussian probability distribution used to determine the size of a particle, m - r ₀ characteristic length of an averaging region, m - V averaging volume, m³ - V volume of the-phase contained in the averaging volume,V, m³ - x position of the centroid of an averaging area, m - x ₀ reference position of the centroid of an averaging area, m - y position vector locating points in the-phase relative to the centroid, m Greek Letters V /V, volume average porosity - _{a i} standard deviation ofa _i - _r standard deviation ofr - intrinsic phase average of      

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.     

Similar solutions and the asymptotics of filtration equations

In this paper we consider the asymptotic behavior of solutions of the quasilinear equation of filtration as t. We prove that similar solutions of the equation u _t = (u )_xx asymptotically represent solutions of the Cauchy problem for the full equation u _t = [(u)]_xx if (u) is close to u for small u.     

Asymptotics of the homogenized moduli for the elastic chess-board composite

We find the asymptotic behavior of the homogenized coefficients of elasticity for the chess-board structure. In the chess board white and black cells are isotropic and have Lamé constants (, ,) and (, ) respectively. We assume that the black cells are soft, so 0. It turns out that the Poisson ratio for this composite tends to zero with .     

Calculation of heat transfer accompanying hypersonic three-dimensional flow of air over slender blunt-nosed cones

The effective length method [1, 2] has been used to make systematic calculations of the heat transfer for laminar and turbulent boundary layers on slender blunt-nosed cones at small angles of attack ( + 5° in a separationless hypersonic air stream dissociating in equilibrium (half-angles of the cones 0 20°, angles of attack 0 15°, Mach numbers 5 M 25). The parameters of the gas at the outer edge of the boundary layer were taken equal to the inviscid parameters on the surface of the cones. Analysis of the results leads to simple approximate dependences for the heat transfer coefficients.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 173–177, September–October, 1981.     

Das Einsetzen der Konvektion in Flüssigkeiten über einer beheizten waagerechten Platte

Zusammenfassung Das instationÄre Temperaturfeld oberhalb einer waagerechten, mit konstanter Wärmestromdichte beheizten Platte in verschiedenen Flüssigkeiten wurde experimentell untersucht. —Die Versuche haben ergeben, da\ sich das Einsetzen der InstabilitÄt (Konvektionsbeginn) durch eine Kenngrö\eK=gq _{0 kr} ² /c _p beschreiben lÄ\t. Hierin istg die lokale Fallbeschleunigung, der isobare thermische Ausdehnungskoeffizient, die dynamische ViskositÄt,c _p die isobare spezifische WärmekapazitÄt,q ₀ die Wärmestromdichte und _kr die kritische Zeit, das ist der Zeitabschnitt zwischen dem Einschalten der Heizung und dem ersten Einsetzen der Konvektion.

---

The onset of convection in a horizontal fluid layer heated from below

The unsteady temperature field on the upper side of a horizontal plate heated with constant heat flux is experimentally studied in different fluids. The beginning of the instability (onset of convection) can be described by a dimensionless number given asK=gq _{0 kr} ² /c _p whereg is the local gravitational acceleration, the expansion coefficient at constant pressure,q ₀ the heat flux, _kr the critical time, namely the time interval between the starting of heating and the onset of convection.

     

Rheology of molten polymers

Summary TheCross equation describes the flow of pseudoplastic liquids in terms of an upper and a lower Newtonian viscosity corresponding to infinite and zero shear, and ₀, and of a third material constant related to the mechanism of rupture of linkages between particles in the intermediate, non-Newtonian flow regime, Calculation of of bulk polymers is important, since it cannot be determined experimentally. The equation was applied to the melt flow data of two low density polyethylenes at three temperatures.Using data in the non-Newtonian region covering 3 decades of shear rate to extrapolate to the zero-shear viscosity resulted in errors amounting to about onethird of the measured ₀ values. The extrapolated upper Newtonian viscosity was found to be independent of temperature within the precision of the data, indicating that it has a small activation energy.The ₀ values were from 100 to 1,400 times larger than the values at the corresponding temperatures.The values of were large compared to the values found for colloidal dispersions and polymer solutions, but decreased with increasing temperature. This shows that shear is the main factor in reducing chain entanglements, but that the contribution of Brownian motion becomes greater at higher temperatures.

---

Zusammenfassung Die Gleichung vonCross beschreibt das Fließverhalten von pseudoplastischen Flüssigkeiten durch drei Konstante: Die obereNewtonsche Viskosität (bei sehr hohen Schergeschwindigkeiten), die untereNewtonsche Viskosität ₀ (bei Scherspannung Null), und eine Materialkonstante, die vom Brechen der Bindungen zwischen Partikeln im nicht-Newtonschen Fließbereich abhängt. Die Berechnung von ist wichtig für unverdünnte Polymere, wo man sie nicht messen kann.Die Gleichung wurde auf das Fließverhalten der Schmelzen von zwei handelsüblichen Hochdruckpolyäthylenen bei drei Temperaturen angewandt. Die Werte von ₀, durch Extrapolation von gemessenen scheinbaren Viskositäten im Schergeschwindigkeitsbereich von 10 bis 4000 sec^–1 errechnet, wichen bis 30% von den gemessenen ₀-Werten ab. Die Aktivierungsenergie der war so klein, daß die-Werte bei den drei Temperaturen innerhalb der Genauigkeit der Extrapolation anscheinend gleich waren. Die ₀-Werte waren 100 bis 1400 mal größer als die-Werte.Im Verhältnis zu kolloidalen Dispersionen und verdünnten Polymerlösungen war das der Schmelzen groß, nahm aber mit steigender Temperatur ab. Deshalb wird die Verhakung der Molekülketten hauptsächlich durch Scherbeanspruchung vermindert, aber der Beitrag derBrownschen Bewegung nimmt mit steigender Temperatur zu.

     

Transport in ordered and disordered porous media I: The cellular average and the use of weighting functions

In this work we consider transport in ordered and disordered porous media using singlephase flow in rigid porous mediaas an example. We defineorder anddisorder in terms of geometrical integrals that arise naturally in the method of volume averaging, and we show that dependent variables for ordered media must generally be defined in terms of thecellular average. The cellular average can be constructed by means of a weighting function, thus transport processes in both ordered and disordered media can be treated with a single theory based on weighted averages. Part I provides some basic ideas associated with ordered and disordered media, weighted averages, and the theory of distributions. In Part II a generalized averaging procedure is presented and in Part III the closure problem is developed and the theory is compared with experiment. Parts IV and V provide some geometrical results for computer generated porous media.Roman Letters A interfacial area of the- interface contained within the macroscopic region, m² - A_e area of entrances and exits for the-phase contained within the macroscopic system, m² - g gravity vector, m/s² - I unit tensor - K traditional Darcy's law permeability tensor, m² - L general characteristic length for volume averaged quantities, m - characteristic length (pore scale) for the-phase - (y) weighting function - m(–y) (y), convolution product weighting function - _v special weighting function associated with the traditional averaging volume - N unit normal vector pointing from the-phase toward the-phase - p pressure in the-phase, N/m² - p0 reference pressure in the-phase, N/m² - p traditional intrinsic volume averaged pressure, N/m² - r0 radius of a spherical averaging volume, m - r position vector, m - r position vector locating points in the-phase, m - averaging volume, m³ - V volume of the-phase contained in the averaging volume, m³ - V _cell volume of a unit cell, m³ - v velocity vector in the-phase, m/s - v traditional superficial volume averaged velocity, m/s - x position vector locating the centroid of the averaging volume or the convolution product weighting function, m - y position vector relative to the centroid, m - y position vector locating points in the-phase relative to the centroid, m Greek Letters indicator function for the-phase - Dirac distribution associated with the- interface - V/V, volume average porosity - mass density of the-phase, kg/m³ - viscosity of the-phase, Ns/m²     

Hyperbolic phenomena in a strongly degenerate parabolic equation

We consider the equation u _t =((u) (u _x )) _x , where >0 and where is a strictly increasing function with lim _s = <. We solve the associated Cauchy problem for an increasing initial function, and discuss to what extent the solution behaves qualitatively like solutions of the first-order conservation law u _t = ((u)) _x . Equations of this type arise, for example, in the theory of phase transitions where the corresponding free-energy functional has a linear growth rate with respect to the gradient.     

On Perturbative Expansions to the Stochastic Flow Problem

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

Gap Estimates for Double-Well Schrödinger Operators and Quantum Stabilization of Anharmonic Crystals

The spectrum of the Schrödinger operator of a one-dimensional quantum anharmonic oscillator of mass m is studied. This spectrum consists of simple (nondegenerate) eigenvalues E _n , $$n\in {\mathbb N}_0$$ such that n ^– E _n + as n + with a certain > 1. The gap parameter =min _n (E _n – E _n-1) is in the center of the study. It is proven that this parameter is a continuous function of m; its small mass and large mass asymptotics are found. The influence of the dependence of on m on the stability of systems of interacting quantum anharmonic oscillators is briefly discussed.     

Uniqueness and the maximum principle for quasilinear elliptic equations with quadratic growth conditions

In this paper, we show that the maximum principle holds for quasilinear elliptic equations with quadratic growth under general structure conditions.Two typical particular cases of our results are the following. On one hand, we prove that the equation (1) {ie77-01} where {ie77-02} and {ie77-03} satisfies the maximum principle for solutions in H ¹()L(), i.e., that two solutions u ₁, u ₂H¹() L() of (1) such that u ₁u₂ on , satisfy u ₁u₂ in . This implies in particular the uniqueness of the solution of (1) in H ₀ ¹ ()L().On the other hand, we prove that the equation (2) {ie77-04} where fH^–1() and g(u)>0, g(0)=0, satisfies the maximum principle for solutions uH¹() such that g(u)¦Du|{²L¹(). Again this implies the uniqueness of the solution of (2) in the class uH ₀ ¹ () with g(u)¦Du|{²L¹().In both cases, the method of proof consists in making a certain change of function u=(v) in equation (1) or (2), and in proving that the transformed equation, which is of the form (3) {ie77-05}satisfies a certain structure condition, which using ((v₁ -v ₂)⁺)ⁿ for some n>0 as a test function, allows us to prove the maximum principle.     

Existence of L 1-Connections Between Equilibria of a Semilinear Parabolic Equation

We give a necessary and sufficient condition for the existence of L ¹-connections between equilibria of a semilinear parabolic equation. By an L ¹-connection from an equilibrium ^– to an equilibrium ⁺ we mean a function u(, t) which is a classical solution on the interval (–, T) for some T and blows up at t = T but continues to exist in the space L ¹ for t [T, ) and satisfies u(, t) ^± (in a suitable sense) as t ±. The main tool in our analysis is the zero number.     

Unsteady convective diffusion in a rotating parallel plate channel

The paper presents an exact analysis of the dispersion of a passive contaminant in a viscous fluid flowing in a parallel plate channel driven by a uniform pressure gradient. The channel rotates about an axis perpendicular to its walls with a uniform angular velocity resulting in a secondary flow. Using a generalized dispersion model which is valid for all time, we evaluate the longitudinal dispersion coefficientsK _i (i=1, 2, ...) as functions of time. It is shown thatK ₁=0 andK ₃,K ₄, ... decay rapidly in comparison withK ₂. ButK ₂ decreases with increasing (the dimensionless rotation parameter) for values of upto approximately =2.2. ThereafterK ₂ increases with further increase in and its value gets saturated for large values of (say, 500) and does not change any further with increase in . A physical explanation of this anomalous behaviour ofK ₂ is given.

---

Instationäre konvektive Diffusion in einem rotierenden Parallelplattenkanal

Zusammenfassung In dieser Untersuchung wird eine exakte Analyse der Ausbreitung eines passiven Kontaminierungsstoffes in einer zähen Flüssigkeit gegeben, die, befördert durch einen gleichförmigen Druckgradienten, in einem Parallelplattenkanal strömt. Der Kanal rotiert mit gleichförmiger Winkelgeschwindigkeit um eine zu seinen Wänden senkrechte Achse, wodurch sich eine Sekundärströmung ausbildet. Unter Verwendung eines generalisierten, für alle Zeiten gültigen Dispersionsmodells werden die longitudinalen DispersionskoeffizientenK _i (i=1, 2, ...) als Funktionen der Zeit ermittelt. Es wird gezeigt, daßK ₁=0 gilt und dieK ₃,K ₄, ... gegenüberK ₂ schnell abnehmen.K ₂ nimmt ab, wenn , der dimensionslose Rotationsparameter, bis etwa zum Wert 2,2 ansteigt. Danach wächstK ₂ mit bis auf einem Endwert an, der etwa ab =500 erreicht wird. Dieses anomale Verhalten vonK ₂ findet eine physikalische Erklärung.

---

List of symbols C solute concentration - D molecular diffusivity - K _i longitudinal dispersion coefficients - 2L depth of the channel - P ₀ dimensionless pressure gradient along main flow - Pe Péclet number - q velocity vector - Q _x,Q _y mass flux along the main flow and the secondary flow directions - dimensionless average velocity along the main flow direction - (x, y, z) Cartesian co-ordinates Greek symbols dimensionless rotation parameter - the inclination of side walls withx-axis - kinematic viscosity - fluid density - dimensionless time - angular velocity of the channel - dimensionless distance along the main flow direction - dimensionless distance along the vertical direction - dimensionless solute concentration - integral of the dispersion coefficientK ₂() over a time interval     

A Strongly Degenerate Quasilinear Equation: the Parabolic Case

We prove the existence and uniqueness of entropy solutions of the Neumann problem for the quasilinear parabolic equation u_t=÷ a(u, Du), where a(z,)=f(z,), and f is a convex function of with linear growth as ||||, satisfying other additional assumptions. In particular, this class includes the case where f(z,)=(z)(), >0, and is a convex function with linear growth as ||||.     

The theory of a vibrating-rod densimeter

The paper presents a theory of a device for the accurate determination of the density of fluids over a wide range of thermodynamic states. The instrument is based upon the measurement of the characteristics of the resonance of a circular section tube, or rod, performing steady, transverse oscillations in the fluid. The theory developed accounts for the fluid motion external to the rod as well as the mechanical motion of the rod and is valid over a defined range of conditions. A complete set of working equations and corrections is obtained for the instrument which, together with the limits of the validity of the theory, prescribe the parameters of a practical design capable of high accuracy.Nomenclature A, B, C, D constants in equation (60) - A _j , B _j constants in equation (18) - a _j ⁺ , a _j ^– wavenumbers given by equation (19) - C _f drag coefficient defined in equation (64) - C _f ^/0 , C _f ^/1 components of C _f in series expansion in powers of - c speed of sound - D _b drag force of fluid b - D ₀ coefficient of internal damping - E extensional modulus - force per unit length - F _j ⁺ , F _j ^– constants in equation (24) - f, g functions of defined in equations (56) - G modulus of rigidity - I second moment of area - K constant in equation (90) - k, k constants defined in equations (9) - L half-length of oscillator - Ma Mach number - m _a mass per unit length of fluid a - m _b added mass per unit length of fluid b - m _s mass per unit length of solid - n _j eigenvalue defined in equation (17) - P power (energy per cycle) - P _a , P _b power in fluids a and b - p pressure - R radius of rod or outer radius of tube - R _c radius of container - R _i inner radius of tube - r radial coordinate - T tension - T _visc temperature rise due to heat generation by viscous dissipation - t time - v _r , v radial and angular velocity components - y lateral displacement - z axial coordinate - dimensionless tension - _a dimensionless mass of fluid a - _b dimensionless added mass of fluid b - _b dimensionless drag of fluid b - dimensionless parameter associated with - ₀ dimensionless coefficient of internal damping - dimensionless half-width of resonance curve - dimensionless frequency difference defined in equation (87) - spatial resolution of amplitude - R, , , _s , increments in R, , , _s , - dimensionless amplitude of oscillation - dimensionless axial coordinate - ratio of to - _a , _b ratios of to for fluids a and b - angular coordinate - parameter arising from distortion of initially plane cross-sections - _f thermal conductivity of fluid - dimensionless parameter associated with - viscosity of fluid - _a , _b viscosity of fluids a and b - dimensionless displacement - _j jth component of - density of fluid - _a , _b density of fluids a and b - _s density of tube or rod material - density of fluid calculated on assumption that * - dimensionless radial coordinate - * dimensionless radius of container - dimensionless times - _{rr rr}, _r radial normal and shear stress components - spatial component of defined in equation (13) - _j jth component of - dimensionless streamfunction - ⁰, ¹ components of in series expansion in powers of - phase angle - _r phase difference - _ra , _rb phase difference for fluids a and b - streamfunction - _j jth component defined in equation (22) - dimensionless frequency (based on ) - _a , _b dimensionless frequency in fluids a and b - _s dimensionless frequency (based on _s ) - angular frequency - ₀ resonant frequency in absence of fluid and internal damping - _r resonant frequency in absence of internal fluid - _ra , _rb resonant frequencies in fluids a and b - dimensionless frequency - dimensionless frequency when _a vanishes - dimensionless frequencies when _a vanishes in fluids a and b - dimensionless resonant frequency when _a , _b, _b and ₀ vanish - dimensionless resonant frequency when _a , _b and _b vanish - dimensionless resonant frequency when _b and _b vanish - dimensionless frequencies at which amplitude is half that at resonance     

Transient laminar free convection in horizontal cylinders

A numerical study of laminar natural convection inside uniformly heated, partially or fully filled horizontal cylinders is made. A coordinate transformation which simplifies the discretization of the equations of motion and energy is utilized. The resulting system of partial differential equations with their boundary conditions is solved using central differences for various Prandtl and Grashof numbers for two different grid sizes. The flow in completely filled cylinders for which experimental data are available is predicted. Close agreement between steady-state predictions and experiments is obtained for temperature and velocity profiles as well as for the streamline contours and isotherms. The technique is further demonstrated by solving the transient natural convection flow inside a partially filled horizontal cylinder with an adiabatic free surface and subjected to uniform wall heating.

---

Laminare freie Konvektion in horizontalen Zylindern

Zusammenfassung Es wurde eine numerische Berechnung der laminaren, freien Konvektion in gleichmäßig beheizten, teilweise oder ganz gefüllten, horizontalen Zylindern durchgeführt. Dabei wird eine Koordinatentransformation benützt, welche die Diskretisierung der Bewegungs- und der Energiegleichung vereinfacht. Das so resultierende System von partiellen Differentialgleichungen wird, zusammen mit seinen Randbedingungen, unter Verwendung einer Differenzenmethode für verschiedene Prandtl und Grashof-Zahlen sowie für zwei verschiedene Gittergrößen gelöst. Für den vollständig gefüllten Zylinder, für den experimentelle Daten verfügbar sind, wird die Strömung vorhergesagt. Dabei wird für stationäre Zustände gute Übereinstimmung zwischen Rechnung und Experiment erzielt. Dies gilt sowohl für den Verlauf der Stromlinien als auch für den der Isothermen. Das Verfahren wird weiterhin am Beispiel der Berechnung instationärer, freier Konvektion in einem partiell gefüllten, horizontalen Zylinder demonstriert, wobei eine adiabate, freie Oberfläche und gleichmäßige Beheizung der Wand angenommen sind.

---

Nomenclature g acceleration due to gravity, m/s² - Gr _R ^* modified Grashof number =gqR⁴/kv² - Gr _R Grashof number =gTR³/v² - H heat function vector, dimensionless - k thermal conductivity, W/mK - L(Y) cord length associated with coordinateY, dimensionless - Pr Prandtl number=v/ - q wall heat flux, W/m² - R radius, m - r(X, Y,Z) distance of a boundary point from the reference axis, dimensionless - S vector derived from the flow field solution, dimensionless - T temperature, K - T _w wall temperature, K - T reference temperature, K - t time, s - u, v velocity components inx, y directions, m/s - U, V dimensionless velocity components inX- and Y-direction normalized withU - U reference velocity=gqR²/k or gTR, m/s - V velocity vector, dimensionless - W vorticity vector, dimensionless - W vorticity, dimensionless - x, y, z cartesian coordinates, m - X, Y, Z cartesian coordinates normalized with a reference length, dimensionless Greek letters thermal diffusivity, m²/s - coefficient of thermal expansion, K^–1 - ,,, non-dimensional coordinates in the transformed domain - non-dimensional temperature =(T–T)_k/q^R or T–T/T_w–T - v kinematic viscosity, m²/s - non-dimensional time=v/R² Gr_Rt or v/R² G _R ^* t - angle measured from the bottom of the cylinder, rads - ^* angle measured from the axis on (– ) plane, rads - heat potential, dimensionless - angle of incidence of the heat flux vector, rads - non-dimensional stream function - vector potential, dimensionless - grid size, dimensionless - ² Laplacian operator - gradient vector     

Laminar forced convection in elliptic ducts

The problem of laminar forced convection heat transfer in short elliptical ducts with (i) uniform wall temperature and (ii) prescribed wall heat flux is examined in detail with the well known Lévêque theory of linear velocity profile near the wall. Moreover, consideration is given to the variation of the slope of the linear velocity profile with the position on the duct wall. A correction factor for the temperature dependent viscosity is included. Expressions for the local and average Nusselt numbers and wall temperatures are obtained. For the case of constant heat flux the Nusselt numbers are higher than for constant wall temperature.The results corresponding to the classical Graetz and Purday problems are deduced as special cases.Nomenclature a, b semiaxes of ellipse, b Graetz number (average), Re Pr D _e/Z - h _i ^o local heat transfer coefficient - J _n(x) Bessel function of order n - K thermal conductivity of the fluid - [X] Laplace transform of X - N _u ^o local Nusselt number, h _i ^o D _e/K - perimeter average Nusselt number - overall average Nusselt number - Nu _w wall Nusselt number - Nu Nusselt number at large distance from the inlet - p Laplace transform parameter - Pr Prandtl number, C _a/K - Re Reynolds number, D _e / _a - T temperature of the fluid - T ₁, T _W inlet and wall temperatures, respectively - u _z local isothermal velocity along the axis of the duct - average fluid velocity - x, y, z Cartesian coordinates, z-axis parallel to the axis of the duct (z=0 at duct inlet) - Z length of the duct - thermal diffusivity, K/C - * correction factor for the temperature dependent viscosity - (x) gamma function - coordinate measured normal to the wall of the duct - _a, _w viscosity of fluid at average and wall temperatures - , , z elliptic cylindrical coordinates - density of fluid - (z) heat flux