Скачать книгу

interesting multi‐scale approach for amine‐based absorbers was proposed by Raynal and Royon‐Lebeaud [6], and reviewed by Raynal et al. [7], who divided the modeling in the three scales illustrated in Figure 2.2. First off, two‐dimensional (2‐D) simulations were used at the small scale to represent the gravity‐driven liquid flow falling down the corrugations of the packing by using an interface tracking numerical technique named Volume‐of‐Fluid (VOF) method. The latter introduces an additional conservation equation for a colour function f, the volume fraction, which varies within the closed interval from 0 through 1. Cells within the mesh with values at both ends of that range correspond to either phase, whereas values in between correspond to the interface. The VOF method has proved a useful tool for the description of a number of free surface interface problems across many engineering fields. It ideally requires however the use of adaptive grids or in the case of fixed grids at least a greater resolution in those cells that are expected to contain the interface [8]. This requirement comes motivated by the discontinuous character of gas–liquid interfaces in Nature, whereas an unnatural transition area between phases appears instead of a sharp interface when using the VOF algorithm. The results can therefore be affected by substantial errors if the mesh is not treated conveniently in the vicinity of the interface. Moreover, the VOF method is characterized by the assumption that the density and viscosity assigned to the interface cells are volume fraction averaged. There is therefore some variability in the density along the computational domain that results in the appearance of additional mass and momentum source terms in the conservation equations. The implications of these additional source terms are yet to be assessed and are minimized when extra mesh refinement is introduced [9]. The average film thickness and the gas–liquid interface velocity were the direct output from the calculations at the small scale in Raynal and Royon‐Lebeaud's work. Then, the authors estimated indirectly the total liquid hold‐up as the product between the average liquid film thickness from the simulation and the specific surface of the packing given by the manufacturer.

Schematic illustration of the flow of information between the three scales proposed in the modeling strategy presented by Raynal and Royon-Lebeaud.

      Source: Adapted from Raynal and Royon‐Lebeaud [6].

      At the intermediate scale, one representative elementary unit (REU, the repeating unit that forms the geometry of the packing) was represented because the pressure drop per unit length in a limited set of REUs matches that of the entire column [10]. Raynal and Royon‐Lebeaud [6] proposed a single‐phase flow approach (with only gas phase) at the REU level, although in their case, information from the small scale would be fed into their REU set‐up in order to account for the effect of the presence of liquid. To do so, two modifications at the REU scale setup were introduced. On the one hand, instead of a non‐slip boundary condition on the packing walls, i.e. velocity equal zero, the authors established a fixed velocity (the gas–liquid interface velocity obtained at the small scale). On the other hand, and in order to compare their data with a set of experiments including gas and liquid flow, the authors corrected the F‐factor by the liquid hold‐up obtained from the average liquid film thickness of the small‐scale simulations. The F‐factor F is a measure of the kinetic energy of the gas phase that enters the packed bed and comes determined by the product between the superficial velocity and the square root of the gas‐phase density. The corrected value of the F‐factor F′ is obtained by dividing it over the volume occupied by the gas phase per unit volume of packing, that is, the complementary value of the liquid hold‐up. A greater F‐factor then, than that of a dry simulation, is obtained for the same value of the gas superficial velocity. This is in accordance with the narrowing effect that the presence of the liquid phase has on the channels through which the gas flows, and therefore, greater gas pressure drops are obtained upon wet conditions. These intermediate scale calculations allowed the authors to obtain the pressure drop coefficients K in the three coordinate directions as the ratio per unit length between the pressure drop ΔP and the superficial dynamic pressure of the gas phase.

Скачать книгу