Next Article in Journal
An Image Encryption Algorithm Based on Complex Network Scrambling and Multi-Directional Diffusion
Previous Article in Journal
Research on Throughput-Guaranteed MAC Scheduling Policies in Wireless Networks
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Neural Network-Based Mesh Quality Indicator for Three-Dimensional Cylinder Modelling

1
Science and Technology on Parallel and Distributed Processing Laboratory, National University of Defense Technology, Changsha 410073, China
2
Laboratory of Software Engineering for Complex System, National University of Defense Technology, Changsha 410073, China
3
China Aerodynamics Research and Development Center, Mianyang 621000, China
*
Authors to whom correspondence should be addressed.
Entropy 2022, 24(9), 1245; https://doi.org/10.3390/e24091245
Submission received: 31 July 2022 / Revised: 26 August 2022 / Accepted: 2 September 2022 / Published: 4 September 2022
(This article belongs to the Section Multidisciplinary Applications)

Abstract

:
Evaluating mesh quality prior to performing the computational fluid dynamics (CFD) simulation is an essential step to ensure the acceptable accuracy of cylinder modelling. However, traditional mesh quality indicators are often insufficient since they only check geometric information on individual distorted elements. To yield more accurate results, the current evaluation process usually requires careful manual re-evaluation for quality properties such as mesh distribution and local refinement, which heavily increase the meshing overhead. In this paper, we introduce an efficient quality indicator for varisized cylinder meshes, consisting of a mesh pre-processing method and a neural network-based indicator, Mesh-Net. We also publish a cylinder mesh benchmark dataset. The proposed indicator is trained to study the role of CFD meshes on the accuracy of numerical simulations. It considers both the effect of element geometry (e.g., orthogonality) and quality properties (e.g., smoothness and distribution). Thereafter, the well-trained indicator is used as a black-box to predict the overall quality of the input mesh automatically. Experimental results demonstrate that the proposed indicator is accurate and can be applied in the mesh quality evaluation process without manual interactions.

1. Introduction

Computational fluid dynamics (CFD) plays a vital role in a broad spectrum of scientific and engineering fields, such as bioengineering, aerospace, energy engineering, and manufacturing [1,2,3]. During the CFD simulation, the quality of the generated mesh directly influences the solution accuracy and error magnitude. Many mesh generation methods have been proposed aiming to generate high-quality meshes [4,5]. Unfortunately, the quality of the initial mesh is usually not acceptable. The minimal mesh quality requirement is seldom achieved except on the most elementary problems [6,7,8]. Therefore, the procedure used to handle high-quality mesh generation is divided into three steps: initial mesh generation, mesh quality evaluation, and mesh optimisation. In this meshing process, an efficient mesh quality indicator is particularly important. The indicator determines the direction of subsequent quality optimization and ensures the accuracy of the desired solution. It serves as a basis for assessing the ability of generated mesh to faithfully represent the physics of the flow.
The high degree of complexity (non-linearity) between mesh quality and numerical accuracy makes the quality evaluation an extremely difficult task. It is hard to precisely define the relationship between mesh qualities and their correlations with numerical error [9]. Starting from the observation that regular or equilateral mesh elements are more pleasing, the traditional evaluation procedures focus on evaluating the shape of each element. Such a perspective leads to the formulation of quality indicators in terms of elemental geometric information such as area ratio, edge length, volume ratio, aspect ratio, skewness, minimum (maximum) angle, and gamma coefficient [10,11,12]. However, the above geometric-based indicators are often insufficient since: (1) They only yield geometric information on individual distorted mesh elements, and are useless for quality properties such as mesh distribution and local refinement. These properties affect the simulation accuracy in regions near boundary layers and wing-body configurations. (2) They may give inconsistent evaluation results for the same mesh (see [13,14] and references therein). Moreover, since geometric-based indicators may not guarantee an accurate result, the issue of mesh quality evaluation usually requires careful manual re-evaluation. This process relies heavily on the empirical, descriptive realm of a priori knowledge. As a result, the frequent human-computer interactions needed in the current evaluation process have become a bottleneck to the fully-automatic meshing process and significantly increase the meshing overhead. In order to ensure the cost-efficiency of meshing, it is essential to build an intelligent mesh quality indicator without manual interactions.
In recent years, artificial neural networks have been proven capable of learning complex mapping and replacing human labour in various applications [15,16,17]. The network utilises multiple layers of neural units to learn important features automatically from high-dimensional parameter spaces. By performing an optimisation procedure based on the loss function, the network model is able to approximate the complex and nonlinear mapping from training samples [18]. Despite the widespread success of neural networks in various physical problems, there have been only limited attempts at neural network-based mesh quality evaluation.
In this paper, we propose a mesh quality indicator for three-dimensional cylinders, resulting in a point-based mesh pre-processing method, a neural network Mesh-Net, and a cylinder mesh benchmark dataset. The proposed indicator takes mesh files as input and learns the potential correlations between the mesh quality and simulation error. Compared with traditional quality indicators, which focus on detecting distorted mesh elements, our indicator is more accurate. It considers both the effect of element geometry, such as orthogonality, and quality properties, such as smoothness and mesh distribution. Experimental results demonstrate that the well-trained indicator is able to predict the overall quality of cylinder meshes and achieves an accuracy of up to 98.05%. Moreover, it can be applied in the automatic mesh quality evaluation process without manual interactions, which significantly reduces the meshing overhead. We hope our work can provide future research directions that contribute to efficient mesh generation technology. The proposed benchmark dataset is publicly available at https://github.com/MeshDataset/3D-Cylinder (accessed on 4 October 2021).
The rest of the paper is organized as follows: Section 2 describes related works about existing mesh quality indicators. Section 3 gives details of the proposed neural network-based mesh quality indicator. The experimental results and discussion are presented in Section 4. The conclusion is finally outlined in Section 5.

2. Related Works

2.1. Traditional Mesh Quality Indicators

It is well known that poorly shaped meshes tend to slow convergence and cause instability during the CFD simulation [19,20]. In order to ensure the accuracy of the numerical solution, many indicators have been proposed to check the mesh quality before simulation.
Starting from the observation that regular or equilateral mesh elements are more pleasing, Strang and Fix [21] discussed the minimum angle condition of mesh elements. They stated that the smallest angle of mesh elements should be bounded away from zero. Berzins [6] supported this view and proved that elements with a relatively small included angle might have a negative effect on the solution of the linear algebra problem. Similar conclusions were proposed by Shewchuk [22], who showed that a prerequisite for high-quality mesh elements is that there should be no large included angles.
To make the conclusion more specific, Liu and Joe [23] proposed a quality indicator Q L to identify ‘sliver’ elements.
Q L = 1 8.48528 V l i 6 3 ,
where V is the volume, and l i is the edge length of the examined mesh element.
Bank [24] presented a geometric indicator Q B for CFD mesh quality control:
Q B = 3 12 A l i 2 ,
where A is the area of the mesh element. Weatherill [25] introduced a similar mesh quality indicator in the evaluation process. It is defined by:
Q W = 1 3 A l i 2 ,
Another indicator referred to as the Scaled Jacobian Quality Indicator was proposed in the CUBIT code for the mesh quality [5]. The Scaled Jacobian first computes the triple product at each node of the element corners using the other mesh nodes. It then computes the average of the corner Jacobians. The value of this indicator varies from minus one to plus one. A positive scaled Jacobian is usually considered the minimum quality for an acceptable computational mesh (called inversion-free). In contrast, the negative values of the Scaled Jacobian indicate the presence of distorted elements.
Quality indicators such as Aspect Ratio, Diagonal Ratio, Edge Ratio, and Equiangle Skewness are widely-used in CAE software as quality metrics for mesh elements [26]. For example, The Diagonal Ratio  Q D R is represented as the maximum ratio of the element diagonals:
Q D R = m a x ( d 1 , d 1 d n ) m i n ( d 1 , d 1 d n ) ,
where d i is the length of the element diagonal. By definition, the higher the metric value, the less regularly shaped the examined element. For equilateral elements (square quadrilateral elements or cubic hexahedra), the Diagonal Ratio Q D R is 1.
The above quality metrics provide shape specifications for mesh elements employing geometric formulas (the value usually ranges between 0 and 1, and 1 for an equilateral element). However, several sets of numerical results in [13] have demonstrated that employing different quality metrics to evaluate the same element may lead to inconsistent results. This conclusion is also confirmed by Gao et al. [14]. They performed a thorough numerical study to analyse widely-used quality indicators and their correlations with the stability and accuracy of the simulations. Nearly twenty quality indicators were tested on hexahedral elements. It was observed that the correlations among indicators are ambiguous. The derivation of some geometric element-based quality indicators applies only to specific applications.
Overall, present-day mesh quality indicators tend to assess geometric imperfections (shape, edge length, included angle, Jacobian) on mesh elements. Other considerations such as mesh density and distribution that ensure desirable simulation accuracy are ignored. This deficiency imposes a burden on careful manual re-evaluation, which significantly increases the meshing overhead [9,27].

2.2. Neural Networks for Mesh Quality Evaluation

In recent years, many researchers attempt to explore new methods for complex physical problems using artificial neural networks (ANNs). The main insight of ANNs is the capability of finding nonlinear approximations to complex functions based on the architecture of interconnected neurons. After suitable training, the ANNs are able to predict the desired output accurately. As a result, ANNs have been successfully applied to various CFD problems to improve efficiency and reduce overhead [9,18].
Multi-layer perceptron (MLP), which consists of several layers of neurons, is a specific type of ANNs [28]. The structure of an MLP is separated into three parts, the input layer, hidden layers, and output layer. A multi-layer perceptron with three hidden layers is shown in Figure 1. The neurons in the first hidden layer receives source signals from the input layer and propagate them to the succeeding layers. The signals are passed between all hidden layers (with activation functions) and finally converted into high-level features. The feature values in the output layer indicate the probability of the input belonging to a particular category. In this forward propagation process, the output of each layer is computed as:
N ( x ) = i = 1 m σ   j = 1 n ( ω i j x j + b i ) ,
where n is the number of neural units in the hidden layer. m denotes the number of input units. ω is the weight, and b is the bias. The activation function is represented by σ . The loss function concludes the partial derivatives of the layer outputs with respect to the variables. After that, the adjustable variables in the neurons (weights and biases) are optimised via backpropagation to approximate the nonlinear mapping.
To better learn the local and contextual information from input data, convolutional neural networks (CNNs) are proposed for complex applications such as image classification, regression, and scene recognition [15,16,29]. CNNs employ shift-invariant filters (kernels) followed by pooling units to extract local and global features from feature maps. By minimising the loss function with many hyperparameters, the network obtains the optimum weights and biases for the solving problem.
Chen et al. [9,30,31] first introduced neural networks to the mesh quality evaluation task. They proposed an automatic quality indicator for 2D NACA0012 airfoil meshes using CNNs. The indicator takes geometric characteristics of each mesh element as input (the edge length x, edge length y, and maximum included angle), then feeds them into the construed neural network to identify poor-quality NACA0012 meshes. However, due to the geometric properties of the input features, the input constructing process is computationally expensive, and the proposed method is only applied to two-dimensional meshes.
In this paper, we propose a neural network-based mesh quality indicator, accompanied by a benchmark dataset for three-dimensional cylinder meshes. In the mesh pre-processing phase, the indicator first splits the cylinder mesh into mesh surfaces and extracts mesh points from each surface. The proposed neural network Mesh-Net directly takes mesh points as input without geometric calculation. During the training, Mesh-Net employs fully-convolutional and global average layers to learn the role of mesh geometry and distribution on the accuracy of CFD simulation. The well-designed architecture makes it attractive as an indicator of variable-sized three-dimensional cylinders. After suitable training, the indicator is able to predict the overall quality of the input mesh precisely. It can also be applied in the automatic mesh quality evaluation process without manual interactions.

3. A Neural Network-Based Mesh Quality Indicator for Three-Dimensional Cylinder Modelling

Neural network-based mesh training is the optimisation process by which the relation between the input mesh and quality prediction is established. This process usually requires a large number of labelled mesh samples to learn accurately. However, since annotating CFD meshes with simulation accuracy can be time-consuming and expensive, there has not yet emerged a public three-dimensional mesh dataset. To support our study and address the problem of available mesh datasets, we developed a cylinder mesh benchmark dataset for neural network-based mesh quality evaluation.

3.1. Three-Dimensional Cylinder Mesh Benchmark Dataset

In this section, we introduce the process of building the mesh benchmark dataset used for training. Each mesh sample generation can be divided into four steps: (1) modelling, (2) transforming, (3) simulation, and (4) annotation.
In the initial modelling step, the geometric model of the three-dimensional cylinder was constructed. Then, we generated meshes that varied in mesh size and deformed them to obtain cylinder meshes with different qualities. To this end, we developed an automatic three-dimensional cylinder mesh generator. The generator takes mesh files as input and transforms the input mesh using point reposition, curve translation, or mesh surface rotation. Figure 2 illustrates some of the deformed cases. We can see that a large degree of variance in geometric transformations can be achieved. Using this generator, we have collected a large dataset with 20,480 cylinder meshes that span different mesh sizes and contain a wide variety of quality properties. Notice that the obtained non-fixed size meshes increase the richness of the proposed dataset and make it useful for mesh training tasks involving multiscale cylinder models.
During the simulation step, we performed numerical simulations for each mesh sample on a classical problem. The problem models the steady laminar flow between rotating and stationary concentric cylinders (see material properties in Table 1) [32]. Considering the inner cylinder has radius r 0 , angular velocity w 0 , and temperature T 0 , while the outer is r 1 and T 1 , we calculate the tangential velocity in the annulus at certain radial locations. The motion equations include velocity component u θ , radius r, T, and p as:
u θ = 0
d p d r = ρ u θ 2 r
d 2 u θ d r 2 + d d r u θ r = 0
k d d d r r d T d r + μ d u θ d r u θ r 2 = 0
with boundary conditions:
r = r 0 : u θ = r 0 ω 0 T = T 0 p = p 0
r = r 1 : u θ = 0 T = T 1
After the simulation, we compared the numerical solution with target results in [32] at four radial locations (20 mm, 25 mm, 30 mm, and 35 mm), and manually divided each sample into one of the following four quality categories:
(1)
High-quality Mesh: is a class of acceptable meshes with a very small error in the numerical solution.
(2)
Non-orthogonal Mesh: occurs when the curves or surfaces of the mesh are not vertically orthogonal. Numerical experiments in [5] show that skewed mesh with poor orthogonality can affect the order of accuracy and error magnitude. Non-orthogonal meshes also have a negative impact on the convergence speed.
(3)
Non-smoothness Mesh: is a class of meshes in which the length ratio is distorted, or elements are overlapped in complex domains. One approach to increase the quality is to smooth a collection of nodes (while preserving mesh connectivity) or to optimize node positions (vertex repositioning) [7,33].
(4)
Poor-quality Mesh: represents meshes with poor orthogonality, smoothness, and distribution. According to the analysis in [34], poorly-shaped meshes can cause the ill-conditioned stiffness matrix problem and seriously affect the solutions of the partial differential equations.
To verify the validity of the annotation procedure, we compared the numerical error of meshes in four quality categories. Figure 3 shows the numerical error of 20,480 meshes with different quality categories in the proposed benchmark dataset. We can learn that all high-quality meshes accurately simulate the fluid flow in the cylinder. For Non-orthogonal Mesh, there are small numerical errors (from −5% to 3%) during the CFD simulation. However, these meshes suffer from a slow convergence speed compared with meshes in the High-quality Mesh category. The numerical error of Non-smoothness Mesh ranges from 4% to 10%, while the poor-quality meshes leave a larger simulation error (up to 24.2%) compared to the target results.
Overall, we seek to construct a large collection of mesh samples with accurate solution-based labels. Such data is useful for supervised learning and neural network-based mesh quality evaluation. To achieve this, we built a three-dimensional cylinder mesh dataset containing a total of 20,480 meshes belonging to four categories, with an average of 512 meshes per size per category. The name and detail description of sizes and quality categories are listed in Table 2 and Table 3. Figure 4 shows some mesh samples in the proposed benchmark dataset. The diversity of the meshing ensures the richness and validity of the proposed dataset. We believe that this benchmark dataset contributes to developing advanced mesh understanding algorithms. It can also stimulate innovative research for CFD mesh quality evaluation tasks.

