Analisi del moto di particelle in sospensione mediante simulazioni di fluidodinamica computazionale

La mobilitazione di polveri radioattive nel caso di un incidente di perdita di vuoto (LOVA) all’interno di ITER (International Thermonuclear Experimental Reactor), è uno dei problemi di sicurezza che si sono posti durante la costruzione di tale reattore. Le polveri vengono generate dalla continua erosione da parte del plasma del materiale di contenimento. Ciò porta ad un accumulo delle stesse all’interno della camera di vuoto. Nel caso di un incidente LOVA il rilascio di tali polveri in atmosfera rappresenta un rischio per la salute di lavoratori e della popolazione circostante. Per raccogliere dati su tale tipo di incidente è stata costruita presso il laboratorio dell’università di Tor Vergata una piccola facility, STARDUST, in cui sono stati eseguiti vari esperimenti con delle condizioni iniziali simili a quelle presenti all’interno di ITER. Uno di questi esperimenti in particolare simula la rottura della camera di vuoto mediante l’ingresso di aria in STARDUST, inizialmente posto a 100 Pa, con un rateo di pressurizzazione di 300 Pa s−1. All’interno del serbatoio sono presenti delle polveri che, in differente percentuale, vengono portate in sospensione dal flusso di fluido entrante. In particolare le polveri sono composte da tungsteno (W), acciaio inossidabile (SS – 316 ) e carbonio ( C ). Scopo del presente lavoro è quello di riprodurre il campo di velocità che si genera all’interno di STARDUST nel caso dell’esperimento appena descritto e valutare il moto delle particelle portate in sospensione e la loro successiva deposizione. Ciò viene fatto mediante l’utilizzo di una geometria bidimensionale creata con Salome. Su tale geometria sono costruite differenti mesh strutturate in base al tipo di simulazione che si vuole implementare. Quest’ultima viene poi esportata nel solutore utilizzato, Code_Saturne. Le simulazioni eseguite possono essere suddivise in tre categorie principali. Nella prima (Mesh A) si è cercato di riprodurre i risultati numerici presentati dagli autori della parte sperimentale, che hanno utilizzato il software commerciale Fluent. Nella seconda si è riprodotto il campo di velocità interno a STARUDST sulla base dei dati sperimentali in possesso. Infine nell’ultima parte delle simulazioni si è riprodotto il moto delle particelle sospese all’interno del fluido in STARDUST, valutandone la deposizione. Il moto del fluido pressurizzante è stato considerato come compressibile. Il regime di moto è stato considerato turbolento viste le elevate velocità che provocavano elevati valori del numero di Reynolds per l’aria. I modelli con cui si è trattata la turbolenza sono stati di tre tipi. Il campo di velocità ottenuto è stato leggermente differente nei tre casi e si è proceduto ad un confronto per valutare le cause che hanno portato a tali differenze. Il moto delle particelle è stato trattato mediante l’utilizzo del modello di tracciamento lagrangiano delle particelle implementato in Code_Saturne. Differenti simulazioni sono state eseguite per tenere in considerazione i vari modelli di turbolenza adottati. Si è dunque proceduto ad analizzare similitudini e differenze dei risultati ottenuti.

Dinamica delle galassie a spirale e delle galassie ellittiche

In questa tesi si vuole fornire una descrizione generale delle dinamiche delle galassie ellittiche e a spirale. Nel primo capitolo si danno informazioni generali sulle grandezze che caratterizzano le galassie e come esse vengono classificate. Nel secondo capitolo si espone il concetto di sistema collisionale, si fa notare come le galassie risultino essere sistemi non collisionali e come questo porti delle semplificazioni nella trattazione di questi oggetti e ne spieghi alcune caratteristiche. Si prosegue andando a considerare le equazioni che descrivono il moto (equazione non collisionale di Boltzmann, equazioni di Jeans, teorema del viriale in forma tensoriale) e le informazioni che si possono ricavare. Nel terzo capitolo ci si concentra sulle galassie ellittiche e sulle principali leggi che le descrivono e dalle quali è possibile ottenere stime riguardo distanza e dimensioni. Il quarto e ultimo capitolo è incentrato sulle galassie a spirale e in particolare sulla dinamica del disco, e come si è giunti all'ipotesi dell'esistenza della materia oscura, e sulla dinamica dei bracci.

Dinamica dei fluidi

