Flow in porous media II: The governing equations for immiscible,two-phase flow

The Stokes flow of two immiscible fluids through a rigid porous medium is analyzed using the method of volume averaging. The volume-averaged momentum equations, in terms of averaged quantities and spatial deviations, are identical in form to that obtained for single phase flow; however, the solution of the closure problem gives rise to additional terms not found in the traditional treatment of two-phase flow. Qualitative arguments suggest that the nontraditional terms may be important when / is of order one, and order of magnitude analysis indicates that they may be significant in terms of the motion of a fluid at very low volume fractions. The theory contains features that could give rise to hysteresis effects, but in the present form it is restricted to static contact line phenomena.Roman Letters (, = , , and ) 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² - A _e ^* area of entrances and exits for the-phase contained within a unit cell, m² - g gravity vector, m²/s - H mean curvature of the- interface, m^–1 - H area average of the mean curvature, m^–1 - H – H , deviation of the mean curvature, m^–1 - I unit tensor - K Darcy's law permeability tensor, m² - K permeability tensor for the-phase, m² - K viscous drag tensor for the-phase equation of motion - K viscous drag tensor for the-phase equation of motion - L characteristic length scale for volume averaged quantities, m - characteristic length scale for the-phase, m - n unit normal vector pointing from the-phase toward the-phase (n = –n ) - p _c p – P , capillary pressure, N/m² - 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, m - t time, s - 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² - surface tension of the- interface, N/m - viscous stress tensor for the-phase, N/m² - / kinematic viscosity, m²/s     

Heat transfer from laminar flow in ducts with elliptic cross-section

Summary The cooling of a hot fluid in laminar Newtonian flow through cooled elliptic tubes has been calculated theoretically. Numerical data have been computed for the two values 1.25 and 4 of the axial ratio of the elliptic cross-section . For =1.25 the influence of non-zero thermal resistance between outmost fluid layer and isothermal surroundings has also been investigated. Special attention has been given to the distribution of heat flux around the perimeter; when increases the flux varies more with the position at the circumference. This positional dependence becomes less pronounced, however, as the (position-independent) thermal resistance of the wall increases.Flattening of the conduit, while maintaining its cross-sectional area constant, improves the cooling. Comparison with rectangular pipes shows that this improvement is not as marked with elliptic as with rectangular pipes.Nomenclature A _k =A _{m, n} coefficients of expansion (6) - a, b half-axes of ellipse, b<a - a _p =a _{r, s} coefficients of representation (V) - D hydraulic diameter, = 4S/P; S = cross-sectional area, P = perimeter - D _e equivalent diameter, according to (13) - n coordinate (outward) normal to the tube wall - T temperature of fluid - T _i temperature of fluid at the inlet - T _s temperature of surroundings - v ₀ mean velocity of fluid - v _z longitudinal velocity of fluid - x, y carthesian coordinates coinciding with axes of ellipse - z coordinate in flow direction - , dimensionless half-axes of ellipse, =a/D and =b/D - _t heat transfer coefficient from fluid at bulk temperature to surroundings; equation (11) - _w heat transfer coefficient at the wall; equation (3) - axial ratio of ellipse, = a/b = / - , , , dimensionless coordinates; =x/D, =y/D, =z/D, =n/D - dimensionless temperature, = (T–T _s)/(T _i–T _s) - ₀ cup-mixing mean value of ; equation (10) - thermal conductivity of fluid - _m,n = _k eigenvalue - c volumetric heat capacity of fluid - _{m, n} = _k = _k eigenfunction; equations (6) and (I) - Nu total Nusselt number, = _t D/ - Nusselt number at large distance from the inlet - Nu _w wall Nusselt number, = _w D/, based on _w - Pé Péclet number, = ₀ Dc/     

An experimental study on effects of solid concentration and ζ-potential on pseudoplastic flow of suspensions

