Spatio-temporal air flow properties in a 3D personalised model of the human lung

Titolo: Proprietà del flusso d'aria spazio-temporali in un modello 3D personalizzato del polmone umano

1. Il Problema

La mucosa bronchiale protegge i polmoni intrappolando e rimuovendo le particelle inalate. In condizioni patologiche dove i meccanismi di clearance (clearance mucociliare e tosse) sono compromessi, si verifica un ristagno di muco che aumenta il rischio di infezioni. La fisioterapia toracica (CP) viene spesso prescritta per facilitare il trasporto del muco generando alti flussi d'aria, basandosi sull'ipotesi che tali flussi interagiscano con il muco promuovendone la mobilizzazione. Tuttavia, l'analisi di questa interazione e la dinamica risultante del muco lungo l'albero bronchiale rimangono largamente empiriche.

I modelli matematici e computazionali esistenti (spesso basati su modelli "Weibel-like") presentano un limite fondamentale: assumono che tutti i percorsi dalla trachea agli acini siano geometricamente e biofisicamente identici. Questa idealizzazione non riesce a catturare le eterogeneità spaziali reali, creando un "divario paradigmatico" nella modellazione delle distribuzioni tridimensionali dei processi fisici sottostanti, specialmente durante condizioni di ventilazione a riposo.

2. Metodologia

Gli autori propongono un modello polmonare multi-scala sintetico e personalizzato che integra geometria reale, meccanica dei tessuti e fluidodinamica.

• Geometria 3D e Ricostruzione:

  • Grandi vie aeree: La geometria delle grandi vie aeree e dell'involucro polmonare è ricostruita direttamente da scansioni TC (Tomografia Computerizzata) di pazienti (utilizzando segmentazioni dal challenge LUNA16).
  • Piccole vie aeree: Poiché le TC non risolvono le vie aeree più piccole, queste sono generate algoritmicamente. L'algoritmo utilizza una tessellazione di Voronoi per decomporre i polmoni in sottodomini e ricorsivamente suddivide le regioni fino a raggiungere il volume tipico di un acino. La generazione segue la legge di Hess-Murray per i raggi e mantiene un rapporto costante tra lunghezza e diametro, garantendo coerenza fisiologica.
  • Involucro toracico: Modellato come il guscio convesso dei polmoni dilatato del 10%, riempito da un mezzo elastico che simula la cavità pleurica.

• Meccanica dei Fluidi (Flusso d'aria):

  • Grandi vie aeree (3D): Il flusso è governato dalle equazioni di Navier-Stokes non lineari per fluidi incomprimibili, risolte con elementi finiti 3D e il metodo delle caratteristiche. Sono imposte condizioni di non scorrimento sulle pareti.
  • Piccole vie aeree (0D): Il flusso nelle vie distali è modellato in 1D/0D, includendo effetti di inerzia. La caduta di pressione dipende dalla resistenza idrodinamica, che varia in base al numero di Reynolds (flusso laminare o turbolento).
  • Accoppiamento: Il flusso è vincolato al tasso di variazione del volume delle regioni polmonari associate, calcolato tramite Lagrange multipliers.

• Meccanica dei Tessuti:

  • I tessuti polmonari sono modellati con elasticità a piccole deformazioni (non lineare), con un modulo di Young dipendente dalla deformazione locale, calibrato su dati di compliance statica.
  • L'equazione dello stress include sia la deformazione del tessuto che la pressione dell'aria locale.
  • Condizioni al contorno: Una pressione temporale applicata sulla parte inferiore dell'involucro simula il diaframma, mentre una pressione ridotta simula i muscoli costali. Il segnale di pressione è un'Heaviside smussata (1.3s inspirazione, 2.6s espirazione).

• Simulazione Numerica:

  • Utilizzo di una pipeline automatizzata che integra FreeFem++ (elementi finiti), MeshLab e GMSH.
  • Schema numerico semi-implicito con adattamento del passo temporale per stabilità (CFL < 1).

3. Contributi Chiave

1. Modello Ibrido Multi-scala: Prima integrazione di geometrie reali da TC per le vie aeree centrali con alberi generati algoritmicamente per le vie periferiche, mantenendo la coerenza spaziale tramite tessellazione di Voronoi.
2. Accoppiamento Fluido-Struttura Completo: Modellazione che lega esplicitamente la deformazione del tessuto polmonare, la pressione dell'aria locale e il flusso dinamico, considerando sia l'inerzia del fluido che la compliance delle vie aeree.
3. Analisi Spazio-Temporale ad Alta Risoluzione: Capacità di mappare distribuzioni non uniformi di velocità dell'aria e stress di taglio (shear stress) su un modello anatomicamente realistico, superando i limiti dei modelli simmetrici ideali.

4. Risultati

La simulazione è stata eseguita in condizioni di ventilazione a riposo (volume corrente: 469 ml) su un modello con 2.154 vie aeree ricostruite (diametri da 0.90 a 8.79 mm).

• Eterogeneità del Flusso: I pattern di flusso e di stress di taglio mostrano una complessità significativa. Regioni spazialmente vicine sperimentano pressioni e flussi diversi a causa della geometria intrinseca e della resistenza idrodinamica variabile dei percorsi.
• Asimmetria Inspiro/Espiro:
  • Gli stress di taglio medi sono più alti durante l'espirazione rispetto all'inspirazione.
  • Cause: 1) Asimmetria temporale del segnale di ventilazione; 2) Durante l'inspirazione, l'espansione delle vie aeree riduce la resistenza delle sotto-aree, mentre durante l'espirazione le vie rimangono quasi invariate.
• Distribuzione degli Stress di Taglio: I valori di shear rate coprono un ampio spettro di valori positivi e negativi, con una variabilità crescente nelle vie aeree più piccole.
• Flussi Retrogradi: Sono stati osservati flussi d'aria retrogradi (in direzione opposta al flusso globale), suggerendo una ridistribuzione locale dell'aria tra i compartimenti polmonari.
• Resistenza Idrodinamica: La resistenza totale calcolata è di $1.09 \times 10^5 \, \text{Pa}\cdot\text{m}^{-3}\cdot\text{s}$ e la compliance è di $0.475 \times 10^{-6} \, \text{m}^3\cdot\text{Pa}^{-1}$, valori coerenti con la fisiologia umana.

5. Significato e Implicazioni

Questo studio dimostra che anche in condizioni di riposo, i modelli polmonari multi-scala possono rivelare pattern spazio-temporali complessi che i modelli idealizzati non riescono a catturare.

• Fisioterapia Toracica (CP): L'applicazione di pressioni toraciche altamente eterogenee (come nella CP) potrebbe amplificare queste eterogeneità di flusso e stress di taglio, potenzialmente aiutando la mobilizzazione del muco. Tuttavia, la complessità dei pattern di flusso e la presenza di flussi retrogradi potrebbero interferire con gli obiettivi terapeutici, suggerendo la necessità di una personalizzazione delle terapie.
• Futuri Sviluppi: Il lavoro apre la strada a studi che includano grandi deformazioni (per simulare sforzi massimali o CP intensa) e alla validazione del modello tramite Tomografia a Impedenza Elettrica (EIT). Il modello offre uno strumento potente per comprendere le interazioni aria-muco e ottimizzare le strategie di trattamento per le malattie polmonari ostruttive.