electrostatic problems with ArcaneFEMFunctions

arcaneframework · Aug 22, 2024 · 3874cc5 · 3874cc5
1 parent 09feed6
commit 3874cc5
Show file tree

Hide file tree

Showing 4 changed files with 277 additions and 470 deletions.
diff --git a/electrostatics/Fem.axl b/electrostatics/Fem.axl
@@ -9,9 +9,6 @@
     <variable field-name="phi" name="Phi" data-type="real" item-kind="node" dim="0">
       <description>FEM variable phi on nodes for electrostatic potential</description>
     </variable>
-    <variable field-name="phi_dirichlet" name="PhiDirichlet" data-type="bool" item-kind="node" dim="0">
-      <description>Boolean which is true if Dirichlet node</description>
-    </variable>
     <variable field-name="E" name="E" data-type="real3" item-kind="cell" dim="0">
       <description>electric field vector on cells</description>
     </variable>
@@ -32,16 +29,6 @@
     <simple name="mesh-type" type="string"  default="TRIA3" optional="true">
       <description>Type of mesh provided to the solver</description>
     </simple>
-    <simple name = "enforce-Dirichlet-method" type = "string" default="Penalty" optional="true">
-      <description>
-        Method via which Dirichlet boundary condition is imposed
-      </description>
-    </simple>
-    <simple name = "penalty" type = "real" default="1.e30" optional="true">
-      <description>
-        Penalty value for enforcing Dirichlet condition
-      </description>
-    </simple>
 
     <!-- - - - - - dirichlet-boundary-condition - - - - -->
     <complex name  = "dirichlet-boundary-condition"
@@ -62,6 +49,47 @@
           Value of the boundary condition
         </description>
       </simple>
+      <simple name = "enforce-Dirichlet-method" type = "string" default="Penalty" optional="true">
+        <description>
+          Method via which Dirichlet boundary condition is imposed
+        </description>
+      </simple>
+      <simple name = "penalty" type = "real" default="1.e30" optional="true">
+        <description>
+          Penalty value for enforcing Dirichlet condition
+        </description>
+      </simple>
+    </complex>
+
+    <!-- - - - - - manufactured-solution - - - - -->
+    <complex name  = "manufactured-solution"
+             type  = "ManufacturedSolution"
+             minOccurs = "0"
+             maxOccurs = "unbounded"
+      >
+      <description>
+        Dirichlet boundary condition
+      </description>
+      <simple name="manufactured-dirichlet" type="bool"  default="false" optional="true">
+        <description>
+          Function for Dirichlet boundary condition
+        </description>
+      </simple>
+      <simple name = "enforce-Dirichlet-method" type = "string" default="Penalty" optional="true">
+        <description>
+          Method via which Dirichlet boundary condition is imposed
+        </description>
+      </simple>
+      <simple name = "penalty" type = "real" default="1.e30" optional="true">
+        <description>
+          Penalty value for enforcing Dirichlet condition
+        </description>
+      </simple>
+      <simple name="manufactured-source" type="bool"  default="false" optional="true">
+        <description>
+          Function for manufactured source term condition
+        </description>
+      </simple>
     </complex>
 
     <!-- - - - - - neumann-boundary-condition - - - - -->
@@ -115,6 +143,16 @@
           Value of the point Dirichlet condition
         </description>
       </simple>
+      <simple name = "enforce-Dirichlet-method" type = "string" default="Penalty" optional="true">
+        <description>
+          Method via which Dirichlet boundary condition is imposed
+        </description>
+      </simple>
+      <simple name = "penalty" type = "real" default="1.e30" optional="true">
+        <description>
+          Penalty value for enforcing Dirichlet condition
+        </description>
+      </simple>
     </complex>
 
     <service-instance name = "linear-system"