Electrification of vehicles is necessary to combat greenhouse gas emissions, whichcauses global warming and climate change. There has been a demand in the use ofBattery Electric Vehicles due to this increased awareness of sustainability. However,they have been beset with issues such as range and conservation of energy. A partof the solution could be an energy-efficient thermal management system. A 1-Dimensional (1-D) complete vehicle thermal model in GT-SUITE is used topredict the energy efficiency of the vehicle at the conceptual stage. The model helpspredict results quickly and is used for complete system-level simulations. This thesisfocuses on optimizing the method for predicting the realistic and accurate energyconsumption of the thermal management system in the vehicle at the 1-D level. TheCompact Modular Architecture (CMA) platform, which is the vehicle platform usedin current production cars by Volvo Cars will be used for the study and all studiesare performed using a standardized drive cycle. Initial sensitivity studies with afully open grille are performed to understand the operating points of the variouscomponents of interest to investigate. The mass flow rate, ambient temperature,battery temperature, and cooling fan speed are varied.The Active Grille Shutters (AGS) which provides aerodynamic benefit at high speedis implemented in the existing thermal model which could previously accommodateonly the fully open grille. This allows the grille shutters to vary at different anglesbased on cooling demand. The previous existing thermal model method also neededto be optimized to accommodate the AGS. The cooling fan control logic needed tobe improved for better accuracy and energy consumption prediction. Furthermore,the grille shutters and cooling fan speed is needed to regulate the amount of air flowthrough the heat exchangers for the model to behave as close to a real productionvehicle. A code was developed which generated fan speed and grille shutter anglesbased on mass flow rate values to input in the model.Further investigations were made with the optimized thermal model with AGS tostudy the influence of additional mass flow rate on the mass power consumptionof the thermal management system components of interest. It was observed in the initial sensitivity studies that the additional mass flow rate saw significant power savings. However, with the implementation of AGS and additional mass airflowinto the system, the power due to the variation of shutters is taken into account.The results indicate that the total power consumption gradually decreases with theincreasing mass flow rate. But, this is up to a certain extent where the energyconsumption due to the shutter opening takes over the overall power consumptionof the vehicle and overcomes the savings seen by other components in the system causing the total power to increase.

Wading is the phenomenon where a vehicle drives through water with a relatively deep water level. Sincea large portion of the vehicle is submerged in water it can affect the driveability and function of individualcomponents. Wading is therefore an important phenomenon to be aware of especially today where society moves towards alternative energy sources. This includes water sensitive components when contact with water can generate major consequences. Previous knowledge and experience of wading has been from performing physical tests, but using Computational Fluid Dynamics (CFD) to examine the phenomenon can accelerate the iterative design process. In this thesis, numerical method of wading simulations on complete heavy-duty vehicles using the software STAR-CCM+ are developed. Furthermore, the results from the numerical methods are validated against results from physical tests performed at Scania’s test facility in Södertälje. The numerical methods are divided into a simplified model of a Battery Electric Vehicle (BEV) and a detailed geometry of a gas-driven vehicle from Scania. Beside dividing the wading scenario into the geometries, two different methods are developed, Wave and Wading. The Wave-method includes the vehicle standing still while a water wave is fed in through the inlet of the domain, i.e. allowed to flush over the vehicle, with a velocity of 3.6 km/h and 8 km/h. This method is implemented for both a generic simplified BEV truck and a detailed real-life Scania truck. For the Wading-method, motion is applied to the vehicle where itis driving with a velocity of 3.6 km/h through a digital twin of the water trench available at the test facility. This method is further divided into two cases, Zero Gap and Floating, where the difference is the distance between the tires of the vehicle and ground of the domain. The Floating-case includes a 10 cm distance and the Zero Gap-case has no gap between the tires and ground. The Wading-method is only implemented for the simplified geometry due to the computational cost and complexity. All methods use the Volume of Fluid (VOF) method for multiphase modelling and the Zero Gap-case uses Overset Mesh for modelling motion. The validation of the simulations focuses on the water behaviour such as water surface topology and water flowing inside the vehicle while wading. The results for the Wave-method with both the simplified and detailed truck at 8 km/h shows similarities in the water surface topology between the numerical model and the physical test. The simulations of the Wading-method is not visualising any similarities since the visible wave pattern are few and unclear in the numerical model. An isosurface is used to visualise the surface of the water which generated a smooth topology since no other options, such as vector fields, are added. It is found that the water movement inside the vehicle will affect water sensitive areas, e.g. on the battery packs. It is concluded that the derived methods are a first draft and should be directed towards future development in optimising the methods to lower the computational cost, but also to improve the capturing of the interface between the two phases. Due to instability and computational cost the detailed geometry is not implemented in the Wading-method. The methods are adapted to use different vehicle types since the simplified and detailed geometry are a BEV and a gas-driven truck respectively.

