CFD simulations involving 3D complex geometry have become the norm, however, this hasn’t lessened the…
This monthly mesh we refer readers to our guest blog “It’s all in the numbering – mesh renumbering may improve simulation speed”, published on Another Fine Mesh.
We all know that the mesh plays a vital role in CFD simulations. Yet, not many realize that renumbering (ordering) of the cells in the Finite Volume Method (FVM) can affect the performance of the linear solver and thus the speed of the simulation. In some circumstances, a poorly numbered mesh may result in solver divergence and failure of the simulation. In the guest blog, we use Pointwise and OpenFOAM® to illustrate the effect of mesh renumbering on the matrix structure.
Why renumber mesh?
The goal of mesh renumbering is to reduce the matrix bandwidth and profile. There are many sparse matrix renumbering strategies, with the most commonly used in CFD being the Cuthill-McKee.
Most CFD solvers either renumber the mesh without user intervention or provide an application for the user to renumber the mesh manually. OpenFOAM has an application called renumberMesh that reduces the bandwidth and profile of the coefficient matrices by renumbering the cell label list.
Pointwise stores the mesh internally in a format and numbering that is efficient for mesh generation. This numbering does not necessarily result in a low bandwidth or profile matrix upon export.
Mesh renumbering structured mesh
This case is a simple structured mesh, as shown in Figure 1. The initial bandwidth and profile reported from renumberMesh (OpenFOAM v2006) were 7 and 300, respectively.
Plotting the matrix sparsity pattern (nonzero entries) allows visual inspection of the effect of matrix bandwidth and profile reduction (see Figure 2). A visualisation of the matrix sparsity for the simple structured mesh example before renumbering exhibits diagonal banding.
Figure 3 shows the sparsity pattern for the renumbered simple structure mesh after the application of the Cuthill-McKee renumbering method. In this case, the matrix bandwidth remains unchanged, but the profile was reduced by 18.3% to 245. This mesh renumbering may improve simulation speed.
We have shown the importance of renumbering using meshes exported from Pointwise. Meshes generated natively in OpenFOAM, from applications such as blockMesh, snappyHexMesh and cfMesh should also be renumbered.
Want to give it a go yourself?
Request a Pointwise free trial see link below.
In case you missed our other blogs, see here.
Thanks for reading, the next edition will be in November. Until then, happy meshing!
18. It’s all in the numbering – mesh renumbering may improve simulation speed