Changes
Page history
boundary conditions
authored
Aug 11, 2019
by
Praetorius, Simon
Show whitespace changes
Inline
Side-by-side
examples/navier-stokes.md
View page @
4aef7ea8
...
...
@@ -155,17 +155,13 @@ prob.addVectorOperator(opForce, _v);
Numerical example
-----------------
We consider a square domain $
`\Omega`
$ with $
`L_x=L_y=1`
$ and Dirichlet boundary conditions
We consider a lid-driven cavity problem in a square domain $
`\Omega`
$ with $
`L_x=L_y=1`
$
no external force, and with Dirichlet boundary conditions
```
math
\mathbf{g} = \begin{pmatrix}0 \\ x_1(1 - x_1)(1 - x_0)\end{pmatrix}
```
Thus, the boundary value is zero on top, right, and bottom boundary and parabolic
on the left boundary.
This is the lid-driven cavity problem, where we do not include any other external force
term.
The boundary condition is implemented, by first setting boundary ids for parts of
on the left boundary. The boundary condition is implemented, by first setting boundary ids for parts of
the grid boundary and second, setting the Dirichlet value:
```
c++
double
vel
=
1.0
;
...
...
...
...