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.     

Internal loops in pseudoelastic behaviour of Ti-Ni shape memory alloys: Experiment and modelling

The stress-strain isothermal hysteresis loops due to the incomplete martensitic transformation are analysed for Ti-Ni shape memory alloys. Experiments show the existence of two distinct yield lines for phase transition; one for the forward transformation austenitemartensite (AM), the other for the reverse transformation MA. The tensile behaviour of single crystals with only one yield line (AM) [1] can be considered as an ideal case. An extension of a thermodynamic model for pseudoelasticity [2] allows these two yield lines to be taken into account.

---

Sommario Per leghe Ti-Ni con memoria di forma vengono analizzati i cicli di isteresi isotermici tensione-deformazione prodotti da una incompleta trasformazione martensitica. Gli esperimenti mostrano l'esistenza di due distinte linee di snervamento per la transizione di fase, una verso la trasformazione austenitemartensite (AM), l'altra per la trasformazione inversa MA. Il comportamento a trazione di un singolo cristallo con una sola linea di snervamento (AM) [1], può essere considerato un caso ideale. L'estensione ad un modello termodinamico pseudo-elastico [2] consente di analizzare queste due linee di snervamento.

     

Transient non-Darcy free convection between parallel vertical plates in a fluid saturated porous medium

Transient non-Darcy free convection between two parallel vertical plates in a fluid saturated porous medium is investigated using the generalized momentum equation proposed by Vafai and Tien. The effects of porous inertia and solid boundary are considered in addition to the Darcy flow resistance. Exact solutions are found for the asymptotic states at small and large times. The large time solutions reveal that the velocity profiles are rather sensitive to the Darcy number Da when Da<1. It has also been found that boundary friction alters the velocity distribution near the wall, considerably. Finite difference calculations have also been carried out to investigate the transient behaviour at the intermediate times in which no similarity solutions are possible. This analytical and numerical study reveals that the transient free convection between the parallel plates may well be described by matching the two distinct asymptotic solutions obtained at small and large times.Nomenclature C empirical constant for the Forchheimer term - f velocity function for the small time solution - F velocity function for the large time solution - g acceleration due to gravity - Gr* micro-scale Grashof number - H a half distance between two infinite plates - K permeability - Nu Nusselt number - Pr Prandtl number - t time - T temperature - u, v Darcian velocity components - x, y Cartesian coordinates - effective thermal diffusivity - coefficient of thermal expansion - porosity - dimensionless time - similarity variable - dimensionless temperature - viscosity - kinematic viscosity - density - the ratio of heat capacities     

Propagation of a spherical wave in nonlinearly compressible and elastoplastic materials

The problem of spherical wave propagation in soil under the action of an intense uniformly decreasing load ₀(t) applied to the boundary of a cavity with radius r₀ is considered. Soil with a high stress level is modeled either by ideally nonlinearly compressible or elastoplastic material, taking account of linear irreversible unloading for the material. In contrast to [1–7], in order to describe material movement use is made of strain theory [8] with determining functions = (), _i=_i(_i), where , _i, , _i are the first and second invariants of strain and stress tensors. During material loading these functions are presented in the form of polynomials ()=(_i+₂¦¦), _ii)=(_i-_2i)_i, in which constant coefficients _i, _i=1, 2) are determined by experiment, taking account of the triaxial stressed state of soil. Solution of the problem is constructed by an analytically reversible method, with prescribed shape for the shock-wave (SW) surface in the form of a second-degree polynomial relating to time t and a numerical method of characteristics for a prescribed arbitrarily decreasing load _i(t). On the basis of the analytical equations obtained, calculations are carried out for material parameters (including loading profile) in a computer and stresses and mass velocity of plastic and elastoplastic materials are compared.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 4, pp. 95–100, July–August, 1986.The authors express their sincere thanks to Kh. A. Rakhmatulin for discussing the results of this work.     

Heat transfer in magnetohydrodynamic flow near a stagnation point

Summary Steady, axisymmetric, magnetohydrodynamic flow with a stagnation point on an infinite plane wall is considered with a magnetic field applied normal to the wall. Solutions are obtained in the form of series for the velocity, magnetic field and temperature when the magnetic field parameter () and the ratio of viscosity to magnetic diffusivity () are small. The case=O(1) is considered briefly when solutions which Meyer³⁾ obtained by physical order-of-magnitude arguments are derived mathematically as expansions in. Some remarks are made on the consistency of extending the results to flow within the boundary layer near the nose of a bluff body.     

Scheme of Investigation of the Spectrum of a Family of Perturbed Operators and Its Application to Spectral Problems in Thick Junctions

We develop a scheme for the investigation of the asymptotic behavior of eigenvalues and eigenvectors of a family of self-adjoint compact operators {A: > 0} that act in different spaces and lose their compactness in the limit case 0. We prove the Hausdorff convergence of the spectrum of the operator A to the spectrum of the limit operator A₀, obtain asymptotic estimates for this convergence both to points of the discrete spectrum and to points of the essential spectrum of the operator A₀, and prove asymptotic estimates for eigenvectors of A. This scheme is applied to the investigation of the asymptotic behavior of eigenvalues and eigenfunctions of the Neumann problem in a thick singularly degenerate junction that consists of two domains connected by an -periodic system of thin rods of fixed length.     

Modèles intégro-différentiels isotropes en viscoélasticité linéaire

Summary In this work, we consider isotropic viscoelastic materials characterized by a linear relation between stress and strain and their time derivatives up to the orderm,n. We observe thatColeman's thermodynamics does not imply the causality of the model but that a stronger requirement ofMeixner does always imply it. The thermodynamical restrictions ofMeixner are stated on the basis of theorems for rational positive functions. The propagation of sinusoidal plane waves in infinite medium is studied whenm 2. In the casesm = 0,n = 2 andm = 1,n = 3, it is shown that, when the dissipation is weak, the logarithmic decrement of the waves is smaller than 2 only in a finite frequency band. In the casem = 2,n = 4, we got < 2 only in disjoined finite frequency bands.

---

Zusammenfassung In dieser Arbeit befassen wir uns mit isotropen viskoelastischen Stoffen, die durch eine lineare Beziehung zwischen der Spannung und der Dehnung sowie deren Zeitableitungen bis zur Ordnungm bzw.n charakterisiert sind. Wir stellen fest, daßColemans Thermodynamik nicht die Kausalität des Modells impliziert, die stärkere Forderung vonMeixner dieses hingegen leistet. Die thermodynamischen Einschränkungen vonMeixner werden auf der Basis von Theoremen für rationale positive Funktionen dargestellt. Die Fortpflanzung von sinusförmigen ebenen Wellen in einem unendlichen Medium wird fürm 2 untersucht. Für die Fällem = 0,n = 2 undm = 1,n = 3 wird gezeigt, daß, solange die Dissipation gering ist, das logartithmische Dekrement der Wellen nur in einem endlichen Frequenzband kleiner als 2 ist. Für den Fallm = 2,n = 4 erhält man < 2 nur in nicht überlappenden endlichen Frequenzbändern.

     

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

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

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

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

Deformation to breaking of deep water gravity waves

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

Application of the Chapman-Enskog method to the case of a binary two-temperature gas mixture

An interesting property of the flows of a binary mixture of neutral gases for which the molecular mass ratio =m/M1 is that within the limits of the applicability of continuum mechanics the components of the mixture may have different temperatures. The process of establishing the Maxwellian equilibrium state in such a mixture divides into several stages, which are characterized by relaxation times _i which differ in order of magnitude. First the state of the light component reaches equilibrium, then the heavy component, after which equilibrium between the components is established [1]. In the simplest case the relaxation times differ from one another by a factor of ^*.Here the mixture component temperature difference relaxation time _T /, where is the relaxation time for the light component. If 1, 1, so that _T ~1, then for the characteristic hydrodynamic time scale t~1 the relative temperature difference will be of order unity. In the absence of strong external force fields the component velocity difference is negligibly small, since its relaxation time _vt1.In the case of a fully ionized plasma the Chapman-Enskog method is quite easily extended to the case of the two-temperature mixture [3], since the Landau collision integral is used, which decomposes directly with respect to . In the Boltzmann cross collision integral, the quantity appears in the formulas relating the velocities before and after collision, which hinders the decomposition of this integral with respect to , which is necessary for calculating the relaxation terms in the equations for temperatures differing from zero in the Euler approximation [4] (the transport coefficients are calculated considerably more simply, since for their determination it is sufficient to account for only the first (Lorentzian [5]) terms of the decomposition of the cross collision integrals with respect to ). This led to the use in [4] for obtaining the equations of the considered continuum mixture of a specially constructed model kinetic equation (of the Bhatnagar-Krook type) which has an undetermined degree of accuracy.In the following we use the Boltzmann equations to obtain the equations of motion of a two-temperature binary gas mixture in an approximation analogous to that of Navier-Stokes (for convenience we shall term this approximation the Navier-Stokes approximation) to determine the transport coefficients and the relaxation terms of the equations for the temperatures. The equations in the Burnett approximation, and so on, may be obtained similarly, although this derivation is not useful in practice.     

