Nous nous intéressons dans cette thèse aux équations de Navier-Stokes pour un fluide visqueux incompressible. Dans la première partie, nous étudions le cas d’un fluide homogène. Rappelons que la grande question de la régularité globale en dimension 3 est plus ouverte que jamais : on ne sait pas si la solution de l’équation correspondant à un état initial suffisamment régulier mais arbitrairement loin du repos, va perdurer indéfiniment dans cet état (régularité globale) ou exploser en temps fini(singularité). Une façon d’aborder le problème est de supposer cette éventuelle rupture de régularité et d’envisager les différents scenarii possibles. Après un rapide survol de la structure propre aux équations de Navier-Stokes et des résultats connus à ce jour (chapitre 1), nous nous intéressons(chapitre 2) à l’existence locale (en temps) de solutions dans des espaces de Sobolev qui ne sont pas invariants d’échelle. Partant d’une donnée initiale qui produit une singularité, on prouve l’existence d’une constante optimale qui minore le temps de vie de la solution. Cette constante, donnée parla méthode rudimentaire du point fixe, fournit ainsi un bon ordre de grandeur sur le temps de vie maximal de la solution. Au chapitre 3, nous poursuivons les investigations sur le comportement de telles solutions explosives à la lumière de la méthode des éléments critiques.Dans le seconde partie de la thèse, nous sommes intéressés à un modèle plus réaliste du point de vue de la physique, celui d’un fluide incompressible à densité variable. Ceci est modélisé par les équations de Navier-Stokes incompressible et inhomogènes. Nous avons étudié le caractère globalement bien posé de ces équations dans la situation d’un fluide évoluant dans un tore de dimension 3, avec des données initiales appartenant à des espaces critiques et sans hypothèse de petitesse sur la densité. / This thesis is concerned with incompressible Navier-Stokes equations for a viscous fluid. In the first part, we study the case of an homogeneous fluid. Let us recall that the big question of the global regularity in dimension 3 is still open : we do not know if the solution associated with a data smooth enough and far from the immobile stage will last over time (global regularity) or on the contrary will stop living in finite time and blow up (singularity). The goal of this thesis is to study this regularity break. One way to deal witht his question is to assume that such a phenomen on occurs and to study differents scenarii. The chapter 1 is devoted to a recollection of well-known results. In chapter 2, we are interesting in the local (in time) existence of a solution in some Sobolev spaces which are not invariant under the natural sclaing of Navier-Stokes. Starting with a data generating a singularity, we can prove there exists an optimal lower boundary of the lifes pan of such a solution. In this way, the lower boundary provided by the elementary procedure of fixed-point, gives the correctorder of magnitude. Then, we keep on investigations about the behaviour of regular solution near the blow up, thanks to the method of critical elements (chapter 3).In the second part, we are concerned with a more relevant model, from a physics point of view : the inhomogeneous Navier-Stokes system. We deal with the global well poseness of such a model for a inhomogeneous fluid, evolving on a tor us in dimension 3, with critical data and without smallnes sassumption on the density.

Investigation of material behaviour in a squeeze flow geometry provides an impor- tant technique in rheology and it is relevant also from the technological point of view (some types of dampers, compression moulding). To our best knowledge, the sque- eze flow has not been solved for fluids-like materials with pressure-dependent material moduli. In the main scope of the present thesis, an incompressible fluid whose visco- sity strongly depends on the pressure is studied in both the perfect-slip and the no-slip squeeze flow. It is shown that such a material model can provide interesting departures compared to the classical model for viscous (Navier-Stokes) fluid even on the level of analytical solutions, which are obtained using some physically relevant simplificati- ons. Numerical simulation of a free boundary problem for the no-slip squeeze flow is then developed in the thesis using body-fitted curvilinear coordinates and spectral collocation method. An interesting behaviour is expected especially in the corners of the computational domain where the stress singularities are normally located. Unfor- tunately, numerical results reveal some fundamental drawbacks related to the physical model and its possible improvement is discussed at the end of the thesis.

Cette thèse est consacrée à l'étude des opérateurs strictement hyperboliques à coefficients peu réguliers, aussi bien qu'à l'étude du système d'Euler incompressible à densité variable. Dans la première partie, on montre des estimations a priori pour des opérateurs strictement hyperboliques dont les coefficients d'ordre le plus grand satisfont une condition de continuité log-Zygmund par rapport au temps et une condition de continuité log-Lipschitz par rapport à la variable d'espace. Ces estimations comportent une perte de dérivées qui croît en temps. Toutefois, elles sont suffisantes pour avoir encore le caractère bien posé du problème de Cauchy associé dans l'espace H^inf (pour des coefficients du deuxième ordre ayant assez de régularité).Dans un premier temps, on considère un opérateur complet en dimension d'espace égale à 1, dont les coefficients du premier ordre étaient supposés hölderiens et celui d'ordre 0 seulement borné. Après, on traite le cas général en dimension d'espace quelconque, en se restreignant à un opérateur de deuxième ordre homogène: le passage à la dimension plus grande exige une approche vraiment différente. Dans la deuxième partie de la thèse, on considère le système d'Euler incompressible à densité variable. On montre son caractère bien posé dans des espaces de Besov limites, qui s'injectent dans la classe des fonctions globalement lipschitziennes, et on établit aussi des bornes inférieures pour le temps de vie de la solution ne dépendant que des données initiales. Cela fait, on prouve la persistance des structures géométriques, comme la régularité stratifiée et conormale, pour les solutions de ce système. À la différence du cas classique de densité constante, même en dimension 2 le tourbillon n'est pas transporté par le champ de vitesses. Donc, a priori on peut s'attendre à obtenir seulement des résultats locaux en temps. Pour la même raison, il faut aussi laisser tomber la structure des poches de tourbillon. La théorie de Littlewood-Paley et le calcul paradifférentiel nous permettent d'aborder ces deux différents problèmes. En plus, on a besoin aussi d'une nouvelle version du calcul paradifférentiel, qui dépend d'un paramètre plus grand que ou égal à 1, pour traiter les opérateurs à coefficients peu réguliers. Le cadre fonctionnel adopté est celui des espaces de Besov, qui comprend en particulier les ensembles de Sobolev et de Hölder. Des classes intermédiaires de fonctions, de type logarithmique, entrent, elles aussi, en jeu / The present thesis is devoted both to the study of strictly hyperbolic operators with low regularity coefficients and of the density-dependent incompressible Euler system. On the one hand, we show a priori estimates for a second order strictly hyperbolic operator whose highest order coefficients satisfy a log-Zygmund continuity condition in time and a log-Lipschitz continuity condition with respect to space. Such an estimate involves a time increasing loss of derivatives. Nevertheless, this is enough to recover well-posedness for the associated Cauchy problem in the space $H^infty$ (for suitably smooth second order coefficients).In a first time, we consider acomplete operator in space dimension $1$, whose first order coefficients were assumed Hölder continuous and that of order $0$only bounded. Then, we deal with the general case of any space dimension, focusing on a homogeneous second order operator: the step to higher dimension requires a really different approach. On the other hand, we consider the density-dependent incompressible Euler system. We show its well-posedness in endpoint Besov spaces embedded in the class of globally Lipschitz functions, producing also lower bounds for the lifespan of the solution in terms of initial data only. This having been done, we prove persistence of geometric structures, such as striated and conormal regularity, for solutions to this system. In contrast with the classical case of constant density, even in dimension $2$ the vorticity is not transported by the velocity field. Hence, a priori one can expect to get only local in time results. For the same reason, we also have to dismiss the vortex patch structure. Littlewood-Paley theory and paradifferential calculus allow us to handle these two different problems .A new version of paradifferential calculus, depending on a parameter $ggeq1$, is also needed in dealing with hyperbolic operators with nonregular coefficients. The general framework is that of Besov spaces, which includes in particular Sobolev and Hölder sets. Intermediate classes of functions, of logaritmic type, come into play as well

