|
[Sponsors] |
October 23, 2008, 12:43 |
yPlusRas does not currently wo
|
#1 |
Member
Dennis Kingsley
Join Date: Mar 2009
Location: USA
Posts: 45
Rep Power: 17 |
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. |
|
October 23, 2008, 14:23 |
Dennis,
The source of your
|
#2 |
New Member
James Criner
Join Date: Mar 2009
Posts: 7
Rep Power: 17 |
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 |
|
October 23, 2008, 20:32 |
James,
Thanks for the advic
|
#3 |
Member
Dennis Kingsley
Join Date: Mar 2009
Location: USA
Posts: 45
Rep Power: 17 |
James,
Thanks for the advice, I have it working now. Dennis |
|
October 23, 2008, 22:25 |
Next question, I would like to
|
#4 |
Member
Dennis Kingsley
Join Date: Mar 2009
Location: USA
Posts: 45
Rep Power: 17 |
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(); |
|
April 28, 2010, 09:08 |
|
#5 |
Senior Member
Claus Meister
Join Date: Aug 2009
Location: Wiesbaden, Germany
Posts: 241
Rep Power: 18 |
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 |
|
Thread Tools | Search this Thread |
Display Modes | |
|
|
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 06:09 |
Outlet bc in interFoam | happysimulant | OpenFOAM Running, Solving & CFD | 0 | November 6, 2008 10:06 |
InterFoam | floooo | OpenFOAM Running, Solving & CFD | 0 | November 3, 2008 11:00 |
About interfoam solver | qiu | OpenFOAM Running, Solving & CFD | 0 | May 6, 2007 22:48 |