J. Szantyr – Lecture No. 11 – Computations of Viscous Flows – Finite Element Method
The analysed flow domain is divided into many small parts, so called finite elements. In the selected points of the elements (nodes) the
values of the desired function, e.g. velocity, pressure etc. are to be determined. The solution is sought in the form of a basis function
approximating the solution. Parameters of the approximating function are determined by means of the Variational Calculus.
Example of the one-dimensional solution – 2D Poiseuille’ flow.
The equation describing the flow:
2 0
2
dx d z
L p
If is the approximate solution, then in general: and the measure of the error is the weighed residual:
~
~ 0L
h h
dx L
Wdx L
R
0 0
~
~
~
The condition for minimization of the functional R has the form:
h h
dx dx d dx d z
dx p L
R
0 0
0
~
~ ~
~
~
This can be written as:
h ddx ddx dx pz
h dx0 0
~ ~
~
The interval of integration [0,h] is divided into 4 equal sections
(finite elements) having lengths of l=h/4. Their endpoints determine 5 nodes. In each of them there is the nodal value of velocity:
1
2,3,4
l In the first and last node
there are boundary values:
5 0
1
The sought solution is approximated by the local linear basis functions, defined within each finite element:
x
x
~ 11 where:
x x l
After substituting the local basis functions the condition for minimalization of the functional R takes the form:
4
1 1
0
4
1 1
0
~ ~
~
x z d
l p x
x d d d x d d
l
which may be developed into:
1
0
4 3
4 2
3 2
4 4
3 4
3 4
1
0
2 3
2 3
2 2
1 1
1
1 1
1 1
1 1
x x
x x
x z x
l p
dx x x
d x d x
d x d
x x d x d x x
d d
x x x
d x d dx x
x d x
d x d x
d d l
0 , 1
x
After integration and re-formulation it leads to the system of linear equations for the unknown nodal values of velocity:
where:
z c p
Solution of the system produces the following results:
2 4 232
3 h
z p
3
232
4 h
z p
2 3 2 2 3 4 3 3 4 4
2 3 42 2v v v v 2v v v v 2v v v v v
cl
The analytical solution: steady laminar flow between two parallel infinite flat plates (Poiseuille flow)
Given: const
x p
Boundary conditions:
u=0 for y=h u=0 for y=-h
The Navier-Stokes equation takes the form:
x p dy
u d
1
2 2
After double integration we obtain:
2 1 22
1 y C y C
x y p
u
The integration constants are determined from the boundary condtions:
2 2
2
1 h y
x y p
u
Transformation to the notation adopted in the Finite Element solution leads to:
2 x H y
2
h H
2 22
2 H
x H
x z
p x H
2 2 2 2
32 3 16
4 2
4
H z p H
H H
H z
p H
H
2 22 2
32 3 16
9 4
3 2
4
3 H
z p H
H H
H z
p H
H
2 22 2
32 4 4
2 2
2
H z p H
H H
H z
p H
H
Hence the analytical
solution is exactly equal to the numerical solution using FEM, but only in the nodes.