# icoFoam's pressure possion equation.

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

July 15, 2013, 03:49
icoFoam's pressure possion equation.
#1
Senior Member

Dongyue Li
Join Date: Jun 2012
Location: Torino, Italy
Posts: 754
Rep Power: 10
Hi

In icoFoam we have:
Code:
fvScalarMatrix pEqn
(
fvm::laplacian(rAU, p) == fvc::div(phiHbyA)
);
But via my equation(see the image), Shouldnt it be

Code:
fvScalarMatrix pEqn
(
fvm::laplacian(rAU, p) == fvc::div(HbyA)
);
Thanks
Attached Images
 1.jpg (48.1 KB, 56 views)

 July 16, 2013, 03:30 #2 Senior Member   Dongyue Li Join Date: Jun 2012 Location: Torino, Italy Posts: 754 Rep Power: 10 div(phiHbyA)? Im confused about this divergence of a scalar thing.

 July 16, 2013, 21:34 #3 Senior Member   Dongyue Li Join Date: Jun 2012 Location: Torino, Italy Posts: 754 Rep Power: 10 I changed it to Code: fvScalarMatrix pEqn ( fvm::laplacian(rAU, p) == fvc::div(HbyA) ); It compiles successfully and works in cavity and elbow case!And the solution is almost the same with Code: fvScalarMatrix pEqn ( fvm::laplacian(rAU, p) == fvc::div(PhiHbyA) ); Whats the problem.. Code: fvc::div(HbyA):dimensions [0 0 -1 0 0 0 0]; internalField nonuniform List 9(0.157469 -0.0674757 -0.487085 -0.960691 0.789192 0.715649 30.1762 0.258981 -30.5822); fvc::div(phiHbyA):dimensions [0 0 -1 0 0 0 0]; internalField nonuniform List 9(-0.435688 -0.0982547 0.152804 -0.0448476 0.720155 -0.152663 28.2072 0.236062 -28.5848);

 July 26, 2013, 12:39 #4 Senior Member     Santiago Marquez Damian Join Date: Aug 2009 Location: Santa Fe, Santa Fe, Argentina Posts: 430 Rep Power: 17 Hi, check the p-U coupling section (2.7) in my Ph.D thesis. It will also give you some pointers to good material on this topic. Link to the PhD thesis is provided in my user page at openfoamwiki.net: http://openfoamwiki.net/index.php/User:Santiagomarquezd Regards. __________________ Santiago MÁRQUEZ DAMIÁN, Ph.D. Research Scientist Research Center for Computational Methods (CIMEC) - CONICET/UNL Tel: 54-342-4511594 Int. 7032 Colectora Ruta Nac. 168 / Paraje El Pozo (3000) Santa Fe - Argentina. http://www.cimec.org.ar Last edited by wyldckat; May 24, 2016 at 17:35. Reason: repaired link and description to link

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post Mohsin FLUENT 36 April 29, 2016 17:16 Don456 Main CFD Forum 1 January 19, 2012 16:00 renyun0511 OpenFOAM Programming & Development 0 November 4, 2010 02:38 Sara Main CFD Forum 22 July 30, 2009 04:46 Rui Main CFD Forum 6 October 15, 2000 04:57