Direct Contact Condensation between steam and water as well as bubble entrainment below the water surface play an important role in different accident scenarios for light water reactors. One example is the emergency core cooling water injection into a two-phase mixture. It has to be considered for example to evaluate potential pressurized thermal shock phenomena. This report documents experiments conducted in flat basin inside the TOPFLOW pressure chamber aiming on the generation of a database useful for CFD model development and validation. It comprises 3 different setups: condensation at a stratified flow of sub-cooled water, condensation at a sub-cooled water jet and a combination of both phenomena with steam bubble entrainment. The documentation includes all details on the experimental set up, on experimental conditions (experimental matrices), on the conduction of the experiments, on measuring techniques used and on data evaluation procedures. In addition, selected results are presented.

Counter-Current Flow Limitation (CCFL) is of importance for PWR safety analyses in several accident scenarios connected with loss of coolant. Basing on the experiences obtained during a first series of hot leg tests now new experiments on counter-current flow limitation were conducted in the TOPFLOW pressure vessel. The test series comprises air-water tests at 1 and 2 bar as well as steam-water tests at 10, 25 and 50 bar. During the experiments the flow structure was observed along the hot leg model using a high-speed camera and web-cams. In addition pressure was measured at several positions along the horizontal part and the water levels in the reactor-simulator and steam-generator-simulator tanks were determined. This report documents the experimental setup including the description of operational and special measuring techniques, the experimental procedure and the data obtained. From these data flooding curves were obtained basing on the Wallis parameter. The results show a slight shift of the curves in dependency of the pressure. In addition a slight decrease of the slope was found with increasing pressure. Additional investigations concern the effects of hysteresis and the frequencies of liquid slugs. The latter ones show a dependency on pressure and the mass flow rate of the injected water. The data are available for CFD-model development and validation.

Der vorliegende Bericht gibt einen zusammenfassenden Überblick der im Vorhaben erreichten Ergebnisse. Ziel war die Qualifikation von CFD-Methoden für Zweiphasenströmungen mit Phasenüber¬gang. Dafür werden neuartige experimentelle Daten benötigt. Diese können an der TOPFLOW-Anlage des HZDR generiert werden, da die Anlage Experimente in für die Reaktorsicher-heits¬forschung relevanten Skalen und Parametern mit innovativen Messtechniken verbindet. Die experimentellen Arbeiten umfassen Untersuchungen zu Strömungen in vertikalen Rohren mit Hilfe der ultraschnellen Röntgentomographie, zu Strömungen mit und ohne Phasenübergang in einem Testbassin sowie zur Gegenstrombegrenzung in einem Heißstrangmodell. Diese werden im vorliegenden Bericht nur kurz dargestellt, da es zu allen 3 Versuchsserien ausführliche Dokumentationen in separaten Berichten gibt. Ein wichtiges Ergebnis der Arbeiten zur CFD-Qualifizierung ist der Erstellung des Baseline-Modellkonzepts sowie die Erstellung des Baseline-Modells für polydisperse Blasenströmungen. Damit wird ein wesentlicher Beitrag zur Erhöhung der Vorhersagefähigkeit von CFD-Codes auf Basis des Zwei- oder Mehr-Fluid-Modells erreicht. Das innovative Generalized Two-Phase Flow Konzept (GENTOP) zielt hingegen auf eine Erweiterung der Einsatzmöglichkeiten der Zweiphasen-CFD. In vielen Strömungen treten unterschiedlicher Morphologien der Phasen bzw. Strömungsformen parallel in einer Strömungsdomäne auf. Außerdem gibt es Übergänge zwischen diesen Morphologien. Mit dem GENTOP-Konzept wurde erstmals ein Rahmen geschaffen der die Simulation solcher Strömungen auf konsistente Art und Weise ermöglicht. Spezielle Modellentwicklungen erfolgten mit dem Ziel einer besseren Modellierung des Phasenübergangs.

<div>Many two-phase flow systems might undergo flow instabilities even if the system is adiabatic but operates near the saturation conditions, especially in vertical flow conditions. Such instabilities are caused by <i>flashing</i> of the fluid in flow. Flashing is a sudden phase change in the fluid caused when local saturation enthalpy falls below the fluid enthalpy and the excess energy is used as latent heat for gas generation.</div><div> In the current analysis, a mathematical model is presented for analysis of such instability analytically. The conservation equations have been obtained by statistical averaging in time and space. Then, the concerned system is divided into various regions based on flow conditions, and these averaged equations are used to describe the flow. For flashing-based instability, two parameters are derived from constitutive relationships for the fluid. These two parameters are <i>Flashing Boundary</i> and <i>Gas Generation due to Flashing</i>. These parameters provide for the closure of the mathematical model. Some simple models for flashing have been developed and discussed.</div><div> The mathematical model is then solved analytically for <i>Uniform Heat</i> and <i>Flat Model</i> for the heater and flashing region respectively. The solution is in terms of the characteristic equation which is used to predict the onset of instability caused by flashing. The results are then plotted on the Subcooling-Phase Change number plane. It is observed that inlet and outlet restrictions in the flow does <b>not</b> affect the onset of flashing induced instability as the flow rate is coupled with the pressure drop of the system. This is important as these restrictions play a major role in other two-phase flow instabilities such as <i>Density Wave Oscillations</i></div><div> Finally, the stability boundary in the stability plane is compared to experimental data present for flashing. The comparison was made with data of S. Shi, A. Dixit, and F. Inada. The stability boundary satisfactorily agrees with the experimental data thus corroborating the present mathematical model and analysis.</div>

<div>With the rapid development of the advanced two-phase flow experimental technologies, more experimental databases with extended measurement ranges have been established to support the two-phase flow model development. The advantage of the Two Fluid model in modeling the complex two-phase flow phenomena over the mixture models stands out. One key aspect in the Two Fluid model development is the accurate modeling of the interfacial area between phases, which is strongly related to the interfacial mass, momentum, and energy transfer. As a closure relation of interfacial area concentration (interfacial area per unit volume) for the Two Fluid model, the Interfacial Area Transport Equation (IATE) provides dynamic predictions on the interfacial area change. It substantially solves the shortcoming of using flow-regime-dependent empirical correlations that can introduce numerical discontinuities between flow regimes. </div><div><br></div><div>The IATE has been extensively developed over the past twenty-five years. Many studies targeted on improving its prediction capability by developing bubble interaction source terms based on their experimental data. </div><div>The existing models are usually based on medium and large flow channels, yet the models may not be physically fit the small flow channels. The major reason is that the wall effect can have a larger influence on the two-phase flow in a small flow channel, as the surface area to volume ratio greatly increases. Therefore, the primary objectives of this study are to physically investigate the wall effect on two-phase flow and develop a generalized IATE by extending the application range of existing IATE from large and medium flow channels to small flow channel.</div><div><br></div><div>To achieve the objective, this study established a rigorous database of air-water two-phase flows in a small diameter pipe with its inner diameter of 12.7 mm, focusing on the bubbly-to-slug transition regime. The experimental analysis was performed on the pipe wall effect on the interfacial characteristics, based on the current experimental database and the existing experimental database collected on vertical pipes of different sizes. It is observed that 1) the pipe wall effect can alter the non-uniform radial two-phase distribution; 2) the bubbly-to-slug flow regime transition in a small diameter pipe happens in a smaller void fraction than in a large diameter pipe; 3) the bubble coalescence phenomenon can be more dominant for small pipe flow, and an intensive intergroup transfer can happen for the two-group interfacial area transport in two-phase flows. </div><div>As the interfacial area transport is directly related to the two-phase geometrical configuration, the two-phase geometrical parameters, void fraction and relative bubble size, are identified as the key parameters for modeling.</div><div><br></div><div>In the modeling of IATE source terms, the high geometrical scalability of the model is realized by properly including the wall effect into the modeling consideration. The following major improvements on the existing models are: 1) the inertia subrange assumption on the turbulent-driven interaction is properly improved; 2) the bubble-induced turbulent-driven interactions such as wake entrainment is revised by considering the wall effect on the wake region. In summary, models of bubble interaction due to random collision, wake entrainment, turbulent impact, and shearing-off are revised based on the existing studies on the IATE source terms development. The newly proposed interfacial area transport models are evaluated against an experimental database with 112 test conditions in total from a wide range of experimental pipe diameters from 12.7 mm to 304.8 mm. The new models can accurately capture the drastic intergroup transfer of void fraction and interfacial area concentration between two groups in transition flows. Overall, the relative error of void fraction and interfacial area concentration comparing with the experimental data are within ±15\% and ±10\%, respectively.</div>

Actuellement, les codes de calcul de composants thermohydrauliques de réacteurs nucléaires du CEA et d'EDF utilisent des modèles physiques diphasiques de mélange à 3 ou 4 équations. Or, il existe un fort besoin industriel pour des modèles physiques plus sophistiqués tels que le modèle diphasique à 6 équations voire des modèles multichamps. Par ailleurs, le code système CATHARE du CEA, un des codes systèmes les plus utilisés aujourd'hui sur le plan international, utilise un modèle physique diphasique à 6 équations et un schéma numérique semi-implicite de type ICE à maillages décalés. Le schéma de CATHARE est connu pour sa robustesse dans une large gamme de configurations d'écoulement. En s'inspirant de l'expérience de CATHARE, on propose de mettre en œuvre un schéma volumes finis colocalisés de type « pressure based ». Le but est d'obtenir un décentrement des flux qui assure la robustesse du schéma tout en gardant une bonne précision. De plus, le fait de pouvoir utiliser des maillages colocalisés (structurés ou non-structurés) permet de traiter de différentes configurations complexes de cœur de réacteurs et de réaliser des calculs fin d'inter assemblage. Le schéma doit conserver exactement la masse et l'énergie et la solution numérique doit converger lorsque l'on raffine le maillage. Le schéma doit être capable de traiter des cas d'apparition et de disparition des phases, par exemple le cas de la colonne bouillante où il y a changement de phase dû au transfert de chaleur, des cas de tuyères avec changement de phase dû à un élargissement ou un rétrécissement brusque, ou de séparation de phase par gravité. En outre le schéma doit être capable de traiter des configurations de calcul à faible nombre de Mach, par exemple le cas du renoyage d'un cœur de réacteur. L'objet de la thèse consistera à développer une méthode Volumes Finis co-localisés (dans l'esprit de Ghidaglia et al.) et la direction d'investigation s'inspirera des travaux de Jeong et al. qui a conduit au code CUPID. Références. Ghidaglia, J. M., Kumbaro, A., & Le Coq, G. (2001). On the numerical solution to two fluid models via a cell centered finite volume method. European Journal of Mechanics-B/Fluids, 20(6), 841-867. Jeong, J. J., Yoon, H. Y., Cho, H. K., Kim, J., & Park, I. K. (2008). A semi-implicit numerical scheme for a transient two-fluid three-field model on an unstructured grid. International Communications in Heat and Mass Transfer, 35(5), 597-605. / Currently, there is a strong industrial need for sophisticated physical models such as the two-phase model with 6 equations or multi-field models for the thermo-hydraulic calculation of nuclear reactor components. In addition, the system code CATHARE of CEA, one of the most utilized system codes at the international level, employs a two-phase 6 equations and a semi-implicit numerical scheme of ICE type on staggered grids. The code CATHARE is known for its robustness in a wide range of flow configurations. Drawing on the experience of CATHARE, we propose to implement a collocated finite volume pressure based scheme. The aim is to obtain a decentering of flux which ensures the robustness of the scheme while keeping good accuracy. In addition, being able to use collocated grids (structured or unstructured) can handle different complex configurations. The scheme should conserve exactly the mass and energy and the numerical solution needs to converge when the mesh is refined. The scheme should be able to handle cases of phase appearance and disappearance, for example in the case of boiling column where phase change is due to heat transfer, the case of nozzles with phase change due to a widening or abrupt narrowing, or phase separation by gravity. In addition the scheme should be capable of calculation configurations at low Mach number. The purpose is to develop a co-located Finite Volume method (in the spirit of Ghidaglia et al.) And direction of investigation is build on work of Jeong et al. which led to CUPID code. References: Ghidaglia, J. M., Kumbaro, A., & Le Coq, G. (2001). On the numerical solution to two fluid models via a cell centered finite volume method. European Journal of Mechanics-B/Fluids, 20(6), 841-867. Jeong, J. J., Yoon, H. Y., Cho, H. K., Kim, J., & Park, I. K. (2008). A semi-implicit numerical scheme for a transient two-fluid three-field model on an unstructured grid. International Communications in Heat and Mass Transfer, 35(5), 597-605.

Pour satisfaire aux dernières réglementations en matière d'émissions, des progrès importants sont encore attendus des moteurs à combustion interne. De plus, améliorer l'efficacité du moteur pour réduire les émissions et la consommation de carburant est devenu plus essentiel qu'auparavant. Mais, de nombreux phénomènes complexes restent mal compris dans ce domaine, tels que le processus d'injection de carburant. Nombreux logiciels pour la dynamique des fluides numérique (CFD) prenant en compte le changement de phase (comme la cavitation) et la modélisation de l’injection ont été développés et utilisés avec succès dans le processus d’injection. Néanmoins, il existe peu de codes CFD capables de simuler avec précision des conditions d’injection transcritiques, à partir d'une condition de température de carburant sous-critique vers un mélange supercritique dans la chambre de combustion. En effet, la plupart des modèles existants peuvent simuler des écoulements à phase unique, éventuellement dans des conditions supercritiques, ou des écoulements diphasiques dans des conditions sous-critiques. Par conséquent, il manque un modèle complet capable de traiter les conditions transcritiques, y compris la transition de phase possible entre les régimes souscritiques et supercritiques, ou entre les écoulements monophasiques et diphasiques, de manière dynamique. Cette thèse a pour objectif de relever ce défi.Pour cela, des modèles d'écoulement diphasique compressible de fluide réel basés sur une approche eulérienne-eulérienne avec prise en compte de l'équilibre de phase ont été développés et discutés dans le présent travail. Plus précisément, un modèle à 6-équation entièrement compressibles incluant les équations de bilan des phases liquide et gazeuse résolues séparément ; et un modèle à 4-équation qui résout les équations des bilans liquide et gazeux en équilibre mécanique et thermique sont proposés dans ce manuscrit. L’équation d’état Peng-Robinson EoS est sélectionné pour fermer les deux systèmes et pour faire face aux éventuels changements de phase et à la transition ou à la séparation des phases. En particulier, un solveur d'équilibre de phase a été développé et validé. Ensuite, une série de tests académiques 1D portant sur les phénomènes d'évaporation et de condensation effectués dans des conditions sous-critiques et supercritiques a été simulée et comparée aux données de la littérature et aux résultats académiques disponibles. Ensuite, les modèles d'écoulement en deux phases entièrement compressibles (systèmes à 6-équation et à 4- équation) ont été utilisés pour simuler les phénomènes de cavitation dans une buse 3D de taille réelle afin d'étudier l'effet de l’azote dissous sur la création et le développement de la cavitation. Le bon accord avec les données expérimentales prouve que le solveur proposé est capable de gérer le comportement complexe du changement de phase dans des conditions sous-critiques. Enfin, la capacité du solveur à traiter l’injection transcritique à des pressions et températures élevées a été validée par la modélisation réussie de l’injecteur Spray A du réseau de combustion moteur (ECN). / To satisfy latest stringent emission regulations, important progress is still be expected from internal combustion engines. In addition, improving engine efficiency to reduce the emission and fuel consumption has become more essential than before. But many complex phenomena remain poorly understood in this field, such as the fuel injection process. Numerous software programs for computational fluid dynamics (CFD) considering phase change (such as cavitation) and injection modelling, have been developed and used successfully in the injection process. Nevertheless, there are few CFD codes able to simulate correctly transcritical conditions starting from a subcritical fuel temperature condition towards a supercritical mixture in the combustion chamber. Indeed, most of the existing models can simulate either single-phase flows possibly in supercritical condition or two-phase flows in subcritical condition; lacking therefore, a comprehensive model which can deal with transcritical condition including possible phase transition from subcritical to supercritical regimes, or from single-phase to two-phase flows, dynamically. This thesis aims at dealing with this challenge. For that, real fluid compressible two-phase flow models based on Eulerian-Eulerian approach with the consideration of phase equilibrium have been developed and discussed in the present work. More precisely, a fully compressible 6-equation model including liquid and gas phases balance equations solved separately; and a 4-equation model which solves the liquid and gas balance equations in mechanical and thermal equilibrium, are proposed in this manuscript. The Peng-Robinson equation of state (EoS) is selected to close both systems and to deal with the eventual phase change or phase transition. Particularly, a phase equilibrium solver has been developed and validated. Then, a series of 1D academic tests involving the evaporation and condensation phenomena performed under subcritical and supercritical conditions have been simulated and compared with available literature data and analytical results. Then the fully compressible two-phase flow models (6-Equation and 4-Equation systems) have been employed to simulate the cavitation phenomena in a real size 3D nozzle to investigate the effect of dissolved N2 on the inception and developing of cavitation. The good agreement with experimental data proves the solver can handle the complex phase change behavior in subcritical condition. Finally, the capability of the solver in dealing with the transcritical injection at high pressure and temperature conditions has been further validated through the successful modelling of the engine combustion network (ECN) Spray A injector.

Thesis (Ph.D.)--Massachusetts Institute of Technology, Dept. of Nuclear Engineering, 1980. / MICROFICHE COPY AVAILABLE IN ARCHIVES AND SCIENCE / Includes bibliographical references. / by John Edward Kelly. / Ph.D.

京都大学 / 0048 / 新制・課程博士 / 博士(エネルギー科学) / 甲第21887号 / エネ博第388号 / 新制||エネ||75(附属図書館) / 京都大学大学院エネルギー科学研究科エネルギー基礎科学専攻 / (主査)教授 齊藤 泰司, 教授 横峯 健彦, 准教授 伊藤 啓 / 学位規則第4条第1項該当 / Doctor of Energy Science / Kyoto University / DFAM

<p>Two-phase jet impingement combines the latent heat absorbed by boiling heat transfer with the strong forced convection of an impinging jet. It is a compact and highly effective heat transfer method that is capable of high heat transfer coefficients and high boiling critical heat flux limits. This makes it a suitable technology for electronics immersion cooling applications when configured as a submerged jet of a dielectric coolant. Previous studies have focused on the heat-flux-controlled nucleate boiling performance of an impingement jet up to the critical heat flux. Exploration of other boiling regimes that occur under temperature-controlled surfaces is of fundamental importance to fully understand the design space. It has been shown for free jets that a high and consistent heat flux can be dissipated over a wide range of surface superheats in the transition boiling regime when the surface is temperature-controlled. This effect is strongest in the stagnation zone, directly beneath the jet. Literature that studies this so-called “shoulder effect”, or heat flux shoulder, is scarce and almost completely focused on applications in metals processing using free jets of water as the coolant. It has been hypothesized previously that the impinging subcooled liquid delays and disrupts the start of film boiling, thereby dissipating heat flux levels comparable to that during nucleate boiling. To exploit operation in this unique transition boiling regime for potential applications in immersion cooling of electronics, the occurrence of this shoulder effect, as well as means for estimating the shoulder heat flux across different operating conditions, must be investigated for submerged jets and dielectric coolants. </p> <p>In this work, temperature-controlled submerged jet impingement is experimentally characterized using HFE-7100. A copper heater sized to be completely covered by the jet stagnation zone is increased in surface temperature throughout the transition boiling regime via a PID controller, which allows for steady-state temperature-controlled data to be acquired in this regime. The boiling curves, including critical heat flux and shoulder heat flux, are measured for jet velocities from 0.5-3 m/s and inlet subcooling from 5-30 K. The shoulder effect is shown to exist in these conditions. High-speed imaging is used to relate the flow behavior to the boiling thermal measurements and shows that the shoulder heat flux effect is an enhanced film heat transfer in the film-like mode of transition boiling. Trends and dependencies on inlet subcooling and jet velocity are measured and used to assess available predictive tools. It is observed that there is a proportionality between the critical heat flux and the shoulder heat flux. This implies a mechanistic similarity between the two effects. With further data to correlate, this similarity can potentially be used to predict the shoulder heat flux leveraging existing correlations for the critical heat flux, widening the design space of two-phase jet impingement systems. </p>