In the last 10-15 years, large-eddy simulation (LES) has become well established for non-reacting flows, and several successful models have been developed for the transfer of momentum and kinetic energy to the subfilter-scales (SFS). However, for reacting flows, LES is still undergoing significant development. In particular, for many premixed combustion applications, the chemical reactions are confined to propagating surfaces that are significantly thinner than the computational grids used in practical LES. In these situations, the chemical kinetics and its interaction with the turbulence are not resolved and must be entirely modelled. There is, therefore, a need for accurate and robust physical modelling of combustion at the subfilter-scales. In this thesis, modelled transport equations for progress variable and flame surface density (FSD) were implemented and coupled to the Favre-filtered Navier-Stokes equations for a compressible reactive thermally perfect mixture. In order to reduce the computational costs and increase the resolution of simulating combusting flows, a parallel adaptive mesh (AMR) refinement finite-volume algorithm was extended and used for the prediction of turbulent premixed flames. The proposed LES methodology was applied to the numerical solution of freely propagating flames in decaying isotropic turbulent flow and Bunsen-type flames. Results for both stoichiometric and lean flames are presented. Comparisons are made between turbulent flame structure predictions for methane, propane, hydrogen fuels, and other available numerical results and experimental data. Details of subfilter-scale modelling, numerical solution scheme, computational results, and capabilities of the methodology for predicting premixed combustion processes are included in the discussions. The current study represents the first application of a full transport equation model for the FSD to LES of a laboratory-scale turbulent premixed flame. The comparisons of the LES results of this thesis to the experimental data provide strong support for the validity of the modelled transport equation for the FSD. While the LES predictions of turbulent burning rate are seemingly correct for flames lying within the wrinkled and corrugated flamelet regimes and for lower turbulence intensities, the findings cast doubt on the validity of the flamelet approximation for flames within the thin reaction zones regime.

This thesis presents a method for optimizing cable configuration inside a large magnetic cylindrical steel casing, from the total ampacity point of view. The method is comprised of two main parts, namely: 1) analytically calculating the electromagnetic losses in the steel casing and sheathed cables, for an arbitrary cables configuration, and 2) implementing an algorithm for determining the optimal cables configuration to obtain the best total ampacity. The first part involves approximating the eddy current and hysteresis losses in the casing and cables. The calculation is based on the theory of images, which this thesis expands to apply to casings having both high magnetic permeability and high electric conductivity at the same time. The method of images, in combination with approximating the cable conductors and sheaths as multiple physical filaments, is used to compute the final current distributions in the cables and pipe and thus the associated losses. The accuracy of this computation is assessed against numerical solutions obtained using the Maxwell finite element program by Ansoft. Next, the optimal cable configuration is determined by applying a proposed two-level optimization algorithm. At the outer level, a combinatorial optimization based on a genetic algorithm explores the different possible configurations. The performance of every configuration is evaluated according to its total ampacity, which is calculated using a convex optimization algorithm. The convex optimization algorithm, which forms the inner level of the overall optimization procedure, is based on the barrier method. This proposed optimization procedure is tested for a duct bank installation containing twelve cables and fifteen ducts, comprising two circuits and two cables per phase, and compared with a brute force method of considering all possible configurations. The optimization process is also applied to an installation consisting of a single circuit inside a large magnetic steel casing.