Questo elaborato presenta i principali contenuti della fluidodinamica. Data la vastità dell'argomento, verranno mostrati i tratti essenziali, senza addentrarsi nei dettagli. Il capitolo iniziale introduce il concetto di fluido e i principali metodi di studio. Segue il corpo centrale dell'elaborato, in cui vengono presentate le equazioni fondamentali della dinamica dei fluidi. Successivamente due brevi capitoli si concentrano su una particolare forma dell'equazione del moto e sul tema della vorticità. L'ultimo capitolo tratta brevemente del collasso gravitazionale, una fra le innumerevoli applicazioni della teoria dei fluidi al campo dell'astrofisica. Vettori e tensori verranno rappresentati in coordinate cartesiane e, dove necessario, saranno richiamate in nota le identità utilizzate nei calcoli.

Risolubilità delle equazioni polinomiali

Scopo di questo elaborato è studiare la risolubilità per radicali di un polinomio a coefficienti in un campo di caratteristica zero attraverso lo studio del gruppo di Galois del suo campo di spezzamento. Dopo aver analizzato alcuni risultati su gruppi risolubili e gruppi semplici, vengono studiate le estensioni radicali e risolubili. Viene inoltre dimostrato su un campo K di caratteristica zero il Teorema di Galois, che caratterizza i polinomi risolubili per radicali f a coefficienti in K attraverso la risolubilità del gruppo di Galois G(L/K), dove L è il campo di spezzamento di f. La tesi contiene anche un'esposizione sintetica del metodo introdotto da Lagrange per la risoluzione di equazioni polinomiali di cui si conosca il gruppo di Galois.

Teorie della gravitazione f(R)

Tra tutti i fenomeni naturali osservabili, ne era presente uno particolarmente interessante e con il quale si aveva diretto contatto quotidianamente: la gravità. Dopo le innumerevoli osservazioni astronomiche effettuate da Galileo, fu Newton nel diciassettesimo secolo a capire che il moto dei pianeti era governato dalle medesime leggi che descrivono la caduta dei gravi sulla Terra e fu quindi lui che ci fornì una prima teoria della gravità con la quale si spiegarono le orbite dei pianeti con ottima precisione. Grazie al contributo di Einstein, la teoria si rinnovò e si arricchì, ma rimase pur sempre lontana dall' essere completa, tant' è che ancora oggi sono presenti molte domande a cui non siamo in grado di rispondere. In questo articolo ci occuperemo di tali quesiti, provando a formulare una teoria che sia in accordo con le attuali evidenze sperimentali. Nella prima parte, tratteremo le ragioni che hanno spinto i ricercatori ad introdurre le nuove teorie della gravità f(R); in particolare vedremo la peculiarità delle curve di rotazione delle galassie e perché ci sia il bisogno di tirare in ballo la materia oscura. Discuteremo anche alcuni problemi derivanti dall' evoluzione cosmica e altre incongruenze riguardanti la stabilità delle stelle di neutroni. In seguito mostreremo come ricavare l' equazione di Einstein partendo dai principi variazionali di Hamilton, e estenderemo tale ragionamento con lo scopo di ottenere un' equazione corrispondente ad una gravità modificata. Infine, verranno introdotte le teorie della gravità f(R), per mezzo delle quali cercheremo di discutere alcune possibili spiegazioni alle problematiche mosse nella parte introduttiva.

Il metodo di Newton e le sue varianti per sistemi di equazioni non lineari

In questa tesi sono stati descritti i principali metodi numerici per la risoluzione di sistemi non lineari. Tali metodi sono stati analizzati sia dal punto di vista teorico (analisi di convergenza locale) che pratico (algoritmo e implementazione).

Resistenze al moto in ambito ferroviario: Stato dell'arte, confronto tra rotabili passeggeri e merci

L'obiettivo della seguente trattazione è quello di illustrare i principali metodi di determinazione delle resistenze al moto per materiali rotabili. Verranno quindi proposti e discussi una serie di diversi approcci al problema corredati da esempi e commenti finalizzati a chiarire, oltre al funzionamento degli strumenti, anche l'ambito di validità di questi, evidenziando eventuali differenze nei risultati e interpretandole razionalmente. Nella prima parte verranno introdotti ed analizzati i singoli contributi che, durante l'avanzamento di un generico veicolo ferroviario, generano resistenze al moto: partendo dalle resistenze di tipo meccanico dovute alle caratteristiche del rotabile e alle sue proprietà meccaniche fino a quelle di tipo aerodinamico che, come si vedrà in seguito, dipendono fortemente non solo dalla tipologia di veicolo ma anche dalle condizioni ambientali di esercizio, ma soprattutto dalla velocità di percorrenza. Si passa quindi ad illustrare i principali metodi di determinazione empirici adottati dalle compagnie ferroviarie e dai produttori di materiale rotabile in vari paesi. Nell'ultima parte di questa trattazione si propone una serie di considerazioni finalizzate ad analizzare eventuali soluzioni pratiche per migliorare l'efficienza e la performance dei materiali rotabili sulla base delle informazioni ottenute dallo studio delle varie esperienze analizzate nei capitoli.

