boundary conditions authored by Praetorius, Simon's avatar Praetorius, Simon
......@@ -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;
......
......