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

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

Flow visualization of compressible vortex structures using density gradient techniques

Mathematical results are derived for the schlieren and shadowgraph contrast variation due to the refraction of light rays passing through two-dimensional compressible vortices with viscous cores. Both standard and small-disturbance solutions are obtained. It is shown that schlieren and shadowgraph produce substantially different contrast profiles. Further, the shadowgraph contrast variation is shown to be very sensitive to the vortex velocity profile and is also dependent on the location of the peak peripheral velocity (viscous core radius). The computed results are compared to actual contrast measurements made for rotor tip vortices using the shadowgraph flow visualization technique. The work helps to clarify the relationships between the observed contrast and the structure of vortical structures in density gradient based flow visualization experiments.Nomenclature a Unobstructed height of schlieren light source in cutoff plane, m - c Blade chord, m - f Focal length of schlieren focusing mirror, m - C _T Rotor thrust coefficient, T/( ² R ⁴) - I Image screen illumination, Lm/m ² - l Distance from vortex to shadowgraph screen, m - n _b Number of blades - p Pressure,N/m ² - p Ambient pressure, N/m ² - r, , z Cylindrical coordinate system - r _c Vortex core radius, m - Non-dimensional radial coordinate, (r/r _c ) - R Rotor radius, m - Tangential velocity, m/s - Specific heat ratio of air - Circulation (strength of vortex), m ²/s - Non-dimensional quantity, ² 8²p r _c ² - Refractive index of fluid medium - ₀ Refractive index of fluid medium at reference conditions - Gladstone-Dale constant, m ³/kg - Density, kg/m ³ - Density at ambient conditions, kg/m ³ - Non-dimensional density, (/ ) - Rotor solidity, (n _b c/ R) - Rotor rotational frequency, rad/s     

Determination of permeability tensors for two-phase flow in homogeneous porous media: Theory

In this paper we continue previous studies of the closure problem for two-phase flow in homogeneous porous media, and we show how the closure problem can be transformed to a pair of Stokes-like boundary-value problems in terms of pressures that have units of length and velocities that have units of length squared. These are essentially geometrical boundary value problems that are used to calculate the four permeability tensors that appear in the volume averaged Stokes' equations. To determine the geometry associated with the closure problem, one needs to solve the physical problem; however, the closure problem can be solved using the same algorithm used to solve the physical problem, thus the entire procedure can be accomplished with a single numerical code.Nomenclature a a vector that maps V onto , m^-1. - A a tensor that maps V onto . - A area of the - interface contained within the macroscopic region, m². - A area of the -phase entrances and exits contained within the macroscopic region, m². - A area of the - interface contained within the averaging volume, m². - A area of the -phase entrances and exits contained within the averaging volume, m². - Bo Bond number (= (=(–)g²/). - Ca capillary number (= v/). - g gravitational acceleration, m/s². - H mean curvature, m^-1. - I unit tensor. - permeability tensor for the -phase, m². - viscous drag tensor that maps V onto V. - ^* dominant permeability tensor that maps onto v , m². - ^* coupling permeability tensor that maps onto v , m². - characteristic length scale for the -phase, m. - l characteristic length scale representing both and , m. - L characteristic length scale for volume averaged quantities, m. - n unit normal vector directed from the -phase toward the -phase. - n unit normal vector representing both n and n . - n unit normal vector representing both n and n . - P pressure in the -phase, N/m². - p superficial average pressure in the -phase, N/m². - p intrinsic average pressure in the -phase, N/m². - p –p , spatial deviation pressure for the -phase, N/m². - r ₀ radius of the averaging volume, m. - r position vector, m. - t time, s. - v fluid velocity in the -phase, m/s. - v superficial average velocity in the -phase, m/s. - v intrinsic average velocity in the -phase, m/s. - v –v , spatial deviation velocity in the -phase, m/s. - V volume of the -phase contained within the averaging volmue, m³. - averaging volume, m³. Greek Symbols V /, volume fraction of the -phase. - viscosity of the -phase, Ns/m². - density of the -phase, kg/m³. - surface tension, N/m. - (v +v ^T ), viscous stress tensor for the -phase, N/m².     

On natural convection from a vertical plate with a prescribed surface heat flux in porous media