The pseudoplastic flow of suspensions, alumina or styrene-acrylamide copolymer particles in water or an aqueous solution of glycerin has been studied by the step-shear-rate method. The relation between the shear rate,D, and the shear stress,, in the step-shear-rate measurements, where the state of dispersion was considered to be constant, was expressed as = AD ^1/2 +CD. The effective solid volume fraction,ø _F, andA were dependent on the shear rate and expressed byø _F =aD ^b andA = D . Combining the above relations, the steady flow curve was expressed by = D ^{1/2 +} + ₀ (1 – a D ^b/0.74)^–1.85 D, where ₀ is the viscosity of the medium.With an increase in solid volume fraction and a decreases in the absolute value of the-potential, the flow behavior of the suspensions changed from Newtonian ( = = b = 0), slightly pseudoplastic ( = b = 0), pseudoplastic ( = 0) to a Bingham-like behavior.The change in viscosity of the medium had an effect on the change in the effective volume fraction.     

Image processing of hydrogen bubble flow visualization for determination of turbulence statistics and bursting characteristics

A technique is described which employs automated image processing of hydrogen-bubble flow visualization pictures to establish local, instantaneous velocity profile information. Hydrogen bubble flow visualization sequences are recorded using a high-speed video system and then digitized, stored, and evaluated by a VAX 11/780 computer. Employing special smoothing and gradient detection algorithms, individual bubble-lines are computer identified, which allows local velocity profiles to be constructed using time-of-flight techniques. It is demonstrated how this techniques may be used to 1) determine local velocity behavior as a function of position and time, 2) evaluate time-averaged turbulence properties, and 3) correlate probe-type turbulent burst detection techniques with the corresponding visualization data.List of symbols Re Reynolds number based on momentum thickness, u / - t ⁺ nondimensional time tu ² / - T VITA variance averaging time period - u shear velocity = - u local instantaneous streamwise velocity,x-direction - u local fluctuating streamwise velocity,x-direction - u ⁺ nondimensional streamwise velocity, /u - local normal velocity,y-direction - w local spanwise velocity,z-direction - x ⁺ nondimensional coordinate in streamwise direction xu /v - y ⁺ nondimensional coordinate normal to wall, yu /v Greek momentum thickness, - kinematic viscosity - _w wall shear stress This paper was presented at the Ninth Symposium on Turbulence, University of Missouri-Rolla, October 1–3, 1984     

Chaos and integrability in a nonlinear wave equation

We consider the parametrized family of equations _tt ,u- _xx u-au+u ₂ ² u=O,x(0,L), with Dirichlet boundary conditions. This equation has finite-dimensional invariant manifolds of solutions. Studying the reduced equation to a four-dimensional manifold, we prove the existence of transversal homoclinic orbits to periodic solutions and of invariant sets with chaotic dynamics, provided that =2, 3, 4,.... For =1 we prove the existence of infinitely many first integrals pairwise in involution.     

Remarks on the bifurcation of solutions of a nonlinear eigenvalue problem

Conclusions We have investigated solutions of equation (3) when ² is an eigenvalue of the linearized operator (13) and when it is not. In Section 4 we have shown that for 0 and ² = _i ² we have exactly two nontrivial solutions which bifurcate to the right of _i ² ; these solutions are shown to exist in an interval ( _i ² , _i ² + ₀). The method of Section 3 may then be used to extend these two solutions to the right of _i ² + ₀ providing that ²= _i ² + ₀ is not an eigenvalue of the linear operator (13) evaluated at = ± ₁. Either a solution can be uniquely extended, or there exists a value of ²where the bifurcation method must be applied again³.While the method used here gives the exact number of solutions bifurcating from _i ² , other problems remain open; for example, it is still not proven that the two bifurcating branches have i zeros, as is the case for Hammerstein operators with oscillation kernels [4]. The conjecture of Odeh and Tadjbakhsh that there are exactly 2(i+1) nontrivial solutions in the interval _i ² < _i ^+1/2 remains un-answered, although it would be proven if one could show that there is no secondary bifurcation as in the cases of Kolodner [7] and Coffman [8].     

Heat transfer in laminar flow in a uniformly porous channel with an applied transverse magnetic field

