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

Confused about pEqn

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

Reply
 
LinkBack Thread Tools Display Modes
Old   December 26, 2014, 03:29
Question Confused about pEqn
  #1
New Member
 
BO
Join Date: Dec 2014
Posts: 21
Rep Power: 4
xianqiejiao is on a distinguished road
Hi everyone,

I am currently doing some work on the caveatingFoam. In which the pEqn is as following:

Code:
fvScalarMatrix pEqn
        (
            fvm::ddt(psi, p)
          - (rhol0 + (psil - psiv)*pSat)*fvc::ddt(alphav) - pSat*fvc::ddt(psi)
          + fvc::div(phiv, rho)
          + fvc::div(phiGradp)
          - fvm::laplacian(rhorAUf, p)
        );
According to the literature (Phd thesis by Fabian), the pEqn should be
Code:
fvm::ddt(psi, p)
          - (rhol0 + (psil - psiv)*pSat)*fvc::ddt(alphav) - pSat*fvc::ddt(psi)
          + fvc::div(phiv, rho)
is you write every term one by one.
so what is the last two terms mean? :
Code:
+ fvc::div(phiGradp)
          - fvm::laplacian(rhorAUf, p)
And I am also confused about the last step in pimple loop:
Code:
if (pimple.finalNonOrthogonalIter())
        {
            phiv += (phiGradp + pEqn.flux())/rhof;
        }
Can anyone give me some hint here? If they are for numerical reasons, where can I find information for them? I didn't see much about this in CFD books I have read. I assume they are very important as I can find the similar term in many other solvers.

Kind regards and Merry Christmas!

Last edited by xianqiejiao; December 26, 2014 at 08:40.
xianqiejiao is offline   Reply With Quote

Old   February 29, 2016, 09:36
Default
  #2
New Member
 
Math
Join Date: Jul 2011
Location: chennai
Posts: 14
Rep Power: 7
mg.mithun is on a distinguished road
Send a message via Yahoo to mg.mithun Send a message via Skype™ to mg.mithun
Hi xianqiejiao,
I hope you would have found answers for this question by now. If so can you please share it in the forum so that it will be useful for people like me who are struggling with similar problems. I appreciate your response

MG
mg.mithun 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
Adding a Source Term to Pressure equation pEqn in SimpleFoam oscar_urri OpenFOAM 10 June 4, 2014 12:01
Different Implementation in solving pEqn nishant_hull OpenFOAM Programming & Development 22 June 25, 2013 08:40
Confused about periodic interface Jeff STAR-CCM+ 3 August 15, 2012 22:30
about the pEqn in PISO loop of icoFoam young_Cao OpenFOAM Running, Solving & CFD 4 June 1, 2012 22:36
confused on messages after solver has terminated mactech001 CFX 1 April 21, 2010 16:58


All times are GMT -4. The time now is 04:32.