Daniel Wei's CFD Blog: (http://lakeat.co.nf/)

I have made significant improvement on the mesh export. The manual building up of OpenFOAM mesh now turns up really fast. I got a feeling that it could possibly beat even OpenFOAM's native blockMesh generator. In my test of a 3,200,000 mesh (for flow around a square cylinder). The building process took 90.51s in total. Here are more details:

...... [1] Done assembling the points. ( 7.09000s)
...... [2] Done assembling all cells & intnl cntrs. ( 40.62000s)...

Brief updates: I have succeeded optimizing some pieces of the code, from O(n^2) to O(n), but I still have another key module having the O(nlogn) complexity. It is unacceptable. I need to get rid of it.

Updates:

Finally, I got the boundary file output working. Now it is able to export to the OpenFOAM mesh format. But I still need to do some optimization (hopefully to make it back to log-complexity with hash tables, thanks to c++11). Currently it takes too much time to build an OpenFOAM mesh manually.

Also I have tested the code on a clean ubuntu 64 environment with virtualbox, it seems to me anyone who wants to use it needs vtk and qt installed. Using synaptic package...

Spent four days in figuring out how to construct an OpenFOAM polyMesh manually, it is more complicated than I thought, the difficult part is how to construct a multi-block mesh. One of the annoying issues in a multi-block polymesh construction is how to handle its internal lines/surfaces, where two lines become one. And then the points indices and the face indices all have to change accordingly. There may be many approaches to address this, but I, after many days' thought, decided to use a brutal...

Trust me, I've read more than 20 structured mesh generator codes already in the last year. Interestingly, none of them satisfy me. So I decided to make it my own. I've been working on this small mesh generator for quite a while. It worked before for a single block mesh smoothing. But starting from last month, as I was thinking about how to make it "multi-block", I decided to rewrite the whole code. Because it is quite different from C++ point of view to design classes for a single block...