Efter introduktionen av Eurokoderna har mängden minimiarmering i betongkonstruktioner ökat. Denna ökning beror på skillnader i metoderna som nu används för att beräkna mängden minimiarmering i olika typer av betongkonstruktioner och de som användes tidigare då BKR var den gällande normen. Minimiarmering används i betongkonstruktioner för att omfördela tvångsdeformationer. I detta arbete undersöks om mängden minimiarmering kan minskas utan att dess huvudsakliga funktion går förlorad. I arbetet har ett antal metoder för att beräkna minimiarmering jämförts för att se hur stora skillnaderna i armeringsmängd blir, varför dessa uppkommer och vilka för- och nackdelar som finns med metoderna. De undersökta metoderna inkluderar den som anges i Eurokod 2, förändringar som gjorts i den tyska nationella bilagan till Eurokod 2 och ett förslag till ändring som presenterats av Hallgren i ett preliminärt arbete.   Resultaten bygger på simuleringar utförda i FEM-programmet Atena där olika typfall undersökts. I typfallen har betongens hållfasthetsklass, tvärsnittsdimension, armeringsdimension, täckande betongskikt och vidhäftningsförutsättning varierats. Armeringsmängden varierades för att motsvara den mängd som krävs vid beräkning enligt de tidigare nämnda metoderna. I alla undersökta fall antas tvärsnittet vara utsatt för ren dragbelastning, vilket är fallet vid krympning. Resultaten från de numeriska simuleringarna visade att en större armeringsmängd ger minskade sprickbredder och ökade möjligheter att omfördela tvångsdeformationer till nya sprickor. Bildandet av nya sprickor sker också vid en mindre krympning. Skillnaderna jämfört med om en mindre armeringsmängd används blir dock i de flesta fall små och mängden armering kan reduceras utan att minimiarmeringens sprickfördelande funktion äventyras. Resultaten från de numeriska simuleringarna användes också för att ta fram ett alternativ som ger en mindre mängd minimiarmering utan att armeringens huvudsakliga funktion försvinner. I detta arbete föreslås en förändring rörande koefficienten k i Eurokod. Förändringen av värdet på koefficienten k gör att mängden erforderlig minimiarmering minskar betydligt vid beräkningar, speciellt för tvärsnitt med en större tvärsnittshöjd. Fler numeriska simuleringar utfördes och det kunde säkerställas att en armeringsmängd motsvarande den som ges med den föreslagna ändringen fungerar för alla undersökta typfall. Endast en liten ökning av sprickbredderna sker.   En undersökning av sprickproblematiken under betongens hållfasthetstillväxtfas har gjorts, och resultaten visar att problem med sprickbildning inte beror på den autogena krympningen eller på uttorkningskrympningen. Endast en liten andel av denna krympning hinner uppnås under de första veckorna, samtidigt som hållfasthetstillväxten sker relativt snabbt. Problem med sprickbildning under denna tidsperiod beror istället på andra faktorer såsom en stor värmeutveckling i betongen, en ojämn uttorkning eller en ojämn temperaturfördelning över tvärsnittet i kombination med tvång. Efter den värmeutveckling som skett i betongen under det första dygnet efter gjutning påbörjas en avsvalningsfas. Under denna kan sprickbildning uppstå i betongen, och denna sprickbildning finns kvar när uttorkningskrympningen senare fortskrider. En jämförelse har även gjorts mellan handberäkningar av karaktäristiska sprickbredder enligt Betongföreningens handbok till Eurokod 2 och sprickbredder utlästa ur resultaten från de numeriska simuleringarna. Resultaten visar att metoden överlag fungerar bra, men att vissa justeringar skulle kunna göras för att anpassa den till mindre värden på betongens krympning. / After the introduction of the Eurocodes, the minimum amount of reinforcement for crack control in concrete structures has increased. This is due to differences in themethods for calculating the minimum amount of reinforcement used in the Eurocodes and in BKR, a standard which was used in Sweden before the Eurocodes. Minimum reinforcement for crack control is used in concrete structures to redistribute the tensile stresses caused by restraint during the shrinkage. This thesis examines if the amount of minimum reinforcement can be reduced without compromising the main functions of the reinforcement. In this thesis a number of methods for calculating the minimum reinforcement has been compared to see how large the differences in the amount of reinforcement are, why they occur and what are the advantages and disadvantages of the different methods. The different methods revised in this thesis are the method in Eurocode 2, the changes made in the German National Annex to the Eurocode 2 method and a method proposed by Hallgren in a preliminary work.   The results are based on simulations performed in the FEM-application Atena where different scenarios were examined. In the investigated scenarios, the strength classes of the concrete, the cross-sectional dimensions, the size of the reinforcement bars, the size of the concrete cover and the bond strength has been varied. The amount of reinforcement in each case has been varied to correspond to the required minimum amount according to the different methods for calculating minimum reinforcement. In all investigated cases, the cross section is assumed to be exposed to pure tensile load, which is the case for shrinkage. The results of the simulations showed that a larger amount of reinforcement reduces the width of the cracks, that a larger amount of small cracks are formed, and that new cracks are formed earlier at a lower shrinkage value. The differences are however in most cases small compared to when a lower amount of reinforcement is used. In other words, a lower amount of reinforcement could be used without compromising the main functions of the minimum reinforcement. The results from the simulations have also been used to obtain an alternative method for calculating the minimum reinforcement that gives a lower amount of reinforcement without compromising the main functions. The proposed change in Eurocode is a modification in the value of the coefficient k. By changing the value of the coefficient k, the amount of required minimum reinforcement is decreased significantly, especially for large crosssection heights. More numerical simulations were performed to ensure that the reinforcement amount according to the proposed change works for all scenarios considered in this thesis. The results showed that the reinforcement were capable of redistributing the tensile stresses caused by restraint. Only a small increase in the crack widths could be observed.   An investigation was conducted to find an explanation to why cracks occur during the first weeks after casting of the concrete. The results show that the problem with cracks occurring during this period is neither due to the autogenous shrinkage nor the drying shrinkage. Only a small share of those types of shrinkage are achieved during the first weeks, and at the same time the concrete strength growth is relatively fast. Instead, the problems with cracks seem to be caused by other factors, such as a large heat development in the concrete during the first day, uneven drying shrinkage or an uneven temperature distribution over the cross section in combination with restraint. Cracks can develop during the first days after casting of the concrete due to the shrinkage that occur when the concrete cools off after the large heat development the first day. These cracks remain when the drying shrinkage continues. A comparison was also made between hand calculations of the characteristic crack widths according to “Betongföreningens handbok till Eurokod 2” and the crack widths read out of the results from the numerical simulations. The results showed that the method generally works well, but some minor adjustments could be made to adapt the method to smaller values of concrete shrinkage.

