CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM Programming & Development


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

LinkBack Thread Tools Display Modes
Old   May 14, 2015, 13:32
Default mixture.correct()?
Join Date: Apr 2015
Location: INDIA
Posts: 50
Rep Power: 2
Saideep is on a distinguished road
hi FOAMers;

Wishing you a very good day.

I am working with the new version(OF-2.3.1) interFOAM and trying to change from default PIMPLE to PISO solver. On the way, I have a question in the original openfoam interFoam.C file, regarding, from where is the mixture.correct() function being called to correct the kinematic viscosity. In older versions it was twoPaseProperties.correct(). Even in older versions of github files i failed to find the object twoPhaseProperties.

Could anyone shed some light over this topic?

Thanks and waiting for your answer.

Saideep is offline   Reply With Quote

Old   May 16, 2015, 06:08
Senior Member
Join Date: Aug 2014
Posts: 147
Rep Power: 2
ssss is on a distinguished road
I think that it is dealt inside the turbulenceModel, so when you create the turbulence model:

    autoPtr<incompressible::turbulenceModel> turbulence
        incompressible::turbulenceModel::New(U, phi, mixture)
It automatically reads from the mixture the mu's, and corrects the turbulence itself, without calling any other function than:

turbulence->divDevRhoReff(rho, U)
This should be the way it should work I think
ssss is offline   Reply With Quote

Old   May 16, 2015, 07:30
Default Update!!
Join Date: Apr 2015
Location: INDIA
Posts: 50
Rep Power: 2
Saideep is on a distinguished road
Thanks a lot buddy!!

Appreciate your time and effort in my case. I guess even you are on the same track as of mine.

As of now I understand a bit of your answer(not completely as I am not working with turbulence as of now, i am just changing the default damBreak code from PIMPLE to PISO.) {Reason: I am dealing with interface motion of 2 phases. So, papers that i had read from Ivan Lunati, Raeini, Duon and so on all use PISO algorithm instead of PIMPLE. I dont know the reason as of now and for basic simulations I compared between my results and their results and they are different, so I thought to have a go on PISO once to see if it really works.}

I would be greatful if anyone could explain how PISO is different from PIMPLE (not theoritical(PIMPLE = PISO+SIMPLE) but with help of code.).

Now coming to our original question between twoPhaseProperties and mixture, I someow feel it is deep hidden somewhere and the latest version from 2.3.x the twoPhaseProperties have been replaced with mixture object.{Still searching where it is defined and I find this usage in net but its definition is not been seen.}
And search for definition continues and if anybody comes across this it is worth to share.

Thanks a lot once again from my side for your intrest!!
hope to find the solution soon.
Saideep is offline   Reply With Quote


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

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