طريقة العناصر المنفصلة الموسعة

تُعدّ طريقة العناصر المنفصلة الموسّعة (XDEM) تقنيةً عدديةً تُوسّع نطاق ديناميكيات المواد الحبيبية أو الجسيمات، كما هو موضح في طريقة العناصر المنفصلة الكلاسيكية (DEM) (كوندال [ 1 ] وآلن [ 2 ] )، بإضافة خصائص أخرى مثل الحالة الديناميكية الحرارية ، والإجهاد / الانفعال ، والمجال الكهرومغناطيسي لكل جسيم. وخلافًا لمفهوم ميكانيكا الأوساط المتصلة ، تهدف طريقة XDEM إلى تحليل الطور الجسيمي بكل العمليات المرتبطة به. فبينما تتنبأ طريقة العناصر المنفصلة بموقع واتجاه كل جسيم في المكان والزمان، تُقدّر طريقة العناصر المنفصلة الموسّعة خصائص إضافية مثل درجة الحرارة الداخلية ، وتوزيع الأنواع ، أو التأثير الميكانيكي على الهياكل.
تاريخ
يمكن اعتبار الديناميكا الجزيئية التي طورها ألدر وآخرون في أواخر الخمسينيات [ 3 ] ورحمان في أوائل الستينيات [ 4 ] بمثابة خطوة أولى نحو طريقة العناصر المنفصلة الموسعة، على الرغم من استبدال القوى الناتجة عن التصادمات بين الجسيمات بجهود الطاقة، على سبيل المثال جهود لينارد جونز للجزيئات والذرات كقوى بعيدة المدى لتحديد التفاعل.
وبالمثل، دُرست التفاعلات الديناميكية الموائعية للجسيمات المعلقة في التدفق. واعتُبرت قوى السحب المؤثرة على الجسيمات بفعل السرعة النسبية بينها وبين التدفق قوىً إضافيةً تؤثر عليها. ولذلك، تُحل ظواهر التدفق متعدد الأطوار هذه ، التي تشمل طورًا صلبًا (مثل الجسيمات) وطورًا غازيًا أو مائعًا، باستخدام طرق منفصلة، بينما يُوصف تدفق الغاز أو السائل باستخدام طرق متصلة، ولذا يُطلق عليه نموذج الوسط المتصل والمنفصل المدمج (CCDM) كما طبقه كاواغوتشي وآخرون [ 5 ] ، وهومانز [ 6 ] ، وشو 1997 [ 7 ] ، وشو 1998 [ 8 ]. وبسبب الوصف المنفصل للطور الصلب، تُحذف العلاقات التكوينية ، مما يؤدي إلى فهم أفضل للأساسيات. وقد خلص إلى ذلك أيضًا تشو 2007 وآخرون [ 9 ] ، وتشو 2008 وآخرون. [ 10 ] خلال مراجعة لتدفقات الجسيمات التي تم نمذجتها باستخدام منهجية CCDM. شهدت هذه المنهجية تطورًا كبيرًا خلال العقدين الماضيين، وتصف حركة الطور الصلب باستخدام طريقة العناصر المنفصلة (DEM) على مستوى الجسيم الفردي، بينما تُعالج الأطوار المتبقية باستخدام معادلات نافيير-ستوكس . وبذلك، تُعتبر هذه المنهجية أداة فعالة لدراسة التفاعل بين الطور الجسيمي والطور السائل، كما استعرضها كل من يو وشو [ 11 ] ، وفينغ ويو [ 12 ] ، ودين وآخرون [ 13 ] . واستنادًا إلى منهجية CCDM، تنبأ غريتشكا وآخرون [ 14 ] بخصائص الأسرّة المميعة والمتدفقة.
طُوِّر الأساس النظري لنموذج XDEM عام 1999 على يد بيترز [ 15 ] ، الذي وصف حرق طبقة خشبية متحركة على شبكة أمامية [ 16 ] . لاحقًا، استخدم سيسميك وآخرون [ 17 ] هذا المفهوم للتنبؤ بعملية الفرن في نظام حرق الشبكة. وقد حاول شونغو وآخرون [ 18 ] تطبيقه على العمليات المعقدة لفرن الصهر. تعتمد العديد من برامج ديناميكا الموائع الحسابية، مثل Simcenter STAR-CCM+ و Ansys وAVL-Fire، على المحاكاة العددية لحقن السوائل في بيئة غازية. وتُعالج قطرات الرذاذ بنهج صفري الأبعاد لحساب انتقال الحرارة والكتلة إلى الطور السائل.
المنهجية