This paper presents a theoretical and numerical investigation of the natural convection boundary-layer along a vertical surface, which is embedded in a porous medium, when the surface heat flux varies as (1 +x ²)), where is a constant andx is the distance along the surface. It is shown that for > -1/2 the solution develops from a similarity solution which is valid for small values ofx to one which is valid for large values ofx. However, when -1/2 no similarity solutions exist for large values ofx and it is found that there are two cases to consider, namely < -1/2 and = -1/2. The wall temperature and the velocity at large distances along the plate are determined for a range of values of .Notation g Gravitational acceleration - k Thermal conductivity of the saturated porous medium - K Permeability of the porous medium - l Typical streamwise length - q _w Uniform heat flux on the wall - Ra Rayleigh number, =gK(q _w /k)l/(v) - T Temperature - Too Temperature far from the plate - u, v Components of seepage velocity in the x and y directions - x, y Cartesian coordinates - Thermal diffusivity of the fluid saturated porous medium - The coefficient of thermal expansion - An undetermined constant - Porosity of the porous medium - Similarity variable, =y(1+x ) ^/3/x ^1/3 - A preassigned constant - Kinematic viscosity - Nondimensional temperature, =(T – T )Ra^1/3 k/qw - Similarity variable, = =y(log_e x)^1/3/x ^2/3 - Similarity variable, =y/x ^2/3 - Stream function     

Heat conduction in a cylinder crossed by an electric current with skin effect

The steady periodic temperature distribution in an infinitely long solid cylinder crossed by an alternating current is evaluated. First, the time dependent and non-uniform power generated per unit volume by Joule effect within the cylinder is determined. Then, the dimensionless temperature distribution is obtained by analytical methods in steady periodic regime. Dimensionless tables which yield the amplitude and the phase of temperature oscillations both on the axis and on the surface of copper or nichrome cylindrical electric resistors are presented.

---

Wärmeleitung in einem stromdurchflossenen Zylinder unter Berücksichtigung des Skin-Effektes

Zusammenfassung Es wird die periodische Temperaturverteilung für den eingeschwungenen Zustand in einem unendlich langen, von Wechselstrom durchflossenen Vollzylinder ermittelt. Zuerst erfolgt die Bestimmung der zeitabhängigen, nichgleichförmigen Energiefreisetzung pro Volumeneinheit des Zylinders infolge Joulescher Wärmeentwicklung und anschließend die Ermittlung der quasistationären Temperaturverteilung auf analytischem Wege. Amplitude und Phasenverzögerung der Temperaturschwingungen werden für die Achse und die Oberfläche eines Kupfer- oder Nickelchromzylinders tabellarisch in dimensionsloser Form mitgeteilt.

---

Nomenclature A integration constant introduced in Eq. (2) - ber, bei Thomson functions of order zero - Bi Biot numberhr ₀/ - c speed of light in empty space - c ₁,c ₂ integration constants introduced in Eq. (46) - c _p specific heat at constant pressure - E electric field - E _z component ofE alongz - E time independent part ofE, defined in Eq. (1) - f function ofs and defined in Eq. (11) - g function ofs and defined in Eq. (37) - h convection heat transfer coefficient - H magnetic field - i imaginary uniti=(–1)^1/2 - I electric current - I _eff effective electric currentI _eff=I/2^1/2 - Im imaginary part of a complex number - J _n Bessel function of first kind and ordern - J electric current density - q _g power generated per unit volume - time average of the power generated per unit volume - time averaged power per unit length - r radial coordinate - R electric resistance per unit length - r ₀ radius of the cylinder - Re real part of a complex number - s dimensionless radial coordinates=r/r ₀ - s, s integration variables - t time - T temperature - time averaged temperature - T _f fluid temperature outside the boundary layer - time average of the surface temperature of the cylinder - u, functions ofs, and defined in Eqs. (47) and (48) - W Wronskian - x position vector - x real variable - Y _n Bessel function of second kind and ordern - z unit vector parallel to the axis of the cylinder - z axial coordinate - · modulus of a complex number - equal by definition Greek symbols amplitude of the dimensionless temperature oscillations - electric permittivity - dimensionless temperature defined in Eq. (16) - ₀, ₁, ₂ functions ofs defined in Eq. (22) - thermal conductivity - dimensionless parameter=(2)^1/2 - magnetic permeability - ₀ magnetic permeability of free space - function of defined in Eq. (59) - dimensionless parameter=c _p/() - mass density - electric conductivity - dimensionless time=t - phase of the dimensionless temperature oscillations - function ofs:= ₁+i ₂ - angular frequency - dimensionless parameter=()^1/2 r ₀     

Flow in porous media III: Deformable media

