Effect of particulate thermophoresis in reducing the fouling rate advantages of effusion-cooling

To predict small particle diffusional mass transfer (deposition), including particle thermophoresis, transpiration cooling and variable properties, the coupled ordinary differential equations governing self-similar laminar boundary layers are solved numerically. Under typical combustion turbine conditions, although diffusional deposition rates can be dramatically reduced by transpiration cooling (eg by some 5-decades for mainstream submicron particles corresponding to a Schmidt number of about 10² (or d_p ≈ 0.7 × 10^?2μm) and a wall transpiration-cooled to $\text{T}{}_{\mathrm{<mtext>w</mtext>}}\text{T}{}_{\mathrm{<mtext>e</mtext>}}\text{= 0.8}$) actual deposition rate reductions will be smaller than previously expected (by about 1 decade for particles with Sc ≈ 10²), owing to thermophoretic particle drift ‘caused’ by the colder wall. Such micro-droplets, small enough to behave like ‘heavy molecules’ in combustion systems, are often important because they can cause adherence of the much larger, supermicron, ash particles which inertially impact on the same surface     

Boundary layer in compressible flow of wet water vapour

In this paper the influence of small droplets, with radius $\text{10}{}^{\mathrm{?8}}\text{m}\text{< r < 10}{}^{\mathrm{?6}}\text{m}$, on laminar and turbulent boundary layer behavior is considered. It is found that the laminar boundary layer in a two-phase flow with strongly dispersed liquid retains dissipation energy and that the recovery factor of enthalpy is greater than unity. In turbulent boundary layers small droplets are transported by turbulent diffusion and this leads to the recovery factor being less than unity. Its value in both cases depends mainly on the nondimensional number $\text{Ds}\text{= C}{}_{\mathrm{Le}}\text{L/(U}{}_{\mathrm{e}}{}^2\text{/2)}$. The laminar boundary layer solution for non-equilibrium two-phase flow is obtained. Profiles of the droplet mass fraction, vapour and droplets temperatures and droplet radius are computed for the case of a steady two-dimensional flow. The turbulent boundary layer is treated using a semi-empirical theory assuming thermodynamic equilibrium.     

Correlation of turbulent flow rate-pressure drop data for non-newtonian solutions and slurries in pipes

Isothermal and non-isothermal flow rate-pressure drop data in turbulent flow through smooth pipes have been obtained for non-Newtonian fluids, including aqueous solutions of polymers and aqueous suspensions of titanium dioxide. It has been found that the friction factor, f, is a function of a new form of Reynolds number, Re_B, based on the parameters A, x and w of Bowen's correlation, viz.

$\text{τ}{}_{\mathrm{w}}\text{D}{}^{\mathrm{x}}\text{=A}\text{u}{}^{\mathrm{w}}$

where τ_w is the wall shear strees, ?u the mean velocity, D the pipe diameter; A, x and w are experimentally derived parameters which characterise the fluid.     

A simple stochastic model of two phase flow pressure drop accompanied by boiling inside circular tubes with in-line static mixers

A simple stochastic model has been developed for boiling pressure drop inside a circular tube with in-line static mixers. This model gives rise to the dimensionless correlating equation of the form:

$\text{[f] = a[Pr}{}_{\mathrm{m}}\text{]}{}^{\mathrm{<mtext>?1</mtext><mtext>2</mtext>}}\text{[Re}{}_{\mathrm{L}}\text{]}{}^{\mathrm{?1}}\text{μm}\text{μ}{}_{\mathrm{L}}\text{ρ}\text{ρ}{}^{\mathrm{<mtext>?1</mtext><mtext>2</mtext>}}\text{H}{}_{\mathrm{se}}\text{+ xH}{}_{\mathrm{LG}}\text{C}{}_{\mathrm{pm}}\text{ΔT}{}^{\mathrm{<mtext>?1</mtext><mtext>2</mtext>}}$

This correlation shows good agreement with the experimental data.     

Non-linear boundary and eigenvalue problems for the emden-fowler equations by group methods

Two pragmatic boundary value and eigenvalue problems of the Emden-Fowler equation $\text{(t}{}^{\mathrm{\alpha }}\text{u′)′ + λt}{}^{\mathrm{\beta }}\text{?(u) = 0,?(u) = u}{}^{\mathrm{\gamma }}$ and e^u are studied using the simple one parameter group properties. In all cases boundary value problems are converted into initial value problems using the property of the invariance group. With $\text{?(u) = u}{}^{\mathrm{\gamma }}$ an eigenvalue problem is detailed and calculations presented.     

