A Dissertation submitted to the Faculty of Science, University of the Witwatersrand, Johannesburg, for the degree of Master of Science. School of Computer Science and Applied Mathematics. November 7, 2016. / The equations for the flow of a thin lm of molten glass on a bath of molten tin are extended to the case in which the viscosity of the molten glass depends on the temperature. The continuity equation for an incompressible fluid, the Navier-Stokes equation and the energy balance equation are written in the lubrication (thin fluid lm) approximation. The kinematic boundary condition and the boundary conditions for the normal and tangential stress and the normal heat flux are derived on the upper and lower surfaces of the glass ribbon. It is found for the lubrication approximation that only one equation is obtained for four unknowns which are the two horizontal velocity components, the absolute temperature difference and the thickness of the molten glass rib- bon. The remaining three equations are obtained by taking the calculation to the next order in the square of the ratio of the thickness to length of the glass ribbon. The kinematic edge condition and the edge conditions for the normal and tangential stress and the normal heat flux are derived. The four edge conditions and the boundary conditions at the inlet and outlet give the boundary conditions for the four partial differential equations. It is not the aim of the dissertation to solve the boundary value problem which has been derived, either numerically or analytically. / LG2017

No description available.

No description available.

A generalized Conjugate Gradient like method is used to solve the linear systems of equations formed at each time-integration step of the unsteady, two-dimensional, compressible Navier-Stokes equations of fluid flow. The Navier-Stokes equations are cast in an implicit, upwind finite-volume, flux split formulation. Preconditioning techniques are employed with the Conjugate Gradient like method to enhance the stability and convergence rate of the overall iterative method. The superiority of the new solver is established by comparisons with a conventional Line GaussSeidel Relaxation (LGSR) solver. Comparisons are based on 'number of iterations required to converge to a steady-state solution' and 'total CPU time required for convergence'. Three test cases representing widely varying flow physics are chosen to investigate the performance of the solvers. Computational test results for very low speed (incompressible flow over a backward facing step at Mach 0.1), transonic flow (trailing edge flow in a transonic turbine cascade) and hypersonic flow (shockon- shock interactions on a cylindrical leading edge at Mach 6.0) are presented. For the 1vfach 0.1 case, speed-up factors of 30 (in terms of iterations) and 20 (in terms of CPU time) are found in favor of the new solver when compared with the LGSR solver. The corresponding speed-up factors for the transonic flow case are 20 and 18, respectively. The hypersonic case shows relatively lower speed-up factors of 5 and 4, respectively. This study reveals that preconditioning can greatly enhance the range of applicability and improve the performance of Conjugate Gradient like methods. / Ph. D.

An optimal shape control problem for the steady-state Navier-Stokes equations is considered from an analytical point of view. We examine a rather specific model problem dealing with 2-dimensional channel flow of incompressible viscous fluid: we wish to determine the shape of a bump on a part of the boundary in order to minimize the energy dissipation. To formulate the problem in a comprehensive manner, we study some properties of the Navier-Stokes equations. The penalty method is applied to relax the difficulty of dealing with incompressibility in conjunction with domain perturbations and regularity requirements for the solutions. The existence of optimal solutions for the penalized problem is presented. The computation of the shape gradient and its treatment plays central role in the shape sensitivity analysis. To describe the domain perturbation and to derive the shape gradient, we study the material derivative method and related shape calculus. The shape sensitivity analysis using the material derivative method and Lagrange multiplier technique is presented. The use of Lagrange multiplier techniques,from which an optimality system is derived, is justified by applying a method from functional analysis. Finite element discretizations for the domain and discretized description of the problem are given. We study finite element approximations for the weak penalized optimality system. To deal with inhomogeneous essential boundary condition, the framework of a Lagrange multiplier technique is applied. The split formulation decoupling the traction force from the velocity is proposed in conjunction with the penalized optimality system and optimal error estimates are derived. / Ph. D.