Stokes flow in a deformable medium is considered in terms of an isotropic, linearly elastic solid matrix. The analysis is restricted to steady forms of the momentum equations and small deformation of the solid phase. Darcy's law can be used to determine the motion of the fluid phase; however, the determination of the Darcy's law permeability tensor represents part of the closure problem in which the position of the fluid-solid interface must be determined.Roman Letters A interfacial area of the- interface contained within the macroscopic system, m² - A interfacial area of the- interface contained within the averaging volume, 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 a unit cell, m² - A _e ^* area of entrances and exits for the-phase contained within a unit cell, m² - E Young's modulus for the-phase, N/m² - e i unit base vectors (i = 1, 2, 3) - g gravity vector, m²/s - H height of elastic, porous bed, m - k unit base vector (=e ₃) - characteristic length scale for the-phase, m - L characteristic length scale for volume-averaged quantities, m - n unit normal vector pointing from the-phase toward the-phase (n = -n ) - p pressure in the-phase, N/m² - P p – g·r, N/m² - r ₀ radius of the averaging volume, m - r position vector, m - t time, s - T total stress tensor in the-phase, N/m² - T ⁰ hydrostatic stress tensor for the-phase, N/m² - u displacement vector for the-phase, m - V averaging volume, m₃ - V volume of the-phase contained within the averaging volume, m³ - v velocity vector for the-phase, m/s Greek Letters V /V, volume fraction of the-phase - mass density of the-phase, kg/m³ - shear coefficient of viscosity for the-phase, Nt/m² - first Lamé coefficient for the-phase, N/m² - second Lamé coefficient for the-phase, N/m² - bulk coefficient of viscosity for the-phase, Nt/m² - T –T ⁰ , a deviatoric stress tensor for the-phase, N/m²     

Transport in ordered and disordered porous media IV: Computer generated porous media for three-dimensional systems

In the method of volume averaging, the difference between ordered and disordered porous media appears at two distinct points in the analysis, i.e. in the process of spatial smoothing and in the closure problem. In theclosure problem, the use of spatially periodic boundary conditions isconsistent with ordered porous media and the fields under consideration when the length-scale constraint,r ₀L is satisfied. For disordered porous media, spatially periodic boundary conditions are an approximation in need of further study.In theprocess of spatial smoothing, average quantities must be removed from area and volume integrals in order to extractlocal transport equations fromnonlocal equations. This leads to a series of geometrical integrals that need to be evaluated. In Part II we indicated that these integrals were constants for ordered porous media provided that the weighting function used in the averaging process contained thecellular average. We also indicated that these integrals were constrained by certain order of magnitude estimates for disordered porous media. In this paper we verify these characteristics of the geometrical integrals, and we examine their values for pseudo-periodic and uniformly random systems through the use of computer generated porous media.

Nomenclature

Roman Letters A interfacial area of the- interface associated with the local closure problem, m² - A _e area of entrances and exits for the-phase contained within the averaging system, m² - a _i i=1, 2, 3 gaussian probability distribution used to locate the position of particles - I unit tensor - L general characteristic length for volume averaged quantities, m - L characteristic length for , m - L characteristic length for , m - 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 the traditional volume average - n _i i=1, 2, 3 integers used to locate the position of particles - n unit normal vector pointing from the-phase toward the-phase - n _e outwardly directed unit normal vector at the entrances and exits of 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 - r position vector, m - r _m support of the weighting functionm, m - averaging volume, m³ - V volume of the-phase contained in the averaging volume,, m³ - x positional vector locating the centroid of an averaging volume, m - x ₀ reference position vector associated with the centroid of an averaging volume, m - y position vector locating points 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 - /L, small parameter in the method of spatial homogenization - standard deviation ofa _i - _r standard deviation ofr - _r intrinsic phase average of      

Pulsatile flow in a tube with wall injection and suction

