تجزئة معادلات نافيير-ستوكس

يُعدّ تجزئة معادلات نافيير -ستوكس لديناميكا الموائع إعادة صياغة للمعادلات بطريقة تُمكّن من تطبيقها على ديناميكا الموائع الحسابية . ويمكن تطبيق عدة طرق للتجزئة:

طريقة الحجم المحدود

التدفق غير القابل للانضغاط

نبدأ بالصيغة غير القابلة للانضغاط لمعادلة الزخم. تم تقسيم المعادلة على الكثافة ( P = p / ρ ) وتم دمج الكثافة في حد قوة الجسم . uأنات+uأناuجxج=-Pxأنا+ν2uأناxجxج+وأنا{\displaystyle {\frac {\partial u_{i}}{\partial t}}+{\frac {\partial u_{i}u_{j}}{\partial x_{j}}}=-{\frac {\partial P}{\partial x_{i}}}+\nu {\frac {\partial ^{2}u_{i}}{\partial x_{j}\partial x_{j}}}+f_{i}}

يتم تكامل المعادلة على حجم التحكم لخلية حسابية. V[uأنات+uأناuجxج]دV=V[-Pxأنا+ν2uأناxجxج+وأنا]دV\displaystyle \iiint _{V}\left[{\frac {\partial u_{i}}{\partial t}}+{\frac {\partial u_{i}u_{j}}{\partial x_{j}}}\right]dV=\iiint _{V}\left[-{\frac {\partial P}{\partial x_{i}}}+\nu {\frac {\partial ^{2}u_{i}}{\partial x_{j}\partial x_{j}}}+f_{i}\right]dV}

يُفترض ثبات الحدّ المتغير مع الزمن وحدّ قوة الجسم على كامل حجم الخلية. وتُطبّق نظرية التباعد على حدود الحمل ، وتدرج الضغط ، والانتشار. uأناتV+أuأناuجنجدأ=-أPنأنادأ+أνuأناxجنجدأ+وأناV// حيث n هو متجه الوحدة العمودي على سطح حجم التحكم، و V هو الحجم. إذا كان حجم التحكم متعدد السطوح، وبافتراض ثبات القيم على كل وجه، فيمكن كتابة تكاملات المساحة كمجموع على كل وجه. uأناتV+نبر(uأناuجنجأ)نبر=-نبر(Pنأناأ)نبر+نبر(νuأناxجنجأ)نبر+وأناV{\displaystyle {\frac {\partial u_{i}}{\partial t}}V+\sum _{nbr}\left(u_{i}u_{j}n_{j}A\right)_{nbr}=-\sum _{nbr}\left(Pn_{i}A\right)_{nbr}+\sum _{nbr}\left(\nu {\frac {\partial u_{i}}{\partial x_{j}}}n_{j}A\right)_{nbr}+f_{i}V} حيث يشير الرمز السفلي nbr إلى القيمة عند أي وجه معين.

شبكة ديكارتية ثنائية الأبعاد ذات مسافات منتظمة

بالنسبة لشبكة ديكارتية ثنائية الأبعاد، يمكن توسيع المعادلة إلى uأناتΔxΔy-(uأناuΔy)w+(uأناuΔy)هـ-(uأناvΔx)s+(uأناvΔx)ن=-(PنأناΔy)w-(PنأناΔy)هـ-(PنأناΔx)s-(PنأناΔx)ن-(νuأناxΔy)w+(νuأناxΔy)هـ-(νuأناyΔx)s+(νuأناyΔx)ن+وأنا\begin{aligned}\frac{\partial u_{i}}{\partial t}\Delta x\Delta y-\left(u_{i}u\Delta y\right)_{w}+\left(u_{i}u\Delta y\right)_{e}-\left(u_{i}v\Delta x\right)_{s}+\left(u_{i}v\Delta x\right)_{n}=\\&-\left(Pn_{i}\Delta y\right)_{w}-\left(Pn_{i}\Delta y\right)_{e}-\left(Pn_{i}\Delta x\right)_{s}-\left(Pn_{i}\Delta x\right)_{n}\\&-\left(\nu {\frac {\partial u_{i}}{\partial x}}\Delta y\right)_{w}+\left(\nu {\frac {\partial u_{i}}{\partial x}}\Delta y\right)_{e}-\left(\nu {\frac {\partial u_{i}}{\partial y}}\Delta x\right)_{s}+\left(\nu {\frac {\partial u_{i}}{\partial y}}\Delta x\right)_{n}+f_{i}\end{aligned}}}