توجد العديد من المشكلات الهندسية التي تتضمن مراحل متصلة ومتقطعة، ولا يمكن محاكاة هذه المشكلات بدقة باستخدام الأساليب المتصلة أو المتقطعة. يوفر برنامج XDEM حلاً لبعض هذه التطبيقات الهندسية.
على الرغم من أن البحث والتطوير في الأساليب العددية في كل من مجالي الحلول المنفصلة والمتصلة لا يزالان في طور التقدم، إلا أن الأدوات البرمجية متوفرة. ولربط الأساليب المنفصلة والمتصلة، يتوفر نهجان رئيسيان:
- النهج المتكامل : يتم حل المعادلات التي تصف الظواهر الفيزيائية المتعددة في وقت واحد بواسطة برنامج حل واحد ينتج عنه حل كامل.
- النهج المجزأ أو المتدرج : يتم حل المعادلات التي تصف الظواهر الفيزيائية المتعددة بشكل متسلسل بواسطة حلول مصممة خصيصًا ومتميزة مع تمرير نتائج تحليل واحد كحمل إلى الآخر.
يتطلب النهج الأول برنامج حلٍّ يُعالج جميع المشكلات الفيزيائية ذات الصلة، ولذلك يتطلب جهدًا أكبر في التنفيذ. مع ذلك، توجد حالات يصعب فيها ترتيب معاملات المعادلات التفاضلية المُجمّعة في مصفوفة واحدة .
يُقدّم النهج الأخير، المُجزّأ، الذي يجمع بين عدد من المُحلّلات التي تُمثّل مجالات فيزيائية مُنفصلة، مزايا مُقارنةً بالنهج المُتكامل. فهو يُوفّر درجةً أكبر من المرونة لقدرته على استخدام العديد من المُحلّلات، كما يُتيح تطوير برمجيات أكثر نمطية. مع ذلك، تتطلّب عمليات المحاكاة المُجزّأة خوارزميات ربط مُستقرة ودقيقة.
في إطار مفهوم XDEM المتداخل، تُوصف الحقول المتصلة بحل معادلات الحفظ المتصلة الخاصة بها. كما تُحدد خصائص الجسيمات الفردية، مثل درجة الحرارة، بحل معادلات الحفظ الخاصة بها، والتي تُنتج توزيعًا مكانيًا وزمنيًا داخليًا للمتغيرات ذات الصلة. يُبين الجدول التالي مبادئ الحفظ الرئيسية، مع معادلاتها والمتغيرات المطلوب حلها، والتي تُطبق على جسيم فردي ضمن XDEM.
| قانون الحفاظ على البيئة | معادلة | عامل |
|---|---|---|
| الكتلة (وسط قابل للانضغاط) | الاستمرارية | الضغط/الكثافة |
| الزخم الخطي | نافيير-ستوكس | سرعة |
| طاقة | طاقة | درجة حرارة |
| كتلة الأنواع | نقل الأنواع | الكسور الكتلية |
| الشحنة، التيار | ماكسويل | المجال الكهربائي، المجال المغناطيسي، مجال الإزاحة الكهربائية |
يُعرّف حلّ هذه المعادلات، من حيث المبدأ، مجالًا ثلاثي الأبعاد وعابرًا للمتغيرات ذات الصلة، مثل درجة الحرارة أو الأنواع الكيميائية. مع ذلك، فإن تطبيق مبادئ الحفظ هذه على عدد كبير من الجسيمات عادةً ما يُقيّد الدقة إلى بُعد واحد تمثيلي وزمن واحد على الأكثر، نظرًا لاستهلاك وقت وحدة المعالجة المركزية. تدعم الأدلة التجريبية، على الأقل في هندسة التفاعلات، فرضية أحادية البُعد، كما أشار مان وبيونغ [ 19 ] ، بينما يُشدّد لي وآخرون [ 20 ] على أهمية السلوك العابر.
التطبيقات

