Tectonic and metamorphic evolution of the Central Himalayan Domain in Southeast Zanskar (Kashmir, India)

La région du Zanskar, étudiée dans le cadre de ce travail, se situe au passage entre deux domaines himalayens fortement contrastés, la Séquence Cristalline du Haut Himalaya (HHCS), composée de roches métamorphiques et l'Himalaya Tethysien (TH), composé de séries sédimentaires. La transition entre ces deux domaines est marquée par une structure tectonique majeure, la Zone de Cisaillement du Zanskar (ZSZ), au sein de laquelle on observe une augmentation extrêmement rapide, mais néanmoins graduelle, du degré du métamorphisme entre le TH et le HHCS. Il a été établi que le HHCS n'est autre que l'équivalent métamorphique des séries sédimentaires de la base du TH. C'est principalement lors d'un épisode de mise en place de nappes à vergence sudouest, entre l'Eocène moyen et l'Oligocène, que les séries sédimentaires de la base du TH ont été entraînées en profondeur où elles ont subi un métamorphisme de type barrovien. Au début du Miocène, le HHCS à été exhumé en direction du sud-ouest sous forme d'une grande nappe, délimitée a sa base par le MCT (principal chevauchement central) et à son sommet par la Zone de Cisaillement du Zanskar. L'ensemble des zones barroviennes, de la zone à biotite jusqu'à la zone à disthène, a été cisaillée par les mouvements en faille normale au sommet du HHCS et se retrouve actuellement sur une épaisseur d'environ 1 kilomètre au sein de la ZSZ. La décompression associée à l'exhumation du HHCS a provoqué la fusion partielle d'une partie du HHCS et a donné naissance à des magmas de composition leucogranitiques. Grâce à la géothermobarometrie, et connaissant la géométrie de la ZSZ, il nous a été possible de déterminer que le rejet le long de cette structure d'extension est d'au moins 35?9 kilomètres. Une série d'arguments nous permet cependant de suggérer que ce rejet aurait pu être encore bien plus important (~100km). Les données géochronologiques nous permettent de contraindre la durée des mouvements d'extension le long de la ZSZ à 2.4?0.2 Ma entre 22.2?0.2 Ma et 19.8?0.1 Ma. Ce travail apporte de nouvelles données sur les processus métamorphiques, magmatiques et tectoniques liés aux phénomènes d'extension syn-orogeniques.<br/><br/>The southeastern part of Zanskar is located at the transition between two major Himalayan domains of contrasting metamorphic grade, the High Himalayan Crystalline Sequence (HHCS) and the Tethyan Himalaya (TH). The transition between the TH and the HHCS is marked by a very rapid, although perfectly gradual, decrease in metamorphic grade, which coincides with a major tectonic structure, the Zanskar Shear Zone (ZSZ). It is now an established fact that the relation between the HHCS and the TH is not one of basement-cover type, but that the metasedimentary series of the HHCS represent the metamorphic equivalent of the lowermost sedimentary series of the TH. This transformation of sedimentary series into metamorphic rocks, and hence the differentiation between the TH and the HHCS, is the consequence of crustal thickening associated to the formation of large scale southwest vergent nappes within the Tethyan Himalaya sedimentary series. This, Middle Eocene to Oligocene, episode of crustal thickening and associated Barrovian metamorphism is followed, shortly after, by the exhumation of the HHCS as a, large scale, south-west vergent, nappe. Foreword The exhumation of the HHCS nappe is marked by the activation of two contemporaneous structures, the Main Central Thrust at its base and the Zanskar Shear Zone at its top. Extensional movements along the ZSZ, caused the Barrovian biotite to the kyanite zones to be sheared and constricted within the ~1 km thick shear zone. Decompression associated with the exhumation of the HHCS induced the formation of leucogranitic magmas through vapour-absent partial melting of the highest-grade rocks. The combination of geothermobarometric data with a geometric model of the ZSZ allowed us to constrain the net slip at the top of the HHCS to be at least 35?9 kilometres. A set of arguments however suggests that these movements might have been much more important (~ 100 km). Geochronological data coupled with structural observations constrain the duration of ductile shearing along the ZSZ to 2.4?0.2 Ma between 22.2?0.2 Ma and 19.8?0.1 Ma. This study also addresses the consequences of synorogenic extension on the metamorphic, tectonic and magmatic evolution of the upper parts of the High Himalayan Crystalline Sequence.