Summary The steady laminar flow of an incompressible, viscous, and electrically conducting fluid between two parallel porous plates with equal permeability has been discussed by Terrill and Shrestha [6]. In this paper, using the solution of [6] for the velocity field, the heat transfer problems of (i) uniform wall temperature and (ii) uniform heat flux at wall are solved.For small suction Reynolds numbers we find that the Nusselt number, with increasing Reynolds number, increases for case (i) and decreases for (ii).Nomenclature stream function - 2h channel width - x, y distances measured parallel, perpendicular to the channel walls - U velocity of fluid in the x direction at x=0 - V constant velocity of suction at the wall - nondimensional distance, y/h - nondimensional distance, x/h - f() function defined in (1) - density - coefficient of kinematic viscosity - R suction Reynolds number, V h/ - Re channel Reynolds number, 4U h/ - B ₀ magnetic induction - electrical conductivity - M Hartmann number, B ₀ h(/)^1/2 - K constant defined in (3) - A constant defined in (5) - 4R/Re - q local heat flux per unit area at the wall - k thermal conductivity - T temperature of the fluid - X –1/ ln(1–) - C _p specific heat at constant pressure - j current density - Pr Prandtl number, C _p/k - P mass transfer Péclet number, R Pr - Pe mass transfer Péclet number, P/ - T ₀ temperature at x=0 - T _H() temperature in the fully developed region - T _h(X, ) temperature in the entrance region - Y _n () eigenfunctions, uniform wall temperature - _n eigenvalues - e() function defined by (24) - B _n 2/3 _n ² - A _n constants defined by (28) - a _2m constants defined by (30) - F _n () eigenfunctions, uniform wall heat flux - a _n , b _n , c _n , d _n , e _n constants defined by (45) and (48) - S a parameter, U ²/q - h ₁ heat transfer coefficient - T _m mean temperature - Nu Nusselt number - Nu _T Nusselt number, uniform wall temperature - Nu _q Nusselt number, uniform wall heat flux     

An analytical approach to coupled heat and moisture transport in soil

The simultaneous diffusion of heat and moisture through soil is described by two coupled partial differential equations in which the diffusion coefficients are highly non-linear functions of the dependent variables. The system has been regarded as analytically intractable for any generality of coupled flow. However, for an asymptotically steady state, the equations show a marked periodic stability. Computer simulation indicates that the behaviour quickly becomes entrained to input boundary periodicity for any initial state, regardless of the detailed functional form of the diffusion coefficients. This property allows an harmonic series solution to be assembled. Factors such as amplitude decay, phase shift and wave form evolution may be evaluated. The solution is adapted to boundary conditions pertaining to arid soils and the results validated against the 1968 field data of Rose and the 1973 experiment by Jackson.Notation gradient operator - divergence operator - _A amplitude of surface moisture content variation - _l volumetric liquid content, m³/m³ - _c value for moisture content, at which vapour diffusivity decays to zero - _M mean of surface moisture content variation - _s saturation value of moisture content - tortuosity factor, m/m - _i eigenvalues of ₀ - hypothetical thermal conductivity, J/m/sec/K - ₀ density of saturated water vapour, kg/m³ - _l density of liquid water, kg/m³ - _v density of water vapour, kg/m³ - surface tension, kg/sec² - matric potential, m - C volumetric heat capacity, J/m³/K - D ^* molecular diffusivity of water vapour in the porous medium, m²/sec - D _atm molecular diffusivity of water vapour in air, m²/sec - D _TV thermally induced vapour diffusivity, m²/sec/K - D _Tl thermally induced liquid diffusivity, m²/sec/K - D _v isothermal vapour diffusivity, m²/sec - D _l isothermal liquid diffusivity, m²/sec - L latent heat of vaporisation, J/kg - P atmospheric pressure at soil surface,Pa - R gas constant of water vapour, J/kg/K - T temperature,K - T _M mean temperature at surface, K - T _A temperature amplitude at surface, K - g acceleration due to gravity, m/sec² - h relative humidity, dimensionless - p partial pressure of water vapour,Pa - q _v water vapour flux, kg/m²/sec - t time, sec - z depth, (measured downwards), m     

Non-gaussian diffusion modeling in composite porous media by homogenization: Tail effect

