Setelah melakukan simulasi untuk perpindahan panas pada aplikasi CFDSOF, saya mendapatkan persamaan Curve Fitting yaitu T(x)=−400×4+800×3−480×2+37.78x+303 dengan grafik seperti berikut
Langkah selanjutnya adalah untuk menghitung flux panas dan daya. Untuk hal ini saya menggunakan Hukum Fourier dan mendapatkan persamaan sebagai berikut, q(x)=1600×3−2400×2+960x−37.78. Persamaan tersebut dapat divisualisasikan sebagai:
Selanjutnya, saya akan mencari kontur untuk batang 2D tersebut. Saya mennggunakan bantuan Grok AI serta Google Collab untuk mendapatkan visualisasi tersebut. Saya mendapatkan pseudocode sebagai berikut:
import numpy as np
import matplotlib.pyplot as plt
x = np.array([0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1])
y = np.array([0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1.0])
T = np.array([
[303, 337.78, 351.28, 357.24, 359.91, 360.68, 359.91, 357.24, 351.28, 337.78, 303],
[303, 323.82, 337.11, 344.78, 348.71, 349.92, 348.71, 344.77, 337.11, 323.82, 303],
[303, 317.42, 328.57, 336.04, 340.24, 341.59, 340.24, 336.04, 328.57, 317.42, 303],
[303, 314.27, 323.73, 330.56, 334.61, 335.95, 334.61, 330.56, 323.73, 314.27, 303],
[303, 312.94, 321.50, 327.86, 331.71, 332.99, 317.13, 327.86, 321.50, 312.94, 303],
[303, 313, 321.47, 327.67, 331.37, 332, 331.37, 327.67, 321.47, 313, 303],
[303, 314.60, 323.73, 329.95, 333.52, 334.67, 333.52, 329.95, 323.73, 314.60, 303],
[303, 318.67, 328.88, 334.91, 338.07, 339.06, 338.07, 334.91, 328.88, 318.67, 303],
[303, 328.22, 338.20, 342.73, 344.80, 345.42, 344.80, 342.73, 338.20, 328.22, 303]
])
X, Y = np.meshgrid(x, y)
plt.figure(figsize=(8, 6))
contour = plt.contourf(X, Y, T, cmap=’hot’, levels=20) # Filled contour
plt.colorbar(contour, label=’Temperatur (Kelvin)’)
plt.contour(X, Y, T, colors=’black’, levels=10, linewidths=0.5) # Garis kontur
plt.xlabel(‘x (meter)’)
plt.ylabel(‘y (meter)’)
plt.title(‘Distribusi Temperatur 2D (Kontur)’)
plt.grid(True)
plt.show()
Pseudocode tersebut akan menghasilkan visualisasi sebagai berikut:
