Analisi distribuzioni degli elementi delle matrici di un controllo h-infinito

Lo studio di tesi che segue analizza un problema di controllo ottimo che ho sviluppato con la collaborazione dell'Ing. Stefano Varisco e della Dott.ssa Francesca Mincigrucci, presso la Ferrari Spa di Maranello. Si è trattato quindi di analizzare i dati di un controllo H-infinito; per eseguire ciò ho utilizzato i programmi di simulazione numerica Matlab e Simulink. Nel primo capitolo è presente la teoria dei sistemi di equazioni differenziali in forma di stato e ho analizzato le loro proprietà. Nel secondo capitolo, invece, ho introdotto la teoria del controllo automatico e in particolare il controllo ottimo. Nel terzo capitolo ho analizzato nello specifico il controllo che ho utilizzato per affrontare il problema richiesto che è il controllo H-infinito. Infine, nel quarto e ultimo capitolo ho specificato il modello che ho utilizzato e ho riportato l'implementazione numerica dell'algoritmo di controllo, e l'analisi dei dati di tale controllo.

Mammografia digitale: confronto tra qualità immagine in relazione alle curve di acquisizione dosimetriche

Il presente lavoro è stato svolto presso il Servizio di Fisica Sanitaria dell’Azienda USL della Romagna, Presidio Ospedaliero di Ravenna e consiste nella validazione del dato dosimetrico visualizzato su due apparecchiature per mammografia digitale e nel confronto tra qualità immagine di diverse curve di acquisizione in funzione della dose e della post-elaborazione. Presupposto per l’acquisizione delle immagini è stata la validazione del dato dosimetrico visualizzato sui mammografi, tramite misura diretta della dose in ingresso con strumentazione idonea, secondo protocolli standard e linee guida europee. A seguire, sono state effettuate prove di acquisizione delle immagini radiografiche su due diversi fantocci, contenenti inserti a diverso contrasto e risoluzione, ottenute in corrispondenza di tre curve dosimetriche e in funzione dei due livelli di post-elaborazione delle immagini grezze. Una volta verificati i vari passaggi si è proceduto con l’analisi qualitativa e quantitativa sulle immagini prodotte: la prima valutazione è stata eseguita su monitor di refertazione mammografica, mentre la seconda è stata effettuata calcolando il contrasto in relazione alla dose ghiandolare media. In particolare è stato studiato l’andamento del contrasto cambiando le modalità del software Premium View e lo spessore interposto tra la sorgente di raggi X ed il fantoccio, in modo da simulare mammelle con grandezze differenti.

Quantum biology. Simulazioni di trasferimento di energia in una struttura dimerica

La quantum biology (QB) è un campo di ricerca emergente che cerca di affronta- re fenomeni quantistici non triviali all’interno dei contesti biologici dotandosi di dati sperimentali di esplorazioni teoriche e tecniche numeriche. I sistemi biologici sono per definizione sistemi aperti, caldi,umidi e rumorosi, e queste condizioni sono per loro imprenscindibili; si pensa sia un sistema soggetto ad una veloce decoerenza che sopprime ogni dinamica quantistica controllata. La QB, tramite i principi di noise assisted transport e di antenna fononica sostiene che la presenza di un adeguato livello di rumore ambientale aumenti l’efficienza di un network di trasporto,inoltre se all’interno dello spettro ambientale vi sono specifici modi vibrazionali persistenti si hanno effetti di risonanza che rigenerano la coerenza quantistica. L’interazione ambiente-sistema è di tipo non Markoviano,non perturbativo e di forte non equi- librio, ed il rumore non è trattato come tradizionale rumore bianco. La tecnica numerica che per prima ha predetto la rigenerazione della coerenza all’interno di questi network proteici è stato il TEBD, Time Evolving Block Decimation, uno schema numerico che permette di simulare sistemi 1-D a molti corpi, caratterizzati da interazioni di primi vicini e leggermente entangled. Tramite gli algoritmi numerici di Orthopol l’hamiltoniana spin-bosone viene proiettata su una catena discreta 1-D, tenendo conto degli effetti di interazione ambiente-sistema contenuti nello spettro(il quale determina la dinamica del sistema).Infine si esegue l’evoluzione dello stato.

Teoria dei giochi ed evoluzione La selezione naturale

