Julia Mikhal
Bernard Geurts
Immersed Boundary
Method for
Pulsatile Flow in
Cerebral Aneurysms
−3
8
x 10
3.4
1
3.2
0.8
7
0.6
3
6
0.4
2.8
5
0.2
2.6
4
0
2.4
−0.2
3
2.2
−0.4
2
2
−0.6
1.8
−0.8
1.6
0
10
20
30
40
50
60
70
80
90
100
1
−1
0
2
4
6
8
10
12
14
16
018
0
20
12
2
4
6
8
10
12
10
8
Incompressible Navier-Stokes
6
∇·u=0
∂u
1 2
+ u · ∇u = −∇p +
∇ u+f
∂t
Re
Volume penalization
1
f = − H(x)u
ε
H(x) = 1
=0
2
ε!1
Masking Function in
Ω = Ωs ∪ Ωf
x ∈ Ωs
x ∈ Ωf
Property u ≈ 0
4
if H(x) = 1
Skew-symmetric discretization
0
0
2
4
6
8
10
12