This paper studies similarity solutions for pulsatile flow in a tube with wall injection and suction. The Navier-Stokes equations are reduced to a system of three ordinary differential equations. Two of the equations represent the effects of suction and injection on the steady flow while the third represents the effects of suction and injection on pulsatile flow. Since the equations for steady flow have been studied previously, the analysis centers on the third equation. This equation is solved numerically and by the method of matched asymptotic expansions. The exact numerical solutions compare well with the asymptotic solutions.The effects of suction and injection on pulsatile flow are the following: a) Small values of suction can cause a resonance-like effect for low frequency pulsatile flow. b) The annular effect still occurs but for large injection or suction the frequency at which this effect becomes dominant depends on the cross-flow Reynolds number. c) The maximum shear stress at the wall is decreased by injection, but may be increased or decreased by suction.Nomenclature a radius of the tube - a ₀ ² i ² - A₀, B₀, C₀, D₀, E₀ constant coefficients appearing in the expression for pressure - b a non-dimensionalized length - b ₀ ² i ^{2 2} - b _k complex coefficients of a power series - B - C ₁, C ₂, D complex constants - d - D _1,2 - f() F(a ^1/2)/aV - f ₀,f ₁,... functions of order one used in asymptotic expansions of f() - F(r) rv _r - g() - G(r) a steady component of velocity in axial direction - h() 4/C₀ a ² H(a ^1/2) - h ₀,h ₁,h ₂,...;l ₀,l ₁,l ₂,... functions of order one used in asymptotic expansions for h() in outer regions - H(r) complex valued function giving unsteady component of velocity - H ₀, H ₁, H ₂, ... K ₀, K ₁, K ₂, ...; L ₀, L ₁, L ₂, ... functions of order one used in asymptotic expansions for h() in inner regions - i - J ₀, J ₁, Y ₀, Y ₁ Bessel functions of first and second kind - k - K Rk/2b ² - O order symbol - p pressure - p ₁(z, t) arbitrary function related to pressure - r radial coordinate - r ₀ (1+16 ^{4 4})^1/4 - R Va/, the crossflow Reynolds number - t time - u() G(r)/V - v _r radial velocity - v _z axial velocity - V constant velocity at which fluid is injected or extracted - z axial coordinate - ² a ²/4 - 4.196 - small parameter; =–2/R (Sect. 4); =–R/2 (Sect. 5); =2/R(Sect. 6) - r ²/a ² - * 0.262 - Arctan (4 ^{2 2}) - , inner variables - kinematic viscosity - b - * zero of g() - density - (r, t) arbitrary function related to axial velocity - frequency     

Stabilitätsrechnung der Rohrströmung einer inkompressiblen viskoelastischen Flüssigkeit zweiter Ordnung

Zusammenfassung Für die Durchsatzströmung im Rohr wird mit Hilfe der klassischen hydrodynamischen Stabilitätstheorie gezeigt, daß die inkompressible Flüssigkeit zweiter Ordnungs = –pI + 2(d + 2t ₁ d ² –t ₀ d) stabil ist gegenüber kleinen rotationssymmetrischen Störungen.

---

Summary For Poiseuille pipe flow it is shown by means of the classical theory of hydrodynamic stability, that the incompressible second-order fluids = –pI + 2(d + 2t ₁ d ² –t ₀ d) is stable with respect to small disturbances of rotational symmetry.

---

Nomenklatur a _n Koeffizienten der Reihenentwicklung - c = /k komplexe Wellengeschwindigkeit - d Deformationsgeschwindigkeitstensor - D, D dimensionsloser Deformationsgeschwindigkeitstensor (Grund- und Störtensor) - e _i kovariante Basis - g Vektor der Erdbeschleunigung - I Einheitstensor - k Wellenzahl - M, O, S, Q, T Funktion vonk, Re, ₀ - p, P, p Gesamt-, Grund-, Stördruck - r, (r, , z) dimensionsloser Ortsvektor (Zylinderkoordinaten) - R Rohrradius - Re =U _M R/ Reynoldszahl - s(s ^*=s*pI) Spannungstensor (Isotroper Anteil des ) - t ₀,t ₁ Stoffzeiten, Parameter der Flüssigkeit zweiter Ordnung - t Zeit - u, U, u Vektor der Gesamt-, Grund-, Störgeschwindigkeit - U _M Maximale Grundgeschwindigkeit - v, V, v Vektor der dimensionslosen Gesamt-, Grund-, Störgeschwindigkeit - w Rotationsgeschwindigkeitstensor - W, W Rotationsgeschwindigkeitstensor, dimensionslos (Grund-, Störtensor) - x (x _r ,x ,x _z ) Ortsvektor (Zylinderkoordinaten) - Viskosität - ₀, ₁ dimensionslose Stoffzeiten - dimensionsloser Druck - Dichte - dimensionslose Zeit - Stromfunktion, dimensionslos - komplexe Frequenz, dimensionslos - = e ⁱ /x _i Nablaoperator (e ⁱ kontravariante Basis) - _* Nablaoperator, dimensionslos - R, I Real-, Imaginärteil Mit 4 Abbildungen     

Transport in ordered and disordered porous media II: Generalized volume averaging