Modelli dinamici per fenomeni collettivi in fisica dei sistemi complessi

In questo lavoro si parla di un particolare comportamento emergente nei sistemi complessi: il flocking. Dopo aver dato una panoramica generale di come sia stato affrontato finora lo studio della formazione degli stormi, vengono portati come esempio un modello matematico e uno studio empirico particolare, il quale fonda le basi del flocking su un’interazione topologica. Il modello matematico, basato su un’interazione metrica, viene dapprima presentato, cercando di darne una parziale spiegazione tramite le proprietà delle matrici laplaciane, per poi essere testato attraverso delle simulazioni numeriche. Lo studio empirico, invece, viene presentato nei dettagli fornendo risultati e ipotesi atte a spiegarli. Infine prendendo spunto da questi due lavori, nell’ultima parte vengono posti a confronto due modelli, uno metrico e uno topologico, tramite simulazioni al calcolatore. L’esito di queste simulazioni conferma le ipotesi avanzate per spiegare i risultati empirici.

Equazioni stocastiche backward nel mercato delle emissioni

In questa tesi viene esposto il modello EU ETS (European Union Emission Trading Scheme) per la riduzione delle emissoni di gas serra, il quale viene formalizzato matematicamente da un sistema di FBSDE (Forward Backward Stochastic Differential Equation). Da questo sistema si ricava un'equazione differenziale non lineare con condizione al tempo finale non continua che viene studiata attraverso la teoria delle soluzioni viscosità. Inoltre il modello viene implementato numericamente per ottenere alcune simulazioni dei processi coinvolti.

Studio del trapping adiabatico per mappe stocastiche

In questa tesi vengono trattati argomenti relativi alla dinamica dei fasci di particelle: in particolare si è preso in considerazione il moto betatronico di una particella carica all'interno di un acceleratore circolare. Vengono quindi discussi alcuni aspetti della dinamica trasversa introducendo il formalismo Hamiltoniano e discutendo il modello presentato da Hénon per il caso bidimensionale. Viene poi introdotta la teoria adiabatica al fine di studiare gli effetti intrappolamento di un ensemble di particelle. Infine vengono presentate alcune simulazioni che permettono di poter osservare come il rumore rappresenti un fattore di rilevante importanza nello studio di tali fenomeni.

Analisi di un modello matematico unificato per lo studio del processo angiogenico

In questa tesi si studia l'angiogenesi tumorale, dapprima descrivendo i fenomeni biologici alla base della dinamica cellulare, e successivamente, dopo aver introdotto gli strumenti matematici necessari, sviluppandone un modello seguendo la letteratura esistente basato sulle equazioni differenziali stocastiche e su quelle di Fokker-Planck. Ne vengono infine realizzate simulazioni numeriche.

Equazioni stocastiche lineari e disuguaglianza di Harnack

In questo elaborato si presentano alcuni risultati relativi alle equazioni differenziali stocastiche (SDE) lineari. La soluzione di un'equazione differenziale stocastica lineare è un processo stocastico con distribuzione multinormale in generale degenere. Al contrario, nel caso in cui la matrice di covarianza è definita positiva, la soluzione ha densità gaussiana Γ. La Γ è inoltre la soluzione fondamentale dell'operatore di Kolmogorov associato alla SDE. Nel primo capitolo vengono presentate alcune condizioni necessarie e sufficienti che assicurano che la matrice di covarianza sia definita positiva nel caso, più semplice, in cui i coefficienti della SDE sono costanti, e nel caso in cui questi sono dipendenti dal tempo. A questo scopo gioca un ruolo fondamentale la teoria del controllo. In particolare la condizione di Kalman fornisce un criterio operativo per controllare se la matrice di covarianza è definita positiva. Nel secondo capitolo viene presentata una dimostrazione diretta della disuguaglianza di Harnack utilizzando una stima del gradiente dovuta a Li e Yau. Le disuguaglianze di Harnack sono strumenti fondamentali nella teoria delle equazioni differenziali a derivate parziali. Nel terzo capitolo viene proposto un esempio di applicazione della disuguaglianza di Harnack in finanza. In particolare si osserva che la disuguaglianza di Harnack fornisce un limite superiore a priori del valore futuro di un portafoglio autofinanziante in funzione del capitale iniziale.

Analisi dei fenomeni franosi lenti nell'abitato di Gaggio Montano mediante interferometria radar

