http://www.openstarts.units.it:80/dspace/handle/10077/2548 2013-06-16T17:09:14Z http://www.openstarts.units.it:80/dspace/handle/10077/8662 Title: Numerical investigation of evaporation and condensation of thin films in conjugated heat transfer systems Authors: Sosnowski, Pawel Abstract: Evaporation and condensation of thin liquid films on solid surfaces are common elements of industrial processes. In many cases they have a significant impact on the physics of the studied case. At the same time, experimental studies can prove to be troublesome, mostly because of the amount of possible setups, complex geometries of interest, numerous materials being used and cost. For that reason it is reasonable to study this phenomena using numerical methods. Having the advantage in speed and cost of performance, computational studies become a valuable tool. For evaporation and condensation process, one has to deal with buoyancy driven fluid flows, conjugated heat transfer between gaseous and solid phases, film thickness modeling, vapor phase behavior, and phase transition of the thin fluid film into vapor phase. The strong conjunction and mutual interaction of mentioned effects is the main focus of presented work. The gas phase behavior is being calculated using incompressible Navier-Stokes equations under Boussinesq approximation. The solutions of the partial differential equations are obtained with numerical methods using Eulerian finite volume discretization (Kundu and Cohen [2002]). Time advancement is being treated with second order implicit discretization. For cases with high Reynolds number, large eddy simulation (LES) techniques are used. Due to the complexity of the geometries of interest a dynamic computation of the Smagorinsky constant is preferred, applying the lagrangian dynamic model proposed by Meneveau et al. [1996]. The liquid film present on the surface of the solids is modeled following Petronio[2010]. Since the film is thin, it is assumed that it can be represented only by its thickness. This also leads to assumption that the heat transfer through the film is instantaneous. The vapor is represented by concentration of this phase in the volume of gas. The concentration is transported by convection and diffusion. The phenomena of evaporation and condensation of the thin films are driven by the presence of concentration gradients next to the surfaces. Phase transition of vapor to fluid, or other way around, acts on the energy balance, id. est latent heat is released into the gas when condensation occurs or the solid is cooled during evaporation. The heat transport is modeled in both solid and fluid domains. The case is split into separate regions with different material properties. These regions are solved one by one in a serial way using numerical techniques consistent with domain decomposition methods described by Quarteroni and Valli [1999]. The energy transport among the regions is performed by applying a heat coupling boundary conditions. The main focus of this work is to provide a reliable model for simulation system with complex physics involving fluid motion, heat transport in multi region domains (fluid-solid), vapor transport, thin film evolution and evaporation and condensation effects on energy balance. Proposed model is validated on simple geometries and later applied to problem of evaporation in vertical channel flow. The reference to the channel case is work of Laaroussi et at. [2009]. Presented study aims in providing comprehensive insights into physical effects that appear when the solid wall is being directly modeled and when latent heat transformations are taken into account. The final test is performed on a vertical channel with forced turbulent flow, directly modeled solid walls and evaporation or condensation happening on the boundary. Having the model working within such complex frame allows for its future usage in elaborate industrial applications. Description: 2011/2012 Type: Tesi di dottorato; Doctoral Thesis 2013-04-23T00:00:00Z http://www.openstarts.units.it:80/dspace/handle/10077/8661 Title: Caratterizzazione della permeabilità in mezzi porosi sintetici e naturali. Authors: Siena, Martina Abstract: La presente tesi ha come principale obiettivo lo studio della variabilità di proprietà idrologiche in mezzi porosi, con particolare attenzione alla permeabilità. A tal fine, ci si avvale di un approccio che combina l'analisi di proprietà statistiche e di scaling applicata a dataset di permeabilità, con lo studio di risultati numerici di simulazioni di flusso alla microscala in mezzi porosi. Con la prima analisi è possibile caratterizzare variazioni di permeabilità alla scala di misura (tipicamente dell'ordine del centimetro), mentre la seconda analisi dà una descrizione dell'eterogeneità di permeabilità ad una scala inferiore (nell'ordine del millimetro), ottenuta risolvendo processi fisici alla scala dei pori e derivando le quantità integrali di interesse. L'analisi statistica e di scaling, effettuata sia su distribuzioni di permeabilità sintetiche, sia su dataset raccolti su campioni reali, avvalora la validità dei modelli truncated fractional Brownian motion (tfBm) e truncated fractional Gaussian noise (tfGn), o di processi random sub-Gaussiani ad essi subordinati, per l'interpretazione della variabilità di proprietà idrologiche. Soluzioni numeriche di campi di flusso (i.e. velocità e pressione) alla scala dei pori sono ottenute sia per campioni sintetici, sia per campioni reali, la cui geometria è ricostruita mediante micro-tomografia a raggi X. Diverse metodologie di applicazione delle condizioni al contorno in corrispondenza dell'interfaccia liquido-solido forniscono risultati qualitativamente simili sia in termini di quantità microscopiche, sia in termini di quantità medie.; The work is aimed at providing some insights on the variability of hydrological properties in porous media, focusing in particular on permeability. We consider an approach which combines scaling and statistical analyses of air-permeability datasets with pore-scale numerical simulations of flow through porous media. The former investigation allows to characterize permeability heterogeneity at the centimeter observation scale; the latter provides a description of heterogeneity on a millimeter scale by resolving physical processes occurring at the microscopic scale and deriving up-scaled quantities. Scaling and statistical analyses performed on synthetic permeability distributions as well as on datasets collected on real media support the identification of truncated fractional Brownian motion (tfBm) or truncated fractional Gaussian noise (tfGn) and of sub-Gaussian random processes subordinated to tfBm (or tfGn) as viable models for the interpretation of hydrological properties variability. Pore-scale numerical solutions of flow (i.e., in terms of velocity and pressure distributions) are performed on both randomly generated samples and real porous media reconstructed via X-ray Micro-Tomography. Different approaches for the enforcement of boundary conditions at the fluid-solid interface provide qualitatively similar results in terms of both microscopic and averaged quantities. Description: 2011/2012 Type: Tesi di dottorato; Doctoral Thesis 2013-04-23T00:00:00Z http://www.openstarts.units.it:80/dspace/handle/10077/8656 Title: Dynamical characteristics of the Adriatic Sea using Lagrangian methods Authors: Hariri, Saeed Abstract: Abstract Lagrangian drifters are especially suitable to study transport processes since they move with the currents following closely the motion of water parcels, at least at the large and mesoscale,considering that transport by ocean currents plays a crucial role in many scientific and applied marine issues. In the coastal zone and semi-enclosed seas, transport by surface currents are even more important due to the variety of flotsam and jetsam, like debris, pollutants, oil spills, persons lost at sea, etc., encountered in the vicinity of populated coasts and highly navigated areas. The Adriatic Sea is an interesting and challenging place for such kind of problems because of strong boundary currents and interior gyres. The objective of this thesis is to simulate and study passive tracer dispersion by statistical properties of the currents in the Adriatic Sea from velocity data produced by numerical model outputs and real drifter data set. Each chapter of thesis consists of the content of manuscript which is ready to be submitted to peer reviewed journals. In chapter 2, Statistics of transit and residence times in the surface Adriatic Sea, a semi- enclosed basin of the Mediterranean, are estimated from drifter data and Lagrangian numerical simulations. The results obtained from the drifters are generally underestimated given their short operating lifetimes (half life of ∼ 40 days) compared to the transit and residence times. This bias can be removed by considering a large amount of numerical particles whose trajectories are integrated over a long time (750 days) with a statistical advection-diffusion model of the Adriatic surface circulation. Numerical particles indicate that the maximum transit time to exit the basin is about 216-260 days for objects released near the northern tip of the Adriatic, and that a particle entering on the eastern Otranto Channel will typically exit on the other side of the Channel after 170-185 days. A value of 150-168 days is estimated for the residence time in the Adriatic basin. In chapter 3, we examine the dispersion characteristics of surface drifter pairs in the Adriatic Sea. Relative dispersion of all surface drifters deployed in the Adriatic Sea between 25 August 1990 and 1 January 2007 has been calculated for different initial separation distance, using both chance and original pairs. Relative dispersion ( D 2 ) is explained in two different regimes. First non local regime which changes as a exponential function of time, the second one is local regime and it can be divided in three separated parts, Richardson ( D 2 ∼ t 3 ), ballistic ( D 2 ∼ t 2 ) and diffusive ( D 2 ∼ t ) regimes. The distribution of finite-scale Lyapunov exponent (FLSE) fields in the Adriatic Sea and Lagrangian Structure Function (LSF) have also been calculated from drifter trajectories, which both describe intrinsic physical properties at a given scale, also values of mean separation angle, displacement kurtosis and absolute dispersion have been found which can help us to know more about the behavior of pair trajectories. The aim of this work in chapter 4 is to quantify the near-surface transport properties and residence times in the Adriatic semi-enclosed basin by using synthetic drifters. We analyzed the simulated trajectories computed from the daily averaged velocity fields obtained by the MIT general circulation model implemented in the Adriatic Sea and integrated for the period from October 2006 till the end of 2008. Each numerical particle trajectory was obtained by integrating and interpolating velocity field between grid points using a fourth-order Runge–Kutta scheme and bilinear interpolation. In particular the surface circulation properties in two contrasting years (2007 with mild winter and cold autumn, 2008 with normal winter and hot summer) are here compared. A comparison between the transport statistics for numerical particles crossing three selected sections located along the Italian coast (Conero and Gargano Promontories and Strait of Otranto) and the similar statistics driven by an existing climatology of the Adriatic surface velocity field (obtained by drifters measurements) has been carried out in order to corroborate the model results. Results indicate that the numerical particles are slower in this simulation when comparing them with the particles simulated by the flow field obtained by real drifters. This is because of the less energetic flow field generated by the MIT general circulation model during the selected years. Lagrangian statistics for the entire Adriatic basin after removing the mean Eulerian circulation for numerical particles have also been calculated and it can be found that the values of mean angular momentum, diffusivity and Lagrangian velocity covariance are less than the real drifter observations, but maximum Lagrangian integral time scale is the same. Because of the weather condition observed in 2007 and 2008, and the different kinetic energy of the mean flow (the yearly averages of MKE are 3.1408e-009 and 3.7907e-009 ( km 2 s −2 ) in 2007 and 2008 respectively ) and the mean eddy kinetic energy (2.5638e-009 (2.8961e-009) ( km 2 s −2 ) in 2007(2008)) during these years, the comparison between transport properties and Lagrangian statistics has been done considering these two different periods. The obtained results showed that the effects of wind driven recirculation in north of the Po River (which would be as a sea response to the Bora wind field) and Po River discharge on surface circulation induce the value of residence time to be similar during two years (182 (185) days in 2007 (2008)). Keywords: Adriatic Sea, Lagrangian Statistics, Transit and Residence Times, Relative Dispersion, Finite-Scale Lyapunov Exponent.; Riassunto Le correnti e i processi di trasporto oceanici a grande scala e a mesoscala svolgono un ruolo cruciale in una varietà di campi, siano essi puramente d’interesse scientifico o applicativo. Tra le diverse strumentazioni oceanografiche a disposizione per lo studio della ircolazione oceanica, le boe di deriva di tipo lagrangiano (o “drifters” lagrangiani) sono particolarmente adatte allo studio dei processi di trasporto in oceano data la loro capacità di seguire il moto delle particelle d’acqua e, in definitiva, delle correnti. Nelle aree costiere e nei bacini semi-chiusi, lo studio dei processi di trasporto associato alle correnti è particolarmente importante a causa dell’impatto di inquinanti, di sostanze rilasciate da relitti, di oil-spill sugli ecosistemi costieri, sulle attività economiche e sulle zone costiere densamente abitate. Il mare Adriatico, un bacino semi-chiuso del mar Mediterraneo, è, da questo punto di vista, un’area sensibile e presenta una serie di complessi fenomeni di circolazione quali intense correnti costiere e strutture vorticose a differente scala spaziale. Lo scopo di questa tesi è quello di simulare e quindi studiare i fenomeni di dispersione di traccianti passivi attraverso le proprietà statistiche delle correnti nel mare Adriatico utilizzando correnti superficiali prodotte da modelli di circolazione e da dati sperimentali derivati da drifters reali. La tesi èsuddivisa in capitoli strutturati come articoli scientifici indipendenti e “auto-consistenti”, pronti ad essere sottomessi a riviste scientifiche di tipo “peer-reviewed”. Il capitolo 2 contiene le statistiche dei tempi di transito e dei tempi di residenza nelle acque superficiali del mare Adriatico, calcolate a partire da dati derivati da drifters reali e da simulazioni numeriche lagrangiane. I risultati delle analisi dei dati derivati da drifters reali mostrano in generale una sottostima delle statistiche, data la loro vita operativa relativamente breve (tempo di emivita di 40 giorni). Il problema del bias introdotto dalla breve vita dei drifters reali può essere risolto mediante l’integrazione delle traiettorie di particelle sintetiche su scale temporali più lunghe (750 giorni) con opportuni modelli di circolazione. Per particelle numeriche rilasciate nella parte più settentrionale del bacino Adriatico, l’utilizzo di un modello statistico di avvezione-diffusione associato alla circolazione superficiale fornisce una scala temporale di 216-260 giorni come tempo massimo di uscita dal bacino stesso. In modo analogo, per una particella in ingresso nel bacino nel lato orientale dello stretto di Otranto, il modello fornisce una stima di 170-185 giorni per il tempo di uscita lungo il lato occidentale del canale di Otranto. I tempi di residenza nel bacino Adriatico sono dell’ordine dei 150-168 giorni. Nel Capitolo 3, si esaminano le caratteristiche della dispersione di coppie di drifters superficiali. La dispersione relativa (D2) è calcolata a partire daidrifters rilasciati nel mare Adriatico nel periodo 25 Agosto 1990 - 1 Gennaio 2007, sulla base di diverse separazioni iniziali (chan). La dispersione relativa viene interpretata sulla base di due distinti regimi. Il primo, è un regime non-locale che cambia secondo una funzione esponenziale del tempo; il secondo è un regime locale che può essere a sua volta diviso in tre regimi distinti: Richardson-type ( D 2 ∼ t 3 ), ballistic-type ( D 2 ∼ t 2 ), ediffusive-type ( D 2 ∼ t ) . A partire dalle traiettorie dei drifters, sono state inoltre ricavate le distribuzioni dei campi degli esponenti di Lyapunov a scala-finita (finite-scale Lyapunov exponent,FLSE) e le LagrangianStructureFunctions (LSF), che descrivono le proprietà fisiche intrinseche ad una data scala per il mare Adriatico. Sono stati inoltre calcolati altri descrittori statistici (angolo medio di separazione, displacementkurtosis, e dispersione assoluta) che permettono una conoscenza più dettagliata del comportamento delle coppie di drifters. Il capitolo 4 descrive quantitativamente le proprietà del trasporto superficiale – quasi superficiale e i tempi di residenza nel bacino Adriatico mediante l’uso di drifters sintetici. Le traiettorie analizzate sono state calcolate a partire dai campi di velocità giornalieri forniti dal modello di circolazione generale MITgcm, per il periodo Ottobre 2006 – Dicembre 2008. La traiettoria di ciascun drifter sintetico è stata ottenuta integrando e interpolando i campi dicorrente utilizzando, rispettivamente, uno schema Runge-Kutta del quarto ordine eun’interpolazione bilineare. In particolare, si confrontano le caratteristiche della circolazione in due anni caratterizzati da condizioni climatiche contrastanti (2007, anno con inverno “mite” e autunno “freddo”; 2008, anno con inverno “normale” ed estate “calda”). Al fine di validare il modello, le statistiche di trasporto delle particelle numeriche attraverso tre sezioni lungo la costa italiana (Conero, Gargano e stretto di Otranto) sono state confrontate con le statistiche di trasporto derivate dai drifter per le stesse aree. Il confronto suggerisce che le particelle numeriche derivate dalla simulazione modellistica sono più lente delle particelle “simulate” a partire da campi di velocità dei drifters, a causa di un minor contenuto energetico delle correnti stimate dal modello numerico nel periodo considerato. Dopo aver filtrato la circolazione euleriana media, per le particelle numeriche sono state inoltre calcolate le statistiche lagrangiane per l’intero bacino adriatico. I valori di momento angolare medio, di diffusività e covarianza sono minori di quanto calcolato dai drifter reali; tuttavia la scala temporale integrale lagrangiana rimane invariata tra i due casi (simulazioni - dati reali). Al fine di tener conto delle diverse condizioni meteorologiche tra il 2007 e 2008, e dei diversi valori di energia cinetica del campo medio (MKE = 3.14e-009 , 3.79e-009 km 2 s −2 per gli anni 2007 e 2008) e di energia cinetica media del campo turbolento (2.56e-009 e 2.90e-009 km 2 s −2 per il 2007 e 2008), e per valutare la loro diversa influenza sulle scale di trasporto, i confronti e le statistiche sono state ripetute separatamente per i due anni. I risultati ottenuti hanno evidenziato che gli effetti della circolazione indotta dal vento a nord del fiume Po, e gli effetti dell’apporto di acque dolci del fiume Po stesso sulla circolazione superficiale determinano valori confrontabili dei tempi di residenza (182 giorni e 185 giorni) per i due anni esaminati (2007 e 2008). Description: 2011/2012 Type: Tesi di dottorato; Doctoral Thesis 2013-04-23T00:00:00Z http://www.openstarts.units.it:80/dspace/handle/10077/8643 Title: Modellazione e analisi non lineare di strutture in muratura e in legno Authors: Rinaldin, Giovanni Abstract: Il lavoro di ricerca svolto ha avuto per oggetto lo studio e la definizione di modelli di calcolo e di verifica non lineari per strutture in muratura e legno. Lo studio ha riguardato principalmente la modellazione fenomenologica del comportamento non lineare di pannelli murari e connettori metallici per legno, e la conseguente implementazione delle leggi isteretiche analizzate in molle utente in diversi solutori commerciali/open-source (ABAQUS e OpenSees). In particolare il lavoro svolto si è articolato nei seguenti punti: 1. strutture in muratura: a. studio del comportamento meccanico del materiale muratura e raccolta dati di letteratura b. studio di modelli isteretici fenomenologici adatti alla rappresentazione degli edifici in muratura c. scrittura di routine in Fortran per implementare il modello fenomenologico studiato (modello di Tomazevic) nel programma ABAQUS d. analisi non lineare di edifici in muratura già studiati in letteratura ai fini della validazione del modello proposto e. studio di metodi semplificati per il calcolo sismico di edifici a schiera e scrittura di programma autonomo per l’analisi di pushover con il metodo RigStripTS f. prove di laboratorio sul comportamento di fasce murarie g. studio del modello in mesoscala incluso nel programma ADAPTIC e costruzione di modelli per la simulazione di test sperimentali tramite un meshatore sviluppato ad-hoc. 2. strutture in legno: a. studio del comportamento meccanico dei connettori metallici per le costruzioni in legno lamellare incrociato (X-lam) e a telaio leggero e raccolta dati di letteratura b. studio del comportamento ortotropo del materiale e scrittura routine per materiale utente per la sua rappresentazione in campo elastico c. formulazione di legame isteretico fenomenologico per connettori metallici (squadrette, hold-down e unioni con viti/chiodi) d. scrittura di routine in Fortran per implementare il modello fenomenologico proposto nei programmi ABAQUS e OpenSees e. studio di prove sperimentali su connettori singoli e analisi e implementazione di un metodo di taratura adatto al modello utilizzato, implementato in un programma scritto ad-hoc f. analisi di modelli non-lineari di test di laboratorio eseguiti su pareti in X-lam e a telaio leggero e di test pseudo-dinamici su edifici in X-lam monopiano, condotte con il modello sviluppato g. estensione del modello ad altre tipologie di edifici con struttura lignea. Entrambe le tipologie strutturali sono state trattate con modelli fenomenologici in modo da fornire uno strumento di calcolo non lineare per analisi statiche e dinamiche compatibile con le capacità di calcolo attuali; è stata indagata inoltre la possibilità di distribuire l’onere computazionale a più computer parallelizzando l’analisi strutturale. Particolare attenzione è stata riservata inoltre alla capacità sismica delle strutture analizzate, stimando in modo rigoroso l’energia dissipata in ogni analisi. A tal fine è stato eseguito anche uno studio sui modelli isteretici più comuni atto a verificare l’attendibilità di metodi di verifica globale in campo statico non lineare, come il metodo N2 introdotto da Fajfar e il metodo dello spettro sovra-smorzato proposto dalla norma ATC-40 statunitense. Tale studio è stato condotto tramite la scrittura di un programma ad-hoc per il calcolo di spettri inelastici a partire da record sismici di qualsiasi natura (registrati, artificiali o generati). I risultati ottenuti con tale strumento sono stati utilizzati per validare i metodi sopracitati e verificare la loro attendibilità al variare del periodo strutturale e della capacità dissipativa delle diverse tipologie strutturali. Strutture in muratura La necessità di uno strumento di calcolo flessibile ma potente ha indotto lo sviluppo di un elemento molla che implementasse un ciclo isteretico rappresentativo del comportamento del materiale. È stato scelto dalla letteratura disponibile la legge proposta da Tomazevic e Lutman nel 1996 formulata sulla scorta di numerose prove sperimentali condotte su maschi murari. Tale legge lineare a tratti si compone di una curva scheletro trilineare e di percorsi di scarico/ricarico, oltre che gli opportuni degradi di rigidezza e resistenza necessari per poter caratterizzare correttamente il comportamento del materiale muratura vicino alle condizioni di collasso. La routine implementata utilizza un set di parametri che regolano le rigidezze del sistema e l’ampiezza dei cicli, mentre la resistenza di ciascun elemento molla è calcolata automaticamente sulla base del valore di sforzo normale letto durante l’analisi. Tale operazione è svolta in modo disaccoppiato (cioè è riferita allo sforzo assiale letto al passo precedente durante l’analisi) e il valore resistente cambia sulla base del meccanismo di rottura, che il modello seleziona automaticamente. Tale strumento è stato utilizzato con successo per l’analisi statica non lineare della Parete D, facente parte del progetto Catania, progetto nato con l’intento di valutare la vulnerabilità sismica di edifici storici, poi con analisi dinamiche incrementali sulla parete stessa e su altre testate da Magenes a Pavia, che hanno dimostrato tutta la potenza del modello sviluppato, permettendo di stimare correttamente la capacità di dissipazione isteretica. Il metodo è stato esteso anche alle fasce di piano. Infine, il lavoro sul lato della modellazione fenomenologica è culminato con la formulazione di un macro-elemento, non ancora implementato, che permette di condensare le diverse molle utilizzate per ogni pannello in un solo elemento asta. I risultati con tali strumenti avanzati sono poi stati traslati in metodi e modelli per l’applicazione professionale, come il metodo RigStripTS (Amadio et al.) per il push-over di edifici in muratura e il metodo delle celle per lo studio sismico semplificato di edifici a schiera. L’analisi è stata approfondita con l’utilizzo del un modello a mesoscala per la modellazione dei pannelli murari, sviluppato all’Imperial College e disponibile nel software ADAPTIC (Macorini e Izzuddin). Tramite un meshatore sviluppato ad-hoc per la scrittura dei modelli, sono stati studiati diversi modelli in regime monotono. Strutture in legno Lo studio delle strutture in legno si è concentrato prevalentemente sul Cross-Lam (X-lam) e sul telaio leggero, mancando nelle normative internazionali delle linee guida per il calcolo sismico di tali strutture. Con lo scopo di stimare il comportamento strutturale in campo elastico e plastico e partendo dall’ipotesi, già documentata in letteratura, che tale tipologia costruttiva dissipa energia per i soli connettori, è stato formulato e implementato un modello isteretico lineare a tratti che rappresenta il comportamento ciclico in condizioni statiche e dinamiche di connettori metallici per pannelli in X-lam e per pareti a telaio leggero. Tale legame è stato implementato in elementi multi-molla per rappresentare ogni grado di libertà associato al connettore. Caratteristiche come degrado di resistenza, di rigidezza, effetto di attrito e dominio di resistenza sono state proposte e implementate. Tale modello necessita di una taratura preliminare da svolgersi su prove di laboratorio cicliche di connettori singoli; è stata proposta una metodologia per la taratura ciclica e implementata in un software autonomo sviluppato ad-hoc. Con tale strumento sono stati riprodotti fedelmente test ciclici di laboratorio condotti su singole pareti in X-lam in diverse configurazioni, giungendo poi alla validazione tridimensionale del modello eseguita sull’edificio monopiano testato pseudo-dinamicamente dal CNR-IVALSA a Trento. È stato sviluppato un programma per la meshatura automatica di modelli tridimensionali in X-lam. Analogo approccio è stato adottato per la tipologia a telaio leggero. Infine, è stata sviluppata una routine per il materiale legno che utilizza il modello ortotropo. La routine lavora in campo elastico ma può essere estesa al campo plastico mediante l’aggiunta del dominio resistente di Tsai e Wu, che è stato studiato si dimostra particolarmente adatto al materiale legno. Description: 2011/2012 Type: Tesi di dottorato; Doctoral Thesis 2013-04-22T00:00:00Z