Past studies (primarily on steady state problems) that have compared the penalty and the velocity-pressure finite element formulations on a variety of problems have concluded that both methods yield solutions of comparable accuracy, and that the choice of one method over the other is dictated by which of the two is more efficient. In this work, we show that the penalty finite element method yields inaccurate solutions at large times on a class of transient problems, while the velocity-pressure formulation yields solutions that are in good agreement with the analytical solution. Numerical studies are conducted on various problems to compare these two formulations on the basis of rates of convergence, total number of equations to be solved and accuracy of results. We found that both formulations give almost the same rates of convergence in all problems, however the penalty formulation involves lesser number of equations than the velocity-pressure formulation due to implicit treatment of pressure field, and hence is more efficient. In some of the problems we have also compared a finite volume method with the penalty and velocity-pressure formulations on the basis of accuracy and computational cost.

La méthode IIM (Immersed Interface Method) permet d'étendre certaines méthodes numériques à des problèmes présentant des discontinuités. Elle est utilisée ici pour étudier un fluide incompressible régi par les équations de Navier-Stokes, dans lequel est immergée une membrane exerçant une force singulière. Nous utilisons une méthode de projection dans une grille de différences finies de type MAC. Une dérivation très complète des conditions de saut dans le cas où la viscosité est continue est présentée en annexe. Deux exemples numériques sont présentés : l'un sans membrane, et l'un où la membrane est immobile. Le cas général d'une membrane mobile est aussi étudié en profondeur. / The Immersed Interface Method allows us to extend the scope of some numerical methods to discontinuous problems. Here we use it in the case of an incompressible fluid governed by the Navier-Stokes equations, in which a membrane is immersed, inducing a singular force. We use a projection method and staggered (MAC-type) finite difference approximations. A very complete derivation for the jump conditions is presented in the Appendix, for the case where the viscosity is continuous. Two numerical examples are shown : one without a membrane, and the other where the membrane is motionless. The general case of a moving membrane is also thoroughly studied.

Made available in DSpace on 2014-06-11T19:24:47Z (GMT). No. of bitstreams: 0 Previous issue date: 2005-02-21Bitstream added on 2014-06-13T19:31:51Z : No. of bitstreams: 1 pereira_vd_me_ilha.pdf: 1970071 bytes, checksum: 74646e1883439e38fa62ff7f34d06488 (MD5) / Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES) / Neste trabalho foram feitas simulações de escoamentos incompressiveis por um método de elementos finitos de mínimos quadrados (LSFEM – Least Squares Finite Element Method), usando as formulações velocidade-pressão-vorticidade e velocidade-pressão-tensão, denominadas na literatura de formulações u − p −ω e u = p −τ respectivamente. Estas formulações são preferidas por resultarem em sistemas de equações diferenciais de primeira ordem, o que é mais conveniente para implementação pelo LSFEM. O objetivo principal deste trabalho é a simulação computacional de escoamentos laminares, transicionais e turbulentos através da aplicação da metodologia de simulação de grandes escalas (LES – Large Eddy Simulation) com o modelo de viscosidade turbulenta de Smagorinky para modelar as tensões submalha. Alguns problemas padrões foram resolvidos para validar um código computacional desenvolvido e os resultados são apresentados e comparados com resultados disponíveis na literatura. / In this work simulations of incompressible fluid flows have been done by a Least Squares Finite Element Method (LSFEM) using the velocity-pressure-vorticity and velocity-pressurestress formulations, named, in the literature, u − p −ω and u = p −τ formulations respectively. These formulations are preferred because the resulting equations are partial differential equations of first order, which is more convenient for implementation by LSFEM. The main purpose of this work are the numerical computations of laminar, transitional and turbulent fluid flows through the application of large eddy simulation (LES) methodology using the LSFEM. The Navier- Stokes equations in u − p −ω and u = p −τ formulations are filtered and the eddy viscosity model of Smagorinsky is used for modeling the sub-grid-scale stresses. Some benchmark problems are solved for validate a developed numerical code and the preliminary results are presented and compared with available results from the literature.

