Numerical algorithms for three dimensional computational fluid dynamic problems

Author

Mora Acosta, Josue

Director

Oliva, Asensio

Codirector

Soria Guerrero, Manel

Date of defense

2001-12-20

ISBN

8468919683

Legal Deposit

B-28695-2005



Department/Institute

Universitat Politècnica de Catalunya. Departament de Màquines i Motors Tèrmics

Abstract

The target of this work is to contribute to the enhancement of numerical methods for the simulation of complex thermal systems. Frequently, the factor that limits the accuracy of the simulations is the computing power: accurate simulations of complex devices require fine three-dimensional discretizations and the solution of large linear equation systems.<br/>Their efficient solution is one of the central aspects of this work. Low-cost parallel computers, for instance, PC clusters, are used to do so. The main bottle-neck of these computers is the notwork, that is too slow compared with their floating-point performance.<br/>Before considering linear solution algorithms, an overview of the mathematical models used and discretization techniques in staggered cartesian and cylindrical meshes is provided. <br/>The governing Navier-Stokes equations are solved using an implicit finite control volume method. Pressure-velocity coupling is solved with segregated approaches such as SIMPLEC.<br/>Different algorithms for the solution of the linear equation systems are reviewed: from incomplete factorizations such as MSIP, Krylov solvers such as BICGSTAB and GMRESR to acceleration techniques such as the Algebraic Multi Grid and the Multi Resolution Analysis with wavelts. Special attention is paid to preconditioned Krylov solvers for their application to parallel CFD problems.<br/>The fundamentals of parallel computing in distributed memory computers as well as implemetation details of these algorithms in combination with the domain decomposition method are given. Two different distributed memory computers, a Cray T3E and a PC cluster are used for several performance measures, including network throughput, performance of algebraic subroutines that affect to the overall efficiency of algorithms, and the solver performance. These measures are addressed to show the capabilities and drawbacks of parallel solvers for several processors and their partitioning configurations for a problem model.<br/>Finally, in order to illustrate the potential of the different techniques presented, a three-dimensional CFD problem is solved using a PC cluster. The numerical results obtained are validated by comparison with other authors. The speedup up to 12 processors is measured. An analysis of the computing time shows that, as expected, most of the computational effort is due to the pressure-correction equation,here solved with BiCGSTAB. The computing time algorithm , for different problem sizes, is compared with Schur-Complement and Multigrid.


El trabajo de tesis se centra en la solución numérica de las ecuaciones de navier-Stokes en regimen transitorio, tridimensional y laminar. Los algoritmos utilizados son del tipo segregado (SIMPLEC)y se basan en el uso de técnicas de volumenes finitos, con mallas estructurales del tipo staggered y discretizaciones temporales implícitas. En este contexto, el pricipal, problema son los elevados tiempos de cálculo de las simulaciones, que en buena parte se deben a la solución de los sistemas de ecuaciones lineales. Se hace una revisión de diferentes métodos utilizados típicamente en ordenadores secuenciales: GMRES, BICGSTAB, ACM, MSPIP. <br/>A fin de reducir los tiempos de cálculo se emplean ordenadores paralelos de memoria distribuida, basados en la agrupacion de ordenadores personales convencionales (PC clusters). Por lo que respecta a la potencia de cálculo por procesador, estos sistemas son comparables a los ordenadores paralelos de memoria distribuida convencionales (como el Cray T3E) siendo, su principal problema la baja capacidad de comunicación (elevada latencia, bajo ancho de banda). Este punto condiciona toda la estrategia computacional, obligando a reducir al máximo el número y el tamaño de los mensajes intercambiados. Este aspecto se cuantifica detalladamente en la tesis, realizando medidas de tiempos de cálculo en ambos ordenadores para diversas operaciones críticas para los algoritmos lineales. Tambien se miden y comparan los tiempos de cálculo y speed ups obtenidos en la solución de los sistemas lineales con diferentes algoritmos paralelos (Jacobi, MSIP, GMRES, BICGSTAB) y para diferentes tamaños de malla. <br/>Finalmente, se utilizan las técnicas anteriores para resolver el caso denominado driven cavity, en situacionies tridimensionales y con numeros de Reynolds de hasta 8000. Los resultados obtenidos se utilizan para validar los códigos desarrollados, en base a resultados de otros códigos y también se basa en la comparación con resultados experimentales procedentes de la bibliografía. Se utilizan hasta 12 procesadores, obteniendose spped ups de hasta 9.7 en el cluster de PCs. Se analizan los tiempos de cálculo de cada fase del código, señalandose areas para futuras mejoras. Se comparan los tiempos de cálculo con los algoritmos implementados en otros trabajos. La conclusión final es que los clusters de PCs son una plataforma de gran potencia en los cálculos de dinámica de fluidos computacional.

Keywords

fluid dynamic problems; numerical algorithms

Subjects

536 - Heat. Thermodynamics. Statistical physics

Knowledge Area

2204. Física de fluïds

Documents

01Jma01de01.pdf

3.817Mb

02Jma02de09.pdf

28.83Kb

03Jma03de09.pdf

81.24Kb

04Jma04de09.pdf

89.82Kb

05Jma05de09.pdf

214.0Kb

06Jma06de09.pdf

721.5Kb

07Jma07de09.pdf

1.223Mb

08Jma08de09.pdf

808.2Kb

09Jma09de09.pdf

286.9Kb

 

Rights

ADVERTIMENT. L'accés als continguts d'aquesta tesi doctoral i la seva utilització ha de respectar els drets de la persona autora. Pot ser utilitzada per a consulta o estudi personal, així com en activitats o materials d'investigació i docència en els termes establerts a l'art. 32 del Text Refós de la Llei de Propietat Intel·lectual (RDL 1/1996). Per altres utilitzacions es requereix l'autorització prèvia i expressa de la persona autora. En qualsevol cas, en la utilització dels seus continguts caldrà indicar de forma clara el nom i cognoms de la persona autora i el títol de la tesi doctoral. No s'autoritza la seva reproducció o altres formes d'explotació efectuades amb finalitats de lucre ni la seva comunicació pública des d'un lloc aliè al servei TDX. Tampoc s'autoritza la presentació del seu contingut en una finestra o marc aliè a TDX (framing). Aquesta reserva de drets afecta tant als continguts de la tesi com als seus resums i índexs.

This item appears in the following Collection(s)