Ce travail de thèse est consacré à deux thèmes de recherche en Calcul Scientifique liés par l'approximation numérique de problèmes en mécanique des fluides. Le premier thème concerne l'approximation numérique des équations de Stokes, modélisant les écoulements de fluides incompressibles à vitesse faible. Ce thème est présent dans plusieurs travaux en Calcum Scientifique. La discrétisation en temps est réalisée à l'aide de la méthode de projection. La discrétisation en espace utilise la méthode des éléments finis hybrides qui permet d'imposer de façon exacte la contrainte d'incompressibilité. Cette approche est originale : la méthode des éléments mixtes hybrides est couplée avec une méthode d'éléments finis standards. L'ordre de convergence des deux méthodes est préservé. Le second thème concerne la mise au point de méthodes numériques de type volumes finis pour la résolution de l'équation Level Set. Ces équations interviennent de manière essentielle dans la résolution des problèmes de propagation d'interfaces. Dans cette partie, nous avons développé une nouvelle méthode d'ordre 2 de type MUSCL pour résoudre le système hyperbolique résultant de l'équation Level Set. Nous illustrons ces propriétés par des applications numériques. En particulier nous avons regardé le cas du problème des deux demi-plans pour lequel notre schéma donne une approximation pour le gradient de la fonction Level Set. Par ailleurs, l'ordre de précision attendu est obtenu avec les normes L1 et Linfini pour des fonctions régulières. Pour finir, il est à noter que notre méthode peut être facilement étendue aux problèmes d'Hamilton-Jacobi du premier et du second ordre

No description available.

A Petrov-Galerkin method for the solution of the compressible Euler and Navier-Stokes equations is presented. The method is based on the introduction of an anisotropic balancing diffusion in the local direction of the propogation of the scalar variables. The direction in which the diffusion is added and its magnitude are automatically calculated locally using a criterion that is optimal for one-dimensional transport equations. Algorithms are developed using bilinear quadrilateral and linear triangular elements. The triangular elements are used in conjunction with an adaptive scheme using unstructured meshes. Several applications are presented that show the exceptional stability and accuracy of the method, including the ram accelerator concept for the acceleration of projectiles to ultrahigh velocities. Both two-dimensional and axisymmetric models are employed to evaluate multiple projectile configurations and flow conditions.

This report describes the development of a software for computing viscous incompressiblemultiphase ows. The software does this with solving the coupled non-linear Navier-Stokes(Fluid) and the Cahn-Hilliard (Phase-Field) equations using a Finite Element Method. Thereason for the development is to produce a simulation tool, which eventually is capable ofsimulating the ow of uids inside the OptusAir aeronator manufactured by the Sorubincompany. The solving software developed is built on the ParMetis, PETSc and OpenMPIframeworks. Our primary benchmark has been a geometry resembling the OptusAir product,a cylinder with a rotating bottom. We have made comparisons between simulation resultsand the theory of a free surface in a uniform rotating ow.This thesis shows that the shape of the interface between two uids coincides with theoryto some extent, while the approximate boundary conditions prevent it from coinciding fully.

For the latter half of the last century renormalization methods have played an important part in tackling problems in fundamental physics and in providing a deeper understanding of systems with many interacting scales or degrees of freedom with strong coupling. The study of turbulence is no exception, and this thesis presents an investigation of renormalization techniques available in the study of the statistical theory of homogeneous and isotropic turbulence. The thesis consists of two parts which assess the two main renormalization approaches available in modeling turbulence. In particular we will be focusing on the renormalization procedures developed by McComb and others. The first part of this thesis will discuss Renormalization Group (RG) approaches to turbulence, with a focus on applications to reduce the degrees of freedom in a large-eddy simulation. The RG methods as applied to classical dynamical systems will be reviewed in the context of the Navier-Stokes equations describing fluid flow. This will be followed by introducing a functional based formalism of a conditional average first introduced by McComb, Roberts and Watt [Phys. Rev A 45, 3507 (1992)] as a tool for averaging out degrees of freedom needed in an RG calculation. This conditional average is then used in a formal RG calculation applied to the Navier-Stokes equations, originally done by McComb and Watt [Phys. Rev. A 46, 4797 (1992)], and later revised by Mc- Comb and Johnston [Physica A 292, 346 (2001)]. A correction to the summing of the time-integral detailed in the latter work is shown to introduce an extra viscous life-time term to the denominator of the increment to the renormalized viscosity and is shown to have a negligible effect in the numerical calculations. We follow this study by outlining some problems with the previous approach. In particular it is shown that a cross-term representing the interaction between high and low wavenumber modes which was neglected in the previous studies on the grounds that it does not contribute to energy dissipation, does in fact contribute significantly. A heuristic method is then put forward to include the effects of this term in the RG calculation. This leads to results which agree qualitatively with numerical calculations of eddy-viscosities. We finish this part of the thesis with an application of the RG method to the modeling of a passive scalar advected by a turbulent velocity field. The second part of this thesis will begin by reviewing Eulerian renormalized perturbation theory attempts in closing the infinite moment hierarchy introduced by averaging the Navier-Stokes equations. This is followed by presenting a new formulation of the local energy transfer theory (LET) of McComb et. al. [J. Fluid Mech. 245, 279 (1992)] which resolves some problems of previous derivations. In particular we show by the introduction of time-ordering that some previous problems with the exponential representation of the correlator can be overcome. Furthermore, we show that the singularity in the LET propagator equation cancels by way of a counter-term. We end this study by introducing a single-time Markovian closure based on LET which, unlike other Markovian closures, does not rely on any arbitrary parameters being introduced in the theory.

