Skip to main content

Fick's law and equimolar counter diffusion

Short answer

Fick's law: JA = –DAB dCA/dx (molar flux proportional to concentration gradient), analogous to Fourier's q'' = –k dT/dx. For steady equimolar counter diffusion through a tube of length L, with no net bulk flow, NA = –DAB dCA/dx constant. Integrating: NA = DAB (CA1–CA2)/L. For ideal gases, CA = pA/RT, thus NA = DAB (pA1–pA2)/(RT L). Concentration profile linear.

State Fick's law of diffusion and compare it with Fourier's law. Derive the expression for the molar flux in steady equimolar counter diffusion of two ideal gases A and B between two large vessels connected by a tube of length L.

Fick's law

Mass transfer by molecular diffusion is the movement of a species because of a concentration difference. For a binary mixture, the molar flux of A in the xx direction (relative to the mixture as a whole) is

JA=−DABdCAdxJ_A = -D_{AB}\frac{dC_A}{dx}

where DABD_{AB} is the diffusion coefficient (m2^2/s), CAC_A the molar concentration (kmol/m3^3). The flux is toward lower concentration.

QuantityHeat conductionMass diffusion
Lawq′′=−k dT/dxq'' = -k\,dT/dxJA=−DAB dCA/dxJ_A = -D_{AB}\,dC_A/dx
Driving forceTemperature gradientConcentration gradient
Propertykk (or α\alpha)DABD_{AB}
Units of propertyW/m K (m2^2/s)m2^2/s

Typical DABD_{AB}: gases about 10−510^{-5} m2^2/s, liquids about 10−910^{-9} m2^2/s, solids 10−1010^{-10} to 10−1410^{-14} m2^2/s.

Equimolar counter diffusion

Two vessels, at the same total pressure PP and temperature TT, are connected by a tube of length LL. Gas A diffuses from vessel 1 (partial pressure pA1p_{A1}) to vessel 2 (pA2p_{A2}), and gas B diffuses in the opposite direction at the same molar rate: NA=−NBN_A = -N_B.

 [ vessel 1 ]=======L=======[ vessel 2 ]
   pA1 high        A -->       pA2 low
   pB1 low         <-- B       pB2 high

Since there is no net bulk flow, the flux is only diffusion:

NA=−DABdCAdxN_A = -D_{AB}\frac{dC_A}{dx}

At steady state NAN_A is constant. Integrate from x=0x = 0 (CA1C_{A1}) to x=Lx = L (CA2C_{A2}):

NA=DABL (CA1−CA2)N_A = \frac{D_{AB}}{L}\,(C_{A1}-C_{A2})

For ideal gases CA=pA/RˉTC_A = p_A/\bar R T, so

NA=DABRˉT L (pA1−pA2)N_A = \frac{D_{AB}}{\bar R T\,L}\,(p_{A1}-p_{A2})

with Rˉ=8.314\bar R = 8.314 kJ/kmol K (or 8314 J/kmol K, with pp in Pa). The concentration profile of A is linear in xx, as for conduction through a plane wall.

✓ Answer checked by Claude Sonnet 5.5

All Heat Transfer chapter 7 questions →