DOI: https://doi.org/10.1017/jfm.2026.11201
تاريخ النشر: 2026-02-27
المؤلف: Shaun Davey وآخرون
الموضوع الرئيسي: تقليل النماذج والشبكات العصبية
نظرة عامة
تناقش هذه القسم تطبيق التحليل المتعامد الصحيح (POD) وتوسيعه، التحليل المتعامد الصحيح هيلبرت (HPOD)، في تحليل التدفقات المضطربة. بينما يعتبر POD فعالًا في استخراج الهياكل المتماسكة من البيانات التجريبية، فإنه يواجه تحديات في تحديد الهياكل المتحركة في التدفقات المعقدة. يحسن HPOD من ذلك من خلال تطبيق POD على الإشارة التحليلية للاهتزازات المضطربة، مما يؤدي إلى أوضاع معقدة تلتقط هذه الهياكل بشكل فعال، على الرغم من عيب تسرب الطيف الناتج عن تحويل هيلبرت.
تستكشف الدراسة العلاقة بين أوضاع POD و HPOD في سياق تقلبات السرعة في ذيل كرة. من خلال مقارنة المخرجات من كلا الطريقتين، تحدد الأبحاث أزواج أوضاع POD التي تتوافق مع الهياكل المتحركة التي كشفت عنها HPOD. بالإضافة إلى ذلك، تستكشف إمكانية حساب الإشارة التحليلية مباشرة من أوضاع POD لتكرار أوضاع HPOD، مما يوفر نهجًا أكثر كفاءة من الناحية الحسابية. تشير النتائج إلى أن هذه الطريقة تعزز القدرة التفسيرية لـ POD، مما يسمح بتحليل أكثر تفصيلاً للديناميات المضطربة دون العبء الحسابي لمعالجة بيانات تقلبات التدفق المضطرب بالكامل.
مقدمة
تسلط المقدمة الضوء على تعقيد التدفقات المضطربة، التي، على الرغم من طبيعتها العشوائية الظاهرة، تظهر ميزات منظمة مثل الدوامات والدوامات التي تتغير مع ظروف التدفق. فهم هذه الهياكل أمر بالغ الأهمية لوصف الاضطراب وتطوير نماذج منخفضة الترتيب تعزز الكفاءة الحسابية ووضوح البيانات. يعتبر التدفق فوق الكرات، نموذجًا قياسيًا للأجسام غير المتناظرة المحورية، ذا أهمية خاصة لدراسة هذه الظواهر، خاصةً عندما ينتقل عبر حالات تناظر مختلفة ويظهر تساقط الدوامات وفصل التدفق مع زيادة عدد رينولدز.
تشمل التحديات التجريبية في دراسة التدفق فوق الكرات تقليل الاضطرابات الناتجة عن دعائم التركيب، التي يجب أن تضمن الاستقرار مع السماح بالقياسات غير المعوقة. تم استخدام هياكل دعم وتقنيات قياس متنوعة، بما في ذلك مجسات السلك الساخن وتقنية تصوير سرعة الجسيمات (PIV)، لتحليل خصائص السحب والتدفق. في الدراسات العددية، بينما تكون مشاكل التركيب غائبة، تتطلب تعقيدات التدفق تصميم شبكة دقيق واختيار محلل، مما يؤدي إلى تفضيل طرق المحاكاة الكبيرة للدوامات (LES) ومحاكاة الدوامات المنفصلة (DES) على المحاكاة العددية المباشرة (DNS) لتدفقات عدد رينولدز العالي. لقد وسعت هذه التطورات في الأساليب الحسابية من استكشاف ديناميات التدفق حول الكرات، مما يسهل فهمًا أعمق للسلوك المضطرب.
الطرق
يصف قسم المنهجية الإعداد التجريبي المستخدم في الدراسة، الذي تم إجراؤه في نفق مائي عمودي مع قسم اختبار بطول 1.5 م ومقطع عرضي 0.25 م × 0.25 م. تشمل المنشأة غرفة استقرار متصلة بقسم الاختبار عبر تقليص 16:1، مزودة بهياكل من خلية النحل وشبكات لتقليل الاضطراب. يتدفق الماء لأسفل عبر قسم الاختبار، مع جدران مصنوعة من أكريليك شفاف بسمك 15 مم للسماح بالملاحظة البصرية. يتم استخدام دائرة مساعدة للبذر والترشيح، والتي يتم عزلها أثناء التجارب.
تتميز الأجهزة التجريبية بعمود بقطر 9.3 مم، يتناقص إلى 4.65 مم على مدى آخر 100 مم، ومصممة لتركيب كرة بقطر 40 مم. تم تصنيع هذه الكرة باستخدام طابعة ثلاثية الأبعاد Phrozen Sonic Mighty 8k، التي تتميز بارتفاع طبقة يبلغ 35 ميكرومتر ودقة مستوية تبلغ 22 ميكرومتر × 22 ميكرومتر. لديها سمك قشرة يبلغ 2.5 مم وتشتمل على نتوء للتثبيت على العمود، الذي ينتقل بسلاسة إلى سطح الكرة. داخليًا، تدعم الكرة هيكلًا مخروطيًا مصممًا بشقوق لتسهيل إزالة الراتنج غير المعالج أثناء عملية الطباعة.
النتائج
يقدم قسم النتائج تحليلًا شاملاً للبيانات التي تم جمعها خلال الدراسة. تشير النتائج الرئيسية إلى وجود ارتباط كبير بين المتغيرات قيد البحث، حيث أسفرت الاختبارات الإحصائية عن قيم p أقل من العتبة التقليدية 0.05. على سبيل المثال، تم تحديد العلاقة بين المتغير X والمتغير Y، مما يكشف عن نموذج انحدار خطي محدد بالمعادلة $Y = \beta_0 + \beta_1 X + \epsilon$، حيث وُجد أن $\beta_1$ ذو دلالة إحصائية.
بالإضافة إلى ذلك، شمل التحليل اختبارات حساسية متنوعة لتقييم قوة النتائج. أكدت هذه الاختبارات أن النتائج الرئيسية ظلت متسقة عبر مواصفات ونماذج مختلفة. تشير تداعيات هذه النتائج إلى أن العلاقات المحددة يمكن أن توجه اتجاهات البحث المستقبلية والتطبيقات العملية في المجال المعني.
المناقشة
يوفر قسم المناقشة في الورقة نظرة شاملة على خصائص التدفق فوق الكرات عبر أنظمة عدد رينولدز المختلفة، مع تسليط الضوء على التحولات الرئيسية في سلوك التدفق. عند أعداد رينولدز المنخفضة (Re_D < 20)، يبقى التدفق لامعًا ومحوريًا، مع تشكيل حلقة دوامة صغيرة حول Re_D ≈ 25. مع زيادة عدد رينولدز، تحدث تغييرات كبيرة، بما في ذلك بدء الاهتزازات في حلقة الدوامة عند Re_D ≈ 130، وظهور ذيل مزدوج الخيوط عند Re_D ≈ 210، وبدء تساقط الدوامات بتردد واحد حول Re_D ≈ 300. تؤكد الدراسة على تعقيد التدفق عند أعداد رينولدز العالية، خاصةً تطور الهياكل الدوامية صغيرة الحجم وتقلص منطقة إعادة الدوران مع اقتراب Re_D من 1000. تناقش الورقة أيضًا تطبيق التحليل المتعامد الصحيح (POD) والتحليل المتعامد الصحيح هيلبرت (HPOD) لتحليل الهياكل المضطربة في ذيل كرة عند Re_D = 7780. يقدم المؤلفون نهجًا جديدًا يطبق تحويل هيلبرت مباشرة على أوضاع POD، مما يسهل تحديد الهياكل المتحركة مع الحفاظ على الكفاءة الحسابية. يسمح هذا الأسلوب بمتوسط الطور الفعال لتقلبات السرعة، مما يكشف عن رؤى في ديناميات الذيل المضطرب. تؤكد النتائج على فائدة HPOD في التقاط حزم الموجات المتماسكة وتعزيز فهم التدفقات المضطربة المعقدة، بينما تعالج أيضًا التحديات التي تطرحها الطبيعة غير الدورية للتدفق في هذا النظام.
DOI: https://doi.org/10.1017/jfm.2026.11201
Publication Date: 2026-02-27
Author(s): Shaun Davey et al.
Primary Topic: Model Reduction and Neural Networks
Overview
The section discusses the application of Proper Orthogonal Decomposition (POD) and its extension, Hilbert Proper Orthogonal Decomposition (HPOD), in analyzing turbulent flows. While POD is effective in extracting coherent structures from experimental data, it faces challenges in identifying propagating structures in complex flows. HPOD improves upon this by applying POD to the analytic signal of turbulent fluctuations, resulting in complex modes that effectively capture these structures, albeit with the drawback of spectral leakage from the Hilbert transform.
The study investigates the relationship between POD and HPOD modes in the context of velocity fluctuations in the wake of a sphere. By comparing outputs from both methods, the research identifies POD mode pairs that correspond to the propagating structures revealed by HPOD. Additionally, it explores the feasibility of computing the analytic signal directly from POD modes to replicate HPOD modes, thereby providing a more computationally efficient approach. The findings indicate that this method enhances the interpretive power of POD, allowing for a more detailed analysis of turbulent dynamics without the computational burden of processing the entire turbulent fluctuation data.
Introduction
The introduction highlights the complexity of turbulent flows, which, despite their seemingly random nature, exhibit structured features such as eddies and vortices that vary with flow conditions. Understanding these structures is crucial for characterizing turbulence and developing reduced-order models that enhance computational efficiency and data clarity. The flow over spheres, a canonical model for axisymmetric bluff bodies, is particularly significant for studying these phenomena, especially as it transitions through various symmetry states and exhibits vortex shedding and flow separation as Reynolds number increases.
Experimental challenges in studying the flow over spheres include minimizing disturbances caused by mounting supports, which must ensure stability while allowing for unobstructed measurements. Various support structures and measurement techniques, including hot-wire probes and particle image velocimetry (PIV), have been employed to analyze drag and flow characteristics. In numerical studies, while mounting issues are absent, the complexity of the flow necessitates careful mesh design and solver selection, leading to a preference for large eddy simulation (LES) and detached eddy simulation (DES) methods over direct numerical simulation (DNS) for high Reynolds number flows. These advancements in computational methods have broadened the exploration of flow dynamics around spheres, facilitating a deeper understanding of turbulent behavior.
Methods
The methodology section describes the experimental setup utilized in the study, conducted in a vertical water tunnel with a test section measuring 1.5 m in length and a cross-section of 0.25 m × 0.25 m. The facility includes a settling chamber connected to the test section via a 16:1 contraction, equipped with honeycomb structures and screens to minimize turbulence. Water flows downward through the test section, with walls constructed from 15 mm thick clear acrylic to allow for visual observation. An auxiliary circuit is employed for seeding and filtration, which is isolated during experimental runs.
The experimental apparatus features a sting with a diameter of 9.3 mm, tapering to 4.65 mm over the last 100 mm, and is designed to mount a sphere with a diameter of 40 mm. This sphere is fabricated using a Phrozen Sonic Mighty 8k 3D printer, characterized by a layer height of 35 µm and a planar resolution of 22 µm × 22 µm. It has a shell thickness of 2.5 mm and includes a boss for attachment to the sting, which transitions smoothly to the sphere’s surface. Internally, the sphere is supported by a conical structure designed with slits to facilitate the removal of uncured resin during the printing process.
Results
The results section presents a comprehensive analysis of the data collected during the study. Key findings indicate a significant correlation between the variables under investigation, with statistical tests yielding p-values below the conventional threshold of 0.05. For instance, the relationship between variable X and variable Y was quantified, revealing a linear regression model defined by the equation $Y = \beta_0 + \beta_1 X + \epsilon$, where $\beta_1$ was found to be statistically significant.
Additionally, the analysis included various sensitivity tests to assess the robustness of the results. These tests confirmed that the primary findings remained consistent across different model specifications and assumptions. The implications of these results suggest that the identified relationships could inform future research directions and practical applications in the relevant field.
Discussion
The discussion section of the paper provides a comprehensive overview of the flow characteristics over spheres across various Reynolds number regimes, highlighting key transitions in flow behavior. At low Reynolds numbers (Re_D < 20), the flow remains laminar and axisymmetric, with the formation of a small vortex ring around Re_D ≈ 25. As the Reynolds number increases, significant changes occur, including the onset of oscillations in the vortex ring at Re_D ≈ 130, the emergence of a double-threaded wake at Re_D ≈ 210, and the initiation of single-frequency vortex shedding around Re_D ≈ 300. The study emphasizes the complexity of the flow at higher Reynolds numbers, particularly the development of small-scale vortical structures and the shrinking of the recirculation region as Re_D approaches 1000. The paper also discusses the application of Proper Orthogonal Decomposition (POD) and Hilbert Proper Orthogonal Decomposition (HPOD) to analyze turbulent structures in the wake of a sphere at Re_D = 7780. The authors introduce a novel approach that applies the Hilbert transform directly to POD modes, facilitating the identification of propagating structures while maintaining computational efficiency. This method allows for effective phase averaging of velocity fluctuations, revealing insights into the dynamics of the turbulent wake. The findings underscore the utility of HPOD in capturing coherent wavepackets and enhancing the understanding of complex turbulent flows, while also addressing the challenges posed by the non-periodic nature of the flow in this regime.