المعادلات التفاضلية · معادلة الحرارة

معادلة الحرارة

تقول معادلة الحرارة ⁦uₜ = uₓₓ⁩ إن معدل تغير درجة الحرارة عند كل نقطة تحدده كيفية انحناء منحنى الحرارة هناك. وعلى القضيب، تنعّم كل نتوء، فتزول أدق النتوءات أولًا.

معادلة الحرارة في هذا التصوّر

يعرض الرسم العلوي درجة الحرارة uu على طول قضيب من الصفر إلى π\pi: بلون باهت في البداية وبالأصفر عند الزمن tt. تبيّن الأسهم الزرقاء سرعة ارتفاع درجة الحرارة أو انخفاضها عند كل نقطة. يلوّن الشريط أسفله القضيب حسب الحرارة، ويعرض الصف السفلي حجم كل نمط جيبي على مقياس لوغاريتمي، أو أنماط جيب التمام عند عزل الطرفين.

معدل التغير مع الزمن هو المشتقة الثانية على طول القضيب
ut=uxxu_t = u_{xx}

لماذا تتغير النقطة أسرع ما يكون حيث يشتد انحناء المنحنى؟ لخطوة صغيرة 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: أي الفرق بين متوسط الجارين وu(x)u(x) مضروبًا في 2/h22/h^2. حيث uxx>0u_{xx} > 0، ينحني المنحنى إلى أعلى؛ فتكون النقطة دون متوسط جيرانها القريبين وترتفع حرارتها. وحيث uxx<0u_{xx} < 0، تنخفض حرارتها. وعلى مقطع مستقيم، ut=0u_t = 0 مهما كان ساخنًا. يمكن للمسبار إظهار ذلك بجارين قريبين ومتوسطهما. تصبح القفزات والزوايا ملساء فورًا: لكل t>0t > 0، تكون درجة الحرارة منحنى أملس.

لماذا تزول التموجات الدقيقة أولًا؟ النمط المنفرد، sin⁡kx\sin kx عند تثبيت الطرفين على الصفر أو cos⁡kx\cos kx عند عزلهما، يحقق uxx=−k2uu_{xx} = -k^2 u؛ لذا يحتفظ بشكله ولا يفعل سوى التناقص بالعامل e−k2te^{-k^2 t}. ثلاثة أضعاف التموجات تعطي مشتقة ثانية أكبر تسع مرات عند الارتفاع نفسه، واضمحلالًا أسرع تسع مرات. كل درجة حرارة ابتدائية هي مجموع هذه الأنماط، لذلك تهبط الأعمدة في الصف السفلي بدءًا من اليمين: ينخفض النمط kk إلى أقل من ⁦10%⁩ من بدايته بعد tk=ln⁡10/k2t_k = \ln 10/k^2. في الساعة اللوغاريتمية، كل زيادة في الزمن بعامل عشرة تمحو أنماطًا أعرض بنحو 10≈3.2\sqrt{10} \approx 3.2 مرات. تعطي المسطرة الحمراء 2t2\sqrt t هذا المقياس: النمط المنفرد الذي لا يتجاوز نصف طوله الموجي طول المسطرة يحتفظ على الأكثر بـe−π2/4≈8.48%e^{-\pi^2/4} \approx 8.48\% من سعته.

ما الذي يغيره الطرفان؟ حين يثبتان عند الصفر، يمثلان الجليد، وتقاس uu بالنسبة إلى درجة حرارة الجليد. مع الزمن، يهيمن أدنى نمط موجود في البداية على الشكل. إذا كانت درجة الحرارة الابتدائية غير سالبة في أي موضع وليست صفرًا في كل المواضع، تخرج الحرارة عبر الطرفين، وتنخفض ∫u dx\int u\,dx، ويصبح الشكل قوسًا جيبيًا واحدًا يهبط وفق e−te^{-t}. لا تعبر أي حرارة الطرفين المعزولين، لذا تبقى ∫u dx\int u\,dx ثابتة. تصبح الأنماط حينئذ جيوب تمام، ولا يضمحل النمط k=0k = 0، أي متوسط درجة الحرارة: يقترب القضيب كله من متوسطه. يعمل كلا الشرطين منذ t=0t = 0؛ وأي حالة ابتدائية لا تطابق الشرط عند الطرف، كقضيب منتظم الحرارة بين قطعتي جليد، تتكيف هناك فورًا.

ما الذي تلاحظه

  • تشترك الأسهم في مقياس واحد تحدده أسرع نقطة تغيرًا في تلك اللحظة، ويعطي المفتاح معدل السهم المرجعي. قارن أطوال الأسهم في اللحظة نفسها، لا بين أزمنة مختلفة.
  • درجات الحرارة الأقل من حرارة الجليد سالبة. عند تثبيت الطرفين على الصفر ووجود إشارتين في الحالة الابتدائية، يُعرض الصافي ∫u dx\int u\,dx بدلًا من نسبة.

أدوات التحكم

  • يترك «انتشار» لقطة باهتة كلما زاد الزمن بعامل عشرة.
  • انقر على القضيب لنقل المسبار إليه.
  • اختر «ارسم» واسحب على الرسم العلوي لتحديد درجة الحرارة الابتدائية؛ يعيد الرسم tt إلى الصفر.

مواضيع ذات صلة