Hello everybody,
I'm currently working on a 3D eddy current problem regarding a conducting cube in homogenous magnetic field in the time domain. While doing this some questions arose about the correct implementation of non-homogeneous Dirichlet boundary conditions.
I tried the following so far:
Code:
coeff_dirichlet = CoefficientFunction((
CoefficientFunction([0,0,0,0,0,-25e-3*y*B0*sin(omega * t)]),
CoefficientFunction([0,0,0,25e-3*B0*x*sin(omega * t),0,0]),
CoefficientFunction([0,0,0,0,0,0])))
uh.Set(coeff_dirichlet, BND)
with the entries in the list being the different surfaces (1: bottom, 2: top, 3: left, 4: right, 5: front, 6: back). Drawing these onto the mesh shows in the GUI that the values are set to the correct surfaces.
Even though there is a lot of well written documentation I don't fully understand how to mark the boundaries correctly as the keyword
BND tells
Set that the CoefficientFunction need only be interpolated at the marked boundaries.
Note the attached code. Any help is very much appreciated as I'm pretty new to FEM-Simulation in general.
Best
Sebastian