In this paper we develop the averaged form of the Stokes equations in terms of weighting functions. The analysis clearly indicates at what point one must choose a media-specific weighting function in order to achieve spatially smoothed transport equations. The form of the weighting function that produces the cellular average is derived, and some important geometrical theorems are presented.Roman Letters A interfacial area of the- interface associated with the local closure problem, m² - A _e area of entrances and exits for the-phase contained within the averaging system, m² - A _p surface area of a particle, m² - d _p 6V _p/A_p, effective particle diameter, m - g gravity vector, m/s² - I unit tensor - K _m permeability tensor for the weighted average form of Darcy's law, m² - L general characteristic length for volume averaged quantities, m - L _p general characteristic length for volume averaged pressure, m - L characteristic length for the porosity, m - L _v characteristic length for the volume averaged velocity, m - l characteristic length (pore scale) for the-phase - l _i i=1, 2, 3 lattice vectors, m - (y) weighting function - m(–y) (y), convolution product weighting function - _v special weighting function associated with the traditional averaging volume - m _v special convolution product weighting function associated with the traditional averaging volume - m _g general convolution product weighting function - m _V unit cell convolution product weighting function - m _C special convolution product weighting function for ordered media which produces the cellular average - m _D special convolution product weighting function for disordered media - m _M master convolution product weighting function for ordered and disordered media - n unit normal vector pointing from the-phase toward the-phase - p pressure in the-phase, N/m² - p_m superficial weighted average pressure, N/m² - p _m intrinsic weighted average pressure, N/m² - p traditional intrinsic volume averaged pressure, N/m² - p p p _m , spatial deviation pressure, N/m² - r ₀ radius of a spherical averaging volume, m - r _m support of the convolution product weighting function, m - r position vector, m - r position vector locating points in the-phase, m - V 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 - vm superficial weighted average velocity, m/s - v _m intrinsic weighted average velocity, m/s - V volume of the-phase contained in the averaging volume, m³ - V _p volume of a particle, m³ - v traditional superficial volume averaged velocity, m/s - v v p _m spatial deviation 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 - _m m * . weighted average porosity - mass density of the-phase, kg/m³ - viscosity of the-phase, Ns/m² - V /V, volume fraction of the-phase     

Modelling of heat transfer by conduction in a transition region between a porous medium and an external fluid

We study the modelling of purely conductive heat transfer between a porous medium and an external fluid within the framework of the volume averaging method. When the temperature field for such a system is classically determined by coupling the macroscopic heat conduction equation in the porous medium domain to the heat conduction equation in the external fluid domain, it is shown that the phase average temperature cannot be predicted without a generally negligible error due to the fact that the boundary conditions at the interface between the two media are specified at the macroscopic level.Afterwards, it is presented an alternative modelling by means of a single equation involving an effective thermal conductivity which is a function of point inside the interfacial region.The theoretical results are illustrated by means of some numerical simulations for a model porous medium. In particular, temperature fields at the microscopic level are presented.Roman Letters _sf interfacial area of thes-f interface contained within the macroscopic system m² - A _sf interfacial area of thes-f interface contained within the averaging volume m² - C _p mass fraction weighted heat capacity, kcal/kg/K - 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 _s,l _f microscopic characteristic length m - L macroscopic characteristic length, m - n _fs outwardly directed unit normal vector for thef-phase at thef-s interface - n outwardly directed unit normal vector at the dividing surface. - R ₀ REV characteristic length, m - T _i macroscopic temperature at the interface, K - error on the external fluid temperature due to the macroscopic boundary condition, K - 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³. - _mp volume of the porous medium domain, m³ - _ex volume of the external fluid domain, m³ - _s , _f volume of the considered phase within the volume of the macroscopic system, m³ - dividing surface, m² - x, z spatial coordinates Greek Letters _s, _f volume fraction - ratio of the effective thermal conductivity to the external fluid thermal conductivity - ^* macroscopic thermal conductivity (single equation model) kcal/m s K - _s, _f microscopic thermal conductivities, kcal/m s K - spatial average density, kg/m³ - microscopic temperature, K - ^* microscopic temperature corresponding toT ^*, K - spatial deviation temperature K - error in the temperature due to the macroscopic boundary conditions, K - ^* _i macroscopic temperature at the interface given by the single equation model, K - spatial average - ^s , ^f intrinsic phase average.     

Eine analytische Näherungslösung für die eingefrorene laminare Grenzschichtströmung eines binären Gemisches

Zusammenfassung Für die eingefrorene laminare Grenzschichtströmung eines teilweise dissoziierten binären Gemisches entlang einer stark gekühlten ebenen Platte wird eine analytische Näherungslösung angegeben. Danach läßt sich die Wandkonzentration als universelle Funktion der Damköhler-Zahl der Oberflächenreaktion angeben. Für das analytisch darstellbare Konzentrationsprofil stellt die Damköhler-Zahl den Formparameter dar. Die Wärmestromdichte an der Wand bestehend aus einem Wärmeleitungs- und einem Diffusionsanteil wird angegeben und diskutiert. Das Verhältnis beider Anteile läßt sich bei gegebenen Randbedingungen als Funktion der Damköhler-Zahl ausdrücken.

