<div>This research has been studied a laser-induced spark in flowing gases. The relationship between the minimum ignition energy (MIE), the turbulence intensity, and the flame kernel propagation speed is considered. Plasma emission, produced by the laser-induced spark, and flame kernel generation by the plasma are investigated. The energy balance equation between an ignition energy and energy losses by heat transfer is studied at laminar flows and turbulent flows. Hydrogen and air mixtures were used in a premixed jet burner for ignition experiments. Particle image velocimetry (PIV) examined the velocity and the turbulence intensity under the turbulent flows. The flame kernel development was visualized using Schlieren imaging and infrared images (IR camera). Flame kernel temperatures were measured through Rayleigh scattering and infrared images (IR camera). Plasma evaluations were captured through an intensified CCD camera (ICCD camera). Minimum ignition energies were measured at the laminar flows and the turbulent flows. The MIE decreases with an increase in the turbulence intensity which changed by ignition locations and perforated plates at the constant bulk velocity. Improved mixing rates due to the ignition locations or the geometry of the perforated plates decrease the MIE at the constant bulk velocity. The turbulence intensity increases wrinkles in the flame kernel surface, thus the contact between the flame kernel and reactants increases due to the wrinkles. Therefore, the flame kernel propagation speed increases as the turbulence intensity is higher since the increased reaction by the wrinkles and the contact. Thus, the MIE decreases as the turbulence intensity increases at the constant ignition condition, including bulk velocities and ignition heights, since the high turbulence intensity increases the flame kernel propagation speed. Laser energy differences affect the plasma expansions by the laser absorption. Laser-supported radiation (LSR) wave speeds were measured and calculated using energy balance equations. Velocity does not affect the flame kernel temperature distribution during the early reaction steps because the plasma generates a flame kernel and determines the flame kernel temperature distribution. The MIE increases with increasing the bulk velocity. The energy losses considering convection, conduction, and radiation were calculated using the flame kernel radius, the flame kernel temperature, mixture properties, and the flame speed. The energy balance equation in the ignition of flowing gases is newly written at the laminar flows and the turbulent flows.</div>

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>

Continuous casting is by far the most widely used casting method for steel production all over the world. The Swedish steel industry has specialized in niche products with high quality. To obtain high-quality steel, flow control is crucial when operating the caster. Computational fluid dynamics is used in the research and development of steel production to gain knowledge and optimize flow control, with meshing being a crucial step. One challenge with meshing is finding the balance between computational cost (time and storage) and accuracy. This study investigated whether polyhedral mesh elements would be suffcient for modeling turbulent flow in continuous casting. Firstly, Ansys Fluent's polyhedral converter was used on an existing mesh containing tetrahedral and hexahedral elements. The mesh quality in terms of skewness, orthogonality, and aspect ratio was studied. The results showed that the polyhedral converter lowered the cell count and slightly improved the quality compared to the original mesh. However, more research is needed with other meshes and actual simulations to determine if it is suffcient for modeling turbulent flow in continuous casting. Following that, two new meshes were generated: one polyhedral and one poly-hexcore mesh, both using Ansys Fluent Meshing, and their mesh quality was analyzed as in the first case. These meshes were then simulated with the Large Eddy Simulation model for steel, slag, and air for 60 seconds. They were compared and evaluated using production data in the form of steel layer fluctuations and nail board measurements. The results from the polyhedral and poly-hexcore meshes showed that both meshes had good, equal quality and cell count. Even so, due to the construction of the poly-hexcore, switching between polyhedral and hexahedral elements, and the occurrence of hanging nodes, the simulation showed signs of discretization errors. This might have also contributed to the lower flow velocities at the steel layer compared to the nail boards for the poly-hexcore mesh. The polyhedral mesh showed more stable simulation results with closer velocity trends compared to the nailboard measurements. The steel layer fluctuations for both meshes were somewhat similar to the production data, but further FFT analysis concluded that fluid flow simulations alone are not enough for an FFT comparison to production data since many factors are missing from the simulations, such as mould oscillations and regulations from the stopper. The computational cost varied between time and storage, with the polyhedral mesh requiring 13 % less simulation time but 32.7 % more storage. In summary, it can be concluded that polyhedral mesh elements can be successfully implemented for modeling turbulent flow in continuous castingwith a reasonable computational cost. / Stränggjutning är den främsta metoden använt för stålproduktion över hela världen. Den svenska stålindustrin har specialiserat sig på nischade produkter med hög kvalitet. För att uppnå högkvalitativt stål är flödeskontroll avgörande vid drift av gjutmaskinen. Numeriska strömningsberäkningar används i forskningen och utvecklingen av stålproduktion för att skaffa kunskap och optimera flödeskontrollen, där meshing är ett avgörande steg. En utmaning med meshing är att hitta balansen mellan beräkningskostnad (tid och lagring) och noggrannhet. Denna studie undersökte om polyhedrala mesh-element skulle vara tillräckliga för att modellera turbulent flöde i stränggjutning. Först användes Ansys Fluents polyhedrala konverterare på en befintlig mesh, innehållande tetrahedrala och hexahedrala element. Meshkvalitén i termer av skevhet, ortogonalitet och aspektförhållande studerades. Resultaten visade att den polyhedrala konverteraren sänkte cellantalet och förbättrade kvaliteten något jämfört med den ursprungliga meshen. Dock behövs mer forskning med andra mesher, och faktiska simuleringar, för att avgöra om den är tillräcklig för att modellera turbulent flöde i stränggjutning. Därefter genererades två nya mesher: ett polyhedralt och en poly-hexcore, båda med hjälp av Ansys Fluent Meshing, och precis som i det första fallet, studerades deras meshkvalitet. Dessa två meshes simulerades sedan med Large Eddy Simulation-modellen för stål, slagg och luft i 60 sekunder. De jämfördes och utvärderades med hjälp av produktionsdata i form av stålnivå, fluktuationer och 'nail board'-mätningar. Resultaten från det polyhedrala och poly-hexcore meshen visade att båda mesherna hade god och likvärdiga kvalitéer samt cellantal. Trots detta, på grund av uppbyggnaden av poly-hexcore meshen, alltså växlingen mellan polyhedrala och hexahedrala element och förekomsten av hängande noder, visade simuleringen tecken på diskretiseringsfel. Detta kan också ha bidragit till de lägre flödeshastigheterna vid stålnivån jämfört med 'nail board'-mätningarna för poly-hexcore meshen. Den polyhedrala meshen visade mer stabila simuleringsresultat med hastighetstrender närmare 'nail board'-mätningarna. Amplituden för stålnivåfluktuationerna för båda näten var något lika produktionsdatan, men en vidare FFT-analys visade att endast flödesimuleringar inte är tillräckliga för en FFT-jämförelse med produktionsdata eftersom många faktorer saknas i simuleringarna, såsom svängningar av gjutformen och reglering av inloppashastigheten. Beräkningskostnaden varierade mellan tid och lagring, där den polyhedrala meshen krävde 13 % mindre simuleringstid men 32,7 % mer lagring. Sammanfattningsvis kan det konstateras att polyhedrala mesh-element framgångsrikt kan implementeras för att modellera turbulent flöde i stränggjutning med en rimlig beräkningskostnad.

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.