CFD Modelling of Direct Contact Condensation in Suppression Pools by Applying Condensation Models of Separated Flow

The condensation rate has to be high in the safety pressure suppression pool systems of Boiling Water Reactors (BWR) in order to fulfill their safety function. The phenomena due to such a high direct contact condensation (DCC) rate turn out to be very challenging to be analysed either with experiments or numerical simulations. In this thesis, the suppression pool experiments carried out in the POOLEX facility of Lappeenranta University of Technology were simulated. Two different condensation modes were modelled by using the 2-phase CFD codes NEPTUNE CFD and TransAT. The DCC models applied were the typical ones to be used for separated flows in channels, and their applicability to the rapidly condensing flow in the condensation pool context had not been tested earlier. A low Reynolds number case was the first to be simulated. The POOLEX experiment STB-31 was operated near the conditions between the ’quasi-steady oscillatory interface condensation’ mode and the ’condensation within the blowdown pipe’ mode. The condensation models of Lakehal et al. and Coste & Lavi´eville predicted the condensation rate quite accurately, while the other tested ones overestimated it. It was possible to get the direct phase change solution to settle near to the measured values, but a very high resolution of calculation grid was needed. Secondly, a high Reynolds number case corresponding to the ’chugging’ mode was simulated. The POOLEX experiment STB-28 was chosen, because various standard and highspeed video samples of bubbles were recorded during it. In order to extract numerical information from the video material, a pattern recognition procedure was programmed. The bubble size distributions and the frequencies of chugging were calculated with this procedure. With the statistical data of the bubble sizes and temporal data of the bubble/jet appearance, it was possible to compare the condensation rates between the experiment and the CFD simulations. In the chugging simulations, a spherically curvilinear calculation grid at the blowdown pipe exit improved the convergence and decreased the required cell count. The compressible flow solver with complete steam-tables was beneficial for the numerical success of the simulations. The Hughes-Duffey model and, to some extent, the Coste & Lavi´eville model produced realistic chugging behavior. The initial level of the steam/water interface was an important factor to determine the initiation of the chugging. If the interface was initialized with a water level high enough inside the blowdown pipe, the vigorous penetration of a water plug into the pool created a turbulent wake which invoked the chugging that was self-sustaining. A 3D simulation with a suitable DCC model produced qualitatively very realistic shapes of the chugging bubbles and jets. The comparative FFT analysis of the bubble size data and the pool bottom pressure data gave useful information to distinguish the eigenmodes of chugging, bubbling, and pool structure oscillations.

Pressure drop coefficients for elliptic and circular sections in one, two and three-row arrangements of plate fin and tube heat exchangers

The objective of the present work is the experimental determination of pressure drop coefficients (loss coefficients) for elliptic and circular sections in one, two and three-row arrangements of plate fin and tube heat exchangers. The experiments permitted to correlate the dimensionless loss coefficient with the flow Reynolds number in the rectangular channel formed by the plate fins. The experimental technique consisted of the measurement of the longitudinal pressure distribution along the flow channel, for several values of air mass flow rate. The total number of data runs, each one characterized by the flow Reynolds number, was 216. The present geometry is used in compact heat exchangers for air conditioning systems, heaters, radiators, and others. Also, it is verified the influence of the utilization of elliptic tubes, instead of circular ones, in the pressure drop. The measurements were performed for Reynolds numbers ranging from 200 to 1900.

Effect of different surface roughnesses on a turbulent boundary layer

The classical treatment of rough wall turbulent boundary layers consists in determining the effect the roughness has on the mean velocity profile. This effect is usually described in terms of the roughness function delta U+. The general implication is that different roughness geometries with the same delta U+ will have similar turbulence characteristics, at least at a sufficient distance from the roughness elements. Measurements over two different surface geometries (a mesh roughness and spanwise circular rods regularly spaced in the streamwise direction) with nominally the same delta U+ indicate significant differences in the Reynolds stresses, especially those involving the wall-normal velocity fluctuation, over the outer region. The differences are such that the Reynolds stress anisotropy is smaller over the mesh roughness than the rod roughness. The Reynolds stress anisotropy is largest for a smooth wall. The small-scale anisotropy and interniittency exhibit much smaller differences when the Taylor microscale Reynolds number and the Kolmogorov-normalized mean shear are nominally the same. There is nonetheless evidence that the small-scale structure over the three-dimensional mesh roughness conforms more closely with isotropy than that over the rod-roughened and smooth walls.

Turbulent heat transfer and pressure drop in pinned annular regions

