How can I refine a subdomain in the PDE Toolbox mesh generation tool?

6 次查看（过去 30 天）
Zoi T.
Zoi T. 2016-9-7
编辑： Sean Littleton 2021-1-15
Hello,
I would like to refine a part of my domain without refining the whole mesh. Say for example that I have a unit square domain R1, which contains a circular sub-domain E1. Can I refine only E1 up to a desired element size, without imposing the same mesh size in the rest of the domain?

请先登录，再回答此问题。

采纳的回答

michio
michio 2016-9-7
I assume your problem is in 2D, then refinemesh function accepts an extra input argument "it", that is interpreted as a list of subdomains to refine, if it is a row vector.
See. http://jp.mathworks.com/help/pde/ug/refinemesh.html
  6 个评论
显示 4更早的评论
MarionJ
MarionJ 2017-6-30
It is not possible. I had the same problem. You have to use assempde, parabolic etc. if you want to refine mesh in certain Domains.
Ali Hammouche
Ali Hammouche 2020-7-23
How to generate a triangular mesh by MATLAB with the line joining the cell centre is orthogonal to the edges of the triangle?

请先登录，再进行评论。

更多回答（2 个）

Devosmita Chatterjee
How to refine a mesh for the poisson equation problem in a unitsquare domain by placing many triangles around the point (0.5,0.5) using Adaptive Finite Element Method and a posteriori error estimation? We need to find adaptive meshes for the problem with solution u=a*exp(-a*(r^2)), a =400 using two refinements, four refinements, six refinements, ten refinements. Can anyone provide me the Matlab code for this problem?
Ian Jentz
Ian Jentz 2020-10-2
编辑：Ian Jentz 2020-10-2
There is now a port of refinemesh() available for PDEModel. Like refinemesh(), it works for 2D linear triangular geometries only.
The new PDEModel compatible function is refinePDEMmesh().
This is likely how would run your problem:
create your model as a pde model
model=createpde;
proceed with your geometry definitions.
When you get to meshing, you instantiate your mesh as usuall, but then can refine the mesh on a particular face using refinePDEMmesh().
generateMesh(model);
FaceID = 1; % the face ID of your E1 face?
model = refinePDEMmesh(model,FaceID);
  1 个评论
Sean Littleton
Sean Littleton 2021-1-15
编辑：Sean Littleton 2021-1-15
Hello,
Thank you for sharing this function. It is almost exactly what I need, except for the fact that it does not work with quadratic triangular elements -- the very elements that I need to use for my problem. Are there any plans to add this functionality in the future? If not, is there an approach you would recommend?
Note: The triangular elements I am using are still geometrically linear. The extra nodes are only used for the quadratic shape functions, since I need the derivatives of the solution to be continuous.

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 Geometry and Mesh 的更多信息

标签

产品

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by