المعادلات الأولية

المعادلات الأولية هي مجموعة من المعادلات التفاضلية الجزئية غير الخطية التي تُستخدم لتقريب تدفق الغلاف الجوي العالمي ، وتُستخدم في معظم نماذج الغلاف الجوي . وهي تتكون من ثلاث مجموعات رئيسية من معادلات التوازن:

  1. معادلة الاستمرارية : تمثل حفظ الكتلة.
  2. قانون حفظ الزخم : يتألف من صيغة لمعادلات نافيير-ستوكس التي تصف التدفق الهيدروديناميكي على سطح كرة بافتراض أن الحركة الرأسية أصغر بكثير من الحركة الأفقية (التوازن الهيدروستاتيكي) وأن عمق طبقة السائل صغير مقارنة بنصف قطر الكرة.
  3. معادلة الطاقة الحرارية : تربط درجة الحرارة الكلية للنظام بمصادر الحرارة ومصارفها

يمكن تحويل المعادلات الأولية إلى معادلات لابلاس المدية ، وهي مسألة القيم الذاتية التي يمكن من خلالها تحديد الحل التحليلي للبنية العرضية للتدفق.

بشكل عام، ترتبط جميع أشكال المعادلات الأولية تقريبًا بالمتغيرات الخمسة u و v و ω و T و W وتطورها عبر المكان والزمان.

تمت كتابة المعادلات لأول مرة بواسطة فيلهلم بيركنز . [ 1 ]

التعريفات

  • u{\displaystyle u}هي السرعة المدارية (السرعة في الاتجاه الشرقي الغربي المماسية للكرة)
  • v{\displaystyle v}هي السرعة الزوالية (السرعة في اتجاه الشمال والجنوب مماسية للكرة)
  • ω{\displaystyle \omega }هي السرعة الرأسية في إحداثيات متساوية الضغط
  • تي{\displaystyle T}درجة الحرارة
  • Φ{\displaystyle \Phi }هو الجهد الجيوديسي
  • و{\displaystyle f}يمثل هذا المصطلح قوة كوريوليس ، وهو يساوي2Ωالخطيئة(ϕ){\displaystyle 2\Omega \sin(\phi )}، أينΩأوميغامعدل الدوران الزاوي للأرض (2π/24{\displaystyle 2\pi /24}راديان لكل ساعة فلكية)، وϕ{\displaystyle \phi }هو خط العرض
  • R{\displaystyle R}ثابت الغاز
  • ص{\displaystyle p}الضغط
  • ρ{\displaystyle \rho }الكثافة
  • جص{\displaystyle c_{p}}هي الحرارة النوعية على سطح ذي ضغط ثابت
  • ج{\displaystyle J}هو معدل تدفق الحرارة لكل وحدة زمنية لكل وحدة كتلة
  • دبليو{\displaystyle W}هو الماء القابل للتكثيف
  • Π{\displaystyle \Pi }هي دالة إكسنر
  • θ{\displaystyle \theta }هي درجة الحرارة الكامنة
  • η{\displaystyle \eta }الدوامة المطلقة

القوى التي تسبب حركة الغلاف الجوي

تشمل القوى التي تسبب حركة الغلاف الجوي قوة تدرج الضغط ، والجاذبية ، والاحتكاك اللزج . وتجتمع هذه القوى معًا لتشكل القوى التي تسرع غلافنا الجوي.

تتسبب قوة تدرج الضغط في تسارع يدفع الهواء من مناطق الضغط العالي إلى مناطق الضغط المنخفض. ويمكن التعبير عن ذلك رياضياً كما يلي:

وم=1ρدصدx.{\displaystyle {\frac {f}{m}}={\frac {1}{\rho }}{\frac {dp}{dx}}.}

تعمل قوة الجاذبية على تسريع الأجسام بمعدل 9.8  م/ث² تقريبًا باتجاه مركز الأرض.

يمكن تقريب القوة الناتجة عن الاحتكاك اللزج على النحو التالي:

ور=وأ1ρμ((μv)+(λv)).{\displaystyle f_{r}={f \over a}{1 \over \rho }\mu \left(\nabla \cdot (\mu \nabla v)+\nabla (\lambda \nabla \cdot v)\right).}

باستخدام قانون نيوتن الثاني، يمكن جمع هذه القوى (المشار إليها في المعادلات أعلاه بالتسارعات الناتجة عنها) للحصول على معادلة حركة تصف هذا النظام. ويمكن كتابة هذه المعادلة على النحو التالي:

دvدت=-(1ρ)ص-ز(رر)+ور{\displaystyle {\frac {dv}{dt}}=-({\frac {1}{\rho }})\nabla pg({\frac {r}{r}})+f_{r}}
ز=زهـ.{\displaystyle g=g_{e}.\,}

لذلك، لإكمال نظام المعادلات والحصول على 6 معادلات و6 متغيرات:

  • دvدت=-(1ρ)ص-ز(رر)+(1ρ)[(μv)+(λv)]{\displaystyle {\frac {dv}{dt}}=-({\frac {1}{\rho }})\nabla pg({\frac {r}{r}})+({\frac {1}{\rho }})\left[\nabla \cdot (\mu \nabla v)+\nabla (\lambda \nabla \cdot v)\right]}
  • جvدتيدت+صدαدت=q+و{\displaystyle c_{v}{\frac {dT}{dt}}+p{\frac {d\alpha }{dt}}=q+f}
  • دρدت+ρv=0{\displaystyle {\frac {d\rho }{dt}}+\rho \nabla \cdot v=0}
  • ص=نتي.{\displaystyle p=nT.}

حيث n هي الكثافة العددية بالمول، و T:=RT هي القيمة المكافئة لدرجة الحرارة بالجول/مول.

أشكال المعادلات الأولية

يعتمد الشكل الدقيق للمعادلات الأولية على نظام الإحداثيات الرأسية المُختار، مثل إحداثيات الضغط ، أو إحداثيات الضغط اللوغاريتمي ، أو إحداثيات سيجما . علاوة على ذلك، يمكن تحليل متغيرات السرعة ودرجة الحرارة والجهد الجيوديناميكي إلى مكونات متوسطة ومكونات اضطراب باستخدام تحليل رينولدز .

إحداثيات الضغط في المستوى المماسي الرأسي الديكارتي

في هذه الصيغة، يُختار الضغط كإحداثي رأسي، وتُكتب الإحداثيات الأفقية للمستوى المماس الديكارتي (أي المستوى المماس لنقطة ما على سطح الأرض). لا تأخذ هذه الصيغة انحناء الأرض في الحسبان، ولكنها مفيدة لتصور بعض العمليات الفيزيائية الداخلة في صياغة المعادلات نظرًا لبساطتها النسبية.

لاحظ أن المشتقات الزمنية ذات الحرف D الكبير هي مشتقات مادية . يتكون النظام من خمس معادلات في خمسة مجاهيل.

دuدت-وv=-Φx{\displaystyle {\frac {Du}{Dt}}-fv=-{\frac {\partial \Phi }{\partial x}}}
دvدت+وu=-Φy{\displaystyle {\frac {Dv}{Dt}}+fu=-{\frac {\partial \Phi }{\partial y}}}
  • معادلة الضغط الهيدروستاتيكي ، وهي حالة خاصة من معادلة الزخم الرأسي حيث يُعتبر التسارع الرأسي مهملاً:
0=-Φص-Rتيص{\displaystyle 0=-{\frac {\partial \Phi }{\partial p}}-{\frac {RT}{p}}}
  • معادلة الاستمرارية ، التي تربط التباعد/التقارب الأفقي بالحركة الرأسية في ظل التقريب الهيدروستاتيكي (دص=-ρدΦ{\displaystyle dp=-\rho \,d\Phi }):
