Bilan d'énergie

Pour ce bilan d'énergie, quatre composantes sont prises en compte :

  • L'énergie cinétique : $ E_k (t) $
  • L'énergie potentielle due aux tensions interfaciales : $ E_p (t) $
  • L'énergie potentielle due à la force gravitationnelle : $ E_g (t) $
  • Le travail de dissipation visqueuse : $ L_f (t) $

La conservation de l'énergie implique que l'énergie totale du système ($ E_k + E_p + E_g $) doit décroître pendant la propagation à cause de la dissipation visqueuse ($L_f $). L'équation du bilan d'énergie global (OEB) est donc :

$$ \frac{d}{dt}  \left[ \underbrace{E_k(t)}_\text{énergie cinétique} + \overbrace{
\underbrace{E_p(t)}_\text{tension de surface}
\underbrace{E_g(t)}_\text{force gravitationnelle} }^\text{énergie potentielle} \right] 
+ \underbrace{\frac{d L_f (t)}{dt}}_\text{dissipation visqueuse} = 0 $$

Chaque terme va être calculé afin de déterminer l'équation finale.

  • L'énergie cinétique : $ E_k (t) $

Afin de calculer l'énergie cinétique, on considère que la goutte est composée de plusieurs segments, qui se propagent tous à la vitesse $ U_m (t)$ du centre de masse définie telle que :

$   U_{mx} (t) = \frac{dx_m}{dt} $ et $   U_{my} (t) = \frac{dy_m}{dt} $

Comme la masse volumique du fluide est constante, le centre de masse et le centre géométrique sont les mêmes, définies par :

$  x_m(t) =\frac{ \int_0^{R(t)} x y(x) dx }{ \int_0^{R(t) y(x) dx} } $ et $ y_m(t) =\frac{1}{2} \frac{ \int _0^{R(t)}  y^2(x) dx }{ \int_0^{R(t)}  y(x) dx } $, où $ y(x) $ représente le profil liquide-fluide, c'est-à-dire la la forme de la goutte.

Au final, l'énergie cinétique vaut :

$$ \boxed{ E_k (t) = \frac{m}{2}  \left[ (\frac{d x_m}{dt})^2 + (\frac{d y_m}{dt})^2 \right] } $$

A t= 0, $ \boxed{E_k(0) = 0} $ car la goutte est immobile initialement.

L'avantage d'utiliser le centre de masse est qu'il n'y a pas besoin de connaître le champ de vitesse à l'intérieur de la goutte, ou de supposer une distribution simple qui satisfait l'équation de continuité et les conditions limites; et c'est quasiment impossible de mesurer la vitesse à l'intérieur de la goutte. 

  • L'énergie potentielle due aux tensions interfaciales : $ E_p (t) $

L'énergie potentielle due aux tensions interfaciales  peut s'écrire : $ E_p (t) = E_{lf} (t) + E_{sf} (t) + E_{sl} (t)$ .

Avec $ E_i (t)  = {\gamma}_i*A_i (t) $, où $ E $ est l'énergie potentielle, $ \gamma $ la tension de surface et $ A $ la surface, et les indices $i = lf, sf, sl $ correspondant à l'interface Liquide-Fluide, Solide-Fluide et Solide-Liquide. En utilisant l'équation d'Young ( $ {\gamma}_{sf} = {\gamma}_{sl} + {\gamma}_{lf} cos {\theta}_e $ ) et en sachant que $ A_{sf} (t) + A_{sl} (t) = A_{total} $ (surface du solide) on a :

$$ \boxed{ E_p (t) = {\gamma}_{lf} \left[ A_{lf} (t) - A_sl(t) cos {\theta}_e \right] + {\gamma}_{sf} A_{total} } $$

Avec $ A_{lf} (t) = \int_{-R(t)}^{R(t)} L \frac{dy(x)}{dx} $ et $ A_{sl} = 2 R(t) L $ (où $L$ est l'épaisseur du domaine, suivant z).

A t=0, $ A_{LF} (0) = 2 \pi r_0 L $ et $ A_{sl} (0) =0 $, donc $ \boxed{ E_p (0) = 2 \pi r_0 L {\gamma}_{lf} + {\gamma}_{sf} A_{total}  }$

  • L'énergie potentielle due à la force gravitationnelle : $ E_g (t) $

On se base à nouveau sur l'approche du centre de masse, et donc l'énergie potentielle gravitationnelle de la goutte se propageant vaut donc :  $$  \boxed{ E_g (t) = m g y_m (t) } $$.

A t= 0 , $ \boxed{ E_g (0) = m g r_0} $.

  • Le travail de dissipation visqueuse : $ L_f (t) $

On utilise la formule fournie dans le papier (Y.Gu et D.Li) qui est issue d'une étude précédente, pour la force visqueuse : $ F_v (t) = \frac{3 {\mu}_l }{ {\theta}_d (t) } \frac{dR(t)}{dt} \ln{(\epsilon_\delta^{-1})} $, où $ \mu$ est la viscosité dynamique du fluide, $ \theta $ est l'angle de contact et $\epsilon $ une constante à modifier pour caler au mieux aux résultats numériques. 

La force visqueuse totale exercée par la ligne de contact sur la surface solide vaut $ 2 L F_v (t) $. Par conséquent, le travail total dû à la force visqueuse vaut : $ L_f(t) = \int_0^{R(t)} 2 L F_v (t) dR(t) $.

Finalement : $$ \boxed{ L_f (t) = A * \int_0^t \frac{ R(t) }{ {\theta}_d} \left[ \frac{dR(t)}{dt} \right]^2 dt } $$.

Avec $ A = 6 \pi {\mu}_l \ln{(\epsilon_\delta^{-1})} $, et $ \boxed{ L_f (0) = 0} $.

 

Finalement, l'équation d'équilibre d'énergie au temps t s'écrit $$ \Delta E_k (t) + \Delta E_p(t) + \Delta E_g (t) + \Delta L_f (t) = 0 $$.

En remplacent chaque terme par son expression calculée juste avant : $$ \Rightarrow \frac{m}{2} \left[ (\frac{d x_m}{dt})^2 + (\frac{d y_m}{dt})^2 \right] + \gamma_{lf} \left[A_{lf} - A_{sl} cos\theta_e - 2\pi r_0 L\right] +  mg (y_m -r_0) + A* \int_0^t \frac{L}{\theta(t)} \left[\frac{dR}{dt}\right]^2 dt =0  $$

Cette équation est une équation générale pour modéliser la propagation d'une goutte liquide cylindrique se propageant sur une surface solide. Elle permet de déterminer quantitativement comment le processus de propagation d'une petite goutte de liquide sur une surface solide dépend de la surface et des propriétés physiques du système liquide-fluide-solide, des conditions initiales et de la forme du profil à chaque instant.