Scribd
Upload
150 views

Introduction To Mesh in Matlab

This document provides an introduction to mesh generation in Matlab using DistMesh. It explains that meshes are needed to discretize physical domains for numerical computations in methods like finite volume and finite element methods. DistMesh is introduced as a free mesh generation software that uses a distance function-based approach. The document outlines DistMesh's algorithm, shows examples of how to generate simple and refined meshes, and how to extract boundary information.

Uploaded by

lino
Copyright
© © All Rights Reserved
Available Formats
Download as PDF, TXT or read online on Scribd
150 views

Introduction To Mesh in Matlab

This document provides an introduction to mesh generation in Matlab using DistMesh. It explains that meshes are needed to discretize physical domains for numerical computations in methods like finite volume and finite element methods. DistMesh is introduced as a free mesh generation software that uses a distance function-based approach. The document outlines DistMesh's algorithm, shows examples of how to generate simple and refined meshes, and how to extract boundary information.

Uploaded by

lino
Copyright
© © All Rights Reserved
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 18

Introduction to mesh generation (in Matlab)

By Allan P. Engsig-Karup

Overview
!! Introduction to mesh generation !! Introduction to DistMesh for Matlab !! Goal: Introduce you to DistMesh for use with DG-FEM based models.

2/18

Why do we need a mesh?


!! We need to represent the (usually finite) physical domain in some way discretely for numerical computations. !! In sub domain methods, e.g. Finite volume or FEM methods, it is possible to independently consider the problem solution procedure and mesh generation as two distinct problems. !! This is very convenient if we want to solve more than one problem governed by the same PDEs!
3/18

What defines a mesh?


!! Here we define a mesh as a discrete representation !h of some spatial domain or topology !. !! A mesh can be sub divided into K smaller non-overlapping sub domains K such that !k k . h !h = ! !h !! Mesh generation can be a demanding and non-trivial task. E.g. for complex geometries or objects. !! Unstructured triangular meshes have good support for representing complex domains (or geometries) and mesh adaption (coarsening/refinement).
k =1
F-15. From: www.USEMe.org

Example. Triangulation.

4/18

What defines a mesh?


!! A mesh can be completely defined in terms of (unique) vertices and a mesh element table (triangulation). !! For the purpose of specifying appropriate boundary conditions we may for convenience use a boundary type table. !! Simple meshes can be created manually by hand. However, automatic mesh generation is generally faster and more efficient, although may require some user input for handling complex meshes. !! Note: Mesh data can conveniently be stored for reuse several times.
5/18

Mesh generators available?


!! Lots of standard mesh generators available! These generators can be used to solve a given mesh generation problem. (but may require a translation script) !! An example of a free software distribution for generating unstructured and triangular meshes is DistMesh (Matlab).

6/18

Introduction to DistMesh for Matlab


!! Persson, P.-O. and Strang, G. 2004 A simple mesh generator in Matlab. SIAM Review. Download scripts at: http://www-math.mit.edu/~persson/mesh/index.html !! A simple algorithm that combines a physical principle of force equilibrium in a truss structure with a mathematical representation of the geometry using signed distance functions.

7/18

Introduction to DistMesh for Matlab


!! Algorithm (Conceptual);
Step 1. Define a domain using signed distance functions. Step 2. Distribute a set of nodes interior to the domain. Step 3. Move interior nodes to obtain force equilibirum. Step 4. Apply terminate criterion when all nodes are fixed in space.

!! Post-processing steps (Preparation);


Step 5. Validate output! Step 6. Reorder element vertices to be defined anti-clockwise for use with DG-FEM. Step 7. Setup boundary table. Step 8. Store mesh for reuse.
8/18

Introduction to DistMesh for Matlab


Signed distance function, d(x);