La tesi fornisce una chiave di lettura del rapporto che si instaura tra la biologia e la matematica, nello specificio nell'analisi della Teoria dell'Evoluzione proposta da Charles Darwin e del processo di Selezione Naturale attraverso l'utilizzo di modelli matematici forniti dalla Teoria dei giochi e la Teoria dei giochi Evolutivi.

Applicazione di tecniche di analisi statistica e texture analysis all'imaging di risonanza magnetica di pazienti affetti da glioblastoma multiforme

Il seguente lavoro di tesi si è concentrato sull'analisi statistica dei dati prodotti dall'imaging di risonanza magnetica di pazienti affetti da tumori di alto grado, in particolare glioblastoma multiforme. Le tipologie di acquisizione d'immagine utilizzate sono state l'imaging pesato in T1 e il Diffusion-Weighted Imaging (DWI). Lo studio è stato suddiviso in due fasi: nella prima è stato considerato un campione di pazienti affetti da glioblastoma multiforme che, dopo il trattamento, avessero manifestato una ricaduta della malattia; per questi pazienti è stato quantificato in che modo la dose erogata durante la terapia si sia distribuita sul target del trattamento, in particolare nella porzione di tessuto in cui andrà a svilupparsi la recidiva. Nella seconda fase, è stato selezionato un campione più ristretto che disponesse, per entrambe le modalità di imaging, di un'acquisizione pre-terapia e di un numero sufficiente di esami di follow up; questo al fine di seguire retrospettivamente l'evoluzione della patologia e analizzare tramite metodi statistici provenienti anche dalla texture analysis, i dati estratti dalle regioni tumorali. Entrambe le operazioni sono state svolte tramite la realizzazione di software dedicati, scritti in linguaggio Matlab. Nel primo capitolo vengono fornite le informazioni di base relative ai tumori cerebrali, con un'attenzione particolare al glioblastoma multiforme e alle sue modalità di trattamento. Nel secondo capitolo viene fatta una panoramica della fisica dell'imaging di risonanza magnetica e delle tecniche di formazione delle immagini, con un'ampia sezione è dedicata all'approfondimento dell'imaging in diffusione. Nel terzo capitolo viene descritto il progetto, i campioni e gli strumenti statistici e di texture analysis utilizzati in questo studio. Il quarto capitolo è dedicato alla descrizione puntuale dei software realizzati durante questo lavoro e nel quinto vengono mostrati i risultati ottenuti dall'applicazione di questi ultimi ai campioni di pazienti esaminati.

La risposta cellulare ai diversi tipi di radiazione tramite espressione genica e radiobiologia sistemica

Questo progetto ha confrontato gli effetti indotti da diversi tipi di radiazioni, diversa intensità delle dosi, diverso rateo di dose su sistemi cellulari differenti. In particolare sono stati seguiti due studi differenti, finalizzati all’indagine degli effetti e dei meccanismi indotti da trattamenti radioterapici su cellule in coltura. Nel primo studio -EXCALIBUR- sono stati investigati i meccanismi di induzione e trasmissione del danno a basse dosi di radiazioni, in funzione della qualità della radiazione (raggi gamma e protoni) e della dose. Cellule di glioblastoma umano (T98G) sono state irraggiate con raggi gamma e protoni a due diverse dosi (0,25 Gy e 2 Gy); in questo studio è stata valutata e analizzata la variazione di espressione genica rilevata utilizzando la tecnologia dei microarray. Per mezzo dell’analisi statistica, con due software diversi, si è osservato come nelle cellule irraggiate si attivino i geni legati alla senescenza cellulare; questo risultato è significativo, visto che potrebbe rappresentare una prospettiva terapeutica interessante per molte neoplasie. Il secondo studio –Plasma Focus- ha lo scopo di ampliare le applicazioni nel settore medicale di una sorgente radiante che produce raggi X ad altissimo rateo di dose (plasma focus). In questo studio, l’attenzione è stata posta sulla preparazione dei campioni biologici per l’irraggiamento. Cellule di adenocarcinoma mammario (MCF7) sono state coltivate in laboratorio e posizionate all’interno di appositi portacampioni pronte per essere irraggiate con raggi X ad alto e a basso rateo di dose. Per mezzo della microscopia ottica e della citometria a flusso in fluorescenza, si è osservato come un rateo di dose elevato provochi danni cellulari superiori. L’analisi quantitativa ha mostrato che, nelle cellule trattate con il plasma focus, il 18% risulti danneggiato rispetto al 7% delle cellule di controllo; con i raggi X convenzionali risulta danneggiato l'8% di cellule, rispetto al 3% delle cellule di controllo.

Studio di un nuovo sistema angiotermografico tridimensionale per la diagnostica del tumore al seno.

L’angiotermografia dinamica (DATG) è una tecnica diagnostica per il tumore al seno non invasiva e che non richiede l'utilizzo di radiazioni ionizzanti. La strumentazione fornisce immagini della vascolarizzazione della mammella rilevate e visualizzate da tre placche a cristalli liquidi. Le immagini, sulla base della teoria della neoangiogenesi, permettono la diagnosi anche precoce del tumore al seno. Inoltre, la DATG può essere impiegata su donne di ogni età e grazie al suo basso costo può coadiuvare le altre tecniche per lo screening. Scopo di questo lavoro è la caratterizzazione della strumentazione e lo studio di un nuovo metodo di elaborazione che permetta di ottenere immagini tridimensionali della vascolarizzazione. Poiché le immagini fornite dalla DATG sono a colori, si è elaborato il calcolo della tonalità (hue) per poter avere un valore da mettere in relazione con la temperatura misurata con una termocamera. Sperimentalmente è stato osservato che le placche rispondono a temperature comprese rispettivamente tra 31,5°C e 34,0°C per quella con cornice bianca, con una risoluzione tra 0,50 e 1,50 mm, tra 32,5oC e 34,5°C per quella con cornice grigia, che è in grado di rivelare dettagli con dimensioni minime dell'ordine di 0,50-2,50 mm e tra 34,0°C e 36,5°C per quella con cornice nera, con una risoluzione compresa tra 0,25 mm e 1,50 mm. Le immagini fornite dalla DATG sono bidimensionali. Per ricavare l’informazione sulla profondità di provenienza del segnale, è stata acquisita una sequenza dinamica angiotermografica in vivo. Per ottenere la distribuzione spaziale della vascolarizzazione in tre dimensioni, sulla base dell’equazione del trasporto del calore di Fourier, è stata ipotizzata una relazione esponenziale tra contributo temporale nella formazione dell’immagine e profondità di provenienza del segnale. I primi test sperimentali mostrano con buona approssimazione che è possibile ottenere immagini tridimensionali angiotermografiche, anche se il metodo andrà ulteriormente affinato. In futuro si cercherà di realizzare dei fantocci dedicati che consentano di verificare ulteriormente la modalità con cui si formano le immagini e di sfruttare i dati raccolti per conoscere con maggior accuratezza la profondità da cui proviene il segnale. Si spera, infine, che questo studio possa contribuire a diffondere la conoscenza di tale metodica e il suo utilizzo sia per la diagnosi che per lo screening del tumore al seno.

Valutazione di un sistema di recupero di immagini dermatologiche

Questa tesi si inserisce in un progetto di ricerca fra il gruppo di Matematica della Visione del Prof. Ferri e CA-MI S.r.l. volto a progettare un sistema di recupero di immagini mediante il quale un dermatologo potrà acquisire l’immagine di una lesione e recuperare da un database classificato le immagini più somiglianti. Il concetto stesso di “somiglianza” è formalmente realizzato da una parte dell’omologia persistente (funzioni di taglia). Questa tesi utilizza tali metodi al fine di ottenere una combinazione ottimale dei diversi classificatori che si ottengono utilizzando la modularità intrinseca nella teoria. A questo scopo vengono impiegati due modelli e diversi metodi numerici.

Relazione tra scavi in sotterraneo e deformazioni di versante in terreni argillosi

