1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
|
#import "@preview/charged-ieee:0.1.3": ieee
#show: ieee.with(
title: [Stokes flow - Particle interaction in low Reynolds Number environments],
abstract: [
],
authors: (
(
name: "Claudius Holeksa",
department: [CSSR],
organization: [NORCE Research AS],
location: [Bergen, Norway],
email: "clho@norceresearch.no"
),
),
index-terms: ("Fluid Mechanics", "Review"),
bibliography: bibliography("refs.bib"),
figure-supplement: [Fig.],
)
= Introduction
Flow in porous subsurface structures often is dominated by extremely low velocity and high viscous fluids with low Reynolds numbers $"Re" << 1$.
In such a condition, the flow enters a state called the Stokes Flow or creeping flow, where inertial forces can be assumed to be zero and viscous
forces dominate the whole connected set.
This kind of flow is found in geological processes, in carbohydrate recovery, filtration systems and CO_2 storage where the pore-scale motion
of particles and their interaction with the surrounding porous medium strongly influence macroscopic properties such as the permeability of
rock layers.
Undestanding particle transport, bridging and in general clogging phenomena is critical for the effective use and management of subsurface
reservoirs.
The tiniest changes can alter the established path in a porous geometry and transport behaviour over time.
Thus, the behaviour of individual particles in Stokes Flow form the
For the understanding of porous media, where particle transport, bridging and clogging phenomena are critical to the use of subsurface
reservoirs. Be it either for carbohydrate resource recovery or storage of CO2. Understanding the behaviour of particles in Stokes flow
is key to analyzing more complex phenomena emerging over time such as permeability reduction over periods of time.\
The study of the Stokes Flow has a long history, beginning with Stokes' solution flow along a sphere@stokes1851.\
Subsequent work expanded to explore a multitude of shapes and the interaction of multiple particles.
More recently, research was performed in the understanding of creeping flow which are essential in water treatment, carbohydrate recovery
and CO2 storage. // CITE
This report delves into the central aspects governing the known behaviour past single and multiple spheres.
Explores the effect of it in porous media and examines the more recent particle bridging behaviour found in porous structures.
== Stokes Flow
// Write about Navier
The incompressible Navier Stokes equations are a set of equations defined with the density $rho$, the dynamic viscosity $mu$ and the
velocity field $u(x,t)$
#math.equation(
block: true,
$ ρ((∂)/(∂t)u + (u · ∇)u) = −∇p + μ∇²u + f,\
∇ · u = 0 $
)
While we are often interested in Navier-Stokes flows, a generalized analytical solution does not exist,
most analytical approaches rely on a set of restrictions or assumptions.
With the Reynolds number defined by
$ "Re" = (rho U L ) / mu
$
Luckily we decided to work on a fluid which tends $"Re" arrow 0$ or at least $"Re" << 1$ where due to the viscous forces
the inertial term is negligable and is assumed to be zero.
Which is why we arrive at this term
// Stokes
#math.equation(
block: true,
$ μ∇²u + f = ∇p,\
∇ · u = 0 $
)
An additional restriction by the creeping flow equations is the selection of boundary criteria where it depends on
an infinite space with a set solid boundary interface.
For simplicity a no-slip boundary condition is selected for the solid boundary.
#math.equation(
block: true,
$ v = U "at particles boundary",\
v arrow 0 "as" |x| arrow infinity
$
)
== Solid-Fluid single sphere interaction with Stokes flow
With the previous assumptions a solid sphere of radius $r$ moving with a relative velocity $v$ in an unbounded creeping flow,
we receive the drag force provided by Stokes
#math.equation(
block: true,
$ F_S = 6 pi mu r v $
)
with our boundary conditions specified as
#math.equation(
block: true,
$ u = v "at" a = r,\
u arrow 0 "as" a arrow infinity
$
)
where according to Proudman and Pearson@proudman_pearson_1957 at a Reynolds number of 0.05 the predicted drag is two percent lower
than the possibly more correct value by Proudman and Pearson@proudman_pearson_1957.
//#math.equation(
//block: true,
//$ F = F_S ( 1 + (3/8)"Re" + (9/40)("Re")^2 * (log "Re" + gamma + (5/3) log 2 - (323/360)) + (27/80)*"Re"^3 log "Re" ) $
//)
//Using spherical coordinates $(r,theta,phi.alt)$ and no external force we receive the velocity components
//#math.equation(
//block: true,
//$ v_r = U cos theta (1- (3r)/(2a) + r^3 / (2a^3)) $
//$ v_theta = - U cos theta (1- (3r)/(4a) + r^3 / (4a^3)) $
//)
== Lubrication Forces for Near Contact
When spheres or a sphere near a straight wall approach each other with a small separation $h<<r$ the flow in the gap is dominated by
lubrication forces. For a sphere approaching a straight wall with velocity $U$ normal to the wall we
have
#math.equation(
block: true,
$ F_(L) tilde.eq (6 pi mu r^2 U) / h (1-(9h)/(16r) + ...) $
)
With two spheres of equal radius r and otherwise the same conditions approaching
#math.equation(
block: true,
$ F_(L) tilde.eq (6 pi mu r^2 U) / h (1-(h)/(5r) + ...) $
)
These forces diverge the smaller $h$ becomes.
This is due to the viscous stress increasing on small films of
fluid which then in turn dominate the interaction.
== Clogging of porous structures
Modifying a poiseulle channel flow by narrowing the channel we can observe mainly three different ways how the porosity is reduced.
By sieving, bridging and aggregation of surface deposition.
The least interesting phenomenon is sieving which results purely from a single particle having a too large of a volume to fit into the
pore. While bridiging occurs when multiple particles arrive at a passage at the exact time, forming a stable arch which blocks the flow.
Aggregation involves the deposition of particles gradually narrowing the channel.
And while sieving and bridging results from purely mechanical forces aggregation results from electrostatic interactions between particles, the fluid and surface boundaries at micrometer scales@laurez2025bridging.
== Discusion
The arches resulting from clogging are often inherently unstable, but with the progression of aggregation arrive at less unstable states.
Which in turn makes a stochastic analysis of particles in different conditions@laurez2025bridging based on different widths feasible.
Laurez@laurez2025bridging goes into detail based on CITE FIG how probable clogging mechanisms are, while these are quite indepth, opportunities
in relation to multiphase fluid interaction with the solid structures.
An interesting case
// Insert png
== Lattice-Boltzmann-Method
=== Immersed Boundary Method
=== Homogenized Lattice-Boltzmann
|
