6
$\begingroup$

While answering one of the questions on meshing process, I encountered a lack of understanding on my end for the comparison of the mesh quality.

First, consider an unstructured mesh created in GMSH for a 2-D shape (I am not interested in this particular shape, just here for the illustration). The elements are controlled only through the characteristic length, which basically dictates, how many elements are created along the geometry boundaries.

Fig. 1 Unstructured mesh

Now, I can detect the regions in the mesh that are "kind of suitable" for a structured mesh. That could be done automatically, say, using marching cubes.

Fig. 1 Hybrid structured/unstructured mesh mesh

Now, both meshes satisfy the original criteria for the size of the element. But the second mesh has a lot of elements that are exactly the same - allowing for significant savings in the required amount of computations while featuring a similar number of triangular elements (even a little bit less).

Is the second mesh bad in general? Does a lot depend on the kind of basis functions that will be used?

Now, I am particularly interested in using the meshes in frequency-domain computational electromagnetics, where the common meshing is done at $\lambda/10$, $\lambda$ is the wavelength. Say, both meshes satisfy $\lambda/10$ requirement (largest triangular edge is smaller than $\lambda/10$). Is there something fundamentally bad by using this type of a hybrid structured/unstructured mesh?

NB: I am usually using integral-equation methods; however, don't mind having an understanding of what happens, in general, and for FEM and others, in particular.

$\endgroup$
2
  • $\begingroup$ I don't think the second is bad per se (except that it might polarise your solution a bit since there is a distinct anisotropy in your mesh orientation). That's a very interesting question though... $\endgroup$
    – BlaB
    May 15, 2018 at 22:49
  • $\begingroup$ I'd actually also be curious if anyone had some ideas regarding mixing quads and tets in this type of setup as well. $\endgroup$ May 16, 2018 at 0:19

2 Answers 2

8
$\begingroup$

In my opinion, it is not a good neither a bad mesh. It clearly depends on the PDE you are considering.

The finite space to which the PDE is projected is your mesh, where your operators, e.g. $\vec{\textrm{grad}}$ (gradients), $\textrm{div}$ (divergences), $\triangle$ (Laplacians)... strongly depend on that mesh and become matrices: $$ \vec{\textrm{grad}}\rightarrow\vec{\textrm{grad}}_h=G\quad \textrm{div}\rightarrow\textrm{div}_h=D \quad \triangle\rightarrow\triangle_h=L$$

Elliptic problems do not suffer from problems associated with the mesh, in terms of isotropy/anisotropy because the nature of the associated operator (Laplacian) is isotropic (invariant under translations and rotations), therefore your mesh would be suitable for elliptic problems.

However, other operators like the $\vec{\textrm{grad}}$ depend strongly on the mesh, because they are anisotropic (depend on the direction). This fact is important if you are going to project a hyperbolic equation (Fluid dynamics) onto your mesh, where a prevailing direction arises and thus the solution is affected.

As an example of these matters try to compute the following problems with these two meshes:

  • Mesh a) Square orthogonal mesh.
  • Mesh b) Split mesh a) into triangles.

Problem 1) Solve, for example, the elliptic problem $-\triangle u=f\; \text{in}\;\Omega$ with $u=0\,\forall x\in\partial \Omega$ and a smooth $f$ in the two meshes. Use for example Q1 finite elements for the mesh a) and P2 FE in your mesh b). Compare them.

Problem 2) Solve, for example, the hyperbolic problem $\partial_t u-(\vec{v}^{T}\vec{\textrm{grad}})u = 0$ with $u=g(t)$ and $u(0)=0$ in the inflow boundary, in the two meshes (with $\vec{v}$ in the orthogonal/parallel direction of the boundaries). Use for example finite volumes of second order. Compare them. Now solve both problems with $\vec{v}$ entering the domain with a certain angle. Compare them again.

In problem 1) you won't see great changes in your solution, while in problem 2, the orthogonal mesh is really suitable when solving for $\vec{v}$ oriented perpendicularly to the inflow boundary but misbehaves for an arbitrary angle (clearly the problem is not isotropic). For the triangular mesh, you shouldn't note such a great difference, because a certain anisotropy is already introduced by the mesh.

$\endgroup$
3
$\begingroup$

I think that both meshes are good. But depending of the problem at hand, one might be better than the other. In the domain that you are mentioning, one thing to consider is that you probably don't want element sizes to be too different to satisfy the sampling criterion.

Regarding mesh quality I would suggest that you check this post:

And, if you check reference 1, you would see that depending on the problem at hand some criteria would be better than others, anisotropy in your might be a good thing , for example. You would like to have a finer discretization in directions with higher derivatives (gradients).

Regarding the meshing of the geometry in your example. I would probably go with something like the following

enter image description here

where I have fine tuning over discretization levels for different regions. Thus, I can refine wherever I find it needed according to the solution for a coarser mesh.

References

  1. Jonathan Richard Shewchuk (2002). "What is a Good Linear Finite Element?".

  2. Pébay, Philippe, and Timothy Baker. "Analysis of triangle quality measures." Mathematics of Computation 72.244 (2003): 1817-1839.

$\endgroup$
4
  • $\begingroup$ That's actually a very pretty mesh! And thanks for the nice references, I haven't encountered [2] yet, for some reason. $\endgroup$
    – Anton Menshov
    May 16, 2018 at 14:22
  • 1
    $\begingroup$ @AntonMenshov, hehehe, thanks. Although I don't make them this way just for the sake of being pretty ;). Isn't the link in reference 2 working? $\endgroup$
    – nicoguaro
    May 16, 2018 at 15:23
  • $\begingroup$ No, it is definitely working. I meant I haven't seen the paper [2] before, while I am very familiar with [1]. $\endgroup$
    – Anton Menshov
    May 16, 2018 at 15:27
  • 1
    $\begingroup$ @AntonMenshov, ohhh, you used encountered in a completely different way from my understanding. I think I encountered that reference first in the Gmsh paper. $\endgroup$
    – nicoguaro
    May 16, 2018 at 15:40

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge that you have read and understand our privacy policy and code of conduct.

Not the answer you're looking for? Browse other questions tagged or ask your own question.