div(phi,U) GammaV ?
Hello,
i read some parts about the Gamma Scheme in Jasak's Thesis (pages 107-110) and wanted to try it now for my openFOAM simulations. My question is now, what the Greek_Gamma value should be if i define for example: div (phi,U) GammaV [Greek_Gamma] In the Documentation there is a recommendation for Greek_Gamma = 1 ( i think for maximal stability). But where is the Greek_beta_m value of the chapter in the thesis? Or is this Greek_Gamma this Greek_beta_m? Or is it not possible to define this Greek_beta_m? :confused: What are the best values for this factor in your simulations to get max accuracy (in the best case 2nd order) and still a good convergence rate??? Thanks a lot, Rupert |
What you call Greek_Gamma is beta_m / 2.0. There is no general rule what value to pick, you'll have to play around and find what works best for you. You should not specify a value smaller 0.2 though.
|
Thanks for your reply.
Jasak is writing in his thesis that beta should be between 0.1 and 0.5. So do you mean that the Greek_Psi (sorry for the Gamma :rolleyes: mistake) is two times beta? this would fit then with the bounds to Jasaks description ( beta [0.1;0.5] would lead to Psi [0.2;1]). Or is it really 1/2* beta??? Thanks |
Oh yes, sorry about that! I did mean to say that the value you specify in the fvSchemes file will be divided by 2, which gives you the proper beta_m for the limiter (see also Gamma.H in the C++ documentation).
|
thank you very much
|
Quote:
we set the the grad, div, and laplacian schemes in this way: "gradSchemes : default cellMDLimited Gauss linear 0.5; " "divSchemes : div(phi,U) bounded Gauss GammaV 1.0;" "laplacianSchemes : default Gauss linear limited 0.5; " my Question is : when we want play around to found that what value is a good limiter for GammaV , e.g we decrease the limiter of GammaV tp 0.5, we also should decrease the limiters for grad, and laplacian scheme? or they have no relate to each other? Thank you:) |
All times are GMT -4. The time now is 12:31. |