معادلة برجر

حلول معادلة برجر بدءًا من شرط ابتدائي غاوسيu(x،0)=هـ-x2/2{\displaystyle u(x,0)=e^{-x^{2}/2}}.
حلول من نوع الموجة N لمعادلة برجر، بدءًا من الشرط الأوليu(x،0)=هـ-(x-1)2/2-هـ-(x+1)2/2{\displaystyle u(x,0)=e^{-(x-1)^{2}/2}-e^{-(x+1)^{2}/2}}.

معادلة برجر، أو معادلة باتمان-برجر، هي معادلة تفاضلية جزئية أساسية ومعادلة انتقال حراري [ 1 ] تظهر في مجالات متنوعة من الرياضيات التطبيقية ، مثل ميكانيكا الموائع [ 2 ] ، والصوتيات غير الخطية [ 3 ] ، وديناميكا الغازات ، وتدفق حركة المرور [ 4 ] ، والفيزياء الرياضية [ 5 ] . وقد طُرحت هذه المعادلة لأول مرة من قِبل هاري باتمان عام 1915 [ 6 ] [ 7 ] ، ثم درسها يوهانس مارتينوس برجر عام 1948 [ 8 ]. بالنسبة لحقل معينu(x،ت){\displaystyle u(x,t)}ومعامل الانتشار (أو اللزوجة الحركية ، كما هو الحال في سياق ميكانيكا الموائع الأصلي)ν{\displaystyle \nu }، الشكل العام لمعادلة برجر (المعروفة أيضًا باسم معادلة برجر اللزجة ) في بُعد مكاني واحد هو النظام المبدد :

uت+uux=ν2ux2.{\displaystyle {\frac {\partial u}{\partial t}}+u{\frac {\partial u}{\partial x}}=\nu {\frac {\partial ^{2}u}{\partial x^{2}}}.}