Lo scavo di un tunnel in ammassi rocciosi comporta una ridistribuzione degli stress nell’intorno del foro e induce spesso deformazioni della superficie che si traducono frequentemente in movimenti verticali verso il basso. In letteratura esistono soluzioni analitiche in grado di quantificare gli stress e la subsidenza indotti da uno scavo in sotterraneo. Tali soluzioni funzionano però solo per semplici casi. Al fine di valutare i risultati delle soluzioni analitiche e poter considerare modelli di materiale che si avvicinano maggiormente al comportamento reale del terreno, gli stress e le deformazioni indotte vengono spesso analizzati attraverso modellazioni numeriche. Tali modellazioni vengono svolte di norma tramite l’utilizzo di codici numerici bidimensionali alle differenze finite o agli elementi finiti. L’utilità di un approccio modellistico è particolarmente evidente nel caso dello studio degli spostamenti indotti in un versante dallo scavo di un tunnel. In questi casi le analisi numeriche mostrano che se si considera un versante ideale costituito da materiale elastico, lo scavo provoca spostamenti che tendono a convergere verso il tunnel. Se si considera invece un versante costituito da materiale elasto-plastico, lo scavo, in alcuni casi, può causare uno spostamento verso valle e verso il basso dell’intero versante. Si è voluto quindi indagare il complesso rapporto tra scavo e deformazioni utilizzando un approccio modellistico. A tale scopo è stato condotto uno studio sul tipo di deformazioni che avvengono in un versante teorico, interessato da scavo, a partire da differenti condizioni di stabilità iniziale e considerando tunnel di diverse dimensioni e contropressioni. I versanti considerati nelle modellazioni riproducono le condizioni tipiche dei versanti appenninici argillosi, caratterizzati da basso angolo di attrito interno, geometrie relativamente poco inclinate e coesione fortemente variabile. I risultati ottenuti sono stati successivamente testati analizzando un caso reale di interazione scavo-deformazione relativo a un versante appenninico costituito da Argille a Palombini.

Analisi di facies e stratigrafia dei depositi tardo-quaternari della Pianura Padana centrale

La ricostruzione dettagliata dell’architettura stratigrafica e delle facies delle successioni tardo-quaternarie sepolte al di sotto delle pianure alluvionali moderne sta assumendo un’importanza crescente come strumento per le ricostruzioni paleoambientali e paleoclimatiche e per la correlazione dei sistemi deposizionali continentali e costieri, finalizzata alla caratterizzazione geometrica di corpi sedimentari che rappresentano potenziali serbatoi. Un impulso significativo verso la caratterizzazione dell’architettura stratigrafica dei depositi tardo-quaternari del Bacino Padano è stato fornito, negli ultimi quindici anni, da un’intensa campagna di perforazioni operata dal Servizio Geologico, Sismico e dei Suoli di Regione Emilia-Romagna, nell’ambito della realizzazione della Carta geologica d’Italia in scala 1:50.000. Tuttavia, per ottenere un numero sufficiente di informazioni stratigrafiche è necessario un dataset molto esteso e le perforazioni (sondaggi a carotaggio continuo) presentano costi estremamente elevati. I limiti operativi delle descrizioni dei sondaggi a carotaggio continuo possono essere superati utilizzando metodi indiretti che si basano sullo sfruttamento delle proprietà geotecniche dei depositi tardo-quaternari. In questo lavoro di tesi l’analisi delle caratteristiche geotecniche delle prove CPTU è stata sfruttata per caratterizzare le associazioni di facies della Pianura Padana centrale. Una dettagliata analisi dei log stratigrafici e, in particolare, delle curve di resistenza delle prove CPTU inserite nelle sezioni, la costruzione di sezioni stratigrafiche e un meticoloso lavoro di correlazione basato su 19 datazioni radiometriche, ha rivelato una serie di paleosuoli (inceptisuoli) impilati nei 50-60 m più superficiali della successione tardo-quaternaria della Pianura Padana centrale; questi paleosuoli rappresentano l’unico marker stratigrafico all’interno di una successione costituita prevalentemente da depositi argillosi uniformi. L’analisi delle caratteristiche geotecniche di alcuni paleosuoli e dei loro rapporti con la channel belt del Fiume Po risulta fondamentale per comprendere l’evoluzione recente del Bacino Padano e le caratteristiche geometriche dei corpi sedimentari finalizzata allo sfruttamento di potenziali risorse fluide sotterranee.

Caratterizzazione geochimica dei sedimenti nel sito Caretti di Ferrara