Experiments were performed to determine average heat transfer coefficients and friction factors for turbulent flow through annular ducts with pin fins. The measurements were carried out by means of a double-pipe heat exchanger. The total number of pins attached to the inner wall of the annular region was 560. The working fluids were air, flowing in the annular channel, and water through the inner circular tube. The average heat transfer coefficients of the pinned air-side were obtained from the experimental determination of the overall heat transfer coefficients of the heat exchanger and from the knowledge of the average heat transfer coefficients of the circular pipe (water-side), which could be found in the pertinent literature. To attain fully developed conditions, the heat exchanger was built with additional lengths before and after the test section. The inner circular duct of the heat exchanger and the pin fins were made of brass. Due to the high thermal conductivity of the brass, the small tube thickness and water temperature variation, the surface of the internal tube was practically isothermal. The external tube was made of an industrial plastic which was insulated from the environment by means of a glass wool batt. In this manner, the outer surface of the annular channel can be considered adiabatic. The results are presented in dimensionless forms, in terms of average Nusselt numbers and friction factors as functions of the flow Reynolds number, ranging from 13,000 to 80,000. The pin fin efficiency, which depends on the heat transfer coefficient, is also determined as a function of dimensionless parameters. A comparison of the present results with those for smooth sections (without pins) is also presented. The purpose of such a comparison is to study the influence of the presence of the pins on the pressure drop and heat transfer rate.

Static pressure and wall shear stress distributions in air flow in a seven wire-wrapped rod bundle

An experimental investigation is performed in a turbulent flow in a seven wire-wrapped rod bundle, mounted in an open air facility. Static pressure distributions are measured on central and peripheral rods. By using a Preston tube, the wall shear stress profiles are experimentally obtained along the perimeter of the rods. The geometric parameters of the test section are P/D=1.20 and H/D=15. The measuring section is located at L/D=40 from the air inlet. It is observed that the dimensionless static pressure and wall shear stress profiles are nearly independent of the Reynolds number and strongly dependent of the wire-spacer position, with abrupt variations of the parameters in the neighborhood of the wires.

Prediction of transients and control reactions in a transonic wind tunnel

This work describes a lumped parameter mathematical model for the prediction of transients in an aerodynamic circuit of a transonic wind tunnel. Control actions to properly handle those perturbations are also assessed. The tunnel circuit technology is up to date and incorporates a novel feature: high-enthalpy air injection to extend the tunnel’s Reynolds number capability. The model solves the equations of continuity, energy and momentum and defines density, internal energy and mass flow as the basic parameters in the aerodynamic study as well as Mach number, stagnation pressure and stagnation temperature, all referred to test section conditions, as the main control variables. The tunnel circuit response to control actions and the stability of the flow are numerically investigated. Initially, for validation purposes, the code was applied to the AWT ("Altitude Wind Tunnel" of NASA-Lewis). In the sequel, the Brazilian transonic wind tunnel was investigated, with all the main control systems modeled, including injection.

Direct Simulation Monte Carlo Method Applied to Aerothermodynamics

An axisymmetric supersonic flow of rarefied gas past a finite cylinder was calculated applying the direct simulation Monte Carlo method. The drag force, the coefficients of pressure, of skin friction, and of heat transfer, the fields of density, of temperature, and of velocity were calculated as function of the Reynolds number for a fixed Mach number. The variation of the Reynolds number is related to the variation of the Knudsen number, which characterizes the gas rarefaction. The present results show that all quantities in the transition regime (Knudsen number is about the unity) are significantly different from those in the hydrodynamic regime, when the Knudsen number is small.

Three-dimensional supersonic flow over a spike-nosed body of revolution

The unsteady, viscous, supersonic flow over a spike-nosed body of revolution is numerically investigated by solving the Navier-Stokes equations. The time-accurate computations are performed employing an implicit algorithm based on the second-order time-accurate LU-SGS scheme with the incorporation of a subiteration procedure to maintain time accuracy. The characteristics of the flow field for a Mach number of 3.0, Reynolds number of 7.87 x 10(6)/m, and angles of attack of 5 and 10 degrees are described. Self-sustained asymmetric shock wave oscillations were observed in the numerical computations for these angles of attack. The main characteristic of the flow field, as well as its influence on drag coefficient is discussed.

Large-eddy simulation of wind flows over complex terrains for wind energy applications

