A numerical study of turbulent Rayleigh-Bénard convection

Il flusso di Rayleigh-Bénard, costituito da un fluido racchiuso fra due pareti a diversa temperatura, rappresenta il paradigma della convezione termica. In natura e nelle applicazioni industriali, il moto convettivo avviene principalmente in regime turbolento, rivelando un fenomeno estremamente complesso. L'obiettivo principale di questo elaborato di tesi consiste nell'isolare e descrivere gli aspetti salienti di un flusso turbolento di Rayleigh-Bénard. L'analisi è applicata a dati ottenuti da tre simulazioni numeriche dirette effettuate allo stesso numero di Rayleigh (Ra=10^5) e a numeri di Prandtl differenti (Pr=0.7,2,7). Sulla base di alcune statistiche a singolo punto, vengono definite nel flusso tre regioni caratteritiche: il bulk al centro della cella, lo strato limite termico e quello viscoso in prossimità delle pareti. Grazie all'analisi dei campi istantanei e delle correlazioni spaziali a due punti, sono state poi individuate due strutture fondamentali della convezione turbolenta: le piume termiche e la circolazione a grande scala. L'equazione generalizzata di Kolmogorov, introdotta nell'ultima parte della trattazione, permette di approcciare il problema nella sua complessità, visualizzando come l'energia cinetica viene immessa, si distribuisce e viene dissipata sia nello spazio fisico, sia in quello delle scale turbolente. L'immagine che emerge dall'analisi complessiva è quella di un flusso del tutto simile a una macchina termica. L'energia cinetica viene prodotta nel bulk, considerato il motore del flusso, e da qui fluisce verso le pareti, dove viene infine dissipata.

Preconditioning strategies for the numerical solution of convection-diffusion partial differential equations

Il trattamento numerico dell'equazione di convezione-diffusione con le relative condizioni al bordo, comporta la risoluzione di sistemi lineari algebrici di grandi dimensioni in cui la matrice dei coefficienti è non simmetrica. Risolutori iterativi basati sul sottospazio di Krylov sono ampiamente utilizzati per questi sistemi lineari la cui risoluzione risulta particolarmente impegnativa nel caso di convezione dominante. In questa tesi vengono analizzate alcune strategie di precondizionamento, atte ad accelerare la convergenza di questi metodi iterativi. Vengono confrontati sperimentalmente precondizionatori molto noti come ILU e iterazioni di tipo inner-outer flessibile. Nel caso in cui i coefficienti del termine di convezione siano a variabili separabili, proponiamo una nuova strategia di precondizionamento basata sull'approssimazione, mediante equazione matriciale, dell'operatore differenziale di convezione-diffusione. L'azione di questo nuovo precondizionatore sfrutta in modo opportuno recenti risolutori efficienti per equazioni matriciali lineari. Vengono riportati numerosi esperimenti numerici per studiare la dipendenza della performance dei diversi risolutori dalla scelta del termine di convezione, e dai parametri di discretizzazione.

Numerical and experimental investigation of structural stiffness influence on ratcheting convection cell in granular soils under cyclic loading

Lateral cyclic loaded structures in granular soils can lead to an accumulation of irreversible strains by changing their mechanical response (densification) and forming a closed convective cell in the upper layer of the bedding. In the present thesis the convective cell dimension, formation and grain migration inside this closed volume have been studied and presented in relation to structural stiffness and different loads. This relation was experimentally investigated by applying a cyclic lateral force to a scaled flexible vertical element embedded in dry granular soil. The model was monitored with a camera in order to derive the displacement field by means of the PIV technique. Modelling large soil deformation turns out to be difficult, using mesh-based methods. Consequently, a mesh-free approach (DEM) was chosen in order to investigate the granular flow with the aim of extracting interesting micromechanical information. In both the numerical and experimental analyses the effect of different loading magnitudes and different dimensions of the vertical element were considered. The main results regarded the different development, shape and dimensions of the convection cell and the surface settlements. Moreover, the Discrete Element Method has proven to give satisfactory results in the modelling of large deformation phenomena such as the ratcheting convective cell.