CFD Domain's Creation & Discretization

The next step for developing the CFD simulation for a case study, after the review of physics and theory behind the problem of interest, is to create the CFD/Computational domain. This CFD domain will first be created and then discretized into small finite elements (2D) or finite volumes (3D). The flow field's governing equations are then solved numerically inside these elements or volumes via an iterative process till a converged solution inside the CFD domain is obtained. This section explains the techniques and process of creation and discretization of the CFD domain for the problem of interest. It should be highlighted that this sections provides the most generic steps and techniques to create the CFD domain. Users can implement these steps/techniques in the mesher of their choice or download the pre-generated * .msh file from the "Simulation_Src" folder.

The CFD domain for the 2D laminar flow over a flat plate is chosen to be a finite two dimensional square with sides of 1 [m] as shown in Fig.1. The shape and dimensions of the domain are chosen such that the required number of mesh elements can be controlled and non-dimensional groups of variables can be estimated fast, while most of variables are set to be equal unity. The users can define the dimensions of the CFD domain according to their own problem definition.

Fig.1 - The geometry of CFD domain for the two dimensional laminar flow over a flat plate.

To create the geometry of the visualized CFD domain in Fig.1, it is recommended to start with four points spaced equally with distance of 1 [m] with respect to each other. Then connect them together and define a unit, closed surface. This procedure might need to be implemented in a slightly different manner according to the used CAD/mesher software of choice.

Once the CFD domain is created, it should be discretized into finite elements (a.k.a meshing procedure). In simple words, meshing is the process of "chopping" the geometry of the CFD domain into smaller finite elements and/or volumes, wherein the flow field's governing equations will be solved numerically. High quality mesh elements, with reasonable aspect ratio, would lead to a smoother and more promising numerical solution. Therefore, it is important to increase the quality and level of controllability in the CFD domain at the first step by proper meshing and topology techniques.

In the current problem of interest to increase the level of controllability on mesh resolution, the vertical edges of the CFD domain are discretized into 100 mesh elements starting with a first length of 1 [mm] from the bottom horizontal edge. The mesh element distribution along these edges follow a predefined distribution as shown in Fig. 2. The predefined distribution would control the mesh resolution along the edge from the bottom to top of the domain. The rule of thumb is to have a high mesh resolution mesh in the regions where the flow field behavior is complex and decrease mesh resolution where the flow field behavior is simple. For very complex flow field such as turbulent boundary layer there are specific techniques and guidelines to set the mesh resolution, which will be introduced in other case studies as needed. Therefore, users need to provide a high mesh resolution close to flat plate boundary, where the physics of the flow are more complex (i.e. boundary layer region) and gradually decrease it to the top of the domain where flow field is uniform and simple. This mesh technique provide users the capability to modify the mesh later for more complex flows with minimum amount of effort and time.

Fig.2 - Controlling mesh element distribution along the vertical edges of the CFD domain.

In the next step the horizontal edges of the CFD domain will be meshed. For these edges a uniform distribution of mesh elements along top and bottom edges are used as visualized in Fig.3. It should be noted that if required the users can define specific mesh element distribution to increase or decrease the mesh resolution at the inlet, middle or outlet of the domain. For the current problem of the interest these specific considerations are not necessary.

Fig.3 - Uniform mesh distribution along the horizontal edges of the CFD domain.

Finally the two dimensional CFD domain can be meshed using the proper meshing algorithm of choice via connecting the previously generated mesh elements to each other. Fig. 4 visualizes the outcome of this meshing progress. As shown in Fig.4 the mesh resolution is higher close to the boundary of the flat plat and it decreases moving up toward the top part of the CFD domain.

Fig.4 - Finalized mesh of the CFD domain.

Once the mesh is generated it's quality should be examined. This can be done using the provided features and functions defined in the meshing software. In general these functions calculate the aspect ratio of each elements and report a statistical report on them. The goal is to keep the average aspect ratio close to 1 for a very high quality mesh. In this stage the CFD domain is completely meshed and is ready to be setup for a CFD simulation.


In [ ]: