No abstract provided. / Thesis / Master of Engineering (MEngr) / Scope and contents: This thesis describes the behaviour of iron ore-water suspensions under turbulent flow conditions. This work is divided into two parts. Part I deals with the regimes of transport under steady state flow conditions in circular and horizontal ducts. The heterogeneous flow regime is extensively analyzed; a sequential discrimination of models with an oriented design of experiments have permitted the determination of the best model to correlate hydraulic gradients for these suspensions. A critical discussion on the limit deposit conditions is also included. Part II describes the behaviour of clear water under oscillatory flow conditions. The study demonstrates that the quasi-steady state hypothesis, i.e., fully developed flow assumption, applied to pulsatile turbulent flow under the conditions studied. Observations on the behaviour of iron ore-water suspensions under pulsatile flow are also included. The experiments were carried out using a new air-pulsing technique.

The finite element methods (FEM) have been developed over decades, and together with the growth of computer engineering, they become more and more important in solving large-scale problems in science and industry. The objective of this thesis is to develop high-performance finite element methods (HP-FEM), with two main applications in mind: computational diffusion magnetic resonance imaging (MRI), and simulation of the turbulent flow past a vertical axis wind turbine (VAWT). In the first application, we develop an efficient high-performance finite element framework HP-PUFEM based on a partition of unity finite element method to solve the Bloch-Torrey equation in heterogeneous domains. The proposed framework overcomes the difficulties that the standard approaches have when imposing the microscopic heterogeneity of the biological tissues. We also propose artificial jump conditions at the external boundaries to approximate the pseudo-periodic boundary conditions which allows for the water exchange at the external boundaries for non-periodic meshes. The framework is of a high level simplicity and efficiency that well facilitates parallelization. It can be straightforwardly implemented in different FEM software packages and it is implemented in FEniCS for moderate-scale simulations and in FEniCS-HPC for the large-scale simulations. The framework is validated against reference solutions, and implementation shows a strong parallel scalability. Since such a high-performance simulation framework is still missing in the field, it can become a powerful tool to uncover diffusion in complex biological tissues. In the second application, we develop an ALE-DFS method which combines advanced techniques developed in recent years to simulate turbulence. We apply a General Galerkin (G2) method which is continuous piecewise linear in both time and space, to solve the Navier-Stokes equations for a rotating turbine in an Arbitrary Lagrangian-Eulerian (ALE) framework. This method is enhanced with dual-based a posterior error control and automated mesh adaptation. Turbulent boundary layers are modeled by a slip boundary condition to avoid a full resolution which is impossible even with the most powerful computers available today. The method is validated against experimental data of parked turbines with good agreements. The thesis presents contributions in the form of both numerical methods for high-performance computing frameworks and efficient, tested software, published open source as part of the FEniCS-HPC platform. / <p>QC 20180411</p>

In this study, two major internal erosion initiation processes, suffusion and concentrated leak mechanisms, which lead to both defect formation in a dam’s body and its foundation and high throughflow in dams subjected to internal erosion were studied. This understanding has the potential to facilitate numerical modelling and expedite dam safety assessment studies.  The throughflow properties of coarse rockfill material were studied by; analysing filed pump test data, performing extensive laboratory experiments with a large-scale apparatus and numerically simulating the three-dimensional flow through coarse rock materials, replicating the material used in the laboratory experiments. Results from the tests demonstrate that the parameters of the nonlinear momentum equation of the flow depend on the Reynolds number for pore Reynolds numbers lower than 60000.  Numerical studies were also carried out to conduct numerical experiments. By applying a Lagrangian particle tracking method, a model for estimating the lengths of the flow channels in the porous media was developed.  The shear forces exerted on the coarse particles in the porous media were found to be significantly dependent on the inertial forces of the flow. Suffusion and concentrated leak mechanisms were also studied by means of laboratory experiments to develop a theoretical framework for continuum-based numerical modelling. An erosion apparatus was designed and constructed with the capability of applying hydraulic and mechanical loading. Results were then used to develop constitutive laws of the soil erosion as a function of the applied hydromechanical load for both suffusion and concentrated leak mechanisms. Both the initiation and mass removal rate of were found to be dependent on the soil in-situ stresses. A three-dimensional electrical-resistivity-based tomography method was also adopted for the internal erosion apparatus and was found to be successful in visualising the porosity evolution due to suffusion. / <p>QC 20161006</p>

