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

add a term in continuity equations for twoPhaseEulerfoam

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

LinkBack Thread Tools Search this Thread Display Modes
Old   March 24, 2020, 13:59
Question add a term in continuity equations for twoPhaseEulerfoam
Join Date: Mar 2019
Posts: 57
Rep Power: 3 is on a distinguished road
Hi guys, I need to add a term into the mass conservation equation for twophaseeulerfoam, as shown in following codes which are in the p-equation. I hope I added it in a right position. I think I should use plus sign for this term however the simulation crashed when I use "plus". While I use "minus". It goes well without problem. Do you have any suggestion?
pEqnComp1 =
- fvc::Sp(fvc::ddt(alpha1) + fvc::div(alphaPhi1), rho1)
+ (alpha1*psi1/rho1)*correction(fvm::ddt(p_rgh)) + fvc::laplacian(phase2.turbulence().nut()/mysigma,alpha1);

pEqnComp2 =
- fvc::Sp(fvc::ddt(alpha2) + fvc::div(alphaPhi2), rho2)
+ (alpha2*psi2/rho2)*correction(fvm::ddt(p_rgh)) + fvc::laplacian(phase2.turbulence().nut()/mysigma,alpha2); is offline   Reply With Quote


twophase fluid

Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On

Similar Threads
Thread Thread Starter Forum Replies Last Post
p_rgh initial residual no change with different settings manuc OpenFOAM Running, Solving & CFD 3 June 26, 2018 15:53
time step continuity problem in VAWT simulation lpz_michele OpenFOAM Running, Solving & CFD 5 February 22, 2018 19:50
Help for the small implementation in turbulence model shipman OpenFOAM Programming & Development 25 March 19, 2014 10:08
AMI interDyMFoam for mixer nu problem danny123 OpenFOAM Programming & Development 8 September 6, 2013 02:34
How to write k and epsilon before the abnormal end xiuying OpenFOAM Running, Solving & CFD 8 August 27, 2013 15:33

All times are GMT -4. The time now is 08:47.