Gtheta Computations using UDF
A User Defined FLUENT script (UDF) was created to extract the post processing result from individual cell of the model.
There are four methods used in the calculation of the Gtheta values.
The first method is to extract the individual G and theta values from the cell and sum up all the cells to obtain the total Gtheta. The total Gtheta is then divided by the total cells to get the average Gtheta.
∑(Gcell*Өcell)/N (1)
The second method was to extract individual G value and take the average of the G value. The theta is calculated by volume divided by inlet flow rate.
Өfloc = Vol./(Vol. flow rate)
∑(Gcell)/N* Өfloc (2)
The third method is multiple the sum of the G-value and the sum of the theta value and divide by the number of cells.
∑(Gcell)/N* ∑(Өfloc)/N (3)
Method | Gtheta Value for k-e realizable 1_5 W one baffle |
|---|---|
Ave of Cell's Product (1) | .6698 |
Ave G* Floc residence time (2) | 45 |
Ave G*Ave N (3) | .617 |
The fourth method is to determine Π-value of the energy dissapation region and pressure drop of the flow:
sqrt(Π-cell*K*b*V/(2ν)) (4)
The parameters for the K-e realizable 1_5 Width one baffle geometry for one turn are shown below:
Π-cell: the ratio of the length of the dissipation zone divided by the baffle spacing (approximately 6)
K: Pressure drop coefficient (approximately 4)
b: baffle width (.1 m)
V: Velocity (.1)
ν: kinematic viscosity of water (approx 1e-6 in region of interest)
Method | Gtheta Value for k-e realizable 1_5 W one baffle |
|---|---|
π-value, k (Monroe) | 350 |
| Wiki Markup |
|---|
Gtheta was also plotted using the custom function in FLUENT. (GӨ = Const. \[ε/υ\]^1/2^ \* \[Cell Volume/Vol. flow rate\] )
\\ !aGtheta plot.png! |
G_theta value from the flocculator. The UDF loops through all of the individual cells of the converged solution, and extracts the energy dissipation^1/2 times the cell area. The sum of this quantity for all cells divided by flow and square of viscoisty results in the G_theta value of the entire flocculator as shown in equation 1.
Gθ_baffle=1/Q*vis^(-1/2)*∑(ε^(1/2)*cell_area) Eq. 1
The above formulation has been determined (with Monroe's help) by weighting the value of Gθ by volume.
The Gθ_baffle values weighted by the area or number of cells were also examined, but these formulations do not produce sensible results.
Below is the G_theta value presented for flocculation tanks of different geometries:
Case Geometry (h-height,b-width,N-baffles) | Gθ_Flocculator |
|---|---|
b=.1, h=1, N=1, Clearance height=.15 | 430 |
b=.1, h=.3, N=5, sym bc | 1330 |
b=.1, h=.2, N=5 | 1060 |
...

