
[Sponsors] 
January 23, 2011, 09:42 
fvc::ddtPhiCorr(rUA, U, phi)

#1 
Senior Member
ata kamyabi
Join Date: Aug 2009
Location: Kerman
Posts: 322
Rep Power: 9 
Hi FOAMers
can any one tell me fvc::ddtPhiCorr(rUA, U, phi) in icoFoam code is for what? best regards 

January 24, 2011, 13:05 

#2 
Senior Member
Santiago Marquez Damian
Join Date: Aug 2009
Location: Santa Fe, Santa Fe, Argentina
Posts: 430
Rep Power: 15 
__________________
Santiago MÁRQUEZ DAMIÁN, Ph.D. Research Scientist Research Center for Computational Methods (CIMEC)  CONICET/UNL Tel: 543424511594 Int. 7032 Colectora Ruta Nac. 168 / Paraje El Pozo (3000) Santa Fe  Argentina. http://www.cimec.org.ar 

January 24, 2011, 23:02 
fvc::ddtPhiCorr(rUA, U, phi)

#3 
Senior Member
ata kamyabi
Join Date: Aug 2009
Location: Kerman
Posts: 322
Rep Power: 9 
Thank you very much


August 19, 2011, 22:29 
http://www.openfoam.com/mantisbt/print_bug_page.php?bug_id=169

#4 
New Member
M. Timaji
Join Date: Sep 2010
Location: Tehran
Posts: 6
Rep Power: 7 
the term Phi in UEqn and PEqn is not exactly the same:
In the UEqn: corrected flux from the previous iteration or timestep, In the PEqn: mass flux without the pressure contribution. The term ddtPhiCorr checks for the dimensional units of U and phi to decide what operation has to be performed, and if phi is defined in terms of mass, a division of phi by rho is performed. This represents a problem in the case of zero density. Such a case cannot happen in singlephase flows, however it might happen in multiphase flows. For example, let's consider a momentum equation in the form ddt(alpha*rho*U) + div(alpha*rho*U*U) = ... where alpha is the phase fraction. The equation is then represented in the code as fvm::ddt(alphaRho, U) + fvm::div(alphaf*phi, U) == ... being phi = rho_f * (U_f \cdot S). In such a case, if my understanding is correct, one should compute surfaceScalarField phi = fvc::interpolate(rho)*(fvc::interpolate(U) & mesh.Sf()) + fvc::ddtPhiCorr(rUA, rhoAlpha, U, phi) At this point, since alpha can be zero, ddtPhiCorr will cause a division by zero. This does not represent a problem in the incompressible case, since the equation can be divided by rho, which leads phi to be a volumetric flux, and ddtPhiCorr would be ddtPhiCorr(rUA, alpha, U, phi). However, if one wants to deal with the compressible case, keeping equations in conservative form, the problem appears. 

August 21, 2011, 03:54 

#5 
Senior Member
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,904
Rep Power: 27 
The problem is known and was reported (I guess you read it since you pasted it :P):
http://www.openfoam.com/mantisbt/view.php?id=169
__________________
Alberto Passalacqua GeekoCFD  A free distribution based on openSUSE 64 bit with CFD tools, including OpenFOAM. Available as in both physical and virtual formats. OpenQBMM  An opensource implementation of quadraturebased moment methods. To obtain more accurate answers, please specify the version of OpenFOAM you are using. 

September 10, 2011, 12:47 

#6  
New Member
M. Timaji
Join Date: Sep 2010
Location: Tehran
Posts: 6
Rep Power: 7 
Quote:
Yes alberto I referenced to the link you've mentiond in the title of my post!!! 

February 7, 2013, 23:10 
Simple

#7 
Member
,...
Join Date: Apr 2011
Posts: 92
Rep Power: 6 
Hi FOAMERS
I have modified SIMPLEFOAM for unsteady flows bu adding ddt(U) to the U matrix. Do I also have to add fvc::ddtPhiCorr in PEqu.H to calculate fluxes? How ignoring this term will affect my results? 

Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
phi = pEqn.flux() vs. linearInterpolate(U) & mesh.Sf()  santiagomarquezd  OpenFOAM Programming & Development  32  June 12, 2014 01:50 
Turbulence Model phi vs phi_  doug  OpenFOAM Running, Solving & CFD  4  November 10, 2009 05:33 
Another phi question  ehsan_vaghefi  OpenFOAM Running, Solving & CFD  0  October 24, 2008 19:56 
What does the fvcddtPhiCorrrUA U phi and fvcddtPhiCorrrUA rho U phi mean Any references  dbxmcf  OpenFOAM Running, Solving & CFD  0  October 1, 2008 21:43 
About phi in icoFoam  kar  OpenFOAM Running, Solving & CFD  3  February 20, 2008 06:20 