A dissertation submitted to the Faculty of Science, University of the Witwatersrand, Johannesburg, in fulfilment of requirements for the degree of Master of Science. February 16, 2015. / The Adomian decomposition method is an effective procedure for the analytical solution of a wide class of dynamical systems without linearization or weak nonlinearity assumptions, closure approximations, perturbation theory, or restrictive assumptions on stochasticity. Our aim here is to apply the Adomian decomposition method to steady two-dimensional blood flow through a constricted artery in the presence of a uniform transverse magnetic field. Blood flow is the study of measuring blood pressure and determining flow through arteries. Blood flow is assumed to be Newtonian and is governed by the equation of continuity and the momentum balanced equation (which are known as the Navier-Stokes equations). This model is consistent with the principles of ferro-hydrodynamics and magnetohydrodynamics and takes into account both magnetization and electrical conductivity of blood. We apply the Adomian decomposition method to the equations governing blood flow through arteries in the presence of an external transverse magnetic field. The results show that the e ect of a uniform external transverse magnetic field applied to blood flow through arteries favors the physiological condition of blood. The motion of blood in stenosed arteries can be regulated by applying a magnetic field externally and increasing/decreasing the intensity of the applied field.

A finite difference based solution method is derived for the velocity-pressure formulation of the two-dimensional incompressible Navier-Stokes equations. The method is proven stable using the energy method, facilitated by SBP operators, for characteristic and Dirichlet boundary condition implemented using the SAT technique. The numerical experiments show the utility of high-order finite difference methods as well as emphasize the problem of pressure boundary conditions. Furthermore, we demonstrate that a discretely divergence free solution can be obtained by use of the projection method.

Uma grande limitação dos métodos de diferenças finitas é que eles estão restritos a malhas e domínios retangulares. Para descrever escoamentos em domínios complexos, como, por exemplo, problemas com superfícies livres, faz-se necessário o uso de técnicas acessórias. O método de interfaces imersas é uma dessas técnicas. Nesse trabalho, primeiramente foi desenvolvido um método de projeção, totalmente livre de pressão, para as equações de Navier-Stokes com variáveis primitivas em malha deslocada. Esse método é baseado em diferenças finitas compactas, possuindo segunda ordem temporal e quarta ordem espacial. Esse método foi combinado com o método de interface imersa de Linnick e Fasel [2] para resolver numericamente as equações de Stokes com quarta ordem de precisão. A verificação do código foi feita por meio do método das soluções manufaturadas e da comparação com resultados de outros autores em problemas clássicos da literatura. / A great limitation of finite differences methods is that they are restricted to retangular meshes and domains. In order to describe flows in complex domains, e.g. free surface problems, it is necessary to use accessory techniques. The immersed interface method is one of such techniques. In the present work, firstly, a projection method was developed, which is completely pressure-free, for the Navier-Stokes equations with primitive variables in a staggered mesh. This method is based on compact finite differences, with temporal second-order precision and spatial foruth-order precision. This method was combined with the immersed interface method from Linnick e Fasel [2] in order to numerically solve the Stokes equations with fourth-order precision. The verification of the code was performed with the manufactured solutions method and by comparing results with other authors for some classical problems in the literature.

In this thesis, we study a two fluid system which describes the motion of two charged particles in a strict neutral incompressible plasma. We study the well-posdness of the system in both space dimensions two and three. Regardless of the size of the initial data, we prove the global well-posedness of the Cauchy problem when the space dimension is two. In space dimension three, we construct global weak-solutions, and we prove the local well-posedness of Kato-type solutions. These solutions turn out to be global when the initial data are sufficiently small. We also study the stability of the solution around zero given that the initial data is small and has sufficient regularity. It turns out that our system is a system of regularity-loss and the L2 norm of lower derivatives of the solution decays. At last, this two fluid system can be used to derive the classic MHD at least formally. Arsenio, Ibrahim and Masmoudi (2015) proved that the two fluid system converges to MHD under some constraints. We showed numerically that the two fluid system converges to MHD with no such constraint and found the approximate converge rate. / Graduate

Este trabalho consiste na simulação de escoamentos aerodinâmicos sobre configurações tipo foguete, para os casos de um corpo isolado e de uma geometria multicorpo, utilizando as equações de Navier-Stokes com média de Reynolds, com aproximação de camada fina. Estas equações são implementadas em diferenças finitas, utilizando o algoritmo implícito de fatorização aproximada de Beam e Warming. O fenômeno de turbulência é modelado com o modelo algébrico de viscosidade de vórtice, de duas camadas, de Baldwin e Lomax. Os resultados obtidos são comparados com dados experimentais e a concordância é muito boa. Contornos de pressão e densidade são também apresentados para se comprovar a validade física da simulação. Finalmente, apresentam-se conclusões e idéias para trabalhos futuros.