Diffusion of moisture in deformable porous media. Anisotropic fibrous materials

This paper presents a study on the deformation of anisotropic fibrous porous media subjected to moistening by water in the liquid phase. The deformation of the medium is studied by applying the concept of effective stress. Given the structure of the medium, the displacement of the solid matrix is not taken into account with respect to the displacement of the liquid phase. The transport equations are derived from the model proposed by Narasimhan. The transport coefficients and the relation between the variation in apparent density and effective stress are obtained by test measurements. A numerical model has been established and applied for studying drip moistening of mineral wool samples capable or incapable of deformation.Nomenclature D mass diffusion coefficient [L²t^–1] - e void fraction - g gravity acceleration [Lt^–2] - J mass transfer density [ML^–2t^–1] - K hydraulic conductivity [Lt^–1] - K _s hydraulic conductivity of the solid phase [Lt^–1] - K ^* hydraulic conductivity of the deformable porous medium [Lt^–1] - P pressure of moistening liquid [ML^–1 t^–2] - S degree of saturation - t time [t] - V speed [Lt^–1] - X horizontal coordinate [L] - Z vertical coordinate measured from the bottom of porous medium [L] - z z-coordinate [L] Greek Letters porosity - ₁ total hydric potential [L] - _g gas density [ML^–3] - ₁ liquid density [ML^–3] - ₀ apparent density [ML^–3] - _s density of the solid phase [ML^–3] - density of the moist porous medium [ML^–3] - external load [ML^–1t^–2] - effective stress [ML^–1t^–2] - bishop's parameter - matrix potential or capillary suction [L] Indices g gas - 1 moistening liquid - p direction perpendicular to fiber planes - s solid matrix - t direction parallel to fiber planes - v pore Exponent * movement of solid particles taken into account     

Two-phase geothermal flows with conduction and the connection with Buckley-Leverett theory

Two-phase flows of boiling water and steam in geothermal reservoirs satisfy a pair of conservation equations for mass and energy which can be combined to yield a hyperbolic wave equation for liquid saturation changes. Recent work has established that in the absence of conduction, the geothermal saturation equation is, under certain conditions, asymptotically identical with the Buckley-Leverett equation of oil recovery theory. Here we summarise this work and show that it may be extended to include conduction. In addition we show that the geothermal saturation wave speed is under all conditions formally identical with the Buckley-Leverett wave speed when the latter is written as the saturation derivative of a volumetric flow.Roman Letters C(P, S,q) geothermal saturation wave speed [ms^–1] (14) - c _t (P, S) two-phase compressibility [Pa^–1] (10) - D(P, S) diffusivity [m s^–2] (8) - E(P, S) energy density accumulation [J m^–3] (3) - g gravitational acceleration (positive downwards) [ms^–2] - h _w (P),h _w (P) specific enthalpies [J kg^–1] - J _M (P, S,P) mass flow [kg m^–2 s^–1] (5) - J _E (P, S,P) energy flow [J m^–2s^–1] (5) - k absolute permeability (constant) [m²] - k _w (S),k _s (S) relative permeabilities of liquid and vapour phases - K formation thermal conductivity (constant) [Wm^–1 K^–1] - L lower sheetC<0 in flow plane - m, c gradient and intercept - M(P, S) mass density accumulation [kg m^–3] (3) - O flow plane origin - P(x,t) pressure (primary dependent variable) [Pa] - q volume flow [ms^–1] (6) - S(x, t) liquid saturation (primary dependent variable) - S ^*(x,t) normalised saturation (Appendix) - t time (primary independent variable) [s] - T temperature (degrees Kelvin) [K] - T _sat(P) saturation line temperature [K] - TdT _sat/dP saturation line temperature derivative [K Pa^–1] (4) - T _c ,T _D convective and diffusive time constants [s] - u _w (P),u _s (P),u _r (P) specific internal energies [J kg^–1] - U upper sheetC > 0 in flow plane - U(x,t) shock velocity [m s^–1] - x spatial position (primary independent variable) [m] - X representative length - x, y flow plane coordinates - z depth variable (+z vertically downwards) [m] Greek Letters _P , _S remainder terms [Pa s^–1], [s^–1] - double-valued saturation region in the flow plane - h =h _s –h _w latent heat [J kg^–1] - = _w – _s density difference [kg m^–3] - line envelope - =D _K /D ₀ diffusivity ratio - porosity (constant) - _w (P), _s (P), _t (P, S) dynamic viscosities [Pa s] - v _w (P),v _s (P) kinematic viscosities [m²s^–1] - v ₀ =kh/KT kinematic viscosity constant [m² s^–1] - ₀ =v ₀ dynamic viscosity constant [m² s^–1] - _w (P), _s (P) density [kg m^–3] Suffixes r rock matrix - s steam (vapour) - w water (liquid) - t total - av average - 0 without conduction - K with conduction     

On the laminar forced convection with axial conduction in a circular tube with exponential wall heat flux

The values of the fully developed Nusselt number for laminar forced convection in a circular tube with axial conduction in the fluid and exponential wall heat flux are determined analytically. Moreover, the distinction between the concepts of bulk temperature and mixing-cup temperature, at low values of the Peclet number, is pointed out. Finally it is shown that, if the Nusselt number is defined with respect to the mixing-cup temperature, then the boundary condition of exponentially varying wall heat flux includes as particular cases the boundary conditions of uniform wall temperature and of convection with an external fluid.

---

Über laminare Zwangskonvektion mit Längswärmeleitung in einem Kreisrohr mit exponentiell veränderlichem Wandwärmefluß

Zusammenfassung Es werden die Endwerte der Nusselt-Zahlen für vollausgebildete laminare Zwangskonvektion in einem Kreisrohr mit Längswärmeleitung und exponentiell veränderlichem Wandwärmefluß analytisch ermittelt. Besondere Betonung liegt auf dem Unterschied zwischen den Konzepten für die Mittel- und die Mischtemperatur bei niedrigen Peclet-Zahlen. Schließlich wird gezeigt, daß bei Definition der Nusselt-Zahl bezüglich der Mischtemperatur die Randbedingung exponentiell veränderlichen Randwärmeflusses die Spezialfälle konstanter Wandtemperatur und konvektiven Wärmeaustausches mit einem umgebenden Fluid einschließt.

---

