Dear all, I can create unstructured mesh by myself, and mark node index on boundary. I use the Poisson example code, to get result file, “poisson000000.vtu”. The figure is show by paraview. You can get node inde (id = 1055) in the figure. [螢幕擷取畫面 2022-01-20 190920] In the example code, we can…

I am a novice Fenics user. Thanks for getting back to me. Do you mean that? class Gammathree(SubDomain): def inside(self, x, on_boundary): return near(x[0], 0, DOLFIN_EPS) and near(x[1], 0, DOLFIN_EPS) and near(x[2], 0, DOLFIN_EPS) class Gammatwo(SubDomain): def inside(self, x, on_…

You cannot use the node index, as the index will not necessary be the same in dolfin (as dolfin might renumber for locallity, ordering and parallel distribution). You could mark the node (for instance by a geometric search as shown above) and then use the pointwise method inside the DirichletBC.

Dear @dokken , I totally understand your commentary. As the red dot on the picture, I want to set the DirichletBC on such unstructured mesh, and the mesh density is not consistent. Is there a good way to provide me? [USMfig]

Again, as you know the coordinate, create a condition similar to: class Gammatwo(SubDomain): def inside(self, x, on_boundary): return near(x[0], 375, DOLFIN_EPS) and near(x[1], -52, DOLFIN_EPS) and near(x[2], 409.051, DOLFIN_EPS) point2 = Gammatwo() bc = DirichletBC(V.sub(0), Constant(0…

How to setting Dirichlet boundary condition by node index?

I/O

dokken January 20, 2022, 1:29pm 2

You can have a look at: Boundary conditions on edges or nodes - #6 by dokken

Topic		Replies	Views
Boundary conditions for specific nodes Uncategorized	1	315	August 11, 2022
Boundary condition based on mesh index mesh	0	398	July 12, 2020
Extra argument in def (x,on_boundary) Uncategorized	2	397	April 26, 2021
Apply a DirichletBC to a single node mesh	11	4629	July 26, 2019
Boundary conditions on edges or nodes Uncategorized	10	4672	February 1, 2021

How to setting Dirichlet boundary condition by node index?

Related topics