The present thesis deals with soil resistivity and especially the suitability of replacing various soil models with a two-layer horizontal model. First, the factors that influence soil resistivity as well as the methodology of its measuring are described, and various soil models are introduced. Then, the Ansys program was used to create several soil models: a vertical model with perpendicular or oblique division, a model with a pronounced inhomogeneity and a three-layer model. In all of these models, a Wenner measurement simulation was performed, followed by a two-layer horizontal replacement model. Then the size of the error that was created by replacing the actual model with a two-layer horizontal model was examined.

<p>U radu je pokazano da frekvencijski zavisna impedansa<br />zemlje kao povratnog provodnika, pri proizvoljnoj<br />učestanosti može veoma tačno da se odredi pomoću<br />fluksa kompleksnog Pointingovog vektora. Za<br />izračunavanje kompleksnog Pointingovog vektora,<br />neophodno je prethodno odrediti raspodelu struje u<br />zemlji, kao i raspodelu magnetskog polja i u zemlji i u<br />vazduhu iznad povr&scaron;i zemlje. Obe ove raspodele takođe<br />su značajne za analizu elektroenergetskih sistema.</p> / <p>The thesis shows that the Earth return impedance at<br />arbitrary low frequency can be accurately determined<br />from the complex Poynting vector flux. For the complex<br />Poynting vector calculation, first it is necessary to<br />determine the current distribution inside the ground, as<br />well as the magnetic field distribution both inside the<br />ground and in the air above the ground surface. Both<br />distributions are also significant for power electrical<br />systems analysis.</p>

The objective of this preliminary study was to determine the potential of diffuse correlation spectroscopy (DCS) for assessment of coagulation depth. Coagulation of tissue can occur due to a number of different reasons such as thermal or electrical burns or radiofrequency ablation. DCS is a non-invasive optical technique which can be used to determine the optical and dynamic properties of tissue by fitting a theoretical model of photon propagation in multiply scattering tissue to experimental data obtained from measurements. The DCS measurements were performed on two-layered phantom models that represent healthy tissue with high flow properties with a layer of coagulated tissue with low flow properties on top. Three different phantom models were prepared using gelatin-Intralipid gels, PDMS, and nylon as an upper layer, and an Intralipid solution was used for the bottom layer for all three phantoms. DCS measurements were performed on all three phantom models with varying thicknesses of the upper layers, and varying source-detector separations. The acquired data from the DCS measurement were analyzed in MATLAB in order to obtain the electric field temporal autocorrelation function. A theoretical model describing photon propagation in a two-layered medium was fitted to the obtained data in order to extract the desired parameters. The results showed that the thickness of the gelatin-Intralipid gels could be extracted within a 0.5 mm certainty and the thickness of the PDMS phantoms could also be extracted within approximately 0.7 mm. For the nylon phantoms, the results obtained were not good because the fitting was not successful and the thickness was not extracted appropriately. There is potential in DCS for assessment of burn wound depth but further research and development has to be done in the field in order to obtain more accurate and reliable results.

