We present an immersed interface algorithm for the incompressible Navier Stokes equations. The interface is represented by cubic splines which are interpolated through a set of Lagrangian control points. The position of the control points is implicitly updated using the fluid velocity. The forces that the interface exerts on the fluid are computed from the constitutive relation of the interface and are applied to the fluid through jumps in the pressure and jumps in the derivatives of pressure and velocity. A projection method is used to time advance the Navier-Stokes equations on a uniform cartesian mesh. The Poisson-like equations required for the implicit solution of the diffusive and pressure terms are solved using a fast Fourier transform algorithm. The position of the interface is updated implicitly using a quasi-Newton method (BFGS) within each timestep. Several examples are presented to illustrate the flexibility of the presented approach. / Singapore-MIT Alliance (SMA)

The incompressible Navier-Stokes equations are among the most important partial differential systems arising from classical physics. They are utilized to model a wide range of fluids, from water moving around a naval vessel to blood flowing through the arteries of the cardiovascular system. Furthermore, the secrets of turbulence are widely believed to be locked within the Navier-Stokes equations. Despite the enormous applicability of the Navier-Stokes equations, the underlying behavior of solutions to the partial differential system remains little understood. Indeed, one of the Clay Mathematics Institute's famed Millenium Prize Problems involves the establishment of existence and smoothness results for Navier-Stokes solutions, and turbulence is considered, in the words of famous physicist Richard Feynman, to be "the last great unsolved problem of classical physics." Numerical simulation has proven to be a very useful tool in the analysis of the Navier-Stokes equations. Simulation of incompressible flows now plays a major role in the industrial design of automobiles and naval ships, and simulation has even been utilized to study the Navier-Stokes existence and smoothness problem. In spite of these successes, state-of-the-art incompressible flow solvers are not without their drawbacks. For example, standard turbulence models which rely on the existence of an energy spectrum often fail in non-trivial settings such as rotating flows. More concerning is the fact that most numerical methods do not respect the fundamental geometric properties of the Navier-Stokes equations. These methods only satisfy the incompressibility constraint in an approximate sense. While this may seem practically harmless, conservative semi-discretizations are typically guaranteed to balance energy if and only if incompressibility is satisfied pointwise. This is especially alarming as both momentum conservation and energy balance play a critical role in flow structure development. Moreover, energy balance is inherently linked to the numerical stability of a method. In this dissertation, novel B-spline discretizations for the generalized Stokes and Navier-Stokes equations are developed. The cornerstone of this development is the construction of smooth generalizations of Raviart-Thomas-Nedelec elements based on the new theory of isogeometric discrete differential forms. The discretizations are (at least) patch-wise continuous and hence can be directly utilized in the Galerkin solution of viscous flows for single-patch configurations. When applied to incompressible flows, the discretizations produce pointwise divergence-free velocity fields. This results in methods which properly balance both momentum and energy at the semi-discrete level. In the presence of multi-patch geometries or no-slip walls, the discontinuous Galerkin framework can be invoked to enforce tangential continuity without upsetting the conservation and stability properties of the method across patch boundaries. This also allows our method to default to a compatible discretization of Darcy or Euler flow in the limit of vanishing viscosity. These attributes in conjunction with the local stability properties and resolution power of B-splines make these discretizations an attractive candidate for reliable numerical simulation of viscous incompressible flows. / text

In this thesis, evaluation of H(div)-conforming finite elements is implemented in a high-performance setting and used to solve the incompressible Navier-Stokes equation, obtaining an exactly point-wise divergence-free velocity field. In particular, the anisotropic Raviart-Thomas tensor-product polynomial space is considered, where the finite element operators are evaluated with quadrature in a matrix-free fashion using sum-factorization on tensor-product meshes. The implementation includes evaluation over elements and faces in two- and three-dimensional space, supporting non-conforming meshes with hanging nodes, and using the contravariant Piola transformation to preserve normal components on element boundaries. In terms of throughput, the implementation achieves up to an order of magnitude faster evaluation of finite element operators compared to a matrix-based evaluation. Correctness is demonstrated with optimal convergence rates for various polynomial degrees, as well as exactly divergence-free solutions for the velocity field.

