O estudo de modelos matemáticos para o cálculo de desempenho de rotores, que começou no início do século 20 e se desenvolveu até meados daquele século, com Theodorsen, estagnou-se a partir da década de 1950, quando o uso de turboreatores se difundiu na propulsão aeronáutica. Recentemente a necessidade de desenvolvimento de sistemas propulsivos de alta eficiência para aeronaves de asas fixas e rotativas motiva um estudo mais detalhado sobre a aerodinâmica de rotores. As ferramentas utilizadas para projeto e análise de hélices, entretanto, ainda datam de décadas atrás. O presente estudo visa a comparar dois métodos de análise aerodinâmica de rotores. O primeiro é a formulação clássica da teoria do elemento de pá desenvolvida por Theodorsen. O segundo é um código baseado num método de singularidades tridimensional, segundo a formulação de Morino. Neste trabalho é desenvolvido um código de painéis tridimensional não estacionário de baixa-ordem, com distribuições constantes de dipolo e fonte sobre painéis quadrilaterais ou triangulares. O código é não-estacionário, e possui um algoritmo de relaxação da esteira livre no tempo, e dispõe de um método integral para cálculo de camada-limite e forças viscosas sobre a superfície dos painéis. O código também utiliza o recurso de processamento paralelo para acelerar a montagem das matrizes de coeficientes de influência e os cálculos de relaxação da esteira. Dois estudos de caso são feitos, simulando experimentos com rotores através dos dois métodos supracitados, e os resultados são comparados com dados experimentais encontrados na literatura. Dois outros estudos de caso são realizados para explorar a capacidade do método de painéis em modelar interferência entre o rotor e outras superfícies e condições de operação impossíveis de simular pela teoria clássica. Os resultados são analisados detalhadamente, e as vantagens e limitações de cada formulação são destacadas.

A previsão acurada de fenômenos aerodinâmicos não estacionários tem cada vez mais dependido do emprego de técnicas de Mecânica dos Fluidos Computacional (CFD). A aplicação de tais técnicas em problemas não estacionários depende crucialmente de tecnologias de geração de malha que, usualmente, são processos bastante lentos. Para contornar esta limitação, várias metodologias têm sido desenvolvidas a fim de calcular a dinâmica geométrica do problema físico sem gerar novamente uma malha computacional por completo. Neste sentido, este trabalho apresenta um estudo de técnicas eficientes de geração e movimentação de malhas computacionais para simulação não estacionária em CFD, com aplicações em Aeroelasticidade. Para este propósito, é implementada uma abordagem de movimentação de malhas baseada em Funções de Base Radial (RBF). Nesta metodologia, os deslocamentos dos pontos de malha na superfície da estrutura (móvel) são propagados para os pontos volumétricos (interiores), sem alterar a topologia da malha. Este trabalho conduz uma análise comparativa do emprego de diversos tipos de RBF na movimentação de malhas de níveis de refinamento diferentes, bidimensionais e tridimensionais, no que diz respeito à qualidade da malha deformada resultante. Além disso, aborda estratégias de movimentação e sua influência na qualidade de malha gerada. Como aplicação, é conduzida simulação aerodinâmica computacional com código BREXP3D. Este código utiliza formulação de mecânica dos fluidos não estacionária compressível e não viscosa (equações de Euler), implementada numericamente via método de diferenças finitas, o que exige malhas hexaédricas estruturadas. Este trabalho implementa condições de contorno para simulação de asa tridimensional em malhas monobloco e conduz simulação aerodinâmica estacionária em regime transônico, comparando com resultados experimentais. Problemas de malha e convergência não permitiram executar com sucesso a simulação aerodinâmica computacional não estacionária.

A presente dissertação trata da questão importante sobre dispersão no cálculo do coeficiente de arrasto associado a malhas computacionais em aplicações aeroespaciais. Este tema tem sido amplamente discutido na comunidade de Dinâmica dos Fluidos Computacional (CFD) na última década. Determinar com precisão o valor do coeficiente de arrasto de uma configuração sempre foi um dos maiores problemas em CFD na indústria aeroespacial. Tão importante quanto conhecer a grandeza que se deseja calcular é também conhecer quão preciso somos em determinar essa quantidade. Na indústria aeroespacial, programas computacionais baseados nas equações de Navier-Stokes com média de Reynolds (RANS) e malhas computacionais não estruturadas do tipo hexaédricas e do tipo tetraédricas com camadas de prisma para a camada limite aparecem como uma abordagem eficiente. Neste contexto, o presente trabalho conduz estudos de maneira a identificar e a entender as dispersões no cálculo do coeficiente de arrasto causadas por modificações na topologia da malha computacional, utilizando um único código de CFD. A fim de atingir tal objetivo, a asa ONERA-M6 foi escolhida como caso teste deste estudo por representar uma geometria simples e, assim, tornando mais clara a identificação do problema em questão. Com o objetivo de simplificar ainda mais o problema, as simulações não possuem arrasto induzido, que neste caso significa que o ângulo de ataque é zero para todos os casos. Dois números de Mach foram escolhidos para serem investigados neste estudo, Mach 0,3 e 0,8, por representarem regimes subsônico e transônico, respectivamente. Os estudos consideram um número de Reynolds igual a 3 milhões, que é tipicamente usado na indústria aeroespacial durante os estudos experimentais em túnel de vento. As dispersões no cálculo de arrasto encontradas neste estudo foram discutidas e investigadas com o objetivo de identificar as regiões da asa de maior contribuição da dispersão total e, também, apontar as diferenças na malha superficial e volumétrica que poderiam contribuir para tais dispersões. Baseado nos resultados obtidos neste trabalho, foi possível identificar a precisão no cálculo do coeficiente de arrasto de uma configuração simplificada, devido às modificações na topologia da malha que podem ser associadas a variações geométricas em um estudo de otimização aerodinâmica.

Este trabalho aborda a utilização de uma geometria representativa de um motor mono cilíndrico, com injeção direta de etanol hidratado e ignição assistida por centelha. A partir da utilização de um código comercial de Dinâmica dos Fluidos Computacional, a malha móvel é criada e as equações governantes são resolvidas. O escopo do trabalho inclui a caracterização do escoamento de ar durante a compressão e a formação de mistura estratificada. É analisado o efeito da geometria da câmara na formação da mesma e sua distribuição ao longo da câmara de combustão. É analisada a influência de parâmetros como: pressão do combustível injetado, início da injeção e massa de combustível injetada. O objetivo final é a obtenção da concentração de combustível ao longo do cilindro, durante o desenvolvimento do jato e o movimento do pistão, além da formação de filme nas paredes da câmara. Tais resultados são de extrema importância para a avaliação de estratégias de injeção que tornem possível a partida a frio de veículos em baixas temperaturas sem a aplicação de sistemas auxiliares. Os resultados, como já eram previstos, mostraram que a pressão e o início da injeção de combustível são parâmetros que influenciam diretamente a evaporação de combustível, assim como a quantidade de filme formado, principalmente no pistão.

Mecanismos de aprimoramentos passivos elevam a taxa transferência de calor em trocadores de calor por meio de modificações na superfície de troca térmica. Estas modificações geram turbulência no escoamento sem a necessidade de utilizar dispositivos que consomem energia externa. Este efeito limita o crescimento da camada limite térmica perto da superfície de transferência de calor. Neste trabalho a melhoria da superfície de troca térmica é exercida utilizando as aletas triangulares com venezianas. Este tipo de aleta promove incremento na troca térmica sem a necessidade de adição de material extra na sua fabricação ou qualquer modificação extrema no projeto do trocador de calor. Em virtude disso, o uso de aletas com venezianas em trocadores de calor utilizados no setor automotivo tem se mostrado eficiente para projetos onde há limitações de espaço e demanda por maior troca térmica. O trocador de calor compacto estudado é um intercooler utilizado para arrefecer o ar de admissão em veículos turboalimentados. Para efeito comparativo, foram então construídos dois modelos de intercoolers com dimensões geométricas idênticas, e com a mesma quantidade de tubo e aletas. Assim, um modelo foi fabricado com aleta triangular e o outro com aleta triangular com venezianas. A avaliação dos dois modelos de intercoolers, foram realizadas através de análise teórica e experimental. Os ensaios foram conduzidos em um banco de ensaio denominado calorímetro. Por meio dele os intercoolers com aletas triangulares e com aletas triangulares com venezianas tiveram suas taxas de dissipação de calor obtidas assim como as suas vazões e quedas de pressões em ambos os escoamentos. Os resultados experimentais foram comparados com os resultados obtidos por meio de correlações analíticas provenientes do método da Efetividade-NUT. As diferenças entre os resultados experimental e o teórico calculado pelo método ? - NUT variou entre 0,22% e 5,42% no escoamento interno e 0,06% e 9,92% no escoamento externo. Portanto, as comparações entre os resultados experimental e o calculado pelo método ? - NUT para ambos os modelos de intercooler apresentaram-se satisfatórios cuja a estimativa de incerteza de medição situou-se entre 1,43% e 4,06% para um nível de confiança de 95%. Os resultados obtidos foram apresentados em fator j de Colburn de modo a fornecer uma representação da transferência de calor em termos de parâmetro adimensional em função do número de Reynolds. Foi constatado que o fator j de Colburn na superfície com aleta triangular com venezianas é 2 vezes maior que o da superfície com aleta triangular. Contudo, o fator de atrito da superfície com aleta triangular com venezianas é 5 vezes maior que a superfície com a aleta triangular. Portanto houve um ganho na dissipação do calor com um aumento significativo no fator de atrito. Com o objetivo de avaliar a viabilidade da aleta com venezianas foram feitas análises simultâneas entre o fator j de Colburn e o fator de átrio de Fanning por meio do critério de avaliação de desempenho. Com tais ganhos obtidos e análises realizadas pode-se afirmar que a técnica desenvolvida por este trabalho de tese apresenta-se muito adequada para o desenvolvimento de trocadores de calor automotivos.

ABE, V.C. Determinação de parâmetros de operação de sistema de distribuição de ar frio pelo piso em ambientes de escritórios. 2007. 132 p. Dissertação (Mestrado) - Departamento de Engenharia de Construção Civil da Escola Politécnica da Universidade de São Paulo, São Paulo, 2007. Atualmente, um dos principais requisitos para que um ambiente de escritório atenda plenamente às necessidades de uso e ocupação é a flexibilidade, mas com conforto. O sistema de distribuição de ar frio pelo piso (Underfloor Air Distribution System - UFAD) está sendo utilizado nos edifícios de escritórios para garantir a flexibilização integrada do sistema com o mobiliário. Trata-se de um sistema que fornece ar frio proveniente de um plenum, através de difusores localizados nos painéis do piso elevado. Os difusores de ar podem ser ajustados, ou reposicionados pelos próprios usuários dos ambientes, possibilitando o controle das condições térmicas no entorno próximo, de acordo com suas preferências individuais. Embora o sistema de distribuição de ar pelo piso apresente diversas vantagens com relação ao sistema de distribuição de ar convencional, ainda há falta de informações objetivas e diretrizes de projeto padronizadas, havendo a necessidade de realização de mais pesquisas na área. A introdução da Dinâmica dos Fluidos Computacional (Computational Fluid Dynamics - CFD) apresentou uma alternativa para as pesquisas experimentais relacionadas ao estudo do escoamento do ar. Assim, para o desenvolvimento do presente estudo foram utilizados recursos de dinâmica dos fluidos computacional, com o uso do código FLUENT. O objetivo do presente trabalho é apresentar uma estratégia para a determinação dos parâmetros de operação de sistemas de distribuição de ar pelo piso aplicados a ambientes de escritórios. A estratégia baseiase no estabelecimento de inter-relações entre os principais parâmetros que definem o escoamento do ar em ambientes com o sistema de distribuição de ar pelo piso: a vazão de ar insuflado, a temperatura do ar insuflado e a diferença de pressão entre o plenum pressurizado e o interior do ambiente. O emprego da estratégia possibilita que o projetista identifique as diversas combinações entre os parâmetros de operação que resultem numa mesma condição desejada, permitindo assim a escolha da opção mais satisfatória. / ABE, V.C. Determination of operation parameters of rooms with underfloor air distribution system. 2007. 132 p. Dissertation (Master Course) - Departamento de Engenharia de Construção Civil da Escola Politécnica da Universidade de São Paulo, São Paulo, 2007. Nowadays, one of the main requirements of an office room to fully attend the use and occupation needs is the flexibility, but with comfort. The underfloor air distribution system is being used in the office buildings to guarantee the integrated flexibility of the system with the furniture layout. The system supplies cooled air to the room from a pressurized plenum, passing through air diffusers located in the raised floor panels. The air diffusers can be adjusted or relocated by the room users themselves, making possible the control of the thermal conditions in their surroundings, according to their individual preferences. Though the underfloor air distribution system presents a lot of advantages with regard to conventional overhead distribution systems, there is still a lack of objective information and standardized guidelines, and the necessity of accomplishment of more researches about the subject. The introduction of Computational Fluid Dynamics presented an alternative for the experimental tests related to the study of airflow in rooms. Thus, for the development of the present study a CFD tool was used, the code FLUENT. The purpose of the present work is to develop a strategy for the determination of the operation parameters of an underfloor air distribution system applied to office rooms. The strategy is based on the establishment of interrelations between the main parameters that define the airflow in rooms with the underfloor air distribution system: the input volume flux, the input air temperature and the pressure difference between underfloor plenum and the room. The use of this strategy enables the designer to identify several combinations of the operation parameters that result in one desired condition, allowing the choice of the most satisfactory option.

O sucesso da conclusão de um poço por perfuração e o custo do projeto estão relacionados às propriedades dos fluidos de perfuração. O estudo experimental foi realizado a partir da determinação e análise das curvas de tensão de cisalhamento versus taxa de deformação de suspensões aquosas de bentonita e carboximetilcelulose com diferentes concentrações de NaCl. Investigou-se, inicialmente, o comportamento reológico de suspensões de carboximetilcelulose 0,5% (em massa) com concentrações de sal, NaCl de 0 a 4% (em massa), para temperaturas variando de 14 a 26ºC. Os resultados experimentais possibilitam a caracterização reológica destas soluções como fluidos pseudoplásticos e indicam que a adição de NaCl em soluções com CMC alteram significativamente o comportamento reológico desse tipo de solução. Numa segunda etapa, estudou-se o comportamento reológico de suspensões de bentonita e carboximetilcelulose em diferentes concentrações de sal, NaCl. Empregaram-se soluções contendo 4,8 % de bentonita, 0,5 % de carboximetilcelulose e concentrações de 0 a 4% de NaCl. Constata-se o comportamento tixotrópico destas suspensões e a forte dependência com a concentração de sal. O estudo foi complementado com uma análise comparativa dos resultados do comportamento reológico obtidos por reômetro de cilindros coaxiais modelo Brookfield e o viscosímetro FANN 35 A, usualmente empregado na análise de fluidos de perfuração. O equacionamento do escoamento do fluido nos aparelhos permitiu uma interpretação mais detalhada dos valores descritos na norma, obtendo-se uma equivalência dos valores medidos em campo com os obtidos no reômetro Brookfield apesar de as faixas de operação de taxas de deformação serem distintas. / The properties of drilling fluids have a very significant effect on the successful of a well completion and the project costs. The experimental study was conducted based on the determination and analysis of the curves of shear stress and shear rate of aqueous suspensions composed of bentonite and carboxymetil cellulose at different concentrations of NaCl. The rheological behaviour of carboxymetil cellulose suspensions 0.5% w/w at NaCl concentrations in the range of 0 to 4% w/w from 14 to 26°C was studied. The experimental results allow the rheological characterization of these solutions as pseudoplastic fluid. It was observed that the addition of NaCl to CMC solutions changes its rheological behaviour significantly. In addition, the rheological behaviour of 4,8% bentonite suspensions with 0,5% carboxymethyl cellulose at different concentrations of NaCl in the range of 0 to 4% was studied. The characterization of these solutions showed thixotropic behaviour and a strong dependence on the salt concentration. The study was complemented with a comparative analysis of the results obtained by using a coaxial cylinders rheometer model Brookfield and the FANN Model 35 A viscometer, usually employed in drilling fluids analysis. The analytical solutions allowed a more detailed interpretation of the values described in the standard resulting in an equivalence between the measured values obtained in the field and in the Brookfield rheometer despite the differents shear rate ranges.

Diversas pesquisas têm demonstrado a importância dos combustíveis biorenováveis como alternativas eficientes na redução da demanda de energia fóssil, aquecimento global e seus efeitos ambientais. Os óleos de pinhão manso e de palma (dendê) com suas diferentes composições de ácidos graxos, estabilidade oxidativa e propriedades fluidodinâmicas, têm chamado a atenção no Peru e Brasil, quando o assunto é a diversificação das matérias-primas para a obtenção de um biodiesel otimizado a partir de misturas desses óleos, com desempenho melhorado e maior resistência oxidativa. Neste trabalho, biodiesel foi produzido pelos métodos de transesterificação convencional e transesterificação assistida por ultrassom. Foram estudadas cinco misturas binárias dos óleos brutos pinhão manso e palma pré-tratados (9,1; 18,2; 28,6; 37,0 e 41,7% de óleo de pinhão manso), utilizando etanol e hidróxido de sódio como catalisador. Para determinar as melhores condições de reação foi utilizado o delineamento experimental proposto pelo método de superfície de resposta tendo como variáveis independentes a relação entre as misturas binárias dos óleos de pinhão manso e palma, a razão molar óleo:etanol e a concentração do catalisador NaOH. As melhores misturas binárias dos óleos brutos foram compostas por 37% de pinhão manso e 63% de palma (1PM1,7P) e pela mistura de 28,6% de pinhão manso e 71,4% de palma (1PM2,5P), ambas propiciando altos rendimentos de ésteres etílicos (maiores que 96,6%) por ambos os métodos, com a diferença que o sistema por ultrassom só precisou de 4 minutos. Além disso, os valores de densidade, viscosidade, ponto de fluidez, ponto de entupimento de filtro a frio (1PM1,7P obtido por ultrassom) e a estabilidade oxidativa (1PM2,5P obtido por ultrassom) dos ésteres ficaram dentro da norma ANP (2014). A mistura binária dos óleos brutos de pinhão manso e palma produziram altos teores de ésteres etílicos com características fluidodinâmicas e de estabilidade oxidativa aprimoradas. / Several researches have demonstrated the importance of bio-renewable fuels as efficient alternatives in reducing the contribution of fossil energy demand, global warming and its environmental effects. Jatropha and palm and biodiesel oils that differ in their compositions of fatty acids, oxidative stability and fluid dynamics properties have caught attention in Brazil and Peru in the diversification of raw materials for production of an optimized biodiesel from mixtures of these oils, with improved performance and higher oxidative resistance. In this study, conventional and ultrasonic assisted transesterification was applied to five binary mixtures of the pretreated jatropha and palm crude oils (9.1, 18.2, 28.6, 37.0 and 41.7% of jatropha oil) using ethanol and sodium hydroxide as the catalyst. In order to determine the best reaction conditions, the experimental design proposed by the response surface method was used, having as independent variables the proportion of jatropha and palm oils in the binary mixtures, the molar oil: ethanol ratio and the NaOH catalyst concentration. The best binary oil mixtures were: 37% of jatropha and 63% of palm (1 PM1,7P) and 28.6% of jatropha and 71.4% of palm (1 PM2,5P), both providing high yields of ethyl esters (greater than 96.6%) by either systems, with the difference that the ultrasound system required only 4 minutes. In addition, density, viscosity, pour point, cold filter plug-in point for 1PM1.7P obtained by ultrasound and oxidative stability for 1PM2.5P by ultrasound of the esters were within the ANP standards (2014). The binary mixture of crude jatropha and palm oils produced high levels of ethyl esters with improved fluid dynamics and oxidative stability characteristics.

Neste trabalho, pretendeu-se obter as características da movimentação do ar no início do tempo de admissão em um motor diesel, partindo de sua entrada no cabeçote, passando pelos pórticos de admissão e válvulas, até a câmara de combustão, utilizando dados reais da geometria e do processo de sobrealimentação de um motor diesel. Para tanto, utilizou-se da abordagem clássica da dinâmica dos fluidos para o desenvolvimento de um modelo matemático constituído por equacionamentos que representassem o consumo de ar, a rotacionalidade do movimento do ar nos pórticos e a movimentação de ar no cilindro, descrita por linhas de corrente obtidas a partir da aplicação da teoria de escoamento potencial plano. Os cálculos foram efetuados usando planilhas eletrônicas geradas em aplicativo Excel e, tanto a avaliação dos efeitos da geometria dos condutos quanto a determinação das linhas de corrente, foram realizadas usando o aplicativo gráfico Rhinoceros. Neste trabalho foram estudadas seis condições associadas a seis ângulos de virabrequim do início do tempo de admissão e que geraram seis conjuntos de linhas de corrente. Os resultados foram comparados com ensaios de PIV - particle image velocimetry e apresentaram semelhanças relevantes nos formatos e nos sentidos das linhas de corrente. Concluiu-se que o procedimento utilizado, por sua simplicidade de implementação, não é prático quando há a necessidade de alteração de uma das variáveis envolvidas. Porém, ao desenvolver o equacionamento para os casos apresentados, foi possível estudar o fenômeno em profundidade satisfatória, mesmo com as hipóteses simplificadoras utilizadas. / In this work, the aim was to obtain the air flow characteristics in the beginning of the intake stroke in a diesel engine, considering the inlet in head, its passiage by intake runners and valves, up to the combustion chamber, using values obtained from geometry and supercharging process observed in a specific diesel engine. It was used a classical fluid dynamics approach to develop a mathematical model which represent the air consumption, the rotation of the air at intake runners and the air flow within the cylinder, described by streamlines obtained by the use of potential flow theory. The calculation used Excel electronic spreadsheets and both the effects from geometry of intake runners evaluation and the streamlines generation were conducted in drawing software Rhinoceros. In this work six different conditions, associated at six crankshaft angles in the beginning of the intake stroke were studied, resulting in six streamlines set. The results were compared to a PIV test - particle image velocimetry - and relevant similarities was observed both in shapes and rotation sense of streamlines. The conclusion was that the methodology, considering its simplicity in implementation, it is not the most practical method when a modification in variables involved is necessary. However, this methodology was enough to study the phenomenon in satisfactory depth, even using the simplifying assumptions.

O Tanque de Provas Numérico (TPN) é um laboratório pioneiro em hidrodinâmica aplicada e fruto de uma colaboração entre a indústria brasileira de petróleo (PETROBRAS S.A.) e as principais instituições de pesquisa do país. Seu principal objetivo é atuar como parceiro da indústria offshore e de petróleo, colaborando para a obtenção da autossuficiência da produção nacional de petróleo como uma poderosa ferramenta para projeto e análise de sistemas flutuantes de produção de óleo e gás. O coração do TPN é um cluster de computadores SMP que é hoje um dos maiores agrupamentos do Brasil para fins de pesquisa. Um dos focos de atenção do TPN é a aplicação do Moving Particles Semi-implicit Method (MPS) na exploração de soluções para muitos problemas de Engenharia. Por trabalhar sem a necessidade do uso de malhas (método tradicional Euleriano), o método tem diversas aplicações na simulação de corpos flutuantes e na hidrodinâmica aplicada, sendo atualmente utilizado para realização de estudos sobre a influência do movimento de ondas em navios; simulações de fenômenos que envolvem fragmentações; superfícies livres; grandes deformações; dinâmica de fluidos em condições extremas, como é o caso em processos prospecção do petróleo onde muitas vezes é difícil e economicamente inviável fazer ensaios físicos. Devido ao grande número de partículas utilizadas na simulação de sistemas complexos pelo método MPS, é necessário aproveitar de forma eficiente os recursos computacionais disponíveis para a análise de modelos com o refinamento adequado às aplicações práticas. Com tera-FLOPS disponíveis na rede cluster do TPN para modelagem computacional, há uma grande necessidade de uma solução computacional paralela altamente escalonável que, além disto, seja fácil de manutenção e extensibilidade. Dentro desta linha de pesquisa, foi desenvolvida uma solução com essas características através do emprego de modernas técnicas de engenharia de software. / The Numerical Offshore Tank (TPN) is a pioneer laboratory in applied hydrodynamics and result of collaboration between the Brazilian oil (Petrobras SA) and the major research institutions in the country. Its main purpose is to act as a partner of industry and offshore oil, contributing to the achievement of self-sufficiency of domestic oil production as a powerful tool for design and analysis of floating production systems for oil and gas. The heart of TPN is a cluster of SMP computers that is now one of the largest groupings of Brazil for research purposes. One focus of attention of TPN is the application of Moving Particles Semi-implicit Method (MPS) in exploring solutions to many engineering problems. By working without the use of mesh (Eulerian traditional method), the method has several applications in the simulation of floating bodies and applied hydrodynamics, currently being used for studies on the influence of the movement of ships in waves; simulations of phenomena involving fragmentation; free surfaces, large deformations; fluid dynamics in extreme conditions, as is the case in processes where petroleum exploration is often difficult and uneconomical to do physical tests. Due to the high number of particles used in the simulation of complex systems by the MPS method, it is necessary to efficiently take advantage of the computational resources available for the analysis of models with the refinement suitable for practical applications. With tera-FLOPS available in the TPN network cluster for computational modeling, there is a great need for a parallel highly scalable solution which, moreover, must be easy maintenance and extensibility. Within this line of research, we developed a solution with these characteristics through the use of modern software engineering techniques.