Code_Saturne boundary conditions
Forum rules
Please read the forum usage recommendations before posting.
Please read the forum usage recommendations before posting.
-
WDIAROM
Code_Saturne boundary conditions
How can I impose a Robin boundary condition in Saturne?
Last edited by WDIAROM on Fri Mar 20, 2015 12:52 pm, edited 1 time in total.
-
Yvan Fournier
- Posts: 4085
- Joined: Mon Feb 20, 2012 3:25 pm
Re: Code_Saturne boundary conditions
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
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