Nomenclature A _n dimensionless coefficients employed in the Appendix - Bi Biot numberBi=h _e r ₀/ - c _n dimensionless coefficients defined in Eq. (17) - c _p specific heat at constant pressure of the fluid within the tube, [J kg^–1 K^–1] - f solution of Eq. (15) - h ₁,h ₂ specific enthalpies employed in Eqs. (2) and (4), [J kg^–1] - h _e convection coefficient with a fluid outside the tube, [W m^–2 K^–1] - rate of mass flow, [kg s^–1] - Nu bulk Nusselt number,2r ₀ q _w /[(T _w –T _b )] - Nu _H fully developed value of the bulk Nusselt number for the boundary condition of uniform wall heat flux - Nu _T fully developed value of the bulk Nusselt number for the boundary condition of uniform wall temperature - Nu ^* mixing Nusselt number,2r ₀ q _w /[(T _w –T _m )] - Nu _C ^* fully developed value of the mixing Nusselt number for the boundary condition of convection with an external fluid - Nu _H ^* fully developed value of the mixing Nusselt number for the boundary condition of uniform wall heat flux - Nu _T ^* fully developed value of the mixing Nusselt number for the boundary condition of uniform wall temperature - Pe Peclet number, 2r ₀/ - q ₀ wall heat flux atx=0, [W m^–2] - q _w wall heat flux, [W m^–2] - r radial coordinate, [m] - r ₀ radius of the tube, [m] - s dimensionless radius,s=r/r ₀ - T temperature, [K] - T ₀ temperature constant employed in Eq. (14), [K] - T reference temperature of the fluid external to the tube, [K] - T _b bulk temperature, [K] - T _m mixing or mixing-cup temperature, [K] - T _w wall temperature, [K] - u velocity component in the axial direction, [m s^–1] - mean value ofu, [m s^–1] - x axial coordinate, [m] Greek symbols thermal diffusivity of the fluid within the tube, [m² s^–1] - exponent in wall heat flux variation, [m^–1] - dimensionless parameter - dimensionless temperature =(T _w –T)/(T _w –T _b ) - ^* dimensionless temperature ^*=(T _w –T)/(T _w –T _m ) - thermal conductivity of the fluid within the tube, [W m^–1 K^–1] - density of the fluid within the tube, [kg m^–3]     

Convective heat transfer in the thermal entrance region of parallel-flow noncircular duct heat exchanger arrays

Convective heat transfer properties of a hydrodynamically fully developed flow, thermally developing flow in a parallel-flow, and noncircular duct heat exchanger passage subject to an insulated boundary condition are analyzed. In fact, due to the complexity of the geometry, this paper investigates in detail heat transfer in a parallel-flow heat exchanger of equilateral-triangular and semicircular ducts. The developing temperature field in each passage in these geometries is obtained seminumerically from solving the energy equation employing the method of lines (MOL). According to this method, the energy equation is reformulated by a system of a first-order differential equation controlling the temperature along each line.Temperature distribution in the thermal entrance region is obtained utilizing sixteen lines or less, in the cross-stream direction of the duct. The grid pattern chosen provides drastic savings in computing time. The representative curves illustrating the isotherms, the variation of the bulk temperature for each passage, and the total Nusselt number with pertinent parameters in the entire thermal entry region are plotted. It is found that the log mean temperature difference (T _LM), the heat exchanger effectiveness, and the number of transfer units (NTU) are 0.247, 0.490, and 1.985 for semicircular ducts, and 0.346, 0.466, and 1.345 for equilateral-triangular ducts.

---

Konvektiver Wärmeübergang im thermischen Einlaufgebiet von Gleichstromwärmetauschern mit nichtkreisförmigen Strömungskanälen

Zusammenfassung Die Untersuchung bezieht sich auf das konvektive Wärmeübertragungsverhalten eines Gleichstromwärmetauschers mit nichtkreisförmigen Strömungskanälen bei hydraulisch ausgebildetet, thermisch einlaufender Strömung unter Aufprägung einer adiabaten Randbedingung. Zwei Fälle komplizierter Geometrie, nämlich Kanäle mit gleichseitig dreieckigen und halbkreisförmigen Querschnitten, werden bezüglich des Wärmeübergangsverhaltens bei Gleichstromführung eingehend analysiert. Das sich entwickelnde Temperaturfeld in jedem Kanal von der eben spezifizierten Querschnittsform wird halbnumerisch durch Lösung der Energiegleichung unter Einsatz der Linienmethode (MOL) erhalten. Dieser Methode entsprechend erfolgt eine Umformung der Energiegleichung in ein System von Differentialgleichungen erster Ordnung, welches die Temperaturverteilung auf jeder Linie bestimmt.Die Temperaturverteilung im Einlaufgebiet wird unter Vorgabe von 16 oder weniger Linien über dem Kanalquerschnitt erhalten, wobei die gewählte Gitteranordnung drastische Einsparung an Rechenzeit ergibt. Repräsentative Kurven für das Isothermalfeld, den Verlauf der Mischtemperatur für jeden Kanal und die Gesamt-Nusseltzahl als Funktion relevanter Parameter im gesamten Einlaufgebiet sind in Diagrammform dargestellt. Es zeigt sich, daß die mittlere logarithmische Temperaturdifferenz (T _LM), der Wärmetauscherwirkungsgrad und die Anzahl der Übertragungseinheiten (NTU) folgende Werte annehmen: 0,247, 0,490 und 1,985 für halbkreisförmige Kanäle sowie 0,346, 0,466 und 1,345 für gleichseitig dreieckige Kanäle.

---

Nomenclature A cross sectional area [m²] - a characteristic length [m] - C _c specific heat of cold fluid [J kg^–1 K^–1] - C _h specific heat of hot fluid [J kg^–1 K^–1] - C _p specific heat [J kg^–1 K^–1] - C _r specific heat ratio,C _r=C _c/C_h - D _h hydraulic diameter of duct [m] - f friction factor - k thermal conductivity of fluid [Wm^–1 K^–1] - L length of duct [m] - m mass flow rate of fluid [kg s^–1] - N factor defined by Eq. (20) - NTU number of transfer units - Nu _{x, T} local Nusselt number, Eq. (19) - P perimeter [m] - p pressure [KN m^–2] - Pe Peclet number,RePr - Pr Prandtl number,/ - Q _T total heat transfer [W], Eq. (13) - Q ideal heat transfer [W], Eq. (14) - Re Reynolds number,D _h/ - T temperature [K] - T _b bulk temperature [K] - T _e entrance temperature [K] - T _w circumferential duct wall temperature [K] - u, U dimensional and dimensionless velocity of fluid,U=u/u - , dimensional and dimensionless mean velocity of fluid - w generalized dependent variable - X dimensionless axial coordinates,X=D _h ² /a ² x* - x, x* dimensional and dimensionless axial coordinate,x*=x/D _hPe - y, Y dimensional and dimensionless transversal coordinates,Y=y/a - z, Z dimensional and dimensionless transversal coordinates,Z=z/a Greek symbols thermal diffusivity of fluid [m² s^–1] - * right triangular angle, Fig. 2 - independent variable - T _LM log mean temperature difference of heat exchanger - effectiveness of heat exchanger - generalized independent variable - dimensionless temperature - _b dimensionless bulk temperature - dynamic viscosity of fluid [kg m^–1 s^–1] - kinematic viscosity of fluid [m² s^–1] - density of fluid [kg m^–3] - heat transfer efficiency, Eq. (14) - generalized dependent variable     