Electromagnetic Design of a Solid Steel Rotor Motor for Demanding Operation Environments

The solid-rotor induction motor provides a mechanically and thermally reliable solution for demanding environments where other rotor solutions are prohibited or questionable. Solid rotors, which are manufactured of single pieces of ferromagnetic material, are commonly used in motors in which the rotationspeeds exceed substantially the conventional speeds of laminated rotors with squirrel-cage. During the operation of a solid-rotor electrical machine, the rotor core forms a conductor for both the magnetic flux and the electrical current. This causes an increase in the rotor resistance and rotor leakage inductance, which essentially decreases the power factor and the efficiency of the machine. The electromagnetic problems related to the solid-rotor induction motor are mostly associated with the low performance of the rotor. Therefore, the main emphasis in this thesis is put on the solid steel rotor designs. The rotor designs studied in thisthesis are based on the fact that the rotor construction should be extremely robust and reliable to withstand the high mechanical stresses caused by the rotational velocity of the rotor. In addition, the demanding operation environment sets requirements for the applied materials because of the high temperatures and oxidizing acids, which may be present in the cooling fluid. Therefore, the solid rotors analyzed in this thesis are made of a single piece of ferromagnetic material without any additional parts, such as copper end-rings or a squirrel-cage. A pure solid rotor construction is rigid and able to keep its balance over a large speed range. It also may tolerate other environmental stresses such as corroding substances or abrasive particles. In this thesis, the main target is to improve the performance of an induction motor equipped with a solid steel rotor by traditional methods: by axial slitting of the rotor, by selecting a proper rotor core material and by coating the rotor with a high-resistive stainless ferromagnetic material. In the solid steel rotor calculation, the rotor end-effects have a significant effect on the rotor characteristics. Thus, the emphasis is also put on the comparison of different rotor endfactors. In addition, a corrective slip-dependent end-factor is proposed. The rotor designs covered in this thesis are the smooth solid rotor, the axially slitted solid rotor and the slitted rotor having a uniform ferromagnetic coating cylinder. The thesis aims at design rules for multi-megawatt machines. Typically, mega-watt-size solidrotor machines find their applications mainly in the field of electric-motor-gas-compression systems, in steam-turbine applications, and in various types of largepower pump applications, where high operational speeds are required. In this thesis, a 120 kW, 10 000 rpm solid-rotor induction motor is usedas a small-scale model for such megawatt-range solid-rotor machines. The performance of the 120 kW solid-rotor induction motors is determined by experimental measurements and finite element calculations.

Mass transfer and particleinteractions in granular systems and suspension flows

The purpose of this study was to investigate some important features of granular flows and suspension flows by computational simulation methods. Granular materials have been considered as an independent state ofmatter because of their complex behaviors. They sometimes behave like a solid, sometimes like a fluid, and sometimes can contain both phases in equilibrium. The computer simulation of dense shear granular flows of monodisperse, spherical particles shows that the collisional model of contacts yields the coexistence of solid and fluid phases while the frictional model represents a uniform flow of fluid phase. However, a comparison between the stress signals from the simulations and experiments revealed that the collisional model would result a proper match with the experimental evidences. Although the effect of gravity is found to beimportant in sedimentation of solid part, the stick-slip behavior associated with the collisional model looks more similar to that of experiments. The mathematical formulations based on the kinetic theory have been derived for the moderatesolid volume fractions with the assumption of the homogeneity of flow. In orderto make some simulations which can provide such an ideal flow, the simulation of unbounded granular shear flows was performed. Therefore, the homogeneous flow properties could be achieved in the moderate solid volume fractions. A new algorithm, namely the nonequilibrium approach was introduced to show the features of self-diffusion in the granular flows. Using this algorithm a one way flow can beextracted from the entire flow, which not only provides a straightforward calculation of self-diffusion coefficient but also can qualitatively determine the deviation of self-diffusion from the linear law at some regions nearby the wall inbounded flows. Anyhow, the average lateral self-diffusion coefficient, which was calculated by the aforementioned method, showed a desirable agreement with thepredictions of kinetic theory formulation. In the continuation of computer simulation of shear granular flows, some numerical and theoretical investigations were carried out on mass transfer and particle interactions in particulate flows. In this context, the boundary element method and its combination with the spectral method using the special capabilities of wavelets have been introduced as theefficient numerical methods to solve the governing equations of mass transfer in particulate flows. A theoretical formulation of fluid dispersivity in suspension flows revealed that the fluid dispersivity depends upon the fluid properties and particle parameters as well as the fluid-particle and particle-particle interactions.

