Modelling and control of a floating oscillating water column

A novel numerical model of a Bent Backwards Duct Buoy (BBDB) Oscillating Water Column (OWC) Wave Energy Converter was created based on existing isolated numerical models of the different energy conversion systems utilised by an OWC. The novel aspect of this numerical model is that it incorporates the interdependencies of the different power conversion systems rather than modelling each system individually. This was achieved by accounting for the dynamic aerodynamic damping caused by the changing turbine rotational velocity by recalculating the turbine damping for each simulation sample and applying it via a feedback loop. The accuracy of the model was validated using experimental data collected during the Components for Ocean Renewable Energy Systems (CORES) EU FP-7 project that was tested in Galway Bay, Ireland. During the verification process, it was discovered that the model could also be applied as a valuable tool when troubleshooting device performance. A new turbine was developed and added to a full scale model after being investigated using Computational Fluid Dynamics. The energy storage capacity of the impulse turbine was investigated by modelling the turbine with both high and low inertia and applying three turbine control theories to the turbine using the full scale model. A single Maximum Power Point Tracking algorithm was applied to the low-inertia turbine, while both a fixed and dynamic control algorithm was applied to the high-inertia turbine. These results suggest that the highinertia turbine could be used as a flywheel energy storage device that could help minimize output power variation despite the low operating speed of the impulse turbine. This research identified the importance of applying dynamic turbine damping to a BBDB OWC numerical model, revealed additional value of the model as a device troubleshooting tool, and found that an impulse turbine could be applied as an energy storage system.

Energy conversion: a comparison of fix- and self-referenced wave energy converters

The paper presents an investigation of fix-referenced and self-referenced wave energy converters and a comparison of their corresponding wave energy conversion capacities from real seas. For conducting the comparisons, two popular wave energy converters, point absorber and oscillating water column, and their power conversion capacities in the fixed-referenced and self-referenced forms have been numerically studied and compared. In the numerical models, the deviceâ s power extractions from seas are maximized using the correspondingly optimized power take-offs in different sea states, thus their power conversion capacities can be calculated and compared. From the comparisons and analyses, it is shown that the energy conversion capacities of the self-referenced devices can be significantly increased if the motions of the device itself can be utilized for wave energy conversion; and the self-referenced devices can be possibly designed to be compliant in long waves, which could be a very beneficial factor for device survivability in the extreme wave conditions (normally long waves). In this regards, the self-referenced WECs (wave energy converters) may be better options in terms of wave energy conversion from the targeted waves in seas (frequently the most occurred), and in terms of the device survivability, especially in the extreme waves when compared to the fix-referenced counterparts.

Dynamic Modeling and Seasonal Performance Evaluation of Air-to-Water Heat Pump Systems

In this Thesis a series of numerical models for the evaluation of the seasonal performance of reversible air-to-water heat pump systems coupled to residential and non-residential buildings are presented. The exploitation of the energy saving potential linked to the adoption of heat pumps is a hard task for designers due to the influence on their energy performance of several factors, like the external climate variability, the heat pump modulation capacity, the system control strategy and the hydronic loop configuration. The aim of this work is to study in detail all these aspects. In the first part of this Thesis a series of models which use a temperature class approach for the prediction of the seasonal performance of reversible air source heat pumps are shown. An innovative methodology for the calculation of the seasonal performance of an air-to-water heat pump has been proposed as an extension of the procedure reported by the European standard EN 14825. This methodology can be applied not only to air-to-water single-stage heat pumps (On-off HPs) but also to multi-stage (MSHPs) and inverter-driven units (IDHPs). In the second part, dynamic simulation has been used with the aim to optimize the control systems of the heat pump and of the HVAC plant. A series of dynamic models, developed by means of TRNSYS, are presented to study the behavior of On-off HPs, MSHPs and IDHPs. The main goal of these dynamic simulations is to show the influence of the heat pump control strategies and of the lay-out of the hydronic loop used to couple the heat pump to the emitters on the seasonal performance of the system. A particular focus is given to the modeling of the energy losses linked to on-off cycling.

Resonant metamaterials for the control of Rayleigh waves

Since their emergence, locally resonant metamaterials have found several applications for the control of surface waves, from micrometer-sized electronic devices to meter-sized seismic barriers. The interaction between Rayleigh-type surface waves and resonant metamaterials has been investigated through the realization of locally resonant metasurfaces, thin elastic interfaces constituted by a cluster of resonant inclusions or oscillators embedded near the surface of an elastic waveguide. When such resonant metasurfaces are embedded in an elastic homogeneous half-space, they can filter out the propagation of Rayleigh waves, creating low-frequency bandgaps at selected frequencies. In the civil engineering context, heavy resonating masses are needed to extend the bandgap frequency width of locally resonant devices, a requirement that limits their practical implementations. In this dissertation, the wave attenuation capabilities of locally resonant metasurfaces have been enriched by proposing (i) tunable metasurfaces to open large frequency bandgaps with small effective inertia, and by developing (ii) an analytical framework aimed at studying the propagation of Rayleigh waves propagation in deep resonant waveguides. In more detail, inertial amplified resonators are exploited to design advanced metasurfaces with a prescribed static and a tunable dynamic response. The modular design of the tunable metasurfaces allows to shift and enlarge low-frequency spectral bandgaps without modifying the total inertia of the metasurface. Besides, an original dispersion law is derived to study the dispersive properties of Rayleigh waves propagating in thick resonant layers made of sub-wavelength resonators. Accordingly, a deep resonant wave barrier of mechanical resonators embedded inside the soil is designed to impede the propagation of seismic surface waves. Numerical models are developed to confirm the analytical dispersion predictions of the tunable metasurface and resonant layer. Finally, a medium-size scale resonant wave barrier is designed according to the soil stratigraphy of a real geophysical scenario to attenuate ground-borne vibration.

Design, development, and analysis of hybrid metal-composite tube for industrial application

In the aerospace, automotive, printing, and sports industries, the development of hybrid Carbon Fiber Reinforced Polymer (CFRP)-metal components is becoming increasingly important. The coupling of metal with CFRP in axial symmetric components results in reduced production costs and increased mechanical properties such as bending, torsional stiffness, mass reduction, damping, and critical speed compared to the single material-built ones. In this thesis, thanks to a novel methodology involving a rubbery/viscoelastic interface layer, several hybrid aluminum-CFRP prototype tubes were produced. Besides, an innovative system for the cure of the CFRP part has been studied, analyzed, tested, and developed in the company that financed these research activities (Reglass SRL, Minerbio BO, Italy). The residual thermal stresses and strains have been investigated with numerical models based on the Finite Element Method (FEM) and compared with experimental tests. Thanks to numerical models, it was also possible to reduce residual thermal stresses by optimizing the lamination sequence of CFRP and determining the influence of the system parameters. A novel software and methodology for evaluating mechanical and damping properties of specimens and tubes made in CFRP were also developed. Moreover, to increase the component's damping properties, rubber nanofibers have been produced and interposed throughout the lamination of specimens. The promising results indicated that the nanofibrous mat could improve the material damping factor over 77% and be adopted in CFRP components with a negligible increment of weight or losing mechanical properties.

Flow of Complex Fluids in Geological Fractures

In this study, the lubrication theory is used to model flow in geological fractures and analyse the compound effect of medium heterogeneity and complex fluid rheology. Such studies are warranted as the Newtonian rheology is adopted in most numerical models because of its ease of use, despite non-Newtonian fluids being ubiquitous in subsurface applications. Past studies on Newtonian and non-Newtonian flow in single rock fractures are summarized in Chapter 1. Chapter 2 presents analytical and semi-analytical conceptual models for flow of a shear-thinning fluid in rock fractures having a simplified geometry, providing a first insight on their permeability. in Chapter 3, a lubrication-based 2-D numerical model is first implemented to solve flow of an Ellis fluid in rough fractures; the finite-volumes model developed is more computationally effective than conducting full 3-D simulations, and introduces an acceptable approximation as long as the flow is laminar and the fracture walls relatively smooth. The compound effect of shear-thinning fluid nature and fracture heterogeneity promotes flow localization, which in turn affects the performance of industrial activities and remediation techniques. In Chapter 4, a Monte Carlo framework is adopted to produce multiple realizations of synthetic fractures, and analyze their ensemble statistics pertaining flow for a variety of real non-Newtonian fluids; the Newtonian case is used as a benchmark. In Chapter 5 and Chapter 6, a conceptual model of the hydro-mechanical aspects of backflow occurring in the last phase of hydraulic fracturing is proposed and experimentally validated, quantifying the effects of the relaxation induced by the flow.

Development and Implementation of the Plasma Focus Technology for Radiation Therapy Applications

A Plasma Focus device can confine in a small region a plasma generated during the pinch phase. When the plasma is in the pinch condition it creates an environment that produces several kinds of radiations. When the filling gas is nitrogen, a self-collimated backwardly emitted electron beam, slightly spread by the coulomb repulsion, can be considered one of the most interesting outputs. That beam can be converted into X-ray pulses able to transfer energy at an Ultra-High Dose-Rate (UH-DR), up to 1 Gy pulse-1, for clinical applications, research, or industrial purposes. The radiation fields have been studied with the PFMA-3 hosted at the University of Bologna, finding the radiation behavior at different operating conditions and working parameters for a proper tuning of this class of devices in clinical applications. The experimental outcomes have been compared with available analytical formalisms as benchmark and the scaling laws have been proposed. A set of Monte Carlo models have been built with direct and adjoint techniques for an accurate X-ray source characterization and for setting fast and reliable irradiation planning for patients. By coupling deterministic and Monte Carlo codes, a focusing lens for the charged particles has been designed for obtaining a beam suitable for applications as external radiotherapy or intra-operative radiation therapy. The radiobiological effectiveness of the UH PF DR, a FLASH source, has been evaluated by coupling different Monte Carlo codes estimating the overall level of DNA damage at the multi-cellular and tissue levels by considering the spatial variation effects as well as the radiation field characteristics. The numerical results have been correlated to the experimental outcomes. Finally, ambient dose measurements have been performed for tuning the numerical models and obtaining doses for radiation protection purposes. The PFMA-3 technology has been fully characterized toward clinical implementation and installation in a medical facility.

Early-stage deformation localisation in thrust systems: new perspectives from carbonates in the Italian Southern Alps and Oman Mountains

This thesis investigates mechanisms and boundary conditions that steer the early localisation of deformation and strain in carbonate multilayers involved in thrust systems, under shallow and mid-crustal conditions. Much is already understood about deformation localisation, but some key points remain loosely constrained. They encompass i) the understanding of which structural domains can preserve evidence of early stages of tectonic shortening, ii) the recognition of which mechanisms assist deformation during these stages and iii) the identification of parameters that actually steer the beginning of localisation. To clarify these points, the thesis presents the results of an integrated, multiscale and multi-technique structural study that relied on field and laboratory data to analyse the structural, architectural, mineralogical and geochemical features that govern deformation during compressional tectonics. By focusing on two case studies, the Eastern Southern Alps (northern Italy), where deformation is mainly brittle, and the Oman Mountains (northeastern Oman), where ductile deformation dominates, the thesis shows that the deformation localisation is steered by several mechanisms that mutually interact at different stages during compression. At shallow crustal conditions, derived conceptual and numerical models show that both inherited (e.g., stratigraphic) and acquired (e.g., structural) features play a key role in steering deformation and differentiating the seismic behaviour of the multilayer succession. At the same time, at deeper crustal conditions, strain localises in narrow domains in which fluids, temperature, shear strain and pressure act together during the development of the internal fabric and the chemical composition of mylonitic shear zones, in which localisation took place under high-pressure (HP) and low-temperature (LT) conditions. In particular, results indicate that those shear zones acted as “sheltering structural capsules” in which peculiar processes can happen and where the results of these processes can be successively preserved even over hundreds of millions of years.

A coupled ocean-atmosphere system for short term forecasting in the southern european seas

This coupled model combines two state-of-the-art numerical models, NEMO for the oceanic component and WRF for the atmospheric component and implements them at an appropriate resolution. The oceanic model has been implemented starting from the Mediterranean Forecasting System with a resolution of 1/24° and the domain was extended to exactly match the grid of a newly implemented atmospheric model for the same area. The uncoupled ocean model has been validated against SST observed data, both in the simulation of an extreme event and in the short-term forecast of two seasonal periods. A new setup of the model was successfully tested in which the downward radiative fluxes were prescribed from atmospheric forecasts. Various physical schemes, domain, boundary, and initial conditions were tested with the atmospheric model to obtain the best representation of medicane Ianos. The heat fluxes calculated by the uncoupled models were compared to determine which setup gave the best energy balance between the components of the coupled model. The coupling strategy used is the traditional one, where the ocean is driven by the surface stress, heat fluxes, and radiative fluxes computed in the atmospheric component, which in turn receives the SST and surface currents. As expected, the overall skills of the coupled model are slightly degraded compared to the uncoupled models, even though the positioning and timing of the cyclone at the time of the landfall is enhanced. The mean heat fluxes do not change compared to the uncoupled model, whereas the pattern of the shortwave radiation and latent heat is changed. Moreover, the two energy fluxes are larger in absolute values than those calculated with the MFS formulas. The fact that they have opposite signs give raise to a compensation error that limits the overall degradation of the coupled simulation.

A Trans-dimensional inversion algorithm to model deformation sources with unconstrained shape in finite element domains

Ground deformation provides valuable insights on subsurface processes with pattens reflecting the characteristics of the source at depth. In active volcanic sites displacements can be observed in unrest phases; therefore, a correct interpretation is essential to assess the hazard potential. Inverse modeling is employed to obtain quantitative estimates of parameters describing the source. However, despite the robustness of the available approaches, a realistic imaging of these reservoirs is still challenging. While analytical models return quick but simplistic results, assuming an isotropic and elastic crust, more sophisticated numerical models, accounting for the effects of topographic loads, crust inelasticity and structural discontinuities, require much higher computational effort and information about the crust rheology may be challenging to infer. All these approaches are based on a-priori source shape constraints, influencing the solution reliability. In this thesis, we present a new approach aimed at overcoming the aforementioned limitations, modeling sources free of a-priori shape constraints with the advantages of FEM simulations, but with a cost-efficient procedure. The source is represented as an assembly of elementary units, consisting in cubic elements of a regular FE mesh loaded with a unitary stress tensors. The surface response due to each of the six stress tensor components is computed and linearly combined to obtain the total displacement field. In this way, the source can assume potentially any shape. Our tests prove the equivalence of the deformation fields due to our assembly and that of corresponding cavities with uniform boundary pressure. Our ability to simulate pressurized cavities in a continuum domain permits to pre-compute surface responses, avoiding remeshing. A Bayesian trans-dimensional inversion algorithm implementing this strategy is developed. 3D Voronoi cells are used to sample the model domain, selecting the elementary units contributing to the source solution and those remaining inactive as part of the crust.

ICE Emissions Estimation through Machine-learning techniques to support real-world R&D experiments

The emissions estimation, both during homologation and standard driving, is one of the new challenges that automotive industries have to face. The new European and American regulation will allow a lower and lower quantity of Carbon Monoxide emission and will require that all the vehicles have to be able to monitor their own pollutants production. Since numerical models are too computationally expensive and approximated, new solutions based on Machine Learning are replacing standard techniques. In this project we considered a real V12 Internal Combustion Engine to propose a novel approach pushing Random Forests to generate meaningful prediction also in extreme cases (extrapolation, very high frequency peaks, noisy instrumentation etc.). The present work proposes also a data preprocessing pipeline for strongly unbalanced datasets and a reinterpretation of the regression problem as a classification problem in a logarithmic quantized domain. Results have been evaluated for two different models representing a pure interpolation scenario (more standard) and an extrapolation scenario, to test the out of bounds robustness of the model. The employed metrics take into account different aspects which can affect the homologation procedure, so the final analysis will focus on combining all the specific performances together to obtain the overall conclusions.

Analysis and numerical investigation of dynamic models for liquid chromatography

Magdeburg, Univ., Fak. für Mathematik, Diss., 2013

Orthotropic active strain models for the numerical simulation of cardiac biomechanics

An active strain formulation for orthotropic constitutive laws arising in cardiac mechanics modeling is introduced and studied. The passive mechanical properties of the tissue are described by the Holzapfel-Ogden relation. In the active strain formulation, the Euler-Lagrange equations for minimizing the total energy are written in terms of active and passive deformation factors, where the active part is assumed to depend, at the cell level, on the electrodynamics and on the specific orientation of the cardiac cells. The well-posedness of the linear system derived from a generic Newton iteration of the original problem is analyzed and different mechanical activation functions are considered. In addition, the active strain formulation is compared with the classical active stress formulation from both numerical and modeling perspectives. Taylor-Hood and MINI finite elements are employed to discretize the mechanical problem. The results of several numerical experiments show that the proposed formulation is mathematically consistent and is able to represent the main key features of the phenomenon, while allowing savings in computational costs.

Stochastic models and numerical algorithms for a class of regulatory gene networks.

Regulatory gene networks contain generic modules, like those involving feedback loops, which are essential for the regulation of many biological functions (Guido et al. in Nature 439:856-860, 2006). We consider a class of self-regulated genes which are the building blocks of many regulatory gene networks, and study the steady-state distribution of the associated Gillespie algorithm by providing efficient numerical algorithms. We also study a regulatory gene network of interest in gene therapy, using mean-field models with time delays. Convergence of the related time-nonhomogeneous Markov chain is established for a class of linear catalytic networks with feedback loops.