The doctoral thesis deals with the analysis of the flow in the circular and not circular cross-section pipes by methods using the distribution of the vorticity density. This analysis is particularly focused on the derivation of the new velocity profiles formulas using the above mentioned method. In this work is presented a historical overview of the derived velocity profiles. This overview of already derived velocity profiles will be a fundamental benchmark for newly derived velocity profiles. These new velocity profiles are derived for the circular and not circular cross-section pipes and the derivation is based on the analogy of electromagnetic induction by using Biot-Savart law. It is necessary to apply this analogy at first on solitary vortex filament. By taking this step is possible to get the value of the induced velocity from one solitary vortex filament. Subsequently it is possible to obtain the value of the induced velocity from the vorticity wall and afterwards from the vorticity density distribution over the cross section. This work contains also the results of the experimental measurements of the velocity profiles, and of the CFD simulations. Experimentally measured results are used besides other for the selecting of the most suitable CFD computational model. Selected CFD model will be subsequently declared as a reference model and the valid velocity profiles for this model will serve with the experimentally measured data as a benchmark for the newly derived velocity profiles.

Denna studie har undersökt om förvärmning av luft är en bra metod för att undvika frostbildning i motströmsplattvärmeväxlare i luftbehandlingsaggregat. Frostbildning i värmeväxlare för ventilation är ett problem i kalla klimat och sänker den återvunna energimängden när potentialen för energiåtervinning är som högst.   Tidigare forskning i området har utförts för att bland annat: 1) Definiera säkra förhållanden utan frostbildning i värmeväxlare med laminär luftströmning.  2) Jämföra olika frostkontrollstrategier. 3) Analysera frostbildningens egenskaper och effekt på värmeväxlare i ventilation. Att forskningen inte är entydig för vilken frostkontrollstrategi som är bäst belyser vikten av detta arbete.   I en fallstudie kvantifieras empiriska värmeöverföringskoefficienter som används i en analytisk modell som förutspår energianvändningen för luftvärmning under olika driftfall för ett luftbehandlingsaggregat. Värmeöverföringskoefficienterna tas fram genom mätningar i ett luftbehandlingsaggregat lokaliserat i ett flerbostadshus. Mätobjektet har ett installerat värmebatteri före värmeväxlaren vilket värmer inkommande uteluft till konstant lufttemperatur.  Tillverkaren proklamerar att det sker turbulent luftströmning i den studerade värmeväxlaren. Vid olika lufthastigheter ändras värmeöverföringskoefficienten mer vid turbulent strömning än vid laminär strömning. Att olika lufthastigheter har betydlig inverkan på värmeöverföringskoefficienten överensstämmer med resultatet från mätningarna. Effektiviteten av den aktuella styrningen av förvärmaren har granskats och brister har upptäckts. En ny optimerad styrning föreslås för frostfri drift och minimerad energianvändning. Den optimerade styrningen jämförs sedan med en vanlig avfrostningsmetod och utförs med hjälp av flödesbalansstatistik från ett aggregat med sektionsavfrostningsstyrning. Resultatet visar att förvärmning av luft för frostfri drift använder mindre energi än vid sektionsavfrostningsstyrning.  Den föreslagna styrningen innebär att bypass-spjället öppnas då förvärmaren går på full effekt. Detta för att förhindra påfrostning vid extrema utetemperaturer, då förvärmaren inte ger tillräcklig effekt för att undvika påfrostning. Den föreslagna styrningen regleras utifrån daggpunkten, temperaturen på tilluften innan värmeväxlaren och temperaturen på avluften. Den framtagna värmeöverföringskoefficienten i studien tar inte i beaktning ökad koefficient under fuktiga förhållanden då daggpunkten i frånluften är över fryspunkten och det sker kondens. En framtida studie kan hitta värmeöverföringskoefficienter som förutspår yttemperaturen under kondensutfällning då daggpunkten är över 0°C vilket ger medel för att ta fram en ännu mera energieffektiv styrning. Detta kan utföras genom mätning i ett luftbehandlingsaggregat med hög temperaturverkningsgrad placerat i en fuktig miljö. / This study has analysed if preheating of air is a good method to prevent frost formation in a counter flow plate heat exchanger used in Air Handling Units. Frosting in heat exchangers used in AHU-systems is a problem in cold climates and lowers the energy recovery when its potential is the highest.   Previous research has been done to 1) Define safe operating conditions without frost formation in heat exchangers with laminar flow. 2) Compare different frost control strategies. 3) Analyse frost formation properties and its effects on heat exchangers in AHU-systems. That the research is not unambiguous for which frost control method is the best highlights the importance of this work. Empirical heat transfer coefficients are quantified in this study and are used in the creation of an analytical model of a counter flow heat exchanger that predicts the energy use for heating the air at different operating modes. The heat transfer coefficients are produced by measurements in an AHU-systems located in an apartment building. The AHU has an installed air heater before the heat exchanger which heats the outdoor inlet air to constant temperature.   The manufacturer of the AHU-system proclaims turbulent flow in the heat exchanger. Change of airflow has greater impact on the heat transfer coefficients during turbulent flow compared to laminar flow conditions. This is also derived from the results of the measurements. The effectiveness of the current control system is analyzed, and flaws are discovered. A new control system is therefore proposed for frost free operation and higher energy efficiency. The proposed control system is compared to a frequently used frost control system which uses bypass-dampers to redirect the inlet air to the reheater, to let the warm air stream melt formed frost on the warm air side. This comparison is accomplished by analyzing the usage of the bypass-dampers during different outdoor temperatures in a comparable AHU-system that have the mentioned frost control method. The results show that preheating of inlet air for frost free operation uses less energy than usage of the bypass-dampers to melt formed frost on the warm air side. The proposed control system needs the bypass-dampers to be used when the preheater is on full operation to prevent frost formation at extreme outdoor temperatures when the preheating may not be enough. The proposed control system is regulated by the dewpoint and the temperatures of the exhaust air stream and the supply air stream. The derived heat transfer coefficient that is used to predict the plate-temperature take no account for condensation on the warm air side during humid conditions when the dew point is above the freezing point. Future studies can derive the heat transfer coefficient during condensation which will improve the prediction of the plate-temperature. This would make the system more energy efficient during humid air conditions. This can be done by measurements in an AHU-system with high temperature efficiency placed in a humid environment.