Quaternary active tectonic structures in the offshore Bajo Segura basin (SE Iberian Peninsula Mediterranean Sea).

The Bajo Segura fault zone (BSFZ) is the northern terminal splay of the Eastern Betic shear zone (EBSZ), a large left-lateral strike-slip fault system of sigmoid geometry stretching more than 450 km from Alicante to Almería. The BSFZ extends from the onshore Bajo Segura basin further into the Mediterranean Sea and shows a moderate instrumental seismic activity characterized by small earthquakes. Nevertheless, the zone was affected by large historical earthquakes of which the largest was the 1829 Torrevieja earthquake (IEMS98 X). The onshore area of the BSFZ is marked by active transpressive structures (faults and folds), whereas the offshore area has been scarcely explored from the tectonic point of view. During the EVENT-SHELF cruise, a total of 10 high-resolution single-channel seismic sparker profiles were obtained along and across the offshore Bajo Segura basin. Analysis of these profiles resulted in (a) the identification of 6 Quaternary seismo-stratigraphic units bounded by five horizons corresponding to regional erosional surfaces related to global sea level lowstands; and (b) the mapping of the active sub-seafloor structures and their correlation with those described onshore. Moreover, the results suggest that the Bajo Segura blind thrust fault or the Torrevieja left-lateral strike-slip fault, with prolongation offshore, could be considered as the source of the 1829 Torrevieja earthquake. These data improve our understanding of present deformation along the BSFZ and provide new insights into the seismic hazard in the area.

The coal bearing Cenozoic As Pontes Basin (NW Spain): Geological influence on its coal's characteristics

Lignite deposits in the Cenozoic As Pontes strike-slip basin (northwestern Spain)were formed as a function of specific paleoclimatic conditions and tectonic evolution of the basin.

Growth, production and quality of pineapple in response to herbicide use

In pineapple fields, weed competition is exacerbated by the fact that the crop is small and has a very slow vegetative development. The objective of this study was to determine the effects of herbicides on growth, yield and quality of pineapple, cultivar 'Pérola'. The experimental design was in randomized blocks with four treatments and four replications. Treatments consisted of weeding by hoe and the herbicides diuron; fluazifop-p-butyl and atrazine + S-metolachlor applied in post-emergence. The characteristics evaluated monthly during the vegetative stage were stem diameter, D-leaf length, number of leaves, number of emitted leaves and percentage of natural floral induction. In the reproductive phase, evaluations were made of average fruit weight (g) with and without crown, fruits length and diameter, number of slip, slip-sucker and sucker type seedlings, determination of soluble solids and pH in the pulp. There was no effect of herbicide treatment on the vegetative growth characteristics. Stem diameter increased until 330 days after planting, showing a decrease after this period. The D-leaf grew over time in all treatments, although phytotoxicity symptoms were observed after the first application of herbicides. The traits evaluated on the reproductive phase showed no significant differences in response to treatments. Therefore, the use of diuron fluazifop-p-butyl and atrazine + S-metolachlor did not affect growth, yield and fruit quality of pineapple, cultivar 'Pérola'.

Growth and production of irrigated vitória pineapple grown in semi-arid conditions

This study aimed to evaluate the growth characteristics of irrigated Vitória pineapple plants grown in semi-arid conditions and determine its developmental stages based on those characteristics. It was used a randomized block design with four replicates. The experimental treatments were: plant harvest at 270, 330, 390, 450, 510, 570, 690, 750, and 810 days after planting (DAP). The following variables were determined: plant height, stem diameter, D-leaf length, D-leaf fresh and dry mass, biomass production of plants and plant parts (organs), and vegetative biomass. Five phenological stages are proposed based on vegetative biomass production: < 20% biomass production (V1); 21-40% (V2); 41-60% (V3); 61-80% (V4); and > 80% (V5). The maximum growth rate for plant height, D-leaf length, and stem diameter was observed at the end of the phenological stage V1 (390-411 DAP), and at the end of stage V5 these plant traits had average values of 106, 82, and 7 cm, respectively. The maximum biomass accumulation rates were observed at stages V4 and V5, resulting in a final fruit yield and total fresh biomass of 72 t ha-1 and 326 t ha-1, respectively. Finally, we estimated that 80% of the accumulated biomass may remain in the field after fruit and slip harvest, and could be incorporated as plant residue into the soil.

