[riddim]

Hey there,

I think I found out that I "simply" need to write my own laplacian scheme with boundary conditions ....

http://www.cfd-online.com/Forums/ope...te-scheme.html <<< this threat was very useful to me but I have to admit I'm still a little bit stuck at the start ...
...maybe you could give me some hints how to implement this ...

My aim is to simply solve

Laplacian(A) == -µ_r*µ_0*J

µ_0 is const. and µ_r can be different from cell to cell

A * n == 0
and (A x curl(A)) * n == 0

where n is the surface normal vector.

do I have to implement these boundary conditions inside the laplacian scheme or is it possible to use the standard laplacian scheme and implement these conditions inside my solver ...?
How do i access the surface normal vector?

Where do you think I should start ...?