numerical problems in the calculation of electric field and charge density profiles for mass impregnated non-draining hvdc cables subjected to fast polarity reversals

The goal of this simulation thesis is to present a tool for studying and eliminating various numerical problems observed while analyzing the behavior of the MIND cable during fast voltage polarity reversal. The tool is built on the MATLAB environment, where several simulations were run to achieve oscillation-free results. This thesis will add to earlier research on HVDC cables subjected to polarity reversals. Initially, the code does numerical simulations to analyze the electric field and charge density behavior of a MIND cable for certain scenarios such as before, during, and after polarity reversal. However, the primary goal is to reduce numerical oscillations from the charge density profile. The generated code is notable for its usage of the Arithmetic Mean Approach and the Non-Uniform Field Approach for filtering and minimizing oscillations even under time and temperature variations.

Numerical simulation of interdigitated back contact hetero-junction solar cells

The goal of this thesis is the application of an opto-electronic numerical simulation to heterojunction silicon solar cells featuring an all back contact architecture (Interdigitated Back Contact Hetero-Junction IBC-HJ). The studied structure exhibits both metal contacts, emitter and base, at the back surface of the cell with the objective to reduce the optical losses due to the shadowing by front contact of conventional photovoltaic devices. Overall, IBC-HJ are promising low-cost alternatives to monocrystalline wafer-based solar cells featuring front and back contact schemes, in fact, for IBC-HJ the high concentration doping diffusions are replaced by low-temperature deposition processes of thin amorphous silicon layers. Furthermore, another advantage of IBC solar cells with reference to conventional architectures is the possibility to enable a low-cost assembling of photovoltaic modules, being all contacts on the same side. A preliminary extensive literature survey has been helpful to highlight the specific critical aspects of IBC-HJ solar cells as well as the state-of-the-art of their modeling, processing and performance of practical devices. In order to perform the analysis of IBC-HJ devices, a two-dimensional (2-D) numerical simulation flow has been set up. A commercial device simulator based on finite-difference method to solve numerically the whole set of equations governing the electrical transport in semiconductor materials (Sentuarus Device by Synopsys) has been adopted. The first activity carried out during this work has been the definition of a 2-D geometry corresponding to the simulation domain and the specification of the electrical and optical properties of materials. In order to calculate the main figures of merit of the investigated solar cells, the spatially resolved photon absorption rate map has been calculated by means of an optical simulator. Optical simulations have been performed by using two different methods depending upon the geometrical features of the front interface of the solar cell: the transfer matrix method (TMM) and the raytracing (RT). The first method allows to model light prop-agation by plane waves within one-dimensional spatial domains under the assumption of devices exhibiting stacks of parallel layers with planar interfaces. In addition, TMM is suitable for the simulation of thin multi-layer anti reflection coating layers for the reduction of the amount of reflected light at the front interface. Raytracing is required for three-dimensional optical simulations of upright pyramidal textured surfaces which are widely adopted to significantly reduce the reflection at the front surface. The optical generation profiles are interpolated onto the electrical grid adopted by the device simulator which solves the carriers transport equations coupled with Poisson and continuity equations in a self-consistent way. The main figures of merit are calculated by means of a postprocessing of the output data from device simulation. After the validation of the simulation methodology by means of comparison of the simulation result with literature data, the ultimate efficiency of the IBC-HJ architecture has been calculated. By accounting for all optical losses, IBC-HJ solar cells result in a theoretical maximum efficiency above 23.5% (without texturing at front interface) higher than that of both standard homojunction crystalline silicon (Homogeneous Emitter HE) and front contact heterojuction (Heterojunction with Intrinsic Thin layer HIT) solar cells. However it is clear that the criticalities of this structure are mainly due to the defects density and to the poor carriers transport mobility in the amorphous silicon layers. Lastly, the influence of the most critical geometrical and physical parameters on the main figures of merit have been investigated by applying the numerical simulation tool set-up during the first part of the present thesis. Simulations have highlighted that carrier mobility and defects level in amorphous silicon may lead to a potentially significant reduction of the conversion efficiency.

Shear behaviour of reinforced cconcrete slab under concentrated load: an investigation through non-linear and sequentially linear analysis

English: The assessment of safety in existing bridges and viaducts led the Ministry of Public Works of the Netherlands to finance a specific campaing aimed at the study of the response of the elements of these infrastructures. Therefore, this activity is focused on the investigation of the behaviour of reinforced concrete slabs under concentrated loads, adopting finite element modeling and comparison with experimental results. These elements are characterized by shear behaviour and crisi, whose modeling is, from a computational point of view, a hard challeng, due to the brittle behavior combined with three-dimensional effects. The numerical modeling of the failure is studied through Sequentially Linear Analysis (SLA), an alternative Finite Element method, with respect to traditional incremental and iterative approaches. The comparison between the two different numerical techniques represents one of the first works and comparisons in a three-dimensional environment. It's carried out adopting one of the experimental test executed on reinforced concrete slabs as well. The advantage of the SLA is to avoid the well known problems of convergence of typical non-linear analysis, by directly specifying a damage increment, in terms of reduction of stiffness and resistance in particular finite element, instead of load or displacement increasing on the whole structure . For the first time, particular attention has been paid to specific aspects of the slabs, like an accurate constraints modeling and sensitivity of the solution with respect to the mesh density. This detailed analysis with respect to the main parameters proofed a strong influence of the tensile fracture energy, mesh density and chosen model on the solution in terms of force-displacement diagram, distribution of the crack patterns and shear failure mode. The SLA showed a great potential, but it requires a further developments for what regards two aspects of modeling: load conditions (constant and proportional loads) and softening behaviour of brittle materials (like concrete) in the three-dimensional field, in order to widen its horizons in these new contexts of study.

Numerical simulations of the effects of residual stresses in adhesively bonded composite specimen

In this work the problem of performing a numerical simulation of quasi-static crack propagation within an adhesive layer of a bonded joint under Mode I loading affected by stress field changes due to thermal-chemical shrinkage induced by cure process is addressed. Secondly, a parametric study on fracture critical energy, cohesive strength and Young's modulus is performed. Finally, a particular case of adhesive layer stiffening is simulated in order to verify qualitatively the major effect.

Analytical and numerical study of polarons

The aim of this thesis is to introduce the polaron concept and to perform a DFT numerical calculation of a small polaron in the rutile phase of TiO2. In the first chapters, we present an analytical study of small and large polarons, based on the Holstein and Fröhlich Hamiltonians. The necessary mathematical formalism and physics fundamentals are briefly reviewed in the first chapter. In the second part of the thesis, Density Functional Theory (DFT) is introduced together with the DFT+U correction and its implementation in the Vienna Ab-Initio Simulation Package (VASP). The calculation of a small polaron in rutile is then described and discussed at a qualitative level. The polaronic solution is compared with the one of a delocalized electron. The calculation showed how the polaron creates a new energy level 0.70 eV below the conduction band. The energy level is visible both in the band structure diagram and in the density of states diagram. The electron is localized on a titanium atom, distorting the surrounding lattice. In particular, the four oxygen atoms closer to the titanium atom are displaced by 0.085 Å outwards, whereas the two further oxygen atoms by 0.023 Å. The results are compatible, at a qualitative level, with the literature. Further developments of this work may try to improve the precision of the results and to quantitatively compare them with the literature.

Determination of local density of states in amorphous oxide semiconductors with Kelvin Probe Force Microscopy

Amorphous semiconductors are important materials as they can be deposited by physical deposition techniques on large areas and even on plastic substrates. Therefore, they are crucial for transistors in large active matrices for imaging and transparent wearable electronics. The most widely applied candidate for amorphous thin film transistors production is Indium Gallium Zinc Oxide (IGZO). It is attracting much interest because of its optical transparency, facile processing by sputtering deposition and notable improved charge carrier mobility with respect to hydrogenated amorphous silicon a-Si:H. Degradation of the device and long-term performance issues have been observed if IGZO thin film transistors are subjected to electrical stress, leading to a modification of IGZO channel properties and subthreshold slope. Therefore, it is of great interest to have a reliable and precise method to study the conduction band tail, and the density of states in amorphous semiconductors. The aim of this thesis is to develop a local technique using Kelvin Probe Force Microscopy to study the evolution of IGZO DOS properties. The work is divided into three main parts. First, solutions to the non-linear Poisson-Boltzmann equation of a metal-insulator-semiconductor junction describing the charge accumulation and its relation to DOS properties are elaborated. Second macroscopic techniques such as capacitance voltage (CV) measurements and photocurrent spectroscopy are applied to obtain a non-local estimate of band-tail DOS properties in thin film transistor samples. The third part of my my thesis is dedicated to the KPFM measurements. By fitting the data to the developed numerical model, important parameters describing the amorphous conduction band tail are obtained. The results are in excellent agreement with the macroscopic characterizations. KPFM result is comparable also with non-local optoelectronic characterizations, such as photocurrent spectroscopy.

Debonding and crack propagation in cracked steel beams strengthened with CFRP and subjected to cyclic loads.

The reinforcement methods used to restore or increase the bearing capacity of metal structures are based on the application of steel plates to be bolted or welded to the original structure, which can cause problems to the integrity of the original structure. These difficulties can be overcome with the introduction of fiber-reinforced composite materials. FRPs are characterized by high strength to weight ratio, and they are very resistant to corrosion. In this dissertation a cracked steel I-beam reinforced with Carbon Fiber-Reinforced Polymer will be studied by performing a numerical evaluation of the structure with the commercial Finite Element Method software ABAQUS. The crack propagation will be computed using XFEM, while the debonding of the reinforcement layer will be found by considering a cohesive contact interface between the beam and the CFRP plate. The results will show the efficiency of the strengthening method in increasing the load carrying capacity of the cracked beam, and in reducing the crack opening of the initial notch.

Numerical study on viscoelastic behavior of polypropylene fiber reinforced concrete subjected to temperature variations using LDPM theory

This thesis is focused on the viscoelastic behavior of macro-synthetic fiber-reinforced concrete (MSFRC) with polypropylene studied numerically when subjected to temperature variations (-30 oC to +60 oC). LDPM (lattice discrete particle model), a meso-scale model for heterogeneous composites, is used. To reproduce the MSFRC structural behavior, an extended version of LDPM that includes fiber effects through fiber-concrete interface micromechanics, called LDPM-F, is applied. Model calibration is performed based on three-point bending, cube, and cylinder test for plain concrete and MSFRC. This is followed by a comprehensive literature study on the variation of mechanical properties with temperature for individual fibers and plain concrete. This literature study and past experimental test results constitute inputs for final numerical simulations. The numerical response of MSFRC three-point bending test is replicated and compared with the previously conducted experimental test results; finally, the conclusions were drawn. LDPM numerical model is successfully calibrated using experimental responses on plain concrete. Fiber-concrete interface micro-mechanical parameters are subsequently fixed and LDPM-F models are calibrated based on MSFRC three-point bending test at room temperature. Number of fibers contributing crack bridging mechanism is computed and found to be in good agreement with experimental counts. Temperature variations model for individual constituents of MSFRC, fibers and plain concrete, are implemented in LDPM-F. The model is validated for MSFRC three-point bending stress-CMOD (crack mouth opening) response reproduced at -30 oC, -15 oC, 0 oC, +20 oC, +40 oC and +60 oC. It is found that the model can well describe the temperature variation behavior of MSFRC. At positive temperatures, simulated responses are in good agreement. Slight disagreement in negative regimes suggests an in-depth study on fiber-matrix interface bond behavior with varying temperatures.

Experimental and numerical characterization of flame instability in hydrogen-containing flames

In the framework of the energy transition, the acquisition of proper knowledge of fundamental aspects characterizing the use of alternative fuels is paramount as well as the development of optimized know-how and technologies. In this sense, the use of hydrogen has been indicated as a promising route for decarbonization at the end-users stage in the energy supply chain. However, the elevated reactivity and the low-density at atmospheric conditions of hydrogen pose new challenges. Among the others, the dilution of hydrogen with carbon dioxide from carbon capture and storage systems represents a possible route. However, the interactions between these species have been poorly studied so far. For these reasons, this thesis, in collaboration between the University of Bologna and Technische Universität Bergakademie of Freiberg in Saxony (Germany), investigates the laminar flame of hydrogen-based premixed gas with the dilution of carbon dioxide. An experimental system, called a heat flux burner, was adopted ad different operating conditions. The presence of the cellularity phenomenon, forming the so-called cellular flame, was observed and analysed. Theoretical and visual methods have allowed for the characterization of the investigated flames, opening new alternatives for sustainable energy production via hydrogen transformation.