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: 39.057 ms
Boundary Conditions: 9.537 µs
Matrix: 15.711 ms
Solver: 491.168 ms

Generate movie 01/51 (1.96 %) 5.63 s
Generate movie 02/51 (3.92 %) 4.68 s
Generate movie 03/51 (5.88 %) 4.54 s
Generate movie 04/51 (7.84 %) 4.46 s
Generate movie 05/51 (9.80 %) 4.36 s
Generate movie 06/51 (11.76 %) 4.41 s
Generate movie 07/51 (13.73 %) 4.22 s
Generate movie 08/51 (15.69 %) 4.12 s
Generate movie 09/51 (17.65 %) 3.99 s
Generate movie 10/51 (19.61 %) 3.99 s
Generate movie 11/51 (21.57 %) 3.88 s
Generate movie 12/51 (23.53 %) 3.95 s
Generate movie 13/51 (25.49 %) 3.66 s
Generate movie 14/51 (27.45 %) 3.63 s
Generate movie 15/51 (29.41 %) 3.58 s
Generate movie 16/51 (31.37 %) 3.33 s
Generate movie 17/51 (33.33 %) 3.28 s
Generate movie 18/51 (35.29 %) 3.38 s
Generate movie 19/51 (37.25 %) 3.17 s
Generate movie 20/51 (39.22 %) 2.97 s
Generate movie 21/51 (41.18 %) 2.93 s
Generate movie 22/51 (43.14 %) 2.93 s
Generate movie 23/51 (45.10 %) 2.85 s
Generate movie 24/51 (47.06 %) 2.61 s
Generate movie 25/51 (49.02 %) 2.50 s
Generate movie 26/51 (50.98 %) 2.44 s
Generate movie 27/51 (52.94 %) 2.34 s
Generate movie 28/51 (54.90 %) 2.30 s
Generate movie 29/51 (56.86 %) 2.15 s
Generate movie 30/51 (58.82 %) 2.04 s
Generate movie 31/51 (60.78 %) 2.00 s
Generate movie 32/51 (62.75 %) 1.87 s
Generate movie 33/51 (64.71 %) 1.78 s
Generate movie 34/51 (66.67 %) 1.71 s
Generate movie 35/51 (68.63 %) 1.57 s
Generate movie 36/51 (70.59 %) 1.45 s
Generate movie 37/51 (72.55 %) 1.49 s
Generate movie 38/51 (74.51 %) 1.28 s
Generate movie 39/51 (76.47 %) 1.17 s
Generate movie 40/51 (78.43 %) 1.06 s
Generate movie 41/51 (80.39 %) 976.42 ms
Generate movie 42/51 (82.35 %) 870.90 ms
Generate movie 43/51 (84.31 %) 778.23 ms
Generate movie 44/51 (86.27 %) 693.15 ms
Generate movie 45/51 (88.24 %) 597.38 ms
Generate movie 46/51 (90.20 %) 501.81 ms
Generate movie 47/51 (92.16 %) 388.55 ms
Generate movie 48/51 (94.12 %) 296.01 ms
Generate movie 49/51 (96.08 %) 199.47 ms
Generate movie 50/51 (98.04 %) 99.68 ms
Generate movie 51/51 (100.00 %) 0.00 µs

from EasyFEA import Display, Folder, Models, 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.Results_Dir() + f"{dim}D"
    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, meshSize=a / 4)

    if dim == 2:
        mesh = domain.Mesh_2D([circle], ElemType.TRI6)
    else:
        mesh = domain.Mesh_Extrude([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, direction=[0, 0, -1])

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

    # ----------------------------------------------
    # Simulation
    # ----------------------------------------------
    thermalModel = Models.Thermal(k=1, c=1, thickness=1)
    simu = Simulations.Thermal(mesh, thermalModel)
    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)