In the paper anomalous diffusion appearing in a porous medium composed of two porous components of considerably different diffusion characteristics is examined. The differences in diffusivities are supposed to result either from two medium types being present or from variations in pore size (double porosity media). The long-tail effect is predicted using the homogenization approach based on the application of multiple scale asymptotic developments. It is shown that, if the ratio of effective diffusion coefficients of two porous media is of the order of magnitude smaller or equal O( ²), where is a homogenization parameter, then the macroscopic behaviour of the composite may be affected by the presence of tail-effect. The results of the theoretical analysis were applied to a problem of diffusion in a bilaminate composite. Analytical calculations were performed to show the presence of the long-tail effect in two particular cases.Notations c _i the concentration of chemical species in water within the medium i - D _i the effective diffusion coefficient for the medium i - D _ij ^eff the macroscopic (or effective) diffusion tensor in the composite - ERV the elementary representative volume - h the thickness of the period - l a chracteristic length of the ERV or the periodic cell - L a characteristic macroscopic length - n the volumetric fraction of the material 2 - 1–n the volumetric fraction of the material 1 - N the unit vector normal to - t the time variable - x the macroscopic (or slow) space variable - y the microscopic (or fast) space variable - c _1c ,C _2c ,D _1c ,D _2c the characteristic quantities - T,T _1L ,T _2L ,T _1l ,T _2l the characteristic times - c ₁ ^* ,c ₂ ^* ,D ₁ ^* ,D ₂ ^* ,t ^* the non-dimensional variables - the homogenization parameter - ₁ the domain occupied by the material 1 - ₂ the domain occupied by the material 2 - the interface between the domains 1 and 2 - the total volume of the periodic cell - /x_i the gradient operator - the gradient operator     

A visualization study of the interaction of a free vortex with the wake behind an airfoil

The two-dimensional interaction of a single vortex with a thin symmetrical airfoil and its vortex wake has been investigated in a low turbulence wind tunnel having velocity of about 2 m/s in the measuring section. The flow Reynolds number based on the airfoil chord length was 4.5 × 10³. The investigation was carried out using a smoke-wire visualization technique with some support of standard hot-wire measurements. The experiment has proved that under certain conditions the vortex-airfoil-wake interaction leads to the formation of new vortices from the part of the wake positioned closely to the vortex. After the formation, the vortices rotate in the direction opposite to that of the incident vortex.List of symbols c test airfoil chord - C vortex generator airfoil chord - TA test airfoil - TE test airfoil trailing edge - TE _G vortex generator airfoil trailing edge - t dimensionless time-interval measured from the vortex passage by the test airfoil trailing edge: gDt=(T-T- _TE)·U/c - T time-interval measured from the start of VGA rotation - U free stream velocity - U vortex induced velocity fluctuation - VGA vortex generator airfoil - y distance in which the vortex passes the test airfoil - Z vortex circulation coefficient: Z=/(U · c/2) - vortex generator airfoil inclination angle - vortex circulation - vortex strength: =/2     

Nonstationary vibration of a rotating shaft with nonlinear spring characteristics during acceleration through a critical speed (A critical speed of a summed-and-differential harmonic oscillation)

