|
[Sponsors] |
kklOmega omega boundary condition at the wall - OF3.0 |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
December 15, 2015, 09:07 |
kklOmega omega boundary condition at the wall - OF3.0
|
#1 |
Senior Member
Artur
Join Date: May 2013
Location: Southampton, UK
Posts: 372
Rep Power: 20 |
Hi All,
I have some concerns about how one should specify the omega BC at the wall when using the kklOmega transition model. In a standard k-omega formulation (low-Re) we would fix the wall and cell centre value of boundary cells using this formulation: http://www.cfd-online.com/Wiki/Near-...k-omega_models which is indeed what happens in omegaWallFunctionFvPatchScalarField.C: Code:
scalar omegaVis = 6.0*nuw[faceI]/(beta1_*sqr(y[faceI])); scalar omegaLog = sqrt(k[cellI])/(Cmu25*kappa_*y[faceI]); omega[cellI] += w*sqrt(sqr(omegaVis) + sqr(omegaLog)); But, the omegaWallFunction BC does not work for the kkl model, as it tries to access the G field which doesn't exist: Code:
--> FOAM FATAL ERROR: request for volScalarField::DimensionedInternalField kkLOmega:G from objectRegistry region0 failed available objects of type volScalarField::DimensionedInternalField are (long list of objects) http://www.engineeringmechanics.cz/pdf/20_5_379.pdf and use zeroGradient the results are very good. Wall values of omega obtained are several orders of magnitudes different than what the wall formulation would suggest, however. Can someone please explain why this model requires a different omega BC than the standard k-omega formulation? Would be great if you can point to the lines in the cause that make the difference. All the best, Artur P.S. my mesh is just an orthogonal, hexahedral, 2D mesh with nearly identical parameters as used in the paper by Furst; tried various 2nd and 1st order schemes, linear solvers, etc., all give consistently the same end result. |
|
Tags |
openfoam 3.0.0, rasmodel, transition boundary layer |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
sliding mesh problem in CFX | Saima | CFX | 46 | September 11, 2021 07:38 |
Natural convection in a closed domain STILL NEEDING help! | Yr0gErG | FLUENT | 4 | December 2, 2019 00:04 |
pressure boundary condition for a rotating wall | barteljoris | CFX | 6 | March 16, 2015 17:21 |
Waterwheel shaped turbine inside a pipe simulation problem | mshahed91 | CFX | 3 | January 10, 2015 11:19 |
Dirichlet boundary condition for additional variable on the wall | ftab | CFX | 13 | January 27, 2013 12:24 |