Near Wall Condition

The first approach developed by Hsu (1962) is widely used to determine ONE in conventional size channels and in micro-channels (Sato and Matsumura 1964 Davis and Anderson 1966 Celata et al. 1997 Qu and Mudawar 2002 Ghiaasiaan and Chedester 2002 Li and Cheng 2004 Liu et al. 2005). These models consider the behavior of a single bubble by solving the one-dimensional heat conduction equation with constant wall temperature as a boundary condition. The temperature distribution inside the surrounding liquid is the same as in the undisturbed near-wall flow, and the temperature of the embryo tip corresponds to the saturation temperature in the bubble 7s,b- The vapor temperature in the bubble can be determined from the Young-Laplace equation and the Clausius-Clapeyron equation (assuming a spherical bubble) ... 

In the absence of diffusion, all hydrodynamic models show infinite variances. This is a consequence of the zero-slip condition of hydrodynamics that forces Vz = 0 at the walls of a vessel. In real systems, molecular diffusion will ultimately remove molecules from the stagnant regions near walls. For real systems, W t) will asymptotically approach an exponential distribution and will have finite moments of all orders. However, molecular diffusivities are low for liquids, and may be large indeed. This fact suggests the general inappropriateness of using to characterize the residence time distribution in a laminar flow system. Turbulent flow is less of a problem due to eddy diffusion that typically results in an exponentially decreasing tail at fairly low multiples of the mean residence time. 

Experiments were conducted in a large (-26 m3) radon/thoron test facility (RTTF) designed for calibration purposes and simulation studies (Bigu, 1984). A number of different materials were exposed in the RTTF to a radon/radon progeny or thoron/thoron progeny atmosphere. Exposure of the materials was carried out under laboratory-controlled conditions of radiation level, aerosol concentration, air moisture content and temperature. The materials used were in the form of circular discs of the same thickness (-0.5 mm) and diameter (-25 mm), and they were placed at different locations on the walls of the RTTF at about 1.6 m above the floor. Other samples were placed on horizontal trays. Samples (discs) of different materials were arranged in sets of 3 to 4 they were placed very close to one another to ensure exposure under identical conditions. Exposure time was at least 24 hours to ensure surface activity equilibrium, or near equilibrium, conditions. 

The standard wall function is of limited applicability, being restricted to cases of near-wall turbulence in local equilibrium. Especially the constant shear stress and the local equilibrium assumptions restrict the universality of the standard wall functions. The local equilibrium assumption states that the turbulence kinetic energy production and dissipation are equal in the wall-bounded control volumes. In cases where there is a strong pressure gradient near the wall (increased shear stress) or the flow does not satisfy the local equilibrium condition an alternate model, the nonequilibrium model, is recommended (Kim and Choudhury, 1995). In the nonequilibrium wall function the heat transfer procedure remains exactly the same, but the mean velocity is made more sensitive to pressure gradient effects. 

The Presumed Probability Density Function method is developed and implemented to study turbulent flame stabilization and combustion control in subsonic combustors with flame holders. The method considers turbulence-chemistry interaction, multiple thermo-chemical variables, variable pressure, near-wall effects, and provides the efficient research tool for studying flame stabilization and blow-off in practical ramjet burners. Nonreflecting multidimensional boundary conditions at open boundaries are derived, and implemented into the current research. The boundary conditions provide transparency to acoustic waves generated in bluff-body stabilized combustion zones, thus avoiding numerically induced oscillations and instabilities. It is shown that predicted flow patterns in a combustor are essentially affected by the boundary conditions. The derived nonreflecting boundary conditions provide the solutions corresponding to experimental findings. 

In general, large industrial fixed beds operate under near-adiabatic conditions, whereas small laboratory-scale fixed beds may approach isothermal operation (Ruthven, 1984). Especially, for most environmental applications, for catalytic, adsorption, and ion-exchange operations, the species to be removed are in such low concentrations that the operarion is nearly isothermal. Thus, the heat transfer to the external fixed-bed wall is often of minimal importance. 

Hot-Wall Reactors. Because of the large mass diffusivities and nearly isothermal conditions (except for the entrance zone) in hot-wall, low-pressure reactors (50 Pa), multicomponent diffusion and chemical reactions are critical... 