3.2. Mesh Pre-Processing

For the CNN-based mesh quality evaluation task, developing a representation scheme applicable to mesh samples is a prerequisite for neural network training. Due to the locally dense nature of CFD meshes, existing three-dimensional quantisation methods (e.g., multi-view or volumetric) do not apply to mesh samples. Point cloud features are able to handle the locally dense areas in underlying meshes. However, traditional point cloud representation ignores the spatial correlation between neighbouring points, which is crucially important for CFD mesh quality evaluation.
In our work, we introduce a point-based pre-processing method for cylinder mesh representation. In order to encode the spatial information of mesh points, we first split the cylinder into two-dimensional surfaces along the rotation axis, and then sequentially extract the mesh points from each mesh surface. After that, we combine the obtained point coordinates to form the three-channel point information matrix. Each channel of the matrix represents one of the dimensional coordinates (x, y, or z). The detail of the pre-processing method is shown in Figure 5.
Since the coordinates of each point are explicitly stored in the mesh source file, we can directly use the source file as training input. Compared to the mesh pre-processing in [9], which represents meshes using specific element features (edge length and included angle), our point-based representation is more efficient. It does not require any additional computation for three-dimensional cylinder meshes, which significantly reduces the pre-processing cost. Moreover, benefiting from the fact that the point information matrix incorporates spatial information, we can easily process input meshes without paying attention to the mesh size.
Normalisation is essentially a linear transformation that proportionally compresses and transforms a vector. This transformation keeps linear combinations and linear relational formulas intact, thus ensuring the robustness of a particular model. After normalising the input data, searching the optimal mapping in CNNs can be smoother (more likely to converge towards the optimal solution). Before training, we apply the standard deviation normalisation to the point information matrix. The normalisation formula is:
x * = x x ¯ σ x
y * = y y ¯ σ y
z * = z z ¯ σ z
where σ is the standard deviation of the original data, while x ¯ , y ¯ , z ¯ are the mean values of the original data, respectively.
Finally, the normalised feature matrix is fed into the proposed network Mesh-Net. The training process stops after converging to a local optimum. Thereafter, the trained network can be used as a black-box to analyse the meshing properties (smoothness, orthogonality, and distribution) from the cylinder point feature and automatically output the quality of the input mesh.

3.3. The Structure of Mesh-Net

We now describe the design of the proposed neural network Mesh-Net. It consists of an input layer, five convolutional layers, and a softmax layer. To keep the training and prediction cost low, we did not consider very deep architectures. The network architecture employs fully convolutional layers with no fully connected layer, which enables the network to take input meshes of arbitrary size and produce fixed size output. Figure 5 shows the architecture of Mesh-Net. As depicted in Figure 5, the number of channels (feature maps) in five convolutional layers is 16, 32, 64, 32, and 4, respectively.
As for the kernel size, we dynamically adjust the kernel size in each layer to obtain different receptive fields, rather than using the fixed kernels. At the beginning of the training, we prefer a relatively large receptive field to obtain more local point information. Inspired by the element dependency in the seven-point difference scheme, we set a 7 × 7 kernel in the first convolutional layer to capture the quality features of 49 adjacent mesh points. In the following layers, we gradually reduce the size of the convolutional kernel to obtain a smaller receptive field in high-level features. The kernel size of the next three layers is 5 × 5, 3 × 3, and 1 × 1, respectively.
There is no max-pooling layer in the proposed architecture. Instead, we set the stride in the first three convolutional layers to two to shrink the dimension of feature maps. It is worth noting that we employ a global average operation to calculate the mean value of the elements across dimensions in the fourth convolutional layers. After global averaging, the compressed feature maps are propagated to the softmax output function.
A loss function is employed during the network training phase to measure the discrepancy between the predicted output and the ground-truth tensor. In this work, we use cross-entropy cost function  L 0 to measure the discrepancy between two probability tensors. L 0 is closely related to the Kullback–Leibler divergence, as given by:
L 0 = 1 n [ y ln y ^ + ( 1 y ) ln ( 1 y ^ ) ]
where y ^ represents the approximation of ground-truth y, and n is the number of samples in the mini-batch. Since each part of the network is differentiable, we can compute the derivatives of L 0 and update the parameters with respect to the input. The optimisation can then proceed via backpropagation (gradient descent). The weights of the network are updated iteratively by:
ω i + 1 = ω i η L ω i
where ω i is the weight in the i-th forward propagation and η > 0 is the learning rate. This process culminates with a vector-valued output that values in [0, 1]. It can be viewed as the neural network approximation of the desired function or the probability that the input mesh falls into one of four quality categories: High-quality Mesh, Non-orthogonal Mesh, Non-smoothness Mesh, and Poor-quality Mesh.

4. Experimental Results and Discussions

4.1. Training

