The project was aimed at describing the mixing phenomena relevant for both safety analysis, particularly in steam line break and boron dilution scenarios, and mixing phenomena of interest for economical operation and the structural integrity. Measurement data from a set of mixing experiments, gained by using advanced measurement techniques with enhanced resolution in time and space help to improve the basic understanding of turbulent mixing and to provide data for Computational Fluid Dynamics (CFD) code validation. Slug mixing tests simulating the start-up of the first main circulation pump are performed with two 1:5 scaled facilities: The Rossendorf coolant mixing model ROCOM and the VATTENFALL test facility, modelling a German Konvoi type and a Westinghouse type three-loop PWR, respectively. Additional data on slug mixing in a VVER-1000 type reactor gained at a 1:5 scaled metal mock-up at EDO Gidropress are provided. Experimental results on mixing of fluids with density differences obtained at ROCOM and the FORTUM PTS test facility are made available. Concerning mixing phenomena of interest for operational issues and thermal fatigue, flow distribution data available from commissioning tests (Sizewell-B for PWRs, Loviisa and Paks for VVERs) are used together with the data from the ROCOM facility as a basis for the flow distribution studies. The test matrix on flow distribution and steady state mixing performed at ROCOM comprises experiments with various combinations of running pumps and various mass flow rates in the working loops. Computational fluid dynamics calculations are accomplished for selected experiments with two different CFD codes (CFX-5, FLUENT). Best practice guidelines (BPG) are applied in all CFD work when choosing computational grid, time step, turbulence models, modelling of internal geometry, boundary conditions, numerical schemes and convergence criteria. The BPG contain a set of systematic procedures for quantifying and reducing numerical errors. The knowledge of these numerical errors is a prerequisite for the proper judgement of model errors. The strategy of code validation based on the BPG and a matrix of CFD code validation calculations have been elaborated. Besides of the benchmark cases, additional experiments were calculated by new partners and observers, joining the project later. Based on the "best practice solutions", conclusions on the applicability of CFD for turbulent mixing problems in PWR were drawn and recommendations on CFD modelling were given. The high importance of proper grid generation was outlined. In general, second order discretization schemes should be used to minimise numerical diffusion. First order schemes can provide physically wrong results. With optimised "production meshes" reasonable results were obtained, but due to the complex geometry of the flow domains, no fully grid independent solutions were achieved. Therefore, with respect to turbulence models, no final conclusions can be given. However, first order turbulence models like K-e or SST K-w are suitable for momentum driven slug mixing. For buoyancy driven mixing (PTS scenarios), Reynolds stress models provide better results.

The project was aimed at describing the mixing phenomena relevant for both safety analysis, particularly in steam line break and boron dilution scenarios, and mixing phenomena of interest for economical operation and the structural integrity. Measurement data from a set of mixing experiments, gained by using advanced measurement techniques with enhanced resolution in time and space help to improve the basic understanding of turbulent mixing and to provide data for Computational Fluid Dynamics (CFD) code validation. Slug mixing tests simulating the start-up of the first main circulation pump are performed with two 1:5 scaled facilities: The Rossendorf coolant mixing model ROCOM and the VATTENFALL test facility, modelling a German Konvoi type and a Westinghouse type three-loop PWR, respectively. Additional data on slug mixing in a VVER-1000 type reactor gained at a 1:5 scaled metal mock-up at EDO Gidropress are provided. Experimental results on mixing of fluids with density differences obtained at ROCOM and the FORTUM PTS test facility are made available. Concerning mixing phenomena of interest for operational issues and thermal fatigue, flow distribution data available from commissioning tests (Sizewell-B for PWRs, Loviisa and Paks for VVERs) are used together with the data from the ROCOM facility as a basis for the flow distribution studies. The test matrix on flow distribution and steady state mixing performed at ROCOM comprises experiments with various combinations of running pumps and various mass flow rates in the working loops. Computational fluid dynamics calculations are accomplished for selected experiments with two different CFD codes (CFX-5, FLUENT). Best practice guidelines (BPG) are applied in all CFD work when choosing computational grid, time step, turbulence models, modelling of internal geometry, boundary conditions, numerical schemes and convergence criteria. The BPG contain a set of systematic procedures for quantifying and reducing numerical errors. The knowledge of these numerical errors is a prerequisite for the proper judgement of model errors. The strategy of code validation based on the BPG and a matrix of CFD code validation calculations have been elaborated. Besides of the benchmark cases, additional experiments were calculated by new partners and observers, joining the project later. Based on the "best practice solutions", conclusions on the applicability of CFD for turbulent mixing problems in PWR were drawn and recommendations on CFD modelling were given. The high importance of proper grid generation was outlined. In general, second order discretization schemes should be used to minimise numerical diffusion. First order schemes can provide physically wrong results. With optimised "production meshes" reasonable results were obtained, but due to the complex geometry of the flow domains, no fully grid independent solutions were achieved. Therefore, with respect to turbulence models, no final conclusions can be given. However, first order turbulence models like K-e or SST K-w are suitable for momentum driven slug mixing. For buoyancy driven mixing (PTS scenarios), Reynolds stress models provide better results.

