<p>A new large eddy simulation (LES) code is developed and used to investigate non-local features of turbulent planetary boundary layers (PBLs). The LES code is based on filtered Navier-Stokes equations, which describe motions of incompressible, Boussinesq fluid at high Reynolds numbers. The code computes directly large-scale, non-universal turbulence in the PBL whereas small-scale, universal turbulence is parameterized by a dynamic mixed subgrid closure. The LES code is thoroughly tested against high quality laboratory and field data. </p><p>This study addresses non-local properties of turbulence which emphasis on the stable stratification. Its basic results are as follows. The flow stability in PBLs is generally caused by two mechanisms: the negative buoyancy force (in the stable density stratification) and the Coriolis force (in the rotating system). The latter stabilizes the flow if the earth’s vorticity and the turbulent vorticity are anti-parallel. The Coriolis force stability suppresses large-scale turbulence and makes large eddies asymmetric. The density stratification suppresses vertical scales of turbulence. Joint actions of the Coriolis and the buoyancy forces result in a more complex behavior of turbulence. Particularly, the layers of vigorous turbulence may appear in the course of development of low-level jets in baroclinic atmosphere. </p><p>Non-local effects determine integral measures of PBLs, first of all the PBL depth. This study clearly demonstrates its pronounced dependences on the Coriolis parameter, the Kazanski-Monin internal stability parameter, and newly introduced imposed-stability and baroclinicity parameters. An LES database is created and used to validate an advanced PBL-depth formulation. LES support the idea that PBLs interact with the stably stratified free flow through the radiation of gravity waves, excited by large turbulent eddies at the interface. </p>

<p>In this study the flux attenuation due to sensor displacement has been investigated over sea using an extensive set of data from the "Ocean Horizontal Array Turbulence Study". All previous investigations of the flux attenuation have been performed over land.</p><p>A function developed for correcting fluxes in the homogenous surface layer was compared to measured flux attenuation. This investigation revealed the possibility to find new functions describing the flux attenuation when measurements are carried out over sea. From the measured flux attenuation studied here a change in the form of correction functions was required to improve the estimated flux loss. The most significant difference found in this report compared to the previous landbased study Horst (2006) is for stable conditions, where significantly less flux loss is found over sea. Two new functions describing the attenuation due to sensor displacement over sea have been constructed.</p><p>One of these expressions has a discontinuity at z/L = 0. This is supported by measured flux attenuation. A reasonable interpretation is; however, that this discontinuity is caused by two separate turbulence regimes near neutrality on the stable and unstable side respectively. The discontinuity is thus not believed to be an effect merely of stability. A second correction function which is continuous over all stabilities has therefore also been constructed. These two functions and the correction function from Horst (2006) have been compared to measured flux loss. Based on this comparison the continuous correction function is recommended for correcting scalar fluxes measured over sea. It should be noted, however, that this expression only describes the mean attenuation and has been constructed from measurements at 5 and 5.5 m above mean sea level.</p><p>The theoretical basis used in the development of the function for flux attenuation over land allows for a direct link between a spectral shape and the attenuation expression. This link has been preserved for the new expressions presented in this report. The spectral shape corresponding to the continuous correction function has been compared to measured mean cospectra and also to the cospectra from Horst (2006) corresponding to crosswind displacements.</p><p>At a height of 10 m and a sensor displacement of 0.2 m the mean flux attenuation is about 1.3-4% in the stability interval −1 < z/L < 1.5 when using the new correction functions presented in this report.</p>

