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.

Panel Method Formulation for Oscillating Airfoils in Sonic Flow

The mathematical model for two-dimensional unsteady sonic flow, based on the classical diffusion equation with imaginary coefficient, is presented and discussed. The main purpose is to develop a rigorous formulation in order to bring into light the correspondence between the sonic, supersonic and subsonic panel method theory. Source and doublet integrals are obtained and Laplace transformation demonstrates that, in fact, the source integral is the solution of the doublet integral equation. It is shown that the doublet-only formulation reduces to a Volterra integral equation of the first kind and a numerical method is proposed in order to solve it. To the authors' knowledge this is the first reported solution to the unsteady sonic thin airfoil problem through the use of doublet singularities. Comparisons with the source-only formulation are shown for the problem of a flat plate in combined harmonic heaving and pitching motion.

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.

Analyzing hydraulic head of fluid flow in underground porous medium

Hydraulic head is distributed through a medium with porous aspect. The analysis of hydraulic head from one point to another is used by the Richard's equation. This equation is equivalent to the groundwater ow equation that predicts the volumetric water contents. COMSOL 3.5 is used for computation applying Richard's equation. A rectangle of 100 meters of length and 10 meters of large (depth) with 0,1 m/s fl ux of inlet as source of our fl uid is simulated. The domain have Richards' equation model in two dimension (2D). Hydraulic head increases proportional with moisture content.

Modeling of Liquid-Solid Flow in Industrial Scale

In the present work, liquid-solid flow in industrial scale is modeled using the commercial software of Computational Fluid Dynamics (CFD) ANSYS Fluent 14.5. In literature, there are few studies on liquid-solid flow in industrial scale, but any information about the particular case with modified geometry cannot be found. The aim of this thesis is to describe the strengths and weaknesses of the multiphase models, when a large-scale application is studied within liquid-solid flow, including the boundary-layer characteristics. The results indicate that the selection of the most appropriate multiphase model depends on the flow regime. Thus, careful estimations of the flow regime are recommended to be done before modeling. The computational tool is developed for this purpose during this thesis. The homogeneous multiphase model is valid only for homogeneous suspension, the discrete phase model (DPM) is recommended for homogeneous and heterogeneous suspension where pipe Froude number is greater than 1.0, while the mixture and Eulerian models are able to predict also flow regimes, where pipe Froude number is smaller than 1.0 and particles tend to settle. With increasing material density ratio and decreasing pipe Froude number, the Eulerian model gives the most accurate results, because it does not include simplifications in Navier-Stokes equations like the other models. In addition, the results indicate that the potential location of erosion in the pipe depends on material density ratio. Possible sedimentation of particles can cause erosion and increase pressure drop as well. In the pipe bend, especially secondary flows, perpendicular to the main flow, affect the location of erosion.

Streamside management zone (SMZ) efficiency in herbicide retention from simulated surface flow

Plot-scale overland flow experiments were conducted to evaluate the efficiency of streamside management zones (SMZs) in retaining herbicides in runoff generated from silvicultural activities. Herbicide retention was evaluated for five different slopes (2, 5, 10, 15, and 20%), two cover conditions (undisturbed O horizon and raked surface), and two periods under contrasting soil moisture conditions (summer dry and winter wet season) and correlated to O horizon and site conditions. Picloram (highly soluble in water) and atrazine (moderately sorbed into soil particles) at concentrations in the range of 55 and 35 µg L-1 and kaolin clay (approximately 5 g L-1) were mixed with 13.000 liters of water and dispersed over the top of 5 x 10 m forested plots. Surface flow was collected 2, 4, 6, and 10 m below the disperser to evaluate the changes in concentration as it moved through the O horizon and surface soil horizon-mixing zone. Results showed that, on average, a 10 m long forested SMZ removed around 25% of the initial concentration of atrazine and was generally ineffective in reducing the more soluble picloram. Retention of picloram was only 6% of the applied quantity. Percentages of mass reduction by infiltration were 36% for atrazine and 20% for picloram. Stronger relationships existed between O horizon depth and atrazine retention than in any other measured variable, suggesting that better solid-solution contact associated with flow through deeper O horizons is more important than either velocity or soil moisture as a determinant of sorption.

The left ventricular contractility of the rat heart is modulated by changes in flow and a1-adrenoceptor stimulation

