The convective Brinkman--Forchheimer equations or the Navier--Stokes equations with damping in bounded or periodic domains $\subset\mathbb{R}^d$, $2\leq d\leq 4$ are considered in this work. The existence and uniqueness of a global weak solution in the Leray-Hopf sense satisfying the energy equality to the system: $$\partial_t\boldsymbol{u}-μÎ\boldsymbol{u}+(\boldsymbol{u}\cdot\nabla)\boldsymbol{u}+α\boldsymbol{u}+β|\boldsymbol{u}|^{r-1}\boldsymbol{u}+\nabla p=\boldsymbol{f},\ \nabla\cdot\boldsymbol{u}=0,$$ (for all values of $β>0$ and $μ>0$, whenever the absorption exponent $r>3$ and $2βμ\geq 1$, for the critical case $r=3$) is proved. We exploit the monotonicity as well as the demicontinuity properties of the linear and nonlinear operators and the Minty-Browder technique in the proofs. Finally, we discuss the existence of global-in-time strong solutions to such systems in periodic domains.
معرّف الورقة : 2412.20940العنوان : حول معادلات Brinkman-Forchheimer الحمليةالمؤلفون : Sagar Gautam, Manil T. Mohan (معهد الهند للتكنولوجيا روركي)التصنيف : math.AP (التحليل الرياضي)تاريخ النشر : 30 ديسمبر 2024رابط الورقة : https://arxiv.org/abs/2412.20940 تدرس هذه الورقة معادلات Brinkman-Forchheimer الحملية (أو معادلات Navier-Stokes مع التخميد) على مجالات محدودة أو دورية R d \mathbb{R}^d R d (2 ≤ d ≤ 4 2\leq d\leq 4 2 ≤ d ≤ 4 ). يتم إثبات وجود وتفرد الحلول الضعيفة العامة التي تحقق معادلة الطاقة بمعنى Leray-Hopf:
∂ t u − μ Δ u + ( u ⋅ ∇ ) u + α u + β ∣ u ∣ r − 1 u + ∇ p = f , ∇ ⋅ u = 0 \partial_t u - \mu\Delta u + (u \cdot \nabla)u + \alpha u + \beta|u|^{r-1}u + \nabla p = f, \quad \nabla \cdot u = 0 ∂ t u − μ Δ u + ( u ⋅ ∇ ) u + αu + β ∣ u ∣ r − 1 u + ∇ p = f , ∇ ⋅ u = 0
وهذا صحيح لجميع β > 0 \beta > 0 β > 0 و μ > 0 \mu > 0 μ > 0 عندما يكون مؤشر الامتصاص r > 3 r > 3 r > 3 ، وفي الحالة الحرجة r = 3 r = 3 r = 3 يتطلب 2 β μ ≥ 1 2\beta\mu \geq 1 2 β μ ≥ 1 . يستخدم الإثبات رتابة المؤثرات الخطية واللاخطية وشبه الاستمرارية وتقنية Minty-Browder. وأخيراً، يتم مناقشة وجود الحلول القوية العامة لمثل هذه الأنظمة على المجالات الدورية.
تحديات معادلات Navier-Stokes الكلاسيكية : وجود وتفرد الحلول القوية العامة لمعادلات Navier-Stokes غير القابلة للضغط ثلاثية الأبعاد يعتبر من أهم المشاكل غير المحلولة في الرياضيات (مسائل الألفية من معهد Clay).دراسة النماذج المعدلة : في السنوات الأخيرة، قام الرياضيون بدراسة القابلية العامة للحل من خلال تعديل معادلات NSE ثلاثية الأبعاد الكلاسيكية، بما في ذلك:معادلات Navier-Stokes مع حد التخميد r ∣ u ∣ r − 1 u r|u|^{r-1}u r ∣ u ∣ r − 1 u معادلات NSE ثلاثية الأبعاد المروضة المسائل المعدلة مع حد الامتصاص ∣ u ∣ r − 2 u |u|^{r-2}u ∣ u ∣ r − 2 u الأهمية الفيزيائية لمعادلات Brinkman-Forchheimer : تصف حركة السوائل غير القابلة للضغط في الوسائط المسامية المشبعة، حيث:μ \mu μ : معامل Brinkman (اللزوجة الفعالة)α \alpha α : معامل Darcy (نفاذية الوسط المسامي)β \beta β : معامل Forchheimer (يتناسب مع معدل مسام المادة)تحسين النظرية : توجد صعوبات تقنية في الأدبيات الموجودة في حالة المجالات المحدودة، خاصة في إنشاء معادلة الطاقةابتكار الطرق : استخدام نظرية الرتابة وتقنية Minty-Browder لتوفير طرق إثبات جديدةتحليل الحالات الحرجة : دراسة متعمقة للحالة الحرجة r = 3 r = 3 r = 3 وشروطهاوجود وتفرد الحلول الضعيفة العامة : إثبات وجود وتفرد الحلول الضعيفة العامة للنظام عند β > 0 \beta > 0 β > 0 , μ > 0 \mu > 0 μ > 0 للحالة r > 3 r > 3 r > 3 تحليل شامل للحالة الحرجة : إنشاء وجود وتفرد الحلول الضعيفة العامة للحالة الحرجة r = 3 r = 3 r = 3 تحت الشرط 2 β μ ≥ 1 2\beta\mu \geq 1 2 β μ ≥ 1 إثبات صارم لمعادلة الطاقة : استخدام تقنيات تقريب خاصة (بناءً على تقريب دوال ذاتية لمؤثر Stokes) لإنشاء معادلة الطاقة التي يحققها الحل الضعيفتطبيق نظرية الرتابة : إنشاء منهجي لرتابة وشبه استمرارية المؤثرات الخطية واللاخطية، وتطبيق تقنية Minty-Browderنتائج الانتظامية للحلول القوية : إنشاء نظرية وجود وانتظامية الحلول القوية على المجالات الدوريةدراسة نظام معادلات Brinkman-Forchheimer الحملية التالي:
{ ∂ u ∂ t − μ Δ u + ( u ⋅ ∇ ) u + α u + β ∣ u ∣ r − 1 u + ∇ p = f في O × ( 0 , T ) ∇ ⋅ u = 0 في O × ( 0 , T ) u = 0 في ∂ O × ( 0 , T ) u ( 0 ) = u 0 في O \begin{cases}
\frac{\partial u}{\partial t} - \mu\Delta u + (u \cdot \nabla)u + \alpha u + \beta|u|^{r-1}u + \nabla p = f & \text{في } O \times (0,T) \\
\nabla \cdot u = 0 & \text{في } O \times (0,T) \\
u = 0 & \text{في } \partial O \times (0,T) \\
u(0) = u_0 & \text{في } O
\end{cases} ⎩ ⎨ ⎧ ∂ t ∂ u − μ Δ u + ( u ⋅ ∇ ) u + αu + β ∣ u ∣ r − 1 u + ∇ p = f ∇ ⋅ u = 0 u = 0 u ( 0 ) = u 0 في O × ( 0 , T ) في O × ( 0 , T ) في ∂ O × ( 0 , T ) في O
حيث O ⊂ R d O \subset \mathbb{R}^d O ⊂ R d (2 ≤ d ≤ 4 2 \leq d \leq 4 2 ≤ d ≤ 4 ) مجال محدود أو دوري.
تعريف فضاءات الدوال الأساسية:
V = { u ∈ C 0 ∞ ( O , R d ) : ∇ ⋅ u = 0 } V = \{u \in C_0^{\infty}(O,\mathbb{R}^d) : \nabla \cdot u = 0\} V = { u ∈ C 0 ∞ ( O , R d ) : ∇ ⋅ u = 0 } H H H : إغلاق V V V في L 2 ( O ) L^2(O) L 2 ( O ) V V V : إغلاق V V V في H 0 1 ( O ) H_0^1(O) H 0 1 ( O ) L ~ p \tilde{L}^p L ~ p : إغلاق V V V في L p ( O ) L^p(O) L p ( O ) مؤثر Stokes: A u = − P Δ u Au = -P\Delta u A u = − P Δ u ، حيث P P P هو الإسقاط المتعامد Helmholtz-Hodge.
الشكل ثلاثي الخطية: b ( u , v , w ) = ∫ O ( u ⋅ ∇ ) v ⋅ w d x b(u,v,w) = \int_O (u \cdot \nabla)v \cdot w dx b ( u , v , w ) = ∫ O ( u ⋅ ∇ ) v ⋅ w d x
المؤثر المقابل: B ( u , v ) = P ( u ⋅ ∇ ) v B(u,v) = P(u \cdot \nabla)v B ( u , v ) = P ( u ⋅ ∇ ) v , B ( u ) = B ( u , u ) B(u) = B(u,u) B ( u ) = B ( u , u )
C ( u ) = P ( ∣ u ∣ r − 1 u ) C(u) = P(|u|^{r-1}u) C ( u ) = P ( ∣ u ∣ r − 1 u ) ، يحقق ⟨ C ( u ) , u ⟩ = ∥ u ∥ L r + 1 r + 1 \langle C(u),u \rangle = \|u\|_{L^{r+1}}^{r+1} ⟨ C ( u ) , u ⟩ = ∥ u ∥ L r + 1 r + 1
النظرية 2.5 : بالنسبة لـ r > 3 r > 3 r > 3 ، المؤثر G ( u ) = μ A u + B ( u ) + β C ( u ) G(u) = \mu Au + B(u) + \beta C(u) G ( u ) = μ A u + B ( u ) + βC ( u ) يحقق:
⟨ G ( u ) − G ( v ) , u − v ⟩ + ϱ ∥ u − v ∥ H 2 ≥ 0 \langle G(u) - G(v), u - v \rangle + \varrho \|u - v\|_H^2 \geq 0 ⟨ G ( u ) − G ( v ) , u − v ⟩ + ϱ ∥ u − v ∥ H 2 ≥ 0
حيث ϱ = r − 3 2 μ ( r − 1 ) ( 2 β μ ( r − 1 ) ) 2 r − 3 \varrho = \frac{r-3}{2\mu(r-1)}\left(\frac{2}{\beta\mu(r-1)}\right)^{\frac{2}{r-3}} ϱ = 2 μ ( r − 1 ) r − 3 ( β μ ( r − 1 ) 2 ) r − 3 2
النظرية 2.6 : بالنسبة للحالة الحرجة r = 3 r = 3 r = 3 و 2 β μ ≥ 1 2\beta\mu \geq 1 2 β μ ≥ 1 ، المؤثر G G G رتيب عام.
استخدام تقنية التنعيم (mollification) وتسلسل تقريب خاص:
u n ( t ) = P 1 / n u ( t ) = ∑ λ j < n 2 e − λ j / n ⟨ u ( t ) , w j ⟩ w j u_n(t) = P_{1/n}u(t) = \sum_{\lambda_j < n^2} e^{-\lambda_j/n}\langle u(t), w_j \rangle w_j u n ( t ) = P 1/ n u ( t ) = ∑ λ j < n 2 e − λ j / n ⟨ u ( t ) , w j ⟩ w j
هذا التقريب يتقارب في كل من فضاء Sobolev وفضاء Lebesgue، مما يتغلب على الصعوبات التقنية في المجالات المحدودة.
من خلال عدم المساواة الرتيبة:
∫ 0 T e − 2 ϱ t ⟨ G ( v ( t ) ) + ϱ v ( t ) , v ( t ) − u ( t ) ⟩ d t ≥ 0 \int_0^T e^{-2\varrho t}\langle G(v(t)) + \varrho v(t), v(t) - u(t) \rangle dt \geq 0 ∫ 0 T e − 2 ϱt ⟨ G ( v ( t )) + ϱ v ( t ) , v ( t ) − u ( t )⟩ d t ≥ 0
مع شبه الاستمرارية، إثبات G 0 ( t ) = G ( u ( t ) ) G_0(t) = G(u(t)) G 0 ( t ) = G ( u ( t )) .
هذه الورقة عبارة عن بحث نظري بحت، يتم التحقق من النتائج بشكل أساسي من خلال الإثبات الرياضي:
تقريب Faedo-Galerkin : بناء تسلسل حلول تقريبية ذات أبعاد محدودةتقديرات الطاقة : إنشاء تقديرات سابقة موحدةعملية الحد : استخدام الانضغاطية ونظرية الرتابة لأخذ الحدإثبات التفرد : من خلال الطريقة الطاقية وعدم المساواة Grönwallعدم المساواة الطاقية (القضية 3.4):
sup t ∈ [ 0 , T ] ∥ u n ( t ) ∥ H 2 + μ ∫ 0 T ∥ u n ( t ) ∥ V 2 d t + 2 β ∫ 0 T ∥ u n ( t ) ∥ L ~ r + 1 r + 1 d t ≤ ∥ u 0 ∥ H 2 + 1 μ ∫ 0 T ∥ f ( t ) ∥ V ′ 2 d t \sup_{t \in [0,T]} \|u^n(t)\|_H^2 + \mu \int_0^T \|u^n(t)\|_V^2 dt + 2\beta \int_0^T \|u^n(t)\|_{\tilde{L}^{r+1}}^{r+1} dt \leq \|u_0\|_H^2 + \frac{1}{\mu}\int_0^T \|f(t)\|_{V'}^2 dt sup t ∈ [ 0 , T ] ∥ u n ( t ) ∥ H 2 + μ ∫ 0 T ∥ u n ( t ) ∥ V 2 d t + 2 β ∫ 0 T ∥ u n ( t ) ∥ L ~ r + 1 r + 1 d t ≤ ∥ u 0 ∥ H 2 + μ 1 ∫ 0 T ∥ f ( t ) ∥ V ′ 2 d t
النظرية 3.5 : بالنسبة لـ 2 ≤ d ≤ 4 2 \leq d \leq 4 2 ≤ d ≤ 4 ، u 0 ∈ H u_0 \in H u 0 ∈ H ، f ∈ L 2 ( 0 , T ; V ′ ) f \in L^2(0,T;V') f ∈ L 2 ( 0 , T ; V ′ ) ، عندما يكون r > 3 r > 3 r > 3 ، يوجد حل ضعيف فريد يحقق:
u ∈ C ( [ 0 , T ] ; H ) ∩ L 2 ( 0 , T ; V ) ∩ L r + 1 ( 0 , T ; L ~ r + 1 ) u \in C([0,T];H) \cap L^2(0,T;V) \cap L^{r+1}(0,T;\tilde{L}^{r+1}) u ∈ C ([ 0 , T ] ; H ) ∩ L 2 ( 0 , T ; V ) ∩ L r + 1 ( 0 , T ; L ~ r + 1 )
النظرية 3.6 : بالنسبة للحالة الحرجة r = 3 r = 3 r = 3 و 2 β μ ≥ 1 2\beta\mu \geq 1 2 β μ ≥ 1 ، يوجد حل ضعيف فريد للنظام.
النظرية 4.2 : على المجالات الدورية، بالنسبة لـ u 0 ∈ V u_0 \in V u 0 ∈ V ، f ∈ L 2 ( 0 , T ; H ) f \in L^2(0,T;H) f ∈ L 2 ( 0 , T ; H ) ، الحل الضعيف له انتظامية:
u ∈ L ∞ ( 0 , T ; V ) ∩ L 2 ( 0 , T ; D ( A ) ) ∩ L r + 1 ( 0 , T ; L ~ p ( r + 1 ) ) u \in L^{\infty}(0,T;V) \cap L^2(0,T;D(A)) \cap L^{r+1}(0,T;\tilde{L}^{p(r+1)}) u ∈ L ∞ ( 0 , T ; V ) ∩ L 2 ( 0 , T ; D ( A )) ∩ L r + 1 ( 0 , T ; L ~ p ( r + 1 ) )
النظرية 4.4 : بالنسبة لـ f ∈ W 1 , 1 ( [ 0 , T ] ; H ) f \in W^{1,1}([0,T];H) f ∈ W 1 , 1 ([ 0 , T ] ; H ) ، u 0 ∈ V u_0 \in V u 0 ∈ V و A u 0 ∈ H Au_0 \in H A u 0 ∈ H ، يوجد حل قوي فريد:
u ∈ W 1 , ∞ ( [ 0 , T ] ; H ) , A u ∈ L ∞ ( 0 , T ; H ) u \in W^{1,\infty}([0,T];H), \quad Au \in L^{\infty}(0,T;H) u ∈ W 1 , ∞ ([ 0 , T ] ; H ) , A u ∈ L ∞ ( 0 , T ; H )
عدم تبادل مؤثرات الإسقاط : P P P و − Δ -\Delta − Δ لا يتبادلان بالضرورةتأثير شروط الحدود : P ( ∣ u ∣ r − 1 u ) ∣ ∂ O ≠ 0 P(|u|^{r-1}u)|_{\partial O} \neq 0 P ( ∣ u ∣ r − 1 u ) ∣ ∂ O = 0 فشل المتطابقة (1.5) : المتطابقة الأساسية لا تنطبق في المجالات المحدودةتقنية تقريب خاصة : استخدام التقريب بناءً على الدوال الذاتية لمؤثر Stokesطريقة التنعيم : معالجة التنعيم في الاتجاه الزمنينظرية الرتابة : الاستخدام المنهجي للخصائص الرتيبة للمؤثراتCai & Jiu (2008): إثبات وجود حلول ضعيفة عامة على الفضاء الكامل لأي r ≥ 1 r \geq 1 r ≥ 1 Zhang et al. (2011, 2012): تحسين شروط وجود الحلول القوية Zhou (2012): إنشاء وجود حلول قوية عامة لـ r ≥ 3 r \geq 3 r ≥ 3 Hajduk & Robinson (2017): إثبات وجود حلول ناعمة عامة على الطارة ثلاثية الأبعاد لـ r > 3 r > 3 r > 3 ، الحالة الحرجة تتطلب 4 β μ ≥ 1 4\beta\mu \geq 1 4 β μ ≥ 1 هذه الورقة: تحسين إلى 2 β μ ≥ 1 2\beta\mu \geq 1 2 β μ ≥ 1 Antontsev & de Oliveira (2010): إنشاء وجود حلول ضعيفة عامة، لكن لم يتم حل مشكلة معادلة الطاقة Fefferman et al. (2022): توفير تقنيات تقريب أساسية هذه الورقة: حل شامل لمشكلة معادلة الطاقة في الحالة الحرجة إنشاء نظرية حل ضعيف شاملة لمعادلات Brinkman-Forchheimer الحملية حل مشكلة معادلة الطاقة في الحالة الحرجة r = 3 r = 3 r = 3 إنشاء نظرية انتظامية الحلول القوية على المجالات الدورية توفير إطار إثبات موحد بناءً على نظرية الرتابة تقييد الأبعاد : النتائج مقيدة بـ 2 ≤ d ≤ 4 2 \leq d \leq 4 2 ≤ d ≤ 4 الحلول القوية للمجالات المحدودة : لم يتم إنشاء نظرية الحلول القوية على المجالات المحدودة بسبب الصعوبات التقنيةالشروط الحرجة : عند r = 3 r = 3 r = 3 لا يزال يتطلب الشرط 2 β μ ≥ 1 2\beta\mu \geq 1 2 β μ ≥ 1 التوسع للحالات عالية الأبعاد (d > 4 d > 4 d > 4 ) تحسين نظرية الحلول القوية للمجالات المحدودة مزيد من تخفيف الشروط الحرجة دراسة النسخة العشوائية ابتكار الطرق : تطبيق منهجي لنظرية الرتابة وتقنية Minty-Browder، توفير أفكار إثبات جديدةالدقة التقنية : خاصة في إنشاء معادلة الطاقة، حل مشكلة تقنية طويلة الأمداكتمال النتائج : من الحلول الضعيفة إلى الحلول القوية، من الوجود إلى التفرد، تشكيل نظام نظري شاملالكتابة الواضحة : هيكل الورقة معقول، التفاصيل التقنية مفصلةحدود التطبيق : بحث نظري بشكل أساسي، يفتقر إلى التحقق العددي أو التطبيقات العمليةالتعقيد التقني : عملية الإثبات معقدة نسبياً، قابلية القراءة تحتاج إلى تحسينتقييد الشروط : الحالة الحرجة لا تزال تتطلب افتراضات شروط قوية نسبياًالمساهمة النظرية : مساهمة مهمة في نظرية معادلات ديناميكا السوائل التفاضلية الجزئيةقيمة الطريقة : قد تنطبق طريقة الرتابة على مشاكل أخرى ذات صلةالأبحاث اللاحقة : توفير أساس صلب للبحث الإضافيمشاكل ديناميكا السوائل في الوسائط المسامية دراسة معادلات Navier-Stokes مع التخميد التحليل النظري للمعادلات التفاضلية الجزئية اللاخطية تطبيقات نظرية المؤثرات الرتيبة تستشهد الورقة بـ 46 مرجعاً ذا صلة، تتضمن بشكل أساسي:
مراجع نظرية Navier-Stokes الكلاسيكية (Temam, Ladyzhenskaya وغيرهم) الأبحاث الحديثة حول معادلات Brinkman-Forchheimer المراجع الأساسية لنظرية المؤثرات الرتيبة (Barbu, Brezis وغيرهم) الكتب المعيارية في التحليل الدالي والمعادلات التفاضلية الجزئية