CFD Applications for Wave Energy Conversion Devices (MoonWEC) and Turbulent Fountains for Environmental Fluid Mechanics

This thesis is based on two studies that are related to floating wave energy conversion (WEC) devices and turbulent fountains. The ability of the open-source CFD software OpenFOAM® has been studied to simulate these phenomena. The CFD model has been compared with the physical experimental results. The first study presents a model of a WEC device, called MoonWEC, which is patented by the University of Bologna. The CFD model of the MoonWEC under the action of waves has been simulated using OpenFOAM and the results are promising. The reliability of the CFD model is confirmed by the laboratory experiments, conducted at the University of Bologna, for which a small-scale prototype of the MoonWEC was made from wood and brass. The second part of the thesis is related to the turbulent fountains which are formed when a heavier source fluid is injected upward into a lighter ambient fluid, or else a lighter source fluid is injected downward into a heavier ambient fluid. For this study, the first case is considered for laboratory experiments and the corresponding CFD model. The vertical releases of the source fluids into a quiescent, uniform ambient fluid, from a circular source, were studied with different densities in the laboratory experiments, conducted at the University of Parma. The CFD model has been set up for these experiments. Favourable results have been observed from the OpenFOAM simulations for the turbulent fountains as well, indicating that it can be a reliable tool for the simulation of such phenomena.

Hidden Schrödinger symmetry in FLRW and Bianchi I minisuperspace models of cosmology and black holes Schwarzschild dynamics

A broad sector of literature focuses on the relationship between fluid dynamics and gravitational systems. This thesis presents results that suggest the existence of a new kind of fluid/gravity duality not based on the holographic principle. The goal is to provide tools that allow us to systematically unearth hidden symmetries for reduced models of cosmology. The focus is on the field space of these models, i.e. the superspace. In fact, conformal isometries of the supermetric leave geodesics in the field space unaltered; this leads to symmetries of the models. An innovative aspect is the use of the Eisenhart-Duval’s lift. Using this method, systems constrained by a potential can be treated as free ones. Moreover, charges explicitly dependent on time, i.e. dynamical, can be found. A detailed analysis is carried out on three basic models of homogenous cosmology: i) flat Friedmann-Lemaître-Robertson-Walker’s isotropic universe filled with a massless scalar field; ii) Schwarzschild’s black hole mechanics and its extension to vacuum (A)dS gravity; iii) Bianchi’s anisotropic type I universe with a massless scalar field. The results show the presence of a hidden Schrödinger’s symmetry which, being intrinsic to both Navier-Stokes’ and Schrödinger’s equations, indicates a correspondence between cosmology and hydrodynamics. Furthermore, the central extension of this algebra explicitly relates two concepts. The first is the number of particles coming from the fluid picture; while the second is the ratio between the IR and UV cutoffs that weighs how much a theory has of “classical” over “quantum”. This suggests a spacetime that emerges from an underlying world which is described by quantum building blocks. These quanta statistically conspire to appear as gravitational phenomena from a macroscopic point of view.

Fluidodinâmica de sementes de brócolos em leito fluizado e leito de jorro

The spouted and fluidized bed technologies are usually employed in operations of drying, coating and granulation of particles by the chemical and pharmaceutical industries. The use of these techniques in agronomy is limited to the treatment of some species of seeds. In this work, the objective was to analyse the fluid-dynamics of fluidized and spouted beds when broccoli (Brassica oleracea L. var. Italica) seeds are used and also to verify the influence on seed germination after 60 min of seed exposition to spouting or fluidization, at room temperature. The fluid-dynamics was defined by the measurements of the bed pressure drop as a function of the air flow rate for different seeds loads. The experimental conditions were based on the physical properties of the seeds and were limited by the apparatus dimensions. The cone-cylindrical bed was constructed in plexyglass to permit flow visualization. The values of the parameters: maximum pressure drop, minimum spouting flow rate and pressure drop, and stable spout pressure drop were experimentally obtained from the fluid-dynamic analysis and were compared with the values calculated by empirical equations found in the literature. The same procedure was carried out with the fluidized bed and the important parameters for this regime were the air velocity and the bed pressure drop at minimum fluidization. The analysis of seed germination indicated that no damage was caused to the seeds by the spout or fluidization processes.

Metamorfose planar via métodos Level Set e Particle Level Set para a reconstrução de superfícies tridimensionais

Evolving interfaces were initially focused on solutions to scientific problems in Fluid Dynamics. With the advent of the more robust modeling provided by Level Set method, their original boundaries of applicability were extended. Specifically to the Geometric Modeling area, works published until then, relating Level Set to tridimensional surface reconstruction, centered themselves on reconstruction from a data cloud dispersed in space; the approach based on parallel planar slices transversal to the object to be reconstructed is still incipient. Based on this fact, the present work proposes to analyse the feasibility of Level Set to tridimensional reconstruction, offering a methodology that simultaneously integrates the proved efficient ideas already published about such approximation and the proposals to process the inherent limitations of the method not satisfactorily treated yet, in particular the excessive smoothing of fine characteristics of contours evolving under Level Set. In relation to this, the application of the variant Particle Level Set is suggested as a solution, for its intrinsic proved capability to preserve mass of dynamic fronts. At the end, synthetic and real data sets are used to evaluate the presented tridimensional surface reconstruction methodology qualitatively.

Metamorfose planar via métodos Level Set e Particle Level Set para a reconstrução de superfícies tridimensionais

Evolving interfaces were initially focused on solutions to scientific problems in Fluid Dynamics. With the advent of the more robust modeling provided by Level Set method, their original boundaries of applicability were extended. Specifically to the Geometric Modeling area, works published until then, relating Level Set to tridimensional surface reconstruction, centered themselves on reconstruction from a data cloud dispersed in space; the approach based on parallel planar slices transversal to the object to be reconstructed is still incipient. Based on this fact, the present work proposes to analyse the feasibility of Level Set to tridimensional reconstruction, offering a methodology that simultaneously integrates the proved efficient ideas already published about such approximation and the proposals to process the inherent limitations of the method not satisfactorily treated yet, in particular the excessive smoothing of fine characteristics of contours evolving under Level Set. In relation to this, the application of the variant Particle Level Set is suggested as a solution, for its intrinsic proved capability to preserve mass of dynamic fronts. At the end, synthetic and real data sets are used to evaluate the presented tridimensional surface reconstruction methodology qualitatively.

Propriedades dinâmicas de fluidos por simulação computacional: métodos híbridos atomístico-contínuo

Computational methods for the calculation of dynamical properties of fluids might consider the system as a continuum or as an assembly of molecules. Molecular dynamics (MD) simulation includes molecular resolution, whereas computational fluid dynamics (CFD) considers the fluid as a continuum. This work provides a review of hybrid methods MD/CFD recently proposed in the literature. Theoretical foundations, basic approaches of computational methods, and dynamical properties typically calculated by MD and CFD are first presented in order to appreciate the similarities and differences between these two methods. Then, methods for coupling MD and CFD, and applications of hybrid simulations MD/CFD, are presented.

Verification of a mixed high-order accurate DNS code for laminar turbulent transition by the method of manufactured solutions

This paper presents results on a verification test of a Direct Numerical Simulation code of mixed high-order of accuracy using the method of manufactured solutions (MMS). This test is based on the formulation of an analytical solution for the Navier-Stokes equations modified by the addition of a source term. The present numerical code was aimed at simulating the temporal evolution of instability waves in a plane Poiseuille flow. The governing equations were solved in a vorticity-velocity formulation for a two-dimensional incompressible flow. The code employed two different numerical schemes. One used mixed high-order compact and non-compact finite-differences from fourth-order to sixth-order of accuracy. The other scheme used spectral methods instead of finite-difference methods for the streamwise direction, which was periodic. In the present test, particular attention was paid to the boundary conditions of the physical problem of interest. Indeed, the verification procedure using MMS can be more demanding than the often used comparison with Linear Stability Theory. That is particularly because in the latter test no attention is paid to the nonlinear terms. For the present verification test, it was possible to manufacture an analytical solution that reproduced some aspects of an instability wave in a nonlinear stage. Although the results of the verification by MMS for this mixed-order numerical scheme had to be interpreted with care, the test was very useful as it gave confidence that the code was free of programming errors. Copyright (C) 2009 John Wiley & Sons, Ltd.

3D CFD Simulation of Vortex-induced Vibration of Cylinder

The 3D flow around a circular cylinder free to oscillate transversely to the free stream was simulated using Computational Fluid Dynamics (CFD) and the Spalart-Allmaras Detached Eddy Simulation (DES) turbulence model for a Reynolds number Re = 10(4). Simulations were carried out for a small mass-damping parameter m*zeta = 0.00858, where m* = 3.3 and zeta = 0.0026. We found good agreement between the numerical results and experimental data. The simulations predicted the high observed amplitudes of the upper branch of vortex-induced vibrations for low mass-damping parameters.

Numerical prediction of gas concentrations and fluctuations above a triangular hill within a turbulent boundary layer

The objective of the present work is to propose a numerical and statistical approach, using computational fluid dynamics, for the study of the atmospheric pollutant dispersion. Modifications in the standard k-epsilon turbulence model and additional equations for the calculation of the variance of concentration are introduced to enhance the prediction of the flow field and scalar quantities. The flow field, the mean concentration and the variance of a flow over a two-dimensional triangular hill, with a finite-size point pollutant source, are calculated by a finite volume code and compared with published experimental results. A modified low Reynolds k-epsilon turbulence model was employed in this work, using the constant of the k-epsilon model C(mu)=0.03 to take into account the inactive atmospheric turbulence. The numerical results for the velocity profiles and the position of the reattachment point are in good agreement with the experimental results. The results for the mean and the variance of the concentration are also in good agreement with experimental results from the literature. (C) 2009 Elsevier Ltd. All rights reserved.

Galerkin finite element method for incompressible thermofluid flows framed within the bond graph theory

