CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

Decomposition Methods: which ones are best for large scale simulations?

Register Blogs Community New Posts Updated Threads Search

Like Tree4Likes
  • 1 Post By olesen
  • 1 Post By rnburne
  • 2 Post By Koushikc

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   January 7, 2010, 13:03
Default Decomposition Methods: which ones are best for large scale simulations?
  #1
New Member
 
Join Date: Nov 2009
Posts: 17
Rep Power: 17
misakagan is on a distinguished road
What decomposition method (Metis, Scotch, simple, hierarchical) you guys favor for large scaled simulations?
misakagan is offline   Reply With Quote

Old   January 8, 2010, 03:22
Default
  #2
Senior Member
 
Mark Olesen
Join Date: Mar 2009
Location: https://olesenm.github.io/
Posts: 1,714
Rep Power: 40
olesen has a spectacular aura aboutolesen has a spectacular aura about
Quote:
Originally Posted by misakagan View Post
What decomposition method (Metis, Scotch, simple, hierarchical) you guys favor for large scaled simulations?
Metis or scotch is what we use, since we generally don't have a better idea ourselves of how to decompose the domain any better (nor the time to figure it out anew for each geometry / cpu count combination).

Simply try out metis/scotch and see if you are happy with the results. You can view the separate processor* domains in paraview (eg, via paraFoam) to see if it matches up with your expectations.

BTW: I believe that future versions of OpenFOAM might be switching to scotch instead of metis due to licensing aspects. I found that the scotch decomposition was fairly similar to metis anyhow.
rajibroy likes this.
olesen is offline   Reply With Quote

Old   January 8, 2010, 05:35
Default
  #3
Senior Member
 
BastiL
Join Date: Mar 2009
Posts: 530
Rep Power: 20
bastil is on a distinguished road
Quote:
Originally Posted by olesen View Post
BTW: I believe that future versions of OpenFOAM might be switching to scotch instead of metis due to licensing aspects. I found that the scotch decomposition was fairly similar to metis anyhow.
Regarding licensing I think you may use Metis from 1.5 for testing purposes only. Otherwise, you need to have a license for that. Afaik this is the main reason why it changed to scotch in 1.6.
What I miss is a parallel partitioning method freeing us from the need to have machines with very large memory.

Regards
bastil is offline   Reply With Quote

Old   January 8, 2010, 05:38
Default
  #4
New Member
 
Join Date: Nov 2009
Posts: 17
Rep Power: 17
misakagan is on a distinguished road
Hi.

I'm using also scotch or metis since they do not require any manual input. I checked the subdomains in paraview and they seem to be convenient. But the problem is that I get a terrible speed-up (scaling) if I increase the processor number. In other words, nearly all computational effort is invested on solving the pressure-poisson equation in PISO loop, which does not scale well with increased processor number.

Are you satisfied with parallelisation of OpenFOAM? I'm asking because considering this computational performance I can not simulate anything big. Maybe I'm doing something wrong.
misakagan is offline   Reply With Quote

Old   January 8, 2010, 05:58
Default
  #5
Senior Member
 
Mark Olesen
Join Date: Mar 2009
Location: https://olesenm.github.io/
Posts: 1,714
Rep Power: 40
olesen has a spectacular aura aboutolesen has a spectacular aura about
Quote:
Originally Posted by bastil View Post
Regarding licensing I think you may use Metis from 1.5 for testing purposes only. Otherwise, you need to have a license for that.
Here https://projects.coin-or.org/BuildTo...Metis?rev=1294 it states
Quote:
For commercial-use, the essential restriction is for the vendors not to be reselling Metis.
The restriction appears to apply to the vendor, not the end-user - http://glaros.dtc.umn.edu/gkhome/met...faq#distribute

Last edited by olesen; January 8, 2010 at 06:07. Reason: added reference to faq
olesen is offline   Reply With Quote

Old   October 29, 2015, 07:13
Default
  #6
