diff --git a/example/maxwell_sibc/sphere.pro b/example/maxwell_sibc/sphere.pro index 51517b2c56a357d16baae4f32928f8ab0f625199..43abd03e6e94fe66377021fcc0acc98a6f1bab94 100644 --- a/example/maxwell_sibc/sphere.pro +++ b/example/maxwell_sibc/sphere.pro @@ -16,17 +16,15 @@ Function{ Mu0 = 4*Pi*1e-7; // [H/m] Eps0 = 1/(Mu0 * C0^2); // [F/m] - // Conductivity (coming from sphere.dat) // - If(isSC == 1) // Superconductor - Sigma[] = Sigma - J[] / (Mu0 * Lambda^2 * Re[Puls[]]); // [S/m] + // Leontovich SIBC (H-Formulation) // + If(isSC == 0) + SL[] = Sqrt[(J[] * Puls[] * Mu0) / Sigma]; // H-Formulation EndIf - If(isSC == 0) // Normal conductor - Sigma[] = Sigma; // [S/m] + If(isSC == 1) + SL[] = 0.5 * Puls[]^2 * Mu0^2 * Lambda^3 * Sigma + + J[] * Puls[] * Mu0 * Lambda; // H-Formulation EndIf - - // Leontovich SIBC (H-Formulation) // - SL[] = Sqrt[(J[] * Re[Puls[]] * Mu0) / Sigma[]]; // H-Formulation - SIbc[] = -1.0 / SL[]; // E-Formulation + SIbc[] = -1.0 / SL[]; // E-Formulation } Integration{