In this thesis, an approach is presented for the generation of grids suitable for the simulation of high Reynolds number viscous flows in complex three-dimensional geometries. The automatic and reliable generation of such grids is today on the biggest bottlenecks in the industrial CFD simulation environment.<p><p>In the proposed approach, unstructured tetrahedral grids are employed for the regions far from the viscous boundaries of the domain, while semi-structured layers of high aspect ratio prismatic and hexahedral elements are used to provide the necessary grid resolution inside the boundary layers and normal to the viscous walls. The definition of the domain model is based on the STEP ISO standard and the topological information contained in the model is used for applying the hierarchical grid generation parameters defined by the user. An efficient, high-quality and robust algorithm is presented for the generation of the unstructured simplicial (triangular of tetrahedral) part of the grid. The algorithm is based on the Delaunay triangulation and the internal grid points are created following a centroid or frontal approach. For the surface grid generation, a hybrid approach is also proposed similar to the volume.<p>Semi-structured grids are generated on the surface grid (both on the edges and faces of the domain) to improve the grid resolution around convex and concave ridges and corners, by aligning the grid elements in the directions of high solution gradients along the surface. A method is also developed for automatically setting the grid generation parameters related to the surface grid generation based on the curvature of the surface in order to obtain an accurate and smooth surface grid. Finally, a semi-structured prismatic/hexahedral grid generation algorithm is presented for the generation of the part of grid close to the viscous walls of the domain. The algorithm is further extended with improvements meant to increase the grid quality around concave and convex ridges of the domain, where the semi-structured grids are known to be inadequate.<p><p>The combined methodology is demonstrated on a variety of complex examples mainly from the automotive and aeronautical industry. / Doctorat en sciences appliquées / info:eu-repo/semantics/nonPublished

L’émergence de bruits auparavant inaudibles dans les réservoirs à carburants automobiles requiert des constructeurs une meilleure compréhension des phénomènes physiques intervenants au sein de leurs produits. Dans cette thèse, différents travaux ont été conduits autour de l’étude du ballottement de fluide dans une cuve rigide rectangulaire partiellement remplie de fluide et soumise à une excitation extérieure. La première partie présente un état de l’art sur le sloshing suivant trois approches complémentaires - approche analytique, approche numérique et approche expérimentale - permettant d’orienter les travaux. Dans une deuxième partie, une étude préliminaire sur le sloshing dans une cuve rectangulaire soumise à une excitation harmonique forcée est réalisée. La confrontation des résultats numériques entre une approche linéaire - basée sur la théorie d’écoulement potentiel tenant compte de la viscosité du fluide [Schotté et Ohayon, 2013] - et une approche non linéaire commerciale – basée sur la résolution des équations de Navier-Stokes - permet de définir un paramètre de linéarité. Ce dernier permet de déterminer les cas de sloshing qui nécessitent une résolution non linéaire et ceux pour lesquels la théorie linéaire suffit pour prédire le phénomène. La troisième partie de ce document présente une étude expérimentale du ballottement de fluide dans une cuve rectangulaire rigide soumise à un freinage automobile. Deux niveaux de remplissage créant deux types d’impacts contre les parois (avec et sans enfermement de poche d’air) ont été analysés. Les essais menés ont permis de mesurer les forces engendrées par le mouvement du fluide, les pressions d’impact en paroi ainsi que le champ de vitesse par méthode Particle Image Velocimetry (PIV). Ce chapitre constitue une importante base de données expérimentales ayant permis d’étudier précisément le phénomène physique. L’étude est complétée par une confrontation des résultats expérimentaux avec des résultats Computational Fluid Dynamics (CFD). Enfin, pour conclure ce mémoire, une étude du sloshing dans un réservoir en tenant compte de la Fluid-Structure Interaction (FSI) est présentée. Le choix du couplage a été porté sur un schéma partitionné itératif faible avec, dans un premier temps, une approche potentielle instationnaire, puis avec une approche Volume Of Fluid (VOF) pour la physique fluide. Les limites d’un tel couplage dans le cas d’étude d’un réservoir partiellement rempli de fluide et attaché de manière flexible en fonction du rapport de masse fluide-réservoir ont été mises en évidence. La correction du schéma de couplage par l’effet de masse ajoutée présentée dans [Song et al., 2013] permet la résolution d’un système couplé quel que soit le rapport de masse en jeu et améliore de manière significative la convergence en réduisant également fortement le temps de calcul. / The present thesis focuses on an investigation of the sloshing phenomenon in a partially filled fuel tank submitted to a harmonic excitation motion. In the first part, the confrontation of numerical results between a linear approach - taking into account viscosity - and a nonlinear approach based on a commercial code leads to define a parameter of linearity. This parameter allows determining cases of sloshing who require non-linear resolution and those who need a linear theory to predict the phenomenon. An experimental study of fluid sloshing in a rectangular tank submitted to an automotive braking is conducted. Tests leaded allow measuring global forces engendered by the motion of the fluid, pressure of fluid impact and velocity field by PIV. This chapter provides an important data base and helps to investigate on the physical phenomenon. This study is completed by CFD results. To conclude, a numerical model for fluid-structure interactions is presented. Limits of this segregated partitioned coupling in case of sloshing in tank flexibly attached are highlighted, depending mostly on the mass ratio between fluid and tank structure. An added-mass term is integrated to the corrected staggered scheme ensuring systematically the convergence of the coupled solution and reducing significantly the iterations required.