---

An analytical approximation for the frozen laminar boundary layer flow of a binary mixture

An analytical approximation is derived for the frozen laminar boundary layer flow of a partially dissociated binary mixture along a strongly cooled flat plate. The concentration at the wall is shown to be a universal function of the Damkohler-number for the wall reaction. The Damkohlernumber also serves as a parameter of shape for the concentration profile which is presented in analytical form. The heat transfer at the wall depending on a conduction and a diffusion flux is derived and discussed. The ratio of these fluxes is expressed as a function of the Damkohler-number if the boundary conditions are known.

---

Formelzeichen A Atom - A₂ Molekül - C Konstante in Gl. (20) - c₁=1/(2C) Konstante in Gl. (35) - c_p spezifische Wärme bei konstantem Druck - D binärer Diffusionskoeffizient - Ec=u ² /(2h_f) Eckert-Zahl - h spezifische Enthalpie - h_t=h+u²/2 totale spezifische Enthalpie - h _A ⁰ spezifische Dissoziationsenthalpie - K_w Reaktionsgeschwindigkeitskonstante der heterogenen Wandreaktion - 1= /( ) Champman-Rubesin-Parameter - Le=Pr/Sc Lewis-Zahl - M Molmasse - p statischer Druck - Pr= c_pf/ Prandtl-Zahl - q_w Wärmestromdichte an der Wand - q_cw, q_dw Wärmeleitungsbzw. Diffusionsanteil der Wärmestromdichte an der Wand - universelle Gaskonstante - R=/(2M_a) individuelle Gaskonstante der molekularen Komponente - Re_x= u x/ Reynolds-Zahl - Sc=/( D) Schmidt-Zahl - T absolute Temperatur - T_d=h _A ⁰ /R charakteristische Dissoziationstemperatur - u, v x- und y-Komponenten der Geschwindigkeit - U=u/u normierte x-Komponente der Geschwindigkeit - x, y Koordinaten parallel und senkrecht zur Platte Griechische Symbole - =A/ Dissoziationsgrad - Grenzschichtdicke - ₂ Impulsverlustdicke - Damköhler-Zahl der Oberflächenreaktion - =T/T normierte Temperatur - =y/ normierter Wandabstand - Wärmeleitfähigkeit - dynamische Viskosität - , ^* Ähnlichkeitskoordinaten - Dichte - Schubspannung Indizes A auf ein Atom bezogen - M auf ein Molekül bezogen - f auf den eingefrorenen Zustand bezogen - w auf die Wand bezogen - auf den Außenrand der Grenzschicht bezogen     

Mass transport and reaction in catalyst pellets

