Hiện tượng vận chuyển

Sự lan truyền của các phân tử của một hỗn hợp ở trong chất lỏng hay chất khí từ nơi chúng đưa vào, trong điều kiện không có các dịch chuyển vĩ mô ở trong chất lỏng hay chất khí, được gọi là *hiện tượng khuếch tán*. Sự xuất hiện các lực ma sát giữa hai lớp chất lỏng hay chất khí chuyển động với các vận tốc khác nhau là biểu hiện của *lực ma sát nội* hay *tính nhớt*. Sự chuyển của năng lượng từ các miền nóng hơn sang các miền lạnh hơn, khi không có các sự dịch chuyển của chất lỏng hay chất khí hay các sự chảy đối lưu trong chúng được gọi là *sự dẫn nhiệt*. Cơ chế bên trong của ba hiện tượng đó là một - đó là chuyển động nhiệt hỗn loạn của các phân tử đã dẫn đến sự dịch chuyển của chúng. Các hiện tượng đó được gọi là các **hiện tượng vận chuyển** vì trong sự khuếch tán có sự vận chuyển chất của hỗn hợp, trong sự ma sát nội có sự vận chuyển xung lượng, còn trong sự dẫn nhiệt có sự vận chuyển nhiệt lượng.

## Lý thuyết

### Đoạn đường tự do của các phân tử
Đoạn đường tự do trung bình của các phân tử là khoảng cách trung bình mà phân tử đi qua giữa các va chạm kế tiếp của nó với các phân tử khác. Trong một đơn vị thời gian, phân tử được tách ra va chạm với tất cả các phân tử, mà tâm của chúng nằm trong hình trụ gẫy có bán kính bằng đường kính d của phân tử, còn độ dài của nó bằng vận tốc trung bình tương đối \bar{u}_{td} của phân tử. Thể tích của hình trụ đó bằng \bar{u}_{td}.\pi.d^2, còn số các phân tử trong đó bằng n.\bar{u}. \pi.d^2. Vì trong mỗi va chạm, có đồng thời hai phân tử tham gia, nên tổng số các va chạm trong một đơn vị thể tích sau một đơn vị thời gian bằng \tfrac{1}{2}\pi.d^2.\bar{u}_{td}.n^2.

Sau một đơn vị thời gian, phân tử va chạm \pi.d^2.n.u_{td} lần trên đoạn đường \bar{u}_{tuyet doi}, do đó đoạn đường tự do trung bình bằng l=\tfrac{u_{tuyetdoi}}{\pi.d^2.n.u_{td}}, hay vì \bar{u}_{td}=\sqrt{2}.\bar{u}_{tuyetdoi} nên l=\tfrac{1}{\sqrt{2}.\pi.d^2.n}

### Hiện tượng khuếch tán
Nếu trong chất khí có phân bố một hỗn hợp với nồng độ n(x), thì khi gradient nồng độ \frac{\operatorname{d}\!n}{\operatorname{d}\!x}\neq0, chuyển động hỗn loạn của các phân tử sẽ tạo điều kiện san bằng nồng độ dọc theo trục Ox.

Dòng các phân tử của hỗn hợp qua diện tích \Delta S đặt vuông góc với trục Ox theo chiều dương bằng \Delta N_+=\tfrac{1}{6}.n_1.\bar{u}.\Delta S.\Delta t, trong đó n_1 là nồng độ của hỗn hợp ở bên trái diện tích, \Delta t là khoảng thời gian.

Dòng các phân tử của hỗn hợp qua cùng diện tích đó theo chiều âm bằng \Delta N_-=\tfrac{1}{6}.n_2.\bar{u}.\Delta S\Delta t, trong đó n_2 là nồng độ các phân tử của hỗn hợp ở bên phải diện tích.

Hệ số \tfrac{1}{6} trong các đẳng thức trên là do tính xác suất của chuyển động hỗn loạn của các phân tử.

Vì sự san bằng các nồng độ xảy ra do các va chạm của những phần tử nên các nồng độ n_1 và n_2 cần được coi là xác định trên các khoảng cách bằng độ dài của đoạn đường tự do l tính từ diện tích \Delta S.

Dòng các phân tử của hỗn hợp qua một đơn vị diện tích trong một đơn vị thời gian bằng:

l=\frac{\Delta N_+-\Delta N_-}{\Delta S\Delta t}=\frac{1}{6}(n_1-n_2)\bar{u}=-\frac{1}{3}\bar{u}l\frac{n_1-n_2}{2l}=-\frac{1}{3}\bar{u}l\tfrac{\Delta n}{\Delta x}

Tích \frac{1}{3}\bar{u}l=D được gọi là *hệ số khuếch tán*.

#### Định luật Fick
Dòng các phân tử của hỗn hợp qua một đơn vị diện tích sau một đơn vị thời gian tỷ lệ thuận với gradient của nồng độ:

I=-D.\frac{\Delta n}{\Delta x}

Dấu trừ chứng tỏ rằng, dòng khuếch tán hướng theo chiều giảm của nồng độ hỗn hợp.