First, we provide new classes of initial data, that grant short time uniqueness of the associated weak Leray-Hopf solutions of the three dimensional Navier-Stokes equations. The main novelty here is the establishment of certain continuity properties near the initial time, for weak Leray-Hopf solutions with initial data in supercritical Besov spaces. The techniques used here build upon related ideas of Calderón. Secondly, we prove local regularity up to the at part of the boundary, for certain classes of solutions to the Navier-Stokes equations, provided that the velocity field belongs to L_∞(-1; 0; L^{3, β}(B(1) ⋂ ℝ³ ₊)) with 3 ≤ β < ∞. What enables us to build upon the work of Escauriaza, Seregin and Šverák [27] and Seregin [100] is the establishment of new scale-invariant estimates, new estimates for the pressure near the boundary and a convenient new ϵ-regularity criterion. Third, we show that if a weak Leray-Hopf solution in ℝ³ ₊×]0,∞[ has a finite blow-up time T, then necessarily lim_t↑T||v(·, t)||_{L^3,β(ℝ³ ₊)} = ∞ with 3 < β < ∞. The proof hinges on a rescaling procedure from Seregin's work [106], a new stability result for singular points on the boundary, suitable a priori estimates and a Liouville type theorem for parabolic operators developed by Escauriaza, Seregin and Šverák [27]. Finally, we investigate a notion of global-in-time solutions to the Navier- Stokes equations in ℝ³, with solenoidal initial data in the critical Besov space ?^-1/4_4,∞(ℝ³), which has certain continuity properties with respect to weak* convergence of the initial data. Such properties are motivated by the strategy used by Seregin [106] to show that if a weak Leray-Hopf solution in ℝ³×]0,∞[ has a finite blow-up time T, then necessarily lim_t↑T ||v(·, t)||_{L₃(ℝ³)} = ∞. We prove new decomposition results for Besov spaces, which are key in the conception and existence theory of such solutions.

In this report the regular finite differences method (FDM) and a least-squares radial basis function-generated finite differences method (RBF-FD-LS) is used to solve the two-dimensional incompressible Navier-Stokes equations for the lid driven cavity problem. The Navier-Stokes equations is solved using stream function-vorticity formulation. The purpose of the report is to compare FDM and RBF-FD-LS with respect to accuracy and computational cost. Both methods were implemented in MATLAB and the problem was solved for Reynolds numbers equal to 100, 400 and 1000. In the report we present the solutions obtained as well as the results from the comparison. The results are discussed and conclusions are drawn. We came to the conclusion that RBF-FD-LS is more accurate when the stepsize of the grids used is held constant, while RBF-FD-LS costs more than FDM for similar accuracy.

Dans cette thèse, nous étudions des modèles mathématiques concernant certains problèmes d'écoulement de fluide à densité variable. Le premier chapitre résume l'ensemble de la thèse et se concentre sur les résultats obtenus, la nouveauté et la comparaison avec la littérature existante. Dans le deuxième chapitre, nous étudions la stabilisation locale des équations non homogènes de Navier-Stokes dans un canal 2d autour du flot de Poiseuille. Nous concevons un contrôle feedback de la vitesse qui agit sur l'entrée du domaine de sorte que la vitesse et la densité du fluide soient stabilisées autour du flot de Poiseuille, à condition que la densité initiale soit donnée par une constante additionnée d'une perturbation dont le support se situe loin du bord latéral du canal. Dans le troisième chapitre, nous étudions un système couplant les équations de Navier-Stokes compressibles à une structure élastique située à la frontière du domaine fluide. Nous prouvons l'existence locale de solutions solides pour ce système couplé. Dans le quatrième chapitre, notre objectif est d'étudier la nulle- contrôlabilité d'un problemè d'interaction fluide-structure linéarisé dans un canal bi dimensional. L'écoulement du fluide est ici modélisé par les équations de Navier-Stokes compressibles. En ce qui concerne la structure, nous considérons une poutre de type Euler-Bernoulli amortie située sur une partie du bord. Dans ce chapitre, nous établissons une inégalité d'observabilité pour le problème considéré d'interaction fluid-structure linéarisé qui constitue le premier pas vers la preuve de la nulle contrôlabilité du système. / In this thesis we study mathematical models concerning some fluid flow problems with variable density. The first chapter is a summary of the entire thesis and focuses on the results obtained, novelty and comparison with the existing literature. In the second chapter we study the local stabilization of the non-homogeneous Navier-Stokes equations in a 2d channel around Poiseuille flow. We design a feedback control of the velocity which acts on the inflow boundary of the domain such that both the fluid velocity and density are stabilized around Poiseuille flow provided the initial density is given by a constant added with a perturbation, such that the perturbation is supported away from the lateral boundary of the channel. In the third chapter we prove the local in time existence of strong solutions for a system coupling the compressible Navier-Stokes equations with an elastic structure located at the boundary of the fluid domain. In the fourth chapter our objective is to study the null controllability of a linearized compressible fluid structure interaction problem in a 2d channel where the structure is elastic and located at the fluid boundary. In this chapter we establish an observability inequality for the linearized fluid structure interaction problem under consideration which is the first step towards the direction of proving the null controllability of the system.

Chorin’s method, Incremental Pressure Correction Scheme (IPCS) and Crank-Nicolson’s method (CN) are three numerical methods that were investigated in this study. These methods were here used for solving the incompressible Navier-Stokes equations, which describe the motion of an incompressible fluid, in three different benchmark problems. The methods were stabilized using residual based artificial viscosity, which was introduced to avoid instability. The methods were compared in terms of accuracy and computational time. Furthermore, a theoretical study of adaptivity was made, based on an a posteriori error estimate and an adjoint problem. The implementation of the adaptivity is left for future studies. In this study we consider the following three well-known benchmark problems: laminar 2D flow around a cylinder, Taylor-Green vortex and lid-driven cavity problem. The difference of the computational time for the three methods were in general relatively small and differed depending on which problem that was investigated. Furthermore the accuracy of the methods also differed in the benchmark problems, but in general Crank-Nicolson’s method gave less accurate results. Moreover the stabilization technique worked well when the kinematic viscosity of the fluid was relatively low, since it managed to stabilize the numerical methods. In general the solution was affected in a negative way when the problem could be solved without stabilization for higher viscosities.

Numerical solution of differential equations having multitude of scales in the solution field is one of the most challenging research areas, but highly demanded in scientific and industrial applications. One of the natural approaches for handling such problems is to separate the scales and approximate the solution of the segregated scales with appropriate numerical method. Variational multiscale method (VMS) is a predominant method in the paradigm of scale separation schemes. In our work we have used the VMS technique to develop a numerical scheme for computations of turbulent flows in time-dependent domains. VMS allows separation of the entire range of scales in the flow field into two or three groups, thereby enabling a different numerical treatment for the different groups. In the context of computational fluid dynamics(CFD), VMS is a significant new improvement over the classical large eddy simulation (LES). VMS does away with the commutation errors arising due to filtering in LES. Further, in a three-scale VMS approach the model for the subgrid scale can be contained to only a part of the resolved scales instead of effecting the entire range of resolved scales. The projection based VMS scheme that we have developed gives a robust and efficient method for solving problems of turbulent fluid flows in deforming domains, governed by incompressible Navier {Stokes equations. In addition to the existing challenges due to turbulence, the computational complexity of the problem increases further when the considered domain is time-dependent. In this work, we have used an arbitrary Lagrangian-Eulerian (ALE) based VMS scheme to account for the domain deformation. In the proposed scheme, the large scales are represented by an additional tensor valued space. The resolved large and small scales are computed in a single unified equation, and the effect of unresolved scales is confined only to the resolved small scales, by using a projection operator. The popular Smagorinsky eddy viscosity model is used to approximate the effects of unresolved scales. The used ALE approach consists of an elastic mesh update technique. Moreover, a computationally efficient scheme is obtained by the choice of orthogonal finite element basis function for the resolved large scales, which allows to reformulate the ALE-VMS system matrix into the standard form of the NSE system matrix. Thus, any existing Navier{Stokes solver can be utilized for this scheme, with modifications. Further, the stability and error estimates of the scheme using a linear model of the NSE are also derived. Finally, the proposed scheme has been validated by a number of numerical examples over a wide range of problems.

No description available.

Dans cette thèse, on s'intéresse à deux problèmes provenant de l'étude mathématique des fluides incompressibles visqueux : la propagation de la régularité tangentielle et le mouvement d'une surface libre.La première question concerne plus particulièrement l'étude qualitative de l'évolution de quantités thermodynamiques telles que la température dans l'équation de Boussinesq sans diffusion et la densité dans le système de Navier-Stokes non homogène. Typiquement, on suppose que ces deux quantités sont, à l'instant initial, discontinues le long d'une interface à régularité h"oldérienne. Comme conséquence de résultats de propagation de régularité tangentielle pour le champ de vitesses, on établit que la régularité des interfaces persiste pour tout temps aussi bien en dimension deux d'espace, qu'en dimension supérieure (avec condition de petitesse). Notre approche suit celle du travail de J.-Y. Chemin dans les années 90 pour le problème des poches de tourbillon dans les fluides incompressiblesparfaits.Dans le cas présent, outre cette hypothèse de régularité tangentielle, nous n'avons besoin que d'une régularité critique sur le champ de vitesses.La démonstration repose sur le calcul para-différentiel et les espaces de multiplicateurs.Dans la dernière partie de la thèse, on considère le problème à frontière libre pour le système de Navier-Stokes incompressible à deux phases. Ce système permet de décrire l'évolution d'un mélange de deux fluides non miscibles tels que l'huile et l'eau par exemple. Différents cas de figure sont étudiés : le cas d'un réservoir borné, d'une goutte ou d'une rivière à profondeur finie.On établit l'existence et l'unicité à temps petit pour ce problème. Notre démonstration repose fortement sur des propriétés de régularité maximale parabolique de type $L_p$-$L_q / This thesis is dedicated to two different problems in the mathematical study of the viscous incompressible fluids: the persistence of tangential regularity and the motion of a free surface.The first problem concerns the study of the qualitative properties of some thermodynamical quantities in incompressible fluid models, such as the temperature for Boussinesq system with no diffusion and the density for the non-homogeneous Navier-Stokes system. Typically, we assume those two quantities to be initially piecewise constant along an interface with H"older regularity.As a consequence of stability of certain directional smoothness of the velocity field, we establish that the regularity of the interfaces persist globally with respect to time both in the two dimensional and higher dimensional cases (under some smallness condition). Our strategy is borrowed from the pioneering works by J.-Y.Chemin in 1990s on the vortex patch problem for ideal fluids.Let us emphasize that, apart from the directional regularity, we only impose rough (critical) regularity on the velocity field. The proof requires tools from para-differential calculus and multiplier space theory.In the last part of this thesis, we are concerned with the free boundary value problem for two-phase density-dependent Navier-Stokes system.This model is used to describe the motion of two immiscible liquids, like the oil and the water. Such mixture may occur in different situations, such as in a fixed bounded container, in a moving bounded droplet or in a river with finite depth. We establish the short time well-posedness for this problem. Our result strongly relies on the $L_p$-$L_q$ maximal regularity theoryfor parabolic equations