The objective of this research is to develop an efficient and accurate methodology to resolve flow non-linearity of fluid-structural interaction. To achieve this purpose, a numerical strategy to apply the detached-eddy simulation (DES) with a fully coupled fluid-structural interaction model is established for the first time. The following novel numerical algorithms are also created: a general sub-domain boundary mapping procedure for parallel computation to reduce wall clock simulation time, an efficient and low diffusion E-CUSP (LDE) scheme used as a Riemann solver to resolve discontinuities with minimal numerical dissipation, and an implicit high order accuracy weighted essentially non-oscillatory (WENO) scheme to capture shock waves. The Detached-Eddy Simulation is based on the model proposed by Spalart in 1997. Near solid walls within wall boundary layers, the Reynolds averaged Navier-Stokes (RANS) equations are solved. Outside of the wall boundary layers, the 3D filtered compressible Navier-Stokes equations are solved based on large eddy simulation(LES). The Spalart-Allmaras one equation turbulence model is solved to provide the Reynolds stresses in the RANS region and the subgrid scale stresses in the LES region. An improved 5th order finite differencing weighted essentially non-oscillatory (WENO) scheme with an optimized epsilon value is employed for the inviscid fluxes. The new LDE scheme used with the WENO scheme is able to capture crisp shock profiles and exact contact surfaces. A set of fully conservative 4th order finite central differencing schemes are used for the viscous terms. The 3D Navier-Stokes equations are discretized based on a conservative finite differencing scheme, which is implemented by shifting the solution points half grid interval in each direction on the computational domain. The solution points are hence located in the center of the grid cells in the computational domain (not physical domain). This makes it possible to use the same code structure as a 2nd order finite volume method. A finite differencing high order WENO scheme is used since a finite differencing WENO scheme is much more efficient than a finite volume WENO scheme. The unfactored line Gauss-Seidel relaxation iteration is employed for time marching. For the time accurate unsteady simulation, the temporal terms are discretized using the 2nd order accuracy backward differencing. A pseudo temporal term is introduced for the unsteady calculation following Jameson's method. Within each physical time step, the solution is iterated until converged based on pseudo time step. A general sub-domain boundary mapping procedure is developed for arbitrary topology multi-block structured grids with grid points matched on sub-domain boundaries. The interface of two adjacent blocks is uniquely defined according to each local mesh index system (MIS) which is specified independently. A pack/unpack procedure based on the definition of the interface is developed to exchange the data in a 1D array to minimize data communication. A secure send/receive procedure is employed to remove the possibility of blocked communication and achieve optimum parallel computation efficiency. Two terms, "Order" and "Orientation", are introduced as the logics defining the relationship of adjacent blocks. The domain partitioning treatment of the implicit matrices is to simply discard the corner matrices so that the implicit Gauss-Seidel iteration can be implemented within each subdomain. This general sub-domain boundary mapping procedure is demonstrated to have high scalability. Extensive numerical experiments are conducted to test the performance of the numerical algorithms. The LDE scheme is compared with the Roe scheme for their behavior with RANS simulation. Both the LDE and the Roe scheme can use high CFL numbers and achieve high convergence rates for the algebraic Baldwin-Lomax turbulence model. For the Spalart-Allmaras one equation turbulence model, the extra equation changes the Jacobian of the Roe scheme and weakens the diagonal dominance. It reduces the maximum CFL number permitted by the Roe scheme and hence decreases the convergence rate. The LDE scheme is only slightly affected by the extra equation and maintains high CFL number and convergence rate. The high stability and convergence rate using the Spalart-Allmaras one equation turbulence model is important since the DES uses the same transport equation for the turbulence stresses closure. The RANS simulation with the Spalart-Allmaras one equation turbulence model is the foundation for DES and is hence validated with other transonic flows including a 2D subsonic flat plate turbulent boundary layer, 2D transonic inlet-diffuser, 2D RAE2822 airfoil, 3D ONERA M6 wing, and a 3D transonic duct with shock boundary layer interaction. The predicted results agree very well with the experiments. The RANS code is then further used to study the slot size effect of a co-flow jet (CFJ) airfoil. The DES solver with fully coupled fluid-structural interaction methodology is validated with vortex induced vibration of a cylinder and a transonic forced pitching airfoil. For the cylinder, the laminar Navier-Stokes equations are solved due to the low Reynolds number. The 3D effects are observed in both stationary and oscillating cylinder simulation because of the flow separations behind the cylinder. For the transonic forced pitching airfoil DES computation, there is no flow separation in the flow field. The DES results agree well with the RANS results. These two cases indicate that the DES is more effective on predicting flow separation. The DES code is used to simulate the limited cycle oscillation of NLR7301 airfoil. For the cases computed in this research, the predicted LCO frequency, amplitudes, averaged lift and moment, all agree excellently with the experiment. The solutions appear to have bifurcation and are dependent on the initial perturbation. The developed methodology is able to capture the LCO with very small amplitudes measured in the experiment. This is attributed to the high order low diffusion schemes, fully coupled FSI model, and the turbulence model used. This research appears to be the first time that a numerical simulation of LCO matches the experiment. The DES code is also used to simulate the CFJ airfoil jet mixing at high angle of attack. In conclusion, the numerical strategy of the high order DES with fully coupled FSI model and parallel computing developed in this research is demonstrated to have high accuracy, robustness, and efficiency. Future work to further maturate the methodology is suggested.

