COMSOL Multiphysics™

meshrefine
clear fem
fem.geom = square2 + move(square2,0,-1) + move(square2,-1,-1);
fem.mesh = meshinit(fem,'hmax',0.4);
subplot(2,2,1), meshplot(fem)
fem.mesh = meshrefine(fem,'tri',[1:50;ones(1,50)]);
subplot(2,2,2), meshplot(fem)
fem.mesh = meshrefine(fem,'subdomain',...
[1 2],'rmethod','longest');
subplot(2,2,3), meshplot(fem)
fem.mesh = meshrefine(fem,'subdomain',1);
subplot(2,2,4), meshplot(fem)
Algorithm
The 2D algorithm is described by the steps below:
1 Pick the initial set of elements to be refined.
2 Either divide all edges of the selected elements in half (regular refinement), or
divide the longest edge in half (longest edge refinement).
3 Divide the longest edge of any element that has a divided edge.
4 Repeat step 3 until no further edges are divided.
5 Introduce new points on all divided edges, and replace all divided entries in e by
two new entries.
6 Form the new elements. If all three sides are divided, new elements are formed
by joining the side midpoints. If two sides are divided, the midpoint of the
longest edge is joined with the opposing corner and with the other midpoint. If
only the longest edge is divided, its midpoint is joined with the opposing corner.
Compatibility
To achieve compatibility with FEMLAB 2.3, the geometry input is automatically
converted to a geometry object using the function geomobject. The geometry
input can be any analyzed geometry. See geomobjectfor details.
The output u is obsolete from FEMLAB 2.2. Use asseminit to interpolate the
solution to a new mesh.
See Also
femmesh, geomcsg, meshinit
292 | CHAPTER 1: COMMAND REFERENCE