The flow field and forces on an isolated oscillating NACA 0012 airfoil in a uniform flow is studied using viscous vortex particle method. The simulations are carried out at very low chord (c) based Reynolds number (Re=1000), motivated by the current interest in development of Micro Air Vehicles (MAV). The airfoil is forced to oscillate in both heave and pitch at different normalized oscillation frequencies (f), which is represented by the non-dimensional reduced frequency fc/U).( From the unsteady loading on the airfoil, the net energy transfer to the airfoil is calculated to determine the propensity for the airfoil to undergo self-induced oscillations or flutter at these very low Reynolds numbers. The simulations are carried out using a viscous vortex particle method that utilizes discrete vortex elements to represent the vorticity in the flow field. After validation of the code against test cases in the literature, simulations are first carried out for the stationary airfoil at different angles of attack, which shows the stall characteristics of the airfoil at this very low Reynolds numbers. For the airfoil oscillating in heave, the airfoil is forced to oscillate at different reduced frequencies at a large angle of attack in the stall regime. The unsteady loading on the blade is obtained at different reduced frequencies. This is used to calculate the net energy transfer to the airfoil from the flow, which is found to be negative in all cases studied. This implies that stall flutter or self-induced oscillations are not possible under the given heave conditions. The wake vorticity dynamics is presented for the different reduced frequencies, which show that the leading edge vortex dynamics is progressively more complex as the reduced frequency is increased from small values. For the airfoil oscillating in pitch, the airfoil is forced to oscillate about a large mean angle of attack corresponding to the stall regime. The unsteady moment on the blade is obtained at different reduced frequencies, and this is used to calculate the net energy transfer to the airfoil from the flow, which is found to be positive in all cases studied. This implies that stall flutter or self-induced oscillations are possible in the pitch mode, unlike in the heave case. The wake vorticity dynamics for this case is found to be relatively simple compared to that in heave. The results of the present simulations are broadly in agreement with earlier stall flutter studies at higher Reynolds numbers that show that stall flutter does not occur in the heave mode, but can occur in the pitch mode. The main difference in the present very low Reynolds number case appears to be the broader extent of the excitation region in the pitch mode compared to large Re cases studied earlier. region in the pitch mode compared to large Re cases studied earlier.

With steadily increasing computational power, computational fluid dynamics (CFD) can be applied to unsteady problems such as seakeeping simulations. Therefore, a good balance between accuracy and computational speed is required. This thesis investigates the application of CFD to seakeeping performance predictions and aims to propose a best-practice procedure for efficient seakeeping simulations. The widely used KVLCC2 research vessel serves as a test case for this thesis and FINEŠ/Marine software package is used for CFD computations. In order to validate the simulations, results are compared to recent experimental data from SSPA as well as predictions with potential ˛ow code SHIPFLOW® Motions. As for the calm water simulations, both inviscid and viscous ˛ow computations are performed in combination with three mesh refinement levels. Seakeeping simulations with regular head waves of different wavelengths are set-up correspondingly. Furthermore, different strategies for time discretization are investigated. With the given computational resources, it is not feasible to complete seakeeping simulations with a ˝ne mesh. However, already the coarse meshes give good agreement to experiments and SHIPFLOW® Motions' predictions. Viscous ˛ow simulations turn out to be more robust than Euler ˛ow computations and thus should be preferred. Regarding the time discretization, a fixed time discretization of 150 steps per wave period has shown the best balance between accuracy and speed. Based on these findings, a best-practice procedure for seakeeping performance predictions in FINEŠ/Marine is established. Taking the most efficient settings obtained from head wave simulations, the vessel is subjected to oblique waves with 160° encounter angle. Under similar wave conditions, CFD predictions of a similar thesis show close agreement in terms of added wave resistance. Compared to the previous head wave conditions of this study, added resistance in 160° oblique waves is found to be significantly higher. This underlines that oblique bow quartering waves represent a relevant case for determining the maximum required power of a ship. CFD and potential ˛ow show similar accuracy with respect to ship motions and added wave resistance, albeit potential ˛ow outperforms CFD in terms of computational speed. Hence, CFD should be applied in cases where viscous effects are known to have large influence on a vessel's seakeeping behavior. This can be the case if motion control and damping devices are to be evaluated, for instance. / Tack vare den stadigt ökande beräkningskraften kan beräkningsuiddynamik (CFD) idag användas på beräkningsintensiva problem som sjöegenskapssimulationer. Den här rapporten undersöker användning av CFD på sjöegenskapsprestanda och syftar till att foreslå ett best-practice förfaringssätt för effektiv sjöegenskapssimulationer. Forskningsskrovet KVLCC2 fungerar som ett testfall för denna rapport och FINE—/Marine-mjukvarupaketet används för CFD-beräkningar. Viktiga parametrar, såsom ödestyp, beräkningsnät och tidssteg varierars systematiskt. Resultaten jämförs med experiment gjorda vid SSPA. Baserat på resultaten förelås en best-practice. Den föreslagna best-practice användas vidare för berökningar av sjöegenskaper i sneda vågor. Jämförelse av resultaten med liknande studier visar god överensstämmelse. Genom att använda det föreslagna förfarandet för best-practice kan CFD-sjöegenskapssimulationer användas på fall där viskösa krafter måste beaktas, till exempel rörelseregleringsanordningar.

The present work deals with the extension of Acoustic Flux Vector Splitting (AFVS) scheme for the Compressible Viscous flow computations. Accurate viscous flow computations require much finer grids with adequate clustering of grid points in certain regions. Viscous flow computations are performed on unstructured triangulated grids. Solving Navier-Stokes equations involves the inviscid Euler part and the viscous part. The inviscid part of the fluxes are computed using the Acoustic Flux Vector Splitting scheme and the viscous part which is diffusive in nature does not require upwinding and is taken care using a central difference type of scheme. For these computations both the cell centered and the cell vertex finite volume methods are used. Higher order accuracy on unstructured meshes is achieved using the reconstruction procedure. Test cases are chosen in such a way that the performance of the scheme can be evaluated for different range of mach numbers. We demonstrate that higher order AFVS scheme in conjunction with a suitable grid adaptation strategy produce results that compare well with other well known schemes and the experimental data. An assessment of the relative performance of the AFVS scheme with the Roe scheme is also presented.

This thesis deals with the dynamics of a thin liquid film falling down a heated plate. The heating yields surface tension gradients that induce thermocapillary stresses on the free surface, thus affecting the stability and the evolution of the film. Accounting for the coherence of the flow due to viscosity, two main approaches that reduce the dimensionality of the original problem are usually considered depending on the flow rate (as measured by the Reynolds number): the `long wave' asymptotic expansion for small Reynolds numbers and the `integral boundary layer' approximation for moderate Reynolds numbers. The former suffers from singularities and the latter from incorrectness of the instability threshold for the occurrence of hydrodynamic waves. Thus, the aim of this thesis is twofold: in a first part, we define quantitatively the validity of the `long wave' evolution equation (Benney equation) for the film thickness h including the thermocapillary effect; and in a second part, we improve the `integral boundary layer' approach by combining a gradient expansion to a weighted residual method. <p>In the first part, we further investigate the Benney equation in its validity domain in the case of periodically inhomogeneous heating in the streamwise direction. It induces steady-state deformations of the free surface with increased transfer rate in regions where the film is thinner, and also in average. The inhomogeneities of the heating also modify the nature of travelling wave solutions at moderate temperature gradients and allows for suppressing wave motion at larger ones.<p>Moreover, large temperature gradients (for instance positive ones) in the streamwise direction produce large local film thickening that may in turn become unstable with respect to transverse disturbances such that the flow may organize in rivulet-like structures. The mechanism of such instability is elucidated via an energy analysis. The main features of the rivulet pattern are described experimentally and recovered by direct numerical simulations.<p>In the second part, various models are obtained, which are valid for larger Reynolds numbers than the Benney equation and account for second-order viscous and inertial effects. We then elaborate a strategy to select the optimal model in terms of linear stability properties and existence of nonlinear solutions (solitary waves), for the widest possible range of parameters. This model -- called reduced model -- is a system of three coupled evolution equations for the local film thickness h, the local flow rate q and the surface temperature Ts. Solutions of this model indicate that the interaction of the hydrodynamic and thermocapillary modes is non-trivial, especially in the region of large-amplitude solitary waves.<p>Finally, the three-dimensional evolution of the solutions of the reduced model in the presence of periodic forcing and noise compares favourably with available experimental data in isothermal conditions and with direct numerical simulations in non-isothermal conditions.<p><p>------------------------------------------------<p><p>Cette thèse analyse la dynamique d'un film mince s'écoulant le long d'une paroi chauffée. Le chauffage crée des gradients de tension superficielle qui induisent des tensions thermocapillaires à la surface libre, altérant ainsi la stabilité et l'évolution du film. Grâce à la cohérence de l'écoulement assurée par la viscosité, deux approches permettant de réduire la dimensionnalité du problème original sont habituellement considérées suivant le débit (mesuré par le nombre de Reynolds): l'approximation asymptotique dite `longues ondes' pour les faibles nombres de Reynolds et l'approximation `intégrale couche limite' pour les nombres de Reynolds modérés. Cependant, la première approximation souffre de singularités et la dernière de prédictions imprécises du seuil de stabilité des ondes hydrodynamiques à la surface du film. Le but de cette thèse est donc double: dans une première partie, il s'agit de déterminer, de manière quantitative, la validité de l'équation d'évolution `longues ondes' (ou équation de Benney) pour l'épaisseur du film h, en y incluant l'effet thermocapillaire; et dans une seconde partie, il s'agit d'améliorer l'approche `intégrale couche limite' en combinant un développement en gradients avec une méthode aux résidus pondérés.<p>Dans la première partie, nous étudions l'équation de Benney, dans son domaine de validité, dans le cas d'un chauffage inhomogène et périodique dans la direction de l'écoulement. Cela induit des déformations permanentes de la surface libre avec un accroissement du transfert de chaleur dans les régions où le film est plus mince, mais aussi en moyenne. Un chauffage inhomogène modifie également la nature des solutions d'ondes progressives pour des gradients de températures modérés et conduit même à leur suppression pour des gradients de températures plus importants. De plus, ceux-ci, lorsqu'ils sont par exemple positifs le long de l'écoulement, produisent des épaississements localisés du film qui peuvent à leur tour devenir instables par rapport à des perturbations suivant la direction transverse à l'écoulement. Ce dernier s'organise alors sous forme d'une structure en rivulets. Le mécanisme de cette instabilité est élucidé via une analyse énergétique des perturbations. Les principales caractéristiques des structures en rivulets sont décrites expérimentalement et retrouvées par l'intermédiaire de simulations numériques. <p>Dans la seconde partie, nous dérivons une famille de modèles valables pour des nombres de Reynolds plus grands que l'équation de Benney, qui prennent en compte les effets visqueux et inertiels du second ordre. Nous élaborons ensuite une stratégie pour sélectionner le modèle optimal en fonction de ses propriétés de stabilité linéaire et de l'existence de solutions non-linéaires (ondes solitaires), et ce pour la gamme de paramètres la plus large possible. Ce modèle -- appelé modèle réduit -- est un système de trois équations d'évolution couplées pour l'épaisseur locale de film h, le débit local q et la température de surface Ts. Les solutions de ce modèle indiquent que l'interaction des modes hydrodynamiques et thermocapillaires n'est pas triviale, spécialement dans le domaine des ondes solitaires de grande amplitude. Finalement, l'évolution tri-dimensionnelle des solutions du modèle réduit en présence d'un forçage périodique ou d'un bruit se compare favorablement aux données expérimentales disponibles en conditions isothermes, ainsi qu'aux simulations numériques directes en conditions non-isothermes<p> / Doctorat en sciences appliquées / info:eu-repo/semantics/nonPublished

La capacité à percevoir des courants dans un fluide s'est développée chez de nombreuses espèces animales, dans des contextes écologiques très variés qui couvrent aussi bien les interactions proies-prédateurs, la sélection sexuelle ou l'orientation dans un environnement. Parmi ces espèces animales, les grillons détectent les courants d'air générés notamment lors de l'attaque de leurs prédateurs à l'aide de deux organes appelés "cerques", situés à l'arrière de leur abdomen et recouverts de poils mécano-sensoriels. Ces senseurs sont considérés comme les détecteurs les plus sensibles du monde animal. Il leur suffit de capter l'énergie d'un dixième d'un photon pour déclencher un potentiel d'action au niveau du neurone sensoriel. Ce manuscrit présente à la fois le développement des outils de mesures sans contact adaptés à ces questions d'écologie sensorielle ainsi que les méthodes numériques simulant les processus physiques à l'oeuvre. L'étude du fonctionnement des senseurs a nécessité l'adaptation des méthodes de mesures non intrusives de très grande précision tel que la Vélocimétrie par Imagerie de Particules (PIV). La couche limite oscillante dans laquelle évoluent les poils a été visualisée et a servi à déterminer la réponse de poils modélisés par des systèmes oscillatoires du second ordre. Le couplage visqueux entre poils a été lui aussi caractérisé en adaptant la PIV à des mesures à très petites échelles sur des poils biomimétiques micro-electro-mécanique (MEMS). Les mesures des perturbations générées lors des attaques d'araignées, principales prédatrices des grillons, nous ont aidé à valider des modélisations numériques, réalisées à l'aide des techniques de dynamique des fluides computationnelles (CFD) par résolution des équations de Navier Stokes via la méthode des éléments finis (FEM). La mise au point et l'utilisation de techniques de métrologie optique en dynamique des fluides semi-visqueux et l'analyse des données nous permettent de revisiter la sensibilité extrême du système sensoriel du grillon et de placer ces mesures dans un contexte plus large, d'écologie sensorielle. En particulier, nous montrons que ces soies sont placées en groupe compact et exercent entre elles un fort couplage aérodynamique visqueux, qui réduit fortement leur sensibilité "de groupe". Ce fort couplage interroge l'intérêt d'avoir des récepteurs aussi performants individuellement, s'ils perdent leur sensibilité lorsqu'ils fonctionnent en réseau. Finalement, les réactions des poils à des mouvements de fluides générés par un piston mimant les attaques réelles d'araignées ont pu être déterminées à l'aide d'une caméra rapide, puis simulées et validées après avoir développé un modèle mécanique du poil répondant à des stimuli transitoires. / Flow sensing is used by a vast number of animals in various ecological contexts, from preypredator interactions to mate selection, and orientation to flow itself. Among these animals, crickets use hundreds of filiform hairs on two cerci as an early warning system to detect remote potential predators. Over the years, the cricket hairs have been described as the most sensitive sensor in the animal kingdom. The energy necessary for the emission of an action potential by its sensory neuron was estimated to be a tenth of the energy of a photon. This PhD thesis aims to describe recent technological advances in the measurement and model of flows around biological and artificial flow sensors in the context of organismal sensory ecology. The study and understanding of the performance of sensory systems requires a high spatial precision of non-intrusive measurement methods. Thus, non-contacting measurement methods such as and Particle Image Velocimetry (PIV), originally developed by aerodynamics and fluid mechanics engineers, have been used to measure flows of biological relevance. The viscous oscillatory boundary layer surrounding filiform hairs has been visualized and used as input to model the mechanical response of these hairs, described as second order mechanical systems. The viscous hydrodynamic coupling occurring within hair canopy was also characterized using PIV measurements on biomimetic micro-electro-mechanical systems (MEMS) hairs, mimicking biological ones. Using PIV, we have also measured the air flow upstream of hunting spiders. We prove that this flow is highly conspicuous aerodynamically, due to substantial air displacement detectable up to several centimeters in front of the running predator. This disturbance of upstream air flows were also assessed by computational fluid dynamics (CFD) with the finite elements method (FEM). The development of non-intrusive measurement and CFD methods and their application to the analysis of the biological flow involved in cricket sensory ecology allowed us to revisit the extreme sensitivity of cricket filiform hairs. We predicted strong hydrodynamic coupling within natural hair canopies and we addressed why hairs are packed together at such high densities, particularly given the exquisite sensitivity of a single hair. We also proposed a new model of hair deflection during the arrival of a predator, by taking into account both the initial and long-term aspects of the flow pattern produced by a lunging predator. We conclude that the length heterogeneity of the hair canopy mirrors the flow complexity of an entire attack, from launch to grasp.

A high-order central essentially non-oscillatory (CENO) finite-volume scheme in combination with a block-based adaptive mesh refinement (AMR) algorithm is proposed for solution of hyperbolic and elliptic systems of conservation laws on body- fitted multi-block mesh. The spatial discretization of the hyperbolic (inviscid) terms is based on a hybrid solution reconstruction procedure that combines an unlimited high-order k-exact least-squares reconstruction technique following from a fixed central stencil with a monotonicity preserving limited piecewise linear reconstruction algorithm. The limited reconstruction is applied to computational cells with under-resolved solution content and the unlimited k-exact reconstruction procedure is used for cells in which the solution is fully resolved. Switching in the hybrid procedure is determined by a solution smoothness indicator. The hybrid approach avoids the complexity associated with other ENO schemes that require reconstruction on multiple stencils and therefore, would seem very well suited for extension to unstructured meshes. The high-order elliptic (viscous) fluxes are computed based on a k-order accurate average gradient derived from a (k+1)-order accurate reconstruction. A novel h-refinement criterion based on the solution smoothness indicator is used to direct the steady and unsteady refinement of the AMR mesh. The predictive capabilities of the proposed high-order AMR scheme are demonstrated for the Euler and Navier-Stokes equations governing two-dimensional compressible gaseous flows as well as for advection-diffusion problems characterized by the full range of Peclet numbers, Pe. The ability of the scheme to accurately represent solutions with smooth extrema and yet robustly handle under-resolved and/or non-smooth solution content (i.e., shocks and other discontinuities) is shown for a range of problems. Moreover, the ability to perform mesh refinement in regions of smooth but under-resolved and/or non-smooth solution content to achieve the desired resolution is also demonstrated.

A high-order central essentially non-oscillatory (CENO) finite-volume scheme in combination with a block-based adaptive mesh refinement (AMR) algorithm is proposed for solution of hyperbolic and elliptic systems of conservation laws on body- fitted multi-block mesh. The spatial discretization of the hyperbolic (inviscid) terms is based on a hybrid solution reconstruction procedure that combines an unlimited high-order k-exact least-squares reconstruction technique following from a fixed central stencil with a monotonicity preserving limited piecewise linear reconstruction algorithm. The limited reconstruction is applied to computational cells with under-resolved solution content and the unlimited k-exact reconstruction procedure is used for cells in which the solution is fully resolved. Switching in the hybrid procedure is determined by a solution smoothness indicator. The hybrid approach avoids the complexity associated with other ENO schemes that require reconstruction on multiple stencils and therefore, would seem very well suited for extension to unstructured meshes. The high-order elliptic (viscous) fluxes are computed based on a k-order accurate average gradient derived from a (k+1)-order accurate reconstruction. A novel h-refinement criterion based on the solution smoothness indicator is used to direct the steady and unsteady refinement of the AMR mesh. The predictive capabilities of the proposed high-order AMR scheme are demonstrated for the Euler and Navier-Stokes equations governing two-dimensional compressible gaseous flows as well as for advection-diffusion problems characterized by the full range of Peclet numbers, Pe. The ability of the scheme to accurately represent solutions with smooth extrema and yet robustly handle under-resolved and/or non-smooth solution content (i.e., shocks and other discontinuities) is shown for a range of problems. Moreover, the ability to perform mesh refinement in regions of smooth but under-resolved and/or non-smooth solution content to achieve the desired resolution is also demonstrated.