River plumes often develop in complex environments, where variable coastal and bottom topography, ambient currents, winds and tides may play important roles in shaping the plume evolution. When all these factors are present, the plume dynamics may become intricate and unclear. The objective of this study is to understand the processes controlling the dynamics of a large river plume that is affected by strong boundary currents, variable winds and complex topography. The Mississippi River (MR) plume is the study case of this dissertation work, and focus is given to the interactions between the plume and the offshore circulation of the Gulf of Mexico (GoM). A series of numerical experiments was designed to investigate the impact of different factors on the development of a large scale river plume in scenarios with variable degrees of complexity. First, a box-like model with an idealized estuary was designed to address the general development of a mid-latitude river plume and assess the variability of the plume with changes in the outflow conditions at the river mouth. The structure and development of the plume in the flat-bottom, receiving basin was highly dependent on the degree of freshwater mixing at the source. Larger freshwater mixing enhanced the estuarine gravitational circulation and modified the dynamical balance at the estuary mouth. Those changes effectively modified the shape of the bulge and length/transport scales of the coastal current. Sloping-bottom conditions further modified the development of the plume. Secondly, a Northern GoM model was designed and numerical experiments were conducted to investigate the specific dynamics of the MR plume, in the presence of both shelf and basin-wide circulation. In particular, buoyancy-driven (due to the MR and all other major Northern GoM rivers) and wind-driven currents were studied on the shelf, while the extension of the Loop Current and associated frontal eddies were considered as major factors in the shelf to offshore interactions; wind-driven, shelfbreak eddies were also considered. Process-oriented experiments demonstrate that westerly and southerly winds promoted the development of a surface Ekman layer that enhances the offshore advection of plume waters. The steep topography in the vicinity of the MR Delta was a favorable condition for that process. When the MR plume was subject to a full-blown scenario (realistically-forced experiment nested within a large-scale model), complex interactions between wind-driven and eddy-driven dynamics determined the fate of the plume waters. Offshore removal is a frequent plume pathway, and the offshore transport can be as large as the wind-driven shelf transport. The offshore pathways depend on the position of the eddies near the shelf edge, their life span and the formation of eddy pairs that generate coherent cross-shelf flows. Strong eddy-plume interactions were observed when the Loop Current (LC) system impinged against the shelfbreak, causing the formation of coherent, narrow low-salinity bands that extended toward the Gulf interior. The offshore transport of MR water is a year-round process, but the interactions between the MR plume and the LC system have large inter-annual variability. Plume to LC interactions are determined by episodic northward intrusions of the LC system in the NGoM. The interactions are dictated by the proximity of the LC system to the MR Delta and by wind effects. On average, plume to LC interactions correspond to ~ 12 % of the year-round, total freshwater transport near the MR Delta, but this percentage can go up to 30 % in individual years. At the time of the plume to LC interactions, an average value of LC freshwater entrainment was estimated to be ~ 4,150 m3 s-1. The findings presented here are a major contribution toward the understanding of the cross-marginal and basin-wide transport of MR waters by a large-scale current system, and the connectivity to remote regions, such as the South Florida region and the Florida Keys.