Investigation of droplet deposition from a turbulent gas stream

Turbulent deposition of particles from two-phase flow onto the smooth wall of a tube has been studied theoretically and experimentally. A model is proposed for the deposition motion of large particles based on turbulent diffusion in the core followed by a free flight towards the wall. The theory shows that within the Stokes regime, the dimensionless deposition velocity $\text{k-}{}_{\mathrm{<mtext>d</mtext>}}\text{/}\text{u}$* depends on Re and $\text{τ}{}^{\mathrm{+}}$ only, where u* is the friction velocity, Re is the tube Reynolds number and $\text{τ}{}^{\mathrm{+}}$ is the dimensionless particle relaxation time. Deposition data are obtained for air-water droplet flow through a 12.7-mm i.d. acrylic tubing at $\text{Re}\text{= 52,500}$ and 94,600. The proposed theory satisfactorily describes the existing deposition data as well as present measurements, covering a wide range of Re and $\text{τ}{}^{\mathrm{+}}$.     

Bifurcations et solutions multiples en cavité 3D différentiellement chaufféeBifurcations and multiple solutions in a differentially heated cubic cavity

In this paper, a differentially heated square/cubic cavity is studied by performing three-dimensional direct numerical simulations. The first bifurcation observed at $\text{R}{}_{\mathrm{a}}\text{≈3.2×10}{}^7$ is due to the 3D vortex structures generated at the end regions of vertical boundary layers near the median plane. The main results of this Note are that the flow returns to a steady state for higher values of the Rayleigh number $\text{R}{}_{\mathrm{a}}$ (7×10⁷ and 10⁸ for example) still exhibiting these 3D vortex structures, and that multiple steady flows which differ by their symmetry properties, are obtained for $\text{R}{}_{\mathrm{a}}\text{=10}{}^8$. However, the flow reverts to unsteadiness for $\text{R}{}_{\mathrm{a}}\text{=3×10}{}^8$. In this latter case, the instability is due to the vertical boundary layers. To cite this article: G. de Gassowski et al., C. R. Mecanique 331 (2003).     

Time-average local thickness measurement in falling liquid film flow

A method for monitoring time-varying local film thickness variation through the detection of laser scattering from suspended latex particles is briefly described. This method was used in conjunction with the Jeffreys theory of drainage from a flat plate to determine time-average local film thickness.Measurements were made at Reynolds numbers (equal to ($\text{4Q/ν}$)) from 145 to 4030 at varying distances along the direction of flow. At the bottom of the flow, 134 cm from the top, average film thickness is given by the expression: where $\text{a}{}_{\mathrm{i}}$ and $\text{n}{}_{\mathrm{i}}$ are constants unique to each of the three Reynolds number regions, wavy laminar, transitional and turbulent.     

Homogenization limit for electrical conduction in biological tissues in the radio-frequency rangeLimite d'homogénéisation pour la conduction électrique dans les tissus biologiques dans le domaine des radiofréquences

We study an evolutive model for electrical conduction in biological tissues, where the conductive intra-cellular and extracellular spaces are separated by insulating cell membranes. The mathematical scheme is an elliptic problem, with dynamical boundary conditions on the cell membranes. The problem is set in a finely mixed periodic medium. We show that the homogenization limit u₀ of the electric potential, obtained as the period of the microscopic structure approaches zero, solves the equation $\text{?}\text{div}\text{(σ}{}_0\text{?}{}_{\mathrm{x}}\text{u}{}_0\text{+A}{}^0\text{?}{}_{\mathrm{x}}\text{u}{}_0\text{+∫}{}_0{}^{\mathrm{t}}\text{A}{}^1\text{(t?τ)?}{}_{\mathrm{x}}\text{u}{}_0\text{(x,τ)}\text{d}\text{τ?}\text{F}\text{(x,t))=0}$ where σ₀>0 and the matrices A⁰, A¹ depend on geometric and material properties, while the vector function $\text{F}$ keeps trace of the initial data of the original problem. Memory effects explicitly appear here, making this elliptic equation of non standard type. To cite this article: M. Amar et al., C. R. Mecanique 331 (2003).     

Linear methods in problems of non-linear differential equations with a small parameter

Methods have been considered for deriving asymptotical formulas for the systems of the type

$\text{ε}{}^{\mathrm{p}}\text{d}\text{x}{}_{\mathrm{k}}\text{d}\text{t}\text{= f}{}_{\mathrm{k}}\text{(x) + ε}\text{f}\text{?}{}_{\mathrm{k}}\text{(x) + …}$

