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

programming & icoFoam

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

Like Tree1Likes
  • 1 Post By sharonyue

Reply
 
LinkBack Thread Tools Display Modes
Old   June 9, 2013, 06:02
Default programming & icoFoam
  #1
Senior Member
 
Sasan Ghomi
Join Date: Sep 2012
Location: Iran - Tehran
Posts: 185
Rep Power: 4
sasanghomi is on a distinguished road
Hi foamers ,
what is the relation between non-orthogonality in the grid and number of solving P.Equ ???
please take a look at this part of icoFoam :

PHP Code:
for (int nonOrth=0nonOrth<=nNonOrthCorrnonOrth++)
            {
                
fvScalarMatrix pEqn
                
(
                    
fvm::laplacian(rUAp) == fvc::div(phi)
                );

                
pEqn.setReference(pRefCellpRefValue);
                
pEqn.solve();

                if (
nonOrth == nNonOrthCorr)
                {
                    
phi -= pEqn.flux();
                }
            } 
why the number of solving P.Equ should be equal nNonOrthogonalCorrectors ??
I think in PISO algorithm the number of solving pressure equation should be equal nCorrectors. I think nNonOrthogonalCorrectors is related to quality of grid generation.
what is the relation between grid's quality and number of solving pressure equation??
Can anyone clarify this issue for me ?

Thanks and best regards,
Sasan.
sasanghomi is offline   Reply With Quote

Old   June 9, 2013, 15:27
Default
  #2
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 8,301
Blog Entries: 34
Rep Power: 84
wyldckat is just really nicewyldckat is just really nicewyldckat is just really nicewyldckat is just really nice
Hi Sasan,

Does this answer your questions: http://openfoamwiki.net/index.php/IcoFoam ?

Best regards,
Bruno
wyldckat is online now   Reply With Quote

Old   June 9, 2013, 16:51
Default
  #3
Senior Member
 
Sasan Ghomi
Join Date: Sep 2012
Location: Iran - Tehran
Posts: 185
Rep Power: 4
sasanghomi is on a distinguished road
Hi Bruno,
Thanks for your reply.
I want to know that why the quantity of solving pressure equation is equal to quantity of nNonOrthogonalCorrectors ?

Thanks and best regards,
Sasan.
sasanghomi is offline   Reply With Quote

Old   June 9, 2013, 17:33
Default
  #4
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 8,301
Blog Entries: 34
Rep Power: 84
wyldckat is just really nicewyldckat is just really nicewyldckat is just really nicewyldckat is just really nice
Hi Sasan,

I'm not sure I understand you correctly... are you asking about the following code?
Code:
if (nonOrth == nNonOrthCorr)
{
    phi -= pEqn.flux();
}
edit: or is it about:
Code:
for (int nonOrth=0; nonOrth<=nNonOrthCorr; nonOrth++)
Best regards,
Bruno

Last edited by wyldckat; June 9, 2013 at 17:34. Reason: see "edit:"
wyldckat is online now   Reply With Quote

Old   June 10, 2013, 02:31
Default
  #5
Senior Member
 
Sasan Ghomi
Join Date: Sep 2012
Location: Iran - Tehran
Posts: 185
Rep Power: 4
sasanghomi is on a distinguished road
Hi Bruno,

I am asking about :

Code:
for (int nonOrth=0; nonOrth<=nNonOrthCorr; nonOrth++)
Thanks and best regards,
Sasan.
sasanghomi is offline   Reply With Quote

Old   June 10, 2013, 07:24
Default
  #6
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 8,301
Blog Entries: 34
Rep Power: 84
wyldckat is just really nicewyldckat is just really nicewyldckat is just really nicewyldckat is just really nice
Hi Sasan,

Well, as the wiki page explains:
Quote:
Originally Posted by http://openfoamwiki.net/index.php/IcoFoam
// iteratively correct for non-orthogonality. The non-orthogonal part of the Laplacian is calculated from the most recent
// solution for pressure, using a deferred-correction approach.
Any explanation beyond this (that I could find) is on this thesis:
Quote:
Originally Posted by http://openfoamwiki.net/index.php/IcoFoam
... which I don't have time to read and study .

Best regards,
Bruno
wyldckat is online now   Reply With Quote

Old   August 25, 2013, 23:55
Default
  #7
Senior Member
 
Dongyue Li
Join Date: Jun 2012
Location: Torino, Italy
Posts: 675
Rep Power: 8
sharonyue is on a distinguished road
Hi ,

Description of flux() method
I think this will make it clear.
It is related to quality of grid generation and the non-orthogonal correction in the diffusion term, you can dig it in page.83 of Jasak's thesis.
wyldckat likes this.
sharonyue 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
Density in icoFoam Densidad en icoFoam manuel OpenFOAM Running, Solving & CFD 8 September 22, 2010 04:10
Programming in OpenFOAM vinu OpenFOAM 2 July 11, 2009 10:16
IcoFoam Not Installed and blockMesh command not found cesarbz OpenFOAM Installation 6 July 4, 2008 04:44
Kubuntu uses dash breaks All scripts in tutorials platopus OpenFOAM Bugs 8 April 15, 2008 07:52
IcoFoam parallel woes msrinath80 OpenFOAM Running, Solving & CFD 9 July 22, 2007 02:58


All times are GMT -4. The time now is 08:46.