Mixed convection flow over a vertical plate with localized heating (cooling), magnetic field and suction (injection)

An analysis is carried out to study the effects of localized heating (cooling), suction (injection), buoyancy forces and magnetic field for the mixed convection flow on a heated vertical plate. The localized heating or cooling introduces a finite discontinuity in the mathematical formulation of the problem and increases its complexity. In order to overcome this difficulty, a non-uniform distribution of wall temperature is taken at finite sections of the plate. The nonlinear coupled parabolic partial differential equations governing the flow have been solved by using an implicit finite-difference scheme. The effect of the localized heating or cooling is found to be very significant on the heat transfer, but its effect on the skin friction is comparatively small. The buoyancy, magnetic and suction parameters increase the skin friction and heat transfer. The positive buoyancy force (beyond a certain value) causes an overshoot in the velocity profiles.A mass transfer constant - B magnetic field - C_fx skin friction coefficient in the x-direction - C_p specific heat at constant pressure, kJ.kg^–1.K - C_v specific heat at constant volume, kJ.kg^–1.K^–1 - E electric field - g acceleration due to gravity, 9.81 m.s^–2 - Gr Grashof number - h heat transfer coefficient, W.m².K^–1 - Ha Hartmann number - k thermal conductivity, W.m^–1.K - L characteristic length, m - M magnetic parameter - Nu_x local Nusselt number - p pressure, Pa, N.m^–2 - Pr Prandtl number - q heat flux, W.m^–2 - Re Reynolds number - Re_m magnetic Reynolds number - T temperature, K - T_o constant plate temperature, K - u,v velocity components, m.s^–1 - V characteristic velocity, m.s^–1 - x,y Cartesian coordinates - thermal diffusivity, m².s^–1 - coefficient of thermal expansion, K^–1 - , transformed similarity variables - dynamic viscosity, kg.m^–1.s^–1 - ₀ magnetic permeability - kinematic viscosity, m².s^–1 - density, kg.m^–3 - buoyancy parameter - electrical conductivity - stream function, m².s^–1 - dimensionless constant - dimensionless temperature, K - w, conditions at the wall and at infinity     

Injection moulding of plastics

Summary In continuation of a previous investigation a simple analytical expression is derived in closed form for the thickness distribution of the freeze-off layer which is vitrified at the (flat) wall of an oblong rectangular cavity. As has been pointed out previously, this layer is marked for amorphous polymers by the molecular orientation (birefringence pattern) in the moulded sample. One can show that a more detailed study with the aid of the coupled equations of energy and of motion will not furnish essential improvements. Problems of polymer physics like glass transition or crystallization kinetics at extreme rates of cooling and shearing must be solved first.

---

Zusammenfassung In Fortsetzung einer früheren Untersuchung wurde ein einfacher analytischer Ausdruck in geschlossener Form für die Dickenverteilung der eingefrorenen Schicht abgeleitet, die an der (flachen) Wand eines langgestreckten rechteckigen Formnestes während des Einspritzvorgangs glasig erstarrt. Wie früher auseinandergesetzt wurde, wird diese Schicht bei amorphen Polymeren durch die Molekülorientierung (Doppelbrechungsmuster) im gespritzten Formteil markiert. Man kann zeigen, daß eine eingehendere Studie mit Hilfe der gekoppelten Energie- und Impulsgleichungen keine essentiellen Verbesserungen bringt. Probleme der Polymerphysik, wie Glasübergang oder Kristallisationskinetik bei extremen Abkühlungs- und Schergeschwindigkeiten, müssen erst gelöst werden.

---

List of Symbols a heat diffusivity of polymer melt (averaged overT) [m²s^–1] - B breadth of mould cavity [m] - Br Brinkman number ( ) - c heat capacity of polymer melt (averaged overT) [J kg^–1 K^–1] - F ₀ Fourier number (at _i/4H ²) - h heat transfer coefficient by melt flow [J K^–1 s^–1 m^–2] - h heat transfer coefficient by layer growth [J K^–1 s^–1 m^–2] - H half height of mould cavity [m] - L length of mould cavity [m] - n exponent in eq. [18] (= 0.417) - Nu Nußelt number (2Hh/) - P pressure gradientdP/dz in mould [N m^–3] - t time [s] - t _i injection time [s] - T _g glass transition temperature of polymer [K] - T _i injection temperature of polymer melt [K] - T _l stagnation temperature [K] - T _m mould wall temperature [K] - speed of flow front during mould filling [m s^–1] - x coordinate perpendicular to mould wall [m] - z coordinate in the injection direction [m] - thickness of stagnant layer (atT _l) [m] - ₀ optically detectable freeze-off thickness [m] - ⁺ apparent layer thickness (atT _i) [m] - dimensionless freeze-off thickness (= ₀/2H) - dimensionless distance from entrance (=z/L) - _m dimensionless coordinate of layer maximum - _g dimensionless temperature (= (T _i –T _l)/(T _g –T _m)) - _i dimensionless temperature (= (T _i –T _l)/(T _i –T _m)) - _l dimensionless temperature (= (T _i –T _l)/(T _l –T _m)) - _i viscosity of polymer atT _i [N s m^–3] - _l viscosity of polymer atT _l [N s m^–3] - heat conductivity of polymer melt (averaged) [J K^–1 s^–1 m^–1] - density of polymer melt (averaged) [kg m^–3] - dimensionless time (eq. [11]) - ⁺ dimensionless parameter (eqs. [19a] and [19b]) - dimensionless layer thickness (eq. [12]) - ⁺ dimensionless parameter (eq. [20a]) - dimensionless parameter (eqs. [11a] and [11b]) Formerly at Delft University of Technology, Delft (The Netherlands).Paper presented at the Conference on Chemical Engineering Rheology, Annual Meeting of the Deutsche Rheologische Gesellschaft in Aachen, March 5–7, 1979.With 3 figures and 1 table     

The effects of span position of winglet vortex generator on local heat/mass transfer over a three-row flat tube bank fin