The Sottunga-Jurmo shear zone : structure and deformation history of a crustal-scale ductile shear zone in SW Finland

The aim of this study is to gain a better understanding of the structure and the deformation history of a NW-SE trending regional, crustal-scale shear structure in the Åland archipelago, SW Finland, called the Sottunga-Jurmo shear zone (SJSZ). Approaches involving e.g. structural geology, geochronology, geochemistry and metamorphic petrology were utilised in order to reconstruct the overall deformation history of the study area. The study therefore describes several features of the shear zone including structures, kinematics and lithologies within the study area, the ages of the different deformation phases (ductile to brittle) within the shear zone, as well as some geothermobarometric results. The results indicate that the SJSZ outlines a major crustal discontinuity between the extensively migmatized rocks NE of the shear zone and the unmigmatised, amphibolite facies rocks SW of the zone. The main SJSZ shows overall dextral lateral kinematics with a SW-side up vertical component and deformation partitioning into pure shear and simple shear dominated deformation styles that was intensified toward later stages of the deformation history. The deformation partitioning resulted in complex folding and refolding against the SW margin of the SJSZ, including conical and sheath folds, and in a formation of several minor strike-slip shear zones both parallel and conjugate to the main SJSZ in order to accommodate the regional transpressive stresses. Different deformation phases within the study area were dated by SIMS (zircon U-Pb), ID-TIMS (titanite U-Pb) and 40Ar/39Ar (pseudotachylyte wholerock) methods. The first deformation phase within the ca. 1.88 Ga rocks of the study area is dated at ca. 1.85 Ga, and the shear zone was reactivated twice within the ductile regime (at ca. 1.83 Ga and 1.79 Ga), during which the strain was successively increasingly partitioned into the main SJSZ and the minor shear zones. The age determinations suggest that the orogenic processes within the study area did not occur in a temporal continuum; instead, the metamorphic zircon rims and titanites show distinct, 10-20 Ma long breaks in deformation between phases of active deformation. The results of this study further imply slow cooling of the rocks through 600-700ºC so that at 1.79 Ga, 2 the temperature was still at least 600ºC. The highest recorded metamorphic pressures are 6.4-7.1 kbar. At the late stages or soon after the last ductile phase (ca. 1.79 Ga), relatively high-T mylonites and ultramylonites were formed, witnessing extreme deformation partitioning and high strain rates. After the rocks reached lower amphibolite facies to amphibolite-greenschist facies transitional conditions (ca. 500-550ºC), they cooled rapidly, probably due to crustal uplift and exhumation. The shear zone was reactivated at least once within the semi-brittle to brittle regime between ca. 1.79 Ga and 1.58 Ga, as evidenced by cataclasites and pseudotachylytes. In summary, the results of this study suggest that the Sottunga-Jurmo shear zone (and the South Finland shear zone) defines a major crustal discontinuity, and played a central role in accommodating the regional stresses during and after the Svecofennian orogeny.

High-Speed Solid-Rotor Induction Machine – Electromagnetic Calculation and Design

Within the latest decade high-speed motor technology has been increasingly commonly applied within the range of medium and large power. More particularly, applications like such involved with gas movement and compression seem to be the most important area in which high-speed machines are used. In manufacturing the induction motor rotor core of one single piece of steel it is possible to achieve an extremely rigid rotor construction for the high-speed motor. In a mechanical sense, the solid rotor may be the best possible rotor construction. Unfortunately, the electromagnetic properties of a solid rotor are poorer than the properties of the traditional laminated rotor of an induction motor. This thesis analyses methods for improving the electromagnetic properties of a solid-rotor induction machine. The slip of the solid rotor is reduced notably if the solid rotor is axially slitted. The slitting patterns of the solid rotor are examined. It is shown how the slitting parameters affect the produced torque. Methods for decreasing the harmonic eddy currents on the surface of the rotor are also examined. The motivation for this is to improve the efficiency of the motor to reach the efficiency standard of a laminated rotor induction motor. To carry out these research tasks the finite element analysis is used. An analytical calculation of solid rotors based on the multi-layer transfer-matrix method is developed especially for the calculation of axially slitted solid rotors equipped with wellconducting end rings. The calculation results are verified by using the finite element analysis and laboratory measurements. The prototype motors of 250 – 300 kW and 140 Hz were tested to verify the results. Utilization factor data are given for several other prototypes the largest of which delivers 1000 kW at 12000 min-1.