by constructing an analog of the Schrödinger perturbation theory of the linear operator

$\text{∑}\text{k}\text{[f}{}_{\mathrm{k}}\text{(x) + ε}\text{f}\text{?}{}_{\mathrm{k}}\text{(x)]}\text{?F}\text{?x}{}_{\mathrm{k}}\text{= A}{}_{\mathrm{o}}\text{F + εA}{}_1\text{F.}$

These methods can be extended to some classes of partial differential equations, in particular, to Whitham's non-linear theory.     

Homogenized model of reaction–diffusion in a porous mediumUn modèle homogénéisé de réaction–diffusion dans un milieu poreux

We study the initial boundary value problem for the reaction–diffusion equation,

$\text{?}{}_{\mathrm{t}}\text{u}{}^{\mathrm{\epsilon }}\text{??·(a}{}^{\mathrm{\epsilon }}\text{?u}{}^{\mathrm{\epsilon }}\text{)+g(u}{}^{\mathrm{\epsilon }}\text{)=h}{}^{\mathrm{\epsilon }}$

in a bounded domain $\text{Ω}$ with periodic microstructure $\text{F}{}^{\mathrm{(\epsilon )}}\text{∪}\text{M}{}^{\mathrm{(\epsilon )}}$, where a^ε(x) is of order 1 in $\text{F}{}^{\mathrm{(\epsilon )}}$ and κ(ε) in $\text{M}{}^{\mathrm{(\epsilon )}}$ with κ(ε)→0 as ε→0. Combining the method of two-scale convergence and the variational homogenization we obtain effective models which depend on the parameter θ=lim_ε→0κ(ε)/ε². In the case of strictly positive finite θ the effective problem is nonlocal in time that corresponds to the memory effect. To cite this article: L. Pankratov et al., C. R. Mecanique 331 (2003).     

An experiment to measure the pressure dependence of the zero-shear-rate viscosity

A simple experimental method, based on Stokes' law for falling spheres, has been devised and used to measure the pressure-dependence of the zero-shear-rate viscosity of a polypropylene melt. The experiment was performed by maintaining three thick-walled test cylinders containing the polymer melt and the falling sphere at the same elevated temperature but different pressures for periods of time ranging from 20 to 48 hours.When compared with experiments using high-pressure capillary or rotational viscometers, this experimental method has the advantages that viscous heating is non-existent and the apparatus and data analysis are relatively simple. The principal disadvantage encountered here, thermal degradation at high temperatures, could probably be reduced by molding specimens under vacuum and by shortening the exposure time. Since the falling-sphere experiment provides data at very low shear rates and the capillary and rotational viscometers generate data at high shear rates, the two experimental methods are complementary.The pressure coefficient b [=d(In η₀/dp] was determined for Hercules Pro-fax 6523 polypropylene in two series of experiments at different temperatures. For seven experiments at 218.3°C and pressures up to 97.9 $\text{MN}\text{m}{}^2$ (14,200 psi), the average value of b ± 95% confidence limits was found to be 14.8 ($\text{GN}\text{m}{}^2$)^?1 ± 2.9.The average b was 12.6 ($\text{GN}\text{m}{}^2$)^?1 ± 1.4 in a series of eight experiments at 232.2°C and pressures up to 123 $\text{MN}\text{m}{}^2$ (17,800 psi).     

The Einstein viscosity correction in n dimensions

The effective viscosity of a dilute suspension of rigid $\text{n-}\text{dimensional}$ hyperspheres in a viscous fluid at small particle Reynolds numbers is determined; the result being

$\text{μ}{}_{\mathrm{eff}}\text{=μ}\text{1+}\text{n+2}\text{2}\text{φ}$

. Expressions are also given for the $\text{n-}\text{dimensional}$ Stokes velocity and pressure fields for a hypersphere in a pure straining flow.     

Two-phase flow in a vertical pipe and the phenomenon of choking: Homogeneous diffusion model—I. Homogeneous flow models

