Hello,
I have a simulation of welding two pipes. I have a heat source within the mesh at the junction of two pipes. The heat source is switched off at a certain time, and cooling of the material occurs.
I have the code:
class Q(Expression):
def __init__(self, T0, t):
self._T0 = T0
self._t = t
def eval(self, values, x):
if(t < 50):
values[0] = ...
else:
values[0] = 0
a = ...
L = ... + Q*v*ds(20)
where Q - described the Neumann condition. This does not work. Maybe because the Q is within the network between the surfaces?
How do I implement this in my code?
Sincerely, Boris.
