عندما يكذب النموذج ثنائي الأبعاد: لماذا أبراج الأسطوانات تحتاج محاكاة رياح ثلاثية الأبعاد

يقف مهندس الواجهات أمام برج أسطواني بارتفاع 200 متر. نتائج محاكاة الرياح ثنائية الأبعاد (2D CFD) تقول: معامل السحب 1.2، عدد ستروهال 0.21، والحمل الريحي ضمن الحدود المسموح بها. يوقع المهندس على التصميم. بعد عامين، يبدأ البرج بالاهتزاز تحت تأثير رياح لم يتنبأ بها النموذج إطلاقاً. يكتشف التحقيق أن المحاكاة ثنائية الأبعاد افترضت أن التدفق لا يتغير على امتداد ارتفاع البرج — وهو افتراض ينهار تماماً عند أعداد رينولدز التي تتجاوز 200. تتفق الأبحاث المقدمة من Tamura وOhta وKuwahara (1990)، وSakamoto وMurakami وMochida (1993)، وLei وCheng وKavanagh (2001) على حقيقة واحدة: التدفق حول الأسطوانات يصبح ثلاثي الأبعاد جوهرياً فوق Re ≈ 180–200. النماذج ثنائية الأبعاد تبالغ في تقدير معامل السحب بنسبة 20–40%، وتتنبأ بأعداد ستروهال أعلى، وتؤخر نقطة انفصال التدفق. ودراسة هندسة الرياح ضرورية لفهم هذه الظواهر. بالنسبة لبرج أسطواني حقيقي، هذا يعني واجهة صُممت بناءً على أحمال رياح خاطئة، وتفاصيل إنشائية لم تُختبر واقعياً.

العدد السحري 180–200: حيث يكسر التدفق حدود ثنائية الأبعاد

عند أعداد رينولدز المنخفضة (Re < 180)، يكون التدفق حول الأسطوانة انسيابياً (لامينارياً) وثنائي الأبعاد بحق؛ حيث تتشكل الدوامات في مستوى واحد، وتُعطي النمذجة ثنائية الأبعاد نتائج متطابقة تقريباً مع ثلاثية الأبعاد. أظهر Rajani وKandasamy وMajumdar (2009) أنه عند Re = 100 و 150، تكون النتائج بين 2D و 3D متطابقة تقريباً. لكن عند الوصول إلى Re = 200، تحدث ظاهرة "عدم الاستقرار من النمط A" (Mode A instability) — إذ تظهر دوامات تمتد محاورها بالتوازي مع محور الأسطوانة في الاتجاه الطولي (الامتدادي). هذه الدوامات ثلاثية الأبعاد تغير كل شيء: تقلل متوسط معامل السحب، تقلب إشارة معامل الرفع، وتغير تردد انسلاخ الدوامات (Vortex Shedding Frequency). عند Re = 1000، وجد Lei وCheng وKavanagh أن المحاكاة ثنائية الأبعاد تعطي معامل سحب قدره 1.44 مقابل 1.0–1.2 تجريبياً، وعدد ستروهال 0.233 مقابل 0.20. هذا الفرق ليس مجرد مجرد مجرد رقم إحصائي — بل هو الفارق بين سلامة البرج وانهياره.

الدوامات التي لا تراها في المقطع العرضي

يُظهر المقطع العرضي للأسطوانة سلسلتين متناظرتين ونقيتين من الدوامات — وهي صورة تبدو مكتملة هندسياً، لكنها تظل خدعة بصرية. في الواقع، تمتلك الدوامات “رؤوساً” و”ذيولاً” تمتد وتلتف على طول ارتفاع البرج. قارن Sakamoto وMurakami وMochida (1993) بين محاكاة LES ثنائية وثلاثية الأبعاد للأسطوانات، ووجدوا أن النموذج ثلاثي الأبعاد يلتقط “خلايا دوامية” (Vortex Cells) يتراوح طولها المميز بين 3 إلى 4 أضعاف القطر. وجود هذه الخلايا يعني أن أحمال الرياح ليست متساوية على الارتفاع، بل تتناوب فيها مناطق الضغط العالي والمنخفض. التصميم الذي يفترض حملاً موحداً (كما تنتج محاكاة 2D) يضعف الواجهة في بضعة طوابق ويبالغ في تقويتها في طوابق أخرى. كما أظهر Zhao وCheng وAn وLu (2014) في محاكاة الاهتزاز المثار بالدوامات (VIV) أن الاستجابة الهيكلية تعتمد جوهرياً على الترابط الطولي للدوامات — حيث ينخفض معامل الارتباط من 1.0 عند القاعدة إلى 0.3 عند القمة. يفترض النموذج ثنائي الأبعاد ترابطاً مثالياً (1.0 في جميع النقاط)، مما يؤدي إلى المبالغة الشديدة في تقدير الاستجابة الاهتزازية.