As with any neural network, the choice of hyperparameters can strongly affect the prediction performance and the rate of convergence. In Section 3.3, we have determined hyperparameters, including the number of layers, the number of layer channels, and the kernel size. In addition, we need to define some training-related parameters, such as activation function in each layer, batch size, and learning rate.
The activation function used in Mesh-Net is responsible for introducing non-linearity into the network. Since the parameter update in each iteration involves the gradient of the activation function, the obtained tiny gradient can lead to a slow convergence or trapping in the local optimum [29]. To accelerate the convergence and avoid vanishing gradient issues, we equipped the convolutional layers with a composite function of the R e L U activation function and batch normalisation. Moreover, we use mini-batches to take a single training step and reshuffle the training set in each epoch after exhausting the entire training set. We find that the stochasticity introduced by shuffling improves the stability and performance in test cases.
The overfitting phenomenon is another problem existing in the training phase (i.e., the trained network ties too closely to the training set and behaves badly during testing). To tackle this problem, we combine the loss function with a regularisation term to avoid overfitting.
L = L 0 + λ 2 n | | ω | | 2
where L 0 is the loss function in Equation (15), ω represents the weights in Mesh-Net, n is the batch size, and λ is the regularisation coefficient.
We use the A d a m optimiser [35] with an exponential learning rate decay, which prevents the training from trapping in a local minimum. A comparison of different settings (batch size and learning rate) is shown in Figure 6. To make the best use of Mesh-Net, we set the initial learning rate to 0.0005, while the batch size was set to 32. The training was performed using the open-source machine learning library TensorFlow [36].

4.2. Prediction

To demonstrate the capability of Mesh-Net when used as a quality indicator, we compared the predictive power of different classifiers on three-dimensional benchmark datasets. During the experiment, we randomly shuffled the samples and employed the first 75% of meshes in each size for training and the latter 25% for testing. For each set, the proportion of different categories of meshes is equal. We ran the training 10 times and took the average accuracy as the final prediction result.
We compared the performance of the Mesh-Net with that of three widely-used machine learning algorithms and one multi-layer perceptron (MLP). The machine learning algorithms are support vector machine (SVM), quadratic discriminant analysis (QDA), and Gaussian Naive Bayes (GNB) [36]. The MLP used in this paper contains five layers, i.e., an input layer, three hidden layers, and an output layer. The number of neural units in three hidden layers is 16, 32, and 64, respectively. The number of neural units in the output layer is 4, which equals the number of quality categories.
Many metrics can be used to measure the performance of neural networks, such as r e c a l l , a c c u r a c y , and F 1 - s c o r e . Since the number of samples in each category balanced (512 × 10 per category), we only use a c c u r a c y to evaluate the performance of different classifiers. The accuracy a c c is defined as:
a c c = T P + T N P + N
where true positive ( T P ) is the number of correctly classified positive instances, true negative ( T N ) is the number of correctly classified negative instances, and P + N represents the total number of instances.
Table 4 reports the accuracy of different methods on the three-dimensional cylinder mesh dataset. Constrained by the classifiers’ limitations on the dimensionality of the input samples, we divided the experiments into two parts. The first part is the mesh training on fixed-size samples (Size 1 test). In this part, all five classifiers are trained and tested. In the second part, Mesh-Net, which accepts meshes of arbitrary size, is trained to test the overall performance across different mesh sizes (Full-size test). As can be seen in Table 4, the accuracy rate of machine learning algorithms is relatively low on the fixed-size test. All machine learning classifiers, SVM, QDA, and GNB, show a prediction accuracy of less than 90%. MLP achieves an accuracy of 95.70% on the mesh quality evaluation task. However, it is clear that the proposed CNN-based indicator is more effective than widely-used machine learning algorithms and MLP. It outperforms other trained classifiers and achieves an accuracy of 98.05% on fixed-size meshes and 96.60% on non-fixed size meshes.
To better understand the prediction across different categories, we present the confusion matrix of the full-size test (see Figure 7). We found that meshes with high-quality (HQ-M) and non-orthogonality (NO-M) achieve accurate quality prediction. Only four meshes (0.31%) were wrongly predicted. For meshes with non-smoothness (NS-M), 45 (3.52%) meshes were misclassified to poor-quality mesh (PQ-M). The category of PQ-M was predicted with the lowest accuracy. The results show that 121 (9.45%) poor-quality testing meshes were wrongly classified. Thirty-four (2.66%) PQ-M samples were misclassified to NO-M, and 87 (6.8%) PQ-M samples to NS-M. The inaccuracy is mainly because mesh quality properties such as non-orthogonality and non-smoothness are easily confused, especially when the point reposition or surface rotation happens. However, the incorrect predictions from Mesh-Net still make sense to the meshing procedure. It identifies part of the quality defects in the input mesh and guides the subsequent mesh optimisation.
We note that the computational complexity increases considerably with the number of neural units. The structure of the CNN-based indicator must be intelligent enough to make the classification task possible and simple enough to keep the training and prediction cost low. Thus, we did not consider very deep architectures. Moreover, the proposed network employs full convolutional layers without fully connected layers, which greatly reduces the number of parameters (orders of magnitude) compared with MLP. The introduction of fully convolutional layers also allows the input of meshes with different sizes.
Overall, we propose a CNN-based quality indicator for three-dimensional cylinder meshes. The proposed network Mesh-Net fully exploits the advantages of receptive field properties of convolutional neural networks. It employs different sizes of kernels to capture the local and global quality features of the preprocessed mesh. During the training, the network learns the relationship between the quality of the cylinder mesh and the error convergence of CFD simulation. Thereafter, the trained indicator can be used as an intelligent quality control model to evaluate mesh quality before CFD simulations.

