Hi team,
I did a convergence study using Code_Saturne v7.0.2. I tested a pseudo-3D Poiseuille flow with inlet boundary (u=4.0*(y-y*y)) and outlet boundary. Also, the domain is (1m * 1m) and uniform mesh with 16, 50, 100 resolutions. In order to ensure accuracy, I set SOLU velocity scheme. However, my results shows about 1st accuracy. Do you think this is reasonable for this case? Here I attached convergence results and setup file.
[img][PureCSaccuracy.jpg]
Best wishes,
Yang
Convergence study in Code_Saturne
Forum rules
Please read the forum usage recommendations before posting.
Please read the forum usage recommendations before posting.
Convergence study in Code_Saturne
- Attachments
-
- setup.xml
- (7.62 KiB) Downloaded 463 times
-
- Posts: 4145
- Joined: Mon Feb 20, 2012 3:25 pm
Re: Convergence study in Code_Saturne
Hello,
In most of our verification cases, we obtain 2nd order convergence, so your results are a bit surprising.
Did you test using the default centered scheme ? SOLU is somewhat more robust, not more precise.
Also, what type of mesh are you using ?
Regards,
Yvan
In most of our verification cases, we obtain 2nd order convergence, so your results are a bit surprising.
Did you test using the default centered scheme ? SOLU is somewhat more robust, not more precise.
Also, what type of mesh are you using ?
Regards,
Yvan
Re: Convergence study in Code_Saturne
Thanks Yvan,
(1) Did you test using the default centered scheme?
I think Code_Saturne default is velocity upwind scheme, and it results is first-order accuracy. Let me test centered scheme then.
(2)what type of mesh are you using ?
I generate mesh using Code_Saturne GUI, as shown below.
Best,
Yang
(1) Did you test using the default centered scheme?
I think Code_Saturne default is velocity upwind scheme, and it results is first-order accuracy. Let me test centered scheme then.
(2)what type of mesh are you using ?
I generate mesh using Code_Saturne GUI, as shown below.
Best,
Yang
Re: Convergence study in Code_Saturne
Hi Yvan,
I tested this case using centered velocity scheme, but it shows similar results.
Best,
Yang
I tested this case using centered velocity scheme, but it shows similar results.
Best,
Yang
-
- Posts: 4145
- Joined: Mon Feb 20, 2012 3:25 pm
Re: Convergence study in Code_Saturne
Hello,
Other important question, how do you measure the error ? Which type of data do you extract ? Do you interpolate it in any way ?
Regards,
Yvan
Other important question, how do you measure the error ? Which type of data do you extract ? Do you interpolate it in any way ?
Regards,
Yvan
Re: Convergence study in Code_Saturne
Hi Yvan,
I used paraview to extract data (velocity distribution) and tested interpolation and no interpolation. The error calculation is through L2-norm:
Best,
Yang
I used paraview to extract data (velocity distribution) and tested interpolation and no interpolation. The error calculation is through L2-norm:
Best,
Yang
Re: Convergence study in Code_Saturne
Hello Yvan,
Thanks for your advice in post-processing. I know the reasons and list here in case someone has similar questions:
(1) I should extract velocity component instead of resultant velocity.
(2)I extract data from the cell centre and avoid any interpolation scheme.
(3)I use velocity centered scheme in Code_Saturne.
Below shows the good results:
Thanks!
Best,
Yang
Thanks for your advice in post-processing. I know the reasons and list here in case someone has similar questions:
(1) I should extract velocity component instead of resultant velocity.
(2)I extract data from the cell centre and avoid any interpolation scheme.
(3)I use velocity centered scheme in Code_Saturne.
Below shows the good results:
Thanks!
Best,
Yang
-
- Posts: 4145
- Joined: Mon Feb 20, 2012 3:25 pm
Re: Convergence study in Code_Saturne
Hello,
That is reassuring.
Now that your postprocessing does not interfere with the convergence order measurement any more, I gues you should also get something close to 2nd order with SOLU, though I do not expect it to be more precise than the centered option (though it may be more robust on distorted meshes, which can be useful in industrial studies).
Best regards,
Yvan
That is reassuring.
Now that your postprocessing does not interfere with the convergence order measurement any more, I gues you should also get something close to 2nd order with SOLU, though I do not expect it to be more precise than the centered option (though it may be more robust on distorted meshes, which can be useful in industrial studies).
Best regards,
Yvan
Re: Convergence study in Code_Saturne
Hi Yvan,
I have tested Code_Saturne three velocity scheme -- automatic, SOLU and second-order scheme.
As you predicted, SOLU and second-order scheme show excellent 2nd-order accuracy, but second-order scheme has lower error. The default setup (automatic) only has first-order accuracy, as shown in below figure.
There is an upwind scheme I didn't test (1st-order accuracy?). If I set automatic, Code_Saturne will adopt which scheme?
Correction: In previous conversation, I showed SOLU scheme has slope -2.53. This result is error because wrong boundary treatment in post-processing.
Best wishes,
Yang
I have tested Code_Saturne three velocity scheme -- automatic, SOLU and second-order scheme.
As you predicted, SOLU and second-order scheme show excellent 2nd-order accuracy, but second-order scheme has lower error. The default setup (automatic) only has first-order accuracy, as shown in below figure.
There is an upwind scheme I didn't test (1st-order accuracy?). If I set automatic, Code_Saturne will adopt which scheme?
Correction: In previous conversation, I showed SOLU scheme has slope -2.53. This result is error because wrong boundary treatment in post-processing.
Best wishes,
Yang
-
- Posts: 4145
- Joined: Mon Feb 20, 2012 3:25 pm
Re: Convergence study in Code_Saturne
Hello Yang,
The default is centered. Or more precisely, centered with slope test. When the slope test kicks in for a given face, it locally switches to pure upwind (for that face and operation), so the convergence order drops to 1. (There is an advanced setting allowing to use blending here, but it is still work in progress).
Best regards,
Yvan
The default is centered. Or more precisely, centered with slope test. When the slope test kicks in for a given face, it locally switches to pure upwind (for that face and operation), so the convergence order drops to 1. (There is an advanced setting allowing to use blending here, but it is still work in progress).
Best regards,
Yvan