CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Programming & Development

implement roughness boundary conditions for Spalart-Allmaras

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   June 25, 2009, 17:34
Default implement roughness boundary conditions for Spalart-Allmaras
  #1
New Member
 
wei wu
Join Date: Mar 2009
Location: us
Posts: 14
Rep Power: 8
wei_wu is on a distinguished road
Hello, all

I'm a new foamer and need help on this problem.

I'm trying to apply the roughness boundary condition on the wall with S-A model.

the wall boundary condition will be : d(nutilda)/dn=nutilda/(d_+0.03*hs)
where d_ is the distance to the wall, which is a variable in the S_A model
hs is the roughness height

my question is howI can implement this new boundary condition.
here is my thought:
there is one basic fixedGradient boundary condition which maybe is a good prototype for me. I will change the evaluate member function in order to get value of (nutilda/(d_+0.03*hs)) for the gradient of nutilda , but I'm so clumsy at C++ and FOAM right now and don't know how to make it happen.

here is the evaluate function:
__________________________________________________ _______________________
template<class Type>
void fixedGradientFvPatchField<Type>::evaluate(const Pstream::commsTypes)
{
if (!this->updated())
{
this->updateCoeffs();
}

Field<Type>:perator=
(
this->patchInternalField() + gradient_/this->patch().deltaCoeffs()
);

fvPatchField<Type>::evaluate();
}
__________________________________________________ ________________________

Thank you all
Wei
wei_wu is offline   Reply With Quote

Reply

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
inlet velocity boundary condition murali CFX 5 August 3, 2012 08:56
Concentric tube heat exchanger (Air-Water) Young CFX 5 October 6, 2008 23:17
Fluent accuracy and boundary conditions Paolo Lampitella FLUENT 0 June 12, 2008 06:25
Outflow boundary conditions Achilleas Tsompanos Main CFD Forum 6 April 2, 2000 16:09
Boundary Conditions Jan Ramboer Main CFD Forum 11 August 16, 1999 08:59


All times are GMT -4. The time now is 17:27.