على المدىuu/x{\displaystyle u\partial u/\partial x}ويمكن إعادة كتابتها أيضًا على النحو التالي(u2/2)/x{\displaystyle \partial (u^{2}/2)/\partial x}عندما يغيب مصطلح الانتشار (أيν=0{\displaystyle \nu =0}تصبح معادلة برجر معادلة برجر غير اللزجة :

uت+uux=0،{\displaystyle {\frac {\partial u}{\partial t}}+u{\frac {\partial u}{\partial x}}=0,}

وهو نموذج أولي لمعادلات الحفظ التي يمكن أن تُحدث انقطاعات ( موجات صدمية ).

سبب تشكل التدرجات الحادة للقيم الصغيرة لـν{\displaystyle \nu }يصبح الأمر واضحًا بشكل بديهي عند فحص الجانب الأيسر من المعادلة. المصطلح/ت+u/x{\displaystyle \partial /\partial t+u\partial /\partial x}من الواضح أن هذا عامل موجي يصف موجة تنتشر في الاتجاه الموجبx{\displaystyle x}- اتجاه بسرعةu{\displaystyle u}بما أن سرعة الموجة هيu{\displaystyle u}المناطق التي تُظهر قيمًا كبيرة منu{\displaystyle u}ستنتشر هذه المناطق نحو اليمين بشكل أسرع من المناطق التي تُظهر قيمًا أصغر لـu{\displaystyle u}بمعنى آخر، إذاu{\displaystyle u}يتناقص فيx{\displaystyle x}- في الاتجاه، في البداية، ثم أكبرu{\displaystyle u}سيلحق من يقبعون في الخلف بالصغارu{\displaystyle u}'s على الجانب الأمامي. دور الحد الانتشارى على الجانب الأيمن هو أساساً منع التدرج من أن يصبح لانهائياً.

معادلة البرغر اللزج

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

uت+uux=0،u(x،0)=و(x){\displaystyle {\frac {\partial u}{\partial t}}+u{\frac {\partial u}{\partial x}}=0,\quad u(x,0)=f(x)}

يمكن بناؤها باستخدام طريقة الخصائص . لنفترضت{\displaystyle t}ليكن المعيار الذي يميز أي خصائص معينة فيx{\displaystyle x}-ت{\displaystyle t}في المستوى، تُعطى المعادلات المميزة بواسطة

دxدت=u،دuدت=0.{\displaystyle {\frac {dx}{dt}}=u,\quad {\frac {du}{dt}}=0.}

يُخبرنا تكامل المعادلة الثانية أنu{\displaystyle u}ثابت على طول الخاصية، ويُظهر تكامل المعادلة الأولى أن الخصائص عبارة عن خطوط مستقيمة، أي

u=ج،x=uت+ξ{\displaystyle u=c,\quad x=ut+\xi }

أينξ{\displaystyle \xi }هي النقطة (أو المعلمة) على المحور السيني ( عند t  =  0) في المستوى x - t الذي يُرسم منه المنحنى المميز.u{\displaystyle u}فيx{\displaystyle x}يمكن معرفة المحور السيني من الشرط الابتدائي ومن حقيقة أنu{\displaystyle u}لا يتغير أثناء تحركنا على طول الخاصية المنبعثة من كل نقطةx=ξ{\displaystyle x=\xi }نكتبu=ج=و(ξ){\displaystyle u=c=f(\xi )}على كل خاصية. لذلك، فإن مجموعة مسارات الخصائص التي يتم تحديدها بواسطةξ{\displaystyle \xi }يكون

x=و(ξ)ت+ξ.{\displaystyle x=f(\xi )t+\xi .}

وبالتالي، يُعطى الحل بواسطة

u(x،ت)=و(ξ)=و(x-uت)،ξ=x-و(ξ)ت.{\displaystyle u(x,t)=f(\xi )=f(x-ut),\quad \xi =xf(\xi )t.}

هذه علاقة ضمنية تحدد حل معادلة برجر غير اللزجة بشرط عدم تقاطع الخصائص. إذا تقاطعت الخصائص، فلن يكون هناك حل كلاسيكي للمعادلة التفاضلية الجزئية، مما يؤدي إلى تشكل موجة صدمية . يعتمد تقاطع الخصائص من عدمه على الشرط الابتدائي. في الواقع، يُعطى زمن التكسر قبل تشكل الموجة الصدمية بالمعادلتين [ 9 ] و[ 10 ].

تب=-1معلوماتx(و(x)).{\displaystyle t_{b}={\frac {-1}{\inf _{x}\left(f^{\prime }(x)\right)}}.}

التكامل الكامل لمعادلة برجر غير اللزجة

الحل الضمني الموصوف أعلاه يحتوي على دالة اختياريةو{\displaystyle f}يُطلق عليه التكامل العام. مع ذلك، فإن معادلة برجر غير اللزجة، كونها معادلة تفاضلية جزئية من الرتبة الأولى ، لها أيضًا تكامل كامل يحتوي على ثابتين اختياريين (للمتغيرين المستقلين). [ 11 ] وقدّم سوبرامانيان تشاندراسيكار التكامل الكامل في عام 1943، [ 12 ] والذي يُعطى بالصيغة التالية:

u(x،ت)=أx+بأت+1.{\displaystyle u(x,t)={\frac {ax+b}{at+1}}.}

أينأ{\displaystyle a}وب{\displaystyle b}هي ثوابت اختيارية. يحقق التكامل الكامل شرطًا ابتدائيًا خطيًا، أيو(x)=أx+ب{\displaystyle f(x)=ax+b}يمكن أيضًا إنشاء التكامل العام باستخدام التكامل الكامل المذكور أعلاه.

معادلة البرغر اللزج

يمكن تحويل معادلة برجر اللزجة إلى معادلة خطية عن طريق تحويل كول-هوف ، [ 13 ] [ 14 ] [ 15 ]

u(x،ت)=-2νxlnφ(x،ت)،{\displaystyle u(x,t)=-2\nu {\frac {\partial }{\partial x}}\ln \varphi (x,t),}

مما يحولها إلى المعادلة

2νx[1φ(φت-ν2φx2)]=0،{\displaystyle 2\nu {\frac {\partial }{\partial x}}\left[{\frac {1}{\varphi }}\left({\frac {\partial \varphi }{\partial t}}-\nu {\frac {\partial ^{2}\varphi }{\partial x^{2}}}\right)\right]=0,}

والتي يمكن دمجها فيما يتعلق بـx{\displaystyle x}للحصول على

φت-ν2φx2=φدو(ت)دت،{\displaystyle {\frac {\partial \varphi }{\partial t}}-\nu {\frac {\partial ^{2}\varphi }{\partial x^{2}}}=\varphi {\frac {df(t)}{dt}},}

أيندو/دت{\displaystyle df/dt}هي دالة اعتباطية للزمن. بإدخال التحويلφφهـو{\displaystyle \varphi \to \varphi e^{f}}(وهذا لا يؤثر على الوظيفة)u(x،ت){\displaystyle u(x,t)})، وتختزل المعادلة المطلوبة إلى معادلة الحرارة [ 16 ]

