How to Model Sparging in COMSOL Multiphysics®

August 22, 2018

Sparging, a mass transfer process between a gas and a liquid, is commonly encountered in industrial applications such as beverage carbonation and photobioreactors, or even at home for aeration in aquariums. In this blog post, we detail how to model a type of sparging, carbonation, using the COMSOL Multiphysics® software.

What Is Sparging?

Sparging is the process of bubbling a gas into a liquid to either dissolve it in the liquid or to strip the liquid off of some dissolved components. As mentioned above, the process entails mass transfer from a gas into a liquid or vice versa. One example of sparging is the process of oxygen stripping, in which nitrogen is bubbled through a liquid to push out or remove oxygen and other dissolved gases. Another example is dissolution, where a gas is dissolved into a liquid. The dissolved species can further react with the other components in the liquid, altering the liquid’s chemical composition. Carbonation and aeration are common examples of such sparging/dissolution.

A close-up photo of carbonated water in a bottle.
Bubbles from carbonation in a bottle of sparkling water.

When we open a bottle of soda, sparkling water, or any carbonated beverage, we experience a “fizz” due to a pressure release accompanied by buoyancy, leading to the rise of carbon dioxide bubbles. Before we open the bottle, however, carbon dioxide is dissolved into the beverage (CO2 (aq)). This dissolution is achieved by passing carbon dioxide bubbles (usually at high pressures for supersaturation) through the beverage. When the carbon dioxide bubbles pass through the beverage, there is mass transfer from the gas to the liquid across the surface of the bubbles and the dissolved carbon dioxide reacts with water to form carbonic acid. This entire process is known as carbonation.

In this blog post, we will take a look at how to model the carbonation process via the dissolution of carbon dioxide in water. In other words, we will make virtual carbonated water.

A Closer Look at Carbonation

Consider a glass of water into which carbon dioxide is introduced from its base, as shown in the schematic below. The carbon dioxide bubbles rise due to buoyancy and move through water, transferring mass and momentum and finally exiting through the top surface. Note that the dissolved carbon dioxide reacts with water to produce carbonic acid. Therefore, in order to model and understand carbonation in this glass, we need to resolve:

  • Mass and momentum transport of water and carbon dioxide bubbles; i.e., multiphase flow
  • Mass transfer of carbon dioxide from bubbles to liquid; i.e., dissolution
  • Mass transport of dissolved carbon dioxide and any subsequent reactions that take place (like the production of carbonic acid)

A diagram explaining the carbonation process in a glass of sparkling water.
Schematic of carbonation in a glass of water, with carbon dioxide bubbles entering from the bottom.

Modeling Sparging in COMSOL Multiphysics®

Let’s take a look at the required interfaces and conditions in COMSOL Multiphysics to resolve the above-mentioned aspects of sparging. The problem, as depicted in the above schematic, is axisymmetric. Therefore, we can set up the model using a 2D axisymmetric component. In terms of modeling multiphase flow, we use the Bubbly Flow, Laminar Flow interface because the volume fraction of bubbles is less than 1% and the density of carbon dioxide gas is negligible compared to the water. When the volume fraction of bubbles is less than 1%, particle tracing can also be used to solve for the transport of bubbles. However, note that bubbly flow can be used to model systems where the bubble volume fraction is greater than 1% (around 10% or less).

In addition, laminar flow has been chosen because the Reynolds number, calculated based on the bubble velocity, is lower than 2000 for the internal flow.

Tip: Read a blog post on the different types of multiphase flow interfaces in the COMSOL® software, including their advantages, disadvantages, and the corresponding underlying assumptions.

The Bubbly Flow, Laminar Flow interface solves for velocity, pressure, and the effective gas density. It also allows for the computation of the number density of gas bubbles and the relevant interfacial area (when the Solve for interfacial area option is chosen in the settings). This enables the computation of mass transfer between the bubbles and liquid using two-film theory, which will be discussed later.

A diagram showing how to set up a sparging model in COMSOL Multiphysics.
Schematic of the problem setup in COMSOL Multiphysics, with the relevant boundary and domain conditions for the Bubbly Flow and Transport of Diluted Species interfaces. Mass transfer of carbon dioxide is specified through a Reactions condition and the carbonic acid reaction is specified through an Equilibrium Reaction condition, both indicated by Reactions in this schematic. Note that the default domain conditions, such as Fluid Properties, Transport Properties, and Initial Values, are not indicated.

While modeling the fluid flow, note that we are assuming that the deformation of the free surface — i.e., the top surface of the water — is negligible. However, as this is a free surface, a Slip wall condition is chosen at the boundary, as shown in the figure above. This allows for the tangential motion of the water but restricts the normal velocity to zero. Gas inlet and outlet conditions are specified under the Wall boundary conditions in this interface. (More details about this process can be found in a previous blog post on modeling bubbly flow in a glass of beer.) The inlet gas flux is ramped from zero to the required value for consistency, a method of load ramping in time that is generally a good practice.

