[Other] Creating large meshes (>10,000,00) in parallel

fontania · April 16, 2012, 14:01

Fellow members and users of OpenFOAM,

I am having problems creating large meshes (greater than 10,000,000) elements. No surprise I am running out of memory on a single computer. I am looking to create a mesh up to 100,000,000 elements. The domain is simply a unit 1x1x1 cube. Does anyone have a tested procedure for parallel discretization or refinement?

Thanks in advance!

lovecraft22 · April 16, 2012, 15:26

How can your domain be a 1x1x1 cube?
Have you tried with decomposePar and then running sHM in parallel with mpirun?

fontania · April 16, 2012, 16:04

1x1x1 cube: I am doing a scaling analysis on a cluster to understand how many processors to run different problem sizes. My base case is the incompressible/icoFOAM/cavity tutorial.

The problem I will eventually be solving is airflow within a large building. After the scaling analysis do I...

1. Create a STL file of the domain
2. Use blockMesh for a basic overlapped representation.
3. Use decomposePar to decompose the initial mesh to the number of processors I plan on running the simulation
4. Use snappyHexMesh in parallel to create the actual domain discretization.

Thanks for the quick reply.

lovecraft22 · April 16, 2012, 16:29

Sorry, I'm missing the point then… The question is…?

wyldckat · April 16, 2012, 17:14

Quick answer: see the tutorial "./incompressible/simpleFoam/turbineSiting"

April 16, 2012, 14:01	Creating large meshes (>10,000,00) in parallel	#1
fontania New Member Join Date: Apr 2012 Location: Ames, IA Posts: 6 Rep Power: 14	Fellow members and users of OpenFOAM, I am having problems creating large meshes (greater than 10,000,000) elements. No surprise I am running out of memory on a single computer. I am looking to create a mesh up to 100,000,000 elements. The domain is simply a unit 1x1x1 cube. Does anyone have a tested procedure for parallel discretization or refinement? Thanks in advance!

April 16, 2012, 17:14		#5
wyldckat Retired Super Moderator Bruno Santos Join Date: Mar 2009 Location: Lisbon, Portugal Posts: 10,981 Blog Entries: 45 Rep Power: 128	Quick answer: see the tutorial "./incompressible/simpleFoam/turbineSiting" __________________ OpenFOAM: FAQ \| Getting started Forum: How to get help, to post code/output and forum guide Read this before sending me PM

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
Error running simpleFoam in parallel	Yuby	OpenFOAM Running, Solving & CFD	14	October 7, 2021 05:38
parallel Grief: BoundaryFields ok in single CPU but NOT in Parallel	JR22	OpenFOAM Running, Solving & CFD	2	April 19, 2013 17:49
Creating large enclosure containing many persons	Michael	CFX	0	January 21, 2009 06:42
[blockMesh] Tutorial for Creating Simple 3D Meshes	jonlemur	OpenFOAM Meshing & Mesh Conversion	4	August 8, 2007 14:20
Large 3D tetrahedral meshes	Aldo Bonfiglioli	Main CFD Forum	4	August 27, 1999 04:33

April 16, 2012, 15:26		#2
lovecraft22 Senior Member lore Join Date: Mar 2010 Location: Italy Posts: 460 Rep Power: 18	How can your domain be a 1x1x1 cube? Have you tried with decomposePar and then running sHM in parallel with mpirun?

April 16, 2012, 16:04		#3
fontania New Member Join Date: Apr 2012 Location: Ames, IA Posts: 6 Rep Power: 14	1x1x1 cube: I am doing a scaling analysis on a cluster to understand how many processors to run different problem sizes. My base case is the incompressible/icoFOAM/cavity tutorial. The problem I will eventually be solving is airflow within a large building. After the scaling analysis do I... 1. Create a STL file of the domain 2. Use blockMesh for a basic overlapped representation. 3. Use decomposePar to decompose the initial mesh to the number of processors I plan on running the simulation 4. Use snappyHexMesh in parallel to create the actual domain discretization. Thanks for the quick reply.

April 16, 2012, 16:29		#4
lovecraft22 Senior Member lore Join Date: Mar 2010 Location: Italy Posts: 460 Rep Power: 18	Sorry, I'm missing the point then… The question is…?