O concreto é um dos materiais de construção civil mais versáteis, podendo se adaptar a formas diversas quando em seu estado fresco e resistindo a grandes cargas de compressão em seu estado rígido. No entanto, em estruturas mais densamente armadas ou com geometrias mais complexas, pode-se apresentar dificuldades para a moldagem, causando falhas no preenchimento da forma, o que reduz a capacidade resistente da peça e sua vida útil. Neste trabalho foi utilizado o método de partículas lagrangeanas Moving Particle Semi-Implicit (MPS) como base para um simulador para o estudo do comportamento do escoamento de pastas e argamassas cujo comportamento pode se aproximado por modelos reológicos como Bingham e Herschel-Bulkley. Foram propostos módulos para a simulação da viscosidade não-newtoniana, variação térmica no processo de cura e modelagem de turbulência. Para modelar a variação de viscosidade de um fluido não-newtoniano foi utilizado o modelo de Herschel-Bulkley. O modelo de Herschel-Bulkley possui uma singularidade para taxas de deformação muito pequenas, que resulta em valores de viscosidade infinitas, dificuldade contornada pela solução proposta por Papanastasiou (1987). Na modelagem térmica foram analisados dois modelos de dissipação, sendo um original do método e outro calculado por meio do divergente do gradiente utilizando os modelos de partículas e que teria melhores resultados para o cálculo da dissipação térmica. Também foi modelada a convecção térmica, utilizando para isso a hipótese de Boussinesq. A reação de hidratação do concreto foi modelada utilizando uma equação do tipo Hill para representar a elevação de temperatura obtida por meio um ensaio adiabático. Para complementar as simulações, foi utilizado o modelo de turbulência Detached Eddy Simulation (DES), baseado no método Large Eddy Simulation (LES), que utiliza um modelo de parede para simular a interação do fluido. Para a implementação deste modelo de turbulência foi proposto um algoritmo para o cálculo da distância da partícula de fluido à parede. Este algoritmo utiliza estruturas de dados já existentes no método de partículas de modo que sua execução requer muito menos recursos que a busca binária. Apesar do trabalho ter se guiado pela simulação do concreto fresco, que é um material particularmente complexo, outros escoamentos encontrados na engenharia civil também podem ser beneficiados pelo método, como os estudos do escoamento em sistemas prediais de água e esgoto, do escoamento e prevenção de erosão em rios e córregos, do escorregamento de encostas, dos reatores para depuração de esgotos, entre outros. / The concrete is one of the more versatile civil construction materials, it can adapt to various forms when in its fresh state while resisting to high compression loads in its rigid state. However, in some cases like densely reinforced concrete structures and complex geometry concrete structures can present issues to the casting, and failure in proper form filling can occur. These failures can reduce the resistance and the lifetime of the structure. This work used a simulator based on the lagrangean particle method called Moving Particle Semi-Implicit (MPS) to study the concrete behavior in its distinctive characteristics. Also, this work proposed modules to simulate the non-Newtonian viscosity, the thermal process of concrete cure and to model the turbulent flow. To model the variation of viscosity of a non-Newtonian fluid, the Herschel-Bulkley model, which relates the shear stress with the strain rate, was applied. The Herschel-Bulkley model has a singularity at low strain rates, which results in infinite viscosities. To overcome this issue, Papanastasiou (1987) proposed a modification in the model in order to eliminate the singularity. For the thermal modeling, two models for thermal dissipation were compared, the original method and other calculated from the divergence of gradient using the differential operators for the particle model and that could present improved results for the thermal dissipation calculation. Also, the thermal convection was modeled using the Boussinesq hypothesis. The hydration reaction of the concrete was modeled using a Hill type equation to reproduce the temperature elevation. In addition, a Detached Eddy Simulation (DES) based turbulence model with a simple wall model in the interaction of wall and fluid was applied in the simulations. To improve the turbulence model, an algorithm to calculate the distance between fluid and the nearest wall particle was proposed. The algorithm uses already calculated information from particles so that the execution requires much less resources than a binary search. Although the work has been focused on to the modeling of fresh concrete simulation, which is a particularly complex material, other flows found in civil engineering can also be benefited by the method, such as studies of drainage in water and sewage systems, drainage and prevention of erosion into rivers and streams, prevention and damage mitigation of landslides, reactors for sewage treatment among many others.