The paper examines the topological structure of all possible solutions which can exist in flows through adiabatic constant-area ducts for which the homogeneous diffusion model has been assumed. The conservation equations are one-dimensional with the single space variable $\text{z.}$ but gravity effects are included. The conservation equations are coupled with three equations of state: a pure substance, a perfect gas with constant specific heats, and a homogeneous two-phase system in thermodynamic equilibrium. The preferred state variables are pressure $\text{P.}$ enthalpy $\text{h.}$ and mass flux $\text{G}{}^2$.The three conservation equations are first-order but nonlinear. They induce a family of solutions which are interpreted as curves in a four-dimensional phase space conceived as a union of three-dimensional spaces ($\text{P, h, G}{}^2\text{, z}$) with $\text{G}{}^2\text{=}\text{const}$ treated as a parameter. It is shown that all points in these spaces are regular, so that no singular solutions need to be considered. The existence and uniqueness theorem leads to the conclusion that through every point in phase space there passes one and only one solution-curve.The set of differential equations, treated as a system of algebraic equations of each point of the phase space, determines the components of a rate-of-change vector which are obtained explicitly by Cramer's rule. This vector is tangent to the solution curve. Each solution curve turns downward in $\text{z}$ at some specific elevation $\text{z}{}^1$, and this determines the condition for choking. Choking occurs always when the exit flow velocity at $\text{L = z}{}^1$ is equal to the local velocity of propagation of small plane disturbances of sufficiently large wavelength, that is when the flow rate $\text{G}$ becomes equal to a specified, critical flow rate, $\text{G}{}^1$. (The possible dependence of the sonic velocity on frequency in a real flow is ignored, because it has not been allowed for in the equations of the model under study.) A criterion, analogous to the Mach number, which indicates the presence or absence of choking in a cross section is the ratio $\text{K = G/G}{}^7$ of the mass-flow rate $\text{G}$ to the local critical mass flow rate. $\text{G}{}^7$, $\text{K = 1}$ denoting choking. The critical parameters depend only on the thermodynamic properties of the fluid and are independent of the gravitational acceleration and shearing stress at the wall.The topological characteristics of the solutions allow us to study all flow patterns which can, and which cannot, occur in a pipe of given length $\text{L}$ into which fluid is discharged through a rounded entrance from a stagnation reservoir and whose back-pressure is slowly lowered. The set of flow patterns is analogous to that which occurs with a perfect gas, except that the characteristic numerical values are different. They must be obtained by numerical integration and the influence of gravity must be allowed for.The preceding conclusions are valid for all assumptions concerning the shearing stress at the wall which make if dependent on the state parameters only, but not on their derivatives with respect to $\text{z}$. However, the study is limited to upward flows for which the shearing stress at the wall and the gravitational acceleration are codirectional.     

The effect of hydrodynamic coupling on the axisymmetric drainage of thin films

The effect of hydrodynamic coupling of adjacent phases on the axisymmetric drainage of thin films is examined using a prototype model of coalescence. For long times, pressure forces in the film dominate flow in all three regions, and finally all move effectively as one, whereas for short times, profiles are sharp and initial flow differences in the three regions can dominate pressure effects. For intermediate times, temporal evolution of velocity profiles depends in a complicated way on the kinematic viscosity ratio and the parameter $\text{R = (?}{}_{\mathrm{A}}\text{μ}{}_{\mathrm{A}}\text{/?}{}_{\mathrm{B}}\text{/μ}{}_{\mathrm{B}}\text{)}{}^{\mathrm{<mtext>1</mtext><mtext>2</mtext>}}$, as well as on initial conditions and pressure gradient. Generally speaking, the initial flows have less of an effect on overall drainage time than the presence of induced circulation in adjacent phases. Analytical solutions are plotted for a range of systems and representative initial conditions and pressure gradients. In a subsequent article, film-thinning equations are solved using this information.     

Dispersion in annular climbing film flow

The dispersion of a tracer injected as a pulse into a climbing liquid film is investigated for a series of water and air flow rates, and for a number of different electrolyte tracers. It is found that at all flow rates the observed concentration distribution depends on the nature of the tracer. This observation is explained in terms of two effects: molecular diffusion in a viscous sub-layer and ion fractionation associated with droplet formation at the gas-liquid interface. The overall dispersive characteristics of the system are described in terms of a mathematical model assuming dispersed plug flow in both the film and entrained droplets with interchange between these phases. This model is fitted to experimental tracer concentration distributions using a non-linear least-squares regression procedure. The parameter values obtained from the fitting procedure are studied to determine trends with flow rates and tracer properties. Values for a film dispersion parameter, $\text{P}{}_{\mathrm{f}}$, are found to correlate significantly with the molecular diffusion coefficients of the tracers. Consistent values for an ion fractionation coefficient, $\text{k}{}_{\mathrm{if}}$, are also obtained.