Example 10 opuseads heatfluid

examples/antora.yml

@@ -1,8 +1,9 @@
 name: cases
 title: Feel++ Toolbox Cases Guide
-version: '0.106'
+version: '0.106-case-10-opuseads-heatfluid'
 nav:
 - modules/ROOT/nav.adoc
+- modules/poisson/nav.adoc
 - modules/csm/nav.adoc
 - modules/cfd/nav.adoc
 - modules/heat/nav.adoc

examples/modules/heatfluid/examples/opus/eads.geo

@@ -0,0 +1,107 @@
+h = 0.0001;
+// h     = DefineNumber[ 0.1, Name "h"     ];
+e_pcb = DefineNumber[ 2e-3, Name "e_pcb" ]; // 0.3
+e_a   = DefineNumber[ 4e-3, Name "e_a"   ]; // 0.7
+e_ic  = DefineNumber[ 2e-3, Name "e_ic"  ]; // 0.3
+h_1   = DefineNumber[ 2e-2, Name "h_1"   ]; // 0.2
+h_2   = DefineNumber[ 7e-2, Name "h_2"   ]; // 0.6
+h_ic  = DefineNumber[ 2e-2, Name "h_ic"  ]; // 0.2
+h_pcb = DefineNumber[13e-2, Name "h_pcb" ]; // 1
+
+// h_pad = DefineNumber[ 0.04, Name "h_pad" ];
+
+
+Point(1)  = {0, 0, 0, h};
+Point(2)  = {e_pcb, 0, 0, h};
+Point(3)  = {e_pcb + e_a, 0, 0, h};
+Point(4)  = {e_pcb, h_1, 0, h};
+Point(5)  = {e_pcb+e_ic, h_1, 0, h};
+Point(6)  = {e_pcb+e_ic, h_1+h_ic, 0, h};
+Point(7)  = {e_pcb, h_1+h_ic, 0, h};
+Point(8)  = {e_pcb, h_2, 0, h};
+Point(9)  = {e_pcb+e_ic, h_2, 0, h};
+Point(10) = {e_pcb+e_ic, h_2+h_ic, 0, h};
+Point(11) = {e_pcb, h_2+h_ic, 0, h};
+Point(12) = {0, h_pcb, 0, h};
+Point(13) = {e_pcb+e_a, h_pcb, 0, h};
+Point(14) = {e_pcb, h_pcb, 0, h};
+Point(15) = {e_pcb+e_ic, 0, 0, h};
+Point(16) = {e_pcb+e_ic, h_pcb, 0, h};
+
+// Point(17) = {e_pcb, h_pcb + h_pad, 0, h};
+// Point(18) = {e_pcb+e_a, h_pcb + h_pad, 0, h};
+
+
+Line(1)  = {4, 7};
+Line(2)  = {8, 11};
+Line(3)  = {1, 12};
+Line(4)  = {3, 13};
+Line(5)  = {2, 4};
+Line(6)  = {7, 8};
+Line(7)  = {11, 14};
+Line(8)  = {4, 5};
+Line(9)  = {5, 6};
+Line(10) = {6, 7};
+Line(11) = {8, 9};
+Line(12) = {9, 10};
+Line(13) = {10, 11};
+Line(14) = {15, 5};
+Line(15) = {6, 9};
+Line(16) = {10, 16};
+Line(17) = {1, 2};
+Line(18) = {12, 14};
+Line(19) = {2, 15};
+Line(20) = {14, 16};
+Line(21) = {15, 3};
+Line(22) = {16, 13};
+
+// Line(37) = {14, 17};
+// Line(38) = {17, 18};
+// Line(39) = {18, 13};
+
+
+Line Loop(23) = {17, 5, 1, 6, 2, 7, -18, -3};
+Plane Surface(24) = {23};
+Line Loop(25) = {13, -2, 11, 12};
+Plane Surface(26) = {25};
+Line Loop(27) = {10, -1, 8, 9};
+Plane Surface(28) = {27};
+Line Loop(29) = {20, -16, 13, 7};
+Plane Surface(30) = {29};
+Line Loop(31) = {11, -15, 10, 6};
+Plane Surface(32) = {31};
+Line Loop(33) = {8, -14, -19, 5};
+Plane Surface(34) = {33};
+Line Loop(35) = {21, 4, -22, -16, -12, -15, -9, -14};
+Plane Surface(36) = {35};
+
+// Line Loop(40) = {37, 38, 39, -22, -20};
+// Plane Surface(41) = {40};
+
+
+Physical Surface("Air")  = {30,32,34,36}; // 41
+Physical Surface("Pcb")  = {24};
+Physical Surface("IC1")  = {28};
+Physical Surface("IC2")  = {26};
+
+// Physical Line("Fixed")   = {21, 19, 17};
+
+Physical Line("Pcb/Input")  = {17};
+Physical Line("Pcb/Left")   = {3};
+Physical Line("Pcb/Output") = {18};
+
+// Physical Line("Upwards") = {18,20,22};
+
+Physical Line("IC1/Pcb") = {1};
+Physical Line("IC2/Pcb") = {2};
+
+// Physical Line("Air/Pcb") = {3,7,6,5,4};
+// Physical Line("Air/IC1") = {10,9,8};
+// Physical Line("Air/IC2") = {13,12,11};
+
+// Physical Line("Air/PostOutput") = {38};
+Physical Line("Air/Input")         = {21,19};
+Physical Line("Air/Walls")         = {7,6,5,10,9,8,13,12,11}; // 37,39
+Physical Line("Air/Right")         = {4};
+Physical Line("Air/Output")        = {22,20};
+

