30.11.2012 Views

ETUDE EXPERIMENTALE ET NUMERIQUE DE LA DISPERSION ...

ETUDE EXPERIMENTALE ET NUMERIQUE DE LA DISPERSION ...

ETUDE EXPERIMENTALE ET NUMERIQUE DE LA DISPERSION ...

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

tel-00476994, version 1 - 27 Apr 2010<br />

24<br />

déflagration-détonation dans un écoulement diphasique dans le cas d’explosifs<br />

condensés, et en 1987 par Khasainov et Veyssière [10] pour modéliser la détonation<br />

de suspensions de particules solides dans des milieux gazeux réactifs. Ce modèle est<br />

utilisé dans le code EFAE [10] qui est à la base de l’étude numérique présentée dans<br />

ce mémoire.<br />

2.2. Les équations<br />

2.2.1. Equations constitutives du problème<br />

Les hypothèses du modèle existant dans le code EFAE (qui a servi de support à nos<br />

simulations numériques, cf. paragraphe 2.3) sont les suivantes : la taille moyenne des<br />

particules est suffisamment grande devant les distances de la dynamique moléculaire, mais<br />

petite devant les autres dimensions caractéristiques du problème. La fraction volumique des<br />

particules ainsi que les interactions inter-particules sont prises en compte. Les gradients de<br />

température à l’intérieur des particules sont supposés nuls et celles-ci sont considérées comme<br />

incompressibles. Le système d’équations permettant de modéliser l’écoulement considéré<br />

peut s’écrire de la manière suivante :<br />

a) Trois équations de bilan pour décrire la phase gazeuse :<br />

∂ρ<br />

Masse + ∇ ⋅(<br />

ρu)<br />

= J<br />

∂t<br />

u<br />

Quantité de mouvement u u<br />

v p f − + Φ −∇ = ⊗ ∇ +<br />

∂ρ<br />

((ρ ) ) P J<br />

∂t<br />

∂E<br />

Energie + ∇ ⋅ u ( E + PΦ))<br />

= J(<br />

e p + Qmélange<br />

) − f ⋅ v p − q<br />

∂t<br />

b) Trois équations de bilan pour la phase solide :<br />

∂σ<br />

Masse + ∇ ⋅(<br />

σv<br />

p ) = −J<br />

∂t<br />

∂σ<br />

v<br />

( 2.1)<br />

( 2.2)<br />

( ( 2.3)<br />

p<br />

Quantité de mouvement (( σ ) ) ( 1 )<br />

∂t<br />

p p p p<br />

( 2.4)<br />

+ ∇ v ⊗ v = −∇Π − Φ − Jv<br />

+ f ( 2.5)

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!