Orientador: João Batista Campos Silva / Banca: João Batista Aparecido / Banca: Luiz Felipe Mendes de Moura / Resumo: Neste trabalho foram feitas simulações de escoamentos incompressiveis por um método de elementos finitos de mínimos quadrados (LSFEM - Least Squares Finite Element Method), usando as formulações velocidade-pressão-vorticidade e velocidade-pressão-tensão, denominadas na literatura de formulações u − p −ω e u = p −τ respectivamente. Estas formulações são preferidas por resultarem em sistemas de equações diferenciais de primeira ordem, o que é mais conveniente para implementação pelo LSFEM. O objetivo principal deste trabalho é a simulação computacional de escoamentos laminares, transicionais e turbulentos através da aplicação da metodologia de simulação de grandes escalas (LES - Large Eddy Simulation) com o modelo de viscosidade turbulenta de Smagorinky para modelar as tensões submalha. Alguns problemas padrões foram resolvidos para validar um código computacional desenvolvido e os resultados são apresentados e comparados com resultados disponíveis na literatura. / Abstract: In this work simulations of incompressible fluid flows have been done by a Least Squares Finite Element Method (LSFEM) using the velocity-pressure-vorticity and velocity-pressurestress formulations, named, in the literature, u − p −ω and u = p −τ formulations respectively. These formulations are preferred because the resulting equations are partial differential equations of first order, which is more convenient for implementation by LSFEM. The main purpose of this work are the numerical computations of laminar, transitional and turbulent fluid flows through the application of large eddy simulation (LES) methodology using the LSFEM. The Navier- Stokes equations in u − p −ω and u = p −τ formulations are filtered and the eddy viscosity model of Smagorinsky is used for modeling the sub-grid-scale stresses. Some benchmark problems are solved for validate a developed numerical code and the preliminary results are presented and compared with available results from the literature. / Mestre

Le schéma Marker-And-Cell (MAC) est un schéma de discrétisation des équations aux dérivées partielles sur maillages cartésiens, très connu en mécanique des fluides. Nous nous intéressons ici à son analyse mathématique dans le cadre des écoulements incompressibles sur des maillages cartésiens non-uniformes en dimension 2 ou 3. Dans un premier temps nous discrétisons les équations de Navier-Stokes pour un écoulement incompressible stationnaire; nous établissons des estimations a priori sur les suites de vitesses et pressions approchées qui permettent d’une part d'établir l’existence d’une solution au schéma, et d’obtenir la compacité de ces suites lorsque le pas d’espace tend vers 0. Nous montrons alors la convergence de ces suites (à une sous-suite près) vers une solution faible du problème continu, ce qui nécessite une analyse fine du terme de convection non linéaire. Nous nous intéressons ensuite aux équations de Navier-Stokes en régime instationnaire avec une discrétisation en temps implicite. Nous démontrons que le schéma préserve les propriétés de stabilité du problème continu et obtenons ainsi l’existence d’une solution au schéma. Puis, grâce à des techniques de compacité et en passant à la limite dans le schéma, nous démontrons qu’une suite de vitesses approchées converge. Si l’on se restreint au problème de Stokes, et en supposant de plus que la condition initiale de la vitesse est dans H 1 , nous obtenons une estimation sur la pression qui permet de montrer la convergence forte des pressions approchées. Enfin nous étendons l’analyse aux écoulements incompressibles à masse volumique variable. On montre la convergence du schéma. / The Marker-And-Cell (MAC) scheme is a discretization scheme for partial derivative equations on Cartesian meshes, which is very well known in fluid mechanics. Here we are concerned with its mathematical analysis in the case of incompressible flows on two or three dimensional non-uniform Cartesian grids. We first discretize the steady-state incompressible Navier-Stokes equations. We show somea priori estimates that allow to show the existence of a solution to the scheme and some compactness and consistency results. By a passage to the limit on the scheme, we show that the approximate solutions obtained with the MAC scheme converge (up to a subsequence) to a weak solution of the Navier-Stokes equations, thanks to a careful analysis of the nonlinear convection term. Then, we analyze the convergence of the unsteady-case Navier-Stokes equations. The algorithm is implicit in time. We first show that the scheme preserves the stability properties of the continuous problem, which yields, the existence of a solution. Then, invoking compactness arguments and passing to the limit in the scheme, we prove that any sequence of solutions (obtained with a sequence of discretizations the space and time step of which tend to zero) converges up to the extraction of a subsequence to a weak solution of the continuous problem. If we restrict ourselves to the Stokes equations and assume that the initial velocity belongs to H 1, then we obtain estimates on the pressure and prove the convergence of the sequences of approximate pressures. Finally, we extend the analysis of the scheme to incompressible variable density flows. we show the convergence of the scheme.