φت=ν2φx2.{\displaystyle {\frac {\partial \varphi }{\partial t}}=\nu {\frac {\partial ^{2}\varphi }{\partial x^{2}}}.}

يمكن حل معادلة الانتشار . أي، إذاφ(x،0)=φ0(x){\displaystyle \varphi (x,0)=\varphi _{0}(x)}، ثم

φ(x،ت)=14πνت-φ0(x)خبرة[-(x-x)24νت]دx.{\displaystyle \varphi (x,t)={\frac {1}{\sqrt {4\pi \nu t}}}\int _{-\infty }^{\infty }\varphi _{0}(x')\exp \left[-{\frac {(x-x')^{2}}{4\nu t}}\right]dx'.}

الدالة الأوليةφ0(x){\displaystyle \varphi _{0}(x)}يرتبط بالوظيفة الأوليةu(x،0)=و(x){\displaystyle u(x,0)=f(x)}بواسطة

lnφ0(x)=-12ν0xو(x)دx،{\displaystyle \ln \varphi _{0}(x)=-{\frac {1}{2\nu }}\int _{0}^{x}f(x')dx',}

حيث يتم اختيار الحد الأدنى بشكل تعسفي. وبعكس تحويل كول-هوبف، نحصل على

u(x،ت)=-2νxln{14πνت-خبرة[-(x-x)24νت-12ν0xو(x")دx"]دx}{\displaystyle u(x,t)=-2\nu {\frac {\partial }{\partial x}}\ln \left\{{\frac {1}{\sqrt {4\pi \nu t}}}\int _{-\infty }^{\infty }\exp \left[-{\frac {(x-x')^{2}}{4\nu t}}-{\frac {1}{2\nu }}\int _{0}^{x'}f(x'')dx''\right]dx'\right\}}

مما يُبسط، عن طريق التخلص من العامل الزمني في وسيط اللوغاريتم ، إلى

u(x،ت)=-2νxln{-خبرة[-(x-x)24νت-12ν0xو(x")دx"]دx}.{\displaystyle u(x,t)=-2\nu {\frac {\partial }{\partial x}}\ln \left\{\int _{-\infty }^{\infty }\exp \left[-{\frac {(x-x')^{2}}{4\nu t}}-{\frac {1}{2\nu }}\int _{0}^{x'}f(x'')dx''\right]dx'\right\}.}

هذا الحل مشتق من حل معادلة الحرارة لـφ{\displaystyle \varphi }يتلاشى ذلك إلى الصفر عندماx±{\displaystyle x\to \pm \infty }حلول أخرى لـu{\displaystyle u}يمكن الحصول عليها انطلاقاً من حلول لـφ{\displaystyle \varphi }التي تستوفي شروطًا حدودية مختلفة.

بعض الحلول الصريحة لمعادلة برجر اللزجة

تتوفر تعابير صريحة لمعادلة برجر اللزجة. ترد أدناه بعض الحلول ذات الصلة الفيزيائية: [ 17 ]

موجة متحركة تنتشر بثبات

لوu(x،0)=و(x){\displaystyle u(x,0)=f(x)}بحيثو(-)=و+{\displaystyle f(-\infty )=f^{+}}وو(+)=و-{\displaystyle f(+\infty )=f^{-}}وو(x)<0{\displaystyle f'(x)<0}إذن لدينا حل الموجة المتحركة (بسرعة ثابتة)ج=(و++و-)/2{\displaystyle c=(f^{+}+f^{-})/2}) معطى بواسطة

u(x،ت)=ج-و+-و-2tanh[و+-و-4ν(x-جت)].{\displaystyle u(x,t)=c-{\frac {f^{+}-f^{-}}{2}}\tanh \left[{\frac {f^{+}-f^{-}}{4\nu }}(x-ct)\right].}

يُستخدم هذا الحل، الذي اشتقه هاري بيتمان في الأصل عام 1915، [ 6 ] لوصف تغير الضغط عبر موجة صدمية ضعيفة [ 16 ] . عندماو+=2{\displaystyle f^{+}=2}وو-=0{\displaystyle f^{-}=0}وهذا يتبسط إلى

u(x،ت)=21+هـx-تν{\displaystyle u(x,t)={\frac {2}{1+e^{\frac {x-t}{\nu }}}}}

معج=1{\displaystyle c=1}.

دالة دلتا كشرط ابتدائي

لوu(x،0)=2νRهـدلتا(x){\displaystyle u(x,0)=2\nu \mathrm {Re} \delta (x)}، أينRهـ{\displaystyle \mathrm {Re} }(على سبيل المثال، رقم رينولدز ) ثابت، ثم لدينا [ 18 ]

u(x،ت)=νπت[(هـRهـ-1)هـ-x2/4νت1+(هـRهـ-1)هـروج(x/4νت)/2].{\displaystyle u(x,t)={\sqrt {\frac {\nu }{\pi t}}}\left[{\frac {(e^{\mathrm {Re} }-1)e^{-x^{2}/4\nu t}}{1+(e^{\mathrm {Re} }-1)\mathrm {erfc} (x/{\sqrt {4\nu t}})/2}}\right].}

في الحدRهـ0{\displaystyle \mathrm {Re} \to 0}، يتمثل السلوك المحدود في انتشار المصدر، وبالتالي يُعطى بواسطة

u(x،ت)=2νRهـ4πνتخبرة(-x24νت).{\displaystyle u(x,t)={\frac {2\nu \mathrm {Re} }{\sqrt {4\pi \nu t}}}\exp \left(-{\frac {x^{2}}{4\nu t}}\right).}

من ناحية أخرى، في الحدRهـ{\displaystyle \mathrm {Re} \to \infty }، يقترب الحل من حل موجة الصدمة المذكور سابقًا لتشاندراسيكار لمعادلة برجر غير اللزجة، ويُعطى بواسطة

u(x،ت)={xت،0<x<2νRهـت،0،خلاف ذلك.{\displaystyle u(x,t)={\begin{cases}{\frac {x}{t}},\quad 0<x<{\sqrt {2\nu \mathrm {Re} \,t}},\\0,\quad {\text{otherwise}}.\end{cases}}}

يُحدد موقع الموجة الصدمية وسرعتها بواسطةx=2νRهـت{\displaystyle x={\sqrt {2\nu \mathrm {Re} \,t}}}وνRهـ/ت.{\displaystyle {\sqrt {\nu \mathrm {Re} /t}}.}

حل الموجة N

يتألف حل الموجة N من موجة انضغاط تليها موجة تخلخل . ويُعطى حل من هذا النوع بالصيغة التالية:

u(x،ت)=xت[1+1هـRهـ0-1تت0خبرة(-Rهـ(ت)x24νRهـ0ت)]-1{\displaystyle u(x,t)={\frac {x}{t}}\left[1+{\frac {1}{e^{\mathrm {Re} _{0}-1}}}{\sqrt {\frac {t}{t_{0}}}}\exp \left(-{\frac {\mathrm {Re} (t)x^{2}}{4\nu \mathrm {Re} _{0}t}}\right)\right]^{-1}}

أينRهـ0{\displaystyle \mathrm {Re} _{0}}يمكن اعتبارها رقم رينولدز الأولي في وقتت=ت0{\displaystyle t=t_{0}}وRهـ(ت)=(1/2ν)0uدx=ln(1+τ/ت){\displaystyle \mathrm {Re} (t)=(1/2\nu )\int _{0}^{\infty }udx=\ln(1+{\sqrt {\tau /t}})}معτ=ت0هـRهـ0-1{\displaystyle \tau =t_{0}{\sqrt {e^{\mathrm {Re} _{0}}-1}}}، يمكن اعتبارها رقم رينولدز المتغير مع الزمن.

أشكال أخرى

معادلة برجر متعددة الأبعاد

في بعدين أو أكثر، تصبح معادلة برجر

uت+uu=ν2u.{\displaystyle {\frac {\partial u}{\partial t}}+u\cdot \nabla u=\nu \nabla ^{2}u.}

يمكن أيضًا توسيع المعادلة لتشمل حقل المتجهاتu{\displaystyle \mathbf {u} }كما في

uت+uu=ν2u.{\displaystyle {\frac {\partial \mathbf {u} }{\partial t}}+\mathbf {u} \cdot \nabla \mathbf {u} =\nu \nabla ^{2}\mathbf {u} .}

معادلة برجر المعممة

تُوسّع معادلة برجر المعممة الحمل الحراري شبه الخطي إلى شكل أكثر عمومية، أي

uت+ج(u)ux=ν2ux2.{\displaystyle {\frac {\partial u}{\partial t}}+c(u){\frac {\partial u}{\partial x}}=\nu {\frac {\partial ^{2}u}{\partial x^{2}}}.}

أينج(u){\displaystyle c(u)}هي أي دالة اختيارية لـ u. غير لزجν=0{\displaystyle \nu =0}لا تزال المعادلة معادلة زائدية شبه خطية لـج(u)>0{\displaystyle c(u)>0}ويمكن بناء حلها باستخدام طريقة الخصائص كما في السابق. [ 19 ]

معادلة برجر العشوائية

ضوضاء إضافية في الزمكانη(x،ت)=دبليو˙(x،ت){\displaystyle \eta (x,t)={\dot {W}}(x,t)}، أيندبليو{\displaystyle W}هول2(R){\displaystyle L^{2}(\mathbb {R} )}تشكل عملية وينر معادلة برجر العشوائية [ 20 ]

uت+uux=ν2ux2-ληx.{\displaystyle {\frac {\partial u}{\partial t}}+u{\frac {\partial u}{\partial x}}=\nu {\frac {\partial ^{2}u}{\partial x^{2}}}-\lambda {\frac {\partial \eta }{\partial x}}.}

هذه المعادلة التفاضلية الجزئية العشوائية هي النسخة أحادية البعد من معادلة كاردار-باريسي-تشانغ في مجالح(x،ت){\displaystyle h(x,t)}عند الاستبدالu(x،ت)=-λح/x{\displaystyle u(x,t)=-\lambda \partial h/\partial x}.

انظر أيضاً

مراجع

  1. ميسرا، سورين؛ راغوراما راو، إس في؛ بوبا، مانوج كومار (2010-09-01). "نمذجة النطاق الفرعي للشبكة القائمة على نظام الاسترخاء لمحاكاة الدوامات الكبيرة لمعادلة برجر" . المجلة الدولية لديناميكا الموائع الحسابية . 24 (8): 303-315 . Bibcode : 2010IJCFD..24..303M . doi : 10.1080/10618562.2010.523518 . ISSN 1061-8562 . S2CID 123001189 .  
  2. يتعلق الأمر بمعادلة نافيير-ستوكس للزخم مع حذف حد الضغط، معادلة برجر (PDF)المتغير هنا هو سرعة التدفق y = u
  3. ينشأ هذا من معادلة ويسترفيلت بافتراض انتشار الموجات للأمام فقط، واستخدام تحويل إحداثيات إلى إطار زمني متأخر: هنا المتغير هو الضغط
  4. موشا، توشيميتسو؛ هيغوتشي، هيديو (1978-05-01). "تذبذب تيار المرور ومعادلة برغر" . المجلة اليابانية للفيزياء التطبيقية . 17 (5): 811. Bibcode : 1978JaJAP..17..811M . doi : 10.1143/JJAP.17.811 . ISSN 1347-4065 . S2CID 121252757 .  
  5. ينشأ هذا كمعادلة تفاضلية تصف مستويات الطاقة لنظرية المجال الكمومي المشوهة بفعل تأثيرتيتي¯(z){\displaystyle T{\bar {T}}(z)}المشغل. زامولوذيكوف، ألكسندر ب. (2004)، القيمة المتوقعة للحقل المركبتيتي¯{\displaystyle T{\bar {T}}}في نظرية الحقل الكمومي ثنائية الأبعاد ، arXiv : hep-th/0401146، كافاجليا، أندريا؛ نيغرو، ستيفانو؛ سيكسيني، إستفان م.؛ تاتيو، روبرتو (2016). "تيتي¯{\displaystyle T{\bar {T}}}"نظريات الحقول الكمومية ثنائية الأبعاد المشوهة". JHEP . 10 : 112. arXiv : 1608.05534 . doi : 10.1007/JHEP10(2016)112 .
  6. 1 2 باتمان، هـ. (1915). "بعض الأبحاث الحديثة حول حركة السوائل" . مجلة الطقس الشهرية . 43 (4): 163-170 . Bibcode : 1915MWRv...43..163B . doi : 10.1175/1520-0493(1915)43 < 163:SRROTM > 2.0.CO ; 2 .
  7. ويذام، جي بي (2011). الموجات الخطية وغير الخطية (المجلد 42). جون وايلي وأولاده.
  8. بورغرز، ج. م. (1948). "نموذج رياضي يوضح نظرية الاضطراب". التقدم في الميكانيكا التطبيقية . 1 : 171-199 . doi : 10.1016/S0065-2156(08)70100-5 . ISBN 9780123745798.{{cite journal}}عدم توافق رقم ISBN / التاريخ ( مساعدة )
  9. أولفر، بيتر ج. (2013). مقدمة في المعادلات التفاضلية الجزئية . نصوص جامعية في الرياضيات. عبر الإنترنت: سبرينغر. ص 37. doi : 10.1007/978-3-319-02099-0 . ISBN  978-3-319-02098-3. S2CID 220617008 . 
  10. كاميرون، ماريا (29 فبراير 2024). "ملاحظات حول معادلة برجر" (ملف PDF) . قسم الرياضيات بجامعة ميريلاند، الموقع الإلكتروني الشخصي لماريا كاميرون . تم الاطلاع عليه بتاريخ 29 فبراير 2024 .
  11. فورسيث، أ. ر. (1903). رسالة في المعادلات التفاضلية . لندن: ماكميلان.
  12. تشاندراسيكار، س. (1943). حول اضمحلال موجات الصدمة المستوية (تقرير). مختبرات أبحاث المقذوفات. التقرير رقم 423.
  13. كول، جوليان (1951). "حول معادلة قطع مكافئ شبه خطية تظهر في الديناميكا الهوائية" . مجلة الرياضيات التطبيقية الفصلية . 9 (3): 225-236 . doi : 10.1090/qam/42889 . JSTOR 43633894 . 
  14. إيبرهارد هوبف (سبتمبر 1950). "المعادلة التفاضلية الجزئية u t + uu x = μu xx ". مجلة الاتصالات في الرياضيات البحتة والتطبيقية . 3 (3): 201-230 . doi : 10.1002/cpa.3160030302 .
  15. كيفوركيان، ج. (1990). المعادلات التفاضلية الجزئية: تقنيات الحل التحليلي . بلمونت: وادزورث. ص 31-35 . ISBN  0-534-12216-7.
  16. 1 2 لاندو، إل دي، وليفشيتز، إي إم (2013). ميكانيكا الموائع: لاندو وليفشيتز: دورة في الفيزياء النظرية، المجلد 6 (المجلد 6). إلسيفير. الصفحات 352-354.
  17. صالح، أ. "معادلة البرجر". المعهد الهندي لعلوم وتكنولوجيا الفضاء، ثيروفانانثابورام (2016).
  18. ويذام، جيرالد بيريسفورد. الموجات الخطية وغير الخطية. جون وايلي وأولاده، 2011.
  19. كوران، ر.، وهيلبرت، د. أساليب الفيزياء الرياضية. المجلد الثاني.
  20. ↑ وانغ، و.؛ روبرتس، أ. ج . (2015). "تقريب الانتشار للتشابه الذاتي للحمل العشوائي في معادلة برغر". الاتصالات في الفيزياء الرياضية . 333 (3): 1287-1316 . arXiv : 1203.0463 . Bibcode : 2015CMaPh.333.1287W . doi : 10.1007/s00220-014-2117-7 . S2CID 119650369 .