Myocardial contractility depends on several mechanisms such as coronary perfusion pressure (CPP) and flow as well as on a1-adrenoceptor stimulation. Both effects occur during the sympathetic stimulation mediated by norepinephrine. Norepinephrine increases force development in the heart and produces vasoconstriction increasing arterial pressure and, in turn, CPP. The contribution of each of these factors to the increase in myocardial performance needs to be clarified. Thus, in the present study we used two protocols: in the first we measured mean arterial pressure, left ventricular pressure and rate of rise of left ventricular pressure development in anesthetized rats (N = 10) submitted to phenylephrine (PE) stimulation before and after propranolol plus atropine treatment. These observations showed that in vivo a1-adrenergic stimulation increases left ventricular-developed pressure (P<0.05) together with arterial blood pressure (P<0.05). In the second protocol, we measured left ventricular isovolumic systolic pressure (ISP) and CPP in Langendorff constant flow-perfused hearts. The hearts (N = 7) were perfused with increasing flow rates under control conditions and PE or PE + nitroprusside (NP). Both CPP and ISP increased (P<0.01) as a function of flow. CPP changes were not affected by drug treatment but ISP increased (P<0.01). The largest ISP increase was obtained with PE + NP treatment (P<0.01). The results suggest that both mechanisms, i.e., direct stimulation of myocardial a1-adrenoceptors and increased flow, increased cardiac performance acting simultaneously and synergistically.

A pulsatile flow model for in vitro quantitative evaluation of prosthetic valve regurgitation

A pulsatile pressure-flow model was developed for in vitro quantitative color Doppler flow mapping studies of valvular regurgitation. The flow through the system was generated by a piston which was driven by stepper motors controlled by a computer. The piston was connected to acrylic chambers designed to simulate "ventricular" and "atrial" heart chambers. Inside the "ventricular" chamber, a prosthetic heart valve was placed at the inflow connection with the "atrial" chamber while another prosthetic valve was positioned at the outflow connection with flexible tubes, elastic balloons and a reservoir arranged to mimic the peripheral circulation. The flow model was filled with a 0.25% corn starch/water suspension to improve Doppler imaging. A continuous flow pump transferred the liquid from the peripheral reservoir to another one connected to the "atrial" chamber. The dimensions of the flow model were designed to permit adequate imaging by Doppler echocardiography. Acoustic windows allowed placement of transducers distal and perpendicular to the valves, so that the ultrasound beam could be positioned parallel to the valvular flow. Strain-gauge and electromagnetic transducers were used for measurements of pressure and flow in different segments of the system. The flow model was also designed to fit different sizes and types of prosthetic valves. This pulsatile flow model was able to generate pressure and flow in the physiological human range, with independent adjustment of pulse duration and rate as well as of stroke volume. This model mimics flow profiles observed in patients with regurgitant prosthetic valves.

Pressurized hot water flow-through extraction of birch wood

Effective processes to fractionate the main compounds in biomass, such as wood, are a prerequisite for an effective biorefinery. Water is environmentally friendly and widely used in industry, which makes it a potential solvent also for forest biomass. At elevated temperatures over 100 °C, water can readily hydrolyse and dissolve hemicelluloses from biomass. In this work, birch sawdust was extracted using pressurized hot water (PHWE) flow-through systems. The hypothesis of the work was that it is possible to obtain polymeric, water-soluble hemicelluloses from birch sawdust using flow-through PHW extractions at both laboratory and large scale. Different extraction temperatures in the range 140–200 °C were evaluated to see the effect of temperature to the xylan yield. The yields and extracted hemicelluloses were analysed to obtain sugar ratios, the amount of acetyl groups, furfurals and the xylan yields. Higher extraction temperatures increased the xylan yield, but decreased the molar mass of the dissolved xylan. As the extraction temperature increased, more acetic acid was released from the hemicelluloses, thus further decreasing the pH of the extract. There were only trace amounts of furfurals present after the extractions, indicating that the treatment was mild enough not to degrade the sugars further. The sawdust extraction density was increased by packing more sawdust in the laboratory scale extraction vessel. The aim was to obtain extracts with higher concentration than in typical extraction densities. The extraction times and water flow rates were kept constant during these extractions. The higher sawdust packing degree decreased the water use in the extractions and the extracts had higher hemicellulose concentrations than extractions with lower sawdust degrees of packing. The molar masses of the hemicelluloses were similar in higher packing degrees and in the degrees of packing that were used in typical PHWE flow-through extractions. The structure of extracted sawdust was investigated using small angle-(SAXS) and wide angle (WAXS) x-ray scattering. The cell wall topography of birch sawdust and extracted sawdust was compared using x-ray tomography. The results showed that the structure of the cell walls of extracted birch sawdust was preserved but the cell walls were thinner after the extractions. Larger pores were opened inside the fibres and cellulose microfibrils were more tightly packed after the extraction. Acetate buffers were used to control the pH of the extracts during the extractions. The pH control prevented excessive xylan hydrolysis and increased the molar masses of the extracted xylans. The yields of buffered extractions were lower than for plain water extractions at 160–170 °C, but at 180 °C yields were similar to those from plain water and pH buffers. The pH can thus be controlled during extraction with acetate buffer to obtain xylan with higher molar mass than those obtainable using plain water. Birch sawdust was extracted both in the laboratory and pilot scale. The performance of the PHWE flow-through system was evaluated in the laboratory and the pilot scale using vessels with the same shape but different volumes, with the same relative water flow through the sawdust bed, and in the same extraction temperature. Pre-steaming improved the extraction efficiency and the water flow through the sawdust bed. The extracted birch sawdust and the extracted xylan were similar in both laboratory and pilot scale. The PHWE system was successfully scaled up by a factor of 6000 from the laboratory to pilot scale and extractions performed equally well in both scales. The results show that a flow-through system can be further scaled up and used to extract water-soluble xylans from birch sawdust. Extracted xylans can be concentrated, purified, and then used in e.g. films and barriers, or as building blocks for novel material applications.

