[shaileshbg]

Hello Foamers,

Newbie to python and pyFoam, I am looking to automate my simulation using pyFoam libraries.
Problem Statement: I have two different bodies which I have meshed but I want to merge and then stitch the two meshes.

The meshing of my two regions muscle and tissue is successful:

# Muscle Meshing
orig = SolutionDirectory(path.expandvars("muscle"), archive=None, paraviewLink=False)
work = orig.cloneCase("muscle1")

blockRun = BasicRunner(argv=["blockMesh", "-case", work.name], silent=True, server=False, writeState=False)
blockRun.start()
if not blockRun.runOK():
    error("There was a problem with blockMesh")

snapRun = BasicRunner(argv=["snappyHexMesh", "-case", work.name], silent=True, server=False, writeState=False)
snapRun.start()
if not snapRun.runOK():
    error("There was a problem with snappyHexMesh")

# Tissue Meshing
orig = SolutionDirectory(path.expandvars("tissue"), archive=None, paraviewLink=False)
work = orig.cloneCase("tissue1")

blockRun = BasicRunner(argv=["blockMesh", "-case", work.name], silent=True, server=False, writeState=False)
blockRun.start()
if not blockRun.runOK():
    error("There was a problem with blockMesh")

snapRun = BasicRunner(argv=["snappyHexMesh", "-case", work.name], silent=True, server=False, writeState=False)
snapRun.start()
if not snapRun.runOK():
    error("There was a problem with snappyHexMesh")

The bash equivalent of what I want to do:

runApplication mergeMeshes -overwrite . ../tissue/              # Merge tissue mesh onto the muscle mesh
runApplication stitchMesh -partial -toleranceDict toleranceDict -overwrite muscle_ext muscle_int   # partial stitch muscle patches per tolerance
cd 0
rm  meshPhi                           # remove meshPhi
cd ..
rm log.*                             # remove current merge and stitch logs
cd constant/polyMesh/                # move to constant polymesh to remove current zones
rm *Zones                            # remove all current zone files
rm meshModifiers                     # remove mesh modifiers

Any help would be greatly appreciated.

[gschaider]

Originally Posted by shaileshbg

Hello Foamers,

Newbie to python and pyFoam, I am looking to automate my simulation using pyFoam libraries.
Problem Statement: I have two different bodies which I have meshed but I want to merge and then stitch the two meshes.

The meshing of my two regions muscle and tissue is successful:

# Muscle Meshing
orig = SolutionDirectory(path.expandvars("muscle"), archive=None, paraviewLink=False)
work = orig.cloneCase("muscle1")

blockRun = BasicRunner(argv=["blockMesh", "-case", work.name], silent=True, server=False, writeState=False)
blockRun.start()
if not blockRun.runOK():
    error("There was a problem with blockMesh")

snapRun = BasicRunner(argv=["snappyHexMesh", "-case", work.name], silent=True, server=False, writeState=False)
snapRun.start()
if not snapRun.runOK():
    error("There was a problem with snappyHexMesh")

# Tissue Meshing
orig = SolutionDirectory(path.expandvars("tissue"), archive=None, paraviewLink=False)
work = orig.cloneCase("tissue1")

blockRun = BasicRunner(argv=["blockMesh", "-case", work.name], silent=True, server=False, writeState=False)
blockRun.start()
if not blockRun.runOK():
    error("There was a problem with blockMesh")

snapRun = BasicRunner(argv=["snappyHexMesh", "-case", work.name], silent=True, server=False, writeState=False)
snapRun.start()
if not snapRun.runOK():
    error("There was a problem with snappyHexMesh")

The bash equivalent of what I want to do:

runApplication mergeMeshes -overwrite . ../tissue/              # Merge tissue mesh onto the muscle mesh
runApplication stitchMesh -partial -toleranceDict toleranceDict -overwrite muscle_ext muscle_int   # partial stitch muscle patches per tolerance
cd 0
rm  meshPhi                           # remove meshPhi
cd ..
rm log.*                             # remove current merge and stitch logs
cd constant/polyMesh/                # move to constant polymesh to remove current zones
rm *Zones                            # remove all current zone files
rm meshModifiers                     # remove mesh modifiers

Any help would be greatly appreciated.

You mean how to go on from the stuff you already did in PyFoam?um First problem I see is that you assigned both SolutionDirectories to the same variable: work. So lets assume you assigned them to muscleWork and tissueWork. Then something like

mergeRun=BasicRunner(argv=['mergeMeshes','-overwrite',muscleWork.name,tissueWork.name])

might do the trick (not sure whether the -case option is needed with this utility). The shell commands you can either encapsulate in os.system or emulate with os.unlink

[shaileshbg]

Hi Bernhard,

Thank you very much for your reply, this is working perfectly.
Alternatively, I am using subprocess to call in the bash functions directly in python.