On Thu, 5 Dec 2019, Prashant Kumar Jha wrote:
> I can create a network mesh consisting of edge elements in 3D space. This > is non-standard mesh in the sense that I start with empty mesh, add couple > of nodes, add edge element between nodes. Each node has 1 pressure dof. On > this network, I am now trying to solve 1D pde for pressure. For boundary > conditions, I want to fix the two ends of the network (i.e. the first dof > and the last dof) to a specific value. > > DirichletBoundary with boundary id 0 and 1 result in error possibly due to > mesh in 3D where as elements are of 1D nature. That shouldn't be the case. What exactly is the error? Did you set the boundary ids when you created the mesh? > To get around this, I was wondering if I could fix the individual dof to a > specific value. Any suggestions to make above work is appreciated. You can manually add a constraint row for each boundary dof, or you can use a penalty method. --- Roy _______________________________________________ Libmesh-users mailing list Libmesh-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-users