The naphthalene sublimation method was used to study the effects of span position of vortex generators (VGs) on local heat transfer on three-row flat tube bank fin. A dimensionless factor of the larger the better characteristics, JF, is used to screen the optimum span position of VGs. In order to get JF, the local heat transfer coefficient obtained in experiments and numerical method are used to obtain the heat transferred from the fin. A new parameter, named as staggered ratio, is introduced to consider the interactions of vortices generated by partial or full periodically staggered arrangement of VGs. The present results reveal that: VGs should be mounted as near as possible to the tube wall; the vortices generated by the upstream VGs converge at wake region of flat tube; the interactions of vortices with counter rotating direction do not effect Nusselt number (Nu) greatly on fin surface mounted with VGs, but reduce Nu greatly on the other fin surface; the real staggered ratio should include the effect of flow convergence; with increasing real staggered ratio, these interactions are intensified, and heat transfer performance decreases; for average Nu and friction factor (f), the effects of interactions of vortices are not significant, f has slightly smaller value when real staggered ratio is about 0.6 than that when VGs are in no staggered arrangement. A cross section area of flow passage [m²] - A _mim minimum cross section area of flow passage [m²] - a width of flat tube [m] - b length of flat tube [m] - B _pT lateral pitch of flat tube: B _pT = S ₁/T _p - d _h hydraulic diameter of flow channel [m] - D _naph diffusion of naphthalene [m²/s] - f friction factor: f = pd _h/(Lu ² _max/2) - h mass transfer coefficient [m/s] - H height of winglet type vortex generators [m] - j Colburn factor [–] - JF a dimensionless ratio, defined in Eq. (23) [–] - L streamwise length of fin [m] - L _PVG longitudinal pitch of vortex generators divided by fin spacing: L _pVG = l _VG/T _p - l _VG pitch of in-line vortex generators [m] - m mass [kg] - m mass sublimation rate of naphthalene [kg/m²·s] - Nu Nusselt number: Nu = d _h/ - P pressure of naphthalene vapor [Pa] - p non-dimensional pitch of in-line vortex generators: p = l _VG/S ₂ - Pr Prandtl number [–] - Q heat transfer rate [W] - R universal gas constant [m²/s²·K] - Re Reynolds number: Re = ·u _max·d _h/ - S ₁ transversal pitch between flat tubes [m] - S ₂ longitudinal pitch between flat tubes [m] - Sc Schmidt number [–] - Sh Sherwood number [–]: Sh = hd _h/D _naph - Sr staggered ratio [–]: Sr = (2Hsin – C)/(2Hsin) - T _p fin spacing [m] - T temperature [K] - u _max maximum velocity [m/s] - u average velocity of air [m/s] - V volume flow rate of air [m³/s] - x,y,z coordinates [m] - z sublimation depth[m] - heat transfer coefficient [W/m²·K] - heat conductivity [W/m·K] - viscosity [kg/m²·s] - density [kg/m³] - attack angle of vortex generator [°] - time interval for naphthalene sublimation [s] - fin thickness, distance between two VGs around the tube [m] - small interval - C distance between the stream direction centerlines of VGs - p pressure drop [Pa] - 0 without VG enhancement - 1, 2, I, II fin surface I, fin surface II, respectively - atm atmosphere - f fluid - fin fin - local local value - m average - naph naphthalene - n,b naphthalene at bulk flow - n,w naphthalene at wall - VG with VG enhancement - w wall or fin surface     

A specific slit die rheometer for extruded starchy products. Design,validation and application to maize starch

A novel in-line rheometer, called Rheopac, has been designed and built in order to study the rheological behaviour of starchy products or, more generally, of products sensitive to a thermomechanical treatment. It is based on the principle of a twin channel, using a balance of feed rate between each of them, in order to make local shear rate vary in the measuring section without changing the flow conditions into the extruder. A wide range of shear rate could be reached and measurements were performed more swiftly than with a classical slit die. The viscous behaviour of maize starch was studied by taking into account the influence of the thermomechanical history, which modified the starch degradation and thus led to important variations in the viscosity. Experimental results were satisfactorily compared to previously published models.Nomenclature E activation energy (J · mol^–1) - h channel depth (m) - h ₁ depth under the piston valve in channel 1 (m) - h ₂ depth under the piston valve in channel 2 (m) - K consistency (Pa·s ⁿ ) - K ₀ reference consistency (Pa·s ⁿ ) - L total channel length (m) - L _p length of the piston valve (m) - MC moisture content (wet basis) - n power law index - N screw rotation speed (rpm) - P ₀ entrance pressure (Pa) - P _e pressure at the entry of the piston valve (Pa) - Q ₁ flow rate in channel 1 (m3 · s^–1) - Q ₂ flow rate in channel 2 m³·s^–1) - Q _T total flow rate (m3 · s^–1) - R constant of perfect gas (8.314 J·mol^–1·K^–1) - SME specific mechanical energy (kWh · t^–1) - T temperature (°C) - T _a absolute temperature (K) - T _b barrel temperature (°C) - T _d die temperature (°C) - T _p product temperature (°C) - w channel width (m) - W energetical term (J·m^–3) - viscosity (Pa · s) - [gh ₀] intrinsic viscosity of native starch (ml·g^–1) - [] intrinsic viscosity (ml·g^–1) - shear rate (s^–1) - shear rate in measuring section (s^–1) - maximum shear rate (s^–1)     

Optimum rib size to enhance forced convection in a horizontal triangular duct with ribbed internal surfaces

