Investigation of carbon dioxide (CO2) capture in a falling film contactor by computer simulation

In this work, mathematical models of carbon dioxide (CO 2 ) absorption by monoethanolamine amine (MEA) in a falling film contactor are developed. The proposed models aim to predict conversion of the gas–liquid reaction along the contactor, gas–liquid interface temperature profile (axial and radial), liquid film thickness along the contactor length, axial and radial concentration profiles of reactants in liquid film, and axial and radial profiles of velocity in the liquid film. A code written in MatLab was used to obtain these profiles based on multi grid method through programming of kinetic and thermodynamic equations and physical properties of the absorption system. The mathematical model is validated by an experimental measurement based on absorption of CO 2 gas by MEA solution. Four parameters are studied as independent variables namely, mole fraction of carbon dioxide in gaseous mixture, molar concentration of absorbent (MEA, volumetric flow rate of MEA, and its temperature. It is found that the entrance effect of the falling film contactor is related to axial distance from the contactor entrance exponentially:           E = B 0 exp(− B 1 y )An optimization technique based on minimization of the sum of the squared error between the experimental and predicted composition of absorption process is used to obtain B 0 and B 1 . It is found that reaction between carbon dioxide and MEA is instantaneous, and the axial conversion of carbon dioxide in the gas phase varies exponentially with the contactor length.


Introduction
Removal of carbon dioxide has been practiced industrially for several decades.Different processes need to have CO 2 removed.In natural gas processing, CO 2 is removed to reduce the costs of compression, transportation and corrosion.In ammonia manufacture, CO 2 needs to be removed from the hydrogen stream, since it poisons the catalyst for the reaction between H 2 and N 2 .Power plant flue gases are a new application of CO 2 removal processes, compared to the first two.In this case, CO 2 is removed only to reduce greenhouse emissions.This issue is of increasing interest, because global warming is an important environmental and political issue.With the Kyoto protocol of 1997, forty-one industrialized countries agreed to cut the carbon dioxide emissions to approximately 5% less than the emissions in 1990, in a five year period going from 2008 to 2012 [1].This situation makes CO 2 capture an important issue in the economies of most countries.It is a fact that a CO 2 removal facility is an expensive plant, and the operating costs are very high.For a power plant the removal of 90% of the CO 2 from the flue gas can use up to 30% of the energy produced by the plant (IEA Greenhouse Gas R&D Programme) [2].This represents a heavy loss, and the focus of research in this field is to reduce the costs of removal, in particular reducing its energy requirements.CO 2 has some uses and can be sold on the market.It is used mostly in Enhanced Oil Recovery (EOR).CO 2 flooding represents one of the main methods for extracting the final amounts of recoverable oil from depleted reservoirs.CO 2 is also used in the food industry for carbonated beverages [2].It finds uses in smaller quantities as inert gas or as supercritical solvent.The problem of the CO 2 market is that many processes produce CO 2 as a by-product, and there is no need to buy from flue gas plants.In addition the amount of CO 2 needed for the aforementioned processes is much smaller than the amount that can be recovered from power plant flue gases.CO 2 capture is typically done by absorption with alkanolamine-water solutions although other processes were also considered such as catalytic hydrogenation [3] and utilization in renewable energy technologies [4].The alkanolamines are bases, and they react with the acid species CO 2 to form different reaction products.They contain alcohol groups in order to become soluble in water.The first alkanolamine to be used industrially was monoethanolamine (MEA).This is the only amine that will be considered in this work, the chemical structure is shown in Figure 1 compared to the structures of other common alkanolamines, such as methyldiethanolamine (MDEA), diethanolamine (DEA), triethanolamine (TEA), diglycolamine (DGA), and piperazine (PZ).The amines are distinguished as primary, secondary and tertiary, according to the number of organic groups attached to the alkaline nitrogens.MEA is a primary amine, DEA and PZ are secondary amines, and MDEA and TEA are tertiary amines.The different categories of amines differ on the type of mechanisms with which they react with CO 2 , as well as the reaction products and the heats of reaction.Typically primary and secondary amines react forming a carbamate species, and the reaction may or may not proceed through an intermediate called the Zwitterions [5].
Tertiary amines cannot form a carbamate species, because they do not have a hydrogen attached to the nitrogen atom.Typically the tertiary amines react according to equation (3).
Primary and secondary amines usually react fast er than tertiary amines, and CO 2 has higher heats of absorption in these amines.Heats of reaction at 25 °C and unloaded conditions are approximately 20.3 kcal/mole for MEA and 14.8 kcal/mole for MDEA [6].For MEA, the equation is particularly simple.The rate is dominated by the zwitterion formation, rather than by its reaction with another base.The advantage of a fast reacting amine is that the size of the separation equipment (absorbers) is smaller.The drawback is that more heat is required to reverse the reaction and regenerate the solvent.Solvent regeneration is the main obstacle to cost reduction of CO 2 capture.The choice of a proper solvent is important.Some work has been done on mixed solvents.Bishnoi [7] and Dang [8] researched the properties of piperazine promoted MDEA and MEA respectively.Cullinane [9] studied the absorption into piperazine promoted potassium carbonate (K 2 CO 3 ).Piperazine has very fast kinetics and can promote the absorption rates.MDEA, MEA and K 2 CO 3 provide CO 2 capacity in the solvent and reduce the heat of regeneration, which would be otherwise too high if piperazine was to be used by itself.These mixed solvents are being studied at bench scale, but they are not yet used in industry.Sharifi and Omidbakhsh Amiri [10] studied the effect of tower type on gas sweetening process.They found that different internals give different efficiency of gas capturing.The first industrial application of the falling film principles was in the form of falling film evaporator [11].The surprising characteristics of this type of evaporator made the investigators of the field increasingly start to study mass and reaction mechanism through the falling film device, since most absorption and reaction mechanisms are accompanying with a heat transfer process.The falling film device offers greater benefits than the conventional one for many industrial processes, namely [12]: -It can be operated at low temperature difference between heating medium and liquor.-Operation is stable at high turndown ratios.
-It can be cleaned more effectively and more efficiently.
-It is more suitable for computer controlled facilities.
-It is more economical to operate and require less operation attention.
Figure 2 shows a sectional view of the falling film contactor, a liquid of the reactant falls by gravity, completely wetting the solid wall and contacting the vapor which is diluted in an inert gas generally air.The solid surface is refrigerated by extremely circulating water.Due to the high surface to volume ratio of liquid in the column, efficient heat elimination takes place.Additional temperature control is achieved by diluting the reactive gas with air or with nitrogen.In some falling film reactors, the control of the film temperature is of considerable importance since by product formation can occur at higher temperature [13].The engineer in the process industries is usually concerned with the operation of existing plants and the development of new process.In the first case the control, improvement, Fig. 1.Chemical structures of most common alkanolamines [1].and optimization of the operation are the engineer's main objectives.In order to achieve this, a quantitative representation of the process, a model is needed which would give the relationship between the various parts of the system.The objective of modeling is to construct, from theoretical and empirical knowledge of a process, a mathematical formulation which can be used to predict the behavior of this process.Thus, mathematical modeling receives high attention in such falling film processes and other sweetening towers [14] to control the process temperature.In order to achieve this goal for CO 2 -MEA system, a model is needed to obtain the relationship between the different operating and design parameters of the system for the purpose of industrial design.To validate the mathematical model, experimental measurements of CO 2 absorption by MEA are conducted.

Mathematical modeling
The mathematical model consists of differential balance equations in the liquid phase with their corresponding boundary conditions.Through these equations, velocity, concentration, and temperature profiles are obtained.The equations are based on steady state operation, common in industrial and laboratory reactors.

The principle assumptions
Modeling of the liquid phase is based on the following assumptions: -The liquid circulates in a laminar flow, and the gas circulates co-currently in a turbulent flow.-The liquid film is symmetric with respect to the contactor axis.
-The film thickness is small compared to the column radius.
-Liquid reactant and liquid product are assumed to be nonvolatile at working temperatures.-The solubility of the reactant gas in the liquid reactant and in the reaction product is ideal according to Henry's law.

Equations of the mathematical model
The column is divided into N number of annular segments as shown in Figure 3. Thus, there are N number of equations for the liquid reactant and the same number of equations is presented for the dissolved gas.To formulate a comprehensive model, the effects of hydrodynamics, mass transfer, kinetics, and heat transfer have to be included.

Hydrodynamics
The liquid shear and velocity profiles are obtained using Navier-Stockes equations, which state that for an element volume of moving fluid, the rate of change of momentum equals the net of the forces acting on the element volume [15] mathematically expressed with respect to Figure 3 as: Rate of change of momentum per unit volume = viscous force on element per unit volume + pressure force on element per unit volume + gravity force on element per unit volume Which is for steady state one-dimensional flow reduces to: Or in the form of shear stress: For a constant axial pressure gradient equation ( 5) becomes: where, V ¼ r L g y À ∂P ∂y ¼ constant.The gas phase exerts an interfacial shear R i on the liquid film, and the velocity varnishes at the solid boundary, so according to Figure 3, we can write the boundary conditions as: Hanratty and Engen [16], Cohen and Hanratty [17], and others have measured interfacial shear in two phase gasliquid flows, reporting the interfacial friction factor f i defined by: Integrating equation (7) twice, and applying boundary conditions, we obtain: Substitute equation ( 8) into (11), to obtain, C 2 = 0 Substitute equation ( 9) into (10), to obtain, Substitute equation (12) into equation (11) yields: or where Equation ( 13) presents velocity profile in the liquid film.
To develop an equation describes the shear profile through the liquid film; equations (10) and (12) are used: At the wall, equation ( 15) becomes, The film thickness is obtained from the volumetric flow rate per unit length of wetted perimeter, Q L , which is related to the velocity distribution by the following equation: Substitute equation (13) into equation (17) yields: Substitute equations (13a) and (13b) into equation (18) yields: Equation ( 19) is cubic equation in d which is readily solved by substituting the numeric values of Q L , R i , V and m L .
To evaluate the interfacial friction factor (f i ), for a turbulent gas, Blasius equation [18] is utilized: E is a correction factor introduced by our model to account for the gas entrance effect, since the gas flow in this study is not a fully developed turbulent flow.It is assumed further that E decreases exponentially with the axial distance from the contactor entrance: The constants B 0 and B 1 are determined later on from the experimental work.

Mass transfer
The general behavior is considered in the mass transfer model, for which the gaseous solute diffuses in the liquid film and then reaction occurred, so the process of diffusion and chemical reaction can still be represented by an extension of the film theory by a method according to Hatta [19].Thus a zone of reaction between gaseous solute (CO 2 ) and liquid reactant (MEA), which moves away from the gas-liquid interface taking up some position towards the bulk of the liquid.The final location of this reaction zone will be such that the rate of diffusion of CO 2 from the gasliquid interface equals the rate of diffusion of MEA from the main body of the liquid.
Figure 4 presents the paths of CO 2 , MEA, and liquid product in the liquid form, so CO 2 diffuses through the gas film as a result of driving force (P G ÀPi) and diffuses to the reaction zone as a result of equilibrium with Ci.The MEA diffuses from the bulk liquid to the reaction zone under a driving force, q, and nonvolatile product diffuses back to the main bulk of liquid under a driving force (mÀn).

Concentration profile
To determine the concentration profile, the column is divided into N sections of annular segments as in Figure 3, and a differential (mass or mole) balance equation is applied for each segments.Therefore two differential balances sufficient to determine the concentration profiles, an equation for CO 2 (A component) and the second equation is used for MEA (B component), both in liquid phase.For the element shown in Figure 3, a partial differential element of liquid with dimensions of dx Â dy Â 1 is considered.The amount of the liquid carried in y direction by bulk flow, plus the amount of liquid diffusing in x direction, should be sufficient to maintain the proceeding of the chemical reaction.Symbolically: Input rate À output rate + consumption rate = accumulation Accumulation = 0 ∴Input rate À Output rate + consumption rate = 0 Similarly for liquid reactant (B), the differential balance equation is as follows: Yih and Liu [20], Yih and Seagrave [21] studied the effect of gas-liquid interfacial shear on mass transfer through the liquid film, and concluded that shear stress induces a greater turbulence in the liquid phase, decreasing its resistance to mass transfer, and eddy diffusion, D E cannot be disregarded with respect to molecular diffusion.Therefore, any eddy diffusion term, D E , as a function of the radial coordinate X, is introduced in the microscopic mass transfer, and diffusivities in equations ( 22) and ( 23) become, Rearranging equations ( 22) and ( 23) and introduce equation ( 24) yields:

Boundary conditions
Boundary conditions are required to solve equations ( 25) and ( 26): where Equation (28) indicates that the components (A and /or B), mathematically and in fact, cannot pass through the wall.Equation (29) indicates that A crosses the interface at a rate depends on mass transfer of A, and that B is not volatile (assumption number 4 in section 2.1) and cannot cross the interface.

Equilibrium conditions (a) Rate of absorption
In the steady state absorption process, the rate of transfer of material through the gas film is the same as that through the liquid film, and the general equation of mass transfer [18] may be written as: where: N CO 2 : molar flux of gaseous reactant, kmol/m 2 .s(CO 2 ) K OG : overall gas phase mass transfer coefficient, m/s R g : gas-liquid interface temperature, K°P G CO 2 : partial pressure of CO 2 in gas phase, N/m 2 P Ã CO 2 : partial pressure of CO 2 in gas phase which is in equilibrium with bulk liquid phase, N/m 2 .
(b) Henry's constant Danckwerts [22] mentioned that in many cases, so long as the concentration of dissolved gas is small and the temperature and pressure are far removed from the critical temperature and pressure of the gas, Henry's law is obeyed and the concentration C CO 2 of dissolved gas is in equilibrium with a partial pressure Pi of the gas as given by where: C CO 2 : concentration of unreacted CO 2 in solution, kmol/m 3 H e : Henry's law constant, (N/m 2 ), m 3 /kmol) If the gas reacts in solution, Henry's law does not apply to the total concentration of unreacted gas.
For the CO 2 -MEA system, at the temperature range of 0-60 °C and 6-30% MEA Yunda et al. [23] reported the following relation: where: T S : is the gas-liquid interface temperature (K).
(c) Conversion of CO 2 When CO 2 is absorbed into aqueous MEA solution, the following overall reaction should be considered : The form of the equilibrium constant which was given by Danckwerts and Sharma [24] is: Savage and Kim [25] presented the following relation for the equilibrium vapor pressure of CO 2 : where: K: dissociation constant Alatiqi et al. [26] stated the following relation between K and Ts is applicable: Substituting the numeric values of C 1 , C 2 and C 3 in equation ( 36) yields: and equation ( 32) is rewritten as Thus Substituting equation ( 40) into equation ( 35) yields: Substituting equation (41) into equation (30) yields Astarita et al. [27] observed that up to 20% amine concentration the Vapor Liquid Equilibrium (VLE) remained essentially unchanged.For our experimental conditions the above analysis is used for VLE since the amine concentration varies from 3% to 18% in this study.
To obtain the molar conversion of CO 2 : where: A S : total area of mass transfer, m 2 M CO2 : Initial molar rate of CO 2 , kmol/s.

Mass transfer coefficients
(a) Gas-Side Mass Transfer Coefficient The gas-side mass transfer coefficient is estimated using McReady and Hanratty equation [28]: where: k G : gas-side mass transfer coefficient, m/s U*: gas-side interfacial velocity, m/s I: a correction term for entrance effect on heat transfer determined later from experimental results.
Yih and Liu [29] defined friction velocity as: where: gc: conversion factor = kg•m N•sec 2 ¼ 1 Substitute equation ( 9) into (45) and rearrange to obtain: Substituting equation ( 46) into (44) yields: (b) Liquid side mass transfer coefficient To estimate the liquid side mass transfer coefficient, the correlation of Vivian and Peaceman [30] is used.They investigated the characteristics of CO 2 -H 2 O system in a wetted wall column and suggested the following equation: where: The overall mass transfer coefficient for the absorbed component CO 2 consists of the gas, and liquid side mass transfer coefficients and is expressed as: Substituting equation ( 49) into equation ( 42): Where R refers to HOCH 2 CH 2 .At very short exposure times of the liquid to the gas encountered in industrial absorber, the effect of reaction 53 can be neglected, and only reaction 52 influences the absorption rate of carbon dioxide [32].Reaction 52 takes place in two steps: Danckwerts [33] reported that the equilibrium constants of reactions 54 and 55 are 10 À5 and 10 10 l/gmol, respectively.The second step, equation (55), is ionic and instantaneous; whereas the first step, equation ( 54) is second order, that is first order with respect to carbon dioxide and MEA, and rate controlling.Thus, the overall reaction 52 is considered as an irreversible, instantaneous and second order reaction with a stoichiometric coefficient of 2 and the equilibrium constant is 10 5 l/gmol.When the absorption of solute gas A is accompanied by an irreversible reaction, the average rate of absorption, NA, is represented by the following equation: For absorption with a second order reaction of the following stoichiometry:, The penetration model solution for the reaction factor is closely approximated by the following equation [34]: h is a dimensionless parameter given by the following equation: where b ∞ is the reaction factor for the absorption accompanied with an instantaneous reaction and expressed in terms of s as follows:

Rate of reaction
Caplow [35] and Danckwerts [22] suggested that CO 2 reacts with amine to form an intermediate called a "Zwitterion" The Zwitterion can be deprotonated by any base may present in the solution producing a carbamate ion and a protonated base according to the following equation: In this study, the homogeneous catalysis mechanism which is developed by Haruo et al. [31] is considered.Blauwhoff et al. [36] stated that according to this mechanism the rate of the forward reaction in equation ( 55) is developed based on the assumption of a pseudo steady state conditions for the Zwitteron concentration [7] as follows: where: k b : reaction rate constant for Zwitterion deprotonated, m 3 /mol.sk -1 : reverse reaction rate constant, s À1 k 2 : second order reaction rate constant, m 3 /mol.sDepending on the relative magnitude of the individual reaction rate constants, different orders and rate constants may obtained with different amines.Values of k b for various amines (bases) were obtained by Blauwhoff et al. [36] and Versteeg and Van Swaaij [37] correlated the rate constant k 2 with liquid temperature and basicity of amines for five different types of amines and obtained the following correlation: T: temperature, K

Heat transfer
Figure 5 shows a heat balance on a liquid element of dimensions d x * d y *1: -Heat input by convection + Heat input by conduction: -Heat generated Heat input À Heat output + Heat generated = 0 Dividing by d x .dy and rearranging yields: Equation ( 68) is used to develop the temperature profile through the column.

Boundary conditions
Boundary conditions are required to solve equation (68): where: T G : temperature of the gas, °K h G : gas heat transfer coefficient, W/m 2 .°KCondition b is based on the assumption that the column is well insulated from ambient while condition c assumes that heat transfers from the liquid to the gas.

Gas-side heat transfer coefficient
Bird et al. [15] stated that the solutions of many mass transfer problems at low mass transfer rates can be obtained by analogy with corresponding problems in heat transfer.Based on this principles, Mc Ready and Hanratty [28] introduced the following equation for gas side heat transfer coefficient: where: P rG : Prandtle number of gas side, C P m k À Á G k: Gas thermal conductivity, W/m.°K, C P : heat capacity of gas, m: dynamic viscosity of gas, Referring to Figure 5, the heat flux q into the liquid is the sum of heat conducted across the gas film (q cond. ) and heat of absorption which is proportional to the molar flux of solute into the liquid at the interface N CO 2 ½ , symbolically speaking: where: DH S : heat of absorption of CO 2 in MEA, kJ/kmol.NCO 2 (J): molar flux of CO 2 into the liquid at segment J of the contactor, kmol/m 2 .s(J): 1, 2,….N N: no. of annular segments of liquid film along the contactor Using finite difference operator to approximate the differential operator, equation (73) becomes: where: h x : grid spacing in x-direction, m.T NGX-1 : Liquid temperature at distance h x from T S , K N GX : grid points in x-direction Substitute equation (NCO 2 ) 51 in equation (74) yields:

See equation 75 below page
where: To obtain TS, applying Newton-Raphson method on equation (75) yields: Derivation of equation (74) w.r.t T S gives: where: T S0 : the initial estimate of T S Equation ( 78) is the general formula of Newton-Raphson iteration method for obtaining the constants of equation (76).Equations ( 76), (77), and (78) are used in the computer main program to estimate the interfacial gasliquid temperature.

Solution of the mathematical model
Equations ( 24), ( 25) and (68) are two dimensional second order nonlinear partial differential equation include: Through the liquid film, where f, F and S are functions of two independent variables (i.e., x and y).The above equations cannot be solved analytically.So this is done numerically.

Concentration profiles
Concentrations profiles of CO 2 are obtained by applying the implicit finite difference multigrid method developed by Wesseling [38].This method implies the solution of N equations (N being the number of annular segments) of CO 2 component and the same number of equations for the amine component.The systems are solved by matrix calculus.For example, to solve equation ( 24) by the multigrid method, one needs to calculate the coefficients of the differential equation and the parameters of its boundary conditions; such as U y , D E , D A , D B , h G , K G , etc.

Velocity profiles
In order to obtain the velocity profiles given by equation (13), the shear at the gas-side interface, R i , is given by equation ( 9).This value is related to the friction factor, w i , by equation (20).The film thickness, d, is calculated from 19 by iteration.

Estimation of model parameters
The developed mathematical model has three parameters, B 0 , B 1 , and I, whose numerical values are unknown.They must be estimated using experimental data.The problem of fitting the model to experimental data is to find the best value of the parameters that will predict absorption rates of CO 2 which are as close as possible to the experimental values.To develop a comprehensive mathematical model, a mass and an energy balance are conducted on the first annular segment and the outputs from those calculations are used as inputs of the second segment.Mass and energy balance are then conducted on the second segment and so on down the contactor.Velocity profile, film thickness, and physical properties of the materials are estimated for the first segment before applying the balances, but are adjusted for the next segment according to the output temperature and mass flow from each previous segment.This method of calculation gives concentration and temperature profiles in the liquid film of the contactor.The mathematical model developed has the predictive capability to calculate the radial profiles versus contactor height, as well as longitudinal profiles of conversion and temperatures.

Experimental work
This part of the study is carried out to obtain experimental data on the absorption reaction of CO 2 in aqueous solutions of MEA.Air stream containing CO 2 enters at the top of contactor and flows downward co-currently with the aqueous MEA solution.The alkaline solution flows downward through the absorber, and leaves at the bottom of the column.The contactor was designed and installed to operate at different ranges of gas and liquid flow rates.All the heat evolved was transferred to the cooling water which flows in a counter current mode through the shell side.Variation in temperature of inlet liquid was small and can be neglected.Furthermore, the temperature of liquor feed was controlled externally to give a temperature within the range of 25-60 °C (± 2 °C).In this study, four parameters are used to obtain concentration and temperature profiles namely, mole fraction of carbon dioxide in gaseous mixture (Y) 0.05, 0.075, and 0.1, molar concentration of absorbent

Materials
Carbon dioxide is supplied as 20 kg cylinder and air is supplied via an air compressor.Monoethanol amine, chemical indicators and reagents are supplied by Sigma Aldrich.The specifications of all gases and chemicals used in this study are presented in Table 1.

Apparatus and procedure
A schematic flow diagram of the experimental apparatus is shown in Figure 5.The main part of this apparatus is the monotubular falling film contactor.The liquid and gas streams are introduced to the reactor at the required temperature.The liquid MEA is introduced at the top of the column section through a pan distributor made of stainless steel.Hobler et al. [39] showed many different types of liquid distributor for wetting vertical tubes internally.He reviewed their dimensions and operating conditions.In this study a distributing device was used for uniform wetting of tube at a low degree of wetting.This distributor was presented by Kharisov and Kogan [40].It is characterized by simplicity in design and wetting intensity is within the liquid loading condition of the present experimental runs.Basically, it consists of a cylindrical pan with 46 holes (1.0 mm diameter) drilled in the base.The gas was introduced 30 mm below the falling section using spherical distributor of 30 mm diameter, which gives good gas distribution.Carbon dioxide is fed through a regulator and rotameter (FI-1) and mixed with air (supplied by an air compressor) in a homogenizing mixer to obtain a mixture at the desired concentration, to flow co-currently with the liquid downwards the contactor.The mixer is a cylindrical stainless steel box filled with a stainless steel gauze packing.The inlet section, before the falling film zone, is long enough to humidify the gas mixture with the amine solution.The temperatures of the gas and liquid are measured at the inlet and the outlet of the contactor.The temperature of the gas feeding stream was kept approximately constant at 35 °C ± 2 °C using a temperature controller.The absorption column is made of stainless steel AISI 316L tube of nominal diameter 1 in (ANSI B36.19); actual inner diameter 0.0223 m, with a wall thickness of 3 mm, the length of the contactor is 2.5 m.The absorber is supplied  with a cooling jacket, for the purpose of exothermic reactions to be performed.The dimensions of the implanted design are shown in Figure 7.The outlet solution goes to a regeneration tank (200 L), where liquid samples are collected for analysis at the entrance of the tank.The total amine concentration in the collected solution is found by titration with 0.1 N HCl to the end point of methyl orange as an indicator.To determine free amine content in presence of carbonate, the latter precipitated with excess of 10% BaCl 2 and titrated with 0.1 N HCl using phenolphthalein as an indicator.To calculate the amount of carbon dioxide absorbed, m, 5 ml of the collected solution is titrated with 0.1 N HCl against methyl orange as an indicator.Consumption of 0.1 N HCl = a ml.A second 5 ml sample of collected solution of a second liquor sample is added to 20 ml BaCl 2 solution.After 10 min, the sample is titrated with 0.1N HCl using phenolphthalein as an indicator.Consumption of 0.1 N HCl = b ml [41].
The amount of CO 2 absorbed is: A computerized UV-spectrophotometer (Jasco 520, manufacturer Jasco, Japan) is used to measure the concentration of MEA solution.Specifications of the devices and control elements are presented Table 2.

Results and discussion
The operating conditions and the experimental CO 2 conversions are fed to a MatLab code programed based on the models developed in section 2. A predicted CO 2 conversion is obtained and plotted versus the experimental conversion with a 3.978 % standard error.This statistical function gives a first cross check for the reliability of the developed model since the limitation of acceptance for the standard error test is 10%.To observe the behavior of the developed model at different operating conditions compared to the experiment results, carbon dioxide conversion is considered as a criteria of comparison.

Effect of MEA temperature
Reaction temperature has a significant effect on conversion.With an increase in MEA reaction temperature, the rate of chemical reaction increases and viscosity of the solution decreases, which leads to an increase in absorption rate and consequently CO 2 conversion.These results are in agreement with Thomas [42], Maddox [43], and Abid [44]. Figure 8 shows the impact of MEA temperature on CO 2 conversion in MEA solution at different liquid flow rates, constant mole fraction of CO 2 (Y = 0.1) and constant absorbent concentration (3 M).The same trend is observed in Figure 9 at different liquid concentrations, at Y = 0.075 and constant liquid flow rate (Q = 15 L/hr).The highest conversion is obtained at 3 M and 20 L/hr MEA compared to the minimum conversion obtained at a 1 M and 10 L/hr.

Effect of mole fraction of carbon dioxide in gas phase
Experimental and predicted results show a proportional relation between mole fraction of carbon dioxide and CO 2 conversion, starting from 0.05 to 0.1, this behavior is shown in Figure 10 and are in agreement with the results of Hikita et al. [45], Tontiwachwuthikul et al. [46] and Abid [44].relevant literature [6,47,48].Studies of Abid [44], Hitchcock and Cadot [49] and Maloney [50] showed that there is a critical of CO 2 conversion appears upon running at a progressive increase of amine concentration, then it falls to an equilibrium conversion.This behavior can be attributed to following fact, as increasing of MEA concentration can give more chances for the reaction to be occurred, also there is an observed increase in MEA viscosity that results in an additional resistance to solute diffusion.

The developed model 4.4.1 Validation of the developed model
To verify the validity of the developed model, standard error has been determined for the results plotted in Figures 8 through 11 at different operating conditions.Table 3 presents the standard errors obtained accordingly at a 95% confidence level.It is seen that all the standard errors obtained are acceptable (less than 10%).Herby, the developed is reliable.

Results of the developed model
The present model is developed to predict radial and axial velocity, concentration, and temperature profiles through the liquid falling-film along the contactor.The axial conversion of the solute in the gas phase and the axial interfacial temperature can be also predicted.A trial experiment (Y = 0.1, Q L = 20 L/hr and T = 55 °C) is used as an input condition to the developed model.

Interfacial temperature profile
Figure 12 shows the interfacial temperature along the contactor predicted by the model.It is seen that the average increase of interfacial temperature is 3 K through the first meter of the contactor while cooling effect of the gas stream becomes obvious after the first meter and diminishes the interfacial temperature lower than the initial input temperature of liquid.

Concentration distribution
The present model has the capability to predict concentration profile of the reactants through the liquid film, axially and radially.Figures 13 and 14 show a linear surface distribution of free CO 2 in the liquid film and the radial distribution of CO 2 along the contactor.The concentration profiles indicate that the reaction between CO 2 and MEA is instantaneous which is consistent with Davis [51] theory.Figure 15 shows a surface distribution of free amine concentration in the liquid film and Figure 16 shows the radial distribution of free amine in the liquid film at different mole fractions of CO 2 .They indicate that the present model described the effect of amine in liquid solutions from the viewpoint of homogeneous catalysis mechanism theory [18].This theory states that the chemical reaction is fast and mass transfer processes are limited in a small region close to the gas-liquid interface due to the bulk flow of amine, named boundary layer.Thus, the CO 2 + MEA reaction is fast enough that CO 2        approaches equilibrium with the other components of the solution at its interface composition.This is consistent with King [52] theory.Also it is seen in Figure 15 that the assumptions of insignificant depletion of MEA at the interface and the reaction products do not build up at the interface are confirmed and the concentrations of every component in the solution, except CO 2 , can be considered constant in the liquid boundary layer.

Temperature distribution
Figure 17 presents the radial distribution of liquid temperature at different heights from the contactor entrance.Figure 18 presents linear surface distribution of liquid temperature along the contactor.The following observations can be noticed from these figures: -The liquid temperature distribution presents a typical bulge shape.This is due to the fact that the exothermic reaction produces heat and the heat is absorbed by the liquid and gas streams.This means that the liquid absorbs most of the heat in the top of the contactor.
-At constant CO 2 content in gas stream, the liquid flow in the column is relatively small, thus the liquid heat capacity is small.As a consequence, the liquid temperature gradient at the top of the absorber is very steep.Even though the reaction occurs throughout the column, the temperature starts to decrease in the middle of the column, because the liquid is cooled by the flue gas flowing up the column till a constant value of temperature.
-The average increase of liquid temperature along the first meter of the contactor length is within 2 K.

Velocity distribution
Figure 19 presents the linear surface distribution of liquid velocity along the contactor.It can be seen from equation ( 6) that the maximum liquid velocity is obtained at the gas-liquid interface and the minimum value (i.e., Uy = 0) is obtained at tube wall, equation ( 6) is utilized the present model to predict velocity profile in the liquid film along the contactor.The parameters influences velocity profile are gas flow rate, liquid flow rate, liquid viscosity, and tube wall roughness.Liquid velocity is proportional to gas and liquid flow rate while it is inversely proportional to liquid viscosity and tube wall roughness.

Equilibrium concentrations
The present model is capable in predicting the equilibrium concentrations of free amine, carbamate, protonated amine, carbonate, hydroxyl and hydrogen ions for each experimental run.Starting with given values of MEA molarity and amine conversion in solution, one can calculate the equilibrium concentrations of ions for hydrogen, hydroxyl, carbonate, bicarbonate, protonated amine and free amine according to the following equations: They defined the carbonate conversion (x 0 ) as moles CO 2 absorbed per mole MEA, chemically expressed as: And the hydroxyl ion concentration was derived as a simple function of the carbonate conversion(x 0 ): equilibrium constant given by Kent and Eisenberg [53] as a function of temperature.Now, if we substitute equation (84) into equation ( 83) and rearranging yields: Substituting equations ( 83) and (84) into equation (81) and rearranging yields: Carbamate hydrolysis is expressed as follows: Apply electro neutrality yields: Apply a material balance yields: Apply equilibrium on equations (81) and (87) and combine yields: Combine equation ( 87) and ( 88) to obtain H + : Combine equation ( 89) and (90) yields: Substitute equations ( 91) and (93) into equation (92) yields: Substitute equations ( 91), ( 93) and (92) into equation (89) yields: where: Simplifying equation (95) and put K 2 K 3 ¼ 0 yields: Equation ( 99) is a quadratic equation in R and can be solved easily as follows: where:W = (A 2 + 1) .
By successive substitution of equation (100) into equations (93), (91) and (92), the numeric values of R -, H + , and R + are estimated as shown in Table 4.The predicted concentrations indicate a proportional relationship between concentration of bicarbonate ion and carbonate conversion.The concentration of hydroxyl ion is inversely proportional to carbonate conversion.Other concentrations of ions like free amine, carbamate, protonated amine, and carbonate are expressed chemically by equations ( 84), (85), (86), (91), ( 93) and (94) respectively.

Dependency of model parameters on experimental operating conditions
The present model has three parameters in equation (21).They are predicted for each experiment and tabulated in Table 5.The model parameter E, which is introduced for compensation of gas entrance effect on gas phase heat transfer coefficient, has the following behavior with operating conditions: -An average increase of carbon dioxide mole fraction by 33.5% results in an average increase of parameter value by 2.24%.This model behavior is due to the fact that the increase in concentration gradient of carbon dioxide gas results in an increase in absorption rate.-An average increase of the liquid flow rate by 46.7% results in an average increase of parameter value 22.6%.This model behavior results from the inversely thickness of the liquid film which means more resistance to heat transfer.-An average increase of liquid temperature by 24.5%, results in an average increase of parameter value by 62.5%.This model behavior is due to the fact that the positive temperature gradient between liquid and gas has predominant effect on heat transfer and in improving gas solubility.-The amine concentration has a positive effect on the model parameter, for the range of (1-3) M MEA, the model parameter has an average increase of 18.6%.

Conclusion
The mathematical model developed for the absorption process of CO 2 gas by MEA solution in a falling film contactor is able to reproduce some plant data, if a correction is applied to the kinetics of the CO 2 +MEA reaction obtained from the experimental results.The model uses kinetics to represent the mass transfer in the absorber.The following findings are concluded: -The adjusted kinetics are based on falling film contactor experiments.The validity of the model is closely related to the quality and accuracy of those data.More tubes in the falling film absorber data would allow for the regression of a more reliable rate constant.-The mass transfer in the absorber can be controlled by gas phase resistance.The correlation of Yunda et al. [54] for K G is empirical, thus it is accurate for a specific flow range.-The entrance effect of the falling film contactor can be related to the axial distance from the reactor entrance exponentially.-The reaction between CO 2 and MEA is instantaneous, since the reaction takes place at the interface [55,56].

Fig. 3 .
Fig. 3. Falling film contactor divided into N annular segments each of length y and thickness.

Figure 11
Figure 11  shows the impact of MEA concentration on the conversion of CO 2 at different MEA flow rates, constant MEA temperature (T = 55 °C) and constant mole fraction of carbon dioxide (Y = 0.1).It shows a good reliability of the predicted conversion with the experimental results.The highest conversion is obtained at the same liquid concentration (i.e., of 3 M MEA), which in line with

Fig. 8 .
Fig. 8. CO 2 conversion vs. liquid temperature for different liquid flow rates at constant mole.

Fig. 9 .
Fig. 9. CO 2 conversion vs. liquid temperature for different liquid concentration at constant liquid.

Fig. 10 .
Fig. 10.CO 2 conversion vs. mole fraction of CO 2 in gas phase for different liquid flow rates at constant liquid concentration (3 M) and at constant liquid temperature (T = 55 °C).

Fig. 12 .
Fig. 12. Interfacial temperature profile along the contactor length as predicted by present model for CO 2 fraction 0.1, QL = 20 L/hr and CMEA = 3 M.

Fig. 13 .
Fig. 13.Surface distribution of free CO 2 concentration in the liquid film for CO 2 fraction 0.1, QL = 20 L/hr and CMEA = 3 M as predicted by present model.

FilmFig. 14 .
Fig. 14.Radial distribution of free CO 2 in liquid film for CO 2 fraction 0.1, QL = 20 L/hr and CMEA = 3 M as predicted by present model.

Fig. 15 .
Fig. 15.Linear surface distribution of free amine concentration through the liquid film along the contactor for CO 2 fraction 0.1, QL = 20 L/hr and CMEA = 3 M.

Fig. 18 .
Fig. 18.Linear surface distribution of temperature through the liquid along the contactor.
solute in liquid, m 2 /s g y : gravitational constant = 9.78 m/sec 2 .L: length of wetted wall, m.G: mass rate of flow of liquid per wetted length, kg/m.sm L : viscosity of liquid, N•s/m 2 .

Table 1 .
Specification of ex perimental apparatus.

Table 2 .
Specification of chemicals used.

Table 3 .
Statistical Comparison between the present model and experimental results.

Table 4 .
Equilibri um concentration of ions as predicted by the present model.Gheni et al.: Oil & Gas Science and Technology -Rev.IFP Energies nouvelles 73, 43 (2018)

Table 5 .
Model parameters estimated by minimization technique.The axial conversion of carbon dioxide in gas phase varies exponentially with reactor length.-The mathematical model describes the effect of the process variables, liquid temperature, liquid concentration, liquid flow rate, and CO 2 mole fraction in gas phase, on absorption rate of CO 2 in MEA solution.-The highest concentration of MEA solution, 3 M, is efficient in absorption than the lower one (1 M) successively, this is consistent with the recent researchers trend line in using high concentration of MEA solution range 2-5 M. Van der Geer, J., Hanraads, J.A.J., Lupton, R.A., 2010.The art of writing a scientific article.J. Sci.Commun.163, 51-59.Gheni et al.: Oil & Gas Science and Technology -Rev.IFP Energies nouvelles 73, 43 (2018) Reynold number, r L u L d e m L R 'e modified Reynold number, r L u L d e m L ð1ÀeÞ R ef modified Reynold number, r L u L d p 6m L ð1ÀeÞS f R e modified Reynold number, r L u L d e m L e LT sc Schmidt number, m D AB r