I movimenti lenti delle colate in terra sono una caratteristica geomorfologica comune nell’Appennino settentrionale e sono uno dei principali agenti di modellazione del paesaggio. Spesso case e piccoli centri abitati sorgono in zone affette da questo tipo di movimento franoso e di conseguenza subiscono danni causati da piccoli spostamenti. In questo lavoro di Tesi vengono presentati i risultati ottenuti dall’interferometria radar ad apertura sintetica (InSAR) mediante elaborazione tramite StaMPS (Stanford Method of Persistent Scatterers), utilizzando la tecnica avanzata Small Baseline Subset (Berardino et al., 2002). Questo metodo informatico è applicato alle acquisizioni rilevate dai satelliti Envisat e COSMO-SkyMed in orbita ascendente e discendente, ottenendo una copertura di dati che va dal 2004 al 2015, oltre ad un rilevamento geologico-geomorfologico in dettaglio eseguito nell’area di studio. Questa tecnica di telerilevamento è estremamente efficace per il monitoraggio dei fenomeni di deformazione millimetrica che persistono sulla superficie terrestre, basata sull'impiego di serie temporali d’immagini radar satellitari (Ferretti et al., 2000). Lo studio è stato realizzato nel paese di Gaggio Montano nell’Appennino bolognese. In questa zona sono stati identificati diversi corpi di frana che si muovono con deformazioni costanti durante il tempo di investigazione e grazie ai risultati ottenuti dai satelliti è possibile confrontare tale risultato. Gli spostamenti misurati con il metodo InSAR sono dello stesso ordine di grandezza dei movimenti registrati dai sondaggi inclinometrici. Le probabili cause dell’instabilità di versante a Gaggio Montano sono di natura antropica, in quanto alti tassi di deformazione sono presenti nelle zone dove sorgono case di recente costruzione e complessi industriali. Un’altra plausibile spiegazione potrebbe essere data dalla ricarica costante d’acqua, proveniente dagli strati dei Flysch verso l’interno del complesso caotico sottostante, tale dinamica causa un aumento della pressione dell’acqua nelle argille e di conseguenza genera condizioni d’instabilità sul versante. Inoltre, i depositi franosi rilevati nell’area di studio non mostrano nessun tipo di variazione dovuta ad influenze idrologiche. Per questo motivo le serie temporali analizzare tendo ad essere abbastanza lineari e costanti nel tempo, non essendo influenzate da cicli stagionali.

Invarianza di gauge in elettromagnetismo e nelle onde gravitazionali

In questa tesi si studia l'uso dell'invarianza di gauge e la sua applicazione alla fisica nello studio dell'Elettromagnetismo e della Gravità. In particolare si fa uso dell'invarianza di gauge per studiare le soluzioni in forma di onde piane delle equazioni di Maxwell e dell'equazione di campo di Einstein. Nella presente tesi si mostra dunque come sia possibile applicare uno stesso procedimento matematico a due fenomeni fisici distinti e trarne conclusioni simili.

Physical characterization and numerical modeling of a WEST-like tokamak divertor plasma

Il presente elaborato è incentrato sulla modellizzazione del plasma di bordo nei dispositivi per la produzione di energia da fusione nucleare noti come tokamak. La tecnologia che nel corso di tutta la seconda metà del XX secolo fino ad oggi è stata sviluppata a questo fine deve necessariamente scontrarsi con alcuni limiti. Nei tokamak il confinamento del plasma è di tipo magnetico e vincola le particelle a muoversi di moto elicoidale all'interno del vessel, tuttavia il confinamento non risulta perfetto e parte dell'energia si scarica sulle pareti della camera, rischiando pertanto di fondere i materiali. Alcune strategie possono essere messe in atto per limitare questo problema, per esempio agendo sulla geometria del tokamak, oppure sulla fisica, inducendo nel plasma una data concentrazione di impurezze che ionizzino irraggiando parte dell'energia di plasma. Proprio tale meccanismo di perdita è stato simulato in un modello monodimensionale di plasma monofluido di bordo. I risultati del codice numerico relativo al modello dimostrano che per concentrazioni di impurezze crescenti è possibile diminuire in modo significativo flusso di calore e temperatura al divertore. Per di più risulta possibile controllare la posizione del fronte di irraggiamento per mezzo di parametri di controllo del plasma quali la pressione. Si osserva inoltre l'insorgere del cosiddetto fenomeno di biforcazione alle basse temperature di divertore, fenomeno in cui il plasma si comporta in modo instabile a causa di fenomeni fisici tipici delle basse energie ("detachment") e a seguito del quale può improvvisamente spegnersi (disruzione). Infine lo stesso modello è stato migliorato inserendo l'ipotesi di plasma bifluido. Anche per gli ioni viene osservato il fenomeno di biforcazione. I risultati numerici evidenziano le dinamiche dello scambio energetico fra le specie gettando le basi di una progettazione efficiente della chimica del plasma finalizzata al raffreddamento del divertore.