The optimum rib size to enhance heat transfer had been proposed through an experimental investigation on the forced convection of a fully developed turbulent flow in an air-cooled horizontal equilateral triangular duct fabricated on its internal surfaces with uniformly spaced square ribs. Five different rib sizes (B) of 5 mm, 6 mm, 7 mm, 7.9 mm and 9 mm, respectively, were used in the present investigation, while the separation (S) between the center lines of two adjacent ribs was kept at a constant of 57 mm. The experimental triangular ducts were of the same axial length (L) of 1050 mm and the same hydraulic diameter (D) of 44 mm. Both the ducts and the ribs were fabricated with duralumin. For every experimental set-up, the entire inner wall of the duct was heated uniformly while the outer wall was thermally insulated. From the experimental results, a maximum average Nusselt number of the triangular duct was observed at the rib size of 7.9 mm (i.e. relative rib size ). Considering the pressure drop along the triangular duct, it was found to increase almost linearly with the rib size. Non-dimensional expressions had been developed for the determination of the average Nusselt number and the average friction factor of the equilateral triangular ducts with ribbed internal surfaces. The developed equations were valid for a wide range of Reynolds numbers of 4,000 < Re _D < 23,000 and relative rib sizes of under steady-state condition. A Inner surface area of the triangular duct [m²] - A _C Cross-sectional area of the triangular duct [m²] - B Side length of the square rib [mm] - C _P Specific heat at constant pressure [kJ·kg^–1·K^–1] - C ₁, C ₂, C ₃ Constant coefficients in Equations (10), (12) and (13), respectively - D Hydraulic diameter of the triangular duct [mm] - Electric power supplied to heat the triangular duct [W] - f Average friction factor - F View factor for thermal radiation from the duct ends to its surroundings - h Average convection heat transfer coefficient at the air/duct interface [W·m^–2 ·K^–1] - k Thermal conductivity of the air [W·m^–1 ·K^–1] - L Axial length of the triangular duct [mm] - Mass flow rate [kg·s^–1] - n ₁, n ₂, n ₃ Power indices in Equations (10), (12) and (13), respectively - Nu _D Average Nusselt number based on hydraulic diameter - P Fluid pressure [Pa] - Pr Prandtl number of the airflow - _c Steady-state forced convection from the triangular duct to the airflow [W] - _l Heat loss from external surfaces of the triangular duct assembly to the surroundings [W] - _r Radiation heat loss from both ends of the triangular duct to the surroundings [W] - Re _D Reynolds number of the airflow based on hydraulic diameter - S Uniform separation between the centre lines of two consecutive ribs [mm] - T Fluid temperature [K] - T _a Mean temperature of the airflow [K] - T _ai Inlet mean temperature of the airflow [K] - T _ao Outlet mean temperature of the airflow [K] - T _s Mean surface temperature of the triangular duct [K] - T Ambient temperature [K] - U Mean air velocity in the triangular duct [m·s^–1] - _r Mean surface-emissivity with respect to thermal radiation - Dynamic viscosity of the fluid [kg·m^–1·s^–1] - Kinematic viscosity of the airflow [m²·s^–1] - Density of the airflow [kg·m^–3] - Stefan-Boltzmann constant [W·m^–2·K^–4]     

Heat transfer and pressure drop for purely viscous non-Newtonian fluids in turbulent flow through rectangular passages

Experimental measurements of friction factor and heat transfer for the turbulent flow of purely viscous non-Newtonian fluids in a 21 rectangular channel are compared with results previously reported for the circular tube geometry. Comparisons are also made with available analytical and empirical predictions.It is found that the rectangular duct fully established friction factor measurements are within ± 5% of the Dodge-Metzner prediction if the Kozicki generalized Reynolds number is used. A modified form of the simpler explicit equation proposed by Yoo, [i.e.f=0.079n ^0.675(Re ^*)^–0.25], is found to yield predictions for both the rectangular duct and the circular tube geometries with approximately the same accuracy as the Dodge-Metzner equation.Fully developed Stanton numbers for the rectangular duct are in good agreement with the circular tube data over a range ofn from 0.37 to 0.88 for a given Prandtl number,Pr _a , when compared at a fixed value of the Reynolds number based on the apparent viscosity evaluated at the wall shear stress. In general, the experimental data are within ± 20% of Yoo's equation,St=0.0152Re _a ^–0.155 Pr _a ^–2/3 . A new equation is proposed to bring the prediction for circular pipes as well as rectangular channels into better agreement with generally accepted Newtonian heat transfer results.

---

Wärmeübergang und Druckverlust für viskose nicht-Newtonsche Fluide in turbulenter Strömung durch rechteckige Kanäle

Zusammenfassung Es werden Messungen des Reibungsfaktors und des Wärmeübergangs bei turbulenter Strömung viskoser nicht-Newtonscher Fluide in einem rechteckigen Kanal mit dem Seitenverhältnis 21 verglichen mit früheren Ergebnissen, die an runden Rohren gewonnen wurden. Weiterhin werden Vergleiche mit aus der Literatur verfügbaren analytischen und empirischen Beziehungen gemacht.Es zeigte sich, daß die Messungen des Reibungsfaktors im rechteckigen Kanal bei vollausgebildeter Strömung auf ± 5% mit der Vorhersage von Dodge-Metzner übereinstimmen, wenn die von Kozicki verallgemeinerte Reynolds-Zahl verwendet wird. Eine modifizierte Form der einfachen von Yoo vorgeschlagenen einfachen Gleichung in explizierter Form (f=0,079n ^0,675(Re ^*)^–0,25) bewies, daß sie sowohl für den rechteckigen Kanal als auch das runde Rohr die Werte mit fast der gleichen Genauigkeit wie die Methode von Dodge-Metzner vorhersagen kann.Die Stanton-Zahlen für den rechteckigen Kanal bei vollausgebildeter Strömung sind in guter Übereinstimmung mit den Werten für das runde Rohr in einem Bereich vonn= 0,37 – 0,88 für eine gegebene Prandtl-Zahl, wenn man den Vergleich bei einem vorgegebenen Wert der Reynolds-Zahl anstellt, die auf die scheinbare Viskosität — abgeleitet aus der Wandschubspannungbezogen ist. Generell läßt sich sagen, daß die Werte auf ± 20% mit der Gleichung von Yoo (St=0,0152Re _a ^–0,155 )Pr _a ^–2/3 ) übereinstimmen. Es wird eine neue Gleichung vorgeschlagen, welche sowohl die Werte für runde Rohre als auch die für rechteckige Kanäle in bessere Übereinstimmung bringt mit den in der Literatur üblichen Ergebnissen für den Wärmeübergang an Newtonsche Fluide.

---

Nomenclature a constant in Eq. (8) - A area of cross-section of channel [m²] - b constant in Eq. (8) - c _p specific heat of test fluid [J kg^–1 K^–1] - d capillary tube diameter [m] - D _h hydraulic diameter, 4A/P[m] - f Fanning friction factor, _w/(g9 V²/2) - h axially local (spanwise averaged) heat transfer coefficient,q _w /(T_wi-T_b) [Wm^–2 K^–1] - k _f thermal conductivity of test fluid [Wm^–1K^–1] - K consistency index of power law fluid - n power law index - Nu fully established, local (spanwise averaged) Nusselt numberh D _h /k _f - P perimeter of channel [m] - Pr _a Prandtl number based on apparent viscosjity, c _p /k _f - Pr ^* defined as (Re _a Pr _a )/Re ^* - q _w wall heat flux [Wm^–2] - Re _a Reynolds number based on apparent viscosity, VD _h/ - Re Metzner's generalized Reynolds number in Eq. (2) - Re ^* Reynolds number defined in Eq. (8) - St Stanton number,h/( V c_p) - T _b local bulk temperature of the fluid [K] - T _wi local inside wall temperature [K] - T _wo local outside wall temperature [K] - V bulk flow velocity [m s^–1] - x distance from the inlet of channel along flow direction [m] Greek symbols shear rate [s^–1] - apparent viscosity [Pa s] - density of test fluid [kg m^–3] - shear stress [Pa] - _w shear stress at the wall [Pa] Dedicated to Prof. Dr.-Ing. U. Grigull's 75th birthday     