(interior) $< 0 , x ( % ! d ( x) = # 0 , x ( '% (boundary) !> 0 , x & % (exterior) "


Define metric using an appropriate norm. E.g. The usual Euclidian metric.

d=0 !" d<0

#!" d>0
9/18

Introduction to DistMesh for Matlab


!! Combine geometries defined by distance functions using the Union, difference and intersection operations (set theory);

Union: Difference: Intersection;

min(d1 ( x), d2 ( x))


max(d1 ( x ), !d2 ( x )) " max(!d1 ( x ), d2 ( x ))

max(d1 ( x), d2 ( x))


10/18

Introduction to DistMesh for Matlab


Example 1. Create a uniform mesh using DistMesh. Square with hole

Visualized distance function

Mesh

Using DistMesh (in Matlab) in only 3 lines of code:


>> fd=inline('ddiff(drectangle(p,-1,1,-1,1),dcircle(p,0,0,0.4))','p'); >> pfix = [-1,-1;-1,1;1,-1;1,1]; >> [p,t] = distmesh2d(fd,@huniform,0.125,[-1,-1;1,1],pfix);
11/18

Introduction to DistMesh for Matlab


DistMesh output; (two tables) p t Unique vertice coordinates Element to Vertice table (not reordered automatically by DistMesh)

From this we can determine, e.g. >> K=size(t,1); >> Nv=size(p,1); >> Nfaces=size(t,2); >> VX = p(:,1); >> VY = p(:,2); >> EToV = t; %Number of elements %Number of vertices in mesh %Number of faces/element %Vertice x-coordinates %Vertice y-coordinates %Element to Vertice table

12/18

Introduction to DistMesh for Matlab


DG-FEM convention for standard element definitions;

!! Vertices are numbered anti-clockwise. !! Faces are numbered anti-clockwise with the first face beeing the one that connects the first two vertices.

13/18

Introduction to DistMesh for Matlab


Example 2. Create a refined mesh using DistMesh.
Element size function

Size*h0

Square with hole

2.5

1.5

1 1 0.5 0 -0.5 y -1 -1 -0.5 x 0

1 0.5

Visualized element size function

Mesh

>> fd = inline('ddiff(drectangle(p,-1,1,-1,1),dcircle(p,0,0,0.4))','p'); >> pfix = [-1,-1;-1,1;1,-1;1,1]; >> fh = inline(['min( sqrt( p(:,1).^2 + p(:,2).^2 ) , 1 )'],'p'); >> [p,t] = distmesh2d(fd,fh,0.125/2.5,[-1,-1;1,1],pfix);
14/18

Introduction to DistMesh for Matlab


Element size function in DistMesh;
Element size function

2.5

Size*h0

1.5

1 1 0.5 0 -0.5 y -1 -1 -0.5 x 0

1 0.5

From former example;

Visualized element size function

Mesh

>> fh = inline(['min( sqrt( p(:,1).^2 + p(:,2).^2 ) , 1 )'],'p');

>> [p,t] = distmesh2d(fd,fh,h0,[-1,-1;1,1],pfix);


!! Function fh Defines relative sizes of elements in final mesh. (fh=constant result in uniform distribution) !! The initial characteristic size of the elements is h0. !! In final distribution, the characteristic size of the smallest elements in the mesh will be approx. h0;

15/18

Introduction to DistMesh for Matlab


Example 3. Selecting boundary nodes. Square with hole
Nodes can be selected using distance functions; |d| = 0 or |d| < tol

Inner boundary nodes

Outer boundary nodes

>> fdInner = inline(dcircle(p,0,0,0.4),p); >> nodesInner = find(abs(fdInner([p]))<1e-3); >> fdOuter = inline(drectangle(p,-1,1,-1,1),p); >> nodesOuter = find(abs(fdOuter([p]))<1e-3); >> nodesB = find(abs(fd([p]))<1e-3);
16/18

Introduction to DistMesh for Matlab


Example 4. Updating boundary table. Square with hole
>> BCcode = 99; >> BCType = zeros(size(EToV')); % empty BCType table >> BCType = CorrectBCTable(K,EToV,BCType,nodesB,BCcode); The BCType boundary table can be used to create different maps (see the script BuildBCMaps2D.m, Section 6.4 in the textbook) for imposing different types of boundary conditions.
17/18

Final remarks
These notes together with example scripts for DistMesh can be found at my webpage: http://www.imm.dtu.dk/~apek/

Feel free to contact me at [email protected]


18/18

You might also like