La technique micro météorologique de covariance de turbulences est utilisée pour estimer les échanges de CO2 entre les écosystèmes et latmosphère. Des centaines de sites instrumentés, répartis dans le monde entier, lutilisent désormais pour étudier une grande variété décosystèmes. Cette technique est cependant entachée dune erreur systématique lorsquelle est appliquée sur des couverts hauts comme des forêts, en conditions atmosphériques stables, cest-à-dire essentiellement pendant les nuits peu venteuses et sans couverture nuageuse. Pendant ces périodes, en effet, le transport turbulent serait concurrencé par un autre mécanisme de transport qui est ladvection. Dans ce travail, la présence dadvection a été testée sur le site expérimental forestier de Vielsalm (Belgique) et son importance a été évaluée. A cette fin, un dispositif expérimental permettant des mesures de vitesse de vent, de concentration de CO2 et de température de lair à lintérieur de la forêt a été installé. Il a permis la mise en évidence, pendant les périodes stables, découlements gravitationnels se réalisant près du sol suite au refroidissement des surfaces et à la présence dune faible pente. Il a été montré que ces écoulements étaient responsables du transport advectif de CO2. Une analyse dincertitude a été menée à laide de campagnes de mesures spécifiques. Elle a conclu à la faisabilité des mesures de gradients verticaux et surtout horizontaux de CO2 sur le site mais a montré que le plus grand facteur dincertitude portait sur les estimations de la composante verticale de la vitesse au-dessus de et dans la forêt. Malgré ces incertitudes, une analyse fine des épisodes gravitationnels a permis de mettre en évidence un mécanisme cohérent liant les écoulements dair et le champ des concentrations de CO2 et permettant de mieux comprendre comment le CO2 pouvait être transporté latéralement et verticalement par les écoulements gravitationnels. Finalement, la faisabilité dune correction basée sur lestimation des termes advectifs a été évaluée. Il a été montré que les incertitudes portant sur ladvection étaient trop importantes pour permettre daméliorer sensiblement le bilan nocturne des flux de CO2 au moyen de cette méthode. Celle-ci savère toutefois intéressante pour mieux comprendre les processus de transport à loeuvre dans un couvert forestier./The micrometeorological technique of eddy-covariance is used for the estimation of the CO2 exchange between the ecosystems and the atmosphere. Hundreds of instrumented sites, spread all over the world, use it henceforth to study a great variety of ecosystems. This technique is however affected by a systematic error when applied above tall canopies like forests, in stable atmospheric conditions, i.e. primarily during non windy nights without cloud cover. Indeed, during these periods the turbulent transport would be competed with by another transport mechanism which is called advection. In this work, the presence of advection has been tested on the experimental forested site of Vielsalm (Belgium) and its importance has been evaluated. For this purpose, an experimental set-up allowing the measurements of wind velocity, CO2 concentration and temperature of the air inside the forest has been installed. It has allowed the description, for the stable periods, of gravitational flows being carried out close to the ground due to the cooling of surfaces and the presence of a weak slope. These flows were shown to be responsible for advective CO2 transport. An uncertainty analysis has been carried out using dedicated measurement campaigns. It has conclude with the feasibility of measurements of vertical and especially horizontal CO2 gradients on the site but has shown that the greatest factor of uncertainty related to the estimate of the vertical velocity component above and in the forest. In spite of these uncertainties, a fine analysis of the gravitational episodes has made it possible to highlight a coherent mechanism linking the flow field and the CO2 concentration field and making it possible to better understand how CO2 could be transported laterally and vertically by the gravitational flows. Finally, the feasibility of a correction based on the estimate of the advective terms has been evaluated. It has been shown that uncertainties relating to advection were too important to make it possible to appreciably improve the night assessment of CO2 fluxes by means of this method. This one proves however interesting for better understanding the processes of transport at work in a forest cover.

When an electrically conducting fluid moves through a magnetic field, fluid mechanics and electromagnetism are coupled. This interaction is the object of magnetohydrodynamics, a discipline which covers a wide range of applications, from electromagnetic processing to plasma- and astro-physics. In this dissertation, the attention is restricted to turbulent liquid metal flows, typically encountered in steel and aluminium industries. Velocity measurements in such flows are extremely challenging because liquid metals are opaque, hot and often corrosive. Therefore, non-intrusive measurement devices are essential. One of them is the Lorentz force flowmeter. Its working principle is based on the generation of a force acting on a charge, which moves in a magnetic field. Recent studies have demonstrated that this technique can measure efficiently the mean velocity of a liquid metal. In the existing devices, however, the measurement depends on the electrical conductivity of the fluid. In this work, a novel version of this technique is developed in order to obtain measurements that are independent of the electrical conductivity. This is particularly appealing for metallurgical applications, where the conductivity often fluctuates in time and space. The study is entirely numerical and uses a flexible computational method, suitable for industrial flows. In this framework, the cost of numerical simulations increases drastically with the level of turbulence and the geometry complexity. Therefore, the simulations are commonly unresolved. Large eddy simulations are then very promising, since they introduce a subgrid model to mimic the dynamics of the unresolved turbulent eddies. The first part of this dissertation focuses on the quality and reliability of unresolved numerical simulations. The attention is drawn on the ambiguity that may arise when interpretating the results. Owing to coarse resolutions, numerical errors affect the performances of the discrete model, which in turn looses its physical meaning. In this work, a novel implementation of the turbulent strain rate appearing in the models is proposed. As opposed to its usual discretisation, the present strain rate is in accordance with the discrete equations of motion. Two types of flow are considered: decaying turbulence located far from boundaries, and turbulent flows between two parallel and infinite walls. Particular attention is given to the balance of resolved kinetic energy, in order to assess the role of the model. The second part of this dissertation deals with a novel version of Lorentz force flowmeters, consisting in one or two coils placed around a circular pipe. The forces acting on each coil are recorded in time as the liquid metal flows through the pipe. It is highlighted that the auto- or cross-correlation of these forces can be used to determine the flowrate. The reliability of the flowmeter is first investigated with a synthetic velocity profile associated to a single vortex ring, which is convected at a constant speed. This configuration is similar to the movement of a solid rod and enables a simple analysis of the flowmeter. Then, the flowmeter is applied to a realistic three-dimensional turbulent flow. In both cases, the influence of the geometrical parameters of the coils is systematically assessed.