Boundary layer separation is usually an unwanted phenomenon in most technical applications as for instance on airplane wings, on ground vehicles and in internal flows such as diffusers. If separation occurs it leads to loss of lift, higher drag and results in energy losses. It is therefore important to be able to find methods to control and if possible avoid separation altogether without introducing a too heavy penalty such as increased drag, energy consuming suction etc. In the present work we study one such control method, namely the use of vortex generators (VGs), which are known to be able to hinder turbulent boundary layer separation. We first study the downstream development of streamwise vortices behind pairs and arrays of vortex generators and how the strength of the vortices is coupled to the relative size of the vortex generators in comparison to the boundary layer size. Both the amplitude and the trajectory of the vortices are tracked in the downstream direction. Also the influences of yaw and free stream turbulence on the vortices are investigated. This part of the study is made with hot-wire anemometry where all three velocity components of the vortex structure are measured. The generation of circulation by the VGs scales excellently with the VG blade height and the velocity at the blade edge. The magnitude of circulation was found to be independent of yaw angle. The second part of the study deals with the control effect of vortex generators on three different cases where the strength of the adverse pressure gradient (APG) in a turbulent boundary layer has been varied. In this case the measurements have been made with particle image velocimetry. It was found that the streamwise position where the VGs are placed is not critical for the control effect. For the three different APG cases approximately the same level of circulation was needed to inhibit separation. In contrast to some previous studies we find no evidence of a universal detachment shape factor H12, that is independent of pressure gradient. / QC 20101119

This work is concerned with the study of flow stability and turbulence control - two old but still open problems of fluid mechanics. The topics are distinct and are (currently) approached from different directions and with different strategies. This thesis reflects this diversity in subject with a difference in geometry and, consequently, flow structure: the first problem is approached in the study of the flow in a toroidal pipe, the second one in an attempt to reduce the drag in a turbulent channel flow. The flow in a toroidal pipe is chosen as it represents the common asymptotic limit between spatially developing and helical pipes. Furthermore, the torus represents the smallest departure from the canonical straight pipe flow, at least for small curvatures. The interest in this geometry is twofold: it allows us to isolate the effect of the curvature on the flow and to approach straight as well as helical pipes. The analysis features a characterisation of the steady solution as a function of curvature and the Reynolds number. The problem of forcing fluid in the pipe is addressed, and the so-called Dean number is shown to be of little use, except for infinitesimally low curvatures. It is found that the flow is modally unstable and undergoes a Hopf bifurcation that leads to a limit cycle. The bifurcation and the corresponding eigenmodes are studied in detail, providing a complete picture of the instability. The second part of the thesis approaches fluid mechanics from a different perspective: the Reynolds number is too high for a deterministic description and the flow is analysed with statistical tools. The objective is to reduce the friction exerted by a turbulent flow on the walls of a channel, and the idea is to employ a control strategy independent of the small, and Reynolds number-dependent, turbulent scales. The method of choice was proposed by Schoppa &amp; Hussain [Phys. Fluids 10:1049-1051 (1998)] and consists in the imposition of streamwise invariant, large-scale vortices. The vortices are re-implemented as a volume force, validated and analysed. Results show that the original method only gave rise to transient drag reduction while the forcing version is capable of sustained drag reduction of up to 18%. An analysis of the method, though, reveals that its effectiveness decreases rapidly as the Reynolds number is increased. / <p>QC 20161004</p>

The aim of the present study is to develop subgrid-scale models that are relevant for complex flows and combustion. A stochastic model based on a stochastic Smagorinsky constant with adjustable variance and time scale is proposed. The stochastic model is shown to provide for backscatter of both kinetic energy and scalar variance without causing numerical instabilities. A new subgrid-scale scalar flux model is developed using the same kind of methodology that leads to the explicit algebraic scalar flux model, EASFM, for RANS. The new model predicts the anisotropy of the subgrid-scales in a more realistic way than the eddy diffusion model. Both new models were tested in rotating homogeneous shear flow with a passive scalar. Rogallo’s method of moving the frame with the mean flow to enable periodic boundary conditions was used to simulate homogeneous shear flow. / QC 20101119

