The Bubbly Flow, k-ε Interface
The Bubbly Flow, k-ε(bf) interface (), found under the Multiphase Flow>Bubbly Flow>Bubbly Flow, Turbulent Flow branch () when adding a physics interface, is used to model the flow of liquids with dispersed bubbles at high Reynolds numbers.
It is assumed that the bubbles only occupy a small volume fraction and that they always travel with their terminal velocity. It is thereby possible to solve only one set of Navier–Stokes equations for the liquid phase and to let the velocity of the bubbles be guided by a slip model. The pressure distribution is calculated from a mixture-averaged continuity equation. The volume fraction of bubbles is tracked by solving a transport equation for the effective gas density. Turbulence effects are modeled using the standard two-equation k-ε model with realizability constraints and bubble-induced turbulence production. The flow near walls is modeled using wall functions.
The physics interface can also model the distribution of the number density, which can be used to calculate the interfacial area, useful when simulating chemical reactions in the mixture.
The main physics node is the Fluid Properties feature, which adds the equations for turbulent bubbly flow and provides an interface for defining the fluid materials for the liquid and the gas and the slip velocity model to use.
When this physics interface is added, the following default physics nodes are also added in the Model BuilderFluid Properties, Wall (the default boundary types are No slip for the liquid and No gas flux for the gas), and Initial Values.
Except where indicated below, the nodes settings for this physics interface are the same as for The Laminar Bubbly Flow Interface.
Turbulence
The default Turbulence model type is RANS. A different turbulence model can be selected under Turbulence model. The default turbulence model is k-ε.
Wall Treatment
Wall treatment for the k-ε model can only be set to Wall functions. More options become available by selecting another option under Turbulence model.
Edit Turbulence Model Parameters
Turbulence model parameters are optimized to fit as many flow types as possible, but for some special cases, better performance can be obtained by tuning the model parameters. For a description of the turbulence model and the included model parameters see Theory for the Turbulent Flow Interfaces.
Dependent Variables
The dependent variables (field variables) are the Velocity field, liquid phase u (SI unit: m/s), the Pressure p (SI unit: Pa), the Effective gas density rhogeff (SI unit: kg/m3), the Turbulent dissipation rate ep (SI unit: m2/s3), the Turbulent kinetic energy k (SI unit: m2/s2), and the Number density, gas phase nd (SI unit: 1/m3).
The names of variables can be changed but the names of fields and dependent variables must be unique within a component.
Consistent and Inconsistent Stabilization
To display this section, click the Show More Options button () and select Stabilization in the Show More Options dialog box. The settings for this section are the same as for The Laminar Bubbly Flow Interface with the addition of this section: stabilization for the turbulence variables in the Turbulence equations area.
When using a turbulence model, streamline and crosswind diffusion are by default applied to the turbulence equations.
Advanced Settings
To display this section, click the Show More Options button () and select Advanced Physics Options in the Show More Options dialog box. The Turbulence variables scale parameters subsection is available when the Turbulence model type is set to RANS.
In addition to the settings described for The Laminar Bubbly Flow Interface, enter a value for Uscale and Lfact under the Turbulence variables scale parameters subsection.
The Uscale and Lfact parameters are used to calculate absolute tolerances for the turbulence variables. The scaling parameters must only contain numerical values, units or parameters defined under Global Definitions. The scaling parameters cannot contain variables. The parameters are used when a new default solver for a transient study step is generated. If you change the parameters, the new values take effect the next time you generate a new default solver.
Pseudo Time Stepping for Laminar Flow Models in this guide and Pseudo Time Stepping in the COMSOL Multiphysics Reference Manual
Flow in an Airlift Loop Reactor: Application Library path CFD_Module/Verification_Examples/airlift_loop_reactor