The transport of aqueous carbon dioxide and carbonic acid that results from the equilibrium reaction of carbon dioxide and water is modeled using the Transport of Diluted Species interface. This interface solves for the concentrations of aqueous carbon dioxide and carbonic acid. Note that we are only considering the reaction leading to the formation of carbonic acid in this case. The relevant boundary conditions for Bubbly Flow (bf) and Transport of Diluted Species (tds) are shown in the above figure.

The mass transfer of carbon dioxide from bubbles to liquid is modeled through a Mass Transfer condition in the Bubbly Flow interface. In this case, we use the Two-film theory option, which uses Henry’s law to compute the required mass transfer from gas to liquid, m_{gl}, given by:

\[{m_{gl}} = k({c^*} – c)M \cdot a\]

where:

  • k is the mass transfer coefficient that has to be specified, derived from experiments or literature
  • c is the concentration of the solute; in this case, the concentration of aqueous carbon dioxide, cCO2, that is computed through the Transport of Diluted Species interface
  • M is the molecular weight of the gas (carbon dioxide)
  • a is the interfacial area per volume, computed from the Bubbly Flow interface
  • {c^*} is the equilibrium concentration given by Henry’s law (c^* = (p+p_{ref})/H), where:
    • p is the relative pressure computed from the Bubbly Flow interface
    • p_{ref} is the reference pressure specified under the Bubbly Flow interface
    • H is the Henry’s constant derived from experiments or literature

The mass transfer from gas to liquid is a source of aqueous carbon dioxide and is accounted for in the Transport of Diluted Species interface through the Reactions term, as shown below. In addition to the coupling between the Transport of Diluted Species and Bubbly Flow interfaces, the velocity of the liquid phase computed from the Bubbly Flow interface is passed under the Transport Properties domain settings of the Transport of Diluted Species interface. This accounts for the convective transport of aqueous carbon dioxide.

An annotated screenshot of the Mass Transfer condition settings in COMSOL Multiphysics.
Settings of the Mass Transfer condition in the Bubbly Flow interface. All parameters specified in this domain feature are defined under the Parameters node, except for c_CO2, which is the concentration of carbon dioxide computed under the Transport of Diluted Species interface.

An annotated screenshot of the Reactions condition settings in COMSOL Multiphysics.
Settings for the Reactions term under the Transport of Diluted Species interface, where the mass transfer value, bf.mgl, is computed in the Bubbly Flow interface and passed as a source/sink term in the Transport of Diluted Species interface.

An annotated screenshot of the Equilibrium Reaction condition settings in COMSOL Multiphysics.
Settings to specify the equilibrium reaction H2O + CO2 (aq) ↔ H2CO3, with water as the solvent in the Transport of Diluted Species interface, where K_eq is defined in the Parameters node.

Results for the Sparging Process Simulation

The relevant inputs for mass transfer of carbon dioxide and the diffusion of aqueous carbon dioxide are taken from literature (Ref. 1–3).

Comparing plots at t = 0 s and t = 10 s in the figure below, it is clear that the carbon dioxide bubbles, quantified by the gas volume fraction, start rising up in the water and impart momentum to the water. The carbon dioxide bubbles exit through the top surface. However, as water does not flow out through the top surface, it starts recirculating in the glass leading to the vortex formations observed in the figures below. The mass transfer takes place in the region where the gas volume fraction is nonzero, so initially, the aqueous carbon dioxide concentration closely follows the gas volume fraction. However, the aqueous carbon dioxide is carried by convection and diffusion to the remaining portion of the glass. The predominance of convective transport is evident when you compare the velocity, including streamlines, and the aqueous carbon dioxide concentration surface plots below.

Sparging process simulation results at four different times.
Velocity magnitude of water along with the streamlines (left), carbon dioxide bubble volume fraction (center) and aqueous carbon dioxide concentration (right) at t = {0, 10, 50, 100} s.

In COMSOL Multiphysics, we can also create and export an animation of the plots above. This facilitates better visualization of the formation of vortices, their evolution in time, and the transport of the diluted species.

 

Time evolution of the velocity magnitude of water along with the streamlines (left), carbon dioxide bubble volume fraction (center), and aqueous carbon dioxide concentration (right).

As this is a 2D axisymmetric simulation, the COMSOL Multiphysics software automatically creates a 3D data set, which is a revolution of the 2D data set from the computed study. We can then create a 3D animation of the time evolution of carbonic acid, where the formation and transport of carbonic acid strongly follow the transport of aqueous carbon dioxide.

 