Hydropower is an efficient renewable energy source able to regulate electrical grid fluctuations. However, many hydropower plants were built decades ago, and now it is the time for a major refurbishment. The turbine's efficiency is essential and needs to be determined before and after refurbishment. To this end, the flow rate needs to be determined. Amongst different discharge measurement methods, the pressure-time method is relatively inexpensive and easy to perform compared to other methods. In this method, the flow rate is estimated by the integration of the measured differential pressure and the pressure loss due to friction between two cross-sections in a conduit during the deceleration of the liquid mass by closing a valve or guide vanes. The pressure-time method's accuracy depends on how accurate the head loss and the integration endpoint are estimated. Furthermore, the pressure-time method has limitations specified in IEC-60041, which make it challenging to apply on low-head turbines due to the short water passages when the flow is developing. The main focus of work is to improve the accuracy of the pressure-time method and extend its validity for low-head turbine conditions. Numerical simulation and experimental study have been acquired. A CFD model is developed to investigate the effects of the endpoint of integration and friction models on the method's accuracy. The effect of different boundary conditions is studied in the CFD model, and the result is validated with available experimental data. Different frictional models used with the pressure-time method are compared with CFD simulation for the developing and developed flows. A new parameter is suggested to improve deviation related to the flow status; developing and developed. Furthermore, a new methodology is presented, where the flow rate is estimated with the pressure-time method function of several endpoints. Then, experimental investigations of the pressure-time method outside IEC-60041 recommendations for conditions similar to low-head hydropower are presented. A laboratory setup is designed and built to test the pressure-time method. The method is applied for cases with shorter length, smaller UxL, pipe with variable cross-section and shorter distance to irregularity than IEC-60041 recommendations. Different assumptions for calculating the pressure loss and dynamic pressure variation are studied. Moreover, the quasi-steady assumption's accuracy on the head loss estimation and the difference in dynamic pressure are compared with constant values for their coefficients. The systematic uncertainty of the pressure-time method is also calculated based on the Monte Carlo Method (MCM).

The interaction between fluid dynamics and biofilm growth plays a key role in both medical and industrial applications. Biofilms, or bacteria that are embedded in a protective matrix of extracellular polymeric substances, settle on interfaces such as on implanted devices or ship hulls. These biofilms canthen cause infectious diseases or significantly increase drag. In this thesis, we investigate the influence of flow, specifically shear stress, on the development of biofilm. The first paper presents a new facility to investigate biofilm growth in laminar flow cells over extended periods of up to several weeks. Optical coherence tomography is used to obtain three-dimensional scans of the biofilm structure at regular intervals. From these time series, we derive a simple model that relates the growth of an individual microcolony to the growth of the full biofilm depending on the wall shear stress. Additionally, we show that biofilm streamers, thin, flexible filaments that extend into the bulk flow, can form on sharp biofilm structures in laminar flow, even if the substratum is a flat surface. The second contribution is a report detailing preliminary studies on biofilm experiments. We investigate the formation of biofilm in the shear layer behinda backward-facing step. The results indicate a maximum shear stress, beyond which biofilm growth is inhibited. We also provide guidelines for the design of experimental setups for the investigation of the influence of fluid dynamics on biofilm and vice-versa. / Samspelet mellan fluiddynamik och biofilmtillväxt spelar en nyckelroll i både medicinska och industriella tillämpningar. Biofilmer, eller bakterier som är inbäddade i en skyddande matris av extracellulära polymera substanser, sätter sig på ytor som på implanterade enheter eller fartygsskrov. Dessa biofilmer kan sedan orsaka infektionssjukdomar eller avsevärt öka vattenmotståndet. I den här avhandlingen undersöker vi hur flöde, speciellt skjuvspänning, påverkar utvecklingen av biofilm. I den första artikeln presenteras en ny uppställning för att undersöka biofilmstillväxt i flödesceller med laminärt flöde under längre perioder på upp till flera veckor. Optisk koherenstomografi används för att få tredimensionella skanningar av biofilmstrukturen vid regelbundna intervall. Från dessa tidsserier härleder vi en enkel modell som relaterar tillväxten av en enskild mikrokoloni till tillväxten av hela biofilmen beroende på väggskjuvspänning. Dessutom visar vi att biofilm filament som sträcker sig in i bulkflödet, kan bildas på skarpa biofilmstrukturer i laminärt flöde, även om substratum är en plan yta.  Det andra bidraget är en rapport som beskriver preliminära studier av biofilmsexperiment. Vi undersöker bildandet av biofilm i skjuvskiktet bakom ett bakåtvänt steg. Resultaten indikerar en maximal skjuvspänning, bortom vilken biofilmstillväxt hämmas. / <p>QC 240314</p>

