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

YPlusRAS and interFoam

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

Reply
 
LinkBack Thread Tools Display Modes
Old   October 23, 2008, 12:43
Default yPlusRas does not currently wo
  #1
Member
 
Dennis Kingsley
Join Date: Mar 2009
Location: USA
Posts: 45
Rep Power: 8
dkingsley is on a distinguished road
yPlusRas does not currently work with interFoam/rasInterFoma etc. The issue is that the transportProperties dictionary contains data for 2 phases.

Has anyone modified yPlusRas to work for this case, I tried to quickly swap in:

//#include "incompressible/singlePhaseTransportModel/singlePhaseTransportModel.H"
#include "incompressible/incompressibleTwoPhaseMixture/twoPhaseMixture.H"

and calling

//singlePhaseTransportModel laminarTransport(U, phi);
twoPhaseMixture laminarTransport(U, phi);

however, it crashes with the following message:

Selecting incompressible transport model Newtonian
Selecting incompressible transport model Newtonian



request for volScalarField alpha1 from objectRegistry region0 failed
available objects of type volScalarField are

3
(
nu2
nu1
yPlus
)

Can anyone point me in the right direction, I would like the tool to ultimately give me the yPlus values for each phase just like it does for a single phase based on the initial definition of gamma.
dkingsley is offline   Reply With Quote

Old   October 23, 2008, 14:23
Default Dennis, The source of your
  #2
New Member
 
James Criner
Join Date: Mar 2009
Posts: 7
Rep Power: 8
jcriner is on a distinguished road
Dennis,

The source of your error is likely the lack of the gamma field (i.e. you need to read in the gamma field and create the twoPhaseMixture properties for mu() and nu()). Have a look at

src/turbulenceModels/RAS/incompressible/RASModel/RASModel.C
src/transportModels/incompressible/incompressibleTwoPhaseMixture/twoPhaseMixture .C

to see what is going on. Then you could follow the source in createFields.H of rasInterFoam for how to create (read in) the gamma field and create the twoPhaseMixture properties.

James
jcriner is offline   Reply With Quote

Old   October 23, 2008, 20:32
Default James, Thanks for the advic
  #3
Member
 
Dennis Kingsley
Join Date: Mar 2009
Location: USA
Posts: 45
Rep Power: 8
dkingsley is on a distinguished road
James,

Thanks for the advice, I have it working now.

Dennis
dkingsley is offline   Reply With Quote

Old   October 23, 2008, 22:25
Default Next question, I would like to
  #4
Member
 
Dennis Kingsley
Join Date: Mar 2009
Location: USA
Posts: 45
Rep Power: 8
dkingsley is on a distinguished road
Next question, I would like to override the RASModel::YPlus method in the SpalartAllmaras class to add the alternative yPlus formulation that I found in this older post http://www.cfd-online.com/cgi-bin/Op...c=1&post=18488.

Is this a reasonable approach? Is there a more suitable formulation?

tmp<scalarfield> tYp =
sqrt(mag(turbulence->nuEff()().boundaryField()[patchI]*
U.boundaryField()[patchI].snGrad()))/
nu.value()*nearWallDist(mesh).y()[patchI];

scalarField& Yp = tYp();
dkingsley is offline   Reply With Quote

Old   April 28, 2010, 09:08
Default
  #5
Senior Member
 
Claus Meister
Join Date: Aug 2009
Location: Wiesbaden, Germany
Posts: 241
Rep Power: 9
idrama is on a distinguished road
Hello dkingsley,

I read your post! I just want to know what have you change in yPlusRAS.C, i.e. have you only made that change:

//singlePhaseTransportModel laminarTransport(U, phi);
twoPhaseMixture laminarTransport(U, phi);

or have you done a couple of things more?

Cheers
idrama 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
Transient interFoam with MRF jaswi OpenFOAM Running, Solving & CFD 18 October 20, 2011 23:29
Writing yPlusRAS Values velan OpenFOAM Running, Solving & CFD 0 December 30, 2008 07:09
Outlet bc in interFoam happysimulant OpenFOAM Running, Solving & CFD 0 November 6, 2008 11:06
InterFoam floooo OpenFOAM Running, Solving & CFD 0 November 3, 2008 12:00
About interfoam solver qiu OpenFOAM Running, Solving & CFD 0 May 6, 2007 22:48


All times are GMT -4. The time now is 07:06.