From 71416d0bcf354a8aec9d1e8be0dc9769e5c4a76d Mon Sep 17 00:00:00 2001 From: Bernd Hahnebach Date: Thu, 21 May 2020 16:53:14 +0200 Subject: [PATCH] FEM: constants, fix unit of permitivity of vacuum --- src/Mod/Fem/femsolver/elmer/writer.py | 3 ++- src/Mod/Fem/femtools/constants.py | 4 +++- 2 files changed, 5 insertions(+), 2 deletions(-) diff --git a/src/Mod/Fem/femsolver/elmer/writer.py b/src/Mod/Fem/femsolver/elmer/writer.py index 16aec0cfd9..f877588e6a 100644 --- a/src/Mod/Fem/femsolver/elmer/writer.py +++ b/src/Mod/Fem/femsolver/elmer/writer.py @@ -340,7 +340,8 @@ class Writer(object): def _handleElectrostaticConstants(self): self._constant( "Permittivity Of Vacuum", - getConstant("PermittivityOfVacuum", "T^4*I^2/(L*M)")) + getConstant("PermittivityOfVacuum", "T^4*I^2/(L^3*M)")) + # https://forum.freecadweb.org/viewtopic.php?f=18&p=400959#p400959 def _handleElectrostaticMaterial(self, bodies): for obj in self._getMember("App::MaterialObject"): diff --git a/src/Mod/Fem/femtools/constants.py b/src/Mod/Fem/femtools/constants.py index b0afb7ae4c..11a138777a 100644 --- a/src/Mod/Fem/femtools/constants.py +++ b/src/Mod/Fem/femtools/constants.py @@ -41,7 +41,9 @@ def stefan_boltzmann(): def permittivity_of_vakuum(): - return "8.8542e-12 s^4*A^2/(m*kg)" + # https://forum.freecadweb.org/viewtopic.php?f=18&p=400959#p400959 + return "8.8542e-12 s^4*A^2 / (m**3 * kg)" + def boltzmann_constant():