Messung der Wärmeleitfähigkeit von Stickstoff und Kohlenmonoxid bei hohen Temperaturen im Stoßwellenrohr

Zusammenfassung Die Meßergebnisse für die Wärmeleitfähigkeit von Stickstoff bei Temperaturen zwischen 1230 und 6000 K und Drückenzwischen 1 und 10 bar und von Kohlenmonoxid zwischen 1150 und 5000 K bei 1 bar werden mitgeteilt. Diese mit dem Stoßwellenrohr gemessenen Werte werden mit jenen verglichen, die sich aus der strengen kinetischen Gastheorie ergeben. Auch verfügbare Daten anderer Autoren werden zum Vergleich herangezogen.

---

Measurement of thermal conductivity of nitrogen and carbon monoxide at high temperatures in a shock tube

The paper presents results of shock-tube measurements of thermal conductivity of nitrogen at temperatures between 1230 and 6000 K and at pressures between 1 and 10 bar and of carbon monoxide at temperatures between 1150 and 5000 K at 1 bar. Experimental results are compared with several variants of theoretical values, computed from rigorous kinetic theory, and with available data of other authors.

---

Bezeichnungen (Einheiten in Klammern) a [m² s^–1] Temperaturleitzahl - C _p[J mol^–1 K^–1] molare Wärmekapazität - k [J K^–1] Boltzmann-Konstante - M [kg mol^–1] molare Masse - p bar Gesamtdruck - R [J mol^–1 K^–1] Gaskonstante - T [K] thermodynamische Temperatur - t [s] Zeit - U [J mol^–1] innere Energie - w [m s^–1] Geschwindigkeit - x [m] Ortskoordinate - x _i [1] Molanteil der Komponentei im Gasgemisch - [Wm^–1 K^–1] Wärmeleitfähigkeit - [mol m^–3] molare Konzentration Indizes i die Komponentei im Gasgemisch - g bezieht sich auf das (kalte) Gas bei der Wandtemperatur - w bezieht sich auf die feste Wand - p bei konstantem Druck Dieser Beitrag wurde auf dem Thermodynamik-Kolloquium des VDI im Oktober 1969 in Zürich vorgetragen.     

A new transient method for analysis of solidification in the continuous casting of metals

Solidification processes involve complex heat and mass transfer phenomena, the modelling of which requires state-of-the art numerical techniques. An efficient and accurate transient numerical method is proposed for the analysis of phase change problems. This method combines both the enthalpy and the enhanced specific heat approaches in incorporating the effects of latent heat released due to phase change. The sensitivity and accuracy of the proposed method to both temporal and spatial discretization is shown together with closed-form solutions and the results from the enhanced specific heat approach. In order to explore the proposed method fully, a non-linear heat release, as is the case for binary alloys, is also examined. The number of operations required for the new transient approach is less than or equal to the enhanced heat capacity method depending on the averaging method adopted. To demonstrate the potential of this new finite-element technique, measurements obtained on operating machines for the casting of zinc, aluminum and steel are compared with the model predictions. The death/birth technique, together with the proper heat-transfer coefficients, were employed in order to model the casting process with minimal error due to the modelling itself.Nomenclature [A] conductance matrix - [B] matrix containing the derivative of the element shape functions - c, C _p specific heat (J kg^–1°C^–1) - effective specific heat (J kg^–1°C^–1) - f(T) local liquid fraction - f thermal load vector - H enthalpy (J kg^–1) - [H] capacitance matrix - h, h _r,h _c heat transfer coefficient (W m^–2°C^–1) - K thermal conductivity (W m^–1°C^–1) - L latent heat of solidification (J kg^–1) - l overall length (m) - N _i shape functions - Q rate of heat generation per unit volume (J m^–3) - q heat flux (W m^–2) - R residual temperature (°C) - T temperature (°C) - T _s solidus temperature (°C) - T _l liquidus temperature (°C) - T _pouring pouring temperature (°C) - T _top temperature at the top of the mould (°C) - T _w temperature of the water spray (°C) - approximated temperature (°C) - T surrounding temperature (°C) - cooling rate (°C/s) - t time (seconds) - x _i,x, y, z spatial variables (m) - t time step (s) - x element size (m) - diffusivity (m²s^–1) - density (kg m^–3) - time marching parameter - _n direction cosines of the unit outward normal to the boundary     

Supersonic viscous perfect gas flow past a circular cylinder

The complete Navier-Stokes equations are used to calculate supersonic perfect gas flow past a circular isothermal cylinder by the method described in [1]. The effects of the Mach number M=2.5–10 and the Reynolds number Re=30-10⁵ on the flowfield structure and heat transfer to the cylinder wall are investigated. Special attention is paid to the study of the near wake and the local characteristics on the leeward side of the cylinder.Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No.6, pp. 107–115, November–December, 1993.     

Shear stress and normal stress measurements of aqueous polymer solutions in a cone-and-plate rheogoniometer

Summary The rheological behaviour of aqueous solutions of Separan AP-30 and Polyox WSR-301 in a concentration range of 10–10000 wppm is investigated by means of a cone-and-plate rheogoniometer. The relation between the shear stress and the shear rate is for lower shear rates characterized by a timet ₀, which is concentration dependent. Both polymers show for 4000 s^–1 < < 10000 s^–1 a behaviour similar to that of a Bingham material, characterized by a dynamic viscosity ₀ and an apparent yield stress ₀, which also depend on the concentration. The inertial forces are measured for water and some other Newtonian liquids. An explanation is given why the theoretical model developed for these forces does not match the experimental values; the shape of the liquid surface is shear rate dependent. To obtain the first normal stress difference, we have to correct for these inertial forces, the surface tension and the buoyancy. The normal forces, measured for Separan AP-30, appear to be a linear function of the shear rate for 350 s^–1 < < 3300 s^–1.

---