Genesis and Evolution of Brittle Structures In Southwestern Finland and Western South Africa. Insights into Fault Reactivation, Fluid Flow and Structural Maturity in Precambrian Cratons

The bedrock of old crystalline cratons is characteristically saturated with brittle structures formed during successive superimposed episodes of deformation and under varying stress regimes. As a result, the crust effectively deforms through the reactivation of pre-existing structures rather than by through the activation, or generation, of new ones, and is said to be in a state of 'structural maturity'. By combining data from Olkiluoto Island, southwestern Finland, which has been investigated as the potential site of a deep geological repository for high-level nuclear waste, with observations from southern Sweden, it can be concluded that the southern part of the Svecofennian shield had already attained structural maturity during the Mesoproterozoic era. This indicates that the phase of activation of the crust, i.e. the time interval during which new fractures were generated, was brief in comparison to the subsequent reactivation phase. Structural maturity of the bedrock was also attained relatively rapidly in Namaqualand, western South Africa, after the formation of first brittle structures during Neoproterozoic time. Subsequent brittle deformation in Namaqualand was controlled by the reactivation of pre-existing strike-slip faults.In such settings, seismic events are likely to occur through reactivation of pre-existing zones that are favourably oriented with respect to prevailing stresses. In Namaqualand, this is shown for present day seismicity by slip tendency analysis, and at Olkiluoto, for a Neoproterozoic earthquake reactivating a Mesoproterozoic fault. By combining detailed field observations with the results of paleostress inversions and relative and absolute time constraints, seven distinctm superimposed paleostress regimes have been recognized in the Olkiluoto region. From oldest to youngest these are: (1) NW-SE to NNW-SSE transpression, which prevailed soon after 1.75 Ga, when the crust had sufficiently cooled down to allow brittle deformation to occur. During this phase conjugate NNW-SSE and NE-SW striking strike-slip faults were active simultaneous with reactivation of SE-dipping low-angle shear zones and foliation planes. This was followed by (2) N-S to NE-SW transpression, which caused partial reactivation of structures formed in the first event; (3) NW-SE extension during the Gothian orogeny and at the time of rapakivi magmatism and intrusion of diabase dikes; (4) NE-SW transtension that occurred between 1.60 and 1.30 Ga and which also formed the NW-SE-trending Satakunta graben located some 20 km north of Olkiluoto. Greisen-type veins also formed during this phase. (5) NE-SW compression that postdates both the formation of the 1.56 Ga rapakivi granites and 1.27 Ga olivine diabases of the region; (6) E-W transpression during the early stages of the Mesoproterozoic Sveconorwegian orogeny and which also predated (7) almost coaxial E-W extension attributed to the collapse of the Sveconorwegian orogeny. The kinematic analysis of fracture systems in crystalline bedrock also provides a robust framework for evaluating fluid-rock interaction in the brittle regime; this is essential in assessment of bedrock integrity for numerous geo-engineering applications, including groundwater management, transient or permanent CO2 storage and site investigations for permanent waste disposal. Investigations at Olkiluoto revealed that fluid flow along fractures is coupled with low normal tractions due to in-situ stresses and thus deviates from the generally accepted critically stressed fracture concept, where fluid flow is concentrated on fractures on the verge of failure. The difference is linked to the shallow conditions of Olkiluoto - due to the low differential stresses inherent at shallow depths, fracture activation and fluid flow is controlled by dilation due to low normal tractions. At deeper settings, however, fluid flow is controlled by fracture criticality caused by large differential stress, which drives shear deformation instead of dilation.

Development of project cash flow forecasting process

Process management refers to improving the key functions of a company. The main functions of the case company - project management, procurement, finance, and human resource - use their own separate systems. The case company is in the process of changing its software. Different functions will use the same system in the future. This software change causes changes in some of the company’s processes. Project cash flow forecasting process is one of the changing processes. Cash flow forecasting ensures the sufficiency of money and prepares for possible changes in the future. This will help to ensure the company’s viability. The purpose of the research is to describe a new project cash flow forecasting process. In addition, the aim is to analyze the impacts of the process change, with regard to the project control department’s workload and resources through the process measurement, and how the impacts take the department’s future operations into account. The research is based on process management. Processes, their descriptions, and the way the process management uses the information, are discussed in the theory part of this research. The theory part is based on literature and articles. Project cash flow and forecasting-related benefits are also discussed. After this, the project cash flow forecasting as-is and to-be processes are described by utilizing information, obtained from the theoretical part, as well as the know-how of the project control department’s personnel. Written descriptions and cross-functional flowcharts are used for descriptions. Process measurement is based on interviews with the personnel – mainly cost controllers and department managers. The process change and the integration of two processes will allow work time for other things, for example, analysis of costs. In addition to the quality of the cash flow information will improve compared to the as-is process. Analyzing the department’s other main processes, department’s roles, and their responsibilities should be checked and redesigned. This way, there will be an opportunity to achieve the best possible efficiency and cost savings.

Trypan blue exclusion assay by flow cytometry

Dye exclusion tests are used to determine the number of live and dead cells. These assays are based on the principle that intact plasma membranes in live cells exclude specific dyes, whereas dead cells do not. Although widely used, the trypan blue (TB) exclusion assay has limitations. The dye can be incorporated by live cells after a short exposure time, and personal reliability, related to the expertise of the analyst, can affect the results. We propose an alternative assay for evaluating cell viability that combines the TB exclusion test and the high sensitivity of the flow cytometry technique. Previous studies have demonstrated the ability of TB to emit fluorescence when complexed with proteins. According to our results, TB/bovine serum albumin and TB/cytoplasmic protein complexes emit fluorescence at 660 nm, which is detectable by flow cytometry using a 650-nm low-pass band filter. TB at 0.002% (w/v) was defined as the optimum concentration for distinguishing unstained living cells from fluorescent dead cells, and fluorescence emission was stable for 30 min after cell treatment. Although previous studies have shown that TB promotes green fluorescence quenching, TB at 0.002% did not interfere with green fluorescence in human live T-cells stained with anti-CD3/fluorescein isothiocyanate (FITC) monoclonal antibody. We observed a high correlation between the percentage of propidium iodide+CD3/FITC+ and TB+CD3/FITC+ cells, as well as similar double-stained cell profiles in flow cytometry dot-plot graphs. Taken together, the results indicate that a TB exclusion assay by flow cytometry can be employed as an alternative tool for quick and reliable cell viability analysis.

Flow in porous media of oil-water systems

Fluid flow behaviour in porous media is a conundrum. Therefore, this research is focused on filtration-volumetric characterisation of fractured-carbonate sediments, coupled with their proper simulation. For this reason, at laboratory rock properties such as pore volume, permeability and porosity are measured, later phase permeabilities and oil recovery in function of flow rate are assessed. Furthermore, the rheological properties of three oils are measured and analysed. Finally based on rock and fluid properties, a model using COMSOL Multiphysics is built in order to compare the experimental and simulated results. The rock analyses show linear relation between flow rate and differential pressure, from which phase permeabilities and pressure gradient are determined, eventually the oil recovery under low and high flow rate is established. In addition, the oils reveal thixotropic properties as well as non-Newtonian behaviour described by Bingham model, consequently Carreau viscosity model for the used oil is given. Given these points, the model for oil and water is built in COMSOL Multiphysics, whereupon successfully the reciprocity between experimental and simulated results is analysed and compared. Finally, a two-phase displacement model is elaborated.

