Meccanismi di formazione dei tubi di lava

Two analytical models are proposed to describe two different mechanisms of lava tubes formation. A first model is introduced to describe the development of a solid crust in the central region of the channel, and the formation of a tube when crust widens until it reaches the leve\'es. The Newtonian assumption is considered and the steady state Navier- Stokes equation in a rectangular conduit is solved. A constant heat flux density assigned at the upper flow surface resumes the combined effects of two thermal processes: radiation and convection into the atmosphere. Advective terms are also included, by the introduction of velocity into the expression of temperature. Velocity is calculated as an average value over the channel width, so that lateral variations of temperature are neglected. As long as the upper flow surface cools, a solid layer develops, described as a plastic body, having a resistance to shear deformation. If the applied shear stress exceeds this resistance, crust breaks, otherwise, solid fragments present at the flow surface can weld together forming a continuous roof, as it happens in the sidewall flow regions. Variations of channel width, ground slope and effusion rate are analyzed, as parameters that strongly affect the shear stress values. Crust growing is favored when the channel widens, and tube formation is possible when the ground slope or the effusion rate reduce. A comparison of results is successfully made with data obtained from the analysis of pictures of actual flows. The second model describes the formation of a stable, well defined crust along both channel sides, their growing towards the center and their welding to form the tube roof. The fluid motion is described as in the model above. Thermal budget takes into account conduction into the atmosphere, and advection is included considering the velocity depending both on depth and channel width. The solidified crust has a non uniform thickness along the channel width. Stresses acting on the crust are calculated using the equations of the elastic thin plate, pinned at its ends. The model allows to calculate the distance where crust thickness is able to resist the drag of the underlying fluid and to sustain its weight by itself, and the level of the fluid can lower below the tube roof. Viscosity and thermal conductivity have been experimentally investigated through the use of a rotational viscosimeter. Analyzing samples coming from Mount Etna (2002) the following results have been obtained: the fluid is Newtonian and the thermal conductivity is constant in a range of temperature above the liquidus. For lower temperature, the fluid becomes non homogeneous, and the used experimental techniques are not able to detect any properties, because measurements are not reproducible.

Effetti idrodinamici prodotti dal costume tecnico sul nuotatore

L’introduzione dei costumi tecnici nel nuoto ha portato miglioramenti senza precedenti sulla prestazione. I miglioramenti nella velocità di nuoto sono stati attribuiti dalla letteratura a riduzioni nelle resistenze idrodinamiche sul nuotatore. Tuttavia, gli effetti specifici dovuti all’utilizzo di questo tipo di costume non sono ancora completamente chiariti. Questa tesi aveva l’obiettivo di indagare gli effetti del costume tecnico sul galleggiamento statico, sulla posizione del corpo e sulla resistenza idrodinamica in avanzamento passivo. Nello studio preliminare sono stati misurati la spinta idrostatica, i volumi polmonari dinamici e la circonferenza toracica di 9 nuotatori che indossavano un costume tradizionale o un costume tecnico in gomma sintetica. Indossare il costume tecnico ha determinato una riduzione significativa del galleggiamento statico, e la compressione toracica causata da questo tipo di costume potrebbe avere una relazione con la significativa riduzione dei volumi polmonari misurati quando il nuotatore indossa questo tipo di costume. Un successiva analisi prevedeva il traino passivo di 14 nuotatori che mantenevano la miglior posizione idrodinamica di scivolamento indossando un costume tradizionale, tecnico in tessuto e tecnico in gomma. La posizione del corpo in avanzamento è stata misurata con un’analisi cinematica. La resistenza passiva indossando i costumi tecnici è risultata significativamente minore per entrambi i costumi tecnici rispetto alla prova con costume tradizionale. L’analisi condotta attraverso modelli di regressione lineari ha mostrato che una parte della riduzione della resistenza passiva era legata a proprietà intrinseche dei costumi tecnici. Tuttavia, anche l’area di impatto frontale determinata dall’inclinazione del tronco del soggetto in scivolamento e l’inclinazione degli arti inferiori hanno mostrato una marcata influenza sulla resistenza idrodinamica passiva. Pertanto, la riduzione di resistenza idrodinamica durante lo scivolamento passivo effettuato con costume tecnico da nuoto è attribuibile, oltre all’effetto del materiale di composizione del costume, ad una variazione della posizione del corpo del nuotatore.

Multilevel Domain Decomposition Algorithms for Monolithic Fluid-Structure Interaction Problems with Application to Haemodynamics

Finite element techniques for solving the problem of fluid-structure interaction of an elastic solid material in a laminar incompressible viscous flow are described. The mathematical problem consists of the Navier-Stokes equations in the Arbitrary Lagrangian-Eulerian formulation coupled with a non-linear structure model, considering the problem as one continuum. The coupling between the structure and the fluid is enforced inside a monolithic framework which computes simultaneously for the fluid and the structure unknowns within a unique solver. We used the well-known Crouzeix-Raviart finite element pair for discretization in space and the method of lines for discretization in time. A stability result using the Backward-Euler time-stepping scheme for both fluid and solid part and the finite element method for the space discretization has been proved. The resulting linear system has been solved by multilevel domain decomposition techniques. Our strategy is to solve several local subproblems over subdomain patches using the Schur-complement or GMRES smoother within a multigrid iterative solver. For validation and evaluation of the accuracy of the proposed methodology, we present corresponding results for a set of two FSI benchmark configurations which describe the self-induced elastic deformation of a beam attached to a cylinder in a laminar channel flow, allowing stationary as well as periodically oscillating deformations, and for a benchmark proposed by COMSOL multiphysics where a narrow vertical structure attached to the bottom wall of a channel bends under the force due to both viscous drag and pressure. Then, as an example of fluid-structure interaction in biomedical problems, we considered the academic numerical test which consists in simulating the pressure wave propagation through a straight compliant vessel. All the tests show the applicability and the numerical efficiency of our approach to both two-dimensional and three-dimensional problems.