Zusammenfassung Das rheologische Verhalten wäßriger Polymerlösungen von Separan AP-30 und Polyox WSR-301 wird in einem Konzentrationsgebiet von 10–10000 wppm in einem Kegel-Platte-Rheogoniometer untersucht. Der Zusammenhang zwischen Schubspannung und Schergeschwindigkeit wird für niedrige Schergeschwindigkeiten durch eine konzentrationsabhängige Zeitt ₀ gekennzeichnet. Für Schergeschwindigkeiten 4000 s^–1 < < 10000 s^–1 zeigen beide Polymere ein genähert binghamsches Verhalten, gekennzeichnet durch eine dynamische Viskosität ₀ und eine scheinbare Fließgrenze ₀, welche ebenfalls konzentrationsabhängig sind. Die Trägheitskräfte werden für Wasser und einige newtonsche Öle bestimmt. Die Abweichung der experimentellen Ergebnisse vom theoretischen Modell wird durch die Abhängigkeit der Gestalt der Flüssigkeitsoberfläche von der Schergeschwindigkeit erklärt. Um die Werte der ersten Normalspannungsdifferenz zu erhalten, muß man bezüglich der Trägheitskräfte, der Oberflächenspannung und der Auftriebskräfte korrigieren. Die Normalspannungen für Separan AP-30, gemessen für 350 s^–1 < < 3300 s^–1, zeigen eine lineare Abhängigkeit von der Schergeschwindigkeit.

---

c concentration (wppm) - g acceleration of gravity (ms^–2) - K force (N) - K _b buoyant force (N) - K _c force, acting on the cone (N) - K ₀ dimensional constant def. by eq. [24] (N) - K _s force, def. by eq. [22] (N) - M dimensional constant def. by eq. [24] (Ns) - P _s pressure def. by eq. [17] (Nm^–2) - P ₀ average pressure in the liquid atr = 0 (Nm^–2) - P _R average pressure in the liquid atr = R (Nm^–2) - r ₁,r ₂ radii of curved liquid surface (m) - R platen radius (m) - R _w radius of wetted platen area (m) - S _x standard deviation ofx - t ₀ characteristic time def. by eq. [1] (s) - T temperature (°C) - V volume of the submerged part of the cone (m³) - v tangential velocity of liquid (ms^–1) - x distance (m) - angle (rad) - ₀ cone angle (rad) - calibration constant (Nm^–3) - shear rate (s^–1) - dynamic viscosity (mPa · s) - ₀ viscosity def. by eq. [1] (mPa · s) - contact angle (rad) - density (kgm^–3) - static surface tension (Nm^–1) - shear stress (Nm^–2) - ₀ yield stress def. by eq. [1] (Nm^–2) - _c, _p angular velocity (c = cone,p = plate) (s^–1) With 8 figures and 3 tables     

Viscous heating and heat transfer in cone-and-plate rheogoniometry

Summary At higher shear rates the relation between shear stress and shear rate appears to deviate from the for Newtonian fluids expected linear behaviour. In cone-and-plate rheogoniometry one of the most important causes of that is the effect of viscous heating. Accurate measurements carried out with a 10 cm diameter cone and plate lead to a semi-logarithmic, linear relationship between temperature increase and time for a Newtonian oil which dynamic viscosity varies approximately linearly with time. A simple model based on a heat balance describes this behaviour quantitatively.

---

Zusammenfassung Bei newtonschen Flüssigkeiten weisen die Experimente eine Abweichung vom linearen Zusammenhang zwischen Schubspannung und Schergeschwindigkeit auf. Im Kegel-Platte-Meßsystem ist die Wärmeproduktion durch innere Reibung die wichtigste Ursache der Abweichung. Bei newtonschen Flüssigkeiten, deren dynamische Viskosität sich ungefähr linear mit der Temperatur verändert, ergeben sorgfältig ausgeführte Messungen mit einem Kegel von 10 cm Durchmesser einen linearen Zusammenhang zwischen der Zeit und dem Logarithmus der Temperaturzunahme. Ein aus der Wärmebilanz abgeleitetes Modell vermag dieses Verhalten quantitativ zu beschreiben.

---

Symbols A platen surface (m²) - B viscosity constant from eq. [1] (Pa s K^–1) - S _B standard deviation ofB (Pa s K^–1) - S _t0 standard deviation oft ₀ (s) - S _t0 standard deviation oft ₀ (s) - S ₀ standard deviation of ₀ (Pa s) - t time (s) - t ₀ time def. by eq. [5] (s) - t ₀ time def. by eq. [11] (s) - T temperature (°C) - T ₀ temperature of the surrounding air (°C) - T highest experimental temperature (°C) - V volume of the fluid between the platen (m³) - W heat capacity of the system (J K^–1) - heat transfer coefficient (W m^–2 K^–1) - shear rate (s^–1) - dynamic viscosity (Pa s) - ₀ dynamic viscosity atT ₀ (Pa s) - dimensionless temperature def. by eq. [4a] (–) - dimensionless time def. by eq. [4b] (–) - dimensionless time def. by eq. [10] (–) With 4 figures and 2 tables     

Calculation of a wall-stabilized argon arc with account for radiative energy transfer

It was shown experimentally in [1, 2] and in a study by E. I. Asinovskii and A. V. Kirillin reported at the Scientific Technical Conference of the High-Temperature Scientific Research Institute held in 1964 that the heat transfer mechanism in a wall-stabilized argon arc was not purely purely conductive at gas temperatures greater than 11 000° K. Asinovskii and Kirillin also showed that radiative energy transfer is the reason why similarity is lost when the current-voltage characteristics are constructed in dimensionless form. The radiation of an argon arc has been studied experimentally by a number of authors [3–5], Dautov [6] calculated an argon arc without allowing for radiation.In this article an argon arc stabilized by the cooled duct walls is calculated with account for radiation using theoretically computed relationships describing the transport properties of argon plasma. A large portion of the radiated energy pertains to spectral lines whose role has been studied by L. M. Biberman. The authors have used I. T. Yakubov's data on argon radiation published in the journal Optics and Spectroscopy. A method of calculation and data on argon plasma radiation are also given in [7].Reference [8] deals with the problem of the role of radiation in an arc burning in nitrogen. In particular, the above-mentioned loss of similarity follows from the results of this work. However, the relationships used in this article to describe the transport properties of nitrogen plasma were obtained experimentally in [9].Notation r₀ arc radius (cm) - r variablesradius (cm) - T temperature (°K) - heat transfer coefficient (ergcm^–1sec^–1deg^–1) - E electric field intensity (g^1/2cm^–1/2sec^–1) - electrical conductivity (sec^–1) - q₁ heat flux density due to conduction - q₂ heat flux density due to radiation - u divergence of radiative energy flux density in the transparent part of the spectrum (ergcm^–3sec^–1) - u₂ same for part of the spectrum where reabsorption is taken taken into account - m₀ atomic mass - m_e electronic mass - Stefan-Boltzmann constant - h Planck constant - k Boltzmann constant - e electronic charge - p pressure - degree of ionization - N_e electron concentration (cm^–3) - n₀ neutral atom concentration - Q_0e electron-neutral collision cross section - Q_ie electron-ion collision cross section (cm²) - ₀ line center frequency (sec^–1) - ₊ line halfwidth (distance from line center to contour for ) due to effects giving dispersion contour - k _v absorption coefficient (cm^–1) - energy radiated by a hemispherical volume - emissivity of hemispherical volume - radius of hemispherical volume - S line intensity The authorS thank I. T. Yakubov for allowing them to use his data on arc plasma radiation.     