This study investigates the complex flow structure generated by riparian emergent vegetation along the edge of floodplain. Detailed velocity and boundary shear stress measurements were carried out for various arrangements of emergent rigid cylindric rods of 3 mm, 6 mm and 9 mm diameters and for three different rod densities. In addition, the impact of foliage on the flow field was assessed during a series of experiments where brushes were used instead of smooth rods. The results of these new experiments are first presented. In addition to the laboratory data, field data was obtained through Acoustic Doppler Current Profiler measurements for two flood events in a stretch of the river Rhône that can be approximated to a straight compound channel with vegetated banks. The analysis of the flow structure highlights the presence of strong secondary circulation and increased vorticity on the river banks. The rods on the edge of the floodplain increase significantly flow resistance, reducing velocity and decreasing boundary shear stress. Flow rate was seen to decrease with increasing vegetative density for all cases except when foliage was added. This suggests that an optimum threshold density, for which a smaller density would lead to an increased flow rate might exist. Wakes trailing downstream of the vegetation stem, planform coherent structures advected between the main channel and the floodplain, and eddying motion in the flow due to enhanced turbulence anisotropy are among the defining patterns observed in the studied compound channel flows with one line of emergent vegetation along the edge of the floodplain. The Shiono and Knight Method (SKM) was modified in order to account for the increased turbulence activity due to the rods. The drag force term was introduced in the same way as in the work of Rameshwaran and Shiono (2007). However, a new term was added to the transverse shear stress term in the form of an Elder formulation, incorporating a friction drag coefficient which can be derived from the experimental data. In this proposed version, the advection term was set to zero. Another version of the SKM, similar to Rameshwaran and Shiono (2007), was also tested with the addition of a local drag friction only applied in the rod region. The proposed SKM version without the advection term was favored as it can be more closely related to the experimental data and to physical processes. Finally, the capabilities of Telemac-2D were tested against the experimental data for various turbulence models. The Large Eddy Simulation turbulence model highlighted some unsteady flow patterns that were observed during experiments, while satisfactorily predicting the lateral velocity and boundary shear stress distributions.

A large eddy simulation (LES) model has been developed and validated for turbulent non-premixed and partially premixed combustion systems. LES based combustion modelling strategy has the ability to capture the detailed structure of turbulent flames and account for the effects of radiation heat loss. Effects of radiation heat loss is modelled by employing an enthalpy-defect based non-adiabatic flamelet model (NAFM) in conjunction with a steady non-adiabatic flamelet approach. The steady laminar flamelet model (SLFM) is used with multiple flamelet solutions through the development of pre-integrated look up tables. The performance of the non-adiabatic model is assessed against experimental measurements of turbulent CH4/H2 bluff-body stabilized and swirl stabilized jet flames carried out by the University of Sydney combustion group. Significant enhancements in the predictions of mean thermal structure have been observed with both bluff body and swirl stabilized flames by the consideration of radiation heat loss through the non-adiabatic flamelet model. In particular, mass fractions of product species like CO2 and H2O have been improved with the consideration of radiation heat loss. From the Sydney University data the HM3e flame was also investigated with SLFM using multiple flamelet strategy and reasonably fair amount of success has been achieved. In this work, unsteady flamelet/progress variable (UFPV) approach based combustion model which has the potential to describe both non-premixed and partially premixed combustion, has been developed and incorporated in an in-house LES code. The probability density function (PDF) for reaction progress variable and scalar dissipation rate is assumed to follow a delta distribution while mixture fraction takes the shape of a beta PDF. The performance of the developed model in predicting the thermal structure of a partially premixed lifted turbulent jet flame in vitiated co-flow has been evaluated. The UFPV model has been found to successfully predict the flame lift-off, in contrast SLFM results in a false attached flame. The mean lift-off height is however over-predicted by UFPV-δ function model by ~20% for methane based flame and under-predicted by ~50% for hydrogen based flame. The form of the PDF for the reaction progress variable and inclusion of a scalar dissipation rate thus seems to have a strong influence on the predictions of gross characteristics of the flame. Inclusion of scalar dissipation rate in the calculations appears to be successful in predicting the flame extinction and re-ignition phenomena. The beta PDF distribution for the reaction progress variable would be a true prospect for extending the current simulation to predict the flame characteristics to a higher degree.

