Modélisation physique et résolution numérique du transport et du dépôt des particules d’aérosols médicaments dans les voies respiratoires extrathoraciques humaines, Physical modelling and numerical resolution of aerosol medicine particles transport and deposition in human extrathoracic airways

Julien Sandeau

Description

Sous la direction de Daniel Isabey
Thèse soutenue le 18 mars 2009: Paris Est
D'abord développées pour traiter les maladies respiratoires (asthme, BPCO...), les thérapies inhalées sont depuis peu l'objet d'un intérêt particulier pour délivrer des molécules dans le sang via la barrière alvéolo-capillaire. Quelle que soit la cible de ces thérapies, leur efficacité dépend de la quantité de médicament qui arrive sur la zone à traiter dans l'arbre respiratoire. La région extrathoracique (appelée région ET et comprenant la bouche, l'oropharynx, le larynx et le début de la trachée) se comporte alors comme un filtre où se dépose une certaine masse de médicament. Afin d'optimiser une telle administration, il est nécessaire de connaître au mieux l'influence des facteurs caractéristiques de l'inhalation (entre autres, le débit d'inhalation, les diamètres des particules inhalées, les propriétés physiques du gaz porteur). Pour cela, un modèle réaliste de région ET a été construit à partir des données issues de l'imagerie médicale. Après la définition d'un maillage et le choix d'un modèle de turbulence et de paramètres adéquats, les équations de Navier-Stokes ont été résolues numériquement pour (i) décrire l'écoulement dans ce modèle et (ii) pour déterminer les trajectoires des particules inhalées. En parallèle, deux études expérimentales ont été menées sur l'exacte reproduction physique du modèle : l'une pour estimer les pertes de pression provoquées par des écoulements d'air et d'oxhel (mélange hélium (78 %) - oxygène (22 %)), l'autre pour mesurer les vitesses par IRM d'un écoulement d'hélium-3 hyperpolarisé. Les résultats obtenus en terme d'écoulement et de quantité de particules déposées dans la région ET sont en accord avec les rares données expérimentales existantes. Ils montrent l'influence prépondérante de la géométrie tant sur l'écoulement (marqué par la présence, vérifiée expérimentalement, du jet laryngeal) que sur le dépôt. Les zones préférentielles de dépôt se situent là où des variations géométriques brusques surviennent (rétrécissement de section et courbures près de l'épiglotte et en amont de la glotte), et là où le caractère turbulent de l'écoulement est prononcé (en fin de larynx notamment). Ces constatations ont permis de caractériser avec précision le dépôt à l'aide de deux nombres adimensionnés : le nombre de Stokes, caractéristique de l'impaction inertielle et le nombre de Reynolds, caractéristique du régime de l'écoulement. L'étude a permis une meilleure appréhension des influences des divers paramètres à considérer pour modéliser une inhalation d'aérosol médicament et, de surcroît, une prévision des quantités de médicament déposée et transmise à l'arbre respiratoire, soit un outil précieux pour l'optimisation des thérapies inhalées. La méthodologie, ainsidéveloppée, est applicable aux morphologies issues d'une étude clinique en cours et permettra d'affiner la caractérisation
-Particules d'aérosols médicaments
-Voies respiratoires extrathoraciques
Inhaled therapies were originally developed to treat respiratory diseases (e.g., asthma and COPD). They are now the subject of research to deliver medicine to the blood through the alveolar-capillary barrier. The efficacy of such therapies lies in the quantity of medicine that reaches the targeted area in the respiratoy tree. The extrathoracic region (called ET and including the mouth, the oropharynx, the larynx and the beginning of the trachea) has the role of a filter limiting the amount of medicine that is deposited in the targeted area. In order to optimize the administration of aerosol medicine, it is necessary to know the influence of the various parameters of an inhalation (among others, the inhalation flow rate, the diameter of inhaled particles, the physical properties of the cariier gas). Therefore, a realistic model of the ET region has been constructed from data obtained from medical imaging. The Navier-Stokes equations were numerically solved (i) to describe the flow in the model and (ii) to determine the trajectories of inhaled particles. The numerical analysis required the development of an appropriate mesh and the choice of turbulence model and parameters. In parallel, two experimental studies were performed on an exact physical reproduction of the model : one to estimate the pressure losses in air and oxhel (a helium (78 %) - oxygen (22 %) mixture) flows and another to measure the velocity field using a MRI of a hyperpolarized helium-3 flow. The flow and deposition results that were obtained in the model of the ET region are in agreement with the few existing experimental data. They show the dominant effect of the geometry not only on the flow (marked by the experimentally verified presence of the laryngeal jet) but also on the particle deposition. The major areas of deposition are located where sudden geometrical changes are observed (like the constrictions and the bends near the epiglottis and upstream the glottis) and where the turbulent character of the flow is pronounced (particulary at the end of the larynx). These observations allowed for the accurate characterization of the deposition with non-dimensional numbers : the Stokes number, characteristic number of the inertial impaction, and the Reynolds number, characteristic number of the flow regime. The study allows for the improved comprehension of the influences of the various parameters that must be considered to model the inhalation of aerosol medicine. Furthermore, predictions in terms of the deposited and transmitted mass of medicine were possible. The methodology explained herein will be applicable to the morphologies that are to be obtained from an on-going clinical study
Source: http://www.theses.fr/2009PEST0003/document

Informations

• Publié par: Thesee
• Nombre de lectures: 60
• Langue: Français
• Poids de l'ouvrage: 6 Mo

Extrait

UNIVERSITÉ PARIS-EST
Sciences et Ingénierie
Thèse de Doctorat en Sciences
Pour l’obtention du grade de
DOCTEUR DE L’UNIVERSITÉ PARIS-EST
Spécialité : Sciences de l’Ingénieur
Soutenue par :
Julien Sandeau
Modélisation Physique et Résolution Numérique
du Transport et du Dépôt des Particules
d’Aérosols Médicaments
dans les Voies Respiratoires Extrathoraciques Humaines
Soutenance le 18 mars 2009 devant le jury composé de :
C. Baroud CC, HDR, LadHyX, Ecole Polytechnique, Palaiseau, France Rapporteur
M. Riethmuller PU, von Karman Institute, Rhode-St-Genèse, Belgique Rapporteur
O. Boiron MCU, HDR, Ecole Centrale Marseille, Marseille, France Examinateur
G. Caillibotte Ing. R&D., Air Liquide, Les Loges-en-Josas, France Examinateur
B. Housset PU-PH, Hôpital Henri Mondor, Créteil, France Examinateur
D. Isabey DR, INSERM U955, Créteil, France Directeur de thèse
G. Apiou-Sbirlea Resp. Groupe R&D, Air Liquide, Les Loges-en-Josas, France Membre invité
X. Maître CR, CNRS UMR8081, Orsay, France Membre invitéKremlin-Bic?tre,
1
cristolienne
Progressiv
e,
emen
(et
t,
onne
ce
qui
m?moire
conseils),
prend
enn,
sa
appris,
forme
se
d?nitiv
!)
e.
bres
Alors
et
que
son
j'y
cidre,
app
Orsa
ose
P
quelques
toine,
derni?res
?
retouc
2
hes,
olescu,
je
ph
v
c
oudrais
biom?canique
remercier
at
:
v
Daniel,
?
bien
des
?videmmen
des
t,
c
directeur
que
de
Cr?teil,
th?se
?
toujours
et
disp
b
onible,
C?cile,
p
Une
our
qui,
m'a
ma
v
c,
oir
o
fait
or
d?couvrir
de
le
sique
monde
t
de
des
la
l'?quip
connaissance
et
et
le
comprendre
R?douane,
certains
jamais
de
en
ses
vier,
enjeux,
p
Georges,
moteur
resp
?riences
onsable
ensateur
de
indues
p
soup
?le
de
scien
et
tique,
crois?s,
de
,
pro
Loges-en-Josas,
jets,
V
de
aris,
th?ses,
et
de
m'on
stages...
p
(et
Carolin,
m?me
erta,
d'in
dam,
terfaces
hance
graphiques
p
!)
et
p
les
our
t
ses
Paris,
conseils
2008
et
fond
son
o
encadremen
et
t
et
constan
Br
ts,
t
Gabriela,
la
p
y-
our
et
son
respiren
r?le
!),
cen
hacun
tral
mem
d'initiatrice
de
de
e
ce
pulmonaire
pro
cellulaire
jet,
(et
les
triumvir
rapp
:
orteurs
Bruno
de
Gabiel,
ce
a
tra
ares
v
sages
ail,
Xa
Mic
actif
hel
l'inni,
Riethm
our
uller
r?le
et
lors
Charles
exp
Baroud,
IRM
rencon
disp
tr?s
?
resp
heures
ectiv
de
emen
de
t
es
lors
hinoises,
d'une
r?sonances...),
soutenance
celles
de
ceux
th?se
j'ai
et
?
lors
y
d'un
?
cours
aux
de
au
biom?canique
?
du
alence,
p
P
oumon,
?
et
alaiseau
les
ailleurs,
examinateurs,
qui
Bruno
t
Housset
un
et
eu,
Charles
eaucoup.
Boiron,
Pierre-An
p
Rob
our
Mathieu,
a
A
v
b
oir
c
accept?
!
d'?tre
derni?re
les
ens?e
jur?s
celles
de
ceux
ce
tous
tra
jours,
v
balladen
ail,
dans
c
t?te.
hacune
le
des
janvier
p
(ave
ersonnes
en
du
sonor
group
L
e
ckwo
m?dico-gazier
d
logeois
Nic
(et
Swe
particuli?remen
Ge
t
gia
celles
own
qui
Julien
fon2extraordinaires.
3
corrige
Si
oin
timide
ten
que
p.
l'on
orne
soit,
oudrait
il
n
faut
et
bien
pas
que
on
l'on
l'exp
in
ph
terp
corrections
ole
de
;
d'?noncer
l'exp
L
?rience
Champs-Flammarion,
ne
ne
nous
ces
donne
eux-m?mes.
qu'un
se
certain
?
nom
on
bre
et
de
qui
p
de
oin
se
ts
vraimen
isol?s,
?rience
il
serait
faut
lois
les
P
r?unir
scienc
par
oth?se
un
de
trait
elle
con
passera
tin
par
u
p
;
ts
c'est
Ainsi
l?
ne
une
b
v
pas
?ritable
g?n?raliser
g?n?ralisation.
?rience,
Mais
la
on
;
fait
le
plus,
ysicien
la
v
courb
s'abstenir
e
ces
que
et
l'on
con
tracera
ter
passera
t
en
l'exp
tre
toute
les
ue
p
forc?
oin
des
ts
bien
observ
Henri
?s
oincar?,
et
a
pr?s
e
de
l'hyp
ces
,
p
?dition
oin
1968,
ts
159.
;4inhalation
5
talemen
Mo
d'h?lium-3
d?lisation
constatations
Ph
pro
ysique
les
et
section
R?solution
et
Num?rique
La
du

T
tit?
ransp
de
ort
?t
et
le
du
de
D?p
une
?t
os?e
des
cours
P
d'air
articules
mesurer
d'A
us
?rosols
ET
M?dicamen
Ils
ts
(marqu?
dans
d?p
les
ariations
V
amon
oies
prononc?
Respiratoires
v
Extrathoraciques
bre
Humaines
r?gime
R?sum?
param?tres
D'ab
pr?vision
ord
outil
d?v
applicable
elopp
p
?es
par
p
(m?lange
our
%)),
traiter
IRM
les
olaris?.
maladies
t
respiratoires
os?es
(asthme,
accord
BPCO...),
?rimen
les
l'inuence
th?rapies
t
inhal?es
v
son
laryngeal)
t
pr?f?ren
depuis
l?
p
t
eu
de
l'ob
glotte),
jet
de
d'un
larynx
in
ermis
t?r?t
d?p
particulier
adimensionn?s
p
caract?ristique
our
de
d?livrer
L'?tude
des
inuences
mol?cules
our
dans
et,
le
de
sang
l'arbre
via
l'optimisation
la
d?v
barri?re
d'une
alv
d'aner
?olo-capillaire.
de
Quelle
o
que
?coulemen
soit
d'o
la
(78
cible
xyg?ne
de
p
ces
vitesses
th?rapies,
?coulemen
leur
yp
ecacit?
r?sultats
d?p
terme
end
de
de
particules
la
la
quan
t
tit?
v
de
donn?es
m?-
existan
dicamen
tren
t
ond?ran
qui
g?om?trie
arriv
l'?coulemen
e
la
sur
exp
la
du
zone
sur
?
Les
traiter
de
dans
situen
l'arbre
des
respiratoire.
brusques
La
t
r?gion
courbures
extrathoracique
et
(app
de
el?e
l?

turbulen
r?gion
t
ET
n

t).
et
t
comprenan
caract?riser
t
pr?cision
la
?
b
nom
ouc
le
he,
Stok
l'oropharynx,
l'impaction
le
nom
larynx
caract?ristique
et
l'?coulemen
le
p
d?but
appr?hension
de
div
la
consid?rer
trac
d?liser
h?e)
m?dicamen
se
surcro?t,
comp
quan
orte
t
alors
transmise
comme
soit
un
p

th?rapies
ltre
dologie,

?e,
o?
morphologies
se
clinique
d?p
p
ose
caract?risation.
une
ertes
certaine
pression
masse
v
de
qu?es
m?dicamen
des
t.
ts
An
et
d'optimiser
xhel
une
h?lium
telle
%)
administration,
o
il
(22
est
l'autre
n?cessaire
our
de
les
conna?tre
par
au
d'un
mieux
t
l'inuence
h
des
erp
fac-
Les
teurs
obten
caract?ristiques
en
de
d'?coulemen
l'inhalation
et
(en
quan
tre
de
autres,
d?p
le
dans
d?bit
r?gion
d'inhal