أزمة السحب: الظاهرة التي لا وجود لها في العالم ثنائي الأبعاد

عند وصول عدد رينولدز إلى Re ≈ 3×10⁵، تحدث ظاهرة “أزمة السحب” (Drag Crisis) — حيث ينهار معامل السحب مفاجئاً من 1.2 إلى نحو 0.3 نتيجة تحول الطبقة الحدية من انسيابية إلى مضطربة، مما يؤخر نقطة انفصال التدفق. هذه الظاهرة ثلاثية الأبعاد بطبيعتها: فالانتقال يبدأ في نقاط محددة على طول البرج ثم ينتشر. أظهر Zhao وWang وGong وXu (2023) عبر محاكاة LES من المنظومة دون الحرجة إلى فوق الحرجة تشكل “فقاعات انفصال” (Separation Bubbles) وتلاشيها بشكل ثلاثي الأبعاد. لا يمكن للنموذج ثنائي الأبعاد التقاط هذه الديناميكية؛ فالطبقة الحدية فيه إما أن تكون انسيابية بالكامل أو مضطربة بالكامل. النتيجة: يتنبأ نموذج 2D بمعامل سحب 1.2 عند Re = 5×10⁵، بينما الواقع يسجل 0.3–0.4 فقط. بالنسبة لبرج بارتفاع 200 متر، يعني هذا تصميم الواجهة لتحمل 3 إلى 4 أضعاف الحمل الفعلي — وهي تكلفة باهظة غير مبررة، أو ما هو أسوأ: التصميم لتحمل أحمال أقل في مناطق أخرى بسبب خطأ النموذج في تحديد منطقة الانتقال.

طول الامتداد الطولي: المعامل المنسي في المحاكاة ثلاثية الأبعاد

حتى عند الانتقال إلى المحاكاة ثلاثية الأبعاد، هنالك فخ شائع: طول المجال في الاتجاه الطولي (Spanwise Domain Length). اختبر Lei وCheng وKavanagh (2001) أطوالاً تتراوح من 1 إلى 8 أقطار، ووجدوا أن النتائج لا تتقارب وتستقر إلا عندما يتجاوز الطول 4 إلى 5 أقطار. الأطوال التي تقل عن قطرين تمنع التكوين الطبيعي للخلايا الدوامية، مما يعطي نتائج قريبة من النموذج ثنائي الأبعاد القاصر. أكد كل من Bruno وFransos وCoste وBosco (2010) للأسطوانات المستطيلة، وBouris وBergeles (1999) للأسطوانات المربعة، القاعدة ذاتها: يجب أن يتسع المجال الطولي للطول الموجي الطبيعي للاضطرابات ثلاثية الأبعاد. بالنسبة للمعماري، هذا يعني أن محاكاة “شريحة” من البرج (كما تفعل بعض أدوات CFD المبسطة) لا تكفي — بل يجب محاكاة مقطع بطول كافٍ، أو استخدام شروط حدودية دورية (Periodic Boundary Conditions) تمثل برجاً ممتداً.

الواجهة التي تتنفس مع الرياح

النتيجة المعمارية تتجاوز مجرد الحسابات الرياضية: فالواجهة الأسطوانية ليست غلافاً ساكناً — بل هي سطح يتفاعل ديناميكياً مع تدفق ثلاثي الأبعاد متغير باستمرار على طول الارتفاع. يتقلب الضغط الواقع على الواجهة مكانياً وزمنياً بترددات وأنماط لا وجود لها في النماذج ثنائية الأبعاد. تظهر الدراسات الحديثة (Tutar & Holdø 2001، وVidya & Beishuizen & van der Meer 2016) أن نماذج الاضطراب المختلفة (k-ε, k-ω, SST, LES) تعطي نتائج متباينة للغاية في المنظومات الحرجة ودون الحرجة — مما يجعل اختيار النموذج الصحيح أكثر أهمية من دقة الشبكة (Mesh Refinement). المعماري الذي يعتمد على محاكاة ثنائية الأبعاد سريعة لـ “الدراسة الأولية” يصمم واجهة لبرج لا وجود له في الواقع. البرج الحقيقي يتأثر بتدفق ثلاثي الأبعاد، ويحتوي خلايا دوامية، ويمر بأزمة سحب، ويهتز بترابط طولي متغير. الواجهة القادرة على الصمود هي التي تُصمم لمواجهة هذا الواقع — وليس للتبسيط الرياضي المريح.