Wall-dominated turbulent dispersed multiphase flows occur in a variety of industrial, biological and environmental applications. The complex nature of the  arrier and the dispersed phase is elevated to a higher level introducing geometrical complexities such as curved walls. Realising such flows and particulate phases poses challenging problems both from computational and also physical point of view. The present thesis tries to address some of these issues Lagrangian computational frame. In the first step, turbulent flow in straight pipes is simulated by means ofdirect numerical simulation with a spectrally accurate code nek5000 to examine the Reynolds number effect on turbulent statistics. Adding the effect of the curvature to these canonical turbulent pipe flows generates Prandtl’s secondary motion of first kind. These configurations, as primary complex geometries in this study, are examined by means of statistical analysis to unfold the evolutionof turbulent characteristics from a straight pipe configuration. A fundamentally different Prandtl’s secondary motion of second kind is also put to test by means of adding the side-walls to a canonical turbulent channel flow and the evolution of various statistical quantities with varying the duct aspect ratios is discussed. After having obtained a characterisation of the turbulent flow in the geometries of bent pipes and ducts, the dispersion of small heavy particles is modelled in the bent pipe by means of point particles which are one-way coupled to the flow. For this purpose a parallel Lagrangian Particle Tracking (LPT) scheme is implemented in the spectral-element code nek5000. Its numerical accuracy, parallel scalability and general performance in realistic situations are scrutinised in various situations. Also, the resulting particle fields are analysed, showing that even a small degree of geometrical curvature has a profound impact on the particle concentration maps. For each of the aforementioned turbulent flow cases new and challenging questions have arisen to be addressed in the present and upcoming research works. Along with an improved understanding of the particle dispersion in the considered complex geometries, the current project is particularly intended to improve the numerical aspects of the current LPT module suitable for largescale computations. / <p>QC 20140226</p>

In this thesis direct numerical simulation is used to investigate the possibilityto delay the transition from laminar to turbulent in boundary layer flows.Furthermore, modal analysis is used to reveal the coherent structures in highdimensional dynamical systems arising in the flow problems.Among different transition scenarios, the classical transition scenario isanalysed. In this scenario, the laminar-turbulent transition occurs when Tollmien-Schlichting waves are triggered inside the boundary layer and grow exponentiallyas they move downstream in the domain. The aim is to attenuate the amplitudeof these waves using active control strategy based on a row of spatiallylocalised sensors and actuators distributed near the wall inside the boundarylayer. To avoid the high dimensional system arises from discretisation of theNavier Stokes equation, a reduced order model (ROM) based on EigensystemRealisation Algorithm (ERA) is obtained and a linear controller is designed.A plasma actuator is modelled and implemented as an external forcing on theflow. To account for the limitation of the plasma actuators and to further reducethe complexity of the controller several control strategies are examinedand compared. The outcomes reveal successful performance in mitigating theenergy of the disturbances inside the boundary layer.To extract coherent features of the wind turbine wakes, modal decompositiontechnique is employed where a large scale dynamical system is reduced toa fewer number of degrees of freedom. Two decomposition techniques are employed:proper orthogonal decomposition and dynamic mode decomposition.In the former procedure, the flow is decomposed into a set of uncorrelated structureswhich are rank according to their energy. In the latter, the eigenvaluesand eigenvectors of the underlying approximate linear operator is computedwhere each mode is associated with a specific frequency and growth rate. Theresults revealed the structures which are dynamically significant to the onsetof instability in the wind turbine wakes. / <p>QC 20130531</p>

In this thesis we apply linear feedback control to spatially evolving flows in order to minimize disturbance growth. The dynamics is assumed to be described by the linearized Navier--Stokes equations. Actuators and sensor are designed and a Kalman filtering technique is used to reconstruct the unknown flow state from noisy measurements. This reconstructed flow state is used to determine the control feedback which is applied to the Navier--Stokes equations through properly designed actuators. Since the control and estimation gains are obtained through an optimization process, and the Navier--Stokes equations typically forms a very high-dimensional system when discretized there is an interest in reducing the complexity of the equations. One possible approach is to perform Fourier decomposition along (almost) homogeneous spatial directions and another is by constructing a reduced order model by Galerkin projection on a suitable set of vectors. The first strategy is used to control the evolution of a range of instabilities in the classical family of Falkner--Skan--Cooke flows whereas the second is applied to a more complex cavity type of geometry. / QC 20101122