Correct way of using symbolicenergy for Neumann boundary condition

NGSolve
1

Hello

I want to use symbolicenergy for a distributed force on a Neumann boundary. For example, if my force is ‘f’ and my boundary is named ‘surface’, it is correct to use:

fes = VectorH1(mesh, order=3)
u, v = fes.TnT()
a= BilinearForm(fes)
a+= SymbolicEnergy(f*specialfc.normal(3)*u, definedon= ‘surface’)

Should I use AssembleLinearization?

Regards
Hamed

2

preferred way is now:

a += Variation(f * n * u * ds("surface"))

but the difference is more syntactic sugar. in symb energy I think you still need to specify BND integration somewhere

best,
Christopher

3

Thanks for your prompt reply!