examples/modules/heatfluid/examples/opus/eads_normal.json

@@ -0,0 +1,99 @@
+// -*- mode: javascript -*-
+{
+    "Name"      : "Eads",
+    "ShortName" : "Eads",
+
+    "Models" : {
+	"use-model-name" : 1,
+        "fluid" : { "equations" : "Navier-Stokes" }
+    },
+
+    "Materials" : {
+        "Pcb" : {
+	    "markers" : "Pcb",
+            "physics" : "heat",
+            "rho"     : "1.",
+            "k"       : "0.2",
+            "Cp"      : "2e6",
+	    "mu"      : "1e12"
+        },
+	"IC1" : {
+	    "markers" : "IC1",
+	    "physics" : "heat",
+	    "name"    : "IC1",
+	    "rho"     : "1.",
+	    "k"       : "2",
+	    "Cp"      : "1.4e6",
+	    "mu"      : "1e12"
+	},
+	"IC2" : {
+	    "markers" : "IC2",
+	    "physics" : "heat",
+	    "rho"     : "1.",
+	    "k"       : "2",
+	    "Cp"      : "1.4e6",
+	    "mu"      : "1e12"
+	},
+	"Air" : {
+	    "markers" : "Air",
+	    "physics" : ["heat","fluid"],
+	    "rho"     : "1.",
+	    "k"       : "3e-2",
+	    "Cp"      : "1100",
+	    "mu"      : "1.8e-5"
+	}
+    },
+
+    "BoundaryConditions" : {
+	"velocity" : {
+	    "Dirichlet" : {
+		"Air/Right"           : { "expr" : "{0,0}" },
+		"Air/Walls"           : { "expr" : "{0,0}" },
+		"Air/Input"           : { "expr" : "{0,(2e-3+2e-3<x)*(x<2e-3+4e-3)*7e-3*(3/2/(4e-3-2e-3))*(1-((x-(4e-3+2e-3)/2-2e-3)/((4e-3-2e-3)/2))^2)*(1-exp(-t/3))}:x:t" }
+	    }
+	},
+	"fluid" : {
+	    "outlet" : {
+		"Air/Output"          : { "expr" : "0." }
+	    }
+	},
+        "temperature" : {
+            "Dirichlet" : {
+                "Air/Input"           : { "expr" : "300" },
+		"Pcb/Input"           : { "expr" : "300" }
+            },
+            "Neumann" : {
+                "Air/Output"          : { "expr" : "0." },
+                "Pcb/Output"          : { "expr" : "0." },
+		"Air/Right"           : { "expr" : "0." },
+		"Pcb/Left"            : { "expr" : "0." }
+            },
+	    "VolumicForces" : {
+		"IC1"                 : { "expr" : "1e6*(1-exp(-t)):t" },
+		"IC2"                 : { "expr" : "1e6*(1-exp(-t)):t" }
+	    }
+	}
+    },
+
+    "InitialConditions" : {
+	"temperature" : {
+	    "" : {
+		"" : { "expr" : "300" }
+	    }
+	}
+    },
+
+    "PostProcess" : {
+	"use-model-name" : 1,
+	"heat-fluid" : {
+	    "Exports" : {
+		"fields" : [
+		    "fluid.velocity","fluid.pressure",
+		    "heat.temperature","fluid.pid"
+		]
+	    }
+	},
+	"fluid" : {}, "heat" : {}
+    }
+}
+

examples/modules/heatfluid/examples/opus/test1.cfg

@@ -0,0 +1,28 @@
+directory=toolboxes/heatfluid/opus/test1
+
+case.dimension=2
+
+[heat-fluid]
+filename=$cfgdir/eads_normal.json
+mesh.filename=$cfgdir/eads.geo
+gmsh.hsize=5e-4
+
+[heat-fluid.heat]
+initial-solution.temperature=300
+stabilization-gls=1
+bdf.order=1
+
+# gmsh.partition=1
+
+[heat-fluid.fluid]
+bdf.order=1
+# solver=Newton
+pc-type=lu
+
+# gmsh.partition=1
+
+[ts]
+time-step=1
+time-final=1500
+# restart.at-last-save=true
+

examples/modules/heatfluid/nav.adoc

@@ -1,6 +1,7 @@
 * xref:README.adoc[Heat & Fluid]
 ** xref:2Dbuilding_NS/README.adoc[2D Building with natural convection]
 ** xref:cabin/README.adoc[Aerothermal in an Airplane Cabin]
+** xref:opus/README.adoc[Cooling of electronic components]
 ** xref:TurekHron/README.adoc[Forced heat convection around a cylinder]
 ** xref:naturalconvection/README.adoc[Natural heat convection in a cavity]
 ** xref:pipestokesmass/README.adoc[Mass transport with a Stokes flow in a pipe]