# 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: 676
Rep Power: 8
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, 41 views)

 July 16, 2013, 03:30 #2 Senior Member   Dongyue Li Join Date: Jun 2012 Location: Torino, Italy Posts: 676 Rep Power: 8 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: 676 Rep Power: 8 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: 418 Rep Power: 15 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. https://docs.google.com/file/d/0B2lp...it?usp=sharing Regards. __________________ Santiago MÁRQUEZ DAMIÁN, Ph.D. Post-doctoral Fellow Research Center for Computational Mechanics (CIMEC) - CONICET/FICH-UNL T.E.: 54-342-4511594 Ext. 1005 Güemes 3450 - (3000) Santa Fe Santa Fe - Argentina http://www.cimec.org.ar

 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 Don456 Main CFD Forum 1 January 19, 2012 16:00 Mohsin FLUENT 34 January 12, 2012 06:45 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

All times are GMT -4. The time now is 23:21.