.. DO NOT EDIT. .. THIS FILE WAS AUTOMATICALLY GENERATED BY SPHINX-GALLERY. .. TO MAKE CHANGES, EDIT THE SOURCE PYTHON FILE: .. "examples/Thermal/Thermal3.py" .. LINE NUMBERS ARE GIVEN BELOW. .. only:: html .. note:: :class: sphx-glr-download-link-note :ref:`Go to the end ` to download the full example code. .. rst-class:: sphx-glr-example-title .. _sphx_glr_examples_Thermal_Thermal3.py: Thermal3 ======== Transient thermal simulation. .. GENERATED FROM PYTHON SOURCE LINES 11-93 .. image-sg:: /examples/Thermal/images/sphx_glr_Thermal3_001.png :alt: Thermal3 :srcset: /examples/Thermal/images/sphx_glr_Thermal3_001.png :class: sphx-glr-single-img .. image-sg:: /examples/Thermal/images/sphx_glr_Thermal3_002.gif :alt: Thermal3 :srcset: /examples/Thermal/images/sphx_glr_Thermal3_002.gif :class: sphx-glr-single-img .. rst-class:: sphx-glr-script-out .. code-block:: none 0.000 s 0.100 s 0.200 s 0.300 s 0.400 s 0.500 s 0.600 s 0.700 s 0.800 s 0.900 s 1.000 s 1.100 s 1.200 s 1.300 s 1.400 s 1.500 s 1.600 s 1.700 s 1.800 s 1.900 s 2.000 s 2.100 s 2.200 s 2.300 s 2.400 s 2.500 s 2.600 s 2.700 s 2.800 s 2.900 s 3.000 s 3.100 s 3.200 s 3.300 s 3.400 s 3.500 s 3.600 s 3.700 s 3.800 s 3.900 s 4.000 s 4.100 s 4.200 s 4.300 s 4.400 s 4.500 s 4.600 s 4.700 s 4.800 s 4.900 s 5.000 s 5.100 s ==================== Mesh ==================== Element type: HEXA8 Ne = 940, Nn = 1584 ==================== Model ==================== Thermal : thermal conductivity (k) : 1 thermal mass capacity (c) : 1 solver : pypardiso ============= Boundary Conditions ============= Unspecified. =================== Results =================== Unspecified. =================== TicTac =================== Mesh : 80.746 ms Boundary Conditions : 29.564 µs Matrix : 24.198 ms Solver : 1.629 s PyVista_Interface : 12.734 s Generate movie 0/51 Generate movie 1/51 (1.96 %) 7.09 s Generate movie 2/51 (3.92 %) 6.91 s Generate movie 3/51 (5.88 %) 6.75 s Generate movie 4/51 (7.84 %) 6.61 s Generate movie 5/51 (9.80 %) 6.27 s Generate movie 6/51 (11.76 %) 6.37 s Generate movie 7/51 (13.73 %) 6.22 s Generate movie 8/51 (15.69 %) 6.05 s Generate movie 9/51 (17.65 %) 5.93 s Generate movie 10/51 (19.61 %) 5.79 s Generate movie 11/51 (21.57 %) 5.63 s Generate movie 12/51 (23.53 %) 5.55 s Generate movie 13/51 (25.49 %) 5.38 s Generate movie 14/51 (27.45 %) 5.22 s Generate movie 15/51 (29.41 %) 5.18 s Generate movie 16/51 (31.37 %) 4.99 s Generate movie 17/51 (33.33 %) 4.86 s Generate movie 18/51 (35.29 %) 4.71 s Generate movie 19/51 (37.25 %) 4.56 s Generate movie 20/51 (39.22 %) 4.44 s Generate movie 21/51 (41.18 %) 4.29 s Generate movie 22/51 (43.14 %) 4.12 s Generate movie 23/51 (45.10 %) 4.03 s Generate movie 24/51 (47.06 %) 3.86 s Generate movie 25/51 (49.02 %) 3.72 s Generate movie 26/51 (50.98 %) 3.56 s Generate movie 27/51 (52.94 %) 3.43 s Generate movie 28/51 (54.90 %) 3.28 s Generate movie 29/51 (56.86 %) 3.14 s Generate movie 30/51 (58.82 %) 2.97 s Generate movie 31/51 (60.78 %) 2.86 s Generate movie 32/51 (62.75 %) 2.71 s Generate movie 33/51 (64.71 %) 2.58 s Generate movie 34/51 (66.67 %) 2.43 s Generate movie 35/51 (68.63 %) 2.28 s Generate movie 36/51 (70.59 %) 2.13 s Generate movie 37/51 (72.55 %) 2.00 s Generate movie 38/51 (74.51 %) 1.80 s Generate movie 39/51 (76.47 %) 1.72 s Generate movie 40/51 (78.43 %) 1.57 s Generate movie 41/51 (80.39 %) 1.43 s Generate movie 42/51 (82.35 %) 1.28 s Generate movie 43/51 (84.31 %) 1.15 s Generate movie 44/51 (86.27 %) 998.23 ms Generate movie 45/51 (88.24 %) 860.91 ms Generate movie 46/51 (90.20 %) 710.71 ms Generate movie 47/51 (92.16 %) 572.73 ms Generate movie 48/51 (94.12 %) 429.82 ms Generate movie 49/51 (96.08 %) 285.30 ms Generate movie 50/51 (98.04 %) 142.56 ms Generate movie 51/51 (100.00 %) 0.00 µs | .. code-block:: Python :lineno-start: 12 from EasyFEA import Display, Folder, Models, np, Mesher, ElemType, Simulations, PyVista from EasyFEA.Geoms import Line, Domain, Point if __name__ == "__main__": Display.Clear() # ---------------------------------------------- # Configuration # ---------------------------------------------- # outputs folder = Folder.Join(Folder.RESULTS_DIR, "Thermal", "Thermal3") makeMovie = True result = "thermal" # geom R = 10 e = 2 h = 10 a = 1 # load Tmax = 5 N = 50 dt = Tmax / N # ---------------------------------------------- # Mesh # ---------------------------------------------- domain = Domain(Point(R), Point(R + e, h), e / 2) axis = Line(Point(), Point(0, 1, 0)) # Generate the mesh based on the specified dimension angle = 360 * 3 / 4 mesh = Mesher().Mesh_Revolve( domain, [], axis, angle, [angle * np.pi / 180 * R / domain.meshSize], elemType=ElemType.HEXA8, isOrganised=True, ) nodesY0 = mesh.Nodes_Conditions(lambda x, y, z: y == 0) nodesYH = mesh.Nodes_Conditions(lambda x, y, z: y == h) # ---------------------------------------------- # Simulation # ---------------------------------------------- thermalModel = Models.Thermal(k=1, c=1) simu = Simulations.ThermalSimu(mesh, thermalModel, False) simu.rho = 1 simu.add_surfLoad(nodesY0, [5], ["t"]) simu.add_surfLoad(nodesYH, [5], ["t"]) # Set the parabolic algorithm for the solver simu.Solver_Set_Parabolic_Algorithm(alpha=0.5, dt=dt) simu._Set_solutions(simu.problemType, np.ones(mesh.Nn) * -10) print() t = -dt # init time while t < Tmax: t += dt simu.Solve() simu.Save_Iter() print(f"{t:.3f} s", end="\r") # ---------------------------------------------- # Results # ---------------------------------------------- print(simu) PyVista.Plot(simu, result, plotMesh=True, nodeValues=True) if makeMovie: PyVista.Movie_simu(simu, result, folder, f"{result}.gif", plotMesh=True) .. rst-class:: sphx-glr-timing **Total running time of the script:** (0 minutes 11.610 seconds) .. _sphx_glr_download_examples_Thermal_Thermal3.py: .. only:: html .. container:: sphx-glr-footer sphx-glr-footer-example .. container:: sphx-glr-download sphx-glr-download-jupyter :download:`Download Jupyter notebook: Thermal3.ipynb ` .. container:: sphx-glr-download sphx-glr-download-python :download:`Download Python source code: Thermal3.py ` .. container:: sphx-glr-download sphx-glr-download-zip :download:`Download zipped: Thermal3.zip ` .. only:: html .. rst-class:: sphx-glr-signature `Gallery generated by Sphinx-Gallery `_