Persamaan diferensial partial domestial (PDEs) adalah fundamental dalam pemodelan berbagai fenomena fisik, seperti transfer panas, propagasi gelombang, dan dinamika fluid. Menggunakan pustaka Python seperti NumPy dan SciPy menyederhanakan proses pemecahan persamaan ini secara numerik. Artikel ini menyediakan pendekatan langkah- demi langkah untuk menyelesaikan PDE dengan alat-alat ini.

Bezakan Masalah

Sebagai contoh, pertimbangkan persamaan panas satu dimensi:

⁇ u/ ⁇ t = ⁇ 2u/ ⁇ x2

Nyatakan kondisi awal, kondisi batas, dan parameter seperti diffusivity termal α.

Memdiskrediksi Domain

Sebagai contoh, membuat sebuah kisi poin:

import numpy as np

[[CALAT:0]]x = np.linspace(0, 1, 100)

Mengeset langkah waktu dan total waktu simulasi untuk mengontrol stabilitas numerik dan akurasi.

Mengimplementasi Metode Numerikal

Untuk persamaan panas, skema eksplisit update suhu pada setiap titik:

[[Oble]u new = u + r * (u[i+1] - 2*u[i] + u[i-1])[

Parameter isbn= yang tidak diketahui akan diabaikan (bantuan) adalah sebuah parameter fixability yang dihitung sebagai ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇

Memburu Simulasi

Initerate dari langkah waktu, memperbarui array solusi pada setiap iterasi. Gunakan operasi NumPy untuk efisiensi:

for n dalam prege(langkah):

  u[1:-1] = u[1:-1] + r * (u[2] - 2*u[1:-1] + u[:-2])

Terapkan kondisi batas setelah setiap pembaruan untuk mempertahankan kendala masalah.

Hasil Visualisasi

Plot keadaan awal dan akhir untuk mengamati proses difusi.

Kode Contoh sebol:

[[ZOLT:0]]impor mattplotlib.pyplot sebagai plt[

[[GALAL:0]]plt.plot(x, u)[