
[Sponsors] 
Some questions about the functions ddtPhiCorr(rUA, U, phi) and adjustPhi(phi, U, p) 

LinkBack  Thread Tools  Display Modes 
March 25, 2013, 08:29 
Some questions about the functions ddtPhiCorr(rUA, U, phi) and adjustPhi(phi, U, p)

#1 
Member
Join Date: Jun 2009
Posts: 34
Rep Power: 9 
Dear all,
I am reading the icoFoam codes, where I find two functions ddtPhiCorr(rUA, U, phi) and adjustPhi(phi, U, p): volScalarField rUA = 1.0/UEqn.A(); U = rUA*UEqn.H(); phi = (fvc::interpolate(U) & mesh.Sf()) + fvc::ddtPhiCorr(rUA, U, phi); adjustPhi(phi, U, p); Could you tell me why the two functions are used here? And where could I find the explanation of these two functions in official documentation? Thank you very much! lwhy 

March 25, 2013, 11:04 

#2 
Senior Member
Fumiya Nozaki
Join Date: Jun 2010
Location: Yokohama, Japan
Posts: 194
Rep Power: 9 
I'm afraid there aren't any official documents regrading this topic, but you can refer to the source code.
1) adjustPhi(phi, U, p)
2)ddtPhiCorr(rUA, U, phi) You might want to read the following thread fvc::ddtPhiCorr(rUA, U, phi)Hope this helps, Fumiya 

June 5, 2013, 10:15 

#3  
New Member
NaiXian Leslie Lu
Join Date: Jun 2009
Location: France
Posts: 26
Rep Power: 9 
Quote:
Why is adjustPhi taken away in, e.g. compressibleInterFoam? Shouldn't the flux still be adjusted to satisfy continuity at the boundary? Best regards, Lesie 

December 25, 2014, 15:16 

#4 
New Member
Abtin Ansari
Join Date: Oct 2014
Posts: 2
Rep Power: 0 
adjustPhi is taken away from compressible solvers because you are no longer solving Poisson equation for getting the pressure. look at post #4 in pRef and adjustPhi
so there are no singularities to be dealt with. for example lets look at pEqn in laminarSMOKE : fvScalarMatrix pEqn ( fvm::ddt(psi, p) + fvc::div(phi)  fvm::laplacian(rho*rUA, p) ); 

Thread Tools  
Display Modes  

