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

Maybe a bug in twoLiquidMixingFoam?

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

LinkBack Thread Tools Display Modes
Old   March 6, 2016, 04:40
Default Maybe a bug in twoLiquidMixingFoam?
Senior Member
Dongyue Li
Join Date: Jun 2012
Location: Torino, Italy
Posts: 756
Rep Power: 10
sharonyue is on a distinguished road
Hi guys,

In twoLiquidMixingFoam, at first we solve alpha by MULES, then rhoPhi is calculated by the following:

rhoPhi = alphaPhi*(rho1 - rho2) + phi*rho2;
This makes sense. Then alpha is solved in alphaDissEqn.H, rhoPhi is updated as following:

rhoPhi += alpha1Eqn.flux()*(rho1 - rho2);

If we put them together, finally, rhoPhi will be:

rhoPhi = (alphaPhi + alpha1Eqn.flux())*(rho1 - rho2) + phi*rho2;

I think we need to delete the old alphaPhi in the code above, so in the alphaDissEqn.H, it should be:

rhoPhi = alpha1Eqn.flux()*(rho1 - rho2) + phi*rho2;
Am I right?
Im the translator of OpenFOAM User Guide Chinese Edition! But always newbie on CFD.
sharonyue 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

Similar Threads
Thread Thread Starter Forum Replies Last Post
Bug in Workbench CFX Pierre1 CFX 6 August 2, 2017 00:18
SU2_MSH: Periodic boundary conditions bug Zef SU2 1 February 18, 2015 14:28
Serious bug in LES interface fs82 OpenFOAM Bugs 21 November 16, 2009 09:15
Bug reports Mattijs Janssens (Mattijs) OpenFOAM 0 January 10, 2005 11:05
Forum y2k Bug Jonas Larsson Main CFD Forum 1 January 5, 2000 11:22

All times are GMT -4. The time now is 03:48.