Rheologische Messungen an Kunststoff-Schmelzen mit einem neuen Rotationsrheometer

Zusammenfassung In der vorliegenden Arbeit wird ein neues Rotationsrheometer vorgestellt und über Messungen an zwei Polymethylmethacrylat-Formmassen berichtet. Bei dem Rheometer handelt es sich um ein Couette-Rheometer mit feststehendem Innenzylinder als Meßkörper. Der Meßkörper ist beidseitig eingespannt. In dem geschlossenen Meßraum können die Schmelzen bis zu einem Druck von 500 bar belastet werden.Der zeitliche Verlauf der Schubspannung in den Schmelzen wird in Abhängigkeit von Temperatur und Druck aufgezeichnet.

---

Summary A new type of rotational rheometer is described, and results for two samples of polymethylmethacrylate are reported. The rheometer consists of a Couette system with fixed inner cylinder, supported at both ends for torque measurements. Pressure may be varied up to 500 bar. Shear stresses have been recorded as a function of time, temperature and pressure.

---

Nomenklatur C [kp cm^–2 s^–1] Steigung der Anlaufkurve im Nullpunkt - D [kp cm rad^–1] Direktionsmoment - E ₀ [kcal mol^–1] Aktivierungsenergie der Newtonschen Viskosität - G [kp cm^–2] Schubmodul - G [—] Griffith-Zahl - l [mm] Länge des Meßkörpers - p [kp cm^–2] Druck - R _i [mm] Radius des Innenzylinders - R _a [mm] Radius des Außenzylinders - t _max [s] Zeit, bei der das Maximum in der Anlaufkurve auftritt - T [°C] Temperatur - ₀ [cm² kp^–1] Druckkoeffizient der Newtonschen Viskosität - [s^–1] Schergeschwindigkeit - ₀ [kp s cm^–2] Newtonsche Viskosität - (g cm²] Trägheitsmoment des Meßkörpers - v ₀ [s^–1] Eigenfrequenz des Meßsystems - _max [kp cm^–2] maximale Schubspannung - _st [kp cm^–2] stationäre Schubspannung Mit 7 Abbildungen und 1 Tabelle     

On the evolutions of the discontinuities of orders ⩾ m in solutions for quasi-linear systems of PDEs of order m,or for the linearizations of these systems

Summary One considers a system L[u]=0 of PDEs, quasi-linear (according to [1]) and of order m, which possesses a bicharacteristic line , as it happens in the hyperbolic case. For v=0, , –m (>0) let u(^v) be a discontinuity wave of order m+v that solves the system above and whose discontinuity hypersurface includes . The corresponding transport equations along are considered. Furthermore some interesting cases are pointed out, in which these equations turn out to be mutually equivalent in a suitable sense. Some theorems are stated to compare the transport equations for the discontinuities of the above kinds, that are connected with the systems d^hL[u]/dt^h=0 (h=0, , –m) and/or the linearization of the system L[u]=0 around any regular solution of it.

---

Sommario Si considera un sistema L[u]=0 di equazioni alle derivate parziali, quasi lineare (secondo [1]) e di ordine m, il quale sia dotato di qualche bicaratteristica , come accade nel caso iperbolico. Per v=0, , –m(>0) sia u(^v) un'onda di discontinuità di ordine m+v risolvente il detto sistema e avente ipersuperficie di discontinuità contenente Si considerano le relative equazioni di trasporto lungo e si determinano casi interessanti in cui queste equazioni sono mutuamente equivalenti in senso opportuno. Si stabiliscono teoremi di confronto per il trasporto delle discontinuità del tipo suddetto, relative ai sistemi d^hL[u]/dt^h=0 (h=0, , –m) e/o alla linearizazione del sistema L[u]=0 attorno a qualche sua soluzione regolare.

     

Dynamic crack resistance of acrylic plastic

The dependence of the stressintensity factor near the tip of a growing crack in an SO120 acrylic plastic on the crackpropagation velocity K_I(:v) within the range of 10^–5 –300 m/sec is determined. Specific features of crack propagation associated with the shape of the curve K_I(v), which has discontinuities and nonuniqueness intervals, are discussed.     

Flows in the initial sections of channels with permeable walls

Calculations of two types of flows in the initial sections of channels with permeable walls are carried out on the basis of semiempirical turbulence theories during fluid injection only through the walls and during interaction of the external flow with the injected fluid. Experimental studies of the first type [1–3] show that at least within the limits of the lengths L/h<30 and L/a< 50 (2h is the distance between permeable walls of a flat channel anda is the tube radius) the velocity distributions in the laminar and turbulent flow regimes differ little and are nearly self-similar for solutions obtained in [4]. For sufficiently large Reynolds numbers, Re₀>100 (Re₀=v₀h/ or Re₀=v₀ a/, where v₀ is the injection velocity), and small fluid compressibility, the axial velocity component is described by the relations for ideal eddying motion: u=(/2)x× cos (y/2) in a flat channel and u=x cos (y²/2) in atube (the characteristic values for the coordinates are, respectively, h anda). Measurements indicate the existence of a segment of laminar flow; its length depends on the Reynolds number of the injection [3]. In the turbulent regime the maximum generation of turbulent energy occurs significantly farther from the wall than in parallel flow. Flows of the second type in tubes were studied in [5–7]. These studies disclosed that for Reynolds numbers of the flow at the entrance to the porous part of the tube Re=u₀ a/<3.10³ fluid injection with v₀/u₀>0.01 leads to suppression of turbu lence in the initial section of the tube. An analogous phenomenon was observed in the boundary layer with v₀/u₀>0.023 [8, 9]. Laminar-turbulent transition in flows with injection was explained in [10, 11] on the basis of hydrodynamic instability theory, taking into account the non-parallel character of these flows. The mechanisms for the development of turbulence and reverse transition in channels with permeable walls are not theoretically explained. Simple semiempirical turbulence theories apparently are insufficient for this purpose. In the present work results are given of calculations with two-parameter turbulence models proposed in [12, 13] for describing complex flows. Due to the sharp changes of turbulent energy along the channel length, a numerical solution of the complete system of equations of motion was carried out by the finite-difference method [14].Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 43–48, September–October, 1976.