Methane (CH4) and carbon dioxide (CO2) are two important greenhouse gases. Recent studies have shown that lakes, although they cover a small area of the globe, can be very important natural sources of atmospheric CH4 and CO2. It is therefore important to monitor the fluxes of these gases between lakes and the atmosphere in order to understand the processes that govern the exchange. By using the eddy covariance method for lake flux studies, the resolution in time and in space of the fluxes is increased, which gives more information on the governing processes. Eddy covariance measurements at a Swedish lake revealed a diel cycle in the fluxes of both CH4 and CO2, with higher fluxes during nighttime than daytime. The high nighttime CO2 fluxes could to a large extent be explained with enhanced transfer velocities due to waterside convection. For the diel cycle of CH4 flux it was suggested that waterside convection could enhance the transfer velocity, transport CH4 rich water to the surface, as well as trigger ebullition. Simultaneous flux measurements of CH4 and CO2 have been presented using both the eddy covariance method and the floating chambers method of which the latter is the traditional measuring method for lake fluxes. For CO2 the two methods agreed well during some periods but differed considerably during others. Disagreement between the methods might be due to horizontal heterogeneity in partial pressure of CO2 in the lake. The methods agreed better for the CH4 flux measurements. However, it is clear that due to the discontinuous nature of the floating chambers, this method will likely miss important high flux events. The main conclusions of this thesis are: 1) the two gas flux methods are not directly comparable and should be seen as supplementary to each other 2) waterside convection enhances the fluxes of both CH4 and CO2 over the water-air surface. If gas flux measurements are not conducted during nighttime, potential high flux periods might be missed and estimates of the total amount of gas released from lakes to the atmosphere may be biased.

The periods of time-varying turbulence in the atmospheric boundary layer, i.e.\ the morning and evening transitions, are often overlooked or highly idealised by dispersion models. These transitions make up a significant portion of the diurnal cycle and are known to affect the spread of pollution due to the properties of turbulence in the residual and stable layers, resulting in phenomena such as lofting, trapping, and fumigation.\\ Two main simulation techniques are presented for the purpose of modelling the dispersion of passive tracers in both convective and evening transitional boundary layers: Lagrangian stochastic (LS) modelling for 1D, inhomogeneous, non-stationary turbulence; and large-eddy simulation (LES) with a particle model tracing pollutant paths using a combination of the resolved flow velocities and a random displacement model to represent sub-grid scale motions.\\ In the convective boundary layer, LS models more accurately representing the state of turbulence, and including the effect of skewness, are shown to produce dispersion results in closer agreement with LES. By considering individual particle trajectories, a reflective top boundary in LS models is shown to produce un-physical, sharp changes in velocity and position. By applying a correction to the vertical velocity variance based on representing the stable potential temperature gradient above the boundary layer, particles are contained within the boundary layer in a physically accurate way. \\ An LS model for predicting dispersion in time-varying, skewed turbulence is developed and tested for various particle releases in transitional boundary layers with different rates of decay, showing an improvement in accuracy compared with previous LS models. Further improvement is made by applying a correction to the vertical velocity variance to represent the effect of a positive potential temperature gradient developing over the course of the transition. Finally, a developing stable boundary layer is shown to have a significant trapping effect on particles released near the surface. \\

