Group Abstract Group Abstract

Message Boards Message Boards

Intersecting boundary problems when solving PDEs

POSTED BY: Marco Thiel
5 Replies
Attachments:

Hi Gregory,

thanks a lot. Indeed, now it works:

mr = Import["~/Desktop/cowfix4.stl", "BoundaryMeshRegion"];
uif = NDSolveValue[{Inactive[Laplacian][
      u[x, y, z], {x, y, z}] == -0.01, 
    DirichletCondition[u[x, y, z] == 1, z <= 0.08]}, 
   u, {x, y, z} \[Element] mr];

Needs["NDSolve`FEM`"]
ElementMeshSurfacePlot3D[uif, Boxed -> False, ViewPoint -> {0, -4, 2}]

enter image description here

I noticed that you had to trim the cow's tail.

Thanks,

Marco

POSTED BY: Marco Thiel

Thank you very much for your help! That is very useful. I will check TetGenDetectIntersectingFacets out and see whether I can make it work. Thanks a lot!

I just wonder one thing. Some time ago I built a little model of cells in a "microfluidic chamber" - the details don't really matter. The idea is that you have "holes" in a rectangular shape. These holes have their own boundary conditions and you want to solve a diffusion equation. The simulation seemed to work fine in the MMA 10 prerelease version (see results here!) and it stopped working in the full version. If you execute the file that you can download using the link, and you execute it in MMA10 you will get an error message that there are intersecting boundaries and the execution fails.

Marco

PS: I am sorry that the notebook is not commented properly, but it should give an idea of the observation.

POSTED BY: Marco Thiel
POSTED BY: Bruce Miller
Reply to this discussion
Community posts can be styled and formatted using the Markdown syntax.
Reply Preview
Attachments
Remove
or Discard