From ac53191210ddfc9f1d8e4158f0b0a6b7b1ce9585 Mon Sep 17 00:00:00 2001
From: Christophe Geuzaine <cgeuzaine@uliege.be>
Date: Thu, 14 Jan 2021 16:48:02 +0100
Subject: [PATCH] option for freq or time domain
---
Magnetodynamics/Lib_Magnetodynamics2D_av_Cir.pro | 4 ++++
Magnetodynamics/transfo.pro | 7 ++++---
2 files changed, 8 insertions(+), 3 deletions(-)
diff --git a/Magnetodynamics/Lib_Magnetodynamics2D_av_Cir.pro b/Magnetodynamics/Lib_Magnetodynamics2D_av_Cir.pro
index 716a045..846d410 100644
--- a/Magnetodynamics/Lib_Magnetodynamics2D_av_Cir.pro
+++ b/Magnetodynamics/Lib_Magnetodynamics2D_av_Cir.pro
@@ -556,6 +556,10 @@ PostProcessing {
Term { [ {d a} ]; In Vol_Mag; Jacobian Vol; }
}
}
+ { Name norm_of_b; Value {
+ Term { [ Norm[{d a}] ]; In Vol_Mag; Jacobian Vol; }
+ }
+ }
{ Name h; Value {
Term { [ nu[] * {d a} ]; In Vol_Mag; Jacobian Vol; }
}
diff --git a/Magnetodynamics/transfo.pro b/Magnetodynamics/transfo.pro
index 3ec7ac7..d8d3ebb 100644
--- a/Magnetodynamics/transfo.pro
+++ b/Magnetodynamics/transfo.pro
@@ -21,7 +21,10 @@ DefineConstant[
// way as far as circuit-coupling is concerned
ConductorType = {2, Choices{1 = "Massive", 2 = "Coil"}, Highlight "Blue",
Name "Parameters/01Conductor type"}
- Freq = {1, Min 0, Max 1e3, Step 1, Name "Parameters/Frequency"}
+ Freq = {1, Min 0, Max 1e3, Step 1,
+ Name "Parameters/Frequency"}
+ Flag_FrequencyDomain = {0, Choices{0, 1},
+ Name "Parameters/Frequency-domain?"}
mur_Core = {1000, Min 1, Max 10000, Step 1,
Name "Parameters/Core relative permeability"}
];
@@ -94,8 +97,6 @@ Function {
// secondary.
Flag_CircuitCoupling = 1;
-Flag_FrequencyDomain = 0;
-
// Note that the voltage will not be equally distributed in the PLUS and MINUS
// parts, which is the reason why we must apply the total voltage through a
// circuit -- and we cannot simply use a current source like in Tutorial 7a.
--
GitLab