Senior Member
 
Join Date: Mar 2015
Posts: 250
Rep Power: 12
KateEisenhower is on a distinguished road
HI Mark,

can you please tell me how to load the separate processor domains into ParaView?

Best regards,

Kate
Quote:
Originally Posted by olesen View Post
Metis or scotch is what we use, since we generally don't have a better idea ourselves of how to decompose the domain any better (nor the time to figure it out anew for each geometry / cpu count combination).

Simply try out metis/scotch and see if you are happy with the results. You can view the separate processor* domains in paraview (eg, via paraFoam) to see if it matches up with your expectations.

BTW: I believe that future versions of OpenFOAM might be switching to scotch instead of metis due to licensing aspects. I found that the scotch decomposition was fairly similar to metis anyhow.
KateEisenhower is offline   Reply With Quote

Old   October 29, 2015, 10:00
Default
  #7
Member
 
Ron Burnett
Join Date: Feb 2013
Posts: 42
Rep Power: 13
rnburne is on a distinguished road
Nothing special is needed Kate. From within each processor file make the
VTK conversion, then open Paraview and load as you normally would.
Does this help?
KateEisenhower likes this.
rnburne is offline   Reply With Quote

Old   October 29, 2015, 11:41
Default
  #8
Senior Member
 
Join Date: Mar 2015
Posts: 250
Rep Power: 12
KateEisenhower is on a distinguished road
Sorry, somehow I overlooked the constant directory in the processor files.

Have a nice day,

Kate
KateEisenhower is offline   Reply With Quote

Old   August 15, 2017, 21:39
Default
  #9
New Member
 
Koushikc's Avatar
 
Koushik C
Join Date: Aug 2017
Posts: 1
Rep Power: 0
Koushikc is on a distinguished road
Quote:
Originally Posted by misakagan View Post
Hi.

I'm using also scotch or metis since they do not require any manual input. I checked the subdomains in paraview and they seem to be convenient. But the problem is that I get a terrible speed-up (scaling) if I increase the processor number. In other words, nearly all computational effort is invested on solving the pressure-poisson equation in PISO loop, which does not scale well with increased processor number.

Are you satisfied with parallelisation of OpenFOAM? I'm asking because considering this computational performance I can not simulate anything big. Maybe I'm doing something wrong.
Hi,
When it comes to parallel programming, there is a sweet spot between the no of processors you can use and your domain. Say your domain has 10000 nodes/points, then you have to find out what is the optimum no of processors you can add such that the communications are minimal. Vtune amplifier will be helpful in this analysis.
Hope this helps.

Regards,
Koushik

Sent from my MI 5 using CFD Online Forum mobile app
blttkgl and Ramzy1990 like this.
Koushikc is offline   Reply With Quote

Old   December 14, 2023, 03:36
Default
  #10
Member
 
MWRS
Join Date: Apr 2018
Location: Islamabad
Posts: 99
Rep Power: 8
waseeqsiddiqui is on a distinguished road
Hi everyone.
Just wanted to share that I was having trouble with parallel processing. I have found that aspect ratio is very crucial for partitioning interfaces.
waseeqsiddiqui is offline   Reply With Quote

Reply

Tags
domain decomposition, parallel computing


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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Any questions about Runge-Kutta methods Runge_Kutta Main CFD Forum 33 September 9, 2019 17:32
RANS modelling and URANS methods in unsteady flow simulations immortality OpenFOAM Running, Solving & CFD 31 January 3, 2014 08:19
OpenFOAM, Courant number and implicit methods fsaltara OpenFOAM 8 December 28, 2012 05:16
comments on FDM, FEM, FVM, SM, SEM, DSEM, BEM kenn Main CFD Forum 2 July 18, 2004 19:28
I just wonder why "SIMPLE" Junseok Kim Main CFD Forum 21 May 20, 2001 09:47


All times are GMT -4. The time now is 20:11.