From 8ddbd0facad1c6ac10a18c2ada371722b6c44c80 Mon Sep 17 00:00:00 2001
From: Nicolas Marsic <marsic@temf.tu-darmstadt.de>
Date: Tue, 12 Sep 2017 12:29:14 +0200
Subject: [PATCH] Real operator induces a mix between source-free and driven
problem: mathematically speaking we need the full complex omega
---
example/maxwell_sibc/sphere.pro | 16 +++++++---------
1 file changed, 7 insertions(+), 9 deletions(-)
diff --git a/example/maxwell_sibc/sphere.pro b/example/maxwell_sibc/sphere.pro
index 51517b2..43abd03 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{
--
GitLab