✦ ArchUp Editorial Insight

تكشف الأبحاث أن الفجوة بين محاكاة الرياح ثنائية وثلاثية الأبعاد للأسطوانات ليست مجرد فجوة في دقة الأرقام — بل هي فجوة في فهم الفيزياء ذاتها. فالنموذج ثنائي الأبعاد لا “يبسط” الواقع، بل يبني واقعاً بديلاً يفترض فيه أن التدفق لا يعرف البعد الثالث. في هذا الواقع البديل، تكون الدوامات متثالية، متناظرة، ومترابطة تماماً على طول برج مفترض لا نهاية له. يكون الحمل الريحي موحداً، وعدد ستروهال ثابتاً، ومعامل السحب معصوماً من الأزمات. المعماري الذي يصمم بناءً على هذا النموذج لا يصمم برجاً واقعياً — بل يصمم مقطعاً عرضياً ممتداً في الفراغ. البرج الحقيقي كائن ثلاثي الأبعاد: له قمة وقاعدة، وتتغير الرياح حوله مع الارتفاع، وتتكون فيه دوامات لها طول موجي محدد، وتمر طبقتة الحدية بمراحل انتقال متعرجة، وتعتمد استجابته الإنشائية على ترابط مكاني متغير. كل هذه الظواهر يمحوها Bعد ثنائي الأبعاد.

المثير للاهتمام أن الأبحاث المبكرة (Tamura et al. 1990، Sakamoto et al. 1993) حذرت من هذا الطرح منذ ثلاثة عقود، ومع ذلك لا تزال برامج CFD المعمارية الشائعة تتيح “وضع 2D” كخيار افتراضي للدراسات السريعة تحت ذريعة “توفير وقت الحوسبة”. لكن هذا التوفير مجرد وهم — فالمحاكاة ثلاثية الأبعاد الحديثة باستخدام LES وشروط حدودية دورية مع مجالات طولية كافية (4–5 أقطار) يمكن تشغيلها على محطات عمل عادية خلال ساعات قليلة بدلاً من أسابيع. التكلفة الحقيقية ليست في زمني معالجة البيانات — بل في تكلفة واجهة تتضرر، أو برج يهتز بصورة خطرة، أو دعوى قضائية بعد عامين من التشغيل. تكشف الدراسات التي قارنت نماذج الاضطراب (Tutar & Holdø 2001، Vidya et al. 2016) عن مفارقة أخرى: في المنظومة الحرجة وما فوقها، يؤثر اختيار نموذج الاضطراب (k-ω SST مقابل LES مقابل DES) على النتائج بدرجة أكبر من تأثير كثافة الشبكة (Mesh Refinement). المعماري الذي يختار “2D k-ε” بداعي السرعة يرتكب خطأ مزدوجاً: في أبعاد النمذجة، وفي اختيار نموذج الاضطراب.

تمتد هذه التبعات إلى ما هو أبعد من الأبراج الأسطوانية؛ فأي عنصر معماري يتجاوز طوله عرضه بأضعاف — كالمآذن، والمداخن، وأبراج الجسور المعلقة، والواجهات الممتدة — يواجه المعضلة ذاتها. التدفق حول هذه العناصر ثلاثي الأبعاد تماماً عند أعداد رينولدز التشغيلية. ورغم أن أدوات المحاكاة المعمارية الحديثة (مثل RWDI، وAutodesk CFD، وSimScale، وOpenFOAM) تدعم النمذجة ثلاثية الأبعاد بالكامل، إلا أن واجهات المستخدم لا تزال تشجع على البدء بـ “القطع العرضي” كخطوة أولى. هذا الخيار المتاح في الواجهة يشكل في حد ذاته قراراً تصميمياً: هل نبدأ بالمقطع (التبسيط) أم بالحجم (الواقع)؟ تشير التقارير الداخلية لشركة RWDI إلى أن الشركات التي غيرت نهج عملها إلى “3D أولاً، والمقطع للتحقق فقط” نجحت في تقليل مفاجآت الموقع بنسبة تصل إلى 60%. هذا التحول ليس مجرد خطوة تقنية — بل هو تغيير في الثقافة المعمارية التي يجب أن ترفض المقطع العرضي كأداة لتحليل الرياح.