When heterogeneous chemical reactions take place in porous catalysts, mass transport can occur by bulk diffusion, Knudsen diffusion, and convective transport. Previous studies of these phenomena have been largely based on Maxwell's dusty gas model with the convective transport or Darcy flow added to the diffusive transport. This is done in order to satisfy one of the limiting conditions encountered in the study of flow in porous media. A more fundamental approach consists of the use of the method of volume averaging and the general form of the species momentum equation. For an N-component system, this leads to N independent flux relations to be used in conjunction with the volume-averaged species continuity equations.Roman Letters _A (t) surface area of a species body, m² - a _v interfacial area per unit volume, m^-1 - A _e area of entrances and exits for the -phase contained within the averaging volume, m² - A _K area of the - surface contained within the averaging volume, m² - b _A species A body force, N/kg - b mass average body force, N/kg - B inverse tortuosity tensor for bulk diffusion - c total molar concentration, moles/m³ - c _A species A molar concentration, moles/m³ - _A surface concentration of species A, moles/m² - C_A² intrinsic phase average molar concentration, moles/m³ - c _A – C_A², spatial deviation concentration, moles/m³ - c _A mean molecular speed for species A, m/s - binary diffusion coefficient, m²/s - D _A ^{K, eff} Knudsen diffusion coefficient for species A, m²/s - f vector that maps P _A into P _A , m - g gravitational vector, m/s² - G second order tensor that maps N _A into N _A for free molecule flow conditions - H inverse tortuosity tensor for Knudsen diffusion - I unit tensor - j _A c _A u _A ^* , molar diffusive flux, moles/m²s - K Darcy's Law permeability tensor, m² - L macroscopic length scale, m - L _D diffusive length, m - l characteristic length for the -phase, m - l _A mean free path for species A, m - M _A molecular weight of species A, kg/mole - n outwardly directed unit normal vector - n _K unit normal vector directed from the -phase toward the -phase - n outwardly directed unit normal vector at the entrances and exits of the -phase contained within the averaging volume - N _A c _A v _A molar flux of species A, moles/m²s - N _A intrinsic phase average of the species A molar flux, moles/m²s - \~N _A spatial deviation of the molar flux of species A, moles/m²s - p total pressure, N/m² - P p + , total pressure over and above the hydrostatic pressure, N/m² - P _A partial pressure of species A, N/m² - p _A intrinsic phase average partial pressure, N/m² - P_A – p _A, spatial deviation partial pressure, N/m² - P _A p_A + _AA partial pressure of species A over and above the hydrostatic pressure of species A, N/m² - p _ab diffusive force exerted by species B on species A, N/m³ - universal gas constant, N m/moles K - R _A molar rate of production of species A owing to homogeneous chemical reaction, moles/m³s - molar rate of production of species A owing to heterogeneous chemical reaction, moles/m²s - r _A mass rate of production of species A owing to homogeneous chemical reaction, kg/m³s - r ₀ radius of the averaging volume, m - r position vector, m - t time, s - t _A species stress vector, N/m² - T _A species stress tensor, N/m² - T total stress tensor, N/m² - T temperature, K - T spatial average temperature, K - u _A v _A – v, mass diffusion velocity, m/su _A ^* v_A – v^*, molar diffusion velocity, m/s - u _o velocity of the rigid, solid phase relative to some inertial frame, m/s - v _A species velocity, m/s - v mass average velocity, m/s - v ^* molar average velocity, m/s - v _A ^* species velocity of those molecules of species A generated by chemical reaction, m/s - _A (t) volume of a species A body, m³ - averaging volume, m³ - V volume of the -phase contained within the averaging volume, m³ - V volume of the -phase contained within the averaging volume, m³ - v phase average, mass average velocity, m/s - w arbitrary velocity vector, m/s - x _A c _A /c mole fraction of species A - X _A intrinsic phase average mole fraction - X _A – X _A , spatial deviation mole fraction Greek Letters V/V volume fraction of the -phase - _A sum of all terms in the species A momentum equation that are small compared to the diffusive force, N/m³ - viscosity of the -phase, Ns/m² - _A mass density of species A, kg/m³ - total mass density, kg/m³ - _a species viscous stress tensor, N/m² - total viscous stress tensor, N/m² - tortuosity factor - total body force potential function, Nm/kg - _a species body force potential function, Nm/kg - 3.1416 - _{a a} / mass fraction of species A     

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.     

Theoretical and experimental investigations of low Mach number turbulent cavity flows

Theoretical and experimental investigations are conducted for rectangular cavities of varying sizes in low Mach number turbulent flows. Emphasis is put on the characterization of the generation of self-sustained oscillations in order to develop methods of active control applied to the aeroacoustics of cavity flows. A linearized stability analysis for low Mach number flows is proposed in which the interface of the cavity is modeled by a vorticity layer. Subsequently, the cavity flow is investigated experimentally in a subsonic wind tunnel, using pressure measurements and a phase-locked particle image velocimetry system. Experimental results indicate that the oscillation process is governed by convective waves, with no definite influence of convected vortical structures. The good agreement between the experimental data and the predictions given by the model allows the identification of the oscillations of the cavity interface via the parameters issued from the theoretical analysis.List of symbols c speed of sound, m/s - f frequency, Hz - G Greens function - h displacement of the vorticity layer, m - K_R Rayleigh conductivity of the aperture, m - k₀ acoustic wavenumber, rad/m - k,n, Rossiter formula parameters - M Mach number of the freestream - pressure, Pa - Q volume flux, m³/s - Re_L Reynolds number Re_L=UL/ - Re_x1 Reynolds number Re_x1=Ux₁/ - Re Reynolds number Re=U/ - S frequency based Strouhal number - T period, s - t time, s - U,U,U^± freestream velocity, m/s - v velocity, m/s - W,L,D model cavity dimensions, m - w,l,d analytical cavity dimensions, m - x,y Cartesian coordinates, m - boundary-layer thickness, m - vorticity thickness, m - * boundary-layer displacement thickness, m - ,^± velocity potential, m²/s - acoustic wavelength, m - kinematic viscosity, m²/s - boundary-layer momentum thickness, m - ₀ density, kg/m³ - pulsation based Strouhal number - angular frequency, rad/s - vorticity, s ^–1 - , non-dimensional coordinates x₁,y₁ - non-dimensional displacement h     

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     

