قانوني فيك للانتشار (بالإنجليزية: Fick's laws of diffusion)‏ عبارة عن قانونين يمثلان انتشار المواد عند غشاء نافذ ويحسبان تدفق الكتل عبر الغشاء، كما ويمثلان الانتشار الذاتي لمادة في محلول أو خليط ما.^[1] ويسمى هذان القانونان باسم العالم الألماني أدولف فيك الذي اشتقهما عام 1855.

يربط هذه القانون بين تدفق الانتشار مادّة في وسط ما وحقل التركيز لهذه المادّة. الفرضية الكامنة في هذا القانون هو أنّ التدفق يكون من المناطق ذات التركيز المرتفع في الحقل إلى المناطق ذات التركيز المنخفض، بمقدار يتناسب مع تدرج التراكيز في الحقل. لذا، فإذا افترضنا أنّ التدرج موجود فقط على امتداد محور x (أي أنّه لا تغيير في التركيز في المحورين المعامدين)، فإنّ تدفق الانتشار يحقّق:

${\displaystyle {\bigg .}J=-D{\frac {\partial \phi }{\partial x}}{\bigg .}}$

بحيث أن:

J هو تدفق الانتشار، بوحدات ${\displaystyle ({\tfrac {\mathrm {mol} }{m^{2}\cdot s}})}$، وهو يقيس كميّة المادّة المتدفّقة عبر مساحة صغيرة خلال فترة قصيرة من الزمن؛

D هو معامل الانتشار، أو الانتشارية، بوحدات ${\displaystyle ({\tfrac {m^{2}}{s}})}$ وهو يتعلّق بنوع المادة، ويعبّر عن المساحة التي تنتشر إليها مادّة معيّنة خلال فترة زمنية قصيرة؛

${\displaystyle \,\phi }$ هو تركيز المادّة، بوحدات ${\displaystyle ({\tfrac {\mathrm {mol} }{m^{3}}})}$؛

x هو إحداثي ويدل على الموقع الذي يتم به حساب التدفق، بوحدات (${\displaystyle \,m}$).

يجدر الذكر أنّ المعامل ${\displaystyle D}$ هو الوحيد في قانون فيك الذي يمت بصلة لنوع المادّة، وهو يتناسب طرديًا مع مربّع سرعة جزيئات المادّة، أي أنّه يتعلّق بدرجة الحرارة بلزوجة الوسط الذي تتحرّك فيه الجزيئات وبحجم الجزيئات، وفق علاقة ستوكس-أينشتاين. بالإضافة إلى ذلك، فإنّ إشارة الناقص تشير إلى أنّ اتجاه التدفق يكون عكس اتجاه التدرج؛ فمثلاً، إذا ارتفع تركيز المادة على امتداد الاتجاه التصاعدي لـ x، فإنّ التدفق سيكون بالاتجاه التنازلي لقيم x.

كذلك، يذكر أنّه بعدم وجود أية قوى أخرى تؤثر على المادّة (بما في ذلك قوى تكون بين جزيئات المادّة نفسها)، فإنّ كل جزيء للمادّة سيتحرك بشكل عشوائي في الفراغ. لكن بسبب وجود تدرج تراكيز، فمن ناحية إحصائية سيكون هنالك عدد أكبر من جزيئات المادة الموجودة في أماكن التركيز المرتفع التي تتحرك نحو أماكن التركيز المنخفض؛ أي أنّه، بالمجمل، سوف تشهد حركة جزيئات نحو أماكن التركيز المنخفض، مع أنّ هذا لا ينفي تحرك جزء من الجزيئات من أماكن التركيز المنخفض نحو أماكن ذات تركيز أكثر ارتفاعًا بالذات. وهذا هو بالضبط المقصود من قانون فيك الأوّل للانتشار - فإنّه ينظر إلى حقل التركيز، وليس إلى حركة كل جزيء على حدة: بالمجمل، هنالك تدفق لجزيئات المادّة بمقدار يتناسب مع تدرج حقل التراكيز.

أمّا إذا كان حقل التراكيز متغيّرًا بأبعاده الثلاثة، فإنّ تدفق الانتشار في هذه الحالة يحسب وفق تعميم رياضي للقانون السابق:

${\displaystyle {\bigg .}J=-D\nabla \phi {\bigg .}}$

بحيث أنّ ${\displaystyle \nabla \phi }$ هو التدرج ثلاثي الأبعاد لحقل التراكيز: ${\displaystyle \nabla \phi ={\frac {\partial \phi }{\partial x}}+{\frac {\partial \phi }{\partial y}}+{\frac {\partial \phi }{\partial z}}}$.

أخيرًا، يذكر هنا أنّ القوانين أعلاه، حيث القوة الدافعة للانتشار هي تدرج التركيز، صحيحة فقط لمحلول أو خليط مثالي. أمّا في منظومات كيميائيّة أخرى، فإنّ القوة الدافعة لانتشار كل مادة في الخليط هي تدرّج الكمون الكيميائي لتلك المادّة. هنا، يتحوّل قانون فيك الأوّل (بصيغته أحادية البعد) إلى:

${\displaystyle {\bigg .}J=-{\frac {Dc_{i}}{RT}}{\frac {\partial \mu _{i}}{\partial x}}{\bigg .}}$

بحيث أنّ:

${\displaystyle c_{i}}$ هو تركيز المادّة i بوحدات ${\displaystyle {\tfrac {mol}{m^{3}}}}$؛

${\displaystyle R}$ هو ثابت الغازات العام بوحدات ${\displaystyle ({\tfrac {J}{Kmol}})}$؛

${\displaystyle T}$ هي درجة الحرارة المطلقة بوحدات ${\displaystyle K}$؛

${\displaystyle \mu }$ هو الجهد الكيميائي (الكمون الكيميائي) للمادّة i بوحدات ${\displaystyle ({\tfrac {J}{mol}})}$.

إذا ما استمرّت عملية الانتشار فترة كافية، أو إذا كانت عملية الانتشار سريعة، فإنّها ستؤدي في نهاية المطاف إلى تغيير حقل التركيز. قانون فيك الثاني يعطي تصورًا لتلك التغييرات التي تطرأ في حقل التركيز نتيجة عملية الانتشار:

${\displaystyle {\bigg .}{\frac {\partial \phi }{\partial t}}=D{\frac {\partial ^{2}\phi }{\partial x^{2}}}{\bigg .}}$

بحيث أنّ t هو الاحداثي الزمني (بوحدات الثانية، مثلاً)، وباقي الرموز هي كالسابق.

لنفرض أنّ التركيز ${\displaystyle \phi }$ يتغيّر فقط باتجاه محور x. إذا ما طبّقنا قانون حفظ المادة على حجم موحل في الصغر ذو مساحة ${\displaystyle A}$ وسمك لامتناهٍ في الصغر ${\displaystyle dx}$، كما في الرسم، نحصل على المعادلة الآتية:

${\displaystyle \left(\phi (t+dt)-\phi (t)\right)Adx=\left(J_{in}-J_{out}\right)Adt}$

أو:

${\displaystyle \left(\phi (t+dt)-\phi (t)\right)Adx=\left(J(x+dx)-J(x)\right)Adt}$

الجهة اليسرى من المعادلة تعطينا الكمية الكلية للمادة في الحجم الذي نتمعّن به حاليًا. أمّا الجهة اليمنى، فهي عبارة عن الفرق بين كميّة المادّة التي تتدفّق بين جهتي الحجم (السطح الأخضر في الرسم والسطح المقابل له). حسب قانون حفظ المادّة، وإذا فرضنا أنّه بداخل الحجم أعلاه لا تتكوّن أو تختفي أي جزيئات من المادة (أي لا يوجد مصدر أو مصرف للمادة في الحجم)، فإنّ الفرق بين الكميات المتدفقة من الجهتين يساوي الكمّيّة الموجودة في داخل الحجم بكل لحظة.

من الواضح أنّه باختزال A، وبتعديلات جبرية طفيفة، بالإمكان إحالة المعادلة إلى الصورة التالية:

${\displaystyle {\frac {\partial \phi }{\partial x}}=-{\frac {\partial J}{\partial x}}}$

لكن العلاقة بين تدفق الانتشار وبين تركيز المادة معروفة من قانون فيك الأوّل:

${\displaystyle {\frac {\partial \phi }{\partial x}}=-{\frac {\partial J}{\partial x}}={\frac {\partial }{\partial x}}\left(D{\frac {\partial \phi }{\partial x}}\right)=D{\frac {\partial ^{2}\phi }{\partial x^{2}}}}$

بالخطوة الأخيرة فرضنا أن معامل الانتشار، D، لا يتغير مع x، ولذلك بالإمكان إخراجه من المشتقة الجزئية والحصول على قانون فيك الثاني.

قانون فيك الثاني العام (ثلاثي الأبعاد) يكتب بالصورة التالية:

${\displaystyle {\frac {\partial \phi }{\partial t}}=D\nabla ^{2}\phi }$

وهو معروف أيضًا باسم معادلة الحرارة، بحيث أنّ:

${\displaystyle \nabla ^{2}\phi ={\frac {\partial ^{2}\phi }{\partial x^{2}}}+{\frac {\partial ^{2}\phi }{\partial y^{2}}}+{\frac {\partial ^{2}\phi }{\partial z^{2}}}}$

أما إذا كان معامل الانتشار يتعلق بالاحداثيات المكانيّة فيكون قانون فيك الثاني على الصورة الآتية:

${\displaystyle {\frac {\partial \phi }{\partial t}}=\nabla \left(D\nabla \phi \right)}$

• إحدى الحالات التي من المهم معاينتها هي الحالة الراسخة التي لا يتغيّر فيها حقل التراكيز مع الزمن، أي أنّ المشتقة الجزئية بالمتغير t تساوي صفر. في هذه الحالة نحصل من قانون فيك الثاني على معادلة لابلاس:

${\displaystyle \nabla ^{2}\phi =0}$

ويطلق الرياضيون على حلول هذه المعادلة (التي تتعلّق بالشروط الحدودية للمسألة) اسم الدوال التوافقية.

• إحدى حالات الانتشار التي من السهل حلّها، نسبيًا، هي انتشار في وسط نصف لا-نهائي، يتغير فيه تركيز المادة المنتشرة باتجاه واحد فقط (انظر الرسم). الشرط الحدودي في المستوى ${\displaystyle x=0}$ يحدّد نوعيّة الحل الناتج لمعادلة فيك.

لأنّ قانون فيك الثاني أحادي البعد يحتوي على ثلاثة عمليات اشتقاق (اثنتان بالمتغيّر x وواحدة بالمتغير t)، سوف نحتاج إلى ثلاثة شروط حدوديّة لحل المسألة بكاملها. إحدى الحالات البسيطة هي:

1. ${\displaystyle n(x,t=0)=n_{i}}$، ويسمى أيضًا شرط ابتدائي أو أوّلي، ويمثل كون المادة بتركيز أوّلي متساوٍ في كل مكان في اللحظة ${\displaystyle t=0}$؛
2. ${\displaystyle n(x\to \infty ,t)=n_{i}}$، وهو يعبّر عن كون تركيز المادّة في النقاط البعيدة بعدًا هائلاً عن الحد ${\displaystyle x=0}$ مساويًا دائمًا للتركيز الأوّلي (أي لكل t) ولن يتأثر أبدًا بعمليّة الانتشار بسبب البعد الزمني والمكاني.

أمّا بالنسبة للشرط الحدودي الثالث والأهم، في المستوى ${\displaystyle x=0}$، فسنمعن النظر في شرطين مختلفين:

ينص الشرط الأوّل على أنّ التركيز على طول المستوى ${\displaystyle x=0}$ يبقى ثابتًا بكل نقطة زمنيّة، أيّ:

${\displaystyle n(x=0,t)=n_{0}}$

في هذه الحالة فإنّ الحل الناتج سيكون:

${\displaystyle {\frac {n(x,t)-n_{i}}{n_{0}-n_{i}}}=1-{\mbox{erf}}({\frac {x}{2{\sqrt {Dt}}}})={\mbox{erfc}}({\frac {x}{2{\sqrt {Dt}}}})}$

بحيث أنّ ${\displaystyle {\mbox{erf}}}$ هي دالة تابع الارتياب، تحقّق:

${\displaystyle {\mbox{erf}}(x)={\frac {2}{\sqrt {\pi }}}\int _{0}^{x}e^{-t^{2}}\,dt}$

و${\displaystyle {\mbox{erfc}}}$ هي الدالة المتممة لها.

الشرط الثاني ينص على أنّ تدفّق الانتشار في المستوى ${\displaystyle x=0}$، والذي نحصل عليه من قانون فيك الأوّل، هو ثابت مع الزمن:

${\displaystyle -D{\frac {\partial n}{\partial x}}{\Big |}_{x=0}=j_{0}}$

في هذه الحالة، يكون حل الانتشار هو:

${\displaystyle n(x,t)-n_{i}=2j_{0}{\sqrt {\frac {t}{\pi D}}}e^{-{\frac {x^{2}}{4Dt}}}-{\frac {j_{0}x}{D}}{\mbox{erfc}}\left({\frac {x}{2{\sqrt {Dt}}}}\right)}$

يشار إلى أنّ الكميّة ${\displaystyle 2{\sqrt {Dt}}}$ ذات وحدات مسافة وتسمّى مسافة الانتشار، لكونها تظهر في الحلول أعلاه دائمًا في مقام كسر بسطه هو ${\displaystyle x}$، وهي تعطي مقياسًا لمدى انتشار جزيئات المادة في زمن t.

• انتشار
• معامل الانتشار
• تناضح
• تبادل الغازات
• رئة

• أسس الانتشار