Lo studio riguarda la caratterizzazione geochimica dei sedimenti del sito Caretti, situato a est del comune di Ferrara, interessato da un forte caso di contaminazione dovuto alla presenza di due discariche municipali non controllate risalenti agli anni '60-'70. L'obiettivo del lavoro è quello di valutare la geochimica dei sedimenti e delle acque prelevate nel sito al fine di individuare eventuali relazioni tra il chimismo delle due matrici. Inoltre, i dati dell' analisi dei sedimenti sono stati confrontati con quelli di letteratura per delinearne aspetti di provenienza. I dati geochimici dei sedimenti sono stati ottenuti mediante fluorescenza a raggi X (XRF) e sono stati elaborati graficamente e statisticamente, mediante software di analisi statistica (GCDkit). La composizione geochimica delle acque è stata ottenuta mediante tecniche di cromatografia ionica e spettroscopia di assorbimento atomico. La geochimica dei sedimenti ha evidenziato che la composizione chimica mostra associazioni di elementi legate a una componente tessiturale del sedimento. In particolare si rilevano concentrazioni elevate di Fe2O3, Al2O3, TiO2, Cr, Ni, e V in sedimenti a granulometria fine, a differenza delle concentrazioni di SiO2 e Na2O elevate in sedimenti sabbiosi. La distribuzione statistica degli elementi nei sedimenti e nelle acque ha permesso di effettuare una valutazione delle caratteristiche geochimiche dei principali corpi idrogeologici. La caratterizzazione geochimica ha fornito, inoltre, indicazioni sull'origine dei sedimenti. Cromo e Nichel, tra gli elementi in traccia, sono particolarmente efficaci nella valutazione della provenienza permettendo, all'interno del sito Caretti, una distinzione tra popolazioni di provenienza appenninica (bassi valori Cr-Ni) e padana (alti valori Cr-Ni).

Petrografia e geochimica delle vulcaniti del Monte dei Porri di Salina (Isole Eolie)

E’ stata indagata dal punto di vista stratigrafico, geochimico e petrografico la successione di vulcaniti appartenenti alle formazioni Rocce di Barcone (Tufi Grigi dei Porri o Grey Porri Tuffs – GPT), Punta Sallustro, Punta di Megna e Serro dello Sciarato (Lucchi et al., 2013a, b), riferibili all’Epoca Eruttiva 5 (~70 - 57 ka) dell’evoluzione vulcanologica dell’isola Salina (isole Eolie), durante la quale si realizza l’edificazione della porzione medio-basale dello stratocono del Monte dei Porri. Attraverso un approccio basato sull’integrazione e correlazione di metodologie stratigrafiche, petrografiche e geochimiche si è giunti ad una ricostruzione dei meccanismi eruttivi, deposizionali e petrogenetici che sottendono la successione di eventi vulcanici responsabili dell’edificazione dello stratocono suddetto. Utilizzando l’analisi di litofacies si è giunti, anche attraverso le correlazioni stratigrafiche estese sull’isola di Lipari, alla definizione di cinque unità eruttive (EU1-EU5), rappresentanti il riferimento stratigrafico per la campionatura necessaria per caratterizzarle sia petrograficamente che geochimicamente. Il significato vulcanologico delle EU1-5 ha permesso di evidenziare che il Monte dei Porri è costituito, in larga parte, da un’articolata alternanza di depositi da caduta e da correnti piroclastiche (formati sia da pomici che da scorie), alla quale si intercalano colate laviche; all’interno dell’unità i prodotti analizzati evidenziano composizione dei magmi da trachi-dacitica (EU1), a basaltica e andesitico basaltica (EU2), a francamente andesitico basaltica (EU3); quest’ultima composizione caratterizza la EU4 (che evolve anche verso termini andesitici), e la EU5. In conclusione, l’integrazione di tutti i dati (stratigrafici, petrografici e geochimici) permette di affermare che l’attività vulcanica responsabile della messa in posto delle EU1-5 cominci con una fase di apertura del condotto eruttivo (fase esplosiva pliniana) accompagnata dall’emissione dei magmi più evoluti (EU1), residenti al top di una camera magmatica zonata. Ad essa presumibilmente segue il coinvolgimento dei livelli più profondi del reservoir magmatico dove risiedono i magmi più mafici (EU2 e EU3). La ripresa dell’attività vulcanica (dopo una stasi durante il quale i magmi mafici evolvono verso composizioni mediamente evolute) vede infatti l’emissione di prodotti andesitici (EU4) seguita da magmi meno evoluti (EU5) che porta ad un progressivo svuotamento del sistema.

Laplace operator on finite graphs and a network diffusion model for the progression of the Alzheimer disease