5. Conclusions

Mesh-based methods have proved extremely useful in computational fluid dynamics simulations. During the cylinder simulation, the quality of the preassigned mesh affects the accuracy of numerical solutions. Poorly shaped meshes tend to slow convergence or cause analysing instability. Many quality indicators have been proposed to serve as quality control by analysing the geometric information of mesh elements. However, these element-based indicators do not necessarily provide reliable guidance for the subsequent optimisation process. They also require frequent human-computer interactions during the evaluation, which significantly increases the meshing overhead. Therefore, it is desirable to develop an intelligent indicator that automatically learns the quality of the mesh.
In this paper, we present an efficient mesh quality indicator by using convolution neural networks (CNNs). To support our study, we also release a three-dimensional cylinder mesh dataset, which contains 20,480 meshes, with different sizes and qualities. The proposed indicator is trained offline and employs a feedforward approximation to learn the mesh quality properties, such as orthogonality, smoothness, and mesh distribution. It takes mesh files as input and outputs the overall quality of the input mesh to determine if it meets the solver’s requirements. Experimental results show that the proposed method is accurate, computationally efficient, and straightforward.
We believe that the applications of deep learning methods to mesh quality problems are expected to address the challenges posed by frequent manual interactions and reduce the meshing cost. We also hope that the release of large-scale datasets can stimulate innovative research on mesh quality evaluation and advance the development of fully automatic mesh generation.

Author Contributions