The flow distribution in the primary circuit of the pressurized water reactor was studied with experiments and Computational Fluid Dynamics (CFD) simulations. The main focus was on the flow field and mixing in the downcomer of the pressure vessel: how the different factors like the orientation of operating loops, the total loop flow rate and the asymmetry of the loop flow rates affect the outcome. In addition to the flow field studies the overall applicability of CFD methods for primary circuit thermal-hydraulic analysis was evaluated based on the CFD simulations of the mixing experiments of the ROCOM (Rossendorf Coolant Mixing Model) test facility and the mixing experiments of the Paks NPP. The experimental part of the work in work package 3 included series of steady state mixing experiments with the ROCOM test facility and the publication of results of Paks VVER-440 NPP thermal mixing experiments. The ROCOM test facility models a 4-loop KONVOI type reactor. In the steady-state mixing experiments the velocity field in the downcomer was measured using laser Doppler anemometry and the concentration of the tracer solution fed from one loop was measured at the downcomer and at the core inlet plane. The varied parameters were the number and orientation of the operating loops, the total flow rate and the (asymmetric) flow rate of individual loops. The Paks NPP thermal mixing experiments took place during commissioning tests of replaced steam generator safety valves in 1987-1989. It was assumed that in the reactor vessels of Paks VVER-440 NPP equipped with six loops the mixing of the coolant is not ideal. For the realistic determination of the active core inlet temperature field for the transients and accidents associated with different level temperature asymmetry a set of mixing factors were determined. Based on data from the online core monitoring system and a separate mathematical model the mixing factors for loop flows at the core inlet were determined. In the numerical simulation part of the work package 3 the detailed measurements of ROCOM tests were used for the validation of CFD methods for primary circuit studies. The selected steady state mixing experiments were simulated with CFD codes CFX-4, CFX-5 and FLUENT. The velocity field in the downcomer and the mixing of the scalar were compared between CFD simulations and experiments. The CFD simulations of full scale PWR included the simulation of Paks VVER-440 mixing experiment and the simulation of Loviisa VVER-440 downcomer flow field. In the simulations of Paks experiments the experimental and simulated concentration field at the core inlet were compared and conclusions made concerning the results overall and the VVER-440 specific geometry modelling aspects like how to model the perforated elliptic bottom plate and what is the effect of the cold leg bends to the flow field entering to the downcomer. With Loviisa simulations the qualitative comparison was made against the original commissioning experiments but the emphasis was on the CFD method validation and testing. The overall conclusion concerning the CFD modelling of the flow field and mixing in the PWR primary circuit could be that the current computation capacity and physical models also in commercial codes is beginning to be sufficient for simulations giving reliable and useful results for many real primary circuit applications. However the misuse of CFD methods is easy, and the general as well as the nuclear power specific modelling guidelines should be followed when the CFD simulations are made.

The goal of the work described in this report was the experimental investigation of the mixing of coolant with different quality (temperature, boron concentration) in nuclear reactors on the way from the cold leg through the downcomer and lower plenum to the core inlet in a systematic way. The obtained data were used for the clarification of the mixing mechanisms and form a data basis for the validation of computational fluid dynamics (CFD) codes. For these purposes, experiments on slug mixing have been performed at two test facilities, modelling different reactor types in scale 1:5, the Rossendorf and Vattenfall test facilities. The corresponding accident scenario is the start-up of first main coolant pump (MCP) after formation of a slug of lower borated water during the reflux-condenser mode phase of a small break loss of coolant accident (LOCA). The matrices for the experiments were elaborated on the basis of the key phenomena, being responsible for the coolant mixing during pump start-up. Slug mixing tests have also been performed at the VVER-1000 facility of EDO Gidropress to meet the specifics of this reactor type. The mixing of slugs of water of different quality is also very important for pre-stressed thermal shock (PTS) situations. In emergency core cooling (ECC) situations after a LOCA, cold ECC water is injected into the hot water in the cold leg and downcomer. Due to the large temperature differences, thermal shocks are induced at the reactor pressure vessel (RPV) wall. Temperature distributions near the wall and temperature gradients in time are important to be known for the assessment of thermal stresses. One of the important phenomena in connection with PTS is thermal stratification, a flow condition with a vertical temperature profile in a horizontal pipe. Due to the fluctuating character of the flow, this may cause thermal fatigue in the pipe. Besides of thermal fatigue, a single thermal shock can also be relevant for structural integrity, if it is large enough, especially in the case, that the brittle fracture temperature of the RPV material is reduced due to radiation embrittlement. Therefore, additional to the investigations of slug mixing during re-start of coolant circulation, the mixing of slugs or streams of water with higher density with the ambient fluid in the RPV was investigated. The aim of these investigations was to study the process of turbulent mixing under the influence of buoyancy forces caused by the temperature differences. Heat transfer to the wall and thermal conductivity in the wall material have not been considered. Experiments on density driven mixing were carried out at the Rossendorf and the Fortum PTS facilities.

The flow distribution in the primary circuit of the pressurized water reactor was studied with experiments and Computational Fluid Dynamics (CFD) simulations. The main focus was on the flow field and mixing in the downcomer of the pressure vessel: how the different factors like the orientation of operating loops, the total loop flow rate and the asymmetry of the loop flow rates affect the outcome. In addition to the flow field studies the overall applicability of CFD methods for primary circuit thermal-hydraulic analysis was evaluated based on the CFD simulations of the mixing experiments of the ROCOM (Rossendorf Coolant Mixing Model) test facility and the mixing experiments of the Paks NPP. The experimental part of the work in work package 3 included series of steady state mixing experiments with the ROCOM test facility and the publication of results of Paks VVER-440 NPP thermal mixing experiments. The ROCOM test facility models a 4-loop KONVOI type reactor. In the steady-state mixing experiments the velocity field in the downcomer was measured using laser Doppler anemometry and the concentration of the tracer solution fed from one loop was measured at the downcomer and at the core inlet plane. The varied parameters were the number and orientation of the operating loops, the total flow rate and the (asymmetric) flow rate of individual loops. The Paks NPP thermal mixing experiments took place during commissioning tests of replaced steam generator safety valves in 1987-1989. It was assumed that in the reactor vessels of Paks VVER-440 NPP equipped with six loops the mixing of the coolant is not ideal. For the realistic determination of the active core inlet temperature field for the transients and accidents associated with different level temperature asymmetry a set of mixing factors were determined. Based on data from the online core monitoring system and a separate mathematical model the mixing factors for loop flows at the core inlet were determined. In the numerical simulation part of the work package 3 the detailed measurements of ROCOM tests were used for the validation of CFD methods for primary circuit studies. The selected steady state mixing experiments were simulated with CFD codes CFX-4, CFX-5 and FLUENT. The velocity field in the downcomer and the mixing of the scalar were compared between CFD simulations and experiments. The CFD simulations of full scale PWR included the simulation of Paks VVER-440 mixing experiment and the simulation of Loviisa VVER-440 downcomer flow field. In the simulations of Paks experiments the experimental and simulated concentration field at the core inlet were compared and conclusions made concerning the results overall and the VVER-440 specific geometry modelling aspects like how to model the perforated elliptic bottom plate and what is the effect of the cold leg bends to the flow field entering to the downcomer. With Loviisa simulations the qualitative comparison was made against the original commissioning experiments but the emphasis was on the CFD method validation and testing. The overall conclusion concerning the CFD modelling of the flow field and mixing in the PWR primary circuit could be that the current computation capacity and physical models also in commercial codes is beginning to be sufficient for simulations giving reliable and useful results for many real primary circuit applications. However the misuse of CFD methods is easy, and the general as well as the nuclear power specific modelling guidelines should be followed when the CFD simulations are made.

The goal of the work described in this report was the experimental investigation of the mixing of coolant with different quality (temperature, boron concentration) in nuclear reactors on the way from the cold leg through the downcomer and lower plenum to the core inlet in a systematic way. The obtained data were used for the clarification of the mixing mechanisms and form a data basis for the validation of computational fluid dynamics (CFD) codes. For these purposes, experiments on slug mixing have been performed at two test facilities, modelling different reactor types in scale 1:5, the Rossendorf and Vattenfall test facilities. The corresponding accident scenario is the start-up of first main coolant pump (MCP) after formation of a slug of lower borated water during the reflux-condenser mode phase of a small break loss of coolant accident (LOCA). The matrices for the experiments were elaborated on the basis of the key phenomena, being responsible for the coolant mixing during pump start-up. Slug mixing tests have also been performed at the VVER-1000 facility of EDO Gidropress to meet the specifics of this reactor type. The mixing of slugs of water of different quality is also very important for pre-stressed thermal shock (PTS) situations. In emergency core cooling (ECC) situations after a LOCA, cold ECC water is injected into the hot water in the cold leg and downcomer. Due to the large temperature differences, thermal shocks are induced at the reactor pressure vessel (RPV) wall. Temperature distributions near the wall and temperature gradients in time are important to be known for the assessment of thermal stresses. One of the important phenomena in connection with PTS is thermal stratification, a flow condition with a vertical temperature profile in a horizontal pipe. Due to the fluctuating character of the flow, this may cause thermal fatigue in the pipe. Besides of thermal fatigue, a single thermal shock can also be relevant for structural integrity, if it is large enough, especially in the case, that the brittle fracture temperature of the RPV material is reduced due to radiation embrittlement. Therefore, additional to the investigations of slug mixing during re-start of coolant circulation, the mixing of slugs or streams of water with higher density with the ambient fluid in the RPV was investigated. The aim of these investigations was to study the process of turbulent mixing under the influence of buoyancy forces caused by the temperature differences. Heat transfer to the wall and thermal conductivity in the wall material have not been considered. Experiments on density driven mixing were carried out at the Rossendorf and the Fortum PTS facilities.

Celui-ci a pour objectif de valider l'ensemble des modèles physiques utilisés dans un code de simulation numérique pour simuler un écoulement de type annulaire dispersé en conduite rencontré lors d'une procédure de lavage à eau utilisé dans les raffineries. Pour ce faire une banque de données expérimentale est mise en place sur des configurations représentatives de celles utilisées en condition industrielle. La géométrie retenue comporte une zone horizontale d'injection rectiligne avec un injecteur central, suivi d'un coude à 90° situé dans un plan vertical. Différentes conditions expérimentales permettent d'étudier l'influence de la vitesse du gaz, de la condition d'injection du brouillard et de la pression sur les différents processus physiques. Ces résultats comprenant des visualisations du brouillard et du film pariétale, des mesures de taille et de distribution de gouttes,des mesures de débit et d'épaisseur de film, sont analysés pour faire ressortir les principaux mécanismes d'interaction entre le gaz et la phase dispersée, le gaz et le film liquide pariétal et la phase dispersée et le film pariétal. En parallèle, des premières simulations, avec une approche RANS, sont réalisées avec le code CEDRE de l'ONERA et les résultats sont confrontés aux mesures. / This work has been realised within a CIFRE contract with TOTAL. Its aim was to validate all the physical models used in a computation, which simulates an annular dispersed flow through a pipe used in a water washing process in refinery plants. That is why, a whole set of data has been gathered using experimental boundary conditions which are representative to those used in industrial configurations. The geometry is made of a horizontal pipe with a centred nozzle followed by a 90º elbow in the vertical plane. Several experimental boundary conditions enable one to study the influence of the gas velocity, the type of the spray injection and the pressure on the different physical phenomena. These results including spray and liquid film visualisations, droplets distribution and size measurements as well as liquid film thickness and mass flow measurements were analysed in order to extract the main interaction mechanism between the gas and the dispersed phase, the gas and the liquid film, and the dispersed phase and the annular liquid film. Meanwhile, simulations using a RANS approach were realized with the ONERA code named CEDRE and its results were compared to the gathered measurements.