Problem concerning a spherical piston in a compressible medium with “dry” friction

We study the self-similar problem concerning the motion of a spherical piston in a medium with dry friction. The piston moves with constant velocity in a nonideal medium.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 1, pp. 136–140, January–February, 1973.The authors thanks T. F. Kryukov for performing the major part of the calculations.     

On an inversion theorem for conical regions in elasticity theory

Summary Two-dimensional stress singularities in wedges have already drawn attention since a long time. An inverse square-root stress singularity (in a 360° wedge) plays an important role in fracture mechanics.Recently some similar three-dimensional singularities in conical regions have been investigated, from which one may be also important in fracture mechanics.Spherical coordinates are r, , . The conical region occupied by the elastic homogeneous body (and possible anisotropic) has its vertex at r=0. The mantle of the cone is described by an arbitrary function f(, )=0. The displacement components be u. For special values of (eigenvalues) there exist states of displacements (eigenstates) % MathType!MTEF!2!1!+-% feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqef0uAJj3BZ9Mz0bYu% H52CGmvzYLMzaerbd9wDYLwzYbItLDharqqr1ngBPrgifHhDYfgasa% acOqpw0xe9v8qqaqFD0xXdHaVhbbf9v8qqaqFr0xc9pk0xbba9q8Wq% Ffea0-yr0RYxir-Jbba9q8aq0-yq-He9q8qqQ8frFve9Fve9Ff0dme% GabaqaaiGacaGaamqadaabaeaafiaakabbaaa6daaahjxzL5gapeqa% aiaadwhadaWgaaWcbaGaeqOVdGhabeaakiabg2da9iaadkhadaahaa% WcbeqaaiabeU7aSbaakiaadAgadaWgaaWcbaGaeqOVdGhabeaakiaa% cIcacqaH7oaBcaGGSaGaeqiUdeNaaiilaiabfA6agjaacMcaaaa!582B!\[u_\xi = r^\lambda f_\xi (\lambda ,\theta ,\Phi )\],which may satisfy rather arbitrary homogeneous boundary conditions along the generators.The paper brings a theorem which expresses that if is an eigenvalue, then also-1- is an eigenvalue. Though the theorem is related to a known theorem in Potential Theory (Kelvin's theorem), the proof has to be given along quite another line.

---

Zusammenfassung Zwei-dimensionale Spannungssingularitäten in keilförmigen Gebieten sind schon längere Zeit untersucht worden und neuerdings auch ähnliche drei-dimensionale Singularitäten in konischen Gebieten.Kugelkoordinaten sind r, , . Das konische Gebiet hat seine Spitze in r=0. Der Mantel des Kegels lässt sich beschreiben mittels einer willkürlichen Funktion f(, )=0. Die Verschiebungskomponenten seien u. Für spezielle Werte von (Eigenwerte) bestehen Verschiebunszustände % MathType!MTEF!2!1!+-% feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqef0uAJj3BZ9Mz0bYu% H52CGmvzYLMzaerbd9wDYLwzYbItLDharqqr1ngBPrgifHhDYfgasa% acOqpw0xe9v8qqaqFD0xXdHaVhbbf9v8qqaqFr0xc9pk0xbba9q8Wq% Ffea0-yr0RYxir-Jbba9q8aq0-yq-He9q8qqQ8frFve9Fve9Ff0dme% GabaqaaiGacaGaamqadaabaeaafiaakabbaaa6daaahjxzL5gapeqa% aiaadwhadaWgaaWcbaGaeqOVdGhabeaakiabg2da9iaadkhadaahaa% WcbeqaaiabeU7aSbaakiaadAgadaWgaaWcbaGaeqOVdGhabeaakiaa% cIcacqaH7oaBcaGGSaGaeqiUdeNaaiilaiabfA6agjaacMcaaaa!582B!\[u_\xi = r^\lambda f_\xi (\lambda ,\theta ,\Phi )\],welche homogene Randwerte der Beschreibenden des Kegels entlang genügen.Das Bericht bringt ein Theorem, welches aussagt, das und =–1– beide Eigenwerte sind.

     

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

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

Some self-similar solutions of the Leibenzon equation

Self-similar one-dimensional solutions of the Leibenzon equation c²_t= _zz ^k (z 0, k 2) are considered. Approximate solutions are constructed for the two cases in which the initial value = ₁ = const > 0 and on the boundary either a constant value = ₂ < ₁ is maintained or the flow (directed outwards) is given. In the first problem the dependence of the boundary flow on the governing parameters is determined. A characteristic property of the types of motion in question is the existence near the boundary of a region, expanding with time, in which the flow is almost independent of the coordinate.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 145–150, September–October, 1991.     

Analytical Modeling of Nonaqueous Phase Liquid Dissolution with Green's Functions

Equilibrium and bicontinuum nonequilibrium formulations of the advection–dispersion equation (ADE) have been widely used to describe subsurface solute transport. The Green's Function Method (GFM) is particularly attractive to solve the ADE because of its flexibility to deal with arbitrary initial and boundary conditions, and its relative simplicity to formulate solutions for multidimensional problems. The Green's functions that are presented can be used for a wide range of problems involving equilibrium and nonequilibrium transport in semiinfinite and infinite media. The GFM is applied to analytically model multidimensional transport from persistent solute sources typical of nonaqueous phase liquids (NAPLs). Specific solutions are derived for transport from a rectangular source (parallel to the flow direction) of persistent contamination using first, second, or thirdtype boundary or source input conditions. Away from the source, the first and thirdtype condition cannot be expected to represent the exact surface condition. The secondtype condition has the disadvantage that the diffusive flux from the source needs to be specified a priori. Near the source, the thirdtype condition appears most suitable to model NAPL dissolution into the medium. The solute flux from the pool, and hence the concentration in the medium, depends strongly on the mass transfer coefficient. For all conditions, the concentration profiles indicate that nonequilibrium conditions tend to reduce the maximum solute concentration and the total amount of solute that enters the porous medium from the source. On the other hand, during nonequilibrium transport the solute may spread over a larger area of the medium compared to equilibrium transport.     

Spin-up in a saturated porous medium

It is shown that, on the Brinkman model, spin-up is confined to boundary layers whose thickness is of order k ^1/2, and the spin-up is established in a time of order k/, where k, , and denote permeability, density, porosity and dynamic viscosity, respectively.     

Messung und Berechnung des örtlichen und mittleren Stoffübergangs an stumpf angeströmten Kreisscheiben bei unterschiedlicher Turbulenz

Zusammenfassung In einer vergleichenden Literaturübersicht zu Umströmung, Druck- bzw. Geschwindigkeitsverteilung sowie Wärme- und Stoffübergang werden bislang vorliegende Angaben zu stumpf angeströmten Kreisscheiben und -Zylindern zusammengefaßt. Wenige und zudem divergierende Ergebnisse zum Wärme- und Stoffübergang machen grundlegende experimentelle und theoretische Untersuchungen notwendig, wie sie in [l, 2] für die Eichung von Stoffübergangsmeßmethoden benötigt werden.Unter Einbeziehung des quer angeströmten Kreiszylinders wird gezeigt, daß genaue Angaben zum Wärme- und Stoffübergang bei zwei- wie dreidimensionalen Staupunktströmungen bislang nur über die Messung möglich sind. Über gemessene Geschwindigkeitsverteilungen berechnete Stoffübergangskoeffizienten werden von der Messung nicht bestätigt. Sie liegen gegenüber dem Experiment zu niedrig.Die Messungen wurden bei Turbulenzintensiten 0,8%Tu6%, Reynolds-Zahlen 2·10³5 und Scheibendurchmessern 9,3mmd73,7mm durchgeführt. Der Einfluß der Turbulenz auf den Stoffübergang im Staupunkt von Kreisscheiben kann nur näherungsweise über den Smith-Kuethe-Parameter Tu · Re/100 erfaßt werden. Differenzen zwischen Theorie nach Smith und Kuethe für Tu· Re<5 und Messung lassen sich über die Stabilitätstheorie erklären. Für eine genauere Erfassung des Stoffübergangs muß den unterschiedlichen Transportvorgängen über Turbulenzballen oder Längswirbeln sowie der Struktur der Turbulenz Rechnung getragen werden.

---

Measuring and computation of local and average mass transfer to disks in cross flow at different turbulence intensities

The results of different publications concerning the flow, pressure and velocity distributions as well as the heat and mass transfer of disks and cylinders in cross flow are compared by a literature review. A few diverging results for heat and mass transfer require new experimental and theoretical approaches. The calibration of recently developed techniques for the determination of mass transfer rates as published in [1, 2] make these investigations expecially necessary. Including the cylinder in cross flow the authors show, that up to now exact data of heat and mass transfer for two- or three-dimensional flow at a forward stagnation region can be obtained by direct measuring only.Mass transfer coefficients computed from measured velocity distributions are not confirmed by the experimental results. Compared to the experimental data they are too low. The measurements were accomplished for turbulence intensities 0.8%Tu6%, Reynolds-numbers 2· 10³5 and disk diameters 9.3 mm d 73.7 mm.The influence of the turbulence on the stagnation point mass transfer of disks can be obtained only approximately by the Smith-Kuethe-parameter Tu·Re/100. Differences between theoretical results of Smith and Kuethe and experimental ones for Tu·Re/100<5 may be explained by the stability theory. For a more accurate determination of the mass transfer the different transport mechanisms of the scale of turbulence or the tree-dimensional flow pattern like Taylor-Görtler-vortices as well as the structure of the turbulence itself have to be regarded.

---

Bezeichnungen a Temperaturleitkoeffizient - C_p Beiwert für den statischen Druck - C₂, C₃ Gradient der bezogenen Geschwindigkeit U⁺ am Staupunkt bei ebener, räumlicher Strömung - D_A Diffusionskoeffizient von Ammoniak in Luft - d Durchmesser - Fr=Sh/Re Frössling-Zahl für den Stoffübergang - Fr=Nu/Re Frössling-Zahl für den Wärmeübergang - Le=a/D_A Lewis-Zahl - L Bezugslänge - M Maschenweite von Turbulenzgittern - Nu=·d/ Nußbelt-Zahl - n Exponent der Prandtl-bzw. Schmidt-Zahl - Pr=/a Prandtl-Zahl - p Druck, Partialdruck - p_x statischer Druck an der Stelle x am Rand der Grenzschicht - Re=U · d/ Reynolds-Zahl - r Radius - r(x) radiale Distanz von der Rotationsachse eines Körpers zu einem Oberflächenelement - Sc=/D_A Schmidt-Zahl - Sh= _A ·d/D_A Sherwood-Zahl - T absolute Temperatur - Tu Turbulenzintensität (Turbulenzgrad) in% - U Strömungsgeschwindigkeit in x-Richtung am Rand der Grenzschicht - U Hauptströmungsgeschwindigkeit im freien Kanalquerschnitt - U⁺=U/U bezogene Geschwindigkeit in x-Richtung am Rand der Grenzschicht - u Strömungsgeschwindigkeit in x-Richtung, tangential zur Oberfläche - mittlere turbulente Geschwindig-keitsschwankung in x-Richtung - v Strömungsgeschwindigkeit in y-Richtung, normal zur Oberfläche - x Koordinate in Strömungsrichtung, tangential zur Oberfläche - x_G Entfernung vom Turbulenzgitter in Strömungsrichtung - x⁺ bezogene Länge x/r - y Koordinate normal zur Oberfläche - Wärmeübergangskoeffizient - _A Stoffübergangskoeffizient (Ammoniak) - dimensionsloses Temperaturgefälle an der Wand - Keilvariable - Wärmeleitkoeffizient - Wirbelweilenlänge (mm) - kinematische Zähigkeit - transformierte bezogene Länge - _A Partialdichte von Ammoniak Indices B mit Korrektur aufgrund der Verengung - m mittel - S bezogen auf die Kreisscheibe - Z bezogen auf den Kreiszylinder Herrn Prof. Dr.-Ing. habil. Josef Ipfelkofer zum 70. Geburtstag am 7. April 1977 gewidmet.