A key step of the refined sodium bicarbonate (BIR ) production by the Solvay process consists in the gas-liquid mass transfer of carbonic gas CO 2 from a gaseous mixture of CO 2 and air to a sodium carbonate Na 2 CO 3 and bicarbonate NaHCO 3 brine. This transfer takes place in 20-m high and 2.5-m wide bubble columns (the BIR columns). This large size induces a large contact time between gas and liquid, in order to increase the amount of transferred CO 2 . Nevertheless, the gas phase leaving the columns contains an important quantity of CO 2 (only half of the CO 2 is transferred). It causes a huge CO 2 emission to the atmosphere (the equivalent of 150 Smart driving at 100 km/h, for a single column). More, the CO 2 is produced by lime calcination. This process requires a large amount of energy. It represents the major part of the energy consumption to produce BIR. In the past decades, several optimizations of the BIR process were performed by an empiric approach. There are today some limits to this approach for applications requiring high levels of purity and a well-defined granulometry. Moreover, the ecologic pressure becomes more and more constraining for the greenhouse effect gas emission. Accordingly, Solvay is seeking for a more fundamental approach. Currently, several studies are realized in the chemical engineering department at the ULB. The final goal is to create a complete model of a BIR bubble column, taking into account all the phenomena taking place in it. This model will be used in order to optimize the CO 2 transfer rate in the column. As a direct consequence of this optimization, the CO 2 emission and the energy consumption of the BIR production will hopefully be reduced.The objective of the work presented in this paper is to develop a model of the mean CO 2 flux density, expressed by unit of time and interfacial area, in a BIR bubble column, that will be integrated in the future complete model of a BIR column. No solid phase is considered in this paper.The CO 2 transfer rate from a bubble of gas to the brine is controlled by the physico-chemical phenomena occurring in the thin layer of liquid near the interface. After the CO 2 absorption in this liquid layer, the CO 2 is transferred by diffusion and is involved in several chemical reactions. These reactions modify (accelerate) the gas-liquid mass transfer rate.The liquid phase flow around the bubbles, rising up in the columns, influences the mass transfer rate too. The scale of this phenomenon is different from the physico-chemical phenomena. A multiscale approach is then followed.A mathematical modelling of the coupling between mass transfer and chemical reaction is first developed. The equations of the model are solved numerically, using COMSOL Multiphysics. In order to validate experimentally this model, a Mach-Zehnder interferometer is set up. The scale-up of this model to a bubble, in a second time, is carried out by completing the model of the coupling with a representation of the liquid phase flow around a bubble. It is called the bubble-li...