## Services on Demand

## Article

## Indicators

## Related links

- Similars in SciELO

## Share

## Brazilian Journal of Chemical Engineering

*Print version* ISSN 0104-6632

### Braz. J. Chem. Eng. vol.29 no.1 São Paulo Jan./Mar. 2012

#### http://dx.doi.org/10.1590/S0104-66322012000100012

**FLUID DYNAMICS; HEAT AND MASS TRANSFER; AND OTHER TOPICS**

**Simple model for CO _{2} absorption in a bubbling water column**

**I. Martínez ^{I,}^{*}; P. A. Casas^{II}**

^{I}Universidad Politécnica de Madrid, ETSIA, Pz. Cardenal Cisneros 3, 28040, Madrid, Spain. E-mail: isidoro.martinez@upm.es

^{II}Ministerio de Defensa, Instituto Tecnológico "La Marañosa", Carretera de San Martín de la Vega, km. 10, 28330, Madrid, Spain. E-mail: pcasalc@fn.mde.es

**ABSTRACT**

A simple analytical model is presented for the time-evolution of dissolved CO_{2} concentration in a batch water column. The application in mind is the removal of the CO_{2} produced by the propulsion plant of a submarine cruising underwater. The model is based on a constant number of constant size CO_{2} bubbles within the water column, with uniform gas injection through a porous membrane at the bottom, and bubbles flowing upwards by their own buoyancy through a short water column (the exhaust gas is recirculated). Analytical results are validated with experimental data found in the literature and the advantages of this simple analytical model pinpointed.

**Keywords:** Carbon dioxide water solutions; CO_{2} absorption; Bubble column; Dissolution; Batch; Absorption.

**INTRODUCTION **

Hydrocarbon-oxygen-based propulsion systems on board submarines generate great quantities of water and CO_{2} during operation; the latter cannot be stored aboard and must be expelled to outside seawater under submerged condition. Among the several options to remove CO_{2} from submarines cruising underwater (see Potter *et al*., 2000), the one selected here is to dissolve the gas in seawater in batch at low pressure in a bubble column (inside the submarine) and then to pump the water outside, aided by a work exchanger device to minimize pumping power.

Bubble generation is by the flow of gas through a porous membrane at the bottom of the water column, while buoyancy provides the bubble dispersion mechanism that enhances the gas absorption rate. Undissolved gas is collected at the top and recirculated.

Mean bubble radius, r_{0}, is governed by two basic phenomena (Parthasarathy *et al*., 1996, Ribeiro *et al*., 2004):

• Surface coalescence at the porous membrane before the bubble detaches. When the membrane pore size becomes small enough, d<<10^{-4} m, it appears that the size of the bubbles formed on the membrane tends to stabilize around a limiting value approximately equal to r0=0.19 mm (Parthasarathy *et al*., 1996).

• Flow coalescence of bubbles in the bulk. For large gas injection speeds, u_{g}>0.1 m/s, collisions between bubbles in the bulk water favor the formation of bubbles larger than when detached from the membrane surface (Moshtari *et al*., 2009) and the flow regime changes from homogeneous to heterogeneous. Flow coalescence increases with bubble number density (n/V).

The model developed here assumes that there is no flow coalescence (u_{g}<<0.1 m/s), and is intended to be applicable to small bubbles of any size (assumed spherical; in the range r_{0}=0.1 mm to 0.5 mm), although our interest is in small sizes, to maximize absorption rate (for the submarine application one has to dispose of at least 0.05 kg/s of CO_{2} , in a tightly constrained engine room). In spite of the many simplifications introduced into the model (uniform constant bubble radius, constant pressure, pure gas, pure water, no chemical reactions...), its predictions compare favourably with more sophisticated analytical and experimental results.

**MATHEMATICAL MODEL OF ABSORPTION RATE**

The diffusion mass flux jg of a dissolved gas in a liquid of density ρl, gas mass fraction w, and gas diffusivity D_{g}, in spherical coordinates, assuming spherical symmetry, is (Fick's law):

Integration in r from the bubble radius r0 to infinity, and application of mass balance of the gas (pure-CO_{2} ), of density ρ_{g}=p_{g}/(R_{g}T), assuming time independence and constant liquid density, as in Darmana *et al*. (2004), yields a simple expression for the total dissolution time, t_{d}:

Because we deal with short columns, we assume that the residence time of the bubbles in the column, t_{r}=L/u_{b}, is much smaller than the total dissolution time, td, and hence there is a nearly constant bubble radius along the column, since only a small amount of gas dissolves in one pass. The short-column approach also justifies the assumption of constant gas density along the column (small change in hydrostatic pressure). Besides, capillary overpressure inside the bubble is also neglected, because only radii r_{0}>0.1 mm are envisaged, as explained above).

The bubble rising speed, ub (terminal velocity) is determined by the buoyancy/friction matching, the first being the weight of liquid displaced by the bubble and the latter being the drag on a sphere given by Stoke's law for low Reynolds number flows (Re<10):