La méthode IIM (Immersed Interface Method) permet d'étendre certaines méthodes numériques à des problèmes présentant des discontinuités. Elle est utilisée ici pour étudier un fluide incompressible régi par les équations de Navier-Stokes, dans lequel est immergée une membrane exerçant une force singulière. Nous utilisons une méthode de projection dans une grille de différences finies de type MAC. Une dérivation très complète des conditions de saut dans le cas où la viscosité est continue est présentée en annexe. Deux exemples numériques sont présentés : l'un sans membrane, et l'un où la membrane est immobile. Le cas général d'une membrane mobile est aussi étudié en profondeur. / The Immersed Interface Method allows us to extend the scope of some numerical methods to discontinuous problems. Here we use it in the case of an incompressible fluid governed by the Navier-Stokes equations, in which a membrane is immersed, inducing a singular force. We use a projection method and staggered (MAC-type) finite difference approximations. A very complete derivation for the jump conditions is presented in the Appendix, for the case where the viscosity is continuous. Two numerical examples are shown : one without a membrane, and the other where the membrane is motionless. The general case of a moving membrane is also thoroughly studied.

Experiments are done using a smooth solid rigid homogeneous acrylic sphere rolling on an inclined plane which is submerged in water. The motivation for these experiments comes from a need to understand a class of solid-fluid interaction problems that include sediment transport, movement of gravel on ocean floor and river bed due to water currents. Experiments are performed in a glass water tank 15 cm wide by 14 cm deep by 61 cm long which can be tilted to desired angle. The sphere is released from rest on the inclined false bottom of the tank in quiescent water. Our experimental study has twofold aim: (1)to study the boundary layer separation, the three-dimensional eddying motion in the wake and the near-wake structure and(2) to establish hydrodynamic force coefficients by analyzing kinematical data of the sphere motion from start to till it attains terminal velocity. Experiments are carried out at moderate Reynolds number Rearound1500. Previous studies on the first problem exist in the literature for Reup to 350. Previous studies on the second problem do not clearly define the added-mass coefficient and the influence of the water tank side-walls on the drag coefficient. In the first study, the characterization of the wake is done using flow visualization methods (fluoresce in dye visualization and particle streak visualization) and Particle Image Velocimetry (PIV). Laser light sheet obtained from an argon ion continuous laser beam is taken in different orientations to illuminate the fluoresce in dye or 14 m silver-coated hollow glass spheres. These experiments show that the wake behind the rolling sphere up to 1.6 diameters (or 1.6D) downstream is confined within height 1.2Dand width1.2D. At about 1.8Ddownstream, the wake sways alternately on either side of the equatorial plane, moving in lateral-vertical direction and moving out of the confining region; this gives zigzag appearance to the wake. Also in these experiments, we observe that the flow separations from the surface of the rolling sphere show three separation zones. The eddies shed from the primary separation surface on the upper hemisphere are symmetrical about the equatorial plane with Strouhal number St=1.0. The primary separation is affected by the symmetrical secondary separations on the rear surface in the piggyback region — it is the region near the upper rear surface of the sphere behind the transverse equatorial plane and below the primary separation surface. The lower eddies below the primary separation zone are shed alternately on either side of the equatorial plane with shedding frequency St=0.5. Our experiments show that there is a viscous blockage of width 0.4Dat the crevice near the point of contact. On either side of the viscous blockage at the crevice, we see weak symmetric eddies. Based on our experimental observations, we proceed to build a simple physical model of the separated flow on the surface of the rolling sphere. In the second study, the motion of the sphere is photographed and paired data of the displacement and time is obtained for the sphere motion from the start of motion till terminal velocity is reached at about 4.5 sphere diameters from the point of release of the sphere. Equation of motion of the sphere is solved numerically treating added-mass coefficient Ca and drag coefficient Cd as parameters. Experimental data is fitted on these solutions and the best fit gives the values of the force coefficients. Theoretical value of Ca equal to 0.621 is confirmed experimentally. Value of Cd is found to be 1.23 at Re=990 and it is 1.06 at Re= 1900. Side-wall effects become important for ratio of diameter of sphere to width of tank greaterthan0.20.