Existence of fast traveling waves for some parabolic equations: A dynamical systems approach

We study semilinear elliptic equationsu + cu _x =f(u,u) and ² u + cu _x =f(u,u, ² u) in infinite cylinders (x,y) × ⁿ⁺¹ using methods from dynamical systems theory. We construct invariant manifolds, which contain the set of bounded solutions and then study a singular limitc, where the equations change type from elliptic to parabolic. In particular we show that on the invariant manifolds, the elliptic equation generates a smooth dynamical system, which converges to the dynamical system generated by the parabolic limit equation. Our results imply the existence of fast traveling waves for equations like a viscous reactive 2d-Burgers equation or the Cahn-Hillard equation in infinite strips.     

Flow through a helically coiled annulus

In this paper the flow is studied of an incompressible viscous fluid through a helically coiled annulus, the torsion of its centre line taken into account. It has been shown that the torsion affects the secondary flow and contributes to the azimuthal component of velocity around the centre line. The symmetry of the secondary flow streamlines in the absence of torsion, is destroyed in its presence. Some stream lines penetrate from the upper half to the lower half, and if is further increased, a complete circulation around the centre line is obtained at low values of for all Reynolds numbers for which the analysis of this paper is valid, being the ratio of the torsion of the centre line to its curvature.Nomenclature A =constant - a outer radius of the annulus - b unit binormal vector to C - C helical centre line of the pipe - D rL - g 1000 - K Dean number=Re² - L 1+r sin - M (L ²+ ² r ²)^1/2 - n unit normal vector to C - P, P pressure and nondimensional pressure - p ₀, p pressures of O(1) and O() - Re Reynolds number=aW ₀/ - (r, , s), (r, , s) coordinates and nondimensional coordinates - nonorthogonal unit vectors along the coordinate directions - r ₀ radius of the projection of C - t unit tangent vector to C - V _r, V , V _s velocity components along the nonorthogonal directions - V_r, V, V _s nondimensional velocity components along - W ₀ average velocity in a straight annulus Greek symbols , curvature and nondimensional curvature of C - U, V, W lowest order terms for small in the velocity components along the orthogonal directions t - _r, , _s first approximations to V _r , V, V _s for small - =/=/ - kinematic viscosity - density of the fluid - , torsion and nondimensional torsion of C - , stream function and nondimensional stream function - nondimensional streamfunction for U, V - a inner radius of the annulus After this paper was accepted for publication, a paper entitled On the low-Reynolds number flow in a helical pipe, by C.Y. Wang, has appeared in J. Fluid. Mech., Vol 108, 1981, pp. 185–194. The results in Wangs paper are particular cases of this paper for =0, and are also contained in [9].     

On the closure problem for Darcy's law

In a previous derivation of Darcy's law, the closure problem was presented in terms of an integro-differential equation for a second-order tensor. In this paper, we show that the closure problem can be transformed to a set of Stokes-like equations and we compare solutions of these equations with experimental data. The computational advantages of the transformed closure problem are considerable.Roman Letters A interfacial area of the- interface contained within the macroscopic system, m² - A _e area of entrances and exits for the-phase contained within the macroscopic system, m² - A interfacial area of the- interface contained within the averaging volume, m² - A _e area of entrances and exits for the-phase contained within the averaging volume, m² - B second-order tensor used to respresent the velocity deviation - b vector used to represent the pressure deviation, m^–1 - C second-order tensor related to the permeability tensor, m^–2 - D second-order tensor used to represent the velocity deviation, m² - d vector used to represent the pressure deviation, m - g gravity vector, m/s² - I unit tensor - K C ^–1,–D, Darcy's law permeability tensor, m² - L characteristic length scale for volume averaged quantities, m - characteristic length scale for the-phase, m - l _i i=1, 2, 3, lattice vectors, m - n unit normal vector pointing from the-phase toward the-phase - n _e outwardly directed unit normal vector at the entrances and exits of the-phase - p pressure in the-phase, N/m ² - p intrinsic phase average pressure, N/m² - p –p , spatial deviation of the pressure in the-phase, N/m² - r position vector locating points in the-phase, m - r ₀ radius of the averaging volume, m - t time, s - v velocity vector in the-phase, m/s - v intrinsic phase average velocity in the-phase, m/s - v phase average or Darcy velocity in the \-phase, m/s - v –v , spatial deviation of the velocity in the-phase m/s - V averaging volume, m³ - V volume of the-phase contained in 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²     

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