تُعدّ المشكلات التي تتضمن طورًا مستمرًا وآخر منفصلًا ذات أهمية بالغة في تطبيقات متنوعة، مثل صناعة الأدوية (على سبيل المثال، إنتاج الأدوية)، والزراعة، وصناعة الأغذية والتصنيع، والتعدين، والبناء، والآلات الزراعية، وتصنيع المعادن، وإنتاج الطاقة، وعلم الأحياء النظمي. ومن الأمثلة البارزة على ذلك: البن، ورقائق الذرة، والمكسرات، والفحم، والرمل، والوقود المتجدد (مثل الكتلة الحيوية لإنتاج الطاقة والأسمدة).
في البداية، اقتصرت هذه الدراسات على تكوينات التدفق البسيطة كما أشار هومانز [ 21 ] ، إلا أن تشو ويو [ 22 ] أثبتا إمكانية تطبيق هذه الطريقة على تكوين تدفق معقد يتألف من طبقة مميعة وحزام ناقل وفاصل إعصاري. وبالمثل، طبق تشو وآخرون [ 23 ] منهجية CCDM على الهندسة المعقدة لموقد غني/فقير بالوقود لاحتراق الفحم المسحوق في محطة، وقام تشو وآخرون [ 24 ] بنمذجة التدفق المعقد للهواء والماء والفحم وجزيئات الماغنيتيت بأحجام مختلفة في فاصل إعصاري ذي وسط كثيف (DMC).
طُبِّقَ نهج CCDM أيضًا على الأسرّة المميعة، كما استعرضه كلٌّ من رو ونينو [ 25 ] وفينغ ويو [ 26 ] ، وطُبِّقَ من قِبَل فينغ ويو [ 27 ] لدراسة الحركة الفوضوية للجسيمات ذات الأحجام المختلفة في سرير مميع غازي. ويصف كافويا وآخرون [ 28 ] نمذجة الجسيمات المنفصلة للسوائل المتصلة في الأسرّة المميعة الغازية الصلبة. وتشمل التطبيقات الأخرى لـ XDEM التحويل الحراري للكتلة الحيوية على شبكة تعمل للأمام والخلف. كما تم حل ودراسة انتقال الحرارة في الأنظمة الجسيمية الحرارية/المتفاعلة، كما استعرضه بنغ وآخرون [ 29 ] بشكل شامل. ويمثل تشوه حزام ناقل نتيجة اصطدام المواد الحبيبية التي تُفرَغ عبر قناة تطبيقًا في مجال تحليل الإجهاد / الانفعال .
مراجع
- ↑ كوندال، ب.أ.؛ ستراك، أ.د.ل. (1979). "نموذج عددي منفصل للتجمعات الحبيبية". جيوتكنيك . 29 : 47-65 . doi : 10.1680/geot.1979.29.1.47 .
- ↑ ألين، إم بي ؛ تيلدسلي، دي جيه (1990). محاكاة السوائل بالحاسوب . مطبعة كلارندون، أكسفورد.
- ↑ ألدر، بي جيه؛ وينرايت، تي إي (1959). "دراسات في الديناميكا الجزيئية. الجزء الأول: الطريقة العامة". مجلة الفيزياء الكيميائية . 31 (2): 459-466 . رمز Bibcode : 1959JChPh..31..459A . doi : 10.1063/1.1730376 .
- ↑ رحمن، أ. (1964). "الترابطات في حركة الذرات في الأرجون السائل". مجلة الفيزياء ، 136 (2أ): A405– A411. رمز Bibcode : 1964PhRv..136..405R . doi : 10.1103/physrev.136.a405 .
- ↑ كاواغوتشي، ت.؛ تسوجي، ي.؛ تاناكا، ت. (1993). "محاكاة الجسيمات المنفصلة لطبقة مميعة ثنائية الأبعاد". تكنولوجيا المساحيق . 77 : 79-87 . doi : 10.1016/0032-5910(93)85010-7 .
- ↑ هومانز، بي بي بي؛ كويبرز، جيه إيه إم؛ بريلز، دبليو جيه؛ فان سوايج، دبليو بي إم (1996). "محاكاة الجسيمات المنفصلة لتكوين الفقاعات والكتل في طبقة مميعة بالغاز ثنائية الأبعاد: نهج الكرة الصلبة". مجلة هندسة العمليات الكيميائية 51 : 99-118 . CiteSeerX 10.1.1.470.6532 . doi : 10.1016/0009-2509(95)00271-5 . S2CID 17460834 .
- ↑ Xu, BH; Yu, AB (1997). "المحاكاة العددية لتدفق الغاز والمواد الصلبة في طبقة مميعة باستخدام طريقة الجسيمات المنفصلة مع ديناميكيات الموائع الحسابية". مجلة علوم الهندسة الكيميائية . 52 (16): 2785-2809 . doi : 10.1016/s0009-2509(97)00081-x .
- ↑ Xu, BH; Yu, AB (1998). "تعليقات على ورقة بحثية بعنوان: المحاكاة العددية لتدفق الغاز والمواد الصلبة في طبقة مميعة باستخدام طريقة الجسيمات المنفصلة مع ديناميكيات الموائع الحسابية". مجلة علوم الهندسة الكيميائية . 53 (14): 2646-2647 . doi : 10.1016/s0009-2509(98)00086-4 .
- ↑ تشو، إتش بي؛ تشو، زد واي؛ يانغ، آر واي؛ يو، إيه بي (2007). "محاكاة الجسيمات المنفصلة للأنظمة الجسيمية: التطورات النظرية". مجلة علوم الهندسة الكيميائية . 62 (13): 3378-3396 . doi : 10.1016/j.ces.2006.12.089 .
- ↑ تشو، إتش بي؛ تشو، زد واي؛ يانغ، آر واي؛ يو، إيه بي (2008). "محاكاة الجسيمات المنفصلة للأنظمة الجسيمية: مراجعة للتطبيقات والنتائج الرئيسية". مجلة علوم الهندسة الكيميائية . 63 (23): 5728-5770 . doi : 10.1016/j.ces.2008.08.006 .
- ↑ Xu, BH; Yu, AB (2003). "نمذجة تدفق الغاز والمواد الصلبة على مستوى الجسيمات في عمليات التمييع". مجلة التكنولوجيا الكيميائية والتكنولوجيا الحيوية . 78 ( 2-3 ): 111-121 . doi : 10.1002/jctb.788 .
- ↑ فينغ، واي كيو؛ يو، إيه بي؛ يو، إيه بي؛ فينس، إيه. (2004). "تقييم صياغات النماذج في محاكاة الجسيمات المنفصلة لتدفق الغاز والمواد الصلبة". مجلة البحوث الكيميائية والهندسية الصناعية . 43 (26): 8378-8390 . doi : 10.1021/ie049387v .
- ^ دين ، ن.ج. أنالاند، MVS. فان دير هوف، MA؛ كويبرز، جام (2007). “مراجعة نمذجة الجسيمات المنفصلة للأسرة المميعة”. علوم الهندسة الكيميائية . 62 ( 1 – 2): 28 – 44. دوى : 10.1016/j.ces.2006.08.014 .
- ↑ غريتشكا، أ.؛ هاينريش، س.؛ دين، ن. س.؛ فان سينت أنالاند، م.؛ كويبرز، ج. أ. م.؛ مورل، م. (2009). "نمذجة ديناميكا الموائع الحسابية لطبقة مميعة موشورية ذات مدخلين غازيين قابلين للتعديل". المجلة الكندية للهندسة الكيميائية . 87 (2): 318-328 . CiteSeerX 10.1.1.335.4108 . doi : 10.1002/cjce.20143 .
- ↑ بيترز، ب. (1999). "تصنيف أنظمة الاحتراق في طبقة معبأة بناءً على المقاييس الزمنية والمكانية ذات الصلة". الاحتراق واللهب . 116 ( 1-2 ): 297-301 . doi : 10.1016/s0010-2180(98)00048-0 .
- ↑ بيترز، ب. (2002). "قياسات وتطبيق نموذج الجسيمات المنفصلة (DPM) لمحاكاة احتراق طبقة معبأة من جسيمات الوقود الفردية". الاحتراق واللهب . 131 ( 1-2 ): 132-146 . doi : 10.1016/s0010-2180(02)00393-0 .
- ↑ سيمسيك، إي.؛ بروش، ب.؛ ويرتز، س.؛ شيرر، ف.؛ كرول، ف. (2009). "المحاكاة العددية لأنظمة إطلاق النار الشبكية باستخدام طريقة العناصر المنفصلة/ديناميكيات الموائع الحسابية المقترنة (DEM)". تكنولوجيا المساحيق . 193 (3): 266-273 . doi : 10.1016/j.powtec.2009.03.011 .
- ↑ ناتسوي، شونغو؛ أويدا، شيغيرو؛ فان، تشنغيون؛ أندرسون، نيلز؛ كانو، جونيا؛ إينوي، ريو؛ أرياما، تاتسورو (2010). "خصائص تدفق المواد الصلبة وتوزيع الإجهاد، بما في ذلك الظواهر غير المتماثلة في الفرن العالي، تم تحليلها باستخدام طريقة العناصر المنفصلة" . المجلة الدولية للحديد والصلب (ISIJ International ). 50 (2): 207-214 . doi : 10.2355/isijinternational.50.207 .
- ↑ مان، واي إتش؛ بيونغ، آر سي (1994). "دراسة عددية حول احتراق جسيم كربوني منفرد محمول في تدفق ثابت". الاحتراق واللهب . 97 : 1-16 . doi : 10.1016/0010-2180(94)90112-0 .
- ↑ لي، جيه سي؛ ييتر، آر إيه؛ دراير، إف إل (1996). "محاكاة عددية لإشعال جسيم كربوني معزول بالليزر في بيئة ساكنة". الاحتراق واللهب . 105 (4): 591-599 . doi : 10.1016/0010-2180(96)00221-0 .
- ↑ هومانز، بي بي بي؛ كويبرز، جيه إيه إم؛ بريلز، دبليو جيه؛ فان سوايج، دبليو بي إم (1996). "محاكاة الجسيمات المنفصلة لتكوين الفقاعات والكتل في طبقة مميعة بالغاز ثنائية الأبعاد: نهج الكرة الصلبة". مجلة هندسة العمليات الكيميائية 51 : 99-118 . CiteSeerX 10.1.1.470.6532 . doi : 10.1016/0009-2509(95)00271-5 . S2CID 17460834 .
- ↑ تشو، ك.و.؛ يو، أ.ب. (2008). "المحاكاة العددية لتدفقات الجسيمات والسوائل المعقدة". تكنولوجيا المساحيق . 179 (3): 104-114 . doi : 10.1016/j.powtec.2007.06.017 .
- ↑ تشو، هـ.؛ مو، ج.؛ تشاو، ج.؛ سين، ك. (2011). "محاكاة DEM-CFD لتشتت الجسيمات في تدفق ثنائي الطور غاز-صلب لموقد غني/فقير بالوقود". الوقود . 90 (4): 1584-1590 . doi : 10.1016/j.fuel.2010.10.017 .
- ↑ تشو، ك.و.؛ وانغ، ب.؛ يو، أ.ب.؛ فينس، أ.؛ بارنيت، ج.د.؛ بارنيت، ب.ج. (2009). "دراسة باستخدام ديناميكيات الموائع الحسابية وطريقة العناصر المنفصلة لتأثير توزيع كثافة الجسيمات على التدفق متعدد الأطوار وأداء الإعصار ذي الوسط الكثيف". هندسة المعادن . 22 (11): 893-909 . doi : 10.1016/j.mineng.2009.04.008 .
- ↑ رو، ب.ن.؛ نينو، أ.و. (1976). "خلط الجسيمات وفصلها في طبقات التمييع الغازي: مراجعة". تكنولوجيا المساحيق . 15 (2): 141-147 . doi : 10.1016/0032-5910(76)80042-3 .
- ↑ فينغ، واي كيو؛ يو، إيه بي؛ يو، إيه بي؛ فينس، إيه. (2004). "تقييم صياغات النماذج في محاكاة الجسيمات المنفصلة لتدفق الغاز والمواد الصلبة". مجلة البحوث الكيميائية والهندسية الصناعية . 43 (26): 8378-8390 . doi : 10.1021/ie049387v .
- ↑ فينغ، واي كيو؛ يو، إيه بي (2008). "تحليل الحركة الفوضوية للجسيمات ذات الأحجام المختلفة في طبقة مميعة غازية". علم الجسيمات . 6 (6): 549-556 . doi : 10.1016/j.partic.2008.07.011 .
- ↑ كافويا، ك.د.؛ ثورنتون، س.؛ آدامز، م.ج. (2002). "نمذجة الموائع المتصلة بالجسيمات المنفصلة لأسرّة التمييع الغازية الصلبة". مجلة علوم الهندسة الكيميائية . 57 (13): 2395-2410 . doi : 10.1016/s0009-2509(02)00140-9 .
- ↑ بنغ، ز.؛ دورودشي، إ.؛ مقتدري، ب. (2020). "نمذجة انتقال الحرارة في عمليات المحاكاة الحرارية باستخدام طريقة العناصر المنفصلة (DEM): النظرية وتطوير النموذج". التقدم في علوم الطاقة والاحتراق . 79، 100847. doi : 10.1016/j.pecs.2020.100847 . S2CID 218967044 .
- المعادلات التفاضلية العددية
- الفيزياء الحاسوبية