Wind energy has obtained outstanding expectations due to risks of global warming and nuclear energy production plant accidents. Nowadays, wind farms are often constructed in areas of complex terrain. A potential wind farm location must have the site thoroughly surveyed and the wind climatology analyzed before installing any hardware. Therefore, modeling of Atmospheric Boundary Layer (ABL) flows over complex terrains containing, e.g. hills, forest, and lakes is of great interest in wind energy applications, as it can help in locating and optimizing the wind farms. Numerical modeling of wind flows using Computational Fluid Dynamics (CFD) has become a popular technique during the last few decades. Due to the inherent flow variability and large-scale unsteadiness typical in ABL flows in general and especially over complex terrains, the flow can be difficult to be predicted accurately enough by using the Reynolds-Averaged Navier-Stokes equations (RANS). Large- Eddy Simulation (LES) resolves the largest and thus most important turbulent eddies and models only the small-scale motions which are more universal than the large eddies and thus easier to model. Therefore, LES is expected to be more suitable for this kind of simulations although it is computationally more expensive than the RANS approach. With the fast development of computers and open-source CFD software during the recent years, the application of LES toward atmospheric flow is becoming increasingly common nowadays. The aim of the work is to simulate atmospheric flows over realistic and complex terrains by means of LES. Evaluation of potential in-land wind park locations will be the main application for these simulations. Development of the LES methodology to simulate the atmospheric flows over realistic terrains is reported in the thesis. The work also aims at validating the LES methodology at a real scale. In the thesis, LES are carried out for flow problems ranging from basic channel flows to real atmospheric flows over one of the most recent real-life complex terrain problems, the Bolund hill. All the simulations reported in the thesis are carried out using a new OpenFOAM® -based LES solver. The solver uses the 4th order time-accurate Runge-Kutta scheme and a fractional step method. Moreover, development of the LES methodology includes special attention to two boundary conditions: the upstream (inflow) and wall boundary conditions. The upstream boundary condition is generated by using the so-called recycling technique, in which the instantaneous flow properties are sampled on aplane downstream of the inlet and mapped back to the inlet at each time step. This technique develops the upstream boundary-layer flow together with the inflow turbulence without using any precursor simulation and thus within a single computational domain. The roughness of the terrain surface is modeled by implementing a new wall function into OpenFOAM® during the thesis work. Both, the recycling method and the newly implemented wall function, are validated for the channel flows at relatively high Reynolds number before applying them to the atmospheric flow applications. After validating the LES model over simple flows, the simulations are carried out for atmospheric boundary-layer flows over two types of hills: first, two-dimensional wind-tunnel hill profiles and second, the Bolund hill located in Roskilde Fjord, Denmark. For the twodimensional wind-tunnel hills, the study focuses on the overall flow behavior as a function of the hill slope. Moreover, the simulations are repeated using another wall function suitable for smooth surfaces, which already existed in OpenFOAM® , in order to study the sensitivity of the flow to the surface roughness in ABL flows. The simulated results obtained using the two wall functions are compared against the wind-tunnel measurements. It is shown that LES using the implemented wall function produces overall satisfactory results on the turbulent flow over the two-dimensional hills. The prediction of the flow separation and reattachment-length for the steeper hill is closer to the measurements than the other numerical studies reported in the past for the same hill geometry. The field measurement campaign performed over the Bolund hill provides the most recent field-experiment dataset for the mean flow and the turbulence properties. A number of research groups have simulated the wind flows over the Bolund hill. Due to the challenging features of the hill such as the almost vertical hill slope, it is considered as an ideal experimental test case for validating micro-scale CFD models for wind energy applications. In this work, the simulated results obtained for two wind directions are compared against the field measurements. It is shown that the present LES can reproduce the complex turbulent wind flow structures over a complicated terrain such as the Bolund hill. Especially, the present LES results show the best prediction of the turbulent kinetic energy with an average error of 24.1%, which is a 43% smaller than any other model results reported in the past for the Bolund case. Finally, the validated LES methodology is demonstrated to simulate the wind flow over the existing Muukko wind farm located in South-Eastern Finland. The simulation is carried out only for one wind direction and the results on the instantaneous and time-averaged wind speeds are briefly reported. The demonstration case is followed by discussions on the practical aspects of LES for the wind resource assessment over a realistic inland wind farm.

MA311exam 1998

Exam questions and solutions in LaTex

MA311exam 2000 pdf

Exam questions and solutions in PDF

MA311exam 1999

Exam questions and solutions in LaTex

MA311exam 1999 pdf

Exam questions and solutions in PDF

MA311exam 1998 pdf

Exam questions and solutions in PDF

MA311exam 2000

Exam questions and solutions in LaTex