[shereez234]

Hello experts in OpenFoam.

I need some assistance in coding a zero equation turbulence model i.e. Baldwin Lomax.

http://www.cfd-online.com/Wiki/Baldwin-Lomax_model

I have tried to modify the two equation turbulence models and one equation models to remove the PDEs and convert this into a fully algebraic zero equation model without much success.

Another Idea I have is to start from Laminar so that I don't have to remove the PDE's. What is the best way to go about it?

Also if any one else has done some work or know more information please help me out.

With kind regards

Shereez

[shereez234]

Hello;

I have made some progress in removing 2 PDE's from a 2-equation model. But Baldwin Lomax uses y+ (yPlus) for it's execution. Wall distance y_ is easily incorporated in to the code but does any one know how I can incorporate yplus in to my source code without having to redefine the whole calculation of yplus again...

Thanks

[shereez234]

Going forward;

now I am stuck here:

I have a function

F= a * b * y_ * C

I would like to return the Ymax which is the corresponding y_ while F = Fmax (i,e, maximum of function F)

I tried

while ( F == Fmax)
{

return y_(Fmax);

}


But it doesnt work. Can somebody please help me


Thanks

[RodriguezFatz]

Hey shereez,
Your code looks horrible
I can probably help you, but I need more information.
What are a, b, C and y_ ? Functions or constants?
Philipp.

[shereez234]

Originally Posted by RodriguezFatz

Hey shereez,
Your code looks horrible
I can probably help you, but I need more information.
What are a, b, C and y_ ? Functions or constants?
Philipp.

tmp<volScalarField> Baldwin_Fake::L() const
 {
     return f_ * y_ *(1 - exp(yPlus/Aplus_));
 }

 tmp<volScalarField> Baldwin_Fake:: Fnum(const volScalarField& G1) const
 {
   return  (y_ * G1 * (1 - exp(yPlus/Aplus_)));
 }



 tmp<volScalarField> Baldwin_Fake::Fmax(const volScalarField &Fnum) const
  {
      tmp<volScalarField> Y1 = max
      (Fnum, 1e-40);
     return Y1;
}

tmp<volScalarField> Baldwin_Fake::Ymax(const volScalarField& Fmax, const volScalarField& G1) const
   {

                 return  (Fmax / (G1 *(1 - exp(yPlus/Aplus_)) ));

    }
//
  tmp<volScalarField> Baldwin_Fake::Fkleb(const volScalarField& Ymax) const
  {
      tmp<volScalarField> H;
       (1 +5.5 * pow(((y_ * Ckleb_)/Ymax),(6)));

       return pow(H,-1);
  }



  tmp<volScalarField> Baldwin_Fake::YmaxFmax
  (
    const volScalarField& Ymax,
      const volScalarField& Fmax
 ) const
  {
      return (Ymax * Fmax);
  }
//
//
 tmp<volScalarField> Baldwin_Fake::Fwake(const volScalarField& YmaxFmax, const volScalarField& Fmax, const volScalarField& Udiff, const volScalarField& Ymax) const
  {

  //   ((YmaxFmax, Cwk_ * Ymax * ((Udiff * Udiff))/Fmax)),1e-40);
        return max ((((sqr(Udiff))/Fmax) * Cwk_ * Ymax), YmaxFmax );

  }

 tmp<volScalarField> Baldwin_Fake::nut_Inner(const volScalarField& L, const volScalarField& G1) const
  {
     return ((sqr(L)) * G1 * scalar(1.2));
  }

 tmp<volScalarField> Baldwin_Fake::nut_Outer(const volScalarField& Fkleb, const volScalarField& Fwake) const
  {
     return (Fwake* Fkleb * Ccp_ * K_ * scalar(1.2));
  }

 tmp<volScalarField> Baldwin_Fake::yCrossOver(const volScalarField& nut_Inner, const volScalarField& nut_Outer) const
  {
    if (nut_Inner = nut_Outer)
    {
        return 0.0;
    }
  }

Yeah that wasnt the real code. Some how I am able to get the Ymax now. At the moment I am stuck in a position with yCrossOver (the last function in the code) nut_Inner = nut_Outer as you can see from above. I need to return thee min(y_) while nut_Inner is equal to nut_Outer. Do you think I need a bool function to do this?

Thanks for the help I am not an expert in C++ coding

[shereez234]

  tmp<volScalarField> Baldwin_Edit2::Ymax(const volScalarField& Fmax) const
   {
        return y_(Fmax);
   }

I have tried several different ways. I would like to know how I can return the value of wall distance y_ when Fmax = maximum.

Is there any one who can help me?

[RodriguezFatz]

Hi again,

Sorry I forgot about this thread...
Well, I don't understand your code or I don't understand the model!
Can you post the whole code? It looks like you use y+ as a constant...