Avanzado Condiciones de Borde Laplace

[Tempa]

Hola. Estoy resolviendo en C la ecuacion de Laplace para modelar el patron de corriente dentro de un canal MHD. Estoy usando la funcion de correinte, lo cual hace el problema bastante simplificado. Aqui tengo un esquema de mi area:



　　　　　　　Electrode
0.--0--0--0--0--0--0--0-==========-1--1--1--
0. 1
(Plasma Flow) --------> 　　 xxxxxxx (x =B Magnetic Field)
0. 　　 xxxxxxx 1
0. 1
--0--0--0--0--0--0--0--==========-1--1--1-- 　　　　　　　　Electrode

Las condiciones de borde son tales que antes del electrodo la funcion valga 0, y despues 1. Las paredes verticales de entrada y salida valen 0 y 1, respectivamente. Y los electrodos tienen una funcion de manera que el valor vaya de 0 a 1 a medida que avance x. Mi pregunta viene aqui. Como puedo hacer para introducir el efecto de un campo magnetico aplicado en el medio del canal, por ejemplo?. Dentro del canal esta la ecuaion clasica de diferencias finitas, y puedo ver la distribucion de correinte cuando grafico. Pero quisiera introducir el efecto del campo magnetico. Estaba oensando en unas condiciones de borde de Neumann en esos untos interiores, de manera que la ecuacion en esos ountos tendria que cambiarla.

Alguna sugerencia?