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

icoFoam's pressure possion equation.

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

Reply
 
LinkBack Thread Tools Display Modes
Old   July 15, 2013, 03:49
Default icoFoam's pressure possion equation.
  #1
Senior Member
 
Dongyue Li
Join Date: Jun 2012
Location: Torino, Italy
Posts: 676
Rep Power: 8
sharonyue is on a distinguished road
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
File Type: jpg 1.jpg (48.1 KB, 41 views)
sharonyue is offline   Reply With Quote

Old   July 16, 2013, 03:30
Default
  #2
Senior Member
 
Dongyue Li
Join Date: Jun 2012
Location: Torino, Italy
Posts: 676
Rep Power: 8
sharonyue is on a distinguished road
div(phiHbyA)?
\nabla \cdot a

Im confused about this divergence of a scalar thing.
sharonyue is offline   Reply With Quote

Old   July 16, 2013, 21:34
Default
  #3
Senior Member
 
Dongyue Li
Join Date: Jun 2012
Location: Torino, Italy
Posts: 676
Rep Power: 8
sharonyue is on a distinguished road
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<scalar> 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<scalar> 9(-0.435688 -0.0982547 0.152804 -0.0448476 0.720155 -0.152663 28.2072 0.236062 -28.5848);
sharonyue is offline   Reply With Quote

Old   July 26, 2013, 12:39
Default
  #4
Senior Member
 
santiagomarquezd's Avatar
 
Santiago Marquez Damian
Join Date: Aug 2009
Location: Santa Fe, Santa Fe, Argentina
Posts: 418
Rep Power: 15
santiagomarquezd will become famous soon enough
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
santiagomarquezd 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
lid-driven cavity in matlab using BiCGStab Don456 Main CFD Forum 1 January 19, 2012 16:00
Pressure Outlet Guage pressure Mohsin FLUENT 34 January 12, 2012 06:45
How to add correcting pressure equation in SIMPLE of MRFSimpleFOAM? renyun0511 OpenFOAM Programming & Development 0 November 4, 2010 02:38
Pressure enigma. Sara Main CFD Forum 22 July 30, 2009 04:46
pressure Poisson equation in fractinoal step metho Rui Main CFD Forum 6 October 15, 2000 04:57


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