Code_Saturne boundary conditions
Posted: Mon Nov 26, 2012 7:10 pm
How can I impose a Robin boundary condition in Saturne?
Hello,
I'll let Martin or others give you more details, but the boundary conditions have changed quite a bit in versions 2.3 and 3.0: we now save coefa/coefb coefficients separatetly for gradient and flux boundary conditions, so handling of those is more systematic and consistent (in previous versions, quite a few extra/intermediate computations were done, to save a bit of memory, but this led to more tricky code).
So I would definitely recommend moving to version 3.0-beta, for this reason in addition to all other improvements.
If you want to look at the code, most boundary condition computations are done in condli.f90 and clptur.f90. Applying boundary conditions to a variable's derivatives may be a more complex issue, though.
Best regards,
Yvan Fournier