[Raffaele10]

Hello everyone!

I'd like to set an automatic parametric-iterative case changing some parameters in the blockMeshDict. Basically I have a set of geometic parameters on which the mesh depens, what I want to do is to change this parameters and run the simulation automatically. This is a conjugate heat transfer case (liquid-solid) and I want to find the best shape of some component in order the minimize the temperture of the solid.

The idea is:
1-generate n sets of geometric parameters;
2-run n simulations changing the n set of parameters;
3-write a file with the results for each case.

I'm pretty new to OpenFOAM, python and C++ and I don't know which one is the best approach, I know that some use pyFoam, other use DAKOTA.

Can someone suggest a good approch to perform this case or maybe provide a tamplate with a similar case?

Thanks !

[Tobermory]

Sounds like a nice idea. The way I'd do it is as follows:

1. - generate a template case folder with a system/blockMeshDict file with the parameters as you described, and an Allrun script to mesh, setup the model and then run the case (check out the tutorials for an example of an Allrun script)
2. - generate a bash script file (this is basically what an Allrun script is) in the folder above; for each of the parameter sets you want to run, get the script to:
   • make a copy of the template folder, with a unique name
   • cd into that folder and run foamDictionary on blockMeshDict to update the parameter values
   • run the Allrun script, piping output into a local log file
   • execute any other postprocessing that you want to do (you should probably include that in the Allrun file)
3. use a separate log file for the bash script to log the actions that the script has completed

etc. - hopefully you get the idea. You could try run from python, but IMO bash scripts are much simpler to use. Good luck!

[dlahaye]

https://github.com/damogranlabs/classy_blocks ?

[HYJLEE]

Using .m4 file to make your parametrized blockMesh template would be helpful.

This page contains information of m4 file.
https://openfoamwiki.net/index.php/H...ckMesh_with_m4
refer to this page, you can make your own .m4 file

And then, modifying only 'value' of define(variable,value) in .m4 file automatically with some code(it can be based on C++, python, etc..)

If you successfully editted your .m4 file, command "m4 blockMeshDict.m4 > blockMeshDict;blockMesh" will construct mesh with changed parameters.

I hope it would be helpful to find your optimal solution using OpenFOAM.

ps. There are some examples on .m4 file in tutorials of foam-extend. Its m4 file is also can be used not only foam-extend but also OpenFOAM, when its boundary conditions are modified properlly.