.. 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 : scipy ============= Boundary Conditions ============= Unspecified. =================== Results =================== Unspecified. =================== TicTac =================== Mesh : 94.740 ms Boundary Conditions : 37.670 µs Matrix : 36.537 ms Solver : 1.247 s PyVista_Interface : 12.044 s Generate movie 01/52 (1.92 %) 7.04 s Generate movie 02/52 (3.85 %) 6.40 s Generate movie 03/52 (5.77 %) 6.28 s Generate movie 04/52 (7.69 %) 6.19 s Generate movie 05/52 (9.62 %) 6.01 s Generate movie 06/52 (11.54 %) 6.02 s Generate movie 07/52 (13.46 %) 5.74 s Generate movie 08/52 (15.38 %) 5.72 s Generate movie 09/52 (17.31 %) 5.59 s Generate movie 10/52 (19.23 %) 5.58 s Generate movie 11/52 (21.15 %) 5.30 s Generate movie 12/52 (23.08 %) 5.13 s Generate movie 13/52 (25.00 %) 5.08 s Generate movie 14/52 (26.92 %) 4.96 s Generate movie 15/52 (28.85 %) 4.81 s Generate movie 16/52 (30.77 %) 4.72 s Generate movie 17/52 (32.69 %) 4.61 s Generate movie 18/52 (34.62 %) 4.45 s Generate movie 19/52 (36.54 %) 4.31 s Generate movie 20/52 (38.46 %) 4.23 s Generate movie 21/52 (40.38 %) 4.17 s Generate movie 22/52 (42.31 %) 3.97 s Generate movie 23/52 (44.23 %) 3.84 s Generate movie 24/52 (46.15 %) 3.74 s Generate movie 25/52 (48.08 %) 3.60 s Generate movie 26/52 (50.00 %) 3.40 s Generate movie 27/52 (51.92 %) 3.31 s Generate movie 28/52 (53.85 %) 3.18 s Generate movie 29/52 (55.77 %) 3.01 s Generate movie 30/52 (57.69 %) 3.05 s Generate movie 31/52 (59.62 %) 2.75 s Generate movie 32/52 (61.54 %) 2.62 s Generate movie 33/52 (63.46 %) 2.48 s Generate movie 34/52 (65.38 %) 2.36 s Generate movie 35/52 (67.31 %) 2.23 s Generate movie 36/52 (69.23 %) 2.14 s Generate movie 37/52 (71.15 %) 2.00 s Generate movie 38/52 (73.08 %) 1.86 s Generate movie 39/52 (75.00 %) 1.75 s Generate movie 40/52 (76.92 %) 1.60 s Generate movie 41/52 (78.85 %) 1.44 s Generate movie 42/52 (80.77 %) 1.32 s Generate movie 43/52 (82.69 %) 1.18 s Generate movie 44/52 (84.62 %) 1.04 s Generate movie 45/52 (86.54 %) 908.28 ms Generate movie 46/52 (88.46 %) 783.62 ms Generate movie 47/52 (90.38 %) 655.70 ms Generate movie 48/52 (92.31 %) 522.28 ms Generate movie 49/52 (94.23 %) 390.51 ms Generate movie 50/52 (96.15 %) 264.28 ms Generate movie 51/52 (98.08 %) 133.43 ms Generate movie 52/52 (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.Results_Dir() 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.Thermal(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.050 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 `_