Contaminant-Removal From Wastewater in A Secondary Clarifier by Sedimentation
Contaminant-Removal From Wastewater in A Secondary Clarifier by Sedimentation
Contaminant-Removal From Wastewater in A Secondary Clarifier by Sedimentation
This model is licensed under the COMSOL Software License Agreement 5.6.
All trademarks are the property of their respective owners. See www.comsol.com/trademarks.
Introduction
Wastewater treatment is a several-step process for removing contaminants. First, large,
solid particles are removed through sedimentation, flotation, and filtration. In a second
step, biological treatment causes the smaller particles to aggregate, forming so-called flocs.
These flocs can more easily be removed, for instance through sedimentation. The present
example studies the separation of flocs from water in a circular secondary clarifier. To
model the turbulent multiphase flow in the tank, this example uses the Mixture Model, k-
ε interface.
Model Definition
The model geometry is shown in Figure 1.
The clarifier has a diameter of 24 m and a gently sloping bottom, which makes the dept
vary between 3.3 and 4 m, attached to a funnel at the center of the tank. The incoming
sludge, consisting of a mixture of solid flocs and water, enters through the inlet in the
middle of the tank. In the clarifier, gravity causes the flocs to settle to the bottom of the
tank. The tank contains two outlets. One is located at the center, at the bottom of the
funnel. The purpose of this outlet is to remove the sedimented flocs from the tank. There
The mixture enters the clarifier in the form of a jet. The Reynolds number based on the
inlet velocity and diameter is 5·105, which means that the flow is turbulent. Upon impact
with the free surface, the mixture spreads out, causing the turbulence production to
diminish with radial distance from the inlet. The turbulent flow in the tank tends to mix
the phases together, and thus has a negative effect on the separation process. The aim of
this example is to study the turbulent multiphase flow within the circular secondary
clarifier.
For simplicity, you can model the flocs as spherical solid particles of equal size. To solve
for the mixture velocity, pressure and the phase volume fractions, you can use the Mixture
Model, k-ε interface. The mixture model is a multiphase-flow model, particularly well
suited for suspensions of solid particles in liquid at low particle volume fractions. For the
slip velocity, you can use the Hadamard-Rybczynski drag law. See the theory in the CFD
Module User’s Guide for more information.
BOUNDARY CONDITIONS
At the inlet, the velocity is fixed to 1.25 m/s and the dispersed phase volume fraction is
0.003. The turbulence intensity and length scale are set to 5% and 0.14rin by default (see
Inlet Values for the Turbulence Length Scale and Turbulent Intensity in the CFD
Module User’s Guide), where rin = 0.2 m is the radius of the inlet. The velocity at the
bottom outlet is set to 0.05 m/s, while a constant pressure is set at the peripheral outlet.
INITIAL CONDITIONS
Initially, the velocity as well as the solid phase volume fraction is zero in the entire clarifier.
Figure 3: Mixture-velocity streamlines and solid phase volume fraction after 12 hours, when
the flow has reached a steady state solution.
As you can see in Figure 3, the maximum volume fraction is less than 1 %. Hence, the
Mixture Model should be sufficiently accurate. See Two-Phase Flow Modeling of a Dense
Figure 4 shows the dispersed phase mass flux at the inlet and the two outlets.
Figure 4: Mass flux of the dispersed phase at the inlet (blue), peripheral outlet (green) and
central outlet (red).
Md = ρd jd, eff ⋅ n dS
where jd,eff is the effective dispersed phase flux
j d, eff = φ d ( j + ( 1 – φ d )u slip ) – D md ∇φ d
Computing the removal rate from the results shows that the clarifier removes 0.52 − 0.10
= 0.42 kg solid particles per second. The separation efficiency is about 81%.
Figure 5: Volume fraction of the dispersed phase and streamlines for the dispersed (black) and
continuous (white) velocity fields.
To further examine the performance of the clarifier, you can easily modify the model in
several ways. You can, for instance, modify the geometry by adding baffles, changing the
inlet and outlet velocities, increasing the dispersed-phase volume fraction in the incoming
sludge, or changing the density and size of the dispersed particles.
NEW
In the New window, click Model Wizard.
MODEL WIZARD
1 In the Model Wizard window, click 2D Axisymmetric.
2 In the Select Physics tree, select Fluid Flow>Multiphase Flow>Mixture Model>
Mixture Model, Turbulent Flow>Mixture Model, k-ε (mm).
3 Click Add.
4 Click Study.
5 In the Select Study tree, select General Studies>Time Dependent.
6 Click Done.
GEOMETRY 1
To simplify the instructions we import a ready made geometry sequence for this model.
The complete step-by-step instructions for the geometry can be found in the Appendix —
Geometry Modeling Instructions.
Use Mesh Control Edges to obtain a mesh which is aligned with the turbulent shear.
GLOBAL DEFINITIONS
Parameters 1
1 In the Model Builder window, under Global Definitions click Parameters 1.
2 In the Settings window for Parameters, locate the Parameters section.
Add a Step function feature to use for implementing a gradual increase of the inlet and
outlet velocities from zero to their constant values.
DEFINITIONS
Step 1 (step1)
1 In the Home toolbar, click Functions and choose Local>Step.
2 In the Settings window for Step, click to expand the Smoothing section.
3 Locate the Parameters section. In the Location text field, type 0.1.
4 Locate the Smoothing section. In the Size of transition zone text field, type 0.2.
5 Click Plot.
Variables 1
1 In the Home toolbar, click Variables and choose Local Variables.
Mixture Properties 1
1 In the Model Builder window, under Component 1 (comp1)>Mixture Model, k-ε (mm) click
Mixture Properties 1.
2 In the Settings window for Mixture Properties, locate the Continuous Phase Properties
section.
3 From the ρc list, choose User defined. In the associated text field, type rho_c.
4 From the μc list, choose User defined. In the associated text field, type mu_c.
5 Locate the Dispersed Phase Properties section. From the ρd list, choose User defined. In
the associated text field, type rho_d.
6 In the dd text field, type d_d.
7 In the Model Builder window, click Mixture Model, k-ε (mm).
8 In the Settings window for Mixture Model, k-ε, locate the Physical Model section.
9 From the Slip model list, choose Hadamard-Rybczynski.
Gravity 1
1 In the Physics toolbar, click Domains and choose Gravity.
2 Select Domain 1 only.
Wall 2
1 In the Physics toolbar, click Boundaries and choose Wall.
2 Select Boundary 7 only.
3 In the Settings window for Wall, locate the Mixture Boundary Condition section.
4 From the Mixture boundary condition list, choose Slip.
Outlet 1
1 In the Physics toolbar, click Boundaries and choose Outlet.
2 Select Boundary 2 only.
3 In the Settings window for Outlet, locate the Mixture Boundary Condition section.
4 From the list, choose Velocity.
5 Locate the Velocity section. In the J0 text field, type v_out.
Outlet 2
1 In the Physics toolbar, click Boundaries and choose Outlet.
2 Select Boundary 18 only.
Initial Values 1
1 In the Model Builder window, click Initial Values 1.
2 In the Settings window for Initial Values, locate the Initial Values section.
3 In the p text field, type -g_const*z*rho_c.
MESH 1
1 In the Model Builder window, under Component 1 (comp1) click Mesh 1.
2 In the Settings window for Mesh, locate the Physics-Controlled Mesh section.
3 From the Element size list, choose Extra fine.
4 Locate the Mesh Settings section. From the Sequence type list, choose User-
controlled mesh.
Size 1
In the Model Builder window, under Component 1 (comp1)>Mesh 1 right-click Size 1 and
choose Build Selected.
Corner Refinement 1
In the Model Builder window, right-click Corner Refinement 1 and choose Disable.
Boundary Layers 1
1 In the Model Builder window, expand the Component 1 (comp1)>Mesh 1>
Boundary Layers 1 node, then click Boundary Layers 1.
2 In the Settings window for Boundary Layers, click to expand the Corner Settings section.
3 From the Handling of sharp corners list, choose Splitting.
STUDY 1
4 Click to expand the Results While Solving section. Select the Plot check box.
Solution 1 (sol1)
1 In the Study toolbar, click Show Default Solver.
Specify an appropriate initial time step.
RESULTS
2D Plot Group 1
In the Home toolbar, click Add Plot Group and choose 2D Plot Group.
Surface 1
1 Right-click 2D Plot Group 1 and choose Surface.
2 In the Settings window for Surface, click Replace Expression in the upper-right corner of
the Expression section. From the menu, choose Component 1 (comp1)>Mixture Model, k-
ε>phid - Volume fraction, dispersed phase.
3 Click to expand the Range section. Select the Manual color range check box.
4 In the Maximum text field, type 0.006.
5 Click the Zoom Extents button in the Graphics toolbar.
2D Plot Group 1
1 In the Model Builder window, click 2D Plot Group 1.
2 In the Settings window for 2D Plot Group, locate the Plot Settings section.
3 Click Go to Source.
STUDY 1
In the Home toolbar, click Compute.
RESULTS
2D Plot Group 1
Calculate the inflow and outflow rates of the dispersed phase. Start with the inflow rate.
Evaluation Group 1
In the Results toolbar, click Evaluation Group.
Line Integration 1
1 Right-click Evaluation Group 1 and choose Integration>Line Integration.
2 Select Boundary 5 only.
3 In the Settings window for Line Integration, locate the Expressions section.
4 In the table, enter the following settings:
5 Locate the Integration Settings section. Clear the Compute surface integral check box.
The negative sign is used since qd_out is defined as the outward flux.
Line Integration 2
1 In the Model Builder window, right-click Evaluation Group 1 and choose Integration>
Line Integration.
2 Select Boundary 18 only.
3 In the Settings window for Line Integration, locate the Expressions section.
5 Locate the Integration Settings section. Clear the Compute surface integral check box.
Line Integration 3
1 Right-click Evaluation Group 1 and choose Integration>Line Integration.
2 Select Boundary 2 only.
3 In the Settings window for Line Integration, locate the Expressions section.
4 In the table, enter the following settings:
5 Locate the Integration Settings section. Clear the Compute surface integral check box.
6 In the Evaluation Group 1 toolbar, click Evaluate.
TABLE
1 Go to the Table window.
2 Click Table Graph in the window toolbar.
RESULTS
1D Plot Group 2
1 In the Model Builder window, under Results click 1D Plot Group 2.
2 In the 1D Plot Group 2 toolbar, click Plot.
3 In the Settings window for 1D Plot Group, locate the Plot Settings section.
4 Select the y-axis label check box.
5 In the associated text field, type Mass flow (kg/s).
Table Graph 1
1 In the Model Builder window, click Table Graph 1.
2 In the Settings window for Table Graph, click to expand the Legends section.
3 Select the Show legends check box.
3D Plot Group 3
In the Results toolbar, click 3D Plot Group.
Surface 1
1 Right-click 3D Plot Group 3 and choose Surface.
2 In the Settings window for Surface, click to expand the Range section.
3 Select the Manual color range check box.
4 In the Maximum text field, type 6e-3.
3D Plot Group 3
Choose a variable defined on the walls, such as δw+ to visualize the surrounding wall.
Surface 2
1 In the Model Builder window, right-click 3D Plot Group 3 and choose Surface.
2 In the Settings window for Surface, click Replace Expression in the upper-right corner of
the Expression section. From the menu, choose Component 1 (comp1)>Mixture Model, k-
ε>Turbulence variables>mm.Delta_wPlus - Wall resolution in viscous units.
3 Locate the Coloring and Style section. From the Coloring list, choose Uniform.
4 From the Color list, choose Gray.
5 Locate the Expression section. Select the Description check box.
6 Clear the associated text field.
7 In the 3D Plot Group 3 toolbar, click Plot.
8 Click to expand the Title section. From the Title type list, choose None.
9 In the 3D Plot Group 3 toolbar, click Plot.
Streamline 1
1 Right-click 3D Plot Group 3 and choose Streamline.
2 In the Settings window for Streamline, locate the Expression section.
3 In the r component text field, type mm.udr.
Streamline 2
1 Right-click 3D Plot Group 3 and choose Streamline.
2 In the Settings window for Streamline, locate the Expression section.
3 In the r component text field, type mm.ucEffr.
4 In the phi component text field, type 0.
5 In the z component text field, type mm.ucEffz.
6 Select the Description check box.
7 In the associated text field, type Continuous phase (white).
8 Locate the Streamline Positioning section. From the Entry method list, choose
Coordinates.
9 In the x text field, type range(0,0.02,0.2) range(0.5,0.5,12).
10 In the y text field, type 0.
11 In the z text field, type -1^range(1,35).
12 Locate the Coloring and Style section. Find the Point style subsection. From the Color
list, choose White.
13 In the 3D Plot Group 3 toolbar, click Plot.
14 Click the Zoom In button in the Graphics toolbar.
GEOMETRY 1
Polygon 1 (pol1)
1 In the Geometry toolbar, click Polygon.
2 In the Settings window for Polygon, locate the Coordinates section.
3 In the table, enter the following settings:
r (m) z (m)
0 0
12 0
12 -3.3
2 -4
0.5 -7
0.5 -7.4
0 -7.4
0 0
Polygon 2 (pol2)
1 In the Geometry toolbar, click Polygon.
2 In the Settings window for Polygon, locate the Object Type section.
3 From the Type list, choose Open curve.
4 Locate the Coordinates section. From the Data source list, choose Vectors.
5 In the r text field, type 0 0.4 0.4 0.4.
6 In the z text field, type -5.4 -5.4 -5.4 -3.4.
Polygon 3 (pol3)
1 In the Geometry toolbar, click Polygon.
2 In the Settings window for Polygon, locate the Object Type section.
3 From the Type list, choose Open curve.
4 Locate the Coordinates section. From the Data source list, choose Vectors.
5 In the r text field, type 0.2 0.2 0.2 0 0 0.
6 In the z text field, type -3.4 -5.2 -5.2 -5.2 -5.2 -5.4.
Rectangle 1 (r1)
1 In the Geometry toolbar, click Rectangle.
2 In the Settings window for Rectangle, locate the Size and Shape section.
3 In the Width text field, type 0.4.
4 In the Height text field, type 0.4.
5 Locate the Position section. In the r text field, type 11.2.
6 In the z text field, type -0.4.
7 Click Build Selected.
Difference 1 (dif1)
1 In the Geometry toolbar, click Booleans and Partitions and choose Difference.
2 Select the object pol1 only to add it to the Objects to add list.
3 In the Settings window for Difference, locate the Difference section.
4 Find the Objects to subtract subsection. Select the Activate Selection toggle button.
5 Select the objects c1, csol1, and r1 only.
Rectangle 2 (r2)
1 In the Geometry toolbar, click Rectangle.
2 In the Settings window for Rectangle, locate the Size and Shape section.
3 In the Height text field, type 0.5.
4 Locate the Position section. In the r text field, type 11.6.
Fillet 1 (fil1)
1 In the Geometry toolbar, click Fillet.
2 On the object r2, select Points 3 and 4 only.
3 In the Settings window for Fillet, locate the Radius section.
4 In the Radius text field, type 0.5.
Difference 2 (dif2)
1 In the Geometry toolbar, click Booleans and Partitions and choose Difference.
2 Select the objects dif1 and fil1 only.
3 In the Settings window for Difference, locate the Difference section.
4 Find the Objects to subtract subsection. Select the Activate Selection toggle button.
5 Select the object c2 only.
6 Clear the Keep interior boundaries check box.
7 Click Build Selected.
Add a few lines to obtain a better control of the mesh in the turbulent shear regions.
Polygon 4 (pol4)
1 In the Geometry toolbar, click Polygon.
2 In the Settings window for Polygon, locate the Coordinates section.
3 In the table, enter the following settings:
r (m) z (m)
0.19 -3.22
0.2 -2.9
Polygon 5 (pol5)
1 In the Geometry toolbar, click Polygon.
2 In the Settings window for Polygon, locate the Coordinates section.
3 In the table, enter the following settings:
r (m) z (m)
0.2 -2.9
0.35 -0.25
r (m) z (m)
0.2 -0.1
7.6 -0.6