O concreto é um dos materiais de construção civil mais versáteis, podendo se adaptar a formas diversas quando em seu estado fresco e resistindo a grandes cargas de compressão em seu estado rígido. No entanto, em estruturas mais densamente armadas ou com geometrias mais complexas, pode-se apresentar dificuldades para a moldagem, causando falhas no preenchimento da forma, o que reduz a capacidade resistente da peça e sua vida útil. Neste trabalho foi utilizado o método de partículas lagrangeanas Moving Particle Semi-Implicit (MPS) como base para um simulador para o estudo do comportamento do escoamento de pastas e argamassas cujo comportamento pode se aproximado por modelos reológicos como Bingham e Herschel-Bulkley. Foram propostos módulos para a simulação da viscosidade não-newtoniana, variação térmica no processo de cura e modelagem de turbulência. Para modelar a variação de viscosidade de um fluido não-newtoniano foi utilizado o modelo de Herschel-Bulkley. O modelo de Herschel-Bulkley possui uma singularidade para taxas de deformação muito pequenas, que resulta em valores de viscosidade infinitas, dificuldade contornada pela solução proposta por Papanastasiou (1987). Na modelagem térmica foram analisados dois modelos de dissipação, sendo um original do método e outro calculado por meio do divergente do gradiente utilizando os modelos de partículas e que teria melhores resultados para o cálculo da dissipação térmica. Também foi modelada a convecção térmica, utilizando para isso a hipótese de Boussinesq. A reação de hidratação do concreto foi modelada utilizando uma equação do tipo Hill para representar a elevação de temperatura obtida por meio um ensaio adiabático. Para complementar as simulações, foi utilizado o modelo de turbulência Detached Eddy Simulation (DES), baseado no método Large Eddy Simulation (LES), que utiliza um modelo de parede para simular a interação do fluido. Para a implementação deste modelo de turbulência foi proposto um algoritmo para o cálculo da distância da partícula de fluido à parede. Este algoritmo utiliza estruturas de dados já existentes no método de partículas de modo que sua execução requer muito menos recursos que a busca binária. Apesar do trabalho ter se guiado pela simulação do concreto fresco, que é um material particularmente complexo, outros escoamentos encontrados na engenharia civil também podem ser beneficiados pelo método, como os estudos do escoamento em sistemas prediais de água e esgoto, do escoamento e prevenção de erosão em rios e córregos, do escorregamento de encostas, dos reatores para depuração de esgotos, entre outros. / The concrete is one of the more versatile civil construction materials, it can adapt to various forms when in its fresh state while resisting to high compression loads in its rigid state. However, in some cases like densely reinforced concrete structures and complex geometry concrete structures can present issues to the casting, and failure in proper form filling can occur. These failures can reduce the resistance and the lifetime of the structure. This work used a simulator based on the lagrangean particle method called Moving Particle Semi-Implicit (MPS) to study the concrete behavior in its distinctive characteristics. Also, this work proposed modules to simulate the non-Newtonian viscosity, the thermal process of concrete cure and to model the turbulent flow. To model the variation of viscosity of a non-Newtonian fluid, the Herschel-Bulkley model, which relates the shear stress with the strain rate, was applied. The Herschel-Bulkley model has a singularity at low strain rates, which results in infinite viscosities. To overcome this issue, Papanastasiou (1987) proposed a modification in the model in order to eliminate the singularity. For the thermal modeling, two models for thermal dissipation were compared, the original method and other calculated from the divergence of gradient using the differential operators for the particle model and that could present improved results for the thermal dissipation calculation. Also, the thermal convection was modeled using the Boussinesq hypothesis. The hydration reaction of the concrete was modeled using a Hill type equation to reproduce the temperature elevation. In addition, a Detached Eddy Simulation (DES) based turbulence model with a simple wall model in the interaction of wall and fluid was applied in the simulations. To improve the turbulence model, an algorithm to calculate the distance between fluid and the nearest wall particle was proposed. The algorithm uses already calculated information from particles so that the execution requires much less resources than a binary search. Although the work has been focused on to the modeling of fresh concrete simulation, which is a particularly complex material, other flows found in civil engineering can also be benefited by the method, such as studies of drainage in water and sewage systems, drainage and prevention of erosion into rivers and streams, prevention and damage mitigation of landslides, reactors for sewage treatment among many others.