Nonstationary vibration of a flexible rotating shaft with nonlinear spring characteristics during acceleration through a critical speed of a summed-and-differential harmonic oscillation was investigated. In numerical simulations, we investigated the influence of the angular acceleration , the initial angular position of the unbalance _n and the initial rotating speed on the maximum amplitude. We also performed experiments with various angular accelerations. The following results were obtained: (1) the maximum amplitude depends not only on but also on _n and : (2) when the initial angular position _n changes. the maximum amplitude varies between two values. The upper and lower bounds of the maximum amplitude do not change monotonously for the angular acceleration: (3) In order to always pass the critical speed with finite amplitude during acceleration. the value of must exceed a certain critical value.Nomenclature O-xyz rectangular coordinate system - , ₁, ₁ inclination angle of rotor and its projections to thexy- andyz-planes - I _r polar moment of inertia of rotor - I diametral moment of inertia of rotor - i _r ratio ofI _r toI - dynamic unbalance of rotor - directional angle of fromx-axis - c damping coefficient - spring constant of shaft - N _nt ,N _nt nonlinear terms in restoring forees in ₁ and ₁ directions - ₄ representative angle - a small quantity - V. V _u .V _N potential energy and its components corresponding to linear and nonlinear terms in the restoring forees - directional angle - _n coefficients of asymmetrical nonlinear terms - _n coefficients of symmetrical nonlinear terms - coefficients of asymmetrical nonlinear terms experessed in polar coordinates - coefficients of symmetrical nonlinear terms expressed in polar coordinates - rotating speed of shaft - t time - _n initial angular position of att=0 - p natural frequency - p ₁.p _t natural frequencies of forward and backward precessions - , ₁, ₁ total phases of harmonic, forward precession and backward precession components in summed-and-differential harmonic oscillation - , ₁, ₁ phases of harmonic, forward precession and backward precession components in summed-and-differential harmonic oscillation - P, R _t ,R _b amplitudes of harmonic, forward precession and backward precession components in summed-and-differential harmonic oscillation - difference between phases ( = _f – _u) - acceleration of rotor - initial rotating speed - t _t ,r _b amplitudes of nonstationary oscillation during acceleration - (r _t )_max, (r _b )_max maximum amplitudes of nonstationary oscillation during acceleration - (r ₁ ¹ )_max, (r _b ¹ )_max maximum value of angular acceleration of non-passable case - ₀ critical value over which the rotor can always pass the critical speed - p ₁,p ₂,p ₃,p ₄ natural frequencies of experimental apparatus     

On the boundary conditions at the macroscopic level

We study the problem of the boundary conditions specified at the boundary of a porous domain in order to solve the macroscopic transfer equations obtained by means of the volume-averaging method. The analysis is limited to the case of conductive transport but the method can be extended to other cases. A numerical study enables us to illustrate the theoretical results in the case of a model porous medium. Roman Letters _sf interfacial area of the s-f interface contained within the macroscopic system m² - A _sf interfacial area of the s-f interface contained within the averaging volume m² - C _p mass fraction weighted heat capacity, kcal/kg/K - d _s , d _f microscopic characteristic length m - g vector that maps to _s, m - h vector that maps to _f , m - K _eff effective thermal conductivity tensor, kcal/m s K - l REV characteristic length, m - L macroscopic characteristic length, m - n _fs outwardly directed unit normal vector for the f-phase at the f-s interface - n _e outwardly directed unit normal vector at the dividing surface - T ^* macroscopic temperature field obtained by solving the macroscopic equation (3), K - V averaging volume, m³ - V _s , V _f volume of the considered phase within the averaging volume, m³ - volume of the macroscopic system, m³ - _s , _f volume of the considered phase within the volume of the macroscopic system, m³ - dividing surface, m² Greek Letters _s , _f volume fraction - ratio of thermal conductivities - _s , _f thermal conductivities, kcal/m s K - spatial average density, kg/m³ - microscopic temperature, K - ^* microscopic temperature corresponding to T ^* , K - spatial deviation temperature K - error on the temperature due to the macroscopic boundary conditions, K - spatial average - ^s , ^f intrinsic phase average     

An experimental study of the lateral migration of a droplet in a creeping flow

The distribution of droplets in a plane Hagen-Poiseuille flow of dilute suspensions has been measured by a special LDA technique. This method assumes a well defined relation between the velocity of the droplets and their lateral position in the channel. The measurements have shown that the droplet distribution is non-uniform and depends on the viscosity ratio between the droplets and the carrier liquid. The results have been compared with a theory by Chan and Leal describing the lateral migration of suspended droplets.List of symbols a particle radius, m - d half width of the channel, m - Re flow Reynolds number, = 2 _m · d · /µ - flow velocity, m/s - _m flow velocity at the channel axis, m/s - We Weber number, = 2 _m ^Emphasis>/2 · d · / - x distance from center line (x = 0) of the channel, m - non-dimensional distance from the channel center line, _x ^d - y distance along the channel (y = 0 at channel inlet), m - non-dimensional distance along the channel, = y/2d - non-dimensional, normalized distance along the channel, = · _m · µ/ - interfacial tension, N/m - viscosity ratio of dispersed (droplet) phase to viscosity of continuous phase - µ viscosity of continuous phase, Pa · s - density of continuous phase, kg/m³ - phase density difference, kg/m³ Experiments were performed at Max-Planck-Institut, Göttingen     

