Numerical simulation methods like finite element or finite volume methods require a mesh, which provides the discretization of the geometry of the domain as well as the internal points, and elements in which the solution variables (velocities, electrical potential, ion concentrations...) will be computed.
The present result concerns a software package, which generates such grids suitable for computations of mass and heat transfer with electrochemically reacting flow solvers. Typical in these applications is the presence of extremely thin diffusion layers near reacting electrodes, which require anisotropic cells with aspect ratios of several thousands. This means that the spacing of the cells normal to the wall needs to be 3 to 4 orders of magnitude smaller than the spacings tangential to the wall. No commercial package exists which is able to achieve these specifications for the kind of complex geometries encountered in electro-deposition processes.
The method followed is a hybrid grid generation approach which combines the strengths of structured grid generators to produce stretched anisotropic layers of elements near the boundary with the capacities of unstructured grid methods to cope with complex geometries. The procedure is hierarchic, beginning with the construction of the edge grids, which constitute the boundaries of the faces. In a second step the face grids are generated starting from the bounding edge grids, and finally in a third step the volume grids are generated starting from the bounding face grids. A strong emphasis is given to efficient interfacing with the CAD system. The geometry and topology of the configuration can be taken directly from the output of the CAD system, provided that the CAD definition is exported under the form of a STEP file. The STEP format is a recently developed ISO standard for the electronic exchange of product data, supported by most CAD packages.
Extremely flexible user control of the mesh spacing is possible, allowing to construct either isotropic or boundary layer grids adjacent to each topological entity. The spacing parameters (e.g. stretching, growth rate) can be controlled by the user for each topological entity separately (vertex, edge, surface or brep). For the isotropic volume meshing, two different point placement algorithms have been implemented, while a Delaunay triangulation method is used to generate the tetrahedra from a given point cloud. For the anisotropic boundary layer part, semi-structured layers of prismatic elements are created normal to the surface, starting from the triangular surface mesh. When the prismatic elements become isotropic the layer generation is stopped and the remaining domain is filled with isotropic tetrahedra (eventually after the creation of some intermediate pyramidal elements).
To simplify the specification of input parameters for the user a novel method was developed allowing the generation of high quality surface meshes using only a few user parameters. The method is based on the surface curvature to accurately resolve the geometry features and uses the topological model provided in the STEP format to create smooth surface grids on complex CAD models with hundreds of surfaces.
The grid generation package has been developed in close collaboration with the software company ElSyCa, who was responsible for Graphical User Interfacing (GUI), CAD interfacing, visualization and interfacing with solvers. Therefore the present grid generator is part of the prototype of the CAD integrated software developed by ELSYCA (see separate result).
Apart from the applications within the context of this project on electrochemical modelling, VKI also intends to exploit the grid generation software for aeronautical applications at high Reynolds number, for which similar though less severe requirements exist, e.g. for the computation of drag around a wing or aircraft body.