ux+vy+ωص=0{\displaystyle {\frac {\partial u}{\partial x}}+{\frac {\partial v}{\partial y}}+{\frac {\partial \omega }{\partial p}}=0}
تيت+uتيx+vتيy+ω(تيص-Rتيصجص)=ججص{\displaystyle {\frac {\partial T}{\partial t}}+u{\frac {\partial T}{\partial x}}+v{\frac {\partial T}{\partial y}}+\omega \left({\frac {\partial T}{\partial p}}-{\frac {RT}{pc_{p}}}\right)={\frac {J}{c_{p}}}}

عند تضمين بيان لحفظ مادة بخار الماء، يمكن لهذه المعادلات الست أن تشكل الأساس لمخطط التنبؤ العددي بالطقس.

المعادلات الأولية باستخدام نظام إحداثيات سيجما، والإسقاط القطبي المجسم

وفقًا لدليل خدمة الأرصاد الجوية الوطنية رقم 1 - منتجات الفاكس ، يمكن تبسيط المعادلات الأولية إلى المعادلات التالية:

  • الرياح المدارية:
uت=ηv-Φx-جصθπx-zuσ-(u2+v22)x{\displaystyle {\frac {\partial u}{\partial t}}=\eta v-{\frac {\partial \Phi }{\partial x}}-c_{p}\theta {\frac {\partial \pi }{\partial x}}-z{\frac {\partial u}{\partial \sigma }}-{\frac {\partial ({\frac {u^{2}+v^{2}}{2}})}{\partial x}}}
  • الرياح المدارية:
vت=-ηuv-Φy-جصθπy-zvσ-(u2+v22)y{\displaystyle {\frac {\partial v}{\partial t}}=-\eta {\frac {u}{v}}-{\frac {\partial \Phi }{\partial y}}-c_{p}\theta {\frac {\partial \pi }{\partial y}}-z{\frac {\partial v}{\partial \sigma }}-{\frac {\partial ({\frac {u^{2}+v^{2}}{2}})}{\partial y}}}
  • درجة حرارة:
تيت=تيت+uتيx+vتيy+wتيz{\displaystyle {\frac {\partial T}{\partial t}}={\frac {\partial T}{\partial t}}+u{\frac {\partial T}{\partial x}}+v{\frac {\partial T}{\partial y}}+w{\frac {\partial T}{\partial z}}}

يمثل الحد الأول التغير في درجة الحرارة الناتج عن الإشعاع الشمسي الوارد والإشعاع الموجي الطويل الصادر، والذي يتغير مع مرور الوقت على مدار اليوم. أما الحدود الثاني والثالث والرابع فتمثل تغيرات درجة الحرارة الناتجة عن انتقال الحرارة. بالإضافة إلى ذلك، يمثل المتغير T ذو الرمز السفلي التغير في درجة الحرارة على ذلك المستوى. كل قيمة T تختلف عن الأخرى وترتبط بالمستوى الذي تمثله. يتم قسمة هذه القيمة على المسافة بين نقاط الشبكة للحصول على التغير في درجة الحرارة مع تغير المسافة. عند ضرب هذه القيمة في سرعة الرياح على ذلك المستوى، فإن وحدتي كلفن لكل متر ومتر لكل ثانية تعطيان كلفن لكل ثانية. مجموع جميع التغيرات في درجة الحرارة الناتجة عن الحركة في الاتجاهات س ، ص ، ع يعطي التغير الكلي في درجة الحرارة مع مرور الوقت.

  • المياه القابلة للترسيب:
دلتادبليوت=uدبليوx+vدبليوy+wدبليوz{\displaystyle {\frac {\delta W}{\partial t}}=u{\frac {\partial W}{\partial x}}+v{\frac {\partial W}{\partial y}}+w{\frac {\partial W}{\partial z}}}

تعمل هذه المعادلة والرموز بطريقة مشابهة لمعادلة درجة الحرارة. تصف هذه المعادلة حركة الماء من مكان إلى آخر عند نقطة معينة دون الأخذ في الاعتبار الماء الذي يتغير شكله. داخل نظام معين، يكون التغير الكلي في الماء مع الزمن صفرًا. ومع ذلك، يُسمح بتحركات التركيز مع الرياح.

  • سمك الضغط:
تصσ=uxxصσ+vyyصσ+wzzصσ{\displaystyle {\frac {\partial }{\partial t}}{\frac {\partial p}{\partial \sigma }}=u{\frac {\partial }{\partial x}}x{\frac {\partial p}{\partial \sigma }}+v{\frac {\partial }{\partial y}}y{\frac {\partial p}{\partial \sigma }}+w{\frac {\partial }{\partial z}}z{\frac {\partial p}{\partial \sigma }}}

تُسهّل هذه التبسيطات فهم ما يحدث في النموذج. فمثلاً، تنتقل عناصر مثل درجة الحرارة (درجة الحرارة الكامنة)، وكمية الماء القابلة للتكثيف، وإلى حد ما سُمك الضغط، من نقطة إلى أخرى على الشبكة مع الرياح. ويتم التنبؤ بالرياح بطريقة مختلفة قليلاً، حيث تستخدم الجهد الجيوبوتنشالي، والحرارة النوعية، ودالة إكسنر π ، والتغير في إحداثي سيجما.

حل المعادلات الأولية الخطية

يتضمن الحل التحليلي للمعادلات الأولية الخطية تذبذبًا جيبيًا في الوقت وخط الطول، يتم تعديله بواسطة معاملات تتعلق بالارتفاع وخط العرض.

{u،v،Φ}={u^،v^،Φ^}هـأنا(sλ+σت){\displaystyle {\begin{Bmatrix}u,v,\Phi \end{Bmatrix}}={\begin{Bmatrix}{\hat {u}},{\hat {v}},{\hat {\Phi }}\end{Bmatrix}}e^{i(s\lambda +\sigma t)}}

حيث s وσ{\displaystyle \sigma }يمثل كل من و عدد الموجة في المنطقة والتردد الزاوي على التوالي. ويمثل الحل الموجات الجوية والمد والجزر .

عندما يتم فصل المعاملات إلى مكونات الارتفاع والعرض، فإن اعتماد الارتفاع يأخذ شكل الموجات المنتشرة أو المتلاشية (حسب الظروف)، بينما يتم إعطاء اعتماد العرض بواسطة دوال هوف .

لا يمكن التوصل إلى هذا الحل التحليلي إلا بعد تبسيط المعادلات الأساسية وتخطيطها. مع ذلك، فإن العديد من هذه التبسيطات (مثل انعدام التبديد، وتجانس درجة حرارة الغلاف الجوي) لا تتوافق مع ظروف الغلاف الجوي الفعلي. ونتيجة لذلك، غالبًا ما يُحسب الحل العددي الذي يأخذ هذه العوامل في الحسبان باستخدام نماذج الدوران العام ونماذج المناخ .

انظر أيضاً

مراجع

  • بينيستون، مارتن. من الاضطراب إلى المناخ: دراسات عددية للغلاف الجوي باستخدام تسلسل هرمي من النماذج. برلين: سبرينغر، 1998. ISBN 3-540-63495-9
  • فيرث، روبرت. بناء شبكة نموذج الأرصاد الجوية متوسطة النطاق وصغيرة النطاق ودقتها. LSMSA، 2006.
  • تومسون، فيليب. التحليل والتنبؤ العددي للطقس. نيويورك: شركة ماكميلان، 1961.
  • بيلكي، روجر أ. نمذجة الأرصاد الجوية متوسطة النطاق. أورلاندو: أكاديميك برس، 1984. ISBN 0-12-554820-6
  • وزارة التجارة الأمريكية، الإدارة الوطنية للمحيطات والغلاف الجوي، دائرة الأرصاد الجوية الوطنية. دليل دائرة الأرصاد الجوية الوطنية رقم 1 - منتجات الفاكس. واشنطن العاصمة: وزارة التجارة، 1979.

خدمة الأرصاد الجوية الوطنية - موقع البحث والتدريب التعاوني بجامعة ولاية كارولينا الشمالية، مراجعة المعادلات الأولية .