The reaction is also influenced by the heat of reaction that develops during the conversion of the reactants, a problem in tubular screening reactors. In micro structures, the heat transport through the walls of the channels is facilitated owing to their small dimensions. The catalysts are deposited on the walls of these micro structures and will thus have the appropriate environment for exothermic reactions by enabling fast quenching of the reaction with near isothermal conditions. Hence also the heat and mass balance in the reactor will be decoupled, which permits the analytical description of the flow in the screening reactor. 

Analysts should review the technical basis for uncertainties in the measurements. They should develop judgments for the uncertainties based on the plant experience and statistical interpretation of plant measurements. The most difficult aspect of establishing the measurement errors is establishing that the measurements are representative of what they purport to be. Internal reactor CSTR conditions are rarely the same as the effluent flow. Thermocouples in catalyst beds may be representative of near-wall instead of bulk conditions. Heat leakage around thermowells results in lower than actual temperature measurements. 

Aqueous solutions vaporize with nearly the same coefficient as pure water if attention is given to boiling-point elevation and if tbe solution does not become saturated and care is taken to avoid dry wall conditions. 

Since there is no radial bulk transport of fluid between the monolith channels, each channel acts basically as a separate reactor. This may be a disadvantage for exothermic reactions. The radial heat transfer occurs only by conduction through the solid walls. Ceramic monoliths are operated at nearly adiabatic conditions due to their low thermal conductivities. However, in gas-liquid reactions, due to the high heat capacity of the liquid, an external heat exchanger will be sufficient to control the reactor temperature. Also, metallic monoliths with high heal conduction in the solid material can exhibit higher radial heat transfer. 

In most high Reynolds number flows, the wall function approach gives reasonable results without excessive demands on computational resources. It is especially useful for modeling turbulent flows in complex industrial reactors. This approach is, however, inadequate in situations where low Reynolds number effects are pervasive and the hypotheses underlying the wall functions are not valid. Such situations require the application of a low Reynolds number model to resolve near-wall flows. For the low Reynolds number version of k-s models, the following boundary conditions are used at the walls ... 

At the outlet, extrapolation of the velocity to the boundary (zero gradient at the outlet boundary) can usually be used. At impermeable walls, the normal velocity is set to zero. The wall shear stress is then included in the source terms. In the case of turbulent flows, wall functions are used near walls instead of resolving gradients near the wall (refer to the discussion in Chapter 3). Careful linearization of source terms arising due to these wall functions is necessary for efficient numerical implementation. Other boundary conditions such as symmetry, periodic or cyclic can be implemented by combining the formulations discussed in Chapter 2 with the ideas of finite volume method discussed here. More details on numerical implementation of boundary conditions may be found in Patankar (1980) and Versteeg and Malalasekara (1995). 

Both equilibrium and non-equilibrium wall boundary implementations are considered. For equilibrium flows the local production rate of turbulence equals the dissipation rate in the near wall grid node. The first set of wall function boundary conditions reported was apparently used for equilibrium flows by Gosman et al. [59]. Denoting the dependent variables in the first point near the wall by a subscript P, an approximate sketch of their approach is given next. 

Let us consider the case of an arbitrary viscoplastic fluid with yield stress To (similar results for nonlinear viscous fluids correspond to to = 0). To obtain the temperature profile, we proceed as follows. First, in the near-wall shear region 0 < < h- ho, where ho = toL/AP, we solve Eq. (6.5.2) with the boundary conditions (6.5.1). Then in the quasisolid region h - ho 2 < h, we solve Eq. (6.5.2) with - 0 under the boundary condition (6.5.3). Finally, we match the two solutions on the common boundary = ho. This procedure results in the following temperature distribution in the channel ... 

As seen from the Eq. (9), accomplishment of the value (9 o)max requires both increasing the near-wall zone permeability Kf, and decreasing the radius of pores involved in the vaporization process. It is known that bidisperse pore stmetures provide maximum heat flux densities at the highest efficiency [25]. Decrease of circumferential channel dimensions (a and b) also causes increasing of ( o)max, however, it also escalating the entire hydraulic resistance inside the LHP. In certain conditions, it plays the major role in heat flow Q and heat flux qo limitations before the value ( o)max is attained. 

---