Conceptualisation, X.C.; methodology, X.C.; software, Z.W.; writing—review & editing, C.G.; writing—original draft preparation, X.C. and C.G.; supervision, J.L. and Y.P.; funding acquisition, J.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research work was supported in part by the National Key Research and Development Program of China (2021YFB0300101).

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Chen, X.; Liu, J.; Li, S.; Xie, P.; Chi, L.; Wang, Q. TAMM: A New Topology-Aware Mapping Method for Parallel Applications on the Tianhe-2A Supercomputer. In Algorithms and Architectures for Parallel Processing; Vaidya, J., Li, J., Eds.; Springer International Publishing: Cham, Switzerland, 2018; pp. 242–256. [Google Scholar]
  2. Garimella, R.V.; Shashkov, M.; Knupp, P.M. Triangular and quadrilateral surface mesh quality optimization using local parametrization. Comput. Methods Appl. Mech. Eng. 2004, 193, 913–928. [Google Scholar] [CrossRef]
  3. Shi, R.; Lin, J.; Yang, H. Numerical Study on the Coagulation and Breakage of Nanoparticles in the Two-Phase Flow around Cylinders. Entropy 2022, 24, 526. [Google Scholar] [CrossRef] [PubMed]
  4. Chen, X.; Li, T.; Wan, Q.; He, X.; Gong, C.; Pang, Y.; Liu, J. MGNet: A novel differential mesh generation method based on unsupervised neural networks. Eng. Comput. 2022, 1–13. [Google Scholar] [CrossRef]
  5. Blacker, T.D.; Owen, S.J.; Staten, M.L.; Quadros, W.R.; Hanks, B.; Clark, B.W.; Meyers, R.J.; Ernst, C.; Merkley, K.; Morris, R.; et al. CUBIT Geometry and Mesh Generation Toolkit 15.1 User Documentation; U.S. Department of Energy Office of Scientific and Technical Information: Washington, DC, USA, 2016. [CrossRef]
  6. Berzins, M. Mesh Quality: A Function of Geometry, Error Estimates or Both? Eng. Comput. 1999, 15, 236–247. [Google Scholar] [CrossRef]
  7. Katz, A.; Sankaran, V. Mesh quality effects on the accuracy of CFD solutions on unstructured meshes. J. Comput. Phys. 2011, 230, 7670–7686. [Google Scholar] [CrossRef]
  8. Lowrie, W.; Lukin, V.S.; Shumlak, U. A Priori Mesh Quality Metric Error Analysis Applied to a High-Order Finite Element Method. J. Comput. Phys. 2011, 230, 5564–5586. [Google Scholar] [CrossRef]
  9. Chen, X.; Liu, J.; Pang, Y.; Chen, J.; Chi, L.; Gong, C. Developing a new mesh quality evaluation method based on convolutional neural network. Eng. Appl. Comput. Fluid Mech. 2020, 14, 391–400. [Google Scholar] [CrossRef]
  10. Diskin, B.; Thomas, J. Effects of mesh regularity on accuracy of finite-volume schemes. In Proceedings of the 50th AIAA Aerospace Sciences Meeting Including the New Horizons Forum and Aerospace Exposition, Nashville, TN, USA, 9–12 January 2012. [Google Scholar] [CrossRef]
  11. Lo, S. Optimization of tetrahedral meshes based on element shape measures. Comput. Struct. 1997, 63, 951–961. [Google Scholar] [CrossRef]
  12. Huang, Y.; Liupeng, W.; Yi, N. Mesh Quality and More Detailed Error Estimates of Finite Element Method. Numer. Math. Theory, Methods Appl. 2017, 10, 420–436. [Google Scholar] [CrossRef]
  13. Nie, C.; Liu, J.; Sun, S. Study on quality measures fpr tetrahedral mesh. Chin. J. Comput. Mech. 2003, 20, 579–582. [Google Scholar]
  14. Gao, X.; Huang, J.; Xu, K.; Pan, Z.; Deng, Z.; Chen, G. Evaluating Hex-mesh Quality Metrics via Correlation Analysis. Comput. Graph. Forum 2017, 36, 105–116. [Google Scholar] [CrossRef]
  15. Hu, J.; Shen, L.; Sun, G. Squeeze-and-excitation networks. In Proceedings of the IEEE Conference on Computer Vision and Pattern Recognition, Salt Lake City, UT, USA, 18–23 June 2018; pp. 7132–7141. [Google Scholar]
  16. Sandler, M.; Howard, A.; Zhu, M.; Zhmoginov, A.; Chen, L.C. Mobilenetv2: Inverted residuals and linear bottlenecks. In Proceedings of the IEEE Conference on Computer Vision and Pattern Recognition, Salt Lake City, UT, USA, 18–23 June 2018; pp. 4510–4520. [Google Scholar]
  17. Li, M.; Li, S.; Zhang, J.; Wu, F.; Zhang, T. Neural Adaptive Funnel Dynamic Surface Control with Disturbance-Observer for the PMSM with Time Delays. Entropy 2022, 24, 1028. [Google Scholar] [CrossRef] [PubMed]
  18. Ray, D.; Hesthaven, J.S. An artificial neural network as a troubled-cell indicator. J. Comput. Phys. 2018, 367, 166–191. [Google Scholar] [CrossRef]
  19. Knupp, P.M. Achieving finite element mesh quality via optimization of the Jacobian matrix norm and associated quantities. Part I—A framework for surface mesh optimization. Int. J. Numer. Methods Eng. 2000, 48, 401–420. [Google Scholar] [CrossRef]
  20. Geuzaine, C.; Remacle, J.F. Gmsh: A three-dimensional finite element mesh generator with built-in pre- and post-processing facilities. Int. J. Numer. Methods Eng. 2009, 11, 1309–1331. [Google Scholar] [CrossRef]
  21. Strang, G.; Fix, G. An Analysis of The Finite Element Method. J. Appl. Mech. 1973, 41, 62. [Google Scholar] [CrossRef]
  22. Shewchuk, J.R. Delaunay refinement algorithms for triangular mesh generation. Comput. Geom. 2002, 22, 21–74. [Google Scholar] [CrossRef]
  23. Liu, A.; Joe, B. Relationship between tetrahedron shape measures. BIT Numer. Math. 1994, 34, 268–287. [Google Scholar] [CrossRef]
  24. Bank, R.E.; Smith, R.K. Mesh Smoothing Using A Posteriori Error Estimates. SIAM J. Numer. Anal. 1997, 34, 979–997. [Google Scholar] [CrossRef] [Green Version]
  25. Weatherill, N.P.; Morgan, K.; Hassan, O.; Jones, J.W. Large-scale aerospace simulations using unstructured grid technology. In Computational Mechanics for the Twenty-First Century; Civil-Comp Press: Edinburgh, UK, 2000. [Google Scholar]
  26. Antiga, L.; Ene-Iordache, B.; Caverni, L.; Cornalba, G.P.; Remuzzi, A. Geometric reconstruction for computational mesh generation of arterial bifurcations from CT angiography. Comput. Med. Imaging Graph. 2002, 26, 227–235. [Google Scholar] [CrossRef]
  27. Chong, C.S.; Lee, H.P.; Kumar, A.S. Genetic algorithms in mesh optimization for visualization and finite element models. Neural Comput. Appl. 2006, 15, 366–372. [Google Scholar] [CrossRef]
  28. Mall, S.; Chakraverty, S. Comparison of Artificial Neural Network Architecture in Solving Ordinary Differential Equations. Adv. Artif. Neural Syst. 2013, 2013, 181895. [Google Scholar] [CrossRef]
  29. Tan, M.; Le, Q.V. Efficientnet: Rethinking model scaling for convolutional neural networks. In Proceedings of the International Conference on Machine Learning, Long Beach, CA, USA, 9–15 June 2019. [Google Scholar]
  30. Chen, X.; Liu, J.; Gong, C.; Pang, Y.; Chen, B. An Airfoil Mesh Quality Criterion using Deep Neural Networks. In Proceedings of the 12th International Conference on Advanced Computational Intelligence, Dali, China, 14–16 August 2020; pp. 536–541. [Google Scholar]
  31. Chen, X.; Gong, C.; Liu, J.; Pang, Y.; Liang, D.; Chi, L.; Li, K. A Novel Neural Network Approach for Airfoil Mesh Quality Evaluation. J. Parallel Distrib. Comput. 2022, 164, 123–132. [Google Scholar] [CrossRef]
  32. Farwig, R. An Lq-analysis of viscous fluid flow past a rotating obstacle. Tohoku Math. J. 2006, 58, 129–147. [Google Scholar] [CrossRef]
  33. Knupp, P.M. Matrix Norms and the Condition Number: A General Framework to Improve Mesh Quality via Node-Movement. Off. Sci. Tech. Inf. Tech. Rep. 1999, 36, 13–22. [Google Scholar]
  34. Zhang, Y.; Jia, Y.; Wang, S.S.Y. Structured mesh generation with smoothness controls. Int. J. Numer. Methods Fluids 2006, 51, 1255–1276. [Google Scholar] [CrossRef]
  35. Qi, C.R.; Su, H.; Mo, K.; Guibas, L.J. Pointnet: Deep learning on point sets for 3d classification and segmentation. In Proceedings of the IEEE Conference on Computer Vision and Pattern Recognition, Honolulu, HI, USA, 21–26 July 2017; pp. 652–660. [Google Scholar]
  36. Abadi, M.; Barham, P.; Chen, J.; Chen, Z.; Davis, A.; Dean, J.; Devin, M.; Ghemawat, S.; Irving, G.; Isard, M.; et al. TensorFlow: A System for Large-Scale Machine Learning. In Proceedings of the 12th USENIX Conference on Operating Systems Design and Implementation, Savannah, GA, USA, 2–4 November 2016; USENIX Association: Berkeley, CA, USA, 2016; pp. 265–283. [Google Scholar]
