Climatology of Vb cyclones, physical mechanisms and their impact on extreme precipitation over Central Europe

Cyclones, which develop over the western Mediterranean and move northeastward are a major source of extreme weather and known to be responsible for heavy precipitation over the northern side of the Alpine range and Central Europe. As the relevant processes triggering these so-called Vb events and their impact on extreme precipitation are not yet fully understood, this study focuses on gaining insight into the dynamics of past events. For this, a cyclone detection and tracking tool is applied to the ERA-Interim reanalysis (1979–2013) to identify prominent Vb situations. Precipitation in the ERA-Interim and the E-OBS data sets is used to evaluate case-to-case precipitation amounts and to assess consistency between the two data sets. Both data sets exhibit high variability in precipitation amounts among different Vb events. While only 23 % of all Vb events are associated with extreme precipitation, around 15 % of all extreme precipitation days (99 percentile) over the northern Alpine region and Central Europe are induced by Vb events, although Vb cyclones are rare events (2.3 per year). To obtain a better understanding of the variability within Vb events, the analysis of the 10 heaviest and lowest precipitation Vb events reveals noticeable differences in the state of the atmosphere. These differences are most pronounced in the geopotential height and potential vorticity field, indicating a much stronger cyclone for heavy precipitation events. The related differences in wind direction are responsible for the moisture transport around the Alps and the orographical lifting along the northern slopes of the Alps. These effects are the main reasons for a disastrous outcome of Vb events, and consequently are absent in the Vb events associated with low precipitation. Hence, our results point out that heavy precipitation related to Vb events is mainly related to large-scale dynamics rather than to thermodynamic processes.

Physical modeling and numerical simulations of degradation mechanisms in devices and insulators for power applications

In this thesis, a TCAD approach for the investigation of charge transport in amorphous silicon dioxide is presented for the first time. The proposed approach is used to investigate high-voltage silicon oxide thick TEOS capacitors embedded in the back-end inter-level dielectric layers for galvanic insulation applications. In the first part of this thesis, a detailed review of the main physical and chemical properties of silicon dioxide and the main physical models for the description of charge transport in insulators are presented. In the second part, the characterization of high-voltage MIM structures at different high-field stress conditions up to the breakdown is presented. The main physical mechanisms responsible of the observed results are then discussed in details. The third part is dedicated to the implementation of a TCAD approach capable of describing charge transport in silicon dioxide layers in order to gain insight into the microscopic physical mechanisms responsible of the leakage current in MIM structures. In particular, I investigated and modeled the role of charge injection at contacts and charge build-up due to trapping and de-trapping mechanisms in the oxide layer to the purpose of understanding its behavior under DC and AC stress conditions. In addition, oxide breakdown due to impact-ionization of carriers has been taken into account in order to have a complete representation of the oxide behavior at very high fields. Numerical simulations have been compared against experiments to quantitatively validate the proposed approach. In the last part of the thesis, the proposed approach has been applied to simulate the breakdown in realistic structures under different stress conditions. The TCAD tool has been used to carry out a detailed analysis of the most relevant physical quantities, in order to gain a detailed understanding on the main mechanisms responsible for breakdown and guide design optimization.

Climatic patterns and physical controls of chlorophyll-a in the Northeast Atlantic

Tese de Doutoramento, Física, 17 de Dezembro de 2013, Universidade dos Açores.

Mechanisms of synchronization and pattern formation in a lattice of pulse-coupled oscillators

We analyze the physical mechanisms leading either to synchronization or to the formation of spatiotemporal patterns in a lattice model of pulse-coupled oscillators. In order to make the system tractable from a mathematical point of view we study a one-dimensional ring with unidirectional coupling. In such a situation, exact results concerning the stability of the fixed of the dynamic evolution of the lattice can be obtained. Furthermore, we show that this stability is the responsible for the different behaviors.

Electrification of precipitating systems over the Amazon: Physical processes of thunderstorm development

