interCondensatingEvaporatingFoam for closed domain
The solver shows a very strange behaviour for
alpha field for closed domain cases. Simply after running the simulation for a while the entire domain reach
alpha = 1. Indeed closed domains pose numerical challenges but usually the problem appears due to
p underspecification not
alpha. I never had a similar experience with
interFoam in closed domain cases.
To reproduce the issue:
- use the existing tutorial for this solver,
- set the
toppatch to wall.
- use the same boundary conditions for all fields as the
- add reference value for pressure in
Note: I am using OpenFOAM-plus updated to commit e6cbe0b1 to make sure I didn't miss any bug fix.
Thanks for your support, please let me know if further testes are required.