Le domaine des problèmes quantiques à N-corps à l'équilibre et hors d'équilibre sont des sujets majeurs de la Physique et de la Physique de la matière condensée en particulier. Les propriétés d'équilibre de nombreux systèmes unidimensionnels en interaction sont bien comprises d'un point de vue théorique, des chaînes de spins aux théories quantiques des champs dans le continue. Ces progrès ont été rendus possibles par le développement de nombreuses techniques puissantes, comme, par exemple, l'ansatz de Bethe, le groupe de renormalisation, la bosonisation, les états produits de matrices ou la théorie des champs invariante conforme. Même si les propriétés à l'équilibre de nombreux modèles soient connues, ceci n'est en général pas suffisant pour décrire leurs comportements hors d'équilibre, et ces derniers restent moins explorés et beaucoup moins bien compris. Les modèles d'impuretés quantiques représentent certains des modèles à N-corps les plus simples. Mais malgré leur apparente simplicité ils peuvent capturer plusieurs phénomènes expérimentaux importants, de l'effet Kondo dans les métaux aux propriétés de transports dans les nanostructures, comme les points quantiques. Dans ce travail nous considérons un modèle d'impureté appelé "modèle de niveau résonnant en interaction" (IRLM). Ce modèle décrit des fermions sans spin se propageant dans deux fils semi-infinis qui sont couplés à un niveau résonant -- appelé point ou impureté quantique -- via un terme de saut et une répulsion Coulombienne. Nous nous intéressons aux situations hors d'équilibre où un courant de particules s'écoule à travers le point quantique, et étudions les propriétés de transport telles que le courant stationnaire (en fonction du voltage), la conductance différentielle, le courant réfléchi, le bruit du courant ou encore l'entropie d'intrication. Nous réalisons des simulations numériques de la dynamique du modèle avec la méthode du groupe de renormalisation de la matrice densité dépendent du temps (tDMRG), qui est basée sur une description des fonctions d'onde en terme d'états produits de matrices. Nous obtenons des résultats de grande précision concernant les courbes courant-voltage ou bruit-voltage de l'IRLM, dans un grand domaine de paramètres du modèle (voltage, force de l'interaction, amplitude de saut vers le dot, etc.). Ces résultats numériques sont analysés à la lumière de résultats exacts de théorie des champs hors d'équilibre qui ont été obtenus pour un modèle similaire à l'IRLM, le modèle de Sine-Gordon avec bord (BSG). Cette analyse est en particulier basée sur l'identification d'une échelle d'énergie Kondo et d'exposants décrivant les régimes de petit et grand voltage. Aux deux points particuliers où les modèles sont connus comme étant équivalents, nos résultats sont en accord parfait avec la solution exacte. En dehors de ces deux points particuliers nous trouvons que les courbes de transport de l'IRLM et du modèle BSG demeurent très proches, ce qui était inattendu et qui reste dans une certaine mesure inexpliqué. / The fields of in- and out-of-equilibrium quantum many-body systems are major topics in Physics, and in condensed-matter Physics in particular. The equilibrium properties of one-dimensional problems are well studied and understood theoretically for a vast amount of interacting models, from lattice spin chains to quantum fields in a continuum. This progress was allowed by the development of diverse powerful techniques, for instance, Bethe ansatz, renormalization group, bosonization, matrix product states and conformal field theory. Although the equilibrium characteristics of many models are known, this is in general not enough to describe their non-equilibrium behaviors, the latter often remain less explored and much less understood. Quantum impurity models represent some of the simplest many-body problems. But despite their apparent simplicity, they can capture several important experimental phenomena, from the Kondo effect in metals to transport in nanostructures such as point contacts or quantum dots. In this thesis consider a classic impurity model - the interacting resonant level model (IRLM). The model describes spinless fermions in two semi-infinite leads that are coupled to a resonant level -- called quantum dot or impurity -- via weak tunneling and Coulomb repulsion. We are interested in out-of-equilibrium situations where some particle current flows through the dot, and study transport characteristics like the steady current (versus voltage), differential conductance, backscattered current, current noise or the entanglement entropy. We perform extensive state-of-the-art computer simulations of model dynamics with the time-dependent density renormalization group method (tDMRG) which is based on a matrix product state description of the wave functions. We obtain highly accurate results concerning the current-voltage and noise-voltage curves of the IRLM in a wide range parameter of the model (voltage bias, interaction strength, tunneling amplitude to the dot, etc.).These numerical results are analyzed in the light of some exact out-of-equilibrium field-theory results that have been obtained for a model similar to the IRLM, the boundary sine-Gordon model (BSG).This analysis is in particular based on identifying an emerging Kondo energy scale and relevant exponents describing the high- and low- voltage regimes. At the two specific points where the models are known to be equivalent our results agree perfectly with the exact solution. Away from these two points, we find that, within the precision of our simulations, the transport curves of the IRLM and BSG remain very similar, which was not expected and which remains somewhat unexplained.

No description available.

No description available.

[ES] La descarga de vapor en forma de chorro en piscinas con agua subenfriada es una forma muy eficiente de condensar vapor y es muy utilizada en diversas aplicaciones industriales, como la industria farmacéutica, metalúrgica, química, sistemas de propulsión sumergida y nuclear. Siendo especialmente esencial en los sistemas de seguridad de las centrales nucleares, concretamente en la piscina de supresión de presión de los Reactores de Agua en Ebullición (BWR, Boiling Water Reactor) y en el tanque de almacenamiento de agua de reabastecimiento en la contención de los reactores avanzados de agua a presión. La eficiencia de esta técnica, denominada Condensación por Contacto Directo (DCC, Direct Contact Condensation), radica en la rápida condensación del vapor proporcionando una alta capacidad de transferencia de calor e intercambio de masa. En las últimas décadas, se han llevado a cabo muchos experimentos sobre los chorros sumergidos de gases no condensables y vapor puro en piscinas, que han proporcionado mucha información de interés, aunque continúan los esfuerzos para adquirir información adicional. Para estudiar en detalle el comportamiento de las descargas de chorros turbulentos se han llevado a cabo una serie de experimentos en una instalación que ha sido diseñada para representar diversos entornos de interés crítico para la seguridad en centrales nucleares. Su singularidad radica en la capacidad de poder reproducir fenomenologías similares a situaciones específicas que podrían surgir en una central nuclear, desde sistemas de despresurización hasta roturas con pérdida de refrigerante. La instalación posibilita el estudio con vapor, aire o una mezcla de ambos utilizando geometrías simples (boquilla) o complejas (sparger) en una piscina con agua en reposo. Se ha utilizado una cámara de alta velocidad con Dispositivo de Carga Acoplada (CCD, Charged-Coupled Device) para capturar imágenes de la descarga del chorro en la piscina y así determinar la interfase entre el gas y el líquido, utilizando técnicas de visualización directa. Para llevar a cabo el procesamiento de dichas imágenes, se han aplicado una serie de procedimientos de procesado, filtrado y post-procesado utilizando una subrutina implementada en MATLAB. Tras la adquisición de las imágenes, se ha caracterizado el comportamiento del chorro de vapor y aire, el efecto de la fracción volumétrica de aire sobre la mezcla y el efecto de la variación del tamaño de la boquilla. Además, se ha medido la longitud del momento del chorro y se ha comprobado su estrecha relación con el diámetro de la boquilla, la velocidad del chorro y el porcentaje de la mezcla. Por otro lado, se ha caracterizado el comportamiento del chorro horizontal en la zona de transición. La transición indica el cambio de un chorro dominado exclusivamente por la fuerza de momento a un régimen en el que las fuerzas de flotación prevalecen. Esta prueba se ha realizado utilizando cinco boquillas intercambiables y para cada caudal de vapor se ha determinado el porcentaje de aire en el que ocurre la transición. Otro de los principales objetivos de esta tesis es la modelación utilizando herramientas de códigos de Dinámica de Fluidos Computacional (CFD), específicamente descargas de chorros de aire en agua subenfriada. Para esto se ha empleado el software ANSYS CFX. El comportamiento del flujo multifásico se ha simulado utilizando el modelo euleriano-euleriano no homogéneo; debido a que es más conveniente para estudiar fluidos interpenetrados. Los resultados CFD se han validado con datos experimentales y resultados teóricos. Las cifras de mérito de la longitud de penetración adimensional y la longitud de flotabilidad adimensional muestran una buena concordancia con las mediciones experimentales. Las correlaciones para estas variables se obtuvieron en función de números adimensionales para dar generalidad utilizando sólo condiciones iniciales de contorno. / [CA] La descàrrega de vapor en forma de xorrat en piscines amb aigua subrefredada és una forma molt eficient de condensar vapor i s'utilitza àmpliament en diverses aplicacions industrials, com la indústria farmacèutica, metal·lúrgica, química, sistemes de propulsió submergida i nuclear. És especialment essencial en els sistemes de seguretat de les centrals nuclears, concretament en la piscina de supressió de pressió dels Reactors d'Aigua en Bulliment (BWR) i en el dipòsit d'aigua de reaprovisionament a la contenció dels reactors avançats d'aigua a pressió. L'eficiència d'aquesta tècnica, anomenada Condensació per Contacte Directe (DCC), rau en la ràpida condensació del vapor proporcionant una alta capacitat de transferència de calor i intercanvi de massa. En les últimes dècades, s'han realitzat molts experiments sobre els raigs submergits de gasos no condensables i vapor pur en piscines, que han proporcionat molta informació d'interès, tot i que continuen els esforços per adquirir informació addicional. Per estudiar amb detall el comportament de les descàrregues de raigs turbulents s'han dut a terme una sèrie d'experiments en una instal·lació que ha estat dissenyada per representar diversos entorns d'interès crític per a la seguretat en centrals nuclears. La seua singularitat rau en la capacitat de poder reproduir fenomenologies semblants a situacions específiques que podrien sorgir en una central nuclear, des de sistemes de despressurització fins a ruptures amb pèrdua de refrigerant. La instal·lació possibilita l'estudi amb vapor, aire o una barreja de tots dos utilitzant geometries simples (busseta) o complexes (sparger) en una piscina amb aigua en repòs. A més, és possible combinar aquestes situacions amb diferents condicions de cabal, pressió o temperatura. S'ha utilitzat una càmera de alta velocitat (CCD) per capturar imatges de la descàrrega del raig a la piscina i així determinar la interfase entre el gas i el líquid, emprant tècniques de visualització directa. Per dur a terme el processament d'aquestes imatges, s'han aplicat una sèrie de procediments complexes de processat, filtrat i post-processat mitjançant una subrutina implementada en MATLAB. Després de l'adquisició de les imatges, s'ha caracteritzat el comportament del raig de vapor i aire, l'efecte de la fracció volumètrica d'aire sobre la barreja i l'efecte de la variació de la grandària de la busseta. A més, s'ha mesurat la longitud del moment del raig i s'ha comprovat la seua estreta relació amb el diàmetre de la busseta, la velocitat del raig i el percentatge de la barreja. D'altra banda, s'ha caracteritzat el comportament del raig horitzontal en la zona de transició. La transició indica el canvi d'un raig dominat exclusivament per la força de moment a un règim en què les forces de flotació prevalen. Aquesta prova s'ha realitzat utilitzant cinc bussetes intercanviables i per a cada cabal de vapor s'ha determinat el percentatge d'aire en què ocorre la transició. Un altre dels principals objectius d'aquesta tesi és la modelització utilitzant eines de codis de Dinàmica de Fluids Computacionals (CFD), específicament descàrregues de xorrades d'aire en aigua subrefredada. Per a això s'ha utilitzat el software ANSYS CFX. El comportament del flux multifàsic s'ha simulat utilitzant el model eulerià-eulerià no homogeni, ja que és més convenient per estudiar fluids interpenetrats. Els resultats CFD s'han validat amb dades experimentals i resultats teòrics. Les figures de mèrit de la longitud de penetració adimensional i la longitud de flotabilitat adimensional mostren una bona concordança amb les mesures experimentals. Les correlacions per a aquestes variables s'han obtingut en funció de nombres adimensionals per donar generalitat utilitzant només condicions inicials de contorn. El model numèric CFD desenvolupat té la capacitat de simular el comportament de gasos no condensables descarregats en aigua. / [EN] The discharge of steam in the form of a jet in pools with subcooled water is a highly efficient way to condense steam and is widely used in various industrial applications, such as the pharmaceutical, metallurgical, chemical, submerged propulsion, and nuclear industries. It is particularly essential in the safety systems of nuclear power plants, specifically in the pressure suppression pool of Boiling Water Reactors (BWR) and in the water replenishment tank in the containment of advanced pressurized water reactors. The efficiency of this technique, known as Direct Contact Condensation (DCC), lies in the rapid condensation of steam, providing a high heat transfer and mass exchange capacity. In recent decades, numerous experiments have been conducted on submerged jets of non-condensable gases and pure steam in pools, yielding valuable information, although efforts continue to gather additional data. To study in detail the behavior of turbulent jet discharges, a series of experiments have been carried out in a facility designed to simulate various critical environments for safety in nuclear power plants. Its uniqueness lies in the ability to reproduce phenomenologies like specific situations that could arise in a nuclear power plant, from depressurization systems to ruptures with coolant loss. The facility enables the study with steam, air, or a mixture of both using simple geometries (nozzle) or complex geometries (sparger) in a pool with quiescent water. Additionally, it is possible to combine these situations with different flow, pressure, or temperature conditions. A high-speed camera with Charge-Couple Device (CCD) has been used to capture images of the jet discharge in the pool, aiming to determine the gas-liquid interface through direct visualization techniques. Following the image acquisition, a series of complex processing, filtering, and post-processing procedures have been applied using a subroutine implemented in MATLAB. After obtaining the images, the behavior of the steam and air jet, the effect of the volumetric fraction of air on the mixture, and the impact of varying the nozzle size have been characterized. Additionally, the length of the jet momentum has been measured, establishing its close relationship with the nozzle diameter, jet velocity, and mixture percentage. Furthermore, the behavior of the horizontal jet in the transition zone has been characterized. The transition denotes the shift from a jet dominated exclusively by momentum force to a regime where buoyancy forces prevail. This test has been conducted using five interchangeable nozzles, and for each steam flow rate, the percentage of air at which the transition occurs has been determined. Another main objective of this thesis is modeling using Computational Fluid Dynamics (CFD) tools, specifically jet discharges of air into subcooled water. The ANSYS CFX software has been employed for this purpose. Multiphase flow behavior has been simulated using the non-homogeneous Eulerian-Eulerian model, as it is more convenient for studying interpenetrated fluids. CFD results have been validated with experimental data and theoretical results. The dimensionless penetration length and dimensionless buoyancy length merit figures show good agreement with experimental measurements. Correlations for these variables were obtained as a function of dimensionless numbers to provide generality using only initial boundary conditions. The developed CFD numerical model has the capability to simulate the behavior of non-condensable gases discharged into water. / Córdova Chávez, Y. (2024). Experimentación y modelación de mezcla de vapor y gases no-condensables [Tesis doctoral]. Universitat Politècnica de València. https://doi.org/10.4995/Thesis/10251/207266

Les simulations numériques sont devenues un outil indispensable dans la conception des chambres de coupure des disjoncteurs électriques haute tension. Elles sont utilisées non seulement dans le dimensionnement des différentes pièces, mais elles fournissent également une aide précieuse dans la compréhension des phénomènes intervenant entre les deux électrodes au moment de la coupure. L’arc électrique généré entre ces deux électrodes rassemble de nombreux domaines de la physique plus ou moins complexes. Tous ces phénomènes ne sont pas encore parfaitement compris. Avec l’évolution de la puissance de calcul, ces simulations peuvent prendre en compte de plus en plus de phénomènes. Mais pour des raisons de temps de développement, la question des phénomènes à prendre en compte dans ces simulations se pose. Le but de telles simulations est de déterminer de manière rapide si une configuration est plus ou moins capable qu’une autre de couper sous une contrainte donnée. Ainsi, il est important de prendre en compte uniquement les phénomènes physiques importants et nécessaires pour avoir une réponse la plus décisive possible et la plus rapide possible, de la réussite ou non à la coupure d’une configuration testée. Dans cette thèse, nous nous sommes particulièrement intéressés aux déséquilibres thermiques et chimiques qui pourraient intervenir dans les disjoncteurs électriques haute tension au moment de la coupure. En effet, pour des raisons de temps et de coût de calcul, la plupart des simulations numériques actuelles sont réalisées en faisant une hypothèse forte : l’hypothèse d’Equilibre Thermodynamique Local (ETL). Cette hypothèse consiste à considérer que dans chaque maille de notre domaine d’étude et à chaque pas de temps, on a un équilibre thermodynamique réalisé. Faire cette hypothèse nous permet d’utiliser les lois de conservation (masse, quantité de mouvement et énergie) en allégeant le problème. Mais en réalité, cette hypothèse est mise à mal dès que l’on est en présence de forts gradients de température ou de densité. Pour réaliser ces simulations, le code numérique CARBUR a été utilisé. Des modules d’arc électrique (effet Joule et rayonnement) et d’électrode mobile ont été implémentés afin de pouvoir simuler au mieux le comportement du gaz présent dans les disjoncteurs électriques haute tension. Six études différentes ont été réalisées et sont présentées. Ces études portent sur les influences de la forme du bout des électrodes, d’une modélisation en Navier-Stokes par rapport à une modélisation en Euler, de la nature du gaz (SF6, CO2 et N2), du déséquilibre thermique dans le cas de l’azote ou encore du positionnement des termes sources de l’arc électrique dans les différentes équations d’évolution des énergies. Dans ce travail, une étude sur différents modèles cinétiques chimiques est proposée. Dans ces modèles, 5 espèces chimiques sont présentes : N2, N, N+, N2+ et e-. En ce qui concerne la température, on en distingue 4 : T, TVib-N2, TVib-N2+ et Te. / The numerical simulations are become a very important tool to design the high voltage circuit breaker (HVCB) chamber. They help for the understanding of the different phenomena which can take place between the 2 electrodes during an interruption process. The electric arc brings together many fields of physics more or less complex and many of these phenomena are still poorly studied. So many aspects remain to be explored to improve simulations. With the increase of the calculation power, these numerical simulations can take into account more phenomena. However, for reasonable simulation times, we need to know which phenomena are preponderant. The aim of these numerical simulations is to rapidly conclude on the capacity of geometry to success an interruption process compared to different other geometries, under a given stress. In this PhD dissertation, we are particularly interested on thermal and chemical non equilibrium that can occur in HVCB during an interruption process. Currently, most simulations are carried out with a strong hypothesis: the hypothesis of Local Thermodynamic Equilibrium (LTE). This assumption allows us to alleviate the problem and to reduce the computing time. But this assumption becomes not valid when high temperature or density gradients occur. To do these simulations, the CARBUR numerical code has been used. In order to simulate flow behaviors in HVCB, an electrical arc (Joule effect and radiation) model and a module of mobile electrode have been added. Six different studies have been done and are presented: influence of the electrode shape, influence of the Navier-Stokes equations compared to the Euler equations, influence of the gas (SF6, CO2 et N2), influence of the thermal non equilibrium in a nitrogen case, influence of the position of the arc source terms in the different energy equations. In this work, a study on different nitrogen chemical kinetics is proposed. In these models, 5 chemical species are distinguished: N2, N, N+, N2+ and e-. Finally, 4 different temperatures are used: T, TVib-N2, TVib-N2+ and Te.

Depuis plusieurs décennies, une importante activité scientifique se concentre sur la description numérique, théorique ou expérimentale de l'hydrodynamique des écoulements multiphasiques. Ces écoulements sont caractérisés par l'existence d'interfaces, et d'une force à l'interface, la tension superficielle, séparant généralement deux fluides non miscibles. Un cas d'étude dans ce contexte est le problème du drainage d'une unique goutte dans une phase continue, l'ensemble étant soumis à la gravité. Ce système fait apparaître des écoulements récemment décrits pour une goutte d'eau dans l'huile de paraffine. Ce système constitue également un modèle simple pour l'étude des propriétés aux interfaces, Mais d'un point de vue numérique, se pose alors le problème de la stabilité des algorithmes pouvant être utilisés. Les effets aux interfaces impliquent en effet des domaines spatiaux très limités dans lesquels les grandeurs physiques entre les deux fluides sont discontinues. D'importants artéfacts numériques peuvent alors être générés dans les simulations et faire perdre la richesse de la physico-chimie du système considéré. Le problème de la simulation d'écoulements multiphasiques intéresse aussi bien le monde académique que le monde industriel. L'objectif de ce travail de thèse est donc d'implémenter les techniques numériques les plus récentes et de développer un code pour permettre la simulation de l'hydrodynamique de systèmes dispersés. Pour parvenir à ce but, il reste encore des problèmes algorithmiques importants à résoudre comme la prise en compte des effets thermocapillaires et thermosolutaux. Ces deux derniers points sont l'objet de cette thèse. / For several decades, an important scientific activity has focused on the numerical, theoretical and experimental hydrodynamics of drops. This work presents numerical results of a single droplet in the gravity field and in non-isothermal conditions. The simulation such a multiphase system is important in both academic and industrial world. This is particularly the case in the field of emulsions, wetting problems and evaporation. To achieve this goal, there are still important algorithmic problems due to the free moving interfaces and the description of capillary effects. Here, a Volume of Fluid technique has been implemented with high order temporal and spatial schemes to preserve the sharpness of the drop interface. The system under consideration is a simplified model consisting in a single water droplet in a continuous paraffin oil phase. These liquids are immiscible and non-compressible and the overall evolution is unsteady. Capillary contributions such as temperature and surfactant dependent surface tension are fully accounted for. This presentation is aimed to show the capabilities of VOF techniques for the simulations of unsteady multiphase systems in non-isothermal configurations. The role of the droplet initial position and temperature field is described with good numerical stability. There are still important problems remaining in the simulation of free interface systems with such a technique. Spurious currents induced by the description of capillarity can in particular come into play. But these latter can be controlled once the droplet average velocity due to drainage becomes large enough.

Les cellules à hétérojonctions de silicium fabriquées par croissance de couches minces de silicium amorphe hydrogéné (a-Si :H) à basse température sur des substrats de silicium cristallin (c-Si) peuvent atteindre des rendements de conversion photovoltaïque élevés (η=23 % démontré). Les efforts de recherche ayant principalement été orientés vers le cristallin de type p jusqu'à présent en France, ce travail s'attache à l'étude du type n pour d'une part déterminer les performances auxquelles s'attendre avec cette nouvelle filière et d'autre part les améliorer. Pour cela, nous avons mis en œuvre des techniques de caractérisation des matériaux composant la structure et de l’interface (a-Si :H/c-Si) couplées à des outils de simulations numériques afin mieux comprendre les phénomènes de transport électronique. Nous nous sommes également intéressés aux cellules à hétérojonctions avec substrats de silicium multicristallin de type n, le silicium multicristallin étant le matériau le plus répandu actuellement dans la fabrication des cellules photovoltaïques. / In this thesis we focus on the silicon heterostructure combining thin films amorphous silicon (a-Si :H) deposited at low temperature on crystalline silicon (c-Si) substrates. We study the different materials and the interface between them through both characterizations, modelling and numerical simulations. The goal is to better understand the influence of the different parameters (doping level, defects density, band offset, ...) on the photovoltaic solar cell's performances in order to get them improved. Structures with multicrystalline silicon substrates are also studied.

Cette thèse porte sur la modélisation de l’écoulement et du transport en milieu poreux ;nous effectuons des simulations numériques et démontrons des résultats de convergence d’algorithmes.Au Chapitre 1, nous appliquons des méthodes de volumes finis pour la simulation d’écoulements à densité variable en milieu poreux ; il vient à résoudre une équation de convection diffusion parabolique pour la concentration couplée à une équation elliptique en pression.Nous nous appuyons sur la méthode des volumes finis standard pour le calcul des solutions de deux problèmes spécifiques : une interface en rotation entre eau salée et eau douce et le problème de Henry. Nous appliquons ensuite la méthode de volumes finis généralisés SUSHI pour la simulation des mêmes problèmes ainsi que celle d’un problème de bassin salé en dimension trois d’espace. Nous nous appuyons sur des maillages adaptatifs, basés sur des éléments de volume carrés ou cubiques.Au Chapitre 2, nous nous appuyons de nouveau sur la méthode de volumes finis généralisés SUSHI pour la discrétisation de l’équation de Richards, une équation elliptique parabolique pour le calcul d’écoulements en milieu poreux. Le terme de diffusion peut être anisotrope et hétérogène. Cette classe de méthodes localement conservatrices s’applique àune grande variété de mailles polyédriques non structurées qui peuvent ne pas se raccorder.La discrétisation en temps est totalement implicite. Nous obtenons un résultat de convergence basé sur des estimations a priori et sur l’application du théorème de compacité de Fréchet-Kolmogorov. Nous présentons aussi des tests numériques.Au Chapitre 3, nous discrétisons le problème de Signorini par un schéma de type gradient,qui s’écrit à l’aide d’une formulation variationnelle discrète et est basé sur des approximations indépendantes des fonctions et des gradients. On montre l’existence et l’unicité de la solution discrète ainsi que sa convergence vers la solution faible du problème continu. Nous présentons ensuite un schéma numérique basé sur la méthode SUSHI.Au Chapitre 4, nous appliquons un schéma semi-implicite en temps combiné avec la méthode SUSHI pour la résolution numérique d’un problème d’écoulements à densité variable ;il s’agit de résoudre des équations paraboliques de convection-diffusion pour la densité de soluté et le transport de la température ainsi que pour la pression. Nous simulons l’avance d’un front d’eau douce assez chaude et le transport de chaleur dans un aquifère captif qui est initialement chargé d’eau froide salée. Nous utilisons des maillages adaptatifs, basés sur des éléments de volume carrés. / This thesis bears on the modelling of groundwater flow and transport in porous media; we perform numerical simulations by means of finite volume methods and prove convergence results. In Chapter 1, we first apply a semi-implicit standard finite volume method and then the generalized finite volume method SUSHI for the numerical simulation of density driven flows in porous media; we solve a nonlinear convection-diffusion parabolic equation for the concentration coupled with an elliptic equation for the pressure. We apply the standard finite volume method to compute the solutions of a problem involving a rotating interface between salt and fresh water and of Henry's problem. We then apply the SUSHI scheme to the same problems as well as to a three dimensional saltpool problem. We use adaptive meshes, based upon square volume elements in space dimension two and cubic volume elements in space dimension three. In Chapter 2, we apply the generalized finite volume method SUSHI to the discretization of Richards equation, an elliptic-parabolic equation modeling groundwater flow, where the diffusion term can be anisotropic and heterogeneous. This class of locally conservative methods can be applied to a wide range of unstructured possibly non-matching polyhedral meshes in arbitrary space dimension. As is needed for Richards equation, the time discretization is fully implicit. We obtain a convergence result based upon a priori estimates and the application of the Fréchet-Kolmogorov compactness theorem. We implement the scheme and present numerical tests. In Chapter 3, we study a gradient scheme for the Signorini problem. Gradient schemes are nonconforming methods written in discrete variational formulation which are based on independent approximations of the functions and the gradients. We prove the existence and uniqueness of the discrete solution as well as its convergence to the weak solution of the Signorini problem. Finally we introduce a numerical scheme based upon the SUSHI discretization and present numerical results. In Chapter 4, we apply a semi-implicit scheme in time together with a generalized finite volume method for the numerical solution of density driven flows in porous media; it comes to solve nonlinear convection-diffusion parabolic equations for the solute and temperature transport as well as for the pressure. We compute the solutions for a specific problem which describes the advance of a warm fresh water front coupled to heat transfer in a confined aquifer which is initially charged with cold salt water. We use adaptive meshes, based upon square volume elements in space dimension two.

Ette thèse traite d’écoulements de fluides à seuil (ou viscoplastiques) en milieu confiné. Les difficultés analytiques et numériques sont dues à la multivaluation du tenseur des contraintes dans les zones plastiques ainsi qu’à la non-différentiabilité du problème de minimisation associé. Cette thèse s’articule en deux parties.Dans un premier temps, des simulations numériques parallèles très précises à l’aide d’algorithmes de dualité ont été effectuées. Elles ont permis de retrouver des résultats observés expérimentalement dont l’existence d’une ligne de glissement pour l’écoulement au dessus d’un obstacle et le caractère quasi-Poiseuille de la vitesse au-delà de cette ligne. Par ailleurs, la théorie de couche limite viscoplastique définie par Oldroyd (1947, à nombre de Bingham asymptotiquement grand) a été revisitée à nombre de Bingham modéré en milieu confiné. L’étude a mis en œuvre des allers-retours entre ces simulations et les expériences physiques de Luu et al. d’IRSTEA ainsi qu’une dérivation théorique. L’approximation de couche limite est vérifiée dans une certaine mesure à l’intérieur de la cavité. Une adaptation de la notion de couche limite viscoplastique est alors exhibée et permet d’étendre les scalings dérivés par Oldroyd (1947) et Balmforth et al. (J. of Fluid Mech, 2017). Ces scalings sont aussi généralisés au cas de la loi d’Herschel-Bulkley. Dans un second temps, on présente une analyse asymptotique des champs de vitesses et de contraintes pour des écoulements en faible épaisseur (ε). Un développement à l’ordre ε2 de la vitesse permet de trouver une équation de Reynolds à la même précision. Cette équation de Reynolds prolonge les résultats déjà existants dans le cadre newtonien, d’une part et dans le cadre fluide à seuil avec une surface libre, d’autre part. / This thesis is devoted to the flow of yield stress (or viscoplastic) fluids in pipes.Analytical and numerical difficulties lie in the multivaluation of the stress tensor in the plastic regions and in the non-differentiability of the associated minimization problem. This manuscript is organized following two main axes.First, very accurate numerical simulations were carried out using duality methods and parallel multifrontal solvers. Thus, experimental observations were recovered, namely the existence of a slip line for the flow over an obstacle and the Poiseuille-like behaviour of the velocity above this line. Moreover, the viscoplastic boundary layer theory defined by Oldroyd (1947 at high Bingham numbers) was revisited at moderate Bingham numbers in confined areas. This study provided an opportunity to go back and forth between these simulations and the physical measures of Luu et al. from IRSTEA and to perform a theoretical derivation. The boundary layer approximation is valid up to a certain extent in the cavity. An adaptation of the viscoplastic boundary layer definition is then given and allows to generalize the scalings shown by Oldroyd (1947) and Balmforth et al. (JFM 2017). These scalings are also generalized to the Herschel-Bulkley case. Then, an asymptotic analysis of the velocity and stress fields for thin layer (ε) flows is presented. A velocity development up to ε2 lets us find a Reynolds equation of same accuracy. This Reynolds equation extends the already existing results, on the one hand in the newtonian case and on the second hand for free surface flows.