This study investigated the physical processes involved in the development of thunderstorms over southwestern Amazon by hypothesizing causalities for the observed cloud-to-ground lightning variability and the local environmental characteristics. Southwestern Amazon experiences every year a large variety of environmental factors, such as the gradual increase in atmospheric moisture, extremely high pollution due to biomass burning, and intense deforestation, which directly affects cloud development by differential surface energy partition. In the end of the dry period it was observed higher percentages of positive cloud-to-ground (+CG) lightning due to a relative increase in +CG dominated thunderstorms (positive thunderstorms). Positive (negative) thunderstorms initiated preferentially over deforested (forest) areas with higher (lower) cloud base heights, shallower (deeper) warm cloud depths, and higher (lower) convective potential available energy. These features characterized the positive (negative) thunderstorms as deeper (relatively shallower) clouds, stronger (relatively weaker) updrafts with enhanced (decreased) mixed and cold vertically integrated liquid. No significant difference between thunderstorms (negative and positive) and nonthunderstorms were observed in terms of atmospheric pollution, once the atmosphere was overwhelmed by pollution leading to an updraft-limited regime. However, in the wet season both negative and positive thunderstorms occurred during periods of relatively higher aerosol concentration and differentiated size distributions, suggesting an aerosol-limited regime where cloud electrification could be dependent on the aerosol concentration to suppress the warm and enhance the ice phase. The suggested causalities are consistent with the invoked hypotheses, but they are not observed facts; they are just hypotheses based on plausible physical mechanisms.

PHYSICAL FACTORS IN B CELL ACTIN DYNAMICS AND ACTIVATION

Cells adapt to their changing world by sensing environmental cues and responding appropriately. This is made possible by complex cascades of biochemical signals that originate at the cell membrane. In the last decade it has become apparent that the origin of these signals can also arise from physical cues in the environment. Our motivation is to investigate the role of physical factors in the cellular response of the B lymphocyte. B cells patrol the body for signs of invading pathogens in the form of antigen on the surface of antigen presenting cells. Binding of antigen with surface proteins initiates biochemical signaling essential to the immune response. Once contact is made, the B cell spreads on the surface of the antigen presenting cell in order to gather as much antigen as possible. The physical mechanisms that govern this process are unexplored. In this research, we examine the role of the physical parameters of antigen mobility and cell surface topography on B cell spreading and activation. Both physical parameters are biologically relevant as immunogens for vaccine design, which can provide laterally mobile and immobile antigens and topographical surfaces. Another physical parameter that influences B cell response and the formation of the cell-cell junction is surface topography. This is biologically relevant as antigen presenting cells have highly convoluted membranes, resulting in variable topography. We found that B cell activation required the formation of antigen-receptor clusters and their translocation within the attachment plane. We showed that cells which failed to achieve these mobile clusters due to prohibited ligand mobility were much less activation competent. To investigate the effect of topography, we use nano- and micro-patterned substrates, on which B cells were allowed to spread and become activated. We found that B cell spreading, actin dynamics, B cell receptor distribution and calcium signaling are dependent on the topographical patterning of the substrate. A quantitative understanding of cellular response to physical parameters is essential to uncover the fundamental mechanisms that drive B cell activation. The results of this research are highly applicable to the field of vaccine development and therapies for autoimmune diseases. Our studies of the physical aspects of lymphocyte activation will reveal the role these factors play in immunity, thus enabling their optimization for biological function and potentially enabling the production of more effective vaccines.

Secondary instabilities in the flow around two circular cylinders in tandem

Direct stability analysis and numerical simulations have been employed to identify and characterize secondary instabilities in the wake of the flow around two identical circular cylinders in tandem arrangements. The centre-to-centre separation was varied from 1.2 to 10 cylinder diameters. Four distinct regimes were identified and salient cases chosen to represent the different scenarios observed, and for each configuration detailed results are presented and compared to those obtained for a flow around an isolated cylinder. It was observed that the early stages of the wake transition changes significantly if the separation is smaller than the drag inversion spacing. The onset of the three-dimensional instabilities were calculated and the unstable modes are fully described. In addition, we assessed the nonlinear character of the bifurcations and physical mechanisms are proposed to explain the instabilities. The dependence of the critical Reynolds number on the centre-to-centre separation is also discussed.

Simulation of the load-unload response ratio and critical sensitivity in the lattice solid model

The Load-Unload Response Ratio (LURR) method is an intermediate-term earthquake prediction approach that has shown considerable promise. It involves calculating the ratio of a specified energy release measure during loading and unloading where loading and unloading periods are determined from the earth tide induced perturbations in the Coulomb Failure Stress on optimally oriented faults. In the lead-up to large earthquakes, high LURR values are frequently observed a few months or years prior to the event. These signals may have a similar origin to the observed accelerating seismic moment release (AMR) prior to many large earthquakes or may be due to critical sensitivity of the crust when a large earthquake is imminent. As a first step towards studying the underlying physical mechanism for the LURR observations, numerical studies are conducted using the particle based lattice solid model (LSM) to determine whether LURR observations can be reproduced. The model is initialized as a heterogeneous 2-D block made up of random-sized particles bonded by elastic-brittle links. The system is subjected to uniaxial compression from rigid driving plates on the upper and lower edges of the model. Experiments are conducted using both strain and stress control to load the plates. A sinusoidal stress perturbation is added to the gradual compressional loading to simulate loading and unloading cycles and LURR is calculated. The results reproduce signals similar to those observed in earthquake prediction practice with a high LURR value followed by a sudden drop prior to macroscopic failure of the sample. The results suggest that LURR provides a good predictor for catastrophic failure in elastic-brittle systems and motivate further research to study the underlying physical mechanisms and statistical properties of high LURR values. The results provide encouragement for earthquake prediction research and the use of advanced simulation models to probe the physics of earthquakes.

Desenvolupament d'un node sensor sense fils autoalimentat. Mòdul piezo-magnètic

En aquest projecte es pretén implementar un dispositiu capaç de ser auto-suficient i no dependre de cap tipus de pila, bateria o fil elèctric que l’abasteixi d’energia elèctrica. El dispositiu recol·lectarà la energia magnètica generada per la corrent elèctrica a un fil i la transformarà en energia elèctrica, que serà emmagatzemada per el seu posterior ús. A demès, aquest projecte s’ha desenvolupat en col·laboració amb un segon projecte, dintre del qual s’implementarà una xarxa de sensors, mitjançant el protocol MIWI. Aquest projecte es divideix en tres grans blocs. El primer bloc del projecte serà una introducció teòrica de tots els coneixements relacionats amb el concepte d’energy harvesting i els mecanismes físic implicats. Al segon bloc podrem veure com s’han realitzat els càlculs, simulacions i posada en marxa, dels diferents elements que formaran el dispositiu recol·lector d’energia. Per últim en el tercer bloc veurem el prototip ja implementat. Es valoraran els resultats obtinguts, i es veuran els temps que necessitarà per alimentar al microcontrolador.

Seismic attenuation and velocity dispersion in heterogeneous partially saturated porous rocks

Using a numerical approach, we explore wave-induced fluid flow effects in partially saturated porous rocks in which the gas-water saturation patterns are governed by mesoscopic heterogeneities associated with the dry frame properties. The link between the dry frame properties and the gas saturation is defined by the assumption of capillary pressure equilibrium, which in the presence of heterogeneity implies that neighbouring regions can exhibit different levels of saturation. To determine the equivalent attenuation and phase velocity of the synthetic rock samples considered in this study, we apply a numerical upscaling procedure, which permits to take into account mesoscopic heterogeneities associated with the dry frame properties as well as spatially continuous variations of the pore fluid properties. The multiscale nature of the fluid saturation is taken into account by locally computing the physical properties of an effective fluid, which are then used for the larger-scale simulations. We consider two sets of numerical experiments to analyse such effects in heterogeneous partially saturated porous media, where the saturation field is determined by variations in porosity and clay content, respectively. In both cases we also evaluate the seismic responses of corresponding binary, patchy-type saturation patterns. Our results indicate that significant attenuation and modest velocity dispersion effects take place in this kind of media for both binary patchy-type and spatially continuous gas saturation patterns and in particular in the presence of relatively small amounts of gas. The numerical experiments also show that the nature of the gas distribution patterns is a critical parameter controlling the seismic responses of these environments, since attenuation and velocity dispersion effects are much more significant and occur over a broader saturation range for binary patchy-type gas-water distributions. This analysis therefore suggests that the physical mechanisms governing partial saturation should be accounted for when analysing seismic data in a poroelastic framework. In this context, heterogeneities associated with the dry frame properties, which do not play important roles in wave-induced fluid flow processes per se, should be taken into account since they may determine the kind of gas distribution pattern taking place in the porous rock.

Seismic attenuation in heterogeneous partially saturated rocks

We study wave-induced fluid flow effects in porous rocks partially saturated with gas and water, where the saturation patterns are governed by mesoscopic heterogeneities associated with the dry frame properties. The link between the dry frame properties and the gas saturation is defined by the assumption of capillary pressure equilibrium, which in the presence of heterogeneity implies that neighboring regions can exhibit different levels of saturation. In order to determine the equivalent attenuation and phase velocity of the synthetic rock samples considered in this study, we apply a numerical upscaling procedure, which permits to take into account mesoscopic heterogeneities associated with the dry frame properties as well as spatially continuous variations of the pore fluid properties. We consider numerical experiments to analyze such effects in heterogeneous partially saturated porous media, where the saturation field is determined by realistic variations in porosity. Our results indicate that the spatially continuous nature of gas saturation inherent to this study is a critical parameter controlling the seismic response of these environments, which in turn suggests that the physical mechanisms governing partial saturation should be accounted for when analyzing seismic data in a poro-elastic context.

Splitting rivers at their seams: bifurcations and avulsion

River bifurcations are critical but poorly understood elements of many geomorphological systems. They are integral elements of alluvial fans, braided rivers, fluvial lowland plains, and deltas and control the partitioning of water and sediment through these systems. Bifurcations are commonly unstable but their lifespan varies greatly. In braided rivers bars and channels migrate, split and merge at annual or shorter timescales, thereby creating and abandoning bifurcations. This behaviour has been studied mainly by geomorphologists and fluid dynamicists. Bifurcations also exist during avulsion, the process of a river changing course on a floodplain or in a delta, which may take 102103 years and has been studied mainly by sedimentologists. This review synthesizes our current understanding of bifurcations and brings together insights from different research communities and different environmental settings. We consider the causes and initiation of bifurcations and avulsion, the physical mechanisms controlling bifurcation and avulsion evolution, mathematical and numerical modelling of these processes, and the possibility of stable bifurcations. We end the review with some open questions. Copyright (c) 2012 John Wiley & Sons, Ltd.

Fundamentals of heat transfer

Following over 170+ pages and additional appendixes are formed based on content of Course: Fundamentals of Heat Transfer. Mainly this summarizes relevant parts on Book of Fundamentals of Heat and Mass Transfer (Incropera), but also other references introducing the same concepts are included. Student’s point of view has been consideredwith following highlights: (1) Relevant topics are presented in a nutshell to provide fast digestion of principles of heat transfer. (2) Appendixes include terminology dictionary. (3) Totally 22 illustrating examples are connecting theory to practical applications and quantifying heat transfer to understandable forms as: temperatures, heat transfer rates, heat fluxes, resistances and etc. (4) Most important Learning outcomes are presented for each topic separately. The Book, Fundamentals of Heat and Mass Transfer (Incropera), is certainly recommended for those going beyond basic knowledge of heat transfer. Lecture Notes consists of four primary content-wise objectives: (1) Give understanding to physical mechanisms of heat transfer, (2)Present basic concepts and terminology relevant for conduction, convection and radiation (3) Introduce thermal performance analysis methods for steady state and transient conduction systems. (4) Provide fast-to-digest phenomenological understanding required for basic design of thermal models

Studies of Solar Activity with Emphasis on Quasi-Periodic Oscillations

In this thesis mainly long quasi-periodic solar oscillations in various solar atmospheric structures are discussed, based on data obtained at several wavelengths, focussing, however, mainly on radio frequencies. Sunspot (Articles II and III) and quiet Sun area (QSA) (Article I) oscillations are investigated along with quasi-periodic pulsations (QPP) in a flaring event with wide-range radio spectra (Article IV). Various oscillation periods are detected; 3–15, 35–70 and 90 minutes (QSA), 10-60 and 80-130 minutes (in sunspots at various radio frequencies), 3-5, 10-23, 220-240, 340 and 470 minutes (in sunspots at photosphere) and 8-12 and 15-17 seconds (in a solar flare at radio frequencies). Some of the oscillation periods are detected for the first time, while some of them have been confirmed earlier by other research groups. Solar oscillations can provide more information on the nature of various solar structures. This thesis presents the physical mechanisms of some solar structure oscillations. Two different theoretical approaches are chosen; magnetohydrodynamics (MHD) and the shallow sunspot model. These two theories can explain a wide range of solar oscillations from a few seconds up to some hours. Various wave modes in loop structures cause solar oscillations (<45 minutes) both in sunspots and quiet Sun areas. Periods lasting more than 45 minutes in the sunspots (and a fraction of the shorter periods) are related to sunspot oscillations as a whole. Sometimes similar oscillation periods are detected both in sunspot area variations and respectively in magnetic field strength changes. This result supports a concept that these oscillations are related to sunspot oscillations as a whole. In addition, a theory behind QPPs at radio frequencies in solar flares is presented. The thesis also covers solar instrumentation and data sources. Additionally, the data processing methods are presented. As the majority of the investigations in this thesis focus on radio frequencies, also the most typical radio emission mechanisms are presented. The main structures of the Sun, which are related to solar oscillations, are also presented. Two separate projects are included in this thesis. Solar cyclicity is studied using the extensively large solar radio map archieve from Metsähovi Radio Observatory (MRO) at 37 GHz, between 1978 and 2011 (Article V) covering two full solar cycles. Also, some new solar instrumentation (Article VI) was developed during this thesis.

Dynamo Magnétohydrodynamique en champ moyen

De nos jours, il est bien accepté que le cycle magnétique de 11 ans du Soleil est l'oeuvre d'une dynamo interne présente dans la zone convective. Bien qu'avec la puissance de calculs des ordinateurs actuels il soit possible, à l'aide de véritables simulations magnétohydrodynamiques, de résoudre le champ magnétique et la vitessse dans toutes les directions spatiales, il n'en reste pas moins que pour étudier l'évolution temporelle et spatiale de la dynamo solaire à grande échelle, il reste avantageux de travailler avec des modèles plus simples. Ainsi, nous avons utilisé un modèle simplifié de la dynamo solaire, nommé modèle de champ moyen, pour mieux comprendre les mécanismes importants à l'origine et au maintien de la dynamo solaire. L'insertion d'un tenseur-alpha complet dans un modèle dynamo de champ moyen, provenant d'un modèle global-MHD [Ghizaru et al., 2010] de la convection solaire, nous a permis d'approfondir le rôle que peut jouer la force électromotrice dans les cycles magnétiques produits par ce modèle global. De cette façon, nous avons pu reproduire certaines caractéristiques observées dans les cycles magnétiques provenant de la simulation de Ghizaru et al., 2010. Tout d'abord, le champ magnétique produit par le modèle de champ moyen présente deux modes dynamo distincts. Ces modes, de périodes similaires, sont présents et localisés sensiblement aux mêmes rayons et latitudes que ceux produits par le modèle global. Le fait que l'on puisse reproduire ces deux modes dynamo est dû à la complexité spatiale du tenseur-alpha. Par contre, le rapport entre les périodes des deux modes présents dans le modèle de champ moyen diffère significativement de celui trouvé dans le modèle global. Par ailleurs, on perd l'accumulation d'un fort champ magnétique sous la zone convective dans un modèle où la rotation différentielle n'est plus présente. Ceci suggère que la présence de rotation différentielle joue un rôle non négligeable dans l'accumulation du champ magnétique à cet endroit. Par ailleurs, le champ magnétique produit dans un modèle de champ moyen incluant un tenseur-alpha sans pompage turbulent global est très différent de celui produit par le tenseur original. Le pompage turbulent joue donc un rôle fondamental au sein de la distribution spatiale du champ magnétique. Il est important de souligner que les modèles dépourvus d'une rotation différentielle, utilisant le tenseur-alpha original ou n'utilisant pas de pompage turbulent, parviennent tous deux à produire une dynamo oscillatoire. Produire une telle dynamo à l'aide d'un modèle de ce type n'est pas évident, a priori. Finalement, l'intensité ainsi que le type de profil de circulation méridienne utilisés sont des facteurs affectant significativement la distribution spatiale de la dynamo produite.