In this paper a bond graph methodology is used to model incompressible fluid flows with viscous and thermal effects. The distinctive characteristic of these flows is the role of pressure, which does not behave as a state variable but as a function that must act in such a way that the resulting velocity field has divergence zero. Velocity and entropy per unit volume are used as independent variables for a single-phase, single-component flow. Time-dependent nodal values and interpolation functions are introduced to represent the flow field, from which nodal vectors of velocity and entropy are defined as state variables. The system for momentum and continuity equations is coincident with the one obtained by using the Galerkin method for the weak formulation of the problem in finite elements. The integral incompressibility constraint is derived based on the integral conservation of mechanical energy. The weak formulation for thermal energy equation is modeled with true bond graph elements in terms of nodal vectors of temperature and entropy rates, resulting a Petrov-Galerkin method. The resulting bond graph shows the coupling between mechanical and thermal energy domains through the viscous dissipation term. All kind of boundary conditions are handled consistently and can be represented as generalized effort or flow sources. A procedure for causality assignment is derived for the resulting graph, satisfying the Second principle of Thermodynamics. (C) 2007 Elsevier B.V. All rights reserved.

Investigation of a method for predicting AUV derivatives

The paper addresses the problem of autonomous underwater vehicle (AUV) modelling and parameter estimation as a means to predict the dynamic performance of underwater vehicles and thus provide solid guidelines during their design phase. The use of analytical and semi-empirical (ASE) methods to estimate the hydrodynamic derivatives of a popular class of AUVs is discussed. A comparison is done with the results obtained by using computational fluid dynamics to evaluate the bare hull lift force distribution around a fully submerged body. An application is made to the estimation of the hydrodynamic derivatives of the MAYA AUV, an autonomous underwater vehicle developed under a joint Indian-Portuguese project. The estimates obtained were used to predict the turning diameter of the vehicle during sea trials. (C) 2008 Elsevier Ltd. All rights reserved.

Gas flow analysis in a Kraft recovery boiler

Demands for optimal boiler performance and increased concerns in lowering emission have always been the driving force in the reevaluation and evolution of the Kraft boiler: specifically the air distribution strategies that are directly related to achieving increased residence time of flue gas combustion inside the furnace which in turn lowers atmosphere emission levels and enhances boiler operation. This paper presents the results of a study that analyzes the interaction of the different multilevel air injections have on flue gas flow patterns including various quaternary air supply arrangements. Additionally, this study assesses the performance of the CFD (Computational Fluid Dynamics) model against data available in literature. Simulations were performed considering isothermal and incompressible flows, and did not take into account thermal phenomena or chemical reactions. The numerical solutions generated proved to be coherently related to the data available in literature, and provided proof of the efficiency of tertiary level air injection, as well as revealed that quaternary air injection ports arranged in a symmetrical configuration is most suitable for optimal equipment operation. (C) 2010 Elsevier B.V. All rights reserved.

DEVELOPMENT OF A MICRO-HEAT EXCHANGER WITH STACKED PLATES USING LTCC TECHNOLOGY

A green ceramic tape micro-heat exchanger was developed using Low Temperature Co-fired Ceramics technology (LTCC). The device was designed by using Computational Aided Design software and simulations were made using a Computational Fluid Dynamics package (COMSOL Multiphysics) to evaluate the homogeneity of fluid distribution in the microchannels. Four geometries were proposed and simulated in two and three dimensions to show that geometric details directly affect the distribution of velocity in the micro-heat exchanger channels. The simulation results were quite useful for the design of the microfluidic device. The micro-heat exchanger was then constructed using the LTCC technology and is composed of five thermal exchange plates in cross-flow arrangement and two connecting plates, with all plates stacked to form a device with external dimensions of 26 x 26 x 6 mm(3).

The effect of flow arrangement on the pressure drop of plate heat exchangers

For the optimal design of plate heat exchangers (PHEs), an accurate thermal-hydraulic model that takes into account the effect of the flow arrangement on the heat load and pressure drop is necessary. In the present study, the effect of the flow arrangement on the pressure drop of a PHE is investigated. Thirty two different arrangements were experimentally tested using a laboratory scale PHE with flat plates. The experimental data was used for (a) determination of an empirical correlation for the effect of the number of passes and number of flow channels per pass on the pressure drop; (b) validation of a friction factor model through parameter estimation; and (c) comparison with the simulation results obtained with a CFD (computational fluid dynamics) model of the PHE. All three approaches resulted in a good agreement between experimental and predicted values of pressure drop. Moreover, the CFD model is used for evaluating the flow maldistribution in a PHE with two channels Per Pass. (c) 2008 Elsevier Ltd. All rights reserved.

Investigation of Normal Force and Moment Coefficients for an AUV at Nonlinear Angle of Attack and Sideslip Range

This paper presents a comparative study of computational fluid dynamics (CFD) and analytical and semiempirical (ASE) methods applied to the prediction of the normal force and moment coefficients of an autonomous underwater vehicle (AUV). Both methods are applied to the. bare hull of the vehicle and to the body-hydroplane combination. The results are validated through experiments in a towing tank. It is shown that the CFD approach allows for a good prediction of the coefficients over the range of angles of attack considered. In contrast with the traditional ASE formulations used in naval and aircraft fields, an improved methodology is introduced that takes advantage of the qualitative information obtained from CFD flow visualizations.