Tectonic-palaeoenvironmental forcing of clay-mineral assemblages in nonmarine settings: the Oligocene-Miocene As Pontes Basin (Spain)

Two small, alluvial-lacustrine subbasins developed during the early restraining overstep stages of the Oligocene-Miocene As Pontes strike-slip Basin (NW Spain). Later, the basin evolved into a restraining bend stage and an alluvial-swamp-dominated depositional framework developed. The palaeobiological record demonstrates that the Oligocene-Miocene palaeoclimate in NW Spain was subtropical, warm and humid to subhumid.

An active zone characterized by slow normal faults, the northwestern margin of the València trough (NE Iberia): a review

The northwestern margin of the Valencia trough is an area of low strain characterized by slow normal faults and low to moderate seismicity. Since the mid 1990s this area has been the subject of a number of studies on active tectonic which have proposed different approaches to the location of active faults and to the calculation of the parameters that describe their seismic cycle. Fifty-six active faults have been found and a classification has been made in accordance with their characteristics: a) faults with clear evidence of large paleo-, historic or instrumental earthquakes (2/56); b) faults with evidence of accumulated activity during the Plio-Quaternary and with associated instrumental seismicity (7/56); c) faults with evidence of accumulated activity during the Plio-Quaternary and without associated instrumental seismicity (17/56); d) faults with associated instrumental seismicity and without evidence of accumulated activity during the Plio-Quaternary (30/56), and e) faults without evidence of activity or inactive faults. The parameters that describe the seismic cycle of these faults have been evaluated by different methods that use the geological data obtained for each fault except when paleoseismological studies were available. This classification can be applied to other areas with low slip faults because of the simplicity of the approaches adopted. This study reviews the different approaches proposed and describes the active faults located, highlighting the need a) to better understand active faults in slow strain zones through paleoseismological studies, and b) to include them in seismic hazard studies.

Evidence for coseismic events of recurrent prehistoric deformation along the Alhama de Murcia fault, southeastern Spain

The Alhama de Murcia fault is a 85 km long oblique-slip fault, and is related to historical and instrumental seismic activity. A paleoseismic analysis of the Lorca-Totana sector of the fault containing MSK I=VIII historical earthquakes was made in order to identify and quantify its seismic potential. We present 1) the results of the neotectonic, structural and geomorphological analyses and, 2) the results of trenching. In the study area, the Alhama de Murcia fault forms a depressed corridor between two strands, the northwestern fault with morphological and structural features of a reverse component of slip, bounding the La Tercia range to the South, and the southeastern fault strand with evidence of sinistral oblique strike-slip movement. The offset along this latter fault trapped the sediments in transit from the La Tercia range towards the Guadalentín depression. The most recent of these sediments are arranged in three generations of alluvial fans and terraces. The first two trenches were dug in the most recent sediments across the southeastern fault strand. The results indicate a coseismic reverse fault deformation that involved the sedimentary sequence up to the intermediate alluvial fan and the Holocene terrace deposits. The sedimentary evolution observed in the trenches suggests an event of temporary damming of the Colmenar creek drainage to the South due to uplifting of the hanging wall during coseismic activation of the fault. Trench, structural and sedimentological features provide evidence of at least three coseismic events, which occurred after 125,000 yr. The minimum vertical slip rate along the fault is 0.06 mm/yr and the average recurrence period should not exceed 40,000 yr in accordance with the results obtained by fan topographic profiling. Further absolute dating is ongoing to constrain these estimates.

Phase behavior and rheological analysis of reverse liquid crystals and W/I2, W/H2 gel emulsions using an amphiphilic block copolymer.