<div>Thermal management in next generation power electronic devices, radar applications and semiconductor packaging architectures is becoming increasingly challenging due to the need to reject localized high heat fluxes as well as large total powers. Air cooling has been considered as a simple and reliable method for thermal management compared to architectures that incorporate liquid cooling. However, air-cooled heat sinks typically require effective heat spreading to provide the requisite level of area enhancement to dissipate high heat fluxes. Compared to solid metallic heat spreaders, advanced heat sinks that incorporate two-phase heat transfer devices such as vapor chambers can significantly enhance the power dissipation capabilities in such configurations. Vapor chambers are devices that utilize evaporation/boiling processes within a sealed cavity to achieve efficient heat spreading. In high-heat-flux applications, boiling can occur within the internal wick structure of the vapor chamber at the location of the heat input (i.e., the evaporator). The maximum dryout heat flux and thermal resistance of the device is dictated by the resulting two-phase flow and heat transfer in the porous evaporator due to boiling. While various works in the literature have introduced new evaporator wick designs to improve the dryout heat flux during boiling, the enhancement is limited to small, millimeter scale hotspots or at a very high thermal resistance. In additixon, the effective design of such evaporator systems requires mechanistic models that can accurately predict the dryout limit and thermal performance. </div><div> This thesis first explores the usage of a novel ‘two-layer’ evaporator wick for passive high heat flux dissipation over large heater areas at a low thermal resistance. Moreover, a new mechanistic (first principles based) model framework is introduced for dryout limit and thermal performance prediction during boiling in capillary fed evaporators, by considering the resulting simultaneous flow of two phases (liquid and vapor) within the microscale porous media.</div><div> The novel two-layer wick concept uses a thick ‘cap’ layer of porous material to feed liquid to a thin ‘base’ layer through an array of vertical liquid-feeding ‘posts’. Vapor ‘vents’ in the cap layer allow for vapor formed during the boiling process (which is constrained to the base layer) to escape out of the wick. This two-layer structure decouples the functions of liquid resupply and capillary-fed boiling heat transfer, making the design realize high heat flux dissipation greater than 500 W/cm2 over large heat input areas of ~1 cm2. A reduced-order model is first developed to demonstrate the performance of a vapor chamber incorporating such a two-layer evaporator wick design. The model comprises simplified hydraulic and thermal resistance networks for predicting the capillary-limited maximum heat flux and the overall thermal resistance, respectively. The reduced-order model is validated against a higher fidelity numerical model and then used to analyze the performance of the vapor chamber with varying two-layer wick geometric feature sizes. The fabrication of the proposed two-layer wick is then presented. The thermal performance of the fabricated wicks is characterized using a boiling test facility that utilizes high speed visualization to identify the characteristic regimes of boiling operation in the wicks. The performance is also benchmarked to conventional single-layer wicks. </div><div> It is observed that single-layer wicks exhibit an unfavorable boiling regime where the center of the heater area dries out locally, leading to a high value of thermal resistance. The two-layer wicks avoid local dryout due to the distributed feeding provided by the posts and enhance the dryout heat flux significantly compared to single-layer wicks. A two-layer design that consists of a 10 × 10 array of liquid feeding posts provided a 400% improvement in the dryout heat flux. Following a parametric analysis of the effect of particle size, two-layer wicks composed of 180 – 212 µm particles and a 15 × 15 array of liquid feeding posts yielded a maximum heat flux dissipation of 485 W/cm2 over a 1 cm2 heat input area while also maintaining a low thermal resistance of only ~0.052 K/W. The effect of vapor venting and liquid-feeding areas is also experimentally studied. By understanding these effects, a parametrically optimized design is fabricated and shown to demonstrate an extremely high dryout limit of 512 W/cm2. We identify that the unique area-scalability of the two-layer wick design allows it to achieve an unprecedented combination of high total power and low-thermal-resistance heat dissipation over larger areas than was previously possible in the literature.</div><div> The results from the characterization of two-layer wicks revealed that the overall performance of the design was limited by the boiling process in the thin base wick layer. A fundamental model-based understanding of the resulting two-phase flow and heat transfer process in such thin capillary-fed porous media was still lacking. This lack of a mechanistic model precluded the accurate prediction of dryout heat flux and thermal performance of the two-layer wick. Moreover, such an understanding is needed for the optimal design of advanced hybrid evaporator wicks that leverage capillary-fed boiling. Despite the existence of various experimental works, there are currently no mechanistic approaches that model this behavior. To fill this unmet need, this thesis presents a new semi-empirical model for prediction of dryout and thermal resistance of capillary-fed evaporator systems. Thermal conduction across the solid and volumetric evaporation within the pores are solved to obtain the temperature distribution in the porous structure. Capillary-driven lateral liquid flow from the outer periphery of the evaporator to its center, with vapor flow across the thickness, is considered to obtain the local liquid and vapor pressures. Experiments are conducted on sintered copper particle evaporators of different particle sizes and heater areas to collect data for model calibration. To demonstrate the wider applicability of the model for other types of porous evaporators, the model is further calibrated against a variety of dryout limit and thermal resistance data collected from the literature. The model is shown to predict the experimentally observed trends in the dryout limit with mean particle/pore size, heater size, and evaporator thicknesses. This physics–based modeling approach is then implemented into a vapor chamber model to predict the thermal performance limits of air-cooled heat sinks with embedded vapor chambers. The governing energy and momentum equations of a low-cost analytical vapor chamber modeling approach is coupled with the evaporator model to capture the effect of boiling in the evaporator wick. An example case study illustrating the usage of the model is demonstrated and compared to a purely evaporation-based modeling approach, for quantifying the differences in dryout limit prediction, signifying the need to account for boiling in the evaporator wick. </div><div> The understanding gained from this thesis can be utilized for the prediction of dryout and thermal performance during boiling in capillary limited evaporator systems. The work also suggests the usage of a universal relative permeability correlation for the two-phase flow configuration studied herein for capillary-fed boiling, based on a wide calibration to experimental data. The modeling framework can also be readily leveraged to find novel and unexplored designs of advanced evaporator wicks. From an application standpoint, the new vapor chamber model developed here can be used for the improved estimation of performance limits specifically when high heat fluxes are encountered by the device. This will enable better and informed design of air-cooled heat sink architectures with embedded vapor chambers for high performance applications. </div><div><br></div>

