i have read an post in this forum about setting up the boundary condition only at a certain nod. Unfortunatly I don’t know how to get the id or name of the nod form a 3d mesh grid. Would aprreciate any help.
mesh = FE.get_CubeMesh(size=0.1, mesh_split=split_mesh)
fes = VectorH1(mesh, order=degree_el, dirichlet_bbnd=name_Boundarienod)
‘3d cube mesh gets defined here’
def get_CubeMesh(size=0.1, mesh_split=2)
return get_Mesh(size, sizex=size, sizey=size, sizez=size, mesh_split=mesh_split)