Thermal2

Transient thermal simulation.

Static Scene

Interactive Scene

Volume: 0.913
0.000 s
0.010 s
0.020 s
0.030 s
0.040 s
0.050 s
0.060 s
0.070 s
0.080 s
0.090 s
0.100 s
0.110 s
0.120 s
0.130 s
0.140 s
0.150 s
0.160 s
0.170 s
0.180 s
0.190 s
0.200 s
0.210 s
0.220 s
0.230 s
0.240 s
0.250 s
0.260 s
0.270 s
0.280 s
0.290 s
0.300 s
0.310 s
0.320 s
0.330 s
0.340 s
0.350 s
0.360 s
0.370 s
0.380 s
0.390 s
0.400 s
0.410 s
0.420 s
0.430 s
0.440 s
0.450 s
0.460 s
0.470 s
0.480 s
0.490 s
0.500 s

==================== Mesh ====================

Element type: PRISM18
Ne = 288, Nn = 1512

==================== Model ====================


Thermal :
thermal conductivity (k)  : 1
thermal mass capacity (c) : 1

solver : scipy

============= Boundary Conditions =============

Unspecified.

=================== Results ===================

Unspecified.

=================== TicTac ===================

Mesh : 45.364 ms
Boundary Conditions : 17.643 µs
Matrix : 18.022 ms
Solver : 1.201 s
PyVista_Interface : 35.707 ms

Generate movie 01/51 (1.96 %) 8.12 s
Generate movie 02/51 (3.92 %) 7.73 s
Generate movie 03/51 (5.88 %) 7.44 s
Generate movie 04/51 (7.84 %) 7.31 s
Generate movie 05/51 (9.80 %) 7.13 s
Generate movie 06/51 (11.76 %) 7.05 s
Generate movie 07/51 (13.73 %) 6.87 s
Generate movie 08/51 (15.69 %) 6.80 s
Generate movie 09/51 (17.65 %) 6.54 s
Generate movie 10/51 (19.61 %) 6.40 s
Generate movie 11/51 (21.57 %) 6.27 s
Generate movie 12/51 (23.53 %) 6.09 s
Generate movie 13/51 (25.49 %) 6.01 s
Generate movie 14/51 (27.45 %) 5.81 s
Generate movie 15/51 (29.41 %) 5.68 s
Generate movie 16/51 (31.37 %) 5.45 s
Generate movie 17/51 (33.33 %) 5.32 s
Generate movie 18/51 (35.29 %) 5.21 s
Generate movie 19/51 (37.25 %) 4.99 s
Generate movie 20/51 (39.22 %) 4.93 s
Generate movie 21/51 (41.18 %) 4.81 s
Generate movie 22/51 (43.14 %) 4.62 s
Generate movie 23/51 (45.10 %) 4.45 s
Generate movie 24/51 (47.06 %) 4.30 s
Generate movie 25/51 (49.02 %) 4.13 s
Generate movie 26/51 (50.98 %) 3.92 s
Generate movie 27/51 (52.94 %) 3.80 s
Generate movie 28/51 (54.90 %) 3.65 s
Generate movie 29/51 (56.86 %) 3.47 s
Generate movie 30/51 (58.82 %) 3.31 s
Generate movie 31/51 (60.78 %) 3.17 s
Generate movie 32/51 (62.75 %) 3.02 s
Generate movie 33/51 (64.71 %) 2.83 s
Generate movie 34/51 (66.67 %) 2.70 s
Generate movie 35/51 (68.63 %) 2.53 s
Generate movie 36/51 (70.59 %) 2.41 s
Generate movie 37/51 (72.55 %) 2.24 s
Generate movie 38/51 (74.51 %) 2.08 s
Generate movie 39/51 (76.47 %) 1.93 s
Generate movie 40/51 (78.43 %) 1.76 s
Generate movie 41/51 (80.39 %) 1.59 s
Generate movie 42/51 (82.35 %) 1.44 s
Generate movie 43/51 (84.31 %) 1.27 s
Generate movie 44/51 (86.27 %) 1.12 s
Generate movie 45/51 (88.24 %) 957.74 ms
Generate movie 46/51 (90.20 %) 809.88 ms
Generate movie 47/51 (92.16 %) 638.07 ms
Generate movie 48/51 (94.12 %) 473.71 ms
Generate movie 49/51 (96.08 %) 322.38 ms
Generate movie 50/51 (98.04 %) 160.30 ms
Generate movie 51/51 (100.00 %) 0.00 µs

from EasyFEA import Display, Folder, Models, Mesher, ElemType, Simulations, PyVista
from EasyFEA.Geoms import Circle, Domain, Point

if __name__ == "__main__":
    Display.Clear()

    # ----------------------------------------------
    # Configuration
    # ----------------------------------------------
    dim = 3  # Set the simulation dimension (2D or 3D)

    # outputs
    folder = Folder.Join(Folder.RESULTS_DIR, "Thermal", "Thermal2")
    makeMovie = True
    result = "thermal"

    # geom
    a = 1

    # load
    Tmax = 0.5
    N = 50
    dt = Tmax / N

    # ----------------------------------------------
    # Mesh
    # ----------------------------------------------
    domain = Domain(Point(), Point(a, a), a / 4)
    circle = Circle(Point(a / 2, a / 2), diam=a / 3, isHollow=True, meshSize=a / 4)

    if dim == 2:
        mesh = Mesher().Mesh_2D(domain, [circle], ElemType.TRI6)
    else:
        mesh = Mesher().Mesh_Extrude(
            domain, [circle], [0, 0, -a], [4], ElemType.PRISM18
        )

    nodesX0 = mesh.Nodes_Conditions(lambda x, y, z: x == 0)
    nodesXa = mesh.Nodes_Conditions(lambda x, y, z: x == a)
    nodesCircle = mesh.Nodes_Cylinder(circle, [0, 0, -1])

    if dim == 3:
        print(f"Volume: {mesh.volume:.3}")

    # ----------------------------------------------
    # Simulation
    # ----------------------------------------------
    thermalModel = Models.Thermal(k=1, c=1, thickness=1)
    simu = Simulations.ThermalSimu(mesh, thermalModel, False)
    simu.rho = 1

    simu.add_dirichlet(nodesX0, [0], ["t"])
    simu.add_dirichlet(nodesXa, [40], ["t"])

    # Set the parabolic algorithm for the solver
    simu.Solver_Set_Parabolic_Algorithm(alpha=0.5, dt=dt)

    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).show()

    if makeMovie:
        PyVista.Movie_simu(simu, result, folder, f"{result}.gif", plotMesh=True)