#### Phương trình khuếch tán
Nếu nhân đẳng thức trên với khối lượng của phân tử, ta có thể tìm được định luật cơ bản Fick cho hiện tượng khuếch tán thông qua qua dòng chất khuếch tán q và nồng độ C của chất dưới dạng q=-D\tfrac{\partial c}{\partial x}

Dòng khuếch tán qua một mặt kín S sau khoảng thời gian \Delta t=t_1-t_2 bằng:

q_1=-\int\limits_{t_1}^{t_2}\int\limits\int\limits_{S}D\frac{\partial c}{\partial n}\operatorname{d}\!S\operatorname{d}\!t=-\int\limits_{t_1}^{t_2}\int\limits\int\limits_{S}D grad c.n\operatorname{d}\!S\operatorname{d}\!t, trong đó n là pháp tuyến ngoài với mặt.

Nếu F(x,y,z,t)là mật độ của nguồn chất của hỗn hợp, nghĩa là lượng chất được tạo thành sau một đơn vị thời gian trong một đơn vị thể tích do các phản ứng hóa học, thì trong thể tích V giới hạn bởi mặt S trong khoảng thời gian \Delta t sẽ xuất hiện khối lượng:

q_2=\int\limits_{t_1}^{t_2}\int\limits\int\limits_{V}\int \operatorname F\operatorname{d}\!V\operatorname{d}\!t

Các quá trình đã mô tả có kèm theo sự biến thiên của nồng độ hỗn hợp tại mỗi điểm, sau khoảng thời gian \Delta t, độ biến thiên của lượng chất hỗn hợp trong toàn bộ thể tích V bằng:

q_3=\int\limits_{t_1}^{t_2}\int\limits\int\limits_{V}\int\limits\frac{\partial c}{\partial t}\operatorname{d}\!V\operatorname{d}\!t

Phương trình cân bằng khối lượng có dạng q_2-q_1=q_3. Biến đổi tích phân mặt theo công thức:\int\limits\int\limits_{S}Dgradc.n.\operatorname{d}\!S=\int\int\limits_{V}\int div(Dgradc)\operatorname{d}\!V

Xét đến sự cân bằng các khối lượng, chúng ta thu được:

\int\limits_{t_1}^{t_2}\int\int\limits_{V}\int[F+div(Dgradc)-\frac{\partial c}{\partial t}]\operatorname{d}\!V\operatorname{d}\!t=0

Từ đó ta suy ra được phương trình khuếch tán trong môi trường đứng yên:

\frac{\partial c}{\partial t}=\frac{\partial}{\partial x}(D.\frac{\partial c}{\partial x})+\frac{\partial}{\partial y}(D.\frac{\partial c}{\partial y})+\frac{\partial}{\partial z}(D.\frac{\partial c}{\partial z})+F(x,y,z,t)

Đó là phương trình vi phân với các đạo hàm riêng loại parabolic, tương tự với phương trình dẫn nhiệt.

### Nội ma sát hay tính nhớt
Nếu một lớp khí dịch chuyển đối với một lớp khác, nghĩa là gradient vận tốc theo hướng vuông góc với mặt của lớp khác không thì do sự dời chuyển của các phân tử từ lớp này sang lớp khác trong chuyển động hỗn loạn, sự vận xung lượng từ lớp này sang lớp khác sẽ được thực hiện.

Xung lượng được chuyển qua diện tích \Delta S ở trên mặt phân giới của các lớp theo một hướng bằng: mv_1\Delta N_+=\tfrac{1}{6}mv_1n\bar{u}\Delta S\Delta t, trong đó v_1 là vận tốc dịch chuyển định hướng của chất khí ở trong lớp có các phân tử thoát qua đến những khoảng cách bằng đoạn đường tự do tính từ mặt \Delta S trên.

Xung lượng được chuyển qua theo chiều ngược lại bằng: mv_2\Delta N_-=\tfrac{1}{6}mv_2n\bar{u}\Delta S\Delta t, trong đó v_2 là vận tốc dịch chuyển định hướng của chất khí đến những khoảng cách bằng đoạn đường tự do về phía kia của diện tích.

Độ biến thiên xung lượng của lớp bằng xung lực ma sát:

\Delta F_{ms}\Delta t=mv_1\Delta N_+-mv_2\Delta N_-

Lực ma sát tác dụng lên một đơn vị diện tích của mặt phân giới giữa các lớp lân cận bằng:

F_{ms}=\frac{1}{3}nm\bar{u}l\frac{v_2-v_1}{2l}=-\frac{1}{3}nm\bar{u}l\frac{\Delta v}{\Delta x}

Tích số \tfrac{1}{3}l\bar{u}nm=\eta được gọi là *hệ số ma sát nội* hay *độ nhớt.*

Lực ma sát nội khi các lớp trong chất khí chuyển động tỷ lệ thuận với gradient vận tốc:

F_{ms}=-\eta\frac{\Delta v}{\Delta x}

