CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Pre-Processing

Problem with decomposePar (and mapFields) for large problem

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   February 7, 2019, 04:48
Default Problem with decomposePar (and mapFields) for large problem
  #1
Senior Member
 
TWB
Join Date: Mar 2009
Posts: 414
Rep Power: 19
quarkz is on a distinguished road
Hi,

I'm having problem using decomposePar (and mapFields) for large problem.

My grid has abt 0.28 billion cells. I intend to run on a cluster with 1000 - 1500 cpus. It is a DES simulation and I need to use resumed data from a smaller problem, mapField, and then decomposePar.

However, I have some problems.

1. I used to run:

mapFields ../path/to/source -consistent -sourceTime latestTime

then

decomposePar -force

But maybe because the problem is too big, mapFields quit after about 0.5 hr.

2. Another way is to 1st create a folder with dummy files in the time step I want to map, then running

decomposePar -force

then

mapFields path/to/source -sourceTime latestTime -parallelSource -parallelTarget -consistent

or

mpirun -n xx mapFieldsPar -consistent -parallel -sourceTime latestTime path/to/source

However decomposePar also quit after 0.5 hr.

So what other solutions do I have?

Thanks
quarkz is offline   Reply With Quote

Old   February 20, 2019, 10:18
Default
  #2
Senior Member
 
Zander Meiring
Join Date: Jul 2018
Posts: 125
Rep Power: 8
yambanshee is on a distinguished road
Why not make use of the cluster to do the map fields in parallel? Split both domains up, and use mpirun mapFieldsPar
yambanshee is offline   Reply With Quote

Old   February 21, 2019, 10:51
Default
  #3
Senior Member
 
TWB
Join Date: Mar 2009
Posts: 414
Rep Power: 19
quarkz is on a distinguished road
Hi yambanshee,

Thanks for your suggestions. I have tried but using this way, I got the error:

[1021]
[1021]
[1021] --> FOAM FATAL ERROR:
[1021] Unable to set source and target faces
[1021]
[1021] From function void Foam::faceAreaWeightAMI<SourcePatch, TargetPatch>::setNextFaces(int &, int &, int &, const Foam::List<bool> &, Foam::List<int> &, const Foam:ynamicList<int, 0U, 2U, 1U> &, bool) const [with SourcePatch = Foam::PrimitivePatch<Foam::face, Foam::SubList, const Foam::Field<Foam::Vector<double>> &, Foam::Vector<double>>, TargetPatch = Foam::PrimitivePatch<Foam::face, Foam::SubList, const Foam::Field<Foam::Vector<double>> &, Foam::Vector<double>>]
[1021] in file lnInclude/faceAreaWeightAMI.C at line 287.

Hence, I'm not able to use this mtd. Or do you know how to solve this problem?

Thanks.
quarkz is offline   Reply With Quote

Reply


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
Problems with decomposePar and mapFields Fischlurch OpenFOAM Running, Solving & CFD 2 October 7, 2016 08:51
implementation of mapFields into parallel transient case simpomann OpenFOAM Pre-Processing 4 August 2, 2016 05:41
transientSimpleDyMFoam, mapFields and decomposePar pad OpenFOAM Running, Solving & CFD 0 December 3, 2010 06:22


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