Time evolution of carbonic acid concentration.

Concluding Remarks

In this blog post, we discussed how to model sparging with an example of the carbonation of water. In general, a similar process can be used to model mass transfer between any two fluids in dispersed multiphase flows. Note that the Mass Transfer condition is available both in the Bubbly Flow and the Mixture Model interfaces. As is usually the case with COMSOL Multiphysics, you can also specify user-defined expressions under the Mass Transfer condition by scrolling down the list under Mass transfer model.

While modeling such multiphysics problems, it is useful to start with a simple model and then increase complexity for an optimal modeling workflow.

Next Steps

Interested in modeling sparging? Download a carbonation example by clicking the button below.

For more information on modeling dispersed and separated multiphase flows, check out these example models:

References

  1. Sander, R., Compilation of Henry’s law constants for inorganic and organic species of potential importance in environmental chemistry, Air Chemistry Dept., Max-Planck Inst. of Chemistry, 1999.
  2. Hogendoorn, J.A., Bhat, R.V., and Versteeg, G.F., “On the Determination of the Diffusivity of CO 2 in Aqueous and NonAqueous Solvents: Investigations with Laminar Jets and Wetted Wall Columns”, Chemical Engineering Communications, vol. 189, no. 8, pp. 1009–1037, 2002.
  3. Han, J., Eimer, D.A., and Melaaen, M.C., “Liquid phase mass transfer coefficient of carbon dioxide absorption by water droplet”, Energy Procedia, vol. 37, pp. 1728–1735, 2013.

Comments (16)

Leave a Comment
Log In | Registration
Loading...
tarkan ortac
tarkan ortac
August 22, 2018

Hi
That was great and could i have the example of this model on comsol?
I have some problem with this modelling ..

Siva Sashank Tholeti
Siva Sashank Tholeti
August 23, 2018

Hi Tarkan Ortac, thank you for your interest in this model. Yes, the model is available through the button “Get The Carbonation In Water Tutorial” towards the end of the blog. For your convenience I have also provided a link below:
https://www.comsol.com/model/carbonation-in-water-67701

mehran jooya
mehran jooya
November 29, 2019

Hi, thank you Mr.Siva Sashank Tholeti
I’m working on a specific airlift reactor. it was so helpful for me.

I have couple of questions and I’ll be thankful if you would help me.
1-If the bubble volume fraction is greater than 10%, can i use Turbulent Bubbly Flow interface or i should use The Mixture Model, Turbulent Flow interface?
2-I have a mass transfer process (carbon dioxide or oxygen) between a gas(air) and a liquid(water) like your tutorial. can i use (Transport of Diluted Species interface) with your setting (with my coefficient and variables)?
Regards

Siva Sashank Tholeti
Siva Sashank Tholeti
December 2, 2019 COMSOL Employee

Hi Mehran, yes, Mixture Model would be more appropriate for such cases, please see this article for more information on this:
https://www.comsol.com/blogs/which-multiphase-flow-interface-should-i-use/

Note that as long as the species is dilute (less than 10 mol%), it’s transport can be simulated with Transport of Diluted Species interface. Use Transport of Concentrated Species for species with concentrations beyond that (10 mol%).

Minh Tri Luu
Minh Tri Luu
December 19, 2019

Hi Siva,

Thanks a lots for this tutorial. I have a similar problem, however instead of co2 is moved from gas to liquid, I have Co2 from liquid transfer to gas. The liquid is initially supersaturated with CO2, when there are bubbles, the excess CO2 in liquid migrate to the bubble. I wonder if I change the initial CO2 concentration (let say 8 mol/m3), can I use your model to simulate the transfer of CO2 from liquid to bubble ? Thanks

Best
Minh

Siva Sashank Tholeti
Siva Sashank Tholeti
December 19, 2019 COMSOL Employee

Hi Minh,
Yes, you could use a similar approach as the mass transfer can be from continuous to dispersed phase as well (the direction of mass transfer is dependent on the the sign of the value you specify). But if the CO2 concentration that is being tracked is significant, i.e. greater than 10 mol%, you would have to use Transport of Concentrated Species interface. Please contact COMSOL Technical support (support@comsol.com) for more information. Thanks.

Christian McHugh
Christian McHugh
February 25, 2020

Hi, thank you for the demo. I want to modify this by turning off the bubbling source after 2 seconds. I would just monitor the dissolved CO2 concentration after this time. Where would I make this adjustment?

Siva Sashank Tholeti
Siva Sashank Tholeti
December 8, 2020 COMSOL Employee