على شبكة متداخلة ، تكون معادلة الزخم في الاتجاه السيني هي uتΔxΔy-(uuΔy)w+(uuΔy)هـ-(uvΔx)s+(uvΔx)ن=+(PΔy)w-(PΔy)هـ-(νuxΔy)w+(νuxΔy)هـ-(νuyΔx)s+(νuyΔx)ن+وx{\displaystyle {\begin{aligned}&{\frac {\partial u}{\partial t}}\Delta x\Delta y-\left(uu\Delta y\right)_{w}+\left(uu\Delta y\right)_{e}-\left(uv\Delta x\right)_{s}+\left(uv\Delta x\right)_{n}=\\&+\left(P\Delta y\right)_{w}-\left(P\Delta y\right)_{e}-\left(\nu {\frac {\partial u}{\partial x}}\Delta y\right)_{w}+\left(\nu {\frac {\partial u}{\partial x}}\Delta y\right)_{e}-\left(\nu {\frac {\partial u}{\partial y}}\Delta x\right)_{s}+\left(\nu {\frac {\partial u}{\partial y}}\Delta x\right)_{n}+f_{x}\end{aligned}}} ومعادلة الزخم في الاتجاه y هي vتΔxΔy-(vuΔy)w+(vuΔy)هـ-(vvΔx)s+(vvΔx)ن=+(PΔx)s-(PΔx)ن-(νvxΔy)w+(νvxΔy)هـ-(νvyΔx)s+(νvyΔx)ن+وy{\displaystyle {\begin{aligned}&{\frac {\partial v}{\partial t}}\Delta x\Delta y-\left(vu\Delta y\right)_{w}+\left(vu\Delta y\right)_{e}-\left(vv\Delta x\right)_{s}+\left(vv\Delta x\right)_{n}=\\&+\left(P\Delta x\right)_{s}-\left(P\Delta x\right)_{n}-\left(\nu {\frac {\partial v}{\partial x}}\Delta y\right)_{w}+\left(\nu {\frac {\partial v}{\partial x}}\Delta y\right)_{e}-\left(\nu {\frac {\partial v}{\partial y}}\Delta x\right)_{s}+\left(\nu {\frac {\partial v}{\partial y}}\Delta x\right)_{n}+f_{y}\end{aligned}}}

الهدف في هذه المرحلة هو تحديد تعابير القيم الظاهرية لـ u و v و P ، وتقريب المشتقات باستخدام تقريبات الفروق المحدودة . في هذا المثال، سنستخدم الفروق الخلفية للمشتقة الزمنية والفروق المركزية للمشتقات المكانية. بالنسبة لمعادلتي الزخم، تصبح المشتقة الزمنية كما يلي: uأنات=uأنان-uأنان-1Δت{\displaystyle {\frac {\partial u_{i}}{\partial t}}={\frac {u_{i}^{n}-u_{i}^{n-1}}{\Delta t}}} حيث n هو مؤشر الوقت الحالي وΔt هي خطوة الزمن. كمثال على المشتقات المكانية، تصبح المشتقة في حد الانتشار على الواجهة الغربية في معادلة زخم x كما يلي: (ux)w=uأنا،ج-uأنا-1،جΔx{\displaystyle \left({\frac {\partial u}{\partial x}}\right)_{w}={\frac {u_{I,J}-u_{I-1,J}}{\Delta x}}} حيث يمثل I و J مؤشري خلية الزخم x ذات الأهمية.

طريقة العناصر المحدودة

طريقة الفروق المحدودة