La recherche en bruit de combustion est de nos jours majoritairement consacrée au développement d'outils de calcul du bruit rayonné par les flammes. Les méthodes actuelles de CFD telles que la LES ou la DNS sont capables de fournir le champ acoustique rayonné par des sources de bruit, mais elles sont cependant limitées à des domaines de faible taille, ceci dû à leur fort coût de calcul. Pour surmonter cette limitation, on a vu l'émergence de méthodes hybrides. Dans cette approche, les sources de bruit sont découplées du son rayonné. Les sources sont alors calculées par les méthodes de DNS et de LES tandis que l'analogie acoustique permet de calculer le son rayonné par des codes acoustiques, moins coûteux en temps de calcul.L'objet de cette étude est de développer un outil numérique sur la base de l'analogie acoustique de Phillips pour de faibles nombres de Mach. Il permet de prendre en compte l'impact des conditions limites sur le champ acoustique résultant. La LES et le code de calcul acoustique développé ont été utilisés pour évaluer le bruit produit par une flamme turbulente confinée. Les deux techniques donnent des résultats en accord tant que les bonnes quantités sont comparées: il a été observé que le signal de pression obtenu directement à partir de la LES contient une quantité non négligeable d'hydrodynamique, laquelle doit être négligée si on cherche à comparer seulement les champs acoustiques issus des deux codes.L'hypothèse d'un nombre de Mach faible est totalement réaliste si l'on considère l'écoulement présent dans une chambre de combustion. Elle conduit à des simplifications significatives lorsque les analogies acoustiques sont considérées. Cependant, cette hypothèse ne peut pas être utilisée pour l'écoulement en amont (entrée d'air, compresseur) ni en aval (turbine, tuyère) des chambres de combustion aéronautiques. Un outil numérique a été développé pour pallier ce problème. Il est basé sur les équations d'Euler Quasi-1D, qui prennent en compte des écoulements convectifs, non isentropiques et non isenthalpiques. Cet outil permet d'estimer les conditions limites acoustiques qui doivent être imposées sur les entrées/sorties d'une chambre de combustion pour prendre en compte la présence d'un écoulement de nombre de Mach non négligeable, alors que les calculs acoustiques sont eux effectués sous cette hypothèse fortement restrictive. / Today, much of the current effort in combustion noise is the development of efficient numerical tools to calculate the noise radiated by flames. Although unsteady CFD methods such as LES or DNS can directly provide the acoustic field radiated by noise sources, this evaluation is limited to small domains due to high computational costs. Hybrid methods have been developed to overcome this limitation. In these schemes, the noise sources are decoupled from the radiated sound. The sources are still calculated by DNS or LES codes whereas the radiated sound is evaluated by acoustic codes using an acoustic analogy.In the present study, a numerical tool based on the Phllips' analogy for low Mach numbers flows has been developed. This tool accounts for the role of the boundary conditions in the resulting acoustic field. Both LES and the acoustic code developed here are used to assess the noise produced by a turbulent confined flame of a turbulent swirled-stabilized staged combustor. Good agreements are obtained between both techniques as long as the good quantities are compared: the pressure signal obtained directly from LES contains a non negligible amount of hydrodynamics that must be removed when a suitable acoustics-acoustics comparison is sought. The low Mach number assumption is completely realistic when considering the flow within a combustion chamber; it also conducts to considerable simplifications when leading with acoustic analogies. However, it cannot be used for the upstream (air-intake, compressors) and downstream (turbines, nozzle) of an aeronautical combustion chamber. A numerical tool is developed based on the quasi-1D Linearized Euler Equations in order to account for convective, non-isentropic and non-isenthalpic flows. By means of this tool, it is possible to estimate the acoustic boundary conditions that should be imposed at the inlet/oultlet of a given combustion chamber when performing low-Mach number acoustic computations.

Le vieillissement du parc aéronautique, mais aussi la volonté de prolonger le cycle de vie des appareils, impose aux techniques de maintenance des exigences de contrôle accrues en termes de fiabilité et de rapidité. Les principaux enjeux étant la détection, et surtout la caractérisation des micro-fissures pouvant apparaitre dans les pièces sensibles des appareils. Cette caractérisation pouvant aller jusqu'à la reconstruction qualitative, voire quantitative du profil des fissures, ce qui suppose la résolution du problème inverse consistant, à partir des signaux mesurés, à remonter à certaines caractéristiques de la pièce inspectée, notamment géométriques, qui en sont en partie la cause.Ce manuscrit présente une approche originale de conception de sonde d'imagerie par courants de Foucault. Cette approche, ici appliquée à l'imagerie de pièces aéronautiques cylindriques, consiste à concevoir la sonde de façon à satisfaire la double contrainte d'une instrumentation performante et d'une résolution possible du problème inverse. Ainsi, la conception de la sonde est-elle faite de manière à ce que d'une part les signaux mesurés, résultant de ses interactions avec la pièce contrôlée, soient d'amplitude la plus élevée possible, et d'autre part, a n que ces interactions soient modélisables au moyen d'un modèle qui se prête bien à l'inversion.Dans un premier temps, une méthodologie de conception de sonde est présentée, qui vise à optimiser le dimensionnement et l'agencement de ses éléments en se fixant des contraintes de sensibilité et de résolution. Un prototype, conçu d'après cette méthodologie, est réalisé, mis en œuvre et caractérisé. Les résultats expérimentaux obtenus, qui s'avèrent en accord avec une modélisation par éléments finis, offrent une validation du fonctionnement de la sonde, celle-ci permettant effectivement de mettre en évidence les défauts types (micro-fissures) recherchés dans la pièce inspectée. Dans un second temps, la configuration des interactions sonde-pièce inspectée, choisie à dessein, est mise à profit pour élaborer un modèle direct. Celui-ci repose sur l'hypothèse selon laquelle, les défauts étant de faibles dimensions, il est possible de considérer que la perturbation de signal qu'ils engendrent est équivalente à celle que génèreraient des sources virtuelles de courants localisées uniquement dans le volume du défaut. Outre sa simplicité, cette modélisation offre l'avantage d'être adaptée à une mise en œuvre au moyen de la méthode des points sources distribués. Or celle-ci permet de formuler le problème direct sous forme matricielle ce qui constitue une base de nature à faciliter la résolution du problème inverse.C'est à ce problème qu'est consacrée la dernière partie du mémoire. Des méthodes d'inversion y sont proposées, visant à traiter progressivement la complexité du problème. Ainsi un algorithme d'inversion mono-fréquence est-il proposé, qui se montre efficace pour reconstruire des défauts sur de faibles profondeurs. Pour une meilleure reconstruction en profondeur, des algorithmes multifréquence faisant appel à des méthodes de régularisation sont ensuite conçus et appliqués à des signaux correspondant à des défauts de géométries diverses. / The ageing of the aeronautical fleet and the will to increase the aircrafts lifetime require the maintenance techniques to be made always more reliable and fast. In this context, the detection and characterization of the microscopic cracks likely to appear in some sensitive parts of the aircrafts is an important issue to be faced.This work deals with an original approach for the design of an eddy current imaging probe dedicated to the non destructive evaluation of cylindrical fastener parts. This approach consists in designing the probe in such a way that it both satisfies the constraint of an efficient sensing and that of enabling the inverse problem to be solved with a view to the defects reconstruction.Firstly, a probe design methodology is presented that optimizes the emission/reception topology, the choice of the probe elements and their arrangement according to sensitivity and resolution constraints. A prototype is built, implemented and characterized and validation of the design is obtained as the researched defects are displayed with experimental performances that agree with finite elements modelling simulations.Secondly, since the probe relies on a uniform eddy current flow interacting with small defects, a rather simple forward model is proposed based on virtual defect current sources (VDCS). The model implementation is carried out using the distributed point source method leading to a matrix formulation that can facilitate the resolution of the inverse problem.Finally, mono-frequency and multi-frequency methods are proposed for inverting the VDCS forward model and promising defect reconstruction results are obtained.

Numerical simulation allows investigation into the influence of separation distance and rotation on the performance of two vertical-axis hydrokinetic turbines. Compu- tational fluid dynamics is applied to calulate the lift and drag coefficients acting upon interacting NACA 0021 turbine blades for a Reynolds number of Red = 10, 000. To understand the effect of separation distance, large-eddy simulation of the flow around side-by-side and staggered cylinders, ReD = 3,000, and airfoils, Rec = 3,000, are also performed. Based upon the simulations, a drag reduction of 11.3% and 19.8% is determined for the downstream cylinder and airfoil, respectively. A reduction in Reynolds stresses is also observed for the staggered configuration compared to the side-by-side configuration. Due to computational resources of large-eddy simulation, the Reynolds averaged Navier-Stokes method is also applied to investigate the influence of separation distance and rotation on two vertical axis hydrokinetic turbines. The numerical simulations show that a drag reduction of 15.5% occurs when the non-dimensional spanwise and streamwise separation distances, based on turbine diameter, reach 1 and 2, respectively. / October 2016

La méthode des courants de Foucault (CF) est très largement utilisée en milieu industriel pour l’évaluation non destructive (END) de pièces ou de structures électriquement conductrices, parce qu’elle est sensible, robuste peu couteuse et non polluante. Toutefois, dans le cas général, l’évaluation quantitative des paramètres caractéristiques d’une structure est un problème difficile, d’une part parce que les données CF disponibles sont généralement incomplètes, d’autre part parce qu’il est nécessaire de faire appel à des modèle numériques élaborés pour rendre compte des interactions sonde / structure, et enfin, parce que le problème inverse consistant à estimer les paramètres de la structure à partir de la connaissance des interactions physiques et des données CF disponibles est « mal posé ».Dans cette thèse, nous nous intéressons au problème particulier de l’évaluation de jeux entre pièces dans un assemblage métallique. Pour contourner les difficultés liées à la mise en oeuvre de l’évaluation quantitative par CF, nous avons choisi de placer nos travaux dans un cadre multifréquence afin d’enrichir les données d’observation fournies par la sonde, et d’adopter une démarche consistant à élaborer des modèles comportementaux issus de l’analyse statistique des interactions sonde / structure. Partant de l’analyse de l’effet d’un jeu apparaissant entre les plaques d’un assemblage de plaques d’aluminium sur la variation de l’impédance normalisée de la sonde, menée expérimentalement et par simulations à l’aide de codes éléments finis, nous avons tout d’abord bâti un modèle comportemental approché d’interactions reposant sur la linéarité de la relation observée dans une bande de fréquences d’examen liée aux dimensions de la structure étudiée. Dans ce cadre, nous avons proposé, implanté et discuté des performances de plusieurs algorithmes de résolution du problème inverse, permettant d’estimer le jeu entre pièces et l’épaisseur de la pièce « cachée » de l’assemblage. Ensuite, nous avons étendu l’approche proposée à l’aide de réseaux de neurones artificiels, puis appliqué la méthode comportementale proposée à un second problème : celui de la caractérisation de fissures de fatigue réelles dans des pièces massives. Les résultats obtenus inclinent à considérer que la démarche proposée est généralisable à d’autres configurations en END par CF. / The eddy current (EC) method is widely used in the industry for the nondestructive evaluation (NDE) of electrically conductive structures, because it is sensitive, robust, cost effective and non polluting. However, in the general case, the quantitative evaluation of the parameters of a structure is a difficult problem, because the available EC data are generally incomplete, because it is necessary to have numerical models developed to account for the probe / structure interactions, and finally, because the inverse problem which consists in estimating the parameters of an inspected structure from the knowledge of those physical interactions and the available EC data is "ill-posed".In this work, we focus on the evaluation of air gaps between the parts of a metallic assembly. In order to cope with the difficulties associated with the implementation of an EC quantitative evaluation, we choose to use a multi-frequency approach so as to increase the data provided by the EC probe, and to build behavioral models from the statistical analysis of the probe / structure interactions. Experimental data and simulated data based on finite elements modeling are analyzed. A first approach behavioral model of the probe/structure interactions is deduced from these analyses, which is based on the linear relationship observed (in a particular excitation an frequency band depending on the dimensions of the studied structure) between the variations of the sensor normalized impedance and the air gap to evaluate. Moreover, we propose, implement and discuss the performances of several algorithms designed to solve the inverse problem dealing with the estimation of both the multilayered structure air gap and "hidden" plate thicknesses. Then, we extend the proposed multi-frequency behavioral approach thanks to artificial neural networks, and we apply the proposed behavioral method to a second problem : that of the characterization of real fatigue cracks in metallic massive parts. The obtained results let the generalization of the proposed approach to other EC NDE configurations to be envisaged.