Hi Christian,
Sorry for the late reply. To turn off the bubble source after 2 seconds, modify the bubble mass flux and number density under bottom wall (“Wall 2”) condition. The following article can be useful towards this end:
https://www.comsol.com/support/knowledgebase/1245

Thank you,
Siva

Emile Sabeti-Mehr
Emile Sabeti-Mehr
December 8, 2020

Hi Prof Tholeti,

I’m sorry to bother you regarding a different subject matter on your blog. I’m a graduate student in the Integrative Biology department at the University of Guelph, and I’ve been watching and following your online lectures on Modeling Nonisothermal Flow and Modeling External Turbulent Flows. I found them quite helpful in learning how to navigate in COMSOL as I am a new user for my research. The pandemic has delayed my experiments, so I have been focusing on COMSOL to measure drag and lift forces on an Ahmed body as a precursor to my work on freshwater mussels.

I started off with the Ahmed body on the Applications libraries under the Varification Examples in COMSOL and followed what you did here https://www.comsol.com/video/modeling-external-turbulent-flows-in-comsol-multiphysics. I reduced the size of the domain so that the computation time would be less. The system can compute and display all the plots at 40 m/s, however, when I changed the Inlet velocity, (i.e. from 40m/s to 10 m/s) I received the following error messages:

(1) Feature: Stationary Solver 1 (sol1/s1) Turbulence variables Undefined value found. – Detail: Undefined value found in the equation residual vector. There are 10 degrees of freedom giving NaN/Inf in the vector for the variable comp1.ep. at coordinates: (0.111272,1.05708,0.209503), (0.117805,1.05225,0.206438)

(2) Feature: Stationary Solver 1 (sol1/s1)
Turbulence variables

Undefined value found.

– Detail: Undefined value found in the equation residual vector.

There are 6 degrees of freedom giving NaN/Inf in the vector for the variable comp1.ep. at coordinates: (0.197968,0.591908,0.310229), (0.198074,0.582103,0.311047), (0.199137,0.596712,0.319952), (0.199308,0.591908,0.310229), (0.199455,0.582103,0.311047),

(3) Feature: Stationary Solver 1 (sol1/s1) Failed to find a solution for all parameters, even when using the minimum parameter step. Maximum number of segregated iterations reached. Not all parameter steps returned.

I followed the threads on the COMSOL site for solutions and implemented those suggestions, but the same errors keep coming up or the computation goes on for days and does not seem to progress at all.

I would really appreciate your help, if you could please point me in the right direction.

Thank you,

Emile

Siva Sashank Tholeti
Siva Sashank Tholeti
December 8, 2020 COMSOL Employee

Hi Emile,
Thank you for the message and happy to hear you found the lecture series useful. As we would have to look into the model to help with this, please forward this question to COMSOL support (support@comsol.com).

Thank you,
Siva

Jose Juan Bolivar Caballero
Jose Juan Bolivar Caballero
January 15, 2021

Hi!

I want to start by saying that I’m really happy that I found this. I’m doing my thesis and I’m working with algae cultivation. The system is very similar to this simulation with the difference that CO2 is consumed by the algae and the control over the pH is important. I have downloaded this simulation and since I’m interested in the pH value I want to add the association of H2CO3 into HCO3- and H+. I’m having problems with this, I understand that you define the rate of reaction from bf.mgl and then you add the equilibrium but I do not understand how to add the order reaction plus its equilibrium.

I also understand that the bubbles consist of 100% CO2. How can I change the composition of the bubbles? For example, 80% N2, 20 % CO2

Thank you beforehand

Siva Sashank Tholeti
Siva Sashank Tholeti
March 1, 2021 COMSOL Employee

Hi Jose, thank you your comment. As this is an involved question, please contact support@comsol.com for help with this.

Thank you,
Siva

Raghuram I
Raghuram I
February 26, 2021

Hi, Can I do this in 5.3a or 5.4 Version?

Siva Sashank Tholeti
Siva Sashank Tholeti
March 1, 2021 COMSOL Employee

Yes on both fronts

Yu Zhang
Yu Zhang
June 22, 2021

Hi, Siva:

Are there any ways to Explicitly simulate the same processs, i.e. bubble rising and dissolution, by using COMSOL? Would you mind giving me some clues about that? Thanks.

Siva Sashank Tholeti
Siva Sashank Tholeti
June 23, 2021 COMSOL Employee

Hi Yu,

Thank you for the question. If you are interested in modeling the rising of a bubble and tracking the fluid-fluid interface explicitly, then here is an example that models that using level-set interface (however mass transfer across the fluid-fluid interface is not modeled here):

https://www.comsol.com/model/rising-bubble-177

For more information on this and any other related questions, please feels free to contact support@comsol.com

Thank you,
Siva

EXPLORE COMSOL BLOG