Quote:
I agree with you about the value of k and epsilon near the wall. If the reynolds is low then the value of k and espilon should be low too. However in the sonicFoam/nacaAirfoil tutorial (compressible case) k and epsilon are have a high fixed value (several thousands). In incompressible/boundaryFoam/boundaryLaunderSharma (incompressible case) k and epsilon have a very low value at the wall (around 1e-10 for both). I just don't understand on what is based this difference, and why this is different for compressible and incompressible fluid. |
Hi malaboss,
Do you find something wich can explain why there is this difference? Valentin |
Well,
I just watch the boundary conditions of both cases (nacaairfoil and boundaryLaunderSharma) which use LaunderSharma. In nacaAirfoil, wall function are used while it should be low Re wall function or nothing at all in a low Re model. This could explain the high values of k at wall, and as a conclusion, Naca Airfoil would not be a good example to understand which boundary conditions we have to set for a Low Re model. See : http://www.cfd-online.com/Wiki/Low-R...els#References for boundary conditions Be careful though, I may be wrong since I don't know how you use a wall function at the wall for a compressible case like naca airfoil. I tried to see how wall fucntion are implemented in Launder and Sharma models, but all I found was this : http://www.cfd-online.com/Forums/ope...k-e-model.html which states that Launder and Sharma was conceived for incompressible cases. About the general implementation of wall function in compressible cases, it seems like there is no big difference with the incompressible. http://jullio.pe.kr/fluent6.1/help/html/ug/node451.htm I did not manage to see the y+ value of the naca airfoil case as I could not run the case. Please tell me if you find something new and more relevant ! |
Thank you for your help malaboss :)
I don't know if launder-sharma is correct in the compressible cases but it is in the list of tubulent models available in compressible. I will give a try to this model in my case and will say if it's ok. Do you know the boundary conditions for nut and alphat that I have to give at walls? I found posts wich say zeroGradient but I'm not sure about this. |
If we keep on reasoning as we did just before, nut should be set to zero (fixed value) on the wall. There is no turbulence so there is no turbulence viscosity, that is, no viscosity due to turbulence.
For alpha t (again i'm not an expert for compressible flows) it represents the thermal diffusivity, proportionnal to thermal conductivity. Near the wall we have nearly no turbulence, so no convection, hence no conductivity. To me, you can fix the value to zero for alpha t. Where did you find the Boundary conditions you are talking about, as it doesn't match my thoughts ? |
I agree with you. According to me zeroGradient doesn't have any sense.
|
Quote:
1) What is the turbulence model in openfoam used? Is it a Low-Re or high-Re model. My suggestion would be Low-Re model as impinging jet involves complex physics that the standard wall function approach might not be able to predict well. From you initial condition settings (epsilonWallFunction, kqRWallFunction, nutSpalartAllmarasWallFunction you are using a high-Re model right?) I tried modelling in fluent impinging jet axisymmetric - use realisable KE and enhanced wall treatment. I attended a course by fluent that taught us that for complex physics it is better to use enhanced wall treatment. In OpenFOAM, there is no equivalent enhanced wall treatment in the Hi-Re Realisable ke model so another way to get around this is to use Low-Re model such as laundersharmaKE which contains damping functions to resolve boundary layer flow. Some people suggest using nut_wall: nutSpalartAllmarasWallFunction or nutSpalartAllmarasStandardWallFunction to make a Hi-Re model looks like low-Re model. I do not think this is right way to do. I have read the thread somewhere on this matter. So strongly recommend to try LaunderSharmaKE. 2) How does your geometry looks like? can you draw out the problem? 3) What are your fvSchemse and fvSolution? Sometimes these are the things that cause continuity to blow up. 4) For nut, is there any reason you use nutSpalartAllmarasStandardWallfunction over nutkWallFunction if using Hi-Re model? Any advantage? |
All times are GMT -4. The time now is 16:39. |