[kk415]

Hello Foamers,


I am working on a project in which I need to develop a multiphase solver using the concept of dual grid (i.e. using two grid simultaneously, one for flow solver and another for volume fraction scalar. The key idea is to use fine mesh for volume fraction scalar).



So I develop the solver in a similar manner to chtMultiRegion solver in which two regions are overlapping each other. The solver works for serial run but for parallel case it is causing problem.



I have the following question to ask:


1. How can I ensure that after decomposition each processor takes the same set of cells in both regions?
2. Is there any way to create a refine mesh from the parent mesh using splitcells or something, inside a solver?

[kk415]

Working on splitting of mesh I found this document very helpful.


http://www.tfd.chalmers.se/~hani/kur...erceReport.pdf


Now I am working in similar manner, and splitting the cells the following manner.


#include "createMesh.H"
#include "createNewMesh.H" //Create object newMesh

const cellList& cells = mesh.cells(); //List to loop through
labelList refCells; //List of cells to refine
refCells.clear(); //Clear the list
refCells.resize(mesh.cells().size());
forAll(cells,cellI)
{
refCells[cellI]=cellI;
}
Info << "Entering Refinement" << endl;
//******** Cell Refinement **********
// Since this is a pretty long function i would in this section
// like to call upon a function file that refines the mesh.
// That file should be similar to refineMesh.C
// Define refinement dictionary
dictionary refineDict; //Declare
dictionary coeffsDict;
coeffsDict.add("tan1", vector(1, 0, 0));
coeffsDict.add("tan2", vector(0, 1, 0));
wordList directions(2);
directions[0] = "tan1";
directions[1] = "tan2";
refineDict.add("directions", directions); // Add directions to the dictionary
refineDict.add("useHexTopology", "false");// Use standard cutter
refineDict.add("coordinateSystem", "global");
refineDict.add("globalCoeffs", coeffsDict);
refineDict.add("geometricCut", "false");
refineDict.add("writeMesh", "false");
// Multi-directional refinement (does multiple iterations)
Info << "Entering Refinement" << endl;
multiDirRefinement multiRef(newMesh, refCells, refineDict);
Info << "Done with Refinement" << endl;
Info << "Writing new mesh with " << refCells.size() << " refinements" << endl;
//mesh.write();
newMesh.write();
Info << newMesh.cells().size() << endl;


The problem is my parent mesh is also getting refined, when I want only newMesh to be refined.