How we can mark interior nodes and impose some condition on those nodes and then call it in solve()

[dokken]

See for instance: Boundary conditions on edges or nodes - #6 by dokken on how to apply BCs to a function, which in turn can be visualized in Paraview