A new large eddy simulation (LES) code is developed and used to investigate non-local features of turbulent planetary boundary layers (PBLs). The LES code is based on filtered Navier-Stokes equations, which describe motions of incompressible, Boussinesq fluid at high Reynolds numbers. The code computes directly large-scale, non-universal turbulence in the PBL whereas small-scale, universal turbulence is parameterized by a dynamic mixed subgrid closure. The LES code is thoroughly tested against high quality laboratory and field data. This study addresses non-local properties of turbulence which emphasis on the stable stratification. Its basic results are as follows. The flow stability in PBLs is generally caused by two mechanisms: the negative buoyancy force (in the stable density stratification) and the Coriolis force (in the rotating system). The latter stabilizes the flow if the earth’s vorticity and the turbulent vorticity are anti-parallel. The Coriolis force stability suppresses large-scale turbulence and makes large eddies asymmetric. The density stratification suppresses vertical scales of turbulence. Joint actions of the Coriolis and the buoyancy forces result in a more complex behavior of turbulence. Particularly, the layers of vigorous turbulence may appear in the course of development of low-level jets in baroclinic atmosphere. Non-local effects determine integral measures of PBLs, first of all the PBL depth. This study clearly demonstrates its pronounced dependences on the Coriolis parameter, the Kazanski-Monin internal stability parameter, and newly introduced imposed-stability and baroclinicity parameters. An LES database is created and used to validate an advanced PBL-depth formulation. LES support the idea that PBLs interact with the stably stratified free flow through the radiation of gravity waves, excited by large turbulent eddies at the interface.

In this study the flux attenuation due to sensor displacement has been investigated over sea using an extensive set of data from the "Ocean Horizontal Array Turbulence Study". All previous investigations of the flux attenuation have been performed over land. A function developed for correcting fluxes in the homogenous surface layer was compared to measured flux attenuation. This investigation revealed the possibility to find new functions describing the flux attenuation when measurements are carried out over sea. From the measured flux attenuation studied here a change in the form of correction functions was required to improve the estimated flux loss. The most significant difference found in this report compared to the previous landbased study Horst (2006) is for stable conditions, where significantly less flux loss is found over sea. Two new functions describing the attenuation due to sensor displacement over sea have been constructed. One of these expressions has a discontinuity at z/L = 0. This is supported by measured flux attenuation. A reasonable interpretation is; however, that this discontinuity is caused by two separate turbulence regimes near neutrality on the stable and unstable side respectively. The discontinuity is thus not believed to be an effect merely of stability. A second correction function which is continuous over all stabilities has therefore also been constructed. These two functions and the correction function from Horst (2006) have been compared to measured flux loss. Based on this comparison the continuous correction function is recommended for correcting scalar fluxes measured over sea. It should be noted, however, that this expression only describes the mean attenuation and has been constructed from measurements at 5 and 5.5 m above mean sea level. The theoretical basis used in the development of the function for flux attenuation over land allows for a direct link between a spectral shape and the attenuation expression. This link has been preserved for the new expressions presented in this report. The spectral shape corresponding to the continuous correction function has been compared to measured mean cospectra and also to the cospectra from Horst (2006) corresponding to crosswind displacements. At a height of 10 m and a sensor displacement of 0.2 m the mean flux attenuation is about 1.3-4% in the stability interval −1 &lt; z/L &lt; 1.5 when using the new correction functions presented in this report.

The water and heat exchange in thesoil-snow-vegetation-atmosphere system was studied in order toimprove the quantitative knowledge of land surface processes.In this study, numerical simulation models and availabledatasets representing arable land, sub-alpine snowpack, andboreal forest were evaluated at both diurnal and seasonaltimescales. Surface heat fluxes, snow depth, soil temperatures andmeteorological conditions were measured at an agriculturalfield in central Sweden during three winters and two summersfrom 1997 to 2000 within the WINTEX project. A one-dimensionalsimulation model (COUP) was used to simulate the water and heatbalance of the field. Comparison of simulated and measured heatfluxes in winter showed that parameter values governing theupper boundary condition were more important for explainingmeasured fluxes than the formulation of the internal mass andheat balance of the snow cover. The assumption of steady stateheat exchange between the surface and the reference height wasinadequate during stable atmospheric conditions. Independentestimates of the soil heat and water balance together with thecomparison of simulated and measured surface heat fluxes showedthat the eddy-correlation estimates of latent heat fluxes fromthe arable field were on average 40 % too low. The ability of a multi-layered snowpack model (SNTHERM) tosimulate the layered nature of a sub-alpine snowpack wasevaluated based on a dataset from Switzerland. The modelsimulated the seasonal development of snow depth and densitywith high accuracy. However, the models ability to reproducethe strong observed snowpack layering was limited by theneglection of the effect of snow microstructure on snowsettling, and a poor representation of water redistributionwithin the snowpack. The representation of boreal forest in the land surfacescheme used within a weather forecast (ECMWF) model was testedwith a three-year dataset from the NOPEX forest site in centralSweden. The new formulation with separate energy balances forvegetation and the soil/snow beneath the tree cover improvedthe simulation of seasonal and diurnal variations in latent andsensible heat flux. Further improvements of simulated latentheat fluxes were obtained when seasonal variation in vegetationproperties was introduced. Application of the COUP model withthe same dataset showed that simulation of evaporation fromintercepted snow contributed to a better agreement with themeasured sensible heat flux above forests, but also indicatedthat the measurements might have underestimated latent heatflux. The winter sensible heat flux above the forest wasfurther improved if an upper limit of the aerodynamicresistance of 500 s m-1 was applied for stable conditions. A comparison of the water and heat balance of arable landand forest confirmed the general knowledge of the differencesbetween these two surface types. The forest contributed withconsiderably more sensible heat flux to the atmosphere than thearable land in spring and summer due to the lower albedo andrelatively less latent heat flux. Latent heat flux from theforest was higher in winter due to the evaporation ofintercepted snow and rain. The net radiation absorbed by theforest was 60 % higher than that absorbed by the arable land,due to the lower surface albedo in winter. Key words:soil; snow; land surface heat exchange;forest; arable land; eddy-correlation. / QC 20100614

Turbocharging the internal combustion (IC) engine is a common technique to increase the power density. If turbocharging is used with the downsizing technique, the fuel consumption and pollution of green house gases can be decreased. In the turbocharger, the energy of the engine exhaust gas is extracted by expanding it through the turbine which drives the compressor by a shaft. If a turbocharged IC engine is compared with a natural aspirated engine, the turbocharged engine will be smaller, lighter and will also have a better efficiency, due to less pump losses, lower inertia of the system and less friction losses. To be able to further increase the efficiency of the IC engine, the understanding of the highly unsteady flow in turbochargers must be improved, which then can be used to increase the efficiency of the turbine and the compressor. The main objective with this thesis has been to enhance the understanding of the unsteady flow in turbocharger and to assess the sensitivity of inflow conditions on the turbocharger performance. The performance and the flow field in a radial turbocharger turbine working under both non-pulsatile and pulsatile flow conditions has been assessed by using Large Eddy Simulation (LES). To assess the effects of different operation conditions on the turbine performance, different cases have been considered with different perturbations and unsteadiness of the inflow conditions. Also different rotational speeds of the turbine wheel were considered. The results show that the turbine cannot be treated as being quasi-stationary; for example,the shaft power varies for different frequencies of the pulses for the same amplitude of mass flow. The results also show that perturbations and unsteadiness that are created in the geometry upstream of the turbine have substantial effects on the performance of the turbocharger. All this can be summarized as that perturbations and unsteadiness in the inflow conditions to the turbine affect the performance. The unsteady flow field in ported shroud compressor has also been assessed by using LES for two different operational points. For an operational point near surge, the flow field in the entire compressor stage is unsteady, where the driving mechanism is an unsteadiness created in the volute. For an operational point far away from surge, the flow field in the compressor is relatively much more steady as compared with the former case. Although the stable operational point exhibits back-flow from the ported shroud channels, which implies that the flow into the compressor wheel is disturbed due to the structures that are created in the shear layer between the bulk flow and the back-flow from the ported shroud channels. / QC20100622