Persamaan Diferensial · Persamaan panas

Persamaan panas

Persamaan panas uₜ = uₓₓ menyatakan bahwa laju perubahan suhu di setiap titik ditentukan oleh lengkungan profil suhu di sana. Pada batang, persamaan ini menghaluskan semua tonjolan; tonjolan dengan skala terkecil hilang lebih dahulu.

Persamaan panas dalam visualisasi ini

Grafik atas menunjukkan suhu uu sepanjang batang dari 0 hingga π\pi: samar pada awal dan kuning pada waktu tt. Panah biru menunjukkan seberapa cepat setiap titik menghangat atau mendingin. Batang di bawahnya diwarnai menurut suhu, dan baris terbawah menunjukkan ukuran setiap mode sinus pada skala logaritmik, atau mode kosinus saat ujung terisolasi.

Laju perubahan terhadap waktu adalah turunan kedua sepanjang batang
ut=uxxu_t = u_{xx}

Mengapa suatu titik berubah paling cepat di tempat kurva paling melengkung? Untuk langkah kecil hh, uxx(x)≈(u(x−h)+u(x+h)−2u(x))/h2u_{xx}(x) \approx \bigl(u(x-h) + u(x+h) - 2u(x)\bigr)/h^2: yakni 2/h22/h^2 kali selisih antara rata-rata kedua tetangga dan u(x)u(x). Saat uxx>0u_{xx} > 0, kurva melengkung ke atas; titik itu berada di bawah rata-rata tetangga dekatnya dan menghangat. Saat uxx<0u_{xx} < 0, titik itu mendingin. Pada bagian lurus, ut=0u_t = 0, sepanas apa pun bagian itu. Probe dapat menunjukkannya dengan dua tetangga dekat dan rata-ratanya. Lompatan dan sudut langsung menjadi halus: untuk setiap t>0t > 0, suhu membentuk kurva halus.

Mengapa gelombang kecil hilang lebih dahulu? Satu mode, sin⁡kx\sin kx saat kedua ujung dijaga pada 0 atau cos⁡kx\cos kx saat ujung terisolasi, memenuhi uxx=−k2uu_{xx} = -k^2 u. Jadi bentuknya tetap dan hanya menyusut dengan faktor e−k2te^{-k^2 t}. Gelombang tiga kali lebih rapat memiliki turunan kedua sembilan kali lebih besar pada tinggi yang sama dan meluruh sembilan kali lebih cepat. Setiap suhu awal adalah jumlah mode-mode ini, sehingga batang di baris bawah turun mulai dari kanan: mode kk turun di bawah 10% nilai awal setelah tk=ln⁡10/k2t_k = \ln 10/k^2. Pada jam logaritmik, setiap kenaikan waktu sepuluh kali lipat menghilangkan mode sekitar 10≈3.2\sqrt{10} \approx 3.2 kali lebih lebar. Penggaris merah 2t2\sqrt t menunjukkan skala itu: mode tunggal dengan setengah panjang gelombang tidak melebihi penggaris menyimpan paling banyak e−π2/4≈8.48%e^{-\pi^2/4} \approx 8.48\% amplitudonya.

Apa yang diubah oleh kedua ujung? Saat dijaga pada 0, ujung-ujung itu mewakili es, dan uu diukur relatif terhadap suhu es. Seiring waktu, mode terendah yang ada dalam kondisi awal mendominasi bentuk. Jika suhu awal tidak pernah negatif dan tidak semuanya nol, panas keluar lewat ujung, ∫u dx\int u\,dx menurun, dan profil menjadi satu lengkungan sinus yang turun seperti e−te^{-t}. Tidak ada panas yang melewati ujung terisolasi, sehingga ∫u dx\int u\,dx tetap. Mode-mode itu kini berupa kosinus, dan mode k=0k = 0, yaitu suhu rata-rata, tidak pernah meluruh: seluruh batang mendekati rata-ratanya. Kedua jenis kondisi ujung bekerja sejak t=0t = 0. Profil awal yang tidak cocok di ujung, seperti batang bersuhu seragam di antara es, langsung menyesuaikan di sana.

Hal yang perlu diperhatikan

  • Semua panah memakai skala yang sama, ditentukan oleh titik yang berubah paling cepat pada saat itu. Legenda menunjukkan laju untuk panah acuan. Bandingkan panjang panah pada saat yang sama, bukan antarsaat.
  • Suhu di bawah suhu es bernilai negatif. Jika kedua ujung dijaga pada 0 dan suhu awal memiliki kedua tanda, yang ditampilkan adalah nilai bersih ∫u dx\int u\,dx, bukan proporsi.

Kontrol

  • Difusikan meninggalkan cuplikan samar setiap kali waktu naik sepuluh kali lipat.
  • Klik batang untuk memindahkan probe ke sana.
  • Pilih Gambar lalu seret pada grafik atas untuk menentukan suhu awal; menggambar mengembalikan tt ke 0.

Terkait