CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM (http://www.cfd-online.com/Forums/openfoam/)
-   -   how to calculate Cf inside a wall function (http://www.cfd-online.com/Forums/openfoam/85042-how-calculate-cf-inside-wall-function.html)

 aerothermal February 15, 2011 20:02

how to calculate Cf inside a wall function

Dear foamers,

Could you tell me how to calculate Cf (friction coefficient) inside a wall function like alphatWallFunctionFvPatchScalarField.C ?

I already have yPlus, rho and uStar=Cmu^(0.25)*k^(0.5).

Code:

```    forAll(mutw, faceI)     {         label faceCellI = patch().faceCells()[faceI];         scalar pr = muw[faceI]/alphaw[faceI];         scalar uStar = Cmu25*sqrt(k[faceCellI]);         scalar KsPlus = uStar*Ks/(muw[faceI]/rhow[faceI]); ....     }```
Since uStar is approx. uTau=sqrt(tau/rho) for equilibrium boundary layers, now I need a boundary-layer edge velocity Ue.

In order to be more accurate, can I calculate uTau with lines below?

Code:

```    const fvPatchVectorField& U =         patch().lookupPatchField<volVectorField, vector>("U");                scalarField magFaceGradU = mag(U.snGrad());```
and inside loop of mutw in all faceI
Code:

```        scalar utau =             sqrt((mutw[faceI] + muw[faceI])*magFaceGradU[faceI]/rhow[faceI]);```
I see two options but they are not easy for me yet since it is my first development with OF:

1) I know that I can recalculate the uPlus with yPlus to get sqrt(Cf/2), since uPlus=1/sqrt(Cf/2)

2) If I use the other definition of uPlus, the problem will be to know Ue in order to find Cf=tau/(0.5*Ue*Ue*rho).
uPlus= Up*Cmu^(0.25)*k^(0.5)/(tau/rho)
Even I adopt such calculation, how to get value of Up? Is it simply U[faceI]?

I would like to know whether or not there is an easier way. Are the methods above acceptable? Are they correct in programming?

Regards,

aerothermal

 aerothermal July 15, 2011 11:59

No clues on that?

 All times are GMT -4. The time now is 14:44.