Hello to all,

I am trying to solve the Non-Linear Reynolds equation on a rectangular domain with Dirichlet (=-pps) pressure all around.

I create the geometry and the coefficient function for h but have no experience with ngsolve and cant figure hot to do the discretization and initialization properly. Is anyone eager to help me out?

Thank you,

George K

Reynolds.ipynb (12.9 KB)

Hi,

after doing integration by parts your bf term should look something like this:

```
a += (h**3/(6*nu) * u * grad(u) * grad(v) - (p * h * CF((U,V)) * grad(v)) * dx
```

for newton convergence it’s better to set the dirichlet values already on the whole domain and not only on the boundary

Thank you very much Christopher.

I’m one step closer I just have to figure out periodic boundaries now.

you already had the identification correctly in the script. Just wrap the H1 space with a Periodic(H1(…)) and remove the dirichlet from left and right

Hello, I have a follow-up question.

Is it possible to impose a constraint within the iterative process so my result does not go beneath a certain value (eg ambient pressure). This is for cavitation modeling in Journal bearings.

Thank you!