Electrostatic solver

.gitignore

@@ -79,4 +79,4 @@ femmt/examples/paper_thermal_validation.py
 femmt/examples/example_results
 
 # pyspelling
-dictionary.dic
+dictionary.dic

docs/wordlist

@@ -452,9 +452,79 @@ NbSteps
 
 
 
-
-
-
+# electrostatic
+Kapton
+Capacitances
+capacitances
+pF
+ADE
+ACvsBDE
+ADvsBCE
+DvsABCE
+CvsABDE
+BvsACDE
+AvsBCDE
+ABEvsCD
+ABCDvsE
+ABvsCDE
+CDE
+NaN
+xlsx
+InsulationMaterial
+CoreMaterial
+ZeroVoltage
+MagD
+MagE
+EleSta
+airgaps
+OutBoundary
+DirResValsVoltage
+DirResValsCapacitance
+DirResValsCapacitanceFromQV
+DirResValsCharge
+QV
+DirResValsTurn
+Tetrafluoroethylene
+ETFE
+Fluorinated
+FEP
+Propylene
+ECTFE
+Ethylenechlorotrifluoroethylene
+PVDF
+Polyvinylidene
+Elastomers
+TPE
+Polyamide
+PBT
+Polybutylene
+Terephthalate
+PUR
+CPE
+TPR
+Polychloroprene
+Butadiene
+SBR
+Styrene
+EPR
+SBR
+CSPE
+Chlorosulfonated
+fluoropolymer
+PFA
+Perfluoroalkoxy
+PTFE
+Polytetrafluoroethylene
+LCP
+DAP
+Diallyl
+Phthalate
+Nomex
+Phenolic
+bOPET
+xlsxwriter
+Thermoset
+Polyphenylene
 
 # logging
 FEMMTLogger

femmt/data.py

@@ -92,6 +92,12 @@ def update_paths(self, working_directory: str, electro_magnetic_folder_path: str
 
         # Setup file paths
         self.e_m_results_log_path = os.path.join(self.results_folder_path, "log_electro_magnetic.json")
+        ###
+        # for electrostatic
+        self.electrostatic_results_log_path = os.path.join(self.results_folder_path, "log_electro_static.json")
+        self.capacitance_result_log_path = os.path.join(self.results_folder_path, "capacitance_result.json")
+        self.capacitance_matrix_path = os.path.join(self.results_folder_path, "log_capacitance_matrix.json")
+        ###
         self.coordinates_description_log_path = os.path.join(self.results_folder_path, "log_coordinates_description.json")
         self.reluctance_log_path = os.path.join(self.results_folder_path, "log_reluctance_and_inductance.json")
         self.material_log_path = os.path.join(self.results_folder_path, "log_material.json")

femmt/dtos.py

@@ -22,6 +22,14 @@ class StackedCoreDimensions:
     window_h_top: float
     window_h_bot: float
 
+@dataclass
+class BobbinDimensions:
+    """Defines the dimensions of a default core."""
+
+    bobbin_inner_diameter: float
+    bobbin_window_w: float
+    bobbin_window_h: float
+    bobbin_h: float
 
 @dataclass
 class ConductorRow:

femmt/electro_magnetic/fields_electrostatic.pro

@@ -0,0 +1,25 @@
+PostOperation Map_local UsingPost EleSta {
+
+  // Potentials for the entire domain
+  ExtGmsh     = ".pos";
+  Print[ u0, OnElementsOf Region[{Domain}], Name "Potential / V", File StrCat[DirResFields, "Potential", ExtGmsh], LastTimeStepOnly ] ;
+  Print[ u0, OnElementsOf Region[{Core}], Name "Average Potential on Core / V", File StrCat[DirResFields, "Voltage_Core_Map", ExtGmsh], LastTimeStepOnly ];
+
+  // Electric Field vector in the entire domain
+  Print[ e, OnElementsOf Region[{Domain}], Name "Electric Field", File StrCat[DirResFields, "Efield", ExtGmsh], LastTimeStepOnly ] ;
+  //Print[ Welocal, OnElementsOf Region[{Domain}], Name "Stored Energy", File StrCat[DirResFields, "We", ExtGmsh], LastTimeStepOnly ] ;
+  Print[ MagE, OnElementsOf Region[{Domain}], Name "Magnitude Electric Field / V/m", File StrCat[DirResFields, "MagE", ExtGmsh], LastTimeStepOnly ] ;
+
+  // Displacement Field vector in the entire domain
+  Print[ d, OnElementsOf Region[{Domain}], Name "Electric Field Density", File StrCat[DirResFields, "Dfield", ExtGmsh], LastTimeStepOnly ] ;
+  Print[ MagD, OnElementsOf Region[{Domain}], Name "Magnitude Electric Field Density / C/m^2", File StrCat[DirResFields, "MagD", ExtGmsh], LastTimeStepOnly ] ;
+
+
+  // Settings for visualization output (optional)
+  Echo[ Str["View[PostProcessing.NbViews-1].Light=0;
+             View[PostProcessing.NbViews-1].LineWidth = 2;
+             View[PostProcessing.NbViews-1].RangeType=3;
+             View[PostProcessing.NbViews-1].IntervalsType=1;
+             View[PostProcessing.NbViews-1].NbIso = 25;"],
+           File OptionPos];
+}