Accurate and efficient numerical algorithms for solving the three dimensional Navier Stokes equations with a generalized thermodynamic and chemistry model and a one equation turbulence model on structured and unstructured mesh topologies are presented. In the thermo-chemical modeling, particular attention is paid to the modeling of the chemical source terms, modeling of equilibrium thermodynamics, and the modeling of the non-equilibrium vibrational energy source terms. In this work, nonequilibrium thermo-chemical models are applied in the unstructured environment for the first time. A three-dimensional, second-order accurate k-exact reconstruction algorithm for the inviscid and viscous fluxes is presented. Several new methods for determining the stencil required for the inviscid and viscous k-exact reconstruction are discussed. A new simplified method for the computation of the viscous fluxes is also presented. Implementation of the one equation Spalart and Allmaras turbulence model is discussed. In particular, an new integral formulation is developed for this model which allows flux splitting to be applied to the resulting convective flux. Solutions for several test cases are presented to verify the solution algorithms discussed. For the thermo-chemical modeling, inviscid solutions to the three dimensional Aeroassist Flight Experiment vehicle and viscous solutions for the axi-symmetric Ram-II C are presented and compared to experimental data and/or published results. For the hypersonic AFE and Ram-II C solutions, focus is placed on the effects of the chemistry model in flows where ionization and dissociation are dominant characteristics of the flow field. Laminar and turbulent solutions over a flat plate are presented and compared to exact solutions and experimental data. Three dimensional higher order solutions using the k-exact reconstruction technique are presented for an analytic forebody. / Ph. D.

The process of enhancing an upwind finite volume, two-dimensional, thin layer Navier Stokes solver to achieve complete Navier Stokes solutions is described. The shear stress and heat flux contributions are identified and transformed to a generalized coordinate system. The metrics which result from the transformation have a geometrical interpretation in the finite volume formulation and are presented as supporting material. The additional terms which are neglected in the thin-layer approximations, are evaluated and discretized consistently with the finite volume method. Implicit linearizations are applied to the second derivatives tangent to the body surface; however, the cross derivatives are not linearized and are treated conservatively. Validation of the Navier Stokes solver is acquired by comparison to existing computational solutions for a double throat nozzle. Additional viscous solutions for the thin layer and the complete forms of the NS equations are provided for a flat plate shock boundary layer interaction. / Master of Science

The development of the patched grid solution methodology for both the Euler and the Navier-Stokes equations in two dimensions is presented. The governing equations are written in the integral form and the basic numerical algorithm is finite volume. The method is capable of first through third order accuracy in space. The flux vectors associated with the Euler equations are split into two sub-vectors (based on the signs of the characteristic speeds) and discretized separately. The viscous and heat flux contributions are treated with central differences. Patched grid results are demonstrated on shock reflection, subsonic boundary layer, and shock-boundary layer interaction flow problems. The results are compared with non-patched or single zone grids. The patched grid approach shows an improvement in resolution while minimizing storage and computer time. / Master of Science

La première partie de cette thèse traite du développement du code de simulation numérique directe YWC, principalement du développement des conditions aux limites. En effet, une forte contribution scientifique a été apportée aux conditions aux limites appelées "Three dimensional Navier-Stokes characteristic boundary condtions" (3D-NSCBC). Premièrement, la formulation de ces conditions aux arêtes et coins a été complétée, ensuite une extension de la formulation a été proposée pour supprimer les déformations observées en sortie dans le cas d'écoulements non-perpendiculaires à la frontière. <p>De plus, ces conditions ont été étendues au cas des gaz réels et une nouvelle définition du facteur de relaxation pour la pression a été proposée. Ce nouveau facteur de relaxation permet de supprimer les déformations observées en sortie pour des écoulements transcritiques. <p>Les résultats obtenus avec le code YWC ont ensuite été utilisés dans la seconde partie de la thèse pour développer une nouvelle méthode de tabulation basée sur l'analyse en composantes principales. Par rapport aux méthodes existante telles que FPI ou SLFM, la technique proposée, permet une identification automatique des variables à transporter et n'est, de plus, pas lié à un régime de combustion spécifique. Cette technique a permis d'effectuer des calculs d'interaction flamme-vortex en ne transportant que 5 espèces à la place des 9 requises pour le calcul en chimie détaillée complète, sans pour autant perdre en précision. <p>Finalement, dans le but de réduire encore le nombre d'espèces transportées, les techniques T-BAKED et HT-BAKED PCA ont été introduites. En utilisant une pondération des points sous-représentés, ces deux techniques permettent d'augmenter la précision de l'analyse par composantes principales dans le cadre des phénomènes de combustion.<p> / Doctorat en Sciences de l'ingénieur / info:eu-repo/semantics/nonPublished

Les équations fondamentales décrivant la dynamique de l'océan sont en théorie les équations de Navier-Stokes sur une sphère en rotation, auxquelles il faut a jouter une équation d'état pour la densité, et des équations de transport-diffusion pour les traceurs. Toutefois, un certain nombre de considérations physiques et de limitations pratiques ont nécessité le développement de modèles plus simples. En effet, un certain nombre d'hypothèses simplificatrices sont pleinement justifiées du point de vue de la physique des mouvements océaniques, dont les principales sont les approximations de couche mince et de Boussinesq. D'autre part, étant donné les dimensions des bassins océaniques (plusieurs centaines à plusieurs milliers de kilomètres), les coûts de calculs sont un facteur pratique extrêmement limitant. On est, à l'heure actuelle, capable de simuler l'océan mondial avec une résolution de l'ordre de dix kilomètres, en utilisant des modèles dits aux équations primitives, dont le coût de calcul est bien inférieur à celui des équations de Navier-Stokes. On est donc bien loin d'une modélisation complète des phénomènes décrits par ces équations, qui nécessiterait en théorie de considérer des échelles de l'ordre du millimètre. Les équations primitives sont issues des équations complètes de la mécanique des fluides en effectuant l'approximation hydrostatique, justifiée par la faible profondeur des domaines considérés au regard de leur dimension horizontale. Dans cette thèse, nous considérerons les équations de Navier-Stokes (NS) qui sont le coeur du modèle complet évoqué ci-dessus, sans prendre en compte les équations de la densité et des traceurs (salinité, température, etc.). Nous utiliserons l'approximation hydrostatique dans le chapitre 10, et le modèle sera naturellement appelé Navier-Stokes hydrostatique (NSH). Il correspond aux équations primitives dans lesquelles on ne prendrait pas en compte la densité et les traceurs. C'est dans ce cadre que se situe le travail présenté dans cette thèse, avec l'objectif à moyen terme de pouvoir coupler de façon rigoureuse et efficace les équations de Navier-Stokes avec les équations primitives. Dans une première partie, on présentera quelques rappels sur les équations de Navier-Stokes, leur discrétisation, ainsi que le cas-test de la cavité entraînée qui sera utilisé dans tout ce document. Dans une deuxième partie, on met en œuvre les méthodes de Schwarz sur les équations de Stokes et Navier-Stokes, en dérivant notamment des conditions absorbantes exactes et approchées pour ces systèmes. Enfin, dans une troisième partie, on proposera des pistes vers le couplage Navier-Stokes/Navier-Stokes hydrostatique décrit ci-dessus. / Fundamental equations describing the ocean dynamic are theoretically Navier-Stokes equations over a rotating sphere, whom need to add a state equation for the fluid density, and advection-diffusion equations for tracers. However, some physical considerations and practical limitations required to developped more simple models. Indeed, some simplifying hypotheses are well justified from a ocean dynamic point of view, whose principal ones are thin layer and Boussinesq approximations. On the other hand, considering the dimensions of oceans (from serveral hundreds to serveral thousands kilometers), computations costs are a very practical limitating factor. We are, by now, able to simulate the global ocean with about ten kilometers large grid mesh. This is very far from a complete modelisation of all phenomenes decribed by the Navier-Stokes equations, which require to consider scales of milimeters order. Primitives equations derive from complete equations describing fluid mecanics, by doing the hydrostatic approximations, which is justified by the low deepness of considered domains with regard to their horizontal dimension. In this thesis, we are considering Navier-Stokes equations (NS) which are the heart of the complete modele mentionned previously, without holding in account density and tracers equations. We will use the hydrostatic approximations, and the resulting equations will be named as hydrostatic Navier-Stokes equations (NSH).The mid term objective is to couple carefully Navier-Stokes equations with primitive equation. In a first part, we will remind few results for Navier-Stokes equations, their discretization, and the lid-driven cavity which wil be used as a test-case. In a second part, we will use Schwarz method with Stokes and Navier-Stokes equations, deriving in particular exact and approched absorbing interface conditions for these systems. Finally, in a third part, we shall propose first results towards coupling Navier-Stokes and hydrostatic Navier-Stokes equations.