The effects of plants(Typha Latifolia)and root-bed medium on the treatment of domestic sewage within a vertical flow constructed wetland /

The effect that plants {Typha latifolia) as well as root-bed medium physical and chemical characteristics have on the treatment of primary treated domestic wastewater within a vertical flow constructed wetland system was investigated. Five sets of cells, with two cells in each set, were used. Each cell was made of concrete and measured 1 .0 m X 1 .0 m and was 1.3 m deep. Four different root-bed media were tested : Queenston Shale, Fonthill Sand, Niagara Shale and a Michigan Sand. Four of the sets contained plants and a single type of root-bed medium. The influence of plants was tested by operating a Queenston Shale set without plants. Due to budget constraints no replicates were constructed. All of the sets were operated independently and identically for twenty-eight months. Twelve months of data are presented here, collected after 16 months of continuous operation. Root-bed medium type did not influence BOD5 removal. All of the sets consistently met Ontario Ministry of Environment (MOE) requirements (<25 mg/L) for BOD5 throughout the year. The 12 month average BOD5 concentration from all sets with plants was below 2.36 mg/L. All of the sets were within MOE discharge requirements (< 25 mg/L) for suspended solids with set effluent concentrations ranging from 1.53 to 14.80 mg/L. The Queenston Shale and Fonthill Sand media removed the most suspended solids while the Niagara Shale set produced suspended solids. The set containing Fonthill Sand was the only series to meet MOE discharge requirements (< Img/L) for total phosphorus year-round with a twelve month mean effluent concentration of 0.23 mg/L. Year-round all of the root-bed media were well below MOE discharge requirements (< 20mg/L in winter and < 10 mg/L in sumnner) for ammonium. The Queenston Shale and Fonthill Sand sets removed the most total nitrogen. Plants had no effect on total nitrogen removal, but did influence how nitrogen was cycled within the system. Plants increased the removal of suspended solids by 14%, BOD5 by 10% and total phosphorus by 22%. Plants also increased the amount of dissolved oxygen that entered the system. During the plant growing season removal of total phosphorus was better in all sets with plants regardless of media type. The sets containing Queenston Shale and Fonthill Sand media achieved the best results and plants in the Queenston Shale set increased treatment efficiency for every parameter except nitrogen. Vertical flow wetland sewage treatment systems can be designed and built to consistently meet MOE discharge requirements year-round for BOD5, suspended solids, total phosphorus and ammonium. This system Is generally superior to the free water systems and sub-surface horizontal flow systems in cold climate situations.

Investigations in flow injection analysis (FIA): conventional FIA determination of chloride in water; sequential injection (SI)-FIA determination of mercury in water

Flow injection analysis (FIA) was applied to the determination of both chloride ion and mercury in water. Conventional FIA was employed for the chloride study. Investigations of the Fe3 +/Hg(SCN)2/CI-,450 nm spectrophotometric system for chloride determination led to the discovery of an absorbance in the 250-260 nm region when Hg(SCN)2 and CI- are combined in solution, in the absence of iron(III). Employing an in-house FIA system, absorbance observed at 254 nm exhibited a linear relation from essentially 0 - 2000 Jlg ml- 1 injected chloride. This linear range spanning three orders of magnitude is superior to the Fe3+/Hg(SCN)2/CI- system currently employed by laboratories worldwide. The detection limit obtainable with the proposed method was determin~d to be 0.16 Jlg ml- 1 and the relative standard deviation was determined to be 3.5 % over the concentration range of 0-200 Jig ml- 1. Other halogen ions were found to interfere with chloride determination at 254 nm whereas cations did not interfere. This system was successfully applied to the determination of chloride ion in laboratory water. Sequential injection (SI)-FIA was employed for mercury determination in water with the PSA Galahad mercury amalgamation, and Merlin mercury fluorescence detection systems. Initial mercury in air determinations involved injections of mercury saturated air directly into the Galahad whereas mercury in water determinations involved solution delivery via peristaltic pump to a gas/liquid separator, after reduction by stannous chloride. A series of changes were made to the internal hardware and valving systems of the Galahad mercury preconcentrator. Sequential injection solution delivery replaced the continuous peristaltic pump system and computer control was implemented to control and integrate all aspects of solution delivery, sample preconcentration and signal processing. Detection limits currently obtainable with this system are 0.1 ng ml-1 HgO.