• matematika
• metodenumerik
• numerik
• pemrograman
• tutorial

Metode Numerik Sederhana untuk Sebaran Panas pada Mesin Kopi

Simulasi termal atau thermal simulation merupakan teknik numerik yang digunakan untuk memprediksi dan mengoptimalkan sebaran panas dalam satu sistem. Metode ini dapat digunakan dalam teknologi proses manufaktur, pemrosesan material, dan produksi energi.

Dalam industri kopi, simulasi termal dapat digunakan untuk memprediksi suhu dalam proses pembuatan kopi menggunakan mesin espresso atau peralatan kopi lainnya. Simulasi termal membantu perancang dan produsen kopi mengoptimalkan tekanan, suhu, dan waktu untuk menghasilkan kopi yang baik.

Simulasi termal juga dapat membantu dalam mengurangi waktu dan biaya pengujian fisik, sehingga mempercepat proses pengembangan produk dan meningkatkan efisiensi produksi.

Artikel ini membahas penerapan sederhana thermal simulation dari model matematika sampai ke kode program dalam bahasa C.

model matematika

Persamaan di bawah ini menjelaskan hubungan perubahan suhu sistem seturut waktu di bawah pengaruh koefisien difusivitas termal bahan yang digunakan.

Model ini dapat dimodifikasi dengan menambahkan kondisi batas atau parameter tambahan untuk memperkirakan suhu dalam situasi yang lebih kompleks.

Dalam simulasi termal, persamaan ini kemudian akan dipecahkan menggunakan metode numerik seperti metode beda hingga atau metode elemen hingga untuk menghasilkan distribusi suhu di dalam sistem.

∂T/∂t = α∇²T

Keterangan: T : suhu dalam sistem t : waktu α : koefisien difusivitas termal ∇² : operator Laplace

algoritma

1. Inisialisasi parameter simulasi termal a. Tentukan dimensi dan ukuran sistem termal b. Tentukan parameter termal seperti konduktivitas termal dan kapasitas panas c. Tentukan kondisi batas, seperti suhu atau kalor yang ditambahkan atau dihilangkan pada sistem termal d. Tentukan kondisi awal sistem termal, seperti suhu awal atau distribusi suhu awal di dalam sistem

2. Tentukan jangka waktu simulasi dan langkah waktu Δt a. Tentukan waktu total simulasi T b. Tentukan jumlah langkah waktu N = T/Δt

3. Tentukan nilai α a. Tentukan nilai α yang sesuai untuk bahan yang membentuk sistem termal

4. Inisialisasi nilai suhu awal untuk setiap titik di dalam sistem termal a. Tentukan suhu awal pada setiap titik dalam sistem termal

5. Untuk setiap langkah waktu dari 1 sampai N, hitung suhu pada setiap titik di dalam sistem termal a. Hitung laju perpindahan panas di setiap titik dalam sistem termal menggunakan persamaan dasar fisika yang digunakan dalam model matematika b. Hitung suhu pada setiap titik di dalam sistem termal menggunakan metode numerik seperti metode beda hingga atau metode elemen hingga c. Simpan nilai suhu pada setiap titik pada setiap langkah waktu

6. Analisis hasil simulasi termal, dan gunakan data yang dihasilkan untuk memperbaiki dan mengoptimalkan sistem termal a. Analisis hasil simulasi termal, seperti distribusi suhu, waktu respon sistem, atau laju perpindahan panas b. Gunakan data yang dihasilkan untuk memperbaiki sistem termal dan mengoptimalkan performa sistem termal pada situasi yang berbeda

7. Selesai

kode C

#include  <  stdio.h >
#include  <  math.h >

// Define parameter simulasi termal
#define DIMENSION 2 // Dimensi sistem
#define SIZE_X 10 // Ukuran sistem pada sumbu X
#define SIZE_Y 10 // Ukuran sistem pada sumbu Y
#define CONDUCTIVITY 1.0 // Konduktivitas termal
#define HEAT_CAPACITY 1.0 // Kapasitas panas
#define BOUNDARY_CONDITION 0.0 // Kondisi batas suhu
#define TIME_STEPS 100 // Jumlah langkah waktu
#define TIME_STEP_SIZE 0.01 // Ukuran langkah waktu

// Persamaan dasar
double alpha = CONDUCTIVITY / HEAT_CAPACITY;

double temperature[SIZE_X][SIZE_Y];
double temperature_new[SIZE_X][SIZE_Y];

// Inisialisasi suhu awal pada setiap titik di dalam sistem termal
void initialize_temperature() {
    int i, j;
    for (i = 0; i  <  SIZE_X; i++) {
        for (j = 0; j  <  SIZE_Y; j++) {
            temperature[i][j] = BOUNDARY_CONDITION;
        }
    }
}

// Hitung suhu pada setiap titik di dalam sistem termal menggunakan metode beda hingga
void calculate_temperature() {
    int i, j, n;
    double laplacian;

    // Iterasi untuk setiap langkah waktu
    for (n = 0; n  <  TIME_STEPS; n++) {
        // Iterasi untuk setiap titik dalam sistem
        for (i = 1; i  <  SIZE_X - 1; i++) {
            for (j = 1; j  <  SIZE_Y - 1; j++) {
                // Hitung titik laplacian
                laplacian = temperature[i+1][j] + temperature[i-1][j] + temperature[i][j+1] + temperature[i][j-1] - 4*temperature[i][j];
                // Hitung suhu baru
                temperature_new[i][j] = temperature[i][j] + alpha * laplacian * TIME_STEP_SIZE;
            }
        }

        // Update nilai suhu
        for (i = 1; i  <  SIZE_X - 1; i++) {
            for (j = 1; j  <  SIZE_Y - 1; j++) {
                temperature[i][j] = temperature_new[i][j];
            }
        }
    }
}

int main() {
    initialize_temperature();
    calculate_temperature();

    return 0;
}