Nella tesi viene descritto il Network Diffusion Model, ovvero il modello di A. Ray, A. Kuceyeski, M. Weiner inerente i meccanismi di progressione della demenza senile. In tale modello si approssima l'encefalo sano con una rete cerebrale (ovvero un grafo pesato), si identifica un generale fattore di malattia e se ne analizza la propagazione che avviene secondo meccanismi analoghi a quelli di un'infezione da prioni. La progressione del fattore di malattia e le conseguenze macroscopiche di tale processo(tra cui principalmente l'atrofia corticale) vengono, poi, descritte mediante approccio matematico. I risultati teoretici vengono confrontati con quanto osservato sperimentalmente in pazienti affetti da demenza senile. Nella tesi, inoltre, si fornisce una panoramica sui recenti studi inerenti i processi neurodegenerativi e si costruisce il contesto matematico di riferimento del modello preso in esame. Si presenta una panoramica sui grafi finiti, si introduce l'operatore di Laplace sui grafi e si forniscono stime dall'alto e dal basso per gli autovalori. Al fine di costruire una cornice matematica completa si analizza la relazione tra caso discreto e continuo: viene descritto l'operatore di Laplace-Beltrami sulle varietà riemanniane compatte e vengono fornite stime dall'alto per gli autovalori dell'operatore di Laplace-Beltrami associato a tali varietà a partire dalle stime dall'alto per gli autovalori del laplaciano sui grafi finiti.

Algoritmi per la stima del rumore nella ricostruzione di immagini

In questa tesi viene analizzato un metodo di regolarizzazione nel quale il problema di regolarizzazione è formulato come un problema di ottimizzazione vincolata. In questa formulazione è fondamentale il termine noto del vincolo, parametro per il quale dipende la bontà della ricostruzione, che corrisponde ad una stima della norma del rumore sul dato. Il contributo della tesi è l'analisi di un algoritmo proposto per stimare il rumore in applicazioni di image deblurring utilizzando la decomposizione in valori singolari dell'operatore che descrive il blur. Sono stati fatti numerosi test per valutare sia i parametri dell'algoritmo di stima che l'efficacia della stima ottenuta quando è inserita nel termine noto del vincolo del problema di ottimizzazione per la ricostruzione dell'immagine.

La valutazione del rischio di incidente rilevante per i depositi di fitosanitari

I prodotti fitosanitari sono sostanze attive e preparati, contenenti una o più sostanze attive, destinati a proteggere i vegetali, a favorire e regolare i loro processi vitali, a conservarli, a eliminare le parti non necessarie, a frenare o evitare un loro indesiderato accrescimento, ad eliminare le piante indesiderate. Il prodotto fitosanitario è composto normalmente da tre tipologie di sostanze: sostanza attiva, coadiuvante e coformulante. Le sostanze attive rappresentano la parte del prodotto fitosanitario che serve a combattere l’avversità che si vuole controllare, cioè la sostanza tossica che, a seconda della pericolosità e della concentrazione presente nella formulazione, concorre a determinare la classe di tossicità e quindi di pericolosità del prodotto. Gli stabilimenti che detengono fitosanitari sono suddivisi in depositi e siti produttivi. Secondo la classificazione dell’ISPRA, i depositi di fitosanitari costituiscono una delle 15 categorie di attività a rischio di incidente rilevante, mentre, i siti produttivi rientrano nella categoria degli stabilimenti chimici e/o petrolchimici. Il presente lavoro di tesi ha riguardato la valutazione del rischio di incidente rilevante per i depositi di fitosanitari, con particolare attenzione al rischio associato alla produzione di composti tossici a seguito di scenari incidentali di incendio. È stata in primo luogo effettuata una ricerca approfondita finalizzata a comprendere che cosa sono i prodotti fitosanitari, a quale scopo servono e quali sono i rischi ad essi associati. Successivamente, è stato effettuato un approfondimento delle norme a cui i siti che producono e detengono fitosanitari sono soggetti. È stato quindi effettuato un censimento per individuare quali e quanti sono i depositi e i siti di produzione di prodotti fitosanitari suscettibili di causare incidenti rilevanti in Italia. Si è poi proceduto all’analisi storica sulle banche dati MHIDAS, E-MARS ed ARIA, al fine di individuare gli eventi incidentali che hanno interessato depositi e siti produttivi di fitosanitari. Successivamente, si è analizzato il Regolamento CPR-15, quale metodo raccomandato dal ministero olandese VROM (Ministerie van Volkshuisvesting, Ruimtelijke Ordening en Milieu) per la valutazione del rischio associato a depositi di fitosanitari. Come applicazione della metodologia per la valutazione del rischio descritta dal Regolamento CPR 15, è stata effettuata la modellazione di uno scenario di incendio in un deposito reale di fitosanitari attraverso il Warehouse Model, implementato nel software PHAST 6.54. La simulazione ha portato a definire l’estensione delle aree di danno associate alla dispersione dei fumi tossici prodotti dall’incendio.