La prédiction des vibrations induites par un écoulement est essentielle dans la conception des conduits de nombreuses installations industrielles, en particulier dans l’industrie du gaz. Notre étude concerne la prévision du bruit et la vibration des conduits soumis à un écoulement turbulent à faible nombre de Mach. Notre objectif est de présenter une étude numérique et expérimentale permettant aux ingénieurs de mieux comprendre le couplage entre l’excitation aléatoire et le conduit pour deux géométries (circulaire ou rectangulaire). Une approche expérimentale est développée et utilisée pour valider les prévisions numériques. Deux cas sont étudiés : (i) un conduit droit sans singularité, où les modes acoustiques du conduit sont excités par une couche limite turbulente (TBL) et (ii) un conduit droit avec un diaphragme inséré en amont qui génère une source acoustique localisée. La contribution acoustique est déterminée soit par des méthodes de mesure d’interspectres, soit à l’aide des outils de mécanique des fluides numérique (CFD) et d’analogies aéroacoustiques. La réponse de la structure est estimée par une approche dite de « couplage faible » qui utilise des fonctions de transfert modale d’un conduit fini simplement appuyé. Les mesures conduiront à évaluer et suggérer des améliorations de modèles empiriques existants de densité interspectrale de puissance (CPSD) dans un contexte d’écoulements internes turbulents. Une analyse modale expérimentale d’un conduit rectangulaire finie est confrontée à des méthodes de calcul pour évaluer l’effet des conditions aux limites, du rayonnement acoustique et de l’amortissement aérodynamique. Le couplage fluide structure est analysé par la fonction de « joint acceptance » à la fois dans le domaine spatial et dans le domaine des nombres d’onde. L’excitation comprend à la fois les contributions acoustiques et hydrodynamiques à l’aide des CPSD exprimées sur la base des fonctions de cohérence de type Corcos, champ diffus et modes acoustiques d’ordre élevé. Enfin, les études numériques et expérimentales de cette thèse ont été utilisées pour développer un cadre d’étude et de modélisation du bruit et des vibrations dans les conduites, qui relie la dynamique des fluides, les modèles analytiques et empiriques à des techniques efficaces d’analyse aléatoire. / Pipeline and duct vibrations can cause a range of issues from unplanned shutdownsto decreased equipment life time. Thus, the prediction of flow-induced vibrations is essential in piping design in many industrial plants, especially, for Gas industry. This study deals with the prediction of pipe flow noise and vibration at low Mach number. We aim to present a numerical and experimental study which can offer engineers a better understanding of the coupling between random excitation and duct section for two geometries (circular or rectangular). An experimental facility and measurement approach is developed and used to validate numerical predictions. Two cases are investigated: (i) a straight duct with no singularity, duct acoustic modes are excited by the Turbulent Boundary Layer (TBL) and (ii) a straight duct with a diaphragm inserted upstream generating a localized acoustic source. The acoustic contribution is either measured via cross-spectra based methods or calculated using Computational Fluid Dynamics (CFD) and aeroacoustic analogies. The response of the structure is estimated via a ‘blocked’ approach using analytical modal Frequency Response Functions (FRFs) of a simply supported finite duct. Measurements will lead to evaluate and suggest improvements to existing Cross Power Spectral Density (CPSD) empirical models in a context of internal turbulent flows. Experimental modalanalysis of a finite rectangular duct are confronted to computational methods to assess the effect of the Boundary Conditions (BCs), the resistive damping from coupling with the internal acoustic medium and aerodynamic damping. The fluid-structure coupling is analyzed through the joint acceptance function both in the spatial and wave number domain. The excitation includes both the acoustic and hydrodynamic contributions using CPSD written on the basis of Corcos, Diffuse Acoustic Field (DAF) and acoustic duct mode coherence functions. Finally, the numerical and experimental studies in this thesis were used to develop a framework for studying and modelling pipe flow noise and vibration which links CFD, analytical and empirical models to efficient random analysis techniques.