with g=9.8 m/s^{2}, and µ_{l} and ν_{l} being the dynamic and kinematic viscosities of the liquid. The time it takes for a detaching bubble to reach terminal speed is really negligible.

The number of bubbles, n, in a liquid column of height L and cross-sectional area A depends on feeding rate, measured by the gas injection speed at the bottom, u_{g}, and on bubble size r_{0}, as computed by the product of number of bubbles injected per unit time, (the volumetric flow rate is u_{g}A), times the residence time of bubbles in the column, t_{r}=L/u_{b}:

Gas diffusion from bubbles to liquid increases the solute concentration, measured by its mass fraction far from the bubble surface, w. The dissolving rate equation for the whole column is obtained from the solute mass balance:

Integrating Eq. (5) from w=0 at t=0 to a desired mass-fraction before venting, w, one gets the timeevolution for w:

For a given gas injection speed, u_{g}, and bubble radius, r_{0}, Eq. (6) allows one to find an appropriate venting time, t_{v}, to get rid of the CO_{2}-loaded batch water, before the bath gets near saturation and the mass transfer rate becomes impossible. The advantage of this simple analytical model is that the influence of the different controlling parameters is explicit in Eq. (6); i.e., we see that the required area (or the venting time) is proportional to the fourth power of r_{0}, and inversely proportional to u_{g}.

**VALIDATION OF RESULTS **

Experimental data presented in Parthasarathy *et al*. (1996) show that there is a minimum bubble radius (about r_{0min}=0.19 mm) attainable when injecting a gas through a porous membrane into a liquid, due to bubble coalescence at the interface before the bubble detaches.

Tokumura *et al*. (2007) present an involved theoretical model of the dynamics of CO_{2} absorption in water (including different species in solution and chemical reactions), and supporting experimental data. In particular, they show the effect of several gas injection rates (Q_{G0}=1 L/min to Q_{G0}=4 L/min in their Fig. 3b, included here), in the time-evolution of CO_{2} molar concentration (their dissolved inorganic carbon-dioxide-equivalent molarity, DIC), for a cross-sectional area A=75.4 cm^{2}. When we plot our result Eq. (6) in terms of molar concentration, DIC=w_{s}ρ_{l}/M_{g}, for the corresponding injection speeds u_{g}=Q_{G0}/A (2.2 mm/s to 8.8 mm/s), and the property values for their conditions (24ºC and 100 kPa, w_{0}=pS/ρ_{l}=(10^{5})(2×10^{-5})/10^{3}=0.002), we get Fig. 1a.

We can see from Fig. 1 that the behaviour for large injection rates is fully predicted: the water becomes nearly saturated in a couple of minutes (the small differences in saturation values are just due to typical uncertainties in material-property data). Our model is not so good for small injection rates because, in this case dissolution and residence times are of the same order of magnitude, which means that, at least in the initial stage of the process, bubbles completely dissolve before reaching the top of the column, invalidating our assumption of constant bubble radius. For large injection rates, however, the CO_{2} mass-fraction in the water column, w, increases rapidly, and bubble dissolution time also increases quickly by reduction of w_{0}-w in Eq. (2), with the consequence that the radius of bubbles detached from the membrane has a practically constant value during most of the absorption process.

With this assurance of the validity of our model Eq. (6) for high gas injection speeds (u_{g}>4 mm/s) and minimum bubble radius (r_{0}=0.19 mm), it is easy to make preliminary design calculations, as presented now for our application to CO_{2} removal in air-independent submarine propulsion.

We approximate seawater properties by typical plain water values: density ρ_{l}=1000 kg/m^{3}, kinematic viscosity ν_{l}=10^{-6} m^{2}/s; we take D_{g}=2.1×10^{-9} m^{2}/s for CO_{2} diffusion in water; gas density at working conditions ρ_{g}=p_{g}/(R_{g}T)=2×10^{5}/(189·288)=3.7 kg/m^{3}, since absorber operation is expected at a constant pressure p=200 kPa. The saturation mass fraction at this pressure is w_{0}=pS/ρ_{l}=(2×10^{5})(2×10^{-5})/10^{3}=0.004 (4 g of CO_{2} per kg of solution, i.e., a molar concentration c_{0}=w_{0}ρ_{l}/M_{g}=(0.004)(1000)/0.044=91 mol/m^{3}), where S=2.0×10^{-5} kg/(m^{3}·Pa) is the solubility of CO_{2} in water at 15ºC (Carrol *et al*. 1991).

We choose as design parameters for our submarine application to remove =0.05 kg/s (the CO_{2} produced by a 200 kW engine), a water column of L=0.50 m in height, and a gas injection speed of u_{g}=10 mm/s. Assuming a mean bubble radius of r_{0}=0.19 mm (the experimental limit of Parthasarathy *et al*., 1996), in the limit of pure water (w=0) we have:

• Bubble dissolution time, from Eq. (2), is t_{d}=7.9 s.