This article reports the phase behavior determi- nation of a system forming reverse liquid crystals and the formation of novel disperse systems in the two-phase region. The studied system is formed by water, cyclohexane, and Pluronic L-121, an amphiphilic block copolymer considered of special interest due to its aggregation and structural proper- ties. This system forms reverse cubic (I2) and reverse hexagonal (H2) phases at high polymer concentrations. These reverse phases are of particular interest since in the two-phase region, stable high internal phase reverse emulsions can be formed. The characterization of the I2 and H2 phases and of the derived gel emulsions was performed with small-angle X-ray scattering (SAXS) and rheometry, and the influence of temperature and water content was studied. TheH2 phase experimented a thermal transition to an I2 phase when temperature was increased, which presented an Fd3m structure. All samples showed a strong shear thinning behavior from low shear rates. The elasticmodulus (G0) in the I2 phase was around 1 order of magnitude higher than in theH2 phase. G0 was predominantly higher than the viscousmodulus (G00). In the gel emulsions,G0 was nearly frequency-independent, indicating their gel type nature. Contrarily to water-in-oil (W/O) normal emulsions, in W/I2 and W/H2 gel emulsions, G0, the complex viscosity (|η*|), and the yield stress (τ0) decreased with increasing water content, since the highly viscous microstructure of the con- tinuous phase was responsible for the high viscosity and elastic behavior of the emulsions, instead of the volumefraction of dispersed phase and droplet size. A rheological analysis, in which the cooperative flow theory, the soft glass rheology model, and the slip plane model were analyzed and compared, was performed to obtain one single model that could describe the non-Maxwellian behavior of both reverse phases and highly concentrated emulsions and to characterize their microstructure with the rheological properties.

Present-day deformation of the Pyrenees revealed by GPS surveying and earthquake focal mechanisms until 2011

The Pyrenean mountain range is a slowly deforming belt with continuous and moderate seismic activity. To quantify its deformation field, we present the velocity field estimated from a GPS survey of the Pyrenees spanning 18 yr. The PotSis and ResPyr networks, including a total of 85 GPS sites, were installed and first measured in 1992 and 1995 1997, respectively, and remeasured in 2008 and 2010. We obtain a deformation field with velocities less than 1 mm yr−1 across the range. The estimated velocities for individual stations do not differ significantly from zero with 95 per cent confidence. Even so, we estimate a maximum extensional horizontal strain rate of 2.0 ± 1.7 nanostrain per year in a N S direction in the western part of the range. We do not interpret the vertical displacements due to their large uncertainties. In order to compare the horizontal strain rates with the seismic activity, we analyse a set of 194 focal mechanisms using three methods: (i) the 'r' factor relating their P and T axes, (ii) the stress tensors obtained by fault slip inversion and (iii) the strain-rate tensors. Stress and strain-rate tensors are estimated for: (i) the whole data set, (ii) the eastern and western parts of the range separately, and (iii) eight zones, which are defined based on the seismicity and the tectonic patterns of the Pyrenees. Each of these analyses reveals a lateral variation of the deformation style from compression and extension in the east to extension and strike-slip in the west of the range. Although the horizontal components of the strain-rate tensors estimated from the seismic data are slightly smaller in magnitude than those computed from the GPS velocity field, they are consistent within the 2σ uncertainties. Furthermore, the orientations of their principal axes agree with the mapped active faults.

Separation of Gas Bubbles from Circulation Waters and Pulp Suspensions in Open Channel Flow

The objective of this thesis was to study the removal of gases from paper mill circulation waters experimentally and to provide data for CFD modeling. Flow and bubble size measurements were carried out in a laboratory scale open gas separation channel. Particle Image Velocimetry (PIV) technique was used to measure the gas and liquid flow fields, while bubble size measurements were conducted using digital imaging technique with back light illumination. Samples of paper machine waters as well as a model solution were used for the experiments. The PIV results show that the gas bubbles near the feed position have the tendency to escape from the circulation channel at a faster rate than those bubbles which are further away from the feed position. This was due to an increased rate of bubble coalescence as a result of the relatively larger bubbles near the feed position. Moreover, a close similarity between the measured slip velocities of the paper mill waters and that of literature values was obtained. It was found that due to dilution of paper mill waters, the observed average bubble size was considerably large as compared to the average bubble sizes in real industrial pulp suspension and circulation waters. Among the studied solutions, the model solution has the highest average drag coefficient value due to its relatively high viscosity. The results were compared to a 2D steady sate CFD simulation model. A standard Euler-Euler k-ε turbulence model was used in the simulations. The channel free surface was modeled as a degassing boundary. From the drag models used in the simulations, the Grace drag model gave velocity fields closest to the experimental values. In general, the results obtained from experiments and CFD simulations are in good qualitative agreement.