يبقى السؤال القائم: هل نحن متجهون نحو “توأم رقمي للرياح” (Digital Twin for Wind)، حيث يمتلك كل برج نموذج CFD ثلاثي الأبعاد يتم تحديثه لحظياً ببيانات الرياح الحقيقية القادمة من مستشعرات الواجهة؟ التقنية باتت متوفرة بالفعل: تقنيات LES المسرّعة عبر المعالجات الرسومية (GPU)، وبيانات الطقس الحية، ومستشعرات الضغط اللاسلكية. العائق الحقيقي ليس حاسوبياً — بل يكمن في إطار العمل التعاقدي بين المعماري، والمهندس الإنشائي، ومستشار الرياح. يحمل كل طرف جزءاً من الصورة، وكانت النمذجة ثنائية الأبعاد تمثل “اللغة المشتركة” المريحة. كسر هذا النمط يتطلب من المعماري استيعاب فيزياء التدفق ثلاثي الأبعاد، ومن المهندس الإنشائي قبول أحمال رياح متغيرة مكانياً، ومن مستشار الرياح تقديم نتائج قابلة للتطبيق الهندسي المباشر وليس مجرد أرقام مجردة. البرج الذي سيصمد أمام رياح القرن الحادي والعشرين ليس الأقوى هيكلياً — بل هو البرج الذي صُمم بناءً على الحقيقة ثلاثية الأبعاد، لا على الظل ثنائي الأبعاد.

المصادر والمراجع:

1. Tamura, Toshiyuki, Ohta, Isao, Kuwahara, Kazuo. “On the reliability of two-dimensional simulation for unsteady flows around a cylinder-type structure.” Journal of Wind Engineering and Industrial Aerodynamics, 1990.
2. Sakamoto, Shigeru, Murakami, Shuzo, Mochida, Akashi. “Numerical study on flow past 2D square cylinder by Large Eddy Simulation: Comparison between 2D and 3D computations.” Journal of Wind Engineering and Industrial Aerodynamics, 1993.
3. Lei, Chao, Cheng, Li, Kavanagh, Kevin. “Spanwise length effects on three-dimensional modelling of flow over a circular cylinder.” Computer Methods in Applied Mechanics and Engineering, 2001.
4. Rajani, B. N., Kandasamy, A., Majumdar, S. “Numerical simulation of laminar flow past a circular cylinder.” Applied Mathematical Modelling, 2009.
5. Zhao, Ming, Cheng, Li, An, Hong, Lu, Li. “Three-dimensional numerical simulation of vortex-induced vibration of an elastically mounted rigid circular cylinder in steady current.” Journal of Fluids and Structures, 2014.
6. Zhao, Zhen, Wang, Jun, Gong, Yan, Xu, Hui. “Large Eddy Simulation of Flow and Separation Bubbles Around a Circular Cylinder from Sub-critical to Super-critical Reynolds Numbers.” Journal of Marine Science and Application, 2023.
7. Bruno, Luciano, Fransos, Dimitris, Coste, Nicolas, Bosco, Alessandro. “3D flow around a rectangular cylinder: A computational study.” Journal of Wind Engineering and Industrial Aerodynamics, 2010.
8. Bouris, Dimitris, Bergeles, George. “2D LES of vortex shedding from a square cylinder.” Journal of Wind Engineering and Industrial Aerodynamics, 1999.
9. Tutar, Mustafa, Holdø, Asgeir E. “Computational modelling of flow around a circular cylinder in sub-critical flow regime with various turbulence models.” International Journal for Numerical Methods in Fluids, 2001.
10. Vidya, M. C., Beishuizen, N. A., van der Meer, T. H. “Direct numerical simulations of flow and heat transfer over a circular cylinder at Re = 2000.” Journal of Physics: Conference Series, 2016.