lthough extensive research work has been carried out on the drag reduction behavior of polymers and surfactants alone, little progress has been made on the synergistic effects of combined polymers and surfactants. A number of studies have demonstrated that certain types of polymers and surfactants interact with each other to form surfactant-polymer complexes. The formation of such complexes can cause changes in the solution properties and may result in better drag reduction characteristics as compared with pure additives. A series of drag-reducing surfactants and polymers were screened for the synergistic studies. The following two widely used polymeric drag reducing agents (DRA) were chosen: a copolymer of acrylamide and sodium acrylate (referred to as PAM) and polyethylene oxide (PEO). Among the different types of surfactants screened, a cationic surfactant octadecyltrimethylammonium chloride (OTAC) and an anionic surfactant Sodium dodecyl sulfate (SDS) were selected for the synergistic study. In the case of the cationic surfactant OTAC, sodium salicylate (NaSal) was used as a counterion. No counterion was used with anionic surfactant SDS. The physical properties such as viscosity, surface tension and electrical conductivity were measured in order to detect any interaction between the polymer and the surfactant. The drag reduction (DR) ability of both pure and mixed additives was investigated in a pipeline flow loop. The effects of different parameters such as additive concentration, type of water (deionized (DI) or tap), temperature, tube diameter, and mechanical degradation were investigated. The addition of OTAC to PAM solution has a significant effect on the properties of the system. The critical micelle concentration (CMC) of the mixed surfactant-polymer system is found to be different from that of the surfactant alone. The anionic PAM chains collapse upon the addition of cationic OTAC and a substantial decrease in the viscosity occurs. The pipeline flow behaviour of PAM/OTAC mixtures is found to be consistent with the bench scale results. The drag reduction ability of PAM is reduced upon the addition of OTAC. At low concentrations of PAM, the effect of OTAC on the drag reduction behavior is more pronounced. The drag reduction behavior of polymer solutions is strongly influenced by the nature of water (de-ionized or tap). The addition of OTAC to PEO solution exhibited a week interaction based on the viscosity and surface tension measurements. However, the pipeline results showed a considerable synergistic effect, that is, the mixed system gave a significantly higher drag reduction (lower friction factors) as compared with the pure additives (pure polymer or pure surfactant). The synergistic effect in the mixed system was stronger at low polymer concentrations and high surfactant concentrations. Also the resistance against mechanical degradation of the additive was improved upon the addition of OTAC to PEO. The mixed PEO/SDS system exhibited a strong interaction between the polymers (PEO) and the surfactant (SDS), Using electrical conductivity and surface tension measurements, the critical aggregation concentration (CAC) and the polymer saturation point (PSP) were determined. As the PEO concentration is increased, the CAC decreases and the PSP increase. The addition of SDS to the PEO solution exhibits a remarkable increase in the relative viscosity compared to the pure PEO solution. This increase is attributed to the changes in the hydrodynamic radius of the polymer coil. The pipeline flow exhibited a considerable increase in DR for the mixed system as compared to the pure PEO solution. The addition of surfactant always improves the extent of DR up to the PSP. Also the mixed PEO/ SDS system shows better resistance against shear degradation of the additive.

Large-eddy simulation (LES) and unsteady Reynolds-averaged Navier-Stokes (URANS) calculations have been performed to investigate the effects of different mathematical models for scalar variance and its dissipation rate as applied to both a non-reacting bluff-body turbulent flow and an extension to a reacting case. In the conserved scalar formalism, the mean value of a thermo-chemical variable is obtained through the PDF-weighted integration of the local description over the conserved scalar, the mixture fraction. The scalar variance, one of the key parameters for the determination of a presumed β-function PDF, is obtained by solving its own transport equation with the unclosed scalar dissipation rate modelled using either an algebraic expression or a transport equation. The proposed approach is first applied to URANS and then extended to LES. Velocity, length and time scales associated with the URANS modelling are determined using the standard two-equation k-ε transport model. In contrast, all three scales required by the LES modelling are based on the Smagorinsky subgrid scale (SGS) algebraic model. The present study proposes a new algebraic and a new transport LES model for the scalar dissipation rate required by the transport equation for scalar variance, with a time scale consistent with the Smagorinsky SGS model.