• Bubble rising speed, from Eq. (3), is u_{b}=0.080 m/s. It would take then a distance of L=u_{b}t_{d}=(0.08)(8)=0.64 m for a bubble to fully dissolve. The Reynolds number is Re=2u_{b}r_{0}/ν_{l}=30, above the Stokes law limit in Eq. (3), but its effect must compensate somehow in view of the good experimental agreement in Fig. 1. It can be further shown that the mean distance between bubble centres is 0.61 mm, and the mass-transfer area density a=4πr_{0}^{2}n/V=2000 m^{2}/m^{3}.

• Bubble residence time is t_{r}=L/u_{b}=0.5/0.080=6.2 s. Note that the radius decreases only 47% in spite of t_{r}/t_{d}78% because it is the bubble area which decreases linearly with time, as can be seen in Eq. (2).

• Mass-fraction of CO_{2} in the column after venting time t_{v} is, from Eq. (6), wv=0.004(1-exp(-Kt_{v})), with K=0.022 s^{-1}, so that, if we scavenge when the water becomes half-saturated (w_{v}=0.002), then the venting time is t_{v}=ln2/K=31 s.

• Required column cross-section area, A, to get rid of =0.05 kg/s:

which seems reasonable to be accommodated within a submarine engine room, although a smaller size would be preferable. To make the batch processing smoother and to enhance redundancy, 4 o 6 small and slender cylinders could be used to distribute that cross-sectional area.

**CONCLUSIONS **

A simple analytical model has been developed here to predict the time-evolution of CO_{2} absorption in water. This model clearly shows the effects of several parameters that influence the absorption rate: bubble size, injection rate, or temperature (through the change in material properties). Although our final goal was to have a simple formulation of the absorption rate, to be integrated in a larger numerical simulation of the whole CO_{2} removal system for a conventional submarine with air-independent propulsion, we believe that it can be of interest in many other applications within the stated assumptions (relatively large gas injection speeds, small bubble size, and relatively short batch water columns).

**NOMENCLATURE**

A | Area of cross-section | m^{2} |

D_{g} | Diffusion coefficient of CO_{2} | m^{2}/s |

gas in water | ||

d | Membrane pore size | m |

F | Drag | N |

g | Gravity acceleration | m/s^{2} |

j | Mass flow per unit surface | (kg/s)/m^{2} |

K | Mass transport coefficient | 1/s |

L | Bubble column height | m |

Mass flow rate | kg/s | |

n | Amount of substance | mol |

Molar flow rate | mol/s | |

p | Pressure | Pa |

R | Gas constant | J/(kg·K) |

r_{0} | Bubble radius | m |

r | Radial coordinate | m |

S | Solubility of gas in liquid | kg/(m^{3}·Pa) |

T | Temperature | K |

t | Time | s |

u | Velocity | m/s |

V | Volume | m^{3} |

w | Mass fraction |

*Greek Letters*

ρ | Density | kg/m^{3} |

µ | Dynamic viscosity of liquid | Pa·s |

ν | Kinematic viscosity of | m^{2}/s |

liquid |

*Subscripts*

b | Bubble |

c | CO_{2} |

d | Dissolving |

g | Gas (bubbles) |

l | Liquid |

r | Residence |

v | Venting |

0 | At the liquid/gas interface |

Far from the liquid/gas | |

interface |

**REFERENCES**

Carroll, J. J., Slupsky, J. D., Mather, A. E., The solubility of carbon dioxide in water at low pressure. Journal of Physical and Chemical Ref. Data, 20, (6), pp. 1201-1209 (1991). [ Links ]

Darmana, D., Deen, N. G., Kuipers, J. A. M., Modelling of mass transfer and chemical reactions in bubble column reactor using a discrete bubble model. Fifth International Conference on Multiphase Flow, Yokohama, Japan (2004). [ Links ]

Moshtari, B., Babakhani, E. G., Moghaddas, J. S., Experimental study of gas hold-up and bubble behavior in gas-liquid bubble column. Petroleum & Coal, 51, (1), pp. 27-32 (2009). [ Links ]

Parthasarathy, R., Ahmed, N., Size distribution of bubbles generated by fine-pore spargers. Journal of Chemical Engineering of Japan, 29, (6), pp. 1030-1034 (1996). [ Links ]

Potter, I. J., Clavelle, E., Reader, G. T., Kady, J., Carl, M., Exhaust gas management systems for underwater heat engines. In Proceedings of the 2000 International Symposium on Underwater Technology, (Cat. No.00EX418) (2000). [ Links ]

Ribeiro, C. P., Lage, P. L. C., Experimental study of bubble size distributions in a direct-contact evaporator. Brazilian Journal of Chemical Engineering, 21, (1), pp. 69-81 (2004). [ Links ]

Tokumura, M., Baba, M., Kawase, Y., Dynamic modeling and simulation of absorption of carbon dioxide into sea water. Chemical Engineering Science, 62, pp. 7305-7311 (2007). [ Links ]

(Submitted: March 10, 2011 ; Revised: July 22, 2011 ; Accepted: October 27, 2011)

* To whom correspondence should be addressed