La concentration d’un jus de fruit est réalisée par nanofiltration1 sur une membrane tubulaire de 6 mm de diamètre et de 1,2 m de longueur. La nanofiltration est réalisé en mode tangentielle (cross-flow filtration) c'est à dire avec une vitesse de circulation dans le tube (parrallèle à la membrane) et une vitesse de perméation (perpendiculaire à la membrane). Le jus de fruit a une concentration massique de 0,2 % et est retenu totalement par la membrane.
Données :
1 ces calculs peuvent aussi être utilisés pour de l'ultrafiltration ou de l'osmose inverse
1. En considérant l’accumulation du jus de fruit à la surface de la membrane, calculer la concentration en jus de fruit à la surface de la membrane, cm, pour un flux de perméation de 10-6 m/s et 2.10-6 m/s et des vitesses de circulation de 0,05 et 0,1 m/s. Commenter les résultats.
$$c=c_b e^{Jx/D}$$Lors de la filtration, il y a un couplage entre la convection due à la perméation, J, qui amène les solutés (ou les colloïdes) vers la membrane et la diffusion, D, dans la couche limite, $\delta$, qui engendre un "retro"-transport depuis la membrane vers la solution (ou la dispersion). On peut démontrer (voir cours) en réalisant un bilan différentiel sur un élément de volume dans la couche limite que le profil de concentration est exponentiel :
$$c_m=c_b e^{Pe}$$La concentration à la membrane, $c_m$, est alors :
$$Sh=1.86(Re~Sc~\frac{d_H}{L})^{0.33}$$où, $Pe=\frac{J\delta}{D}$, est le nombre de Péclet qui caractérise le rapport entre $\frac{le~flux~advectif}{le~flux~diffusif}$. Pour calculer le nombre de Péclet il est nécessaire d'estimer l'épaisseur de la couche limite. Il est possible de l'estimer à partir d'une corrélation entre nombre sans dimensions pour le transfert dans un tube.
pour Re<2100 $$Sh=0.023Re^{0.8}Sc^{0.33}$$ pour Re>2100
Le code suivant permet de :
- calculer l'épaisseur de couche limite
- caculer le nombre de Péclet
- calculer la concentration à la membrane
import numpy as np
import matplotlib.pyplot as plt
from scipy.optimize import newton
#DATA
#Solution
c0=0.2 #% massique
D=7.e-10 #m2/s
ro=1000 #kg m-3
mu=0.001 #Pa.s
#membrane
d_H=6e-3 #m
L=1.2 #m
# conditions hydrodynamiques
u=[0.05, 0.1] #m/s
J=[2.e-6,5.e-6] #m/s
#Calcul de l'épaisseur de couche limite, du nombre de Péclet et de la concentration à la membrane
deltal=np.zeros(len(u))
Re=np.zeros(len(u))
Sh=np.zeros(len(u))
Pe=np.zeros(len(u)*len(J))
cm=np.zeros(len(u)*len(J))
Sc=mu/(ro*D)
for i in range(len(u)):
Re[i]=ro*u[i]*d_H/mu
if Re[i]<2100:
Sh[i]=1.86*(Re[i]*Sc*d_H/L)**0.33
else:
Sh[i]=0.023*Re[i]**0.8*Sc**0.33
deltal[i]=d_H/Sh[i]
print ('L\'épaisseur de couche limite est de', round(deltal[i]*1e6, 1),'micromètres pour une vitesse tangentielle de ', u[i],'m/s' )
for j in range(len(J)):
Pe[i+len(u)*j]=J[j]*deltal[i]/D
cm[i+len(u)*j]=c0*np.exp(Pe[i+len(u)*j])
print()
print (' J \ u | ', round(u[0],3),' ', round(u[1],3))
print ('------------------------------------------------------------------------------')
print (' | Pe= ', round(Pe[0],3), ' ', round(Pe[1],3))
print (' ', "%.2e"%J[0],' | ')
print (' | cm= ', round(cm[0],3), ' ', round(cm[1],3))
print ('------------------------------------------------------------------------------')
print (' | Pe= ', round(Pe[2],3), ' ', round(Pe[3],3))
print (' ', "%.2e"%J[1],' | ')
print (' | cm= ', round(cm[2],3), ' ', round(cm[3],3))
#Tracé du profil de concentration près de la membrane
x0=np.linspace(-deltal[0],0,100)
x1=np.linspace(-deltal[1],0,100)
c_0=c0*np.exp(Pe[0]*(1+x0/deltal[0]))
c_1=c0*np.exp(Pe[1]*(1+x1/deltal[1]))
c_2=c0*np.exp(Pe[2]*(1+x0/deltal[0]))
c_3=c0*np.exp(Pe[3]*(1+x1/deltal[1]))
plt.plot(x0,c_0, label='u=0.05 m/s J=2 um/s')
plt.plot(x1,c_1, label='u=0.1 m/s J=2 um/s')
plt.plot(x0,c_2, label='u=0.05 m/s J=5 um/s')
plt.plot(x1,c_3, label='u=0.1 m/s J=5 um/s')
plt.legend(loc='upper center')
plt.title('Evolution de la concentration au voisinage de la membrane')
plt.xlabel('Distance à la surface de la membrane en mètre')
plt.ylabel('Concentration du jus de fruit en % massique')
plt.show()
L'épaisseur de couche limite est de 256.7 micromètres pour une vitesse tangentielle de 0.05 m/s L'épaisseur de couche limite est de 204.2 micromètres pour une vitesse tangentielle de 0.1 m/s J \ u | 0.05 0.1 ------------------------------------------------------------------------------ | Pe= 0.733 0.583 2.00e-06 | | cm= 0.416 0.358 ------------------------------------------------------------------------------ | Pe= 1.834 1.459 5.00e-06 | | cm= 1.251 0.86
On note que l'augmentation du flux de perméation engendre une accumulation à la membrane plus importante. Par contre, l'augmentation de la vitesse tangentielle permet de réduire l'accumulation à la membrane. On appelle d'ailleurs cette vitesse la vitesse de balayage tangentiel car elle permet de réduire l'accumulation à la membrane : l'épaisseur de couche limite est plus petite ce qui favorise le transfert de matière par diffusion vers la solution. Mais il faut noter que l'augmentation de la vitesse tangentielle nécessite de l'énergie : il faut donc trouver un compromis.
Le flux de solvant à travers la membrane est directement relié à la concentration à la surface du filtre par la relation suivante : $$J=\frac{\Delta P-\Delta\Pi}{\mu R_m}$$ où
2. Avec les relations précédentes, comment peut-on calculer la concentration à la membrane et le flux de perméation pour une différence de pression donnée ?
$$c_m=c_b e^{Pe}$$L'analyse du transfert de matière a permis de relier la concentration à la membrane, $c_m$ aux conditions opératoires, Pe :
$$J=\frac{\Delta P-\Pi(c_m)}{\mu R_m}$$L'analyse du transfert de quantité de mouvement (mécanique des fluides) permet de traduire l'effet de la concentration, $c_m$, à la membrane (par le biais de la pression osmotique, $\Pi(c_m)$) à travers la loi de Darcy modifiée :
$$\frac{D}{\delta}ln(\frac{c_m}{c_0})-\frac{\Delta P-10^{+5}\frac{133.75c_m}{100-c_m}}{\mu R_m}=0$$On dispose donc d'un système de deux équations à deux inconnues. Si la pression trans-membranaire est fixée, les inconnus sont le flux de perméation et la concentration à la membrane. Si on connait le flux à travers la membrane, les inconnus sont la pression transmembranaire et la concentration à la membrane. Il n' y a pas de solution analytique pour ces équations car le système est non linéaire. Il est par contre possible de combiner les deux équations afin d'avoir une équation à une inconnue ici par exemple en fonction de $c_m$ :
3. Réaliser ce calcul pour des pressions de 2, 5 et 10 bars et pour des vitesses de 0.05 m/s et 0.1 m/s.
Le code suivant permet de résoudre le système d'équations avec la séquence de calculs suivantes :