Công thức Stokes (Gabriel Stokes, 1819 - 1903) đúng cho quả cầu bán kính r chuyển động với vận tốc v trong chất lỏng có độ nhớt \eta, theo công thức này lực cản bằng:

F_{can}=-6\pi\eta r v

### Tính dẫn nhiệt
Cho một dòng nhiệt truyền qua chất khí đựng giữa các thành song song có các nhiệt độ T_1 và T_2 theo phương của trục Ox vuông góc với các thành. Các phân tử đi qua diện tích \Delta S vuông góc với trục mang theo chiều dương của trục một năng lượng bằng:

\omega _1\Delta N_+=\frac{i}{2}kT_1.\frac{1}{6}\bar{u}n\Delta S\Delta t

Năng lượng được mang theo chiều ngược lại bằng:

\omega _2\Delta N_-=\tfrac{i}{2}kT_2.\tfrac{1}{6}\bar{u}n\Delta S\Delta t

Ở đây i là số bậc tự do của các phân tử khí, còn T_1 và T_2 lần lượt là các giá trị nhiệt độ trong các mặt phẳng x-\lambda và x+\lambda. Dòng nhiệt đi qua một đơn vị diện tích trong một đơn vị thời gian bằng:

q=\frac{\omega _1\Delta N_++\omega _2\Delta N_-}{\Delta S\Delta t}=-\frac{1}{3}\bar{u}nl\frac{T_2-T_1}{2l}.\frac{i}{2}k

Đại lượng \lambda=\frac{1}{3}\bar{u}nl\frac{i}{2}k được gọi là *hệ số dẫn nhiệt.*

#### Định luật Fourier
Dòng nhiệt đi qua một đơn vị diện tích trong một đơn vị thời gian tỷ lệ thuận với gradient nhiệt độ:

q=-\lambda\frac{\Delta T}{\Delta x}

Dấu trừ chủng tỏ rằng dòng nhiệt hướng từ lớp có nhiệt độ cao hơn sang lớp có nhiệt độ thấp hơn.

#### Phương trình dẫn nhiệt
Theo định luật Fourier, nhiệt lượng toàn phần đi qua mặt kín S trong khoảng thời gian \Delta t=t_2-t_1 bằng:

q_1=-\int\limits_{t_1}^{t_2}\int\limits_{S}\int\lambda \frac{\partial T}{\partial n}\operatorname{d}\!S\operatorname{d}\!t=-\int\limits_{t_1}^{t_2}\int\limits_{S}\int\limits\lambda gradTn\operatorname{d}\!S\operatorname{d}\!t

Nếu F(x,y,z,t) là mật độ của các nguồn nhiệt thì trong thể tích V sau khoảng thời gian \Delta t, lượng nhiệt tỏa ra bằng:

q_2=\int\limits_{t_1}^{t_2}\iint\limits_V\int\limits F\operatorname{d}\!V\operatorname{d}\!t

Các quá trình nêu trên có kèm theo sự biến thiên nhiệt độ tại mỗi điểm của vật và độ biến thiên toàn phần của lượng nhiệt trong vật sau khoảng thời gian \Delta tbằng:q_3=\int\limits_{t_1}^{t_2}\iint\limits_V\int c_\rho \frac{\delta T}{\delta t}\operatorname{d}\!V\operatorname{d}\!t, trong đó c là nhiệt dung riêng, còn \rho là mật độ.

Nếu dùng phương trình cân bằng nhiệt q_2-q_1=q_3 và biến đổi tích phân mặt theo công thức Ostrogradski:

\iint\limits_SgradT.n\operatorname{d}\!S=\int\int\limits_{V}\int div(\lambda gradT)\operatorname{d}\!V, chúng ta sẽ thu được:

\int\limits_{t_1}^{t_2}\int\int\limits_{V}\int\limits[F+div(\lambda gradT)-c_\rho\frac{\partial T}{\partial t}]\operatorname{d}\!V\operatorname{d}\!t=0, từ đó suy ra phương trình dẫn nhiệt:

c_\rho\frac{\partial T}{\partial t}=\frac{\partial}{\partial x}(\lambda\frac{\partial T}{\partial x})+\frac{\partial}{\partial y}(\lambda\frac{\partial T}{\partial y})+\frac{\partial}{\partial z}(\lambda\frac{\partial T}{\partial z})+F(x,y,z,t)

Đó là phương trình vi phân với các đạo hàm riêng loại parabolic. Đôí với môi trường đồng nhất \lambda là một hằng số, phương trình sẽ đơn giản hơn và có dạng:

\frac{\partial T}{\partial t}=a^2\Delta T+f, trong đó \Delta là toán tử Laplace, còn a^2=\frac{\lambda}{c_\rho} là hệ số dẫn nhiệt, f=\frac{F}{c_\rho}.

Nếu cho các điều kiện giới hạn vá các điều kiện ban đầu, chúng ta có thể chọn được nghiệm duy nhất của phương trình truyền nhiệt.

## Tham khảo

* Nhiệt động lực học