Figure 1. The architecture of MLP.
Figure 1. The architecture of MLP.
Entropy 24 01245 g001
Figure 2. Mesh transforming used in the automatic three-dimensional cylinder mesh generator. (a) Initial mesh, (b) curve translation, (c) mesh surface rotation, (d) point reposition.
Figure 2. Mesh transforming used in the automatic three-dimensional cylinder mesh generator. (a) Initial mesh, (b) curve translation, (c) mesh surface rotation, (d) point reposition.
Entropy 24 01245 g002
Figure 3. The numerical error of different quality meshes in the proposed benchmark dataset.
Figure 3. The numerical error of different quality meshes in the proposed benchmark dataset.
Entropy 24 01245 g003
Figure 4. Mesh samples in the proposed mesh benchmark dataset. The dataset consists of 20,480 labelled cylinder meshes with different sizes and qualities. Each sample falls into one of four categories: High-Quality Mesh (HQ-M), Non-orthogonal Mesh (NO-M), Non-smoothness Mesh (NS-M), or Poor-quality Mesh (PQ-M).
Figure 4. Mesh samples in the proposed mesh benchmark dataset. The dataset consists of 20,480 labelled cylinder meshes with different sizes and qualities. Each sample falls into one of four categories: High-Quality Mesh (HQ-M), Non-orthogonal Mesh (NO-M), Non-smoothness Mesh (NS-M), or Poor-quality Mesh (PQ-M).
Entropy 24 01245 g004
Figure 5. The proposed mesh pre-processing method and the architecture of Mesh-Net.
Figure 5. The proposed mesh pre-processing method and the architecture of Mesh-Net.
Entropy 24 01245 g005
Figure 6. The performance of different batch sizes (left) and learning rates (right).
Figure 6. The performance of different batch sizes (left) and learning rates (right).
Entropy 24 01245 g006
Figure 7. The confusion matrix of the full-size test.
Figure 7. The confusion matrix of the full-size test.
Entropy 24 01245 g007
Table 1. The material properties of the three-dimensional cylinder.
Table 1. The material properties of the three-dimensional cylinder.
ParameterValue
Density1 kg/m 3
Viscosity0.0002 kg/m 5
Radius of the inner cylinder17.8 mm
Radius of the outer cylinder46.28 mm
Angular velocity of the inner wall1 rad/s
Table 2. Ten different sizes of samples in the proposed mesh dataset.
Table 2. Ten different sizes of samples in the proposed mesh dataset.
CaseMesh SizeNumber of Meshes
Size 196 × 30 × 202048
Size 2100 × 30 × 202048
Size 3104 × 30 × 202048
Size 4108 × 30 × 202048
Size 5112 × 31 × 202048
Size 6116 × 31 × 202048
Size 7120 × 31 × 202048
Size 8124 × 31 × 202048
Size 9128 × 32 × 212048
Size 10132 × 32 × 212048
Table 3. The name and detail description of four quality categories.
Table 3. The name and detail description of four quality categories.
LabelQuality CategoriesNumber of Meshes
1 (HQ-M)High-quality Mesh512 × 10
2 (NO-M)Non-orthogonal Mesh512 × 10
3 (NS-M)Non-smoothness Mesh512 × 10
4 (PQ-M)Poor-quality Mesh512 × 10
Table 4. The performance of different classifiers on the cylinder mesh dataset.
Table 4. The performance of different classifiers on the cylinder mesh dataset.
CaseModelAccuracy (%)
Size 1 testSVM89.06%
QDA87.30%
GNB79.49%
MLP95.70%
Mesh-Net98.05%
Full-size testMesh-Net96.60%
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Chen, X.; Wang, Z.; Liu, J.; Gong, C.; Pang, Y. A Neural Network-Based Mesh Quality Indicator for Three-Dimensional Cylinder Modelling. Entropy 2022, 24, 1245. https://doi.org/10.3390/e24091245

AMA Style

Chen X, Wang Z, Liu J, Gong C, Pang Y. A Neural Network-Based Mesh Quality Indicator for Three-Dimensional Cylinder Modelling. Entropy. 2022; 24(9):1245. https://doi.org/10.3390/e24091245

Chicago/Turabian Style

Chen, Xinhai, Zhichao Wang, Jie Liu, Chunye Gong, and Yufei Pang. 2022. "A Neural Network-Based Mesh Quality Indicator for Three-Dimensional Cylinder Modelling" Entropy 24, no. 9: 1245. https://doi.org/10.3390/e24091245

APA Style

Chen, X., Wang, Z., Liu, J., Gong, C., & Pang, Y. (2022). A Neural Network-Based Mesh Quality Indicator for Three-Dimensional Cylinder Modelling. Entropy, 24(9), 1245. https://doi.org/10.3390/e24091245

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop