Simulation of viscoelastic fluid flows in complex geometry at high Weissenburg (Wi) number is still a challenging problem in computational rheology. In this thesis, parallel computing toolbox available in OpenFOAM has been analysed in details. The scalability of parallel viscoelastic flow solver has been critically evaluated under benchmark flow problems, including 2D and 3D 4:1 contraction flow, 2D flow past a cylinder with 50% blockage, using up to 4096 CPU cores and 55 million computational grids. Areas for further improvements in parallel computational rheology are discussed. A new monitoring and preserving molecular conformation method is proposed to overcome the unphysical artefact problem in simulation of the FENE-CD-JS model under benchmark flow problems. It greatly enhances the stability of parallel viscoelastic flow solver in simulation of high Wi number flows and, for the first time, successfully captures elastic turbulence in flow past a cylinder problem. A new constitutive model, named as FENE-CD2-JS model, is proposed to overcome the existing shortcomings of the original FENE dumbbell model. The model accounts for high order interactions between non-equilibrium polymer chain molecules and could reproduce the asymmetric oscillatory vortex dynamics in the 16:1 contraction flow geometry at high Wi number (up to 48) flow observed in the experiments. For the first time, the analysis of the spatial distribution of non-equilibrium polymer conformation, through the conformational tensor, in strong complex flow and the results of their power-law scaling are also presented.

Please read the abstract in the section 00front of this document / Dissertation (M Eng (Mechanical Engineering))--University of Pretoria, 2007. / Mechanical and Aeronautical Engineering / unrestricted

The Richtmyer-Meshkov instability occurs when density-stratified interfaces are impulsively accelerated, typically by a shock wave. We present a numerical method to simulate the Richtmyer-Meshkov instability in cylindrical geometry. The ideal MHD equations are linearized about a time-dependent base state to yield linear partial differential equations governing the perturbed quantities. Convergence tests demonstrate that second order accuracy is achieved for smooth flows, and the order of accuracy is between first and second order for flows with discontinuities. Numerical results are presented for cases of interfaces with positive Atwood number and purely azimuthal perturbations. In hydrodynamics, the Richtmyer-Meshkov instability growth of perturbations is followed by a Rayleigh-Taylor growth phase. In MHD, numerical results indicate that the perturbations can be suppressed for sufficiently large perturbation wavenumbers and magnetic fields.

Orientador: Mauro Augusto Demarzo / Tese (doutorado) - Universidade Estadual de Campinas, Faculdade de Engenharia Civil, Arquitetura e Urbanismo / Made available in DSpace on 2018-08-07T08:31:00Z (GMT). No. of bitstreams: 1 Tacitano_Marcelo_D.pdf: 8426940 bytes, checksum: 2052cac2233a599ef85b5133ca5ceec4 (MD5) Previous issue date: 2006 / Resumo: As estruturas de contenção de valas, normalmente de madeira, aço ou concreto, podem ser constituídas por estacas-prancha, estacas com pranchões, estacões, paredes-diafragma entre outras e precisam ser dimensionadas, de que maneira que atinjam adequadamente os requisitos de funcionalidade, exeqüibilidade, segurança e economia a que se propõem. Dentre os métodos para o seu dimensionamento é possível enquadrá-los em três grandes grupos. Os Métodos Empíricos que se baseiam em resultados de medidas experimentais, os Métodos Semi-Empíricos que admitem como carregamento um diagrama de pressões para ambos os lados da parede, em cada fase de escavação, pressupondo o tipo de grandeza dos deslocamentos e considerando as estroncas e tirantes como apoios fixos e, finalmente, os Métodos Analíticos, que levam em conta as características de resistência e rigidez da estrutura e do maciço e possibilitam o cálculo evolutivo em que os esforços e deslocamentos das fases anteriores são efetivamente levados em conta nos cálculos das fases seguintes. Inicialmente uma ampla revisão bibliográfica sobre os métodos de cálculo de paredes de contenções é apresentada. Após, este trabalho adota como foco de estudo o Modelo de Winkler, através de Método Analítico Unidimensional, que tem sua aplicação prática pela construção do programa CEDEVE (Cálculo Evolutivo de Deslocamentos e Esforços em Valas Escoradas). Este método assimila a parede como uma viga de largura unitária, sendo o solo modelado como molas de comportamento elasto-plástico perfeito incluindo histerese. Estroncas e tirantes, de comportamento elástico, com ou sem esforços iniciais, podem ser introduzidos na estrutura. As ações sobre a estrutura advêm dos empuxos de solo, de água e das eventuais sobrecargas presentes na superfície. Os cálculos são conduzidos de acordo com as fases de escavação, retirando-se as ações (empuxos) e molas correspondentes ao solo escavado e introduzindo-se as estroncas e/ou tirantes a serem instalados, sendo que os esforços e deslocamentos ocorridos nas fases anteriores são devidamente considerados nos cálculos das fases seguintes. O reaterro, quando existir, também é considerado. Um diferencial importante do método de cálculo proposto com relação a outros similares é a possibilidade da inclusão dos efeitos de temperatura nos cálculos dos deslocamentos e esforços, sendo que tais efeitos térmicos podem ser considerados sobre as estroncas, o que promove a tendência de seu alongamento, que sendo parcialmente impedido, gera conseqüentes esforços de compressão e também na própria parede de contenção através de gradientes que induzem flexão, e, portanto também influindo nos esforços sobre o sistema de estroncas. Um estudo prático e numérico é conduzido com o intuito de se validar e verificar o programa CEDEVE, primeiramente testando seus resultados com o programa SAP 2000 e, após, comparando os resultados por ele gerados com vários outros programas disponíveis (SPW2003, DEEP e ESTWIN). O efeito da temperatura sobre as estroncas, calculado pelo CEDEVE, é comparado com alguns resultados de instrumentações disponíveis na bibliografia consultada. Além disso, um estudo comparativo com alguns Métodos Empíricos e Semi-Empíricos é conduzido. De uma forma geral, é possível concluir que o Modelo de Winkler utilizado na modelagem do problema gera resultados satisfatórios e sua relação custo benefício é bastante atraente na análise de paredes de contenção / Abstract: The retaining structures, usually of wood, steel or concrete, can be constituted by sheet pile wall, piles with lagging, ¿in cast¿ pile walls and diaphragm walls among others and need to be designed, so that they reach the requirements about functionality, execution, safety and economy an appropriate way. It is possible divide the design methods in three great groups. The Empirical Methods that has been based on results of experimental measures, the Semi-Empirical Methods that admit as loading a diagram of earth pressures for both sides of the wall, in each excavation phase, presupposing the displacements and considering the struts and anchorage as fixed supports and, finally, the Analytical Methods that take into account the characteristics of strength and stiffness of the structure and soil and they make possible the evolutionary calculation of internal efforts (strut forces, bending moments and shear forces) as well as the displacements, so the previous phases are taken indeed into account in the calculations of the following phases. Initially a wide bibliographical revision on the methods of calculation of retaining structures is presented. After that, this work adopts as focus the Winkler¿s Model, through One-Dimensional Analytical Method that bases the development of the CEDEVE program (Evolutionary Calculation of Displacements and Efforts in Braced Trenches). This program assimilates the wall as a beam of unitary width, being the soil modeled as springs with linear perfectly elastic-plastic behavior including histeresis. Struts and anchorages, of elastic behavior, with or without initials forces, can be introduced in the structure. The actions on the structure occur by the soil pressures, water pressures and eventually overloads in the soil surface. The calculations are performed in agreement with the excavation phases, leaving the actions (soil pressures) and springs corresponding to the dug soil as well as introducing the struts/anchorages, so that the efforts and displacements happened in the previous phases are properly considered in the calculations of the following phases. The process of cover the trench with earth, when it exists, is also considered. A important differential of the CEDEVE program regarding other similar ones is the possibility of the inclusion of the temperature effects in the calculations of the displacements and efforts, and such thermal effects can be considered on the struts, what promotes the tendency of its stretch, that being impeded partially, generates consequent compression and also in the own retaining wall through gradients that induce bending and, therefore, also influencing on the loads on the struts system. A practical and numerical study is done with the intention of validate and verify the CEDEVE program, firstly testing their results with the program SAP 2000 and, after that, comparing the results generated with several other available programs (SPW2003, DEEP and ESTWIN). The effect of the temperature on the struts, calculated by CEDEVE, is compared with some results of available instrumentations in the consulted bibliography. Besides this, a comparative study with some Empirical and Semi-Empirical Methods was done. In general, it is possible to conclude that the Winkler¿s Model used in the modeling of the problem generates satisfactory results and its relationship cost benefit is quite attractive in the analysis of retaining structures / Doutorado / Estruturas / Mestre em Engenharia Civil