Untersuchung des lokalen Stoffüberganges am laminaren und turbulenten Rieselfilm

Zusammenfassung Der lokale Stoffübergang wurde in Abhängigkeit von der Meßlänge, dem Startort und der Zulaufhöhe gemessen. Der Gültigkeitsbereich der Theorie von Nusselt wird ermittelt. Die Reynolds-Zahl nahm Werte zwischen 3,86 und 2496 an. Die örtlich wirkende Hydrodynamik ist entscheidend für das Anwachsen der örtlichen Sherwood-Zahl. Die Genauigkeit aller Versuchsergebnisse kann auf ± 5% abgeschätzt werden.

---

Investigation of the local mass transfer of a laminar and turbulent falling liquid film

The local mass transfer was measured as a function of the measuring length, the starting point and the liquid height above the ring-slot. The range of the Reynolds number was 3,86 Re 2496. The validity of the Nusselt theory and the range of it is shown. The local hydrodynamic is the most important factor of the increase of the local Sherwood number. The accuracy of the measurements is ± 5%.

---

Bezeichnungen a Temperaturleitfähigkeit m²/s=/(c_p) - c Konzentration, c=¯c + c kmol/m³ - c_i0 Konzentration im Flüssigkeitskern kmol/m³ - D Diffusionskoeffizient m²/s - EL-NR Elektrodennummer - Fa Faraday-Konstante A s/kgäq=96,5·10⁶ - g Erdbeschleunigung m/s² - i_G Grenzstromdichte A/m² - u Geschwindigkeit in x-Richtung, u= + u - U Umfang des Rohres m - v Geschwindigkeit in y-Rich- m/stung, v=¯v + v - V^* Volumenstrom m³/s - x Lauflänge, Koordinate in m Strömungsrichtung - x_M Meßlänge für den Stoff-Übergang m - x_ST Startort für den Stoff-Übergang m - y Wegkoordinate senkrecht zur Rohroberfläche m - z Wertigkeit der Elektro-denreaktion kgäq/kmol - ZH Zulaufhöhe m - Wärmeübergangskoeffizient W/m²C - Stoffübergangskoeffizient m/s - Filmdicke m - Wärmeleitfähigkeit W/(mC) - kinematische Viskosität m²/s - Re=u/=V^*/U Reynolds-Zahl - Pr=/a=c_p/ Prandtl-Zahl - Sc=/D Schmidt-Zahl - Nu= / Nusselt-Zahl - Sh= /D Sherwood-Zahl - SHL lokale Sherwood-Zahl - SHM mittlere Sherwood-Zahl - - zeitlich gemittelt - örtlich gemittelt Die Durchführung der Arbeit am Institut für Verfahrens — und Kältetechnik der ETH Zürich bei Prof. Dr. P. Grassmann wurde ermöglicht durch Zuschüsse der Kommission zur Förderung der wissenschaftlichen Forschung und meiner Eltern.     

Time-resolved, 3D, laser-induced fluorescence measurements of fine-structure passive scalar mixing in a tubular reactor

A three-dimensional, time-resolved, laser-induced fluorescence (3D-LIF) technique was developed to measure the turbulent (liquid-liquid) mixing of a conserved passive scalar in the wake of an injector inserted perpendicularly into a tubular reactor with Re=4,000. In this technique, a horizontal laser sheet was traversed in its normal direction through the measurement section. Three-dimensional scalar fields were reconstructed from the 2D images captured at consecutive, closely spaced levels by means of a high-speed CCD camera. The ultimate goal of the measurements was to assess the downstream development of the 3D scalar fields (in terms of the full scalar gradient vector field and its associated scalar energy dissipation rate) in an industrial flow with significant advection velocity. As a result of this advection velocity, the measured 3D scalar field is artificially skewed during a scan period. A method to correct for this skewing was developed, tested and applied. Analysis of the results show consistent physical behaviour.List of symbols  A  Deformation tensor - D_t, D_f  Reactor and injector diameter - L_x, L_y, L_z  Dimensions of the 3D-LIF measurement volume - N_x, N_y, N_z  Number of data samples per measurement volume - Re_m  Reynolds number based on mean velocity - Sc  Schmidt number - f  Focal length - f_c,lens, f_c,array  Cut-off frequency for camera lens and sensor array - f, f  Marginal probability density function for and - f  Joint probability density function of and -  Temporal separation of the 2D data planes -  Temporal resolution of the measurement volume -  Spatial resolution of the measurement volume - ,  Deformation angle and deformation, where =tan -  Fluid energy dissipation rate - ,  Strain limited vorticity and scalar diffusion layers -  Scalar concentration - , _B Kolmogorov and Batchelor length scale - ,  Spherical angles of the scalar gradient vector, -  Kinematic viscosity - _e^–2 Half-thickness (1/e²) of the laser sheet - , _a Kolmogorov and Kolmogorov advection time scales -  Scalar energy dissipation rate -  Scalar diffusivity - 2D, 3D Two- and three-dimensional - DNS Direct numerical simulation - LIF Laser-induced fluorescence - SED Scalar energy dissipation rate - TR Tubular reactor

---

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     

Modellvorstellung zur turbulenten Filmkondensation strömenden Dampfes

Zusammenfassung Der Wärmeübergang bei turbulenter Film kondensation strömenden Dampfes an einer waagerechten ebenen Platte wurde mit Hilfe der Analogie zwischen Impuls-und Wärmeaustausch untersucht. Zur Beschreibung des Impulsaustausches im Film wurde ein Vierbereichmodell vorgestellt. Nach diesem Modell wird die wellige Phasengrenze als starre rauhe Wand angesehen. Die Abhängigkeit einer Schubspannungs-Nusseltzahl von der Film-Reynoldszahl und Prandtlzahl wurde berechnet und dargestellt.

---

A model for turbulent film condensation of flowing vapour

The heat transfer in turbulent film condensation of flowing vapour on a horizontal flat plate was investigated by means of the analogy between momentum and heat transfer. To describe the momentum transfer in the film a four-region model was presented. With this model the wavy interfacial surface is treated as a stiff rough wall. A shear Nusselt number has been calculated and represented as a function of film Reynolds number and Prandtl number.

---

Formelzeichen a Temperaturleitkoeffizient - k Mischungswegkonstante - k _s äquivalente Sandkornrauhigkeit - Nu _x lokale Schubspannungs-Nusseltzahl,Nu _x=x_xv/uw - Pr Prandtlzahl,Pr=v/a - Pr _t turbulente Prandtlzahl,Pr _t =_m/_q - q Wärmestromdichte _q - R Wärmeübergangswiderstand - R_f Wärmeübergangswiderstand des Films - Re _F Reynoldszahl der Filmströmung - T Temperatur - U, V Geschwindigkeitskomponenten des Dampfes in waagerechter und senkrechter Richtung - u, Geschwindigkeitskomponenten des Kondensats in waagerechter und senkrechter Richtung - V Querschwankungsgeschwindigkeit des Kondensats und des Dampfes - u _/gtD Schubspannungsgeschwindigkeit an der Phasengrenze für die Dampfgrenzschicht, uD =(/)^1/2 - u _F Schubspannungsgeschwindigkeit an der Phasengrenze für den Kondensatfilm,u _F =(/)^1/2 - u _w Schubspannungsgeschwindigkeit an der Wand der Kühlplatte,u _w =(w/)^1/2 - y Wandabstand - x Wärmeübergangskoeffizient - gemittelte Kondensatfilmdicke - _s Dicke der zähen Schicht der Filmströmung an der welligen Phasengrenze - ₄ Dicke der zähen Schicht der Filmströmung an der gemittelten glatten Phasengrenze - Wärmeleitzahl - dynamische Viskosität - v kinematische Viskosität - Dichte - Oberflächenspannung - _w Wandschubspannung - Schubspannung an der Phasengrenzfläche - _m turbulente Impulsaustauschgröße - _q turbulente Wärmeaustauschgröße Indizes d Wert des Dampfes - w Wert an der Wand - x lokaler Wert inx - Wert an der Phasengrenze Stoffgrößen ohne Index gelten für das Kondensat     

Über Differentialgleichungen für Wärme- und Stoffaustausch von hygroskopischen Stoffen

Zusammenfassung Die beiden Differentialgleichungssysteme vonKrischer undLykow werden miteinander verglichen. Dabei ergibt sich, daß die in der deutschen und russischen Literatur angewandten mathematischen Modelle der Trocknung von kapillarporösen Körpern praktisch übereinstimmen. Es werden die Transformationsgleichungen der dimensionslosen Kenngrößen angegeben, die die Beziehungen zwischen den beiden Systemen herstellen.

---

The differential equations ofKrischer andLuikow for unsteady internal heat and mass transfer in the porous medium are compared. It is shown, that the mathematical models for drying in the German and Russian literature are equivalent. The transform relations of the non-dimensional parameters between the two models are given.

---

Formelzeichen nach Krischer z laufende Koordinate in Strömungsrichtung in m - R kennzeichnende Abmessung des Körpers in m - t Zeit in h - Raumgewicht bei mittlerer Feuchtigkeit in kg/m³ - _w Teilgewicht des Wassers in 1 m³ Trockengut in kg/m³ - _wa Anfangsfeuchtigkeit in kg/m³ - _D Dampfdichte in kg/m³ - _L Luftvolumen je m³ Trocknungsgut in m³/m³ - Temperatur in °C - _u Umgebungstemperatur in °C - _a Anfangstemperatur in °C - r Verdampfungswärme in kcal/kg - q _E Wärmeentwicklung in kcal/m³ h - c spezifische Wärmekapazität des Trockengutes in kcal/kg grd - Wärmeleitfähigkeit in kcal/m h grd - Feuchtigkeitsleitzahl des Trockengutes in m²/h - wirksame Diffusionszahl von Wasserdampf in Luft in m²/h - Diffusionswiderstandszahl des Trockengutes — - Konstante — - Konstante in kg/m³ grd Formelzeichen nach Lykow X=r/R dimensionslose Koordinate des Körpers;r laufende Koordinate in m;R kennzeichnende Abmessung in m; - Fo=a/R ² Fourier-Zahl;a Temperaturleitzahl in m²/h; Zeit in h - T(X, Fo)=t(r, )– 0/t dimensionslose Temperatur des Körpers im Punkt mit KoordinateX für den ZeitpunktFo;t(r, ) Temperatur in °C; 0 mittlere Anfangstemperatur in °C; t ein vorher angenommener Temperaturunterschied in grd - (X, Fo)= ₀–u(r, )/u dimensionsloses Potential des Stoffübergangs im Punkt mit KoordinateX für den ZeitpunktFo;u(r, ) Feuchtigkeitsgehalt des Trockengutes in kg/kg; ₀ mittlerer Anfangsfeuchtigkeitsgehalt in kg/kg; u ein vorher angenommener Unterschied des Feuchtigkeitsgehalts in kg/kg - Ko= u/c t Kosowitsch-Zahl; Verdampfungswärme in kcal/kg;c spezifische Wärmekapazität des Trockengutes in kcal/kg - Ko*=Ko modifizierte Kosowitsch Zhal; Kenngröße der Zustandsänderung - Pn= t/u Posnowsche Zahl;=a _T ^m /a _m Thermogradientkoeffizient in 1/grd;a _T ^m thermische Stoffübergangszahl (charakterisiert den Stoffstrom unter der Einwirkung von Temperaturgradienten) in m²/h grd;a _m Stoffübergangszahl (charakterisiert den Stoffstrom unter der Einwirkung von Feuchtigkeitsgradienten) in m²/h - Lu=a _m/a Lykowsche Zahl - Ki _q=q _q ()·R/ _q t dimensionsloser Wärmestrom (Kirpitschew-Zahl);q _q() Wärmestrom durch die Körperoberfläche in kcal/m²; _q Wärmeleitfähigkeit in kcal/m² h grd - Ki _m=q _m ()·R/a _{m 0} u dimensionsloser Stoffstrom;q _m() Stoffstrom durch die Körperoberfläche in kg/m² h; ₀ Wichte des Trockengutes in kg/m³