The code presented is a parallel solver for levelset topology optimization method with adaptive mesh refinement
A level set topology optimization solver without adaptive mesh refinement is also uploaded. Users could run this solver on both OpenFOAM 5.0 and 6.0.
Before running this solver, following softwares are needed.
(a) OpenFOAM 6.0
Note that Openfoam only support adaptive mesh refinement for 3D mesh !!! For 2D optimization cases, users should compile the code in folder "meshRefinement2D_of6" and link the dynamic library "libdynamicFvMeshUser.so" (see app/system/controlDict).
Thanks for Luca Cornolti who provided the source code for 2D adaptive mesh refinement solver(https://www.cfd-online.com/Forums/openfoam-community-contributions/118870-2d-adaptive-mesh-refinement-2.html#post703369).
After finishing above works, users should run wmake in the src folder, and then run blockMesh and decomposePar in the app folder. Finally, topology optimization solver can run in the app folder by mpirun –n 16 MTO_levelset, where 16 is the default number of cores for parallel computing.
After the optimization, users should run reconstructPar in the app folder and then run paraFoam to view the results with Paraview.
