In the fourth of a series of posts on technical tips, I will look at mesh generation using The NAG Library. Mesh generation is essential when using the finite element method to solve a system of partial differential equation (PDE). PDEs, where this method is appropriate, can be used, for example, to solve problems in solid mechanics, fluid mechanics and thermal modelling. Using a meshing algorithm is indispensable when the geometry of the numerical domain is not trivial, i.e. not rectangular. Meshing ensures that the spatial discretization of the domain reflects the geometry of the numerical domain.
The NAG Library has several routines in the d06 Chapter to define the boundary, generate the mesh, and the solve the PDEs on the mesh. This technical report illustrates how these routines can be used in conjunction with a sparse iterative solver from Chapter f11 to solve the PDEs using the finite element method.
NAG would like to acknowledge the Institut National de Recherche en Informatique et Automatique, who’s developed the MODULEF software from which the NAG derived the routines in the Chapter.