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

Governing Equations and Rheological Models

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

⁇ (v/ ⁇ t + v) ⁇ v) = - ⁇ p + ⁇ ⁇ ; + / + /

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

  • Bingham model:] hmm = Y0 + mip UM ⁇ ⁇ ⁇ ⁇ 0, where mip is the plastic viscosity and ga ⁇ is the shear rate. This model captures the plug flow behavior of lava — a rigid cap moving over a sheared basal layer.
  • ()Herschel-Bulkley model:] ng = Y0 + K ⁇ gan، حيث K هو مؤشر الاتساق و n هو مؤشر التدفق (n < 1 for shear-thinning, n >) 1 لغاز الصدر - العصيان، وهذا النموذج الأكثر مرونة يمثل السلوك الذي يعتمد على درجة الحرارة والدرجات المتطاولة للمغنمة المعقدة.

القناع نفسه يتطور كبريد وبلورات، وكثيرا ما يستخدم معادلة روسكو - إينستين لتحديث القناع الفعال كوظيفة من وظائف الكسر الكريستالي: ميكرو - ميكرو - ليكسيد - ليكسيد - ليتسن، وهى مشكلة الحرق السطحي - الحرق المكثف - السلوك المكثيف -

عدد من الأساليب المحاكاة لتدفق لافا

Finite Volume Method (FVM)

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

طريقة العناصر النهائية

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

Meshless Methods: Smoothed Particle Hydrodynamics (SPH)

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

النهج المقارنة والنُهج الهجينة

وتتوفر لكل طريقة عددية مواطن القوة والضعف، وتحافظ على البرمجيات وتقويتها بالنسبة للتنبؤات التشغيلية الواسعة النطاق؛ وتتفوق حركة FEM على المرونة والقابلية للتكييف؛ وتُدرج هذه المواد في شكلها الطبيعي لتدفقات سطحية، وتجمع بعض المدونات الحديثة بين هذه الأساليب: حجم محدود أو مذيب لعنصر محدود بالنسبة لغالبية التدفق، مقرونة بنهج معالجة الجسيمات في السطح أو في التكييف.

العمليات المادية الرئيسية في لافا Flow Dynamics

التبريد، والكريستال، وثورة الفيسكوسيتي

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

فتيل الغاز وتدفقه

وهناك تدفقات كثيرة ملاحظــة من الحمم تتضمن تقلبات ملتوية )أساسا H2O، CO2، SO2( تزول كضغط يقترب من الفتحة، ويمكن أن تؤدي البثورات إلى خفض كبير في كثافة الرغاوي ودرجة الارتداد، وأن تعزز سرعة تدفق التدفق وقطعة الارتداد، وكثيرا ما تكون هناك حاجة إلى نماذج تدفقية )بديل + غاز( أو نماذج مختلطة لا تتناسب مع هذا السلوك.

التفاعل الطبوغرافية وتشكيل ليفي

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

بيانات المدخلات وشروط الحدود

وتتطلب محاكاة موثوقة من قبل مركز تنمية الغابات بيانات عالية الجودة، وتتمثل أهم المعايير في ما يلي:

  • (أ) معدل الإهدار (معدل التفريغ): تصريف الحمم في كل وقت من أيام الوحدة، يقاس عادة بال m3/s. ويمكن تقدير ذلك من الصور الحرارية الساتلية أو القياسات الميدانية أو المتوسطات التاريخية.() ويحدد معدل التعبئة كثافة التدفقية ويخضع للمراقبة الأولية على طول التدفق.
  • Initial lava temperature:] For basaltic magmas, the liquidus temperature is:1200°C; for andesitic or dacitic lavas, it is lower (~900-1100°C). The initial temperature affects the cooling rate and prolization kinetics.
  • Topography (DEM): At least 10 m horizontal resolution is recommended for accurate path prediction. Steeper slopes (30°) can cause rapid acceleration and transition from laminar to turbulent flow.
  • Material properties:] density (typically 2500-2800 kg/m3), heat capacity, thermal conductivity, emissivity for radiative cooling, and the rheological parameters (yield stress, consistency index, flow index). Laboratory measurements on natural samples are the gold standard, but in practice these valuesup is tuned against past erld.

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

تقييم ودراسات الحالات الإفرادية

Kīea 2018 منطقة خط الشرق الأدنى

The 2018% eruption of Kīea Volcano on the Island of Hawai `i produced one of the most destructive lava flows in recent history, destroying over 700 structures. Extensive field observations, drone imagery, and satellite data (e.g., from the Hawaii Volcano Observatory and the USGS) provide an excellent benchmarks for CFD models.

جبل إيتنا 2021 باروكسيزمال إيبيسو

وخلال الفترة من شباط/فبراير إلى آذار/مارس ٢٠٢١، لم يكن بوسع جبل إيتنا في صقلية أن يُعد سلسلة من نافورات الحمم المضللة، وكل منها يولد تدفقا قصير الأجل ولكن سريع الحركة، مما أدى إلى حدوث عدة كيلومترات في أسفل سفينة " فالي ديل بوفي " ، وقد وفرت الملاحظات التي أجراها المعهد الوطني الإيطالي للفيزياء وعلم البركان بيانات عالية عن الحل الزمني بشأن معدلات تدفق المياه وتدفقها.

التحديات والحدود

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

  • Compputational cost:] High-resolution 3D simulations of a full eruption can take days or weeks on super computeruters, making real-time forecasting impractical. Even 2D depth-averaged models (e.g., shall require careful mesh design to avoid prohibitive run times.
  • Rheological uncertainty:] The rheology of natural lava is not fully understood, especially for Belgian-rich and bubbly magmas. Laboratory experiments on remelted samples may not represent the behavior of rapidly cooling, gas-charged flows. Uncertainty in yield stress and viscosity can lead to large variations in predicted flow length.
  • Multiscale physics:] Bubble dynamics and Belgian growth occur at scales of micrometers to millimeters, while the flow advances over kilometers. Bridging these scales in a single simulation (multiscale modeling) is an active research topic.
  • Lack of real-time data:] During an ongoing eruption, effusion rate and lava temperature can change rapidly. CFD models that assume constant input fail to capture the waxing and waning phases. Incorporating real-time thermal and positional data from satellites (e.g., VIIRS, Sentinel-2) into a data-way.
  • Topography change:] The lava flow itself modifies the topography by building up new land (e.g., delta formation if lava enters the sea). Most CFD models assume a fixed topography, which may become inaccurate for long-duration eruptions.

الاتجاهات المستقبلية

محاكاة محسنة للتعلم

ويمكن أيضاً تدريب نماذج بديلة تستند إلى التعلم العميق (مثل الشبكات العصبية أو الشبكات العصبية للرسوم البيانية) على مكتبة من المحاكاة السابقة للمؤسسة من أجل إنتاج توقعات قريبة من الحالات، ويمكن لهذه " المصممين " التنبؤ بمسار التدفق والسمكة نظراً لمعدل التعبئة والتضاريس دون حل معادلات الدفن الكاملة - Stokes (20).

استيعاب البيانات والتنبؤات التشغيلية

ويمكن لاستيعاب البيانات في الوقت الحقيقي - ضخ الملاحظات الساتلية مع التنبؤات النموذجية لنموذج CFD - أن يقلل من عدم اليقين ويحسن من المهارات المتوقعة، وقد طُبقت مقادير كلغمان ومرشحات الجسيمات لتكييف البارامترات الارتفاعية مع تقدم التدفق.

نماذج مقترنة لتقييم الأثر

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

خاتمة

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