Ce travail est consacré à la modélisation mathématique et numérique des écoulements eau-air en conduite qui interviennent notamment dans les centrales de production d’électricité ou les réseaux d’eaux usées. On s’intéresse particulièrement aux écoulements mixtes caractérisés par la présence de régimes stratifiés pilotés par des ondes gravitaires lentes, de régimes en charge ou secs (conduite remplie d’eau ou d’air) pilotés par des ondes acoustiques rapides, et de poches d’air piégées. Une modélisation précise de ces écoulements est nécessaire afin de garantir le bon fonctionnement du circuit hydraulique sous-jacent. Alors que la plupart des modèles disponibles dans la littérature se concentrent sur la phase eau en négligeant la présence de l’air, un modèle bicouche compressible prenant en compte les interactions eau-air est proposé dans cette thèse. Sa construction réside dans l’intégration des équations d’Euler barotropes sur la hauteur de chaque phase et dans l’application de la contrainte hydrostatique sur le gradient de pression de l’eau. Le modèle obtenu est hyperbolique et satisfait une inégalité d’entropie en plus d’autres propriétés mathématiques notables, telles que l’unicité des relations de saut ou la positivité des hauteurs et densités de chaque phase. Au niveau discret, la simulation d’écoulements mixtes avec le modèle bicouche compressible soulève plusieurs défis en raison de la disparité des vitesses d’ondes caractérisant chaque régime, des processus de relaxation rapide sous-jacents, et de la disparition de l’une des phases dans les régimes en charge ou sec. Une méthode à pas fractionnaires implicite-explicite est alors développée en s’appuyant sur la relaxation rapide en pression et sur le mimétisme avec les équations de Saint-Venant pour la dynamique lente de la phase eau. En particulier, une approche par relaxation permet d’obtenir une stabilisation du schéma en fonction du régime d’écoulement. Plusieurs cas tests sont traités et démontrent la capacité du modèle proposé à gérer des écoulements mixtes incluant la présence de poches d’air piégées. / The present work is dedicated to the mathematical and numerical modelling of transient air-water flows in pipes which occur in piping systems of several industrial areas such as nuclear or hydroelectric power plants or sewage pipelines. It deals more specifically with the so-called mixed flows which involve stratified regimes driven by slow gravity waves, pressurized or dry regimes (pipe full of water or air) driven by fast acoustic waves and entrapped air pockets. An accurate modelling of these flows is necessary to guarantee the operability of the related hydraulic system. While most of available models in the literature focus on the water phase neglecting the air phase, a compressible two-layer model which accounts for air-water interactions is proposed herein. The derivation process relies on a depth averaging of the isentropic Euler set of equations for both phases where the hydrostatic constraint is applied on the water pressure gradient. The resulting system is hyperbolic and satisfies an entropy inequality in addition to other significant mathematical properties, including the uniqueness of jump conditions and the positivity of heights and densities for each layer. Regarding the discrete level, the simulation of mixed flows with the compressible two-layer model raises key challenges due to the discrepancy of wave speeds characterizing each regime combined with the fast underlying relaxation processes and with phase vanishing when the flow becomes pressurized or dry. Thus, an implicit-explicit fractional step method is derived. It relies on the fast pressure relaxation in addition to a mimetic approach with the shallow water equations for the slow dynamics of the water phase. In particular, a relaxation method provides stabilization terms activated according to the flow regime. Several test cases are performed and attest the ability of the compressible two-layer model to deal with mixed flows in pipes involving air pocket entrapment.

Denomina-se Sistema Corrente do Brasil (CB) ao sistema de correntes de contorno oeste formados pela CB, fluindo para sul-sudoeste, e Corrente de Contorno Intermediária (CCI), fluindo para norte-nordeste, associado aos meandros e vórtices. O sistema bordeja a margem continental brasileira sudeste ao longo de 1200-1800 m. A investigação da dinâmica das ondas baroclínicas de vorticidade superpostas ao Sistema Corrente do Brasil é conduzida através exclusivamente de análise de dados hidrográficos de dois cruzeiros de meso-escala (verão e inverno de 1993) do projeto Circulação Oceânica da Região Oeste do Atlântico Sul - COROAS, a componente brasileira do World Ocean Circulation Experiment - WOCE. Estes cruzeiros foram realizados na porção central do Embaiamento de São Paulo (\'23,5 GRAUS\' - \'27 GRAUS\'S). Objetivando-se a aplicação do Método Dinâmico, estimativa quantitativa de um Nível de Referência (NR) é conduzida comparando-se duas metodologias diferentes e independentes. A primeira utiliza o fato de que a CB transporta Água Tropical (AT) e Água Central do Atlântico Sul (ACAS) e de que a CCI transporta dominantemente Água Intermediária Antártica (AIA). Aplica-se então um Critério Termodinâmico, no qual a interface que separa ACAS e AIA é estimada para os dois cruzeiros via aplicação dos Teoremas de Shtokman. A outra metodologia consiste na obtenção de padrões verticais de velocidade baroclínica absoluta a partir de simulações numéricas com a versão seccional do Princenton Ocean Model - POM inicializados com 14 campos termohalinos interpolados a partir dos dados - é o Critério Dinâmico. Os dois critérios apresentam resultado médio que difere apenas cerca de 7 dbar um do outro. Assim, adota-se valor de 480 dbar como NR para os cálculos geostróficos como nível isobárico que demarca a interface média entre CB e CCI. Mapas de função de corrente geostrófica são gerados através de mapeamento objetivo. Tais mapas apresentam evidências inequívocas de que os ciclones e anticiclones, observados nos experimentos de verão e inverno de 1993, são estruturas de vórtices associadas a meandros do Sistema CB. Este aparecem como estruturas vorticais postadas em lados diametralmente opostos ao eixo da CB: ciclones no lado costeiro, e anticiclones no lado oceânico. A variação vertical das estruturas dos meandros e a existência de uma única inversão de sentido permite comprovar e concluir que são ondas baroclínicas de primeiro modo, corroborando especulações presentes na literatura. A construção de um modelo quase-geostrófico de duas camadas para um oceano não-viscoso no plano \'beta\' e de fundo plano permite a realização de análise de vorticidade potencial. Tal análise permite concluir que o campo de vorticidade básico devido à CB suplanta o planetário e que estas ondas baroclínicas são ondas que devem sua existência ao cisalhamento vertical e horizontal da corrente. A variação frontal em vorticidade potencial baroclínica do Sistema CB é de \'1,7x10 POT.-5 s POT. -1\'. A dominãncia do termo da vorticidade de estiramento, que responde por 60% da variação de vorticidade na frente, permite estabelecer que as ondas capturadas na malha hidrográfica são ondas longas dentro da classe de meso-escala. O princípio de conservação de vorticidade potencial é invocado numa análise que envolve a superposição dos campos de vorticidade potencial e função de corrente para a primeira camada. Esta análise apresenta tanto evidências robustas de propagação das ondas baroclínicas quanto indícios de crescimento do meandro ciclônico, ou seja, de processo de instabilidade. Face aos achados das análises anteriores, teoria linear e método das perturbações são utilizados para obter uma relação de dispersão para as ondas baroclínicas da CB. Esta relação é dependente da velocidade da CB e do gradiente de vorticidade potencial através da corrente / The Brazil Current System is formed by two western boundary currents that flow along the Southeast Brazil continental margin from surface to 1200-1800 m deep. These two currents are the southward-flowing Brazil Current (BC) and the northward-flowing Intermediate Western Boundary Current (IWBC). The investigation of the dynamics of the barociinic vorticity waves superimposed to the BC system is conducted in this work solely from hydrographic data anaiysis and manipulation from two oceanographic meso-scale surveys that were part of the COROAS Project, the Brazilian arm of the Worid Ocean Circulation Experiment. These two cruises sampled the central portion of the São Paulo Bight (23,5°- 27°S) in the summer and winter seasons of 1993. As it was intended to appiy the classical Dynamic Method to the data, a quantitative estimate of the reference levei (RL) is conducted through the comparison of two different and independent methodoiogies. The first uses the previous knowiedge that the BC transports Tropical Water and South Atlantic Central Water (SACW) as well as that the IWBC transports mainly Antarctic Intermediate Water (AAIW). This method, referred here as the Thermodynamic Criterion applies the Shtokman theorems to estimate the interface depth between SACW and AAIW. The second methodology, designated here as the Dynamic Criterion, consists of modeiing absolute baroclinic velocities for all 14 hydrographic transects using the sectional version of the Princeton Ocean Model. The results of the two methods differ oniy in 7 dbar and a RL of 480 dbar is adopted as to represent the average interface isobaric level between BC and IWBC in the following geostrophic calculations. Horizontal distributions of geostrophic stream function are generated using objective mapping. These distributions present unequivocal evidences of cyclonic and anticyclonic structures in both Summer and Winter 1993 cruises. These features are associated to meanders of the BC System. They are depicted in opposing sides of the current core. The cyclone is seen in the BC coastal side as the anticyclone is placed in its oceanic side. The vertical variation of these pattems with a single flow direction inversion confirms that these vortical features are part of a first baroclinic mode wave, which corroborates previous speculations found in the literature. In order to pursue a potential vorticity analysis, a data-derived two-layer quasi-geostrophic model is built assuming an invicid flat-bottomed ocean in the 3-plane. This analysis allowed to conclude that the basic vorticity field associated with both horizontal and vertical shear of the boundary currents are responsible for the baroclinic wave existence. The planetary vorticity gradient is one order of magnitude lower. The frontal variation in baroclinic potential vorticity is 1,7 x i0 s1. The dominance of the stretching vorticity, which accounts for 60% of the gradient variation, leads to classify these oscillatory motions as long meso-scale waves. The potential vorticity conservation principle is invoked on an analysis that consists of superimposing the first layer quasi-gesotrophic potential vorticity and geostrophic stream funtion maps. This analysis revealed that the baroclinic waves are propagating as well as evidences of meander growth, an indication of a possible geophysical instability mechanism, are seen. Given the findings of the previous analyses, linear theory and the perturbation method are used to derive a dispersion relation for the BC System first mode baroclinic waves. The wave frequency is function of the BC velocity as well as the potential vorticity cross-stream gradient.

Denomina-se Sistema Corrente do Brasil (CB) ao sistema de correntes de contorno oeste formados pela CB, fluindo para sul-sudoeste, e Corrente de Contorno Intermediária (CCI), fluindo para norte-nordeste, associado aos meandros e vórtices. O sistema bordeja a margem continental brasileira sudeste ao longo de 1200-1800 m. A investigação da dinâmica das ondas baroclínicas de vorticidade superpostas ao Sistema Corrente do Brasil é conduzida através exclusivamente de análise de dados hidrográficos de dois cruzeiros de meso-escala (verão e inverno de 1993) do projeto Circulação Oceânica da Região Oeste do Atlântico Sul - COROAS, a componente brasileira do World Ocean Circulation Experiment - WOCE. Estes cruzeiros foram realizados na porção central do Embaiamento de São Paulo (\'23,5 GRAUS\' - \'27 GRAUS\'S). Objetivando-se a aplicação do Método Dinâmico, estimativa quantitativa de um Nível de Referência (NR) é conduzida comparando-se duas metodologias diferentes e independentes. A primeira utiliza o fato de que a CB transporta Água Tropical (AT) e Água Central do Atlântico Sul (ACAS) e de que a CCI transporta dominantemente Água Intermediária Antártica (AIA). Aplica-se então um Critério Termodinâmico, no qual a interface que separa ACAS e AIA é estimada para os dois cruzeiros via aplicação dos Teoremas de Shtokman. A outra metodologia consiste na obtenção de padrões verticais de velocidade baroclínica absoluta a partir de simulações numéricas com a versão seccional do Princenton Ocean Model - POM inicializados com 14 campos termohalinos interpolados a partir dos dados - é o Critério Dinâmico. Os dois critérios apresentam resultado médio que difere apenas cerca de 7 dbar um do outro. Assim, adota-se valor de 480 dbar como NR para os cálculos geostróficos como nível isobárico que demarca a interface média entre CB e CCI. Mapas de função de corrente geostrófica são gerados através de mapeamento objetivo. Tais mapas apresentam evidências inequívocas de que os ciclones e anticiclones, observados nos experimentos de verão e inverno de 1993, são estruturas de vórtices associadas a meandros do Sistema CB. Este aparecem como estruturas vorticais postadas em lados diametralmente opostos ao eixo da CB: ciclones no lado costeiro, e anticiclones no lado oceânico. A variação vertical das estruturas dos meandros e a existência de uma única inversão de sentido permite comprovar e concluir que são ondas baroclínicas de primeiro modo, corroborando especulações presentes na literatura. A construção de um modelo quase-geostrófico de duas camadas para um oceano não-viscoso no plano \'beta\' e de fundo plano permite a realização de análise de vorticidade potencial. Tal análise permite concluir que o campo de vorticidade básico devido à CB suplanta o planetário e que estas ondas baroclínicas são ondas que devem sua existência ao cisalhamento vertical e horizontal da corrente. A variação frontal em vorticidade potencial baroclínica do Sistema CB é de \'1,7x10 POT.-5 s POT. -1\'. A dominãncia do termo da vorticidade de estiramento, que responde por 60% da variação de vorticidade na frente, permite estabelecer que as ondas capturadas na malha hidrográfica são ondas longas dentro da classe de meso-escala. O princípio de conservação de vorticidade potencial é invocado numa análise que envolve a superposição dos campos de vorticidade potencial e função de corrente para a primeira camada. Esta análise apresenta tanto evidências robustas de propagação das ondas baroclínicas quanto indícios de crescimento do meandro ciclônico, ou seja, de processo de instabilidade. Face aos achados das análises anteriores, teoria linear e método das perturbações são utilizados para obter uma relação de dispersão para as ondas baroclínicas da CB. Esta relação é dependente da velocidade da CB e do gradiente de vorticidade potencial através da corrente / The Brazil Current System is formed by two western boundary currents that flow along the Southeast Brazil continental margin from surface to 1200-1800 m deep. These two currents are the southward-flowing Brazil Current (BC) and the northward-flowing Intermediate Western Boundary Current (IWBC). The investigation of the dynamics of the barociinic vorticity waves superimposed to the BC system is conducted in this work solely from hydrographic data anaiysis and manipulation from two oceanographic meso-scale surveys that were part of the COROAS Project, the Brazilian arm of the Worid Ocean Circulation Experiment. These two cruises sampled the central portion of the São Paulo Bight (23,5°- 27°S) in the summer and winter seasons of 1993. As it was intended to appiy the classical Dynamic Method to the data, a quantitative estimate of the reference levei (RL) is conducted through the comparison of two different and independent methodoiogies. The first uses the previous knowiedge that the BC transports Tropical Water and South Atlantic Central Water (SACW) as well as that the IWBC transports mainly Antarctic Intermediate Water (AAIW). This method, referred here as the Thermodynamic Criterion applies the Shtokman theorems to estimate the interface depth between SACW and AAIW. The second methodology, designated here as the Dynamic Criterion, consists of modeiing absolute baroclinic velocities for all 14 hydrographic transects using the sectional version of the Princeton Ocean Model. The results of the two methods differ oniy in 7 dbar and a RL of 480 dbar is adopted as to represent the average interface isobaric level between BC and IWBC in the following geostrophic calculations. Horizontal distributions of geostrophic stream function are generated using objective mapping. These distributions present unequivocal evidences of cyclonic and anticyclonic structures in both Summer and Winter 1993 cruises. These features are associated to meanders of the BC System. They are depicted in opposing sides of the current core. The cyclone is seen in the BC coastal side as the anticyclone is placed in its oceanic side. The vertical variation of these pattems with a single flow direction inversion confirms that these vortical features are part of a first baroclinic mode wave, which corroborates previous speculations found in the literature. In order to pursue a potential vorticity analysis, a data-derived two-layer quasi-geostrophic model is built assuming an invicid flat-bottomed ocean in the 3-plane. This analysis allowed to conclude that the basic vorticity field associated with both horizontal and vertical shear of the boundary currents are responsible for the baroclinic wave existence. The planetary vorticity gradient is one order of magnitude lower. The frontal variation in baroclinic potential vorticity is 1,7 x i0 s1. The dominance of the stretching vorticity, which accounts for 60% of the gradient variation, leads to classify these oscillatory motions as long meso-scale waves. The potential vorticity conservation principle is invoked on an analysis that consists of superimposing the first layer quasi-gesotrophic potential vorticity and geostrophic stream funtion maps. This analysis revealed that the baroclinic waves are propagating as well as evidences of meander growth, an indication of a possible geophysical instability mechanism, are seen. Given the findings of the previous analyses, linear theory and the perturbation method are used to derive a dispersion relation for the BC System first mode baroclinic waves. The wave frequency is function of the BC velocity as well as the potential vorticity cross-stream gradient.

Traffic light control (TLC) with transit signal priority (TSP) is an effective way to deal with urban congestion and travel delay. The growing amount of available connected vehicle data offers opportunities for signal control with transit priority, but the conventional control algorithms fall short in fully exploiting those datasets. This paper proposes a novel approach for dynamic TLC with TSP at an urban intersection. We propose a deep reinforcement learning based framework JenaRL to deal with the complex real-world intersections. The optimisation focuses on TSP while balancing the delay of all vehicles. A two-layer state space is defined to capture the real-time traffic information, i.e. vehicle position, type and incoming lane. The discrete action space includes the optimal phase and phase duration based on the real-time traffic situation. An intersection in the inner city of Jena is constructed in an open-source microscopic traffic simulator SUMO. A time-varying traffic demand of motorised individual traffic (MIT), the current TLC controller of the city, as well as the original timetables of the public transport (PT) are implemented in simulation to construct a realistic traffic environment. The results of the simulation with the proposed framework indicate a significant enhancement in the performance of traffic light controller by reducing the delay of all vehicles, and especially minimising the loss time of PT.

Islet transplantation has found its niche in diabetes treatment. It has contributed to a better quality of life and better glycemic control of patients with diabetes suffering from severe hypoglycemia that are not eligible for vascularized pancreas transplantation. Islet isolation is a technically challenging procedure. The different studies within this doctoral thesis aim to improve and standardize different steps in the isolation procedure. They are in particular looking to improve human pancreas preservation during cold storage, to optimize islet release from the exocrine tissue and to assess whether the isolated islet yield can be predicted from a biopsy. We found that pancreas preservation with pre-oxygenated perfluorodecalin (two-layer method) did not improve the ischemic tolerance of the human pancreas as compared to cold storage with the University of Wisconsin (UW) solution. Furthermore, in pancreas with long cold ischemia time (CIT) (&gt;10 hours), Histidine-Tryptophan-Ketoglutarate (HTK) had a limited preservation capacity as compared with the UW solution with respect to isolation outcome. We also found that during enzymatic pancreas digestion, Vitacyte HA was able to provide a similar islet yield and quality as Serva NB1 with less collagenase activity and shorter digestion time. We further describe the first experience with a new GMP manufactured enzyme called Liberase MTF-S for successful human islet isolation. Finally, we found that the isolated islet yield could not be predicted from a biopsy taken from the head of the pancreas concerning solely morphological parameters of the islets tissue. The improvement of pancreas preservation will allow for marginal organs with prolonged cold ischemia time to expand the donor pool. Better knowledge of how the pancreatic extracellular matrix is digested by collagenase will lead to a fast and predictable islet release from the exocrine tissue. By standardizing the isolation procedure and improving organ selection we will increase the success rate in human islet isolation, thereby making islet transplantation available for more patients.

Ce manuscrit porte sur la modélisation et la simulation numérique d’écoulements diphasiques compressibles. Dans ce contexte, les méthodes d’interfaces diffuses sont aujourd’hui bien acceptées. Cependant, un progrès est encore attendu en ce qui concerne la précision de la capture numérique de ces interfaces. Une nouvelle méthode est développée et permet de réduire significativement cette zone de capture. Cette méthode se place dans le contexte des méthodes numériques de type “MUSCL”, très employées dans les codes de production, et sur maillages non-structurés. Ces interfaces pouvant être le lieu où une transition de phase s’opère, celle-ci est considérée au travers d’un processus de relaxation des énergies libres de Gibbs. Un nouveau solveur de relaxation à thermodynamique rapide est développé et s’avère précis, rapide et robuste y compris lors du passage vers les limites monophasiques. En outre, par rapport aux applications industrielles envisagées, une extension de la thermodynamique des phases et du mélange est nécessaire. Une nouvelle équation d’état est développée en conséquence. La formulation est convexe et est basée sur l’équation d’état “Noble-Abel-Stiffened-Gas”. Enfin, sur un autre plan la dispersion de fluides non-miscibles sous l’effet de la gravité est également abordée. Cette problématique fait apparaître de larges échelles de temps et d’espace et motive le développement d’un nouveau modèle multi-fluide de type “shallow water bi-couche”. Sa résolution numérique est également traitée / This manuscript addresses the theoretical modeling and numerical simulation of compressible two-phase flows. In this context, diffuse interface methods are now well-accepted but progress is still needed at the level of numerical accuracy regarding their capture. A new method is developed in this research work, that allows significant sharpening. This method can be placed in the framework of MUSCL-type schemes, widely used in production codes and on unstructured grids. Phase transition is addressed as well through a relaxation process relying on Gibbs free energies. A new instantaneous relaxation solver is developed and happens to be accurate, fast and robust. Moreover, in view of the intended industrial applications, an extension of the thermodynamics of the phases and of the mixture is necessary. A new equation of state is consequently developed. The formulation is convex and based on the “Noble-Abel-Stiffened-Gas” equation of state. In another context, the dispersion of non-miscible fluids under gravity effects is considered as well. This problematic involves large time and space scales and has motivated the development of a new multi-fluid model for “two-layer shallow water” flows. Its numerical resolution is treated as well