Modelos cinéticos da equação linearizada de Boltzmann e um problema de transferência de calor em microescala

Neste trabalho, um problema de transferência de calor da dinâmica de gases rarefeitos, causado pela diferença de temperaturas nas superfícies de um canal, é abordado. O problema é formulado através dos modelos cinéticos BGK, S e Gross-Jackson da equação linearizada de Boltzmann e resolvido, de forma unificada, pelo método analítico de ordenadas discretas (método ADO). Resultados numéricos para as perturbações de densidade e temperatura e também para o fluxo de calor são apresentados e comparados, mostrando que não se pode dizer que algum dos três modelos seja uma melhor aproximação da solução aos resultados da equação linearizada de Boltzmann.

A influência da organização do trabalho na comunicação em um escritório panorâmico : um estudo de caso

Com a evolução tecnológica e a crescente disputa no mercado, as empresas passaram a investir em novas tecnologias e métodos organizacionais que buscassem um aumento da produtividade e qualidade. Esta pesquisa, por meio de um estudo de caso, teve como objetivo levantar e analisar as condições ergonômicas na organização do escritório de uma empresa segundo tendências que contemplam redução de custo e rapidez na comunicação. Neste modelo, o escritório, identificado como aberto ou panorâmico, característica dos sistemas de gerenciamento corporativos, se, por um lado, apresentou vantagens, como a facilidade de comunicação e melhor controle da chefia sobre os trabalhadores, por outro lado, expôs também uma face negativa, porque normalmente tais atividades laborais exigem o exercício da concentração mental, e, em um escritório assim disposto, o trabalhador percebe os murmúrios desse ambiente como uma incômoda sensação que resulta em comprometimento da capacidade de concentração e queda tanto no rendimento quanto na produtividade. Tomou-se como ponto de partida uma queixa de perturbação que se instalou no setor sob a forma de ruído. Este estudo centrou sua investigação em três pontos: a percepção do trabalhador no ambiente organizacional, a sistemática da troca de informações e as características físicas do modelo de escritório. Os resultados indicaram que a troca de informações intensificada, a forma de distribuição do layout e o modelo de escritório adotado afetam negativamente o desempenho do indivíduo, comprometendo-lhe o rendimento laboral e a saúde; no entanto, um vez modificada, não apenas a disposição dos postos de trabalho, mas também as características do mobiliário, e, ainda, levado a efeito um processo educacional quanto ao controle do tom de voz tais mudanças podem tornar o ambiente, para o desenvolvimento mais eficiente das atividades, sem comprometimento do fluxo de comunicação.

Interação heparina-antitrombina : reconhecimento molecular caracterizado por ferramentas de modelagem molecular

A heparina foi isolada no início do século XX e permanece, até os dias atuais, como um dos mais importantes e eficientes agentes terapêuticos de ação antitrombótica. Sua atividade anticoagulante deve-se à ativação da antitrombina (AT), uma serpina responsável pela inibição fisiológica de serino-proteinases plasmáticas, tais como fIIa e fXa. Os esforços no sentido da elucidação dos aspectos estruturais e dinâmicos associados ao reconhecimento molecular da heparina pela AT, em nível atômico, vêm encontrando diversas dificuldades, principalmente associadas aos compostos sacarídicos. Em decorrência de sua elevada polaridade e flexibilidade, glicosaminoglicanos como a heparina são difíceis de estudar e modelar. Soma-se a isto o fato de que os resíduos de iduronato presentes na heparina (IdoA) apresentam um incomum equilíbrio conformacional entre estados de cadeira (1C4) e bote-torcido (2SO), sendo esta última estrutura postulada como a possível conformação bioativa. Sendo assim, este trabalho apresenta um estudo de modelagem molecular do perfil conformacional da heparina, tanto em solução quanto complexada à AT, utilizando cálculos ab initio e simulações de dinâmica molecular (DM) Em decorrência da ausência de parâmetros capazes de descrever polissacarídeos nos campos de força atualmente disponíveis, cargas atômicas foram geradas, utilizando-se os esquemas de Mulliken, Löwdin e cargas ajustadas ao Potencial Eletrostático, através de cálculos quantum-mecânicos na base 6-31G** e testadas na simulação de DM da heparina. Diversas condições de simulação, tais como modelos de água, concentrações de sais e as conformações 1C4 e 2SO do IdoA, foram avaliadas de forma a identificar as melhores condições para a descrição conformacional da heparina em solução. O protocolo de DM obtido foi então utilizado no estudo do complexo AT-heparina. Os resultados obtidos estão de acordo com os dados experimentais atualmente disponíveis acerca da conformação da heparina e da contribuição energética de cada resíduo de aminoácido da AT para a formação do complexo com a heparina. A partir dos cálculos ab initio realizados, foi proposto um refinamento na estrutura da heparina determinada por RMN, enquanto que as simulações de DM permitiram reinterpretar aspectos da estrutura tridimensional da AT determinada por cristalografia de raios-X. Adicionalmente, os dados obtidos sugerem que não há requerimento conformacional para a interação do IdoA com a AT Globalmente, os dados indicam que simulações de DM podem ser utilizadas para representar adequadamente a conformação da heparina, assim como para caracterizar e quantificar suas interações com a AT. Assim sendo, propomos o uso de simulações de DM do complexo entre a AT e a heparina, ou entre a AT e compostos derivados da heparina, como uma ferramenta útil no processo de desenvolvimento de novos agentes antitrombóticos e anticoagulantes.

Produção de quarkonium : aspectos perturbativos e não-perturbativos da QCD

Neste trabalho de tese investigamos o papel de dinâmica perturbativa e não-perturbativa da Cromodinâmica Quântica, a teoria das interações fortes, em processos de produção de quarks pesados e de quarkonium-estados ligados de um par de quarks pesados. Um aspecto importante na produção de quarks pesados consiste no tratamento de ordens mais altas em QCD perturbativa, que abordamos por meio de elementos de matriz QCD em segunda ordem dominante (NLO) e através de um gerador de eventos Monte Carlo, mais útil fenomenologicamente, onde a produção perturbativa de pares Q Q e obtida utilizando elementos de matriz em ordem dominante e a aproximação de chuveiros partônicos de processos em ordens mais altas. Os processos suaves formando estados ligados de quarkonium são descritos em termos do modelo de evaporação de cor (CEM), ou alternativamente através do modelo de interações suaves de cor (SCI) e do modelo da lei das áreas generalizado (GAL). Neste trabalho, calculamos as distribuições em xF e p? para o J= e 0 em hadroprodução em alvo xo e no colisionador p p do Tevatron. Outros observáveis como a seção de choque total para J= , 0 e charme aberto tamb em são reproduzidos. Além disso, extrapolamos os modelos para descrever a produção de J= e no futuro colisionador LHC, onde as taxas de produção de J= estão at e uma ordem de magnitude acima de outra predição da literatura, o qual pode implicar em J= ser um ru do não negligenciável para estudos de violação da simetria CP no LHC. Além disso, com o objetivo de descrever as taxas de produção relativas entre os vários estados de charmonium, desenvolvemos um modelo para o mapeamento do espectro contínuo de massas do par c c produzido perturbativamente, nas ressonâncias de charmonium, onde introduzimos uma correlação entre a massa invariante do par produzido perturbativamente e a massa física do estado de charmonium. Outra abordagem importante ao estudo dos aspectos perturbativos e não-perturbativos da QCD na produção de quarks pesados e o formalismo de fatorização k?, o qual investigamos em processos de fotoprodução de charme e bottom, com ênfase em resultados de um modelo de saturação. Efeitos de evolução DGLAP tamb em são estudados, considerando a derivada da distribuição de glíuons. Analisamos em detalhe seções de choque totais e distribuições em pT , mostrando as regiões de validade de cada descrição. Através do estudo de vários aspectos perturbativos e não-perturbativos da QCD, este trabalho de tese contribui para um melhor entendimento da conexão entre essas duas dinâmicas.

Dinâmica molecular investigada com pulsos de femto-segundos

Técnicas experimentais de óptica não-linear resolvidas no tempo são capazes, em diversas situações, de fornecer mais informações sobre a dinâmica e estrutura molecular do que técnicas sem nenhuma resolução temporal. Nesta tese investigou-se três sistemas com técnicas espectroscópicas resolvidas no tempo, utilizando pulsos ultra-curtos de um sistema LASER amplificado, baseado em cristal de Ti:Safira, e de amplificadores ópticos paramétricos. Para este estudo foram implementadas e analisadas quatro diferentes técnicas espectroscópicas resolvidas no tempo: Espalhamento Raman anti-Stokes Coerente (CARS), Espalhamento hiper-Rayleigh Resolvido no Tempo (TRHRS), Bombeio-Prova e Bombeio- Depleção-Prova. O três sistemas investigados apresentam um grau crescente de complexidade nas suas dinâmicas e interações com o meio. O primeiro dos sistemas estudados foi a dinâmica da molécula de H2 no regime de impacto, na presença de moléculas e átomos perturbadores (N2 e He). Foi possível determinar, pela primeira vez, coeficientes de alargamento para o ramo Q devido à colisões entre H2+H2 utilizando a técnica de CARS resolvido temporalmente, bem como novos coeficientes para o deslocamento das linhas. O emprego da mesma técnica em sistemas binários lançou mais luz sobre a possibilidade de existir alguma não-homogeneidade em tais tipos de colisões. O segundo experimento investigou a dinâmica de relaxação orientacional da acetonitrila em solução utilizando a nova técnica de TRHRS; esta foi desenvolvida e testada com sucesso pela primeira vez neste trabalho. Um modelamento teórico mostra que dentro do modelo de Debye para a difusão rotacional, tempos associados ao terceiro momento do tempo de difusão (τ3) e ao primeiro momento (τ1) devem ser observados. Os resultados experimentais para a acetonitrila apresentam um decaimento consistente com este modelo (τ3~850 fs), concordando também com valores derivados de resultados experimentais e teóricos obtidos por outros métodos. Além disso, também foi observado experimentalmente uma segunda componente mais rápida (~50 fs), que não pôde ser explicada dentro do modelo de difusão rotacional livre. Esta componente foi atribuída a efeitos reorientacionais coletivos em conjunto com efeitos de superposição temporal de pulsos. A interpretação do sinal desta nova técnica foi testada com sucesso realizando o experimento em tetracloreto de carbono.O terceiro sistema investigado foi a dinâmica de relaxação de energia intramolecular de carotenóides. A principal questão, da vii existência ou não de estados eletrônicos adicionais, foi abordada utilizando uma modificação da técnica de Bombeio e Prova, denominada de Bombeio-Depleção-Prova. A técnica de Bombeio-Depleção-Prova foi pela primeira vez utilizada em carotenóides, possibilitando observar dinâmicas nunca vistas por qualquer outra técnica experimental. Dentro da resolução temporal dos experimentos realizados, a presença de novos estados eletrônicos em carotenóides livres em solução é descartada. Além disso, uma nova ferramenta de análise de dados foi desenvolvida. A análise de alvos espectrais baseada em um algoritmo evolutivo trouxe mais argumentos para o modelo clássico de três níveis. A dinâmica de carotenóides livres em solução foi comparada com sua dinâmica quando inseridos no complexo fotossintético de coleta de luz (LH2) da bactéria púrpura, utilizando para isso a técnica de Bombeio e Prova. Os resultados obtidos mostram um sistema de níveis de energia mais complexo do que quando o carotenóide está livre em solução. Esta modificação devido ao ambiente é possivelmente gerada pela quebra de simetria do carotenóide quando inserido no complexo LH2, possibilitando o cruzamento entre sistemas (singlete-triplete).

Estudo da difusão translacional em um modelo para água

Apesar da água ser o líquido mais comum na natureza, suas características ainda não estão totalmente explicadas. A relação entre a forma do potencial intermolecular efetivo que representa as interações presentes na água, as várias anomalias existentes na água e a possível existência de dupla criticalidade ainda é uma questão em aberto. Recentemente descobriu-se que a água apresenta, além de anomalias termodinâmicas, anomalia na difusão translacional e rotacional. Mostrou-se que para um modelo SPC/E para água estas anomalias dinâmicas estão conectadas com a temperatura de máxima densidade (TMD) e que as anomalias dinâmicas ocupam uma região maior que a TMD no diagrama de fases p vs. T. Nesta dissertação investigamos a relação entre a anomalia na densidade e a anomalia na difusão translacional da água em um modelo onde as moléculas interagem através de um potencial que contém duas fontes de interação: uma isotrópica atrativa do tipo van der Waals, v, e uma direcional, u, que pode ser nula ou repulsiva. Reproduzimos o diagrama de fases pressão vs. temperatura aonde a anomalia na densidade e a dupla criticalidade se fazem presente. Mostramos que este sistema apresenta anomalia na difusão translacional e que esta se encontra na mesma região do diagrama de fases pressão vs. temperatura em que a TMD está presente. Diferentemente do que ocorre para água SPC/E, a região de anomalia na difusão translacional ocupa um espaço menor, encontrando-se na parte interna da TMD. Sugerimos que a discretização do potencial afeta mais fortemente a anomalia dinâmica do que a termodinâmica e que leva a um encolhimento na região de anomalia na difusão.

Teorias de calibre supersimétricas formuladas num espaço-tempo não-comutativo tridimensional

Teorias de calibre formuladas em um espaço-tempo não cumulativo têm sido intensamente estudadas nos últimos anos. O interesse nesse assunto possui motivações provenientes da teoria de cordas. Uma das propriedades mais notáveis das teorias das não-cumulativas consiste de uma estrutura não usual de divergências,a chamada mistura UV/IR, que pode levar ao aparecimento de divergências infravermelhas não integráveis. A eliminação de tais divergências é crucial já que elas podem provocar o colapso da série pertubativa. Modelos não- cumulativos supersimétricos tem um lugar proeminente entre as teorias de campo fisicamente interessantes, uma vez que a supersimetria favorece o cancelamento das divergências perigosas. Eles são os melhores candidatos num programa para definir teorias de campo não-cumulativas consistentes. Neste trabalho investigamos a QED e Yang-Mills não-cumulativos supersimétricos em três dimensões usando o formalismo de supercampos. Para o caso abeliano provamos que a mistura UV/IR não é fonte de divergências infravermelhas não integráveis. Além disso, o modelo resulta ser finito na aproximação de um laço. O mesmo se aplica ao caso abeliano porém apenas na representação fundamental do grupo de calibre.

O gás de elétrons degenerado não-comutativo

A dinâmica quântica de sistemas de partícula única não-relativística envolvendo coordenadas não-comutativas, usualmente chamada mecânica quântica não-comutativa, tem sido objeto de numerosas investigações. Neste trabalho estendemos esses estudos para o caso de sistema de muitas partículas. Usamos como protótipo o modelo do gás de elé- trons degenerado cuja dinâmica é bem conhecida no limite comutativo. Nosso objetivo principal aqui é entender qualitativamente, em vez de quantitativamente, as principais modi cações induzidas pela presença de coordenadas não-comutativas. Primeiro veremos que a não-comutatividade modi ca a energia de correlação de troca enquanto preserva a neutralidade elétrica do modelo. Logo, através da teoria de perturbações independente do tempo juntamente com o mapa de Seiberg-Witten mostramos que o potencial de ioniza ção é modi cado pela não-comutatividade e, também, que o parâmetro não-comutativo atua como uma temperatura de referência. Sendo assim, a não-comutatividade levanta a degenerescência do gás de elétrons a temperatura zero.

Ferramentas para análise dinâmica e estrutural de um reboque de linha leve

Neste trabalho é desenvolvida uma metodologia de projeto para identificar as regiões críticas da estrutura de um reboque de linha leve sendo tracionado em pavimentos do tipo rodovia de baixa qualidade e estrada secundária de muito baixa qualidade. Para tanto, são levantados alguns dados experimentais da estrutura, necessários para a aproximação e simulação dinâmica de um modelo simplificado. A excitação da base é realizada por atuadores que simulam as oscilações verticais de um perfil de estrada, a qual é definida de acordo com os estudos realizados por Dodds e Robson (1973). Isto permite a determinação de um histórico de carregamentos das regiões da estrutura do chassi sob a ação das molas da suspensão. Em seguida, é gerado um modelo estrutural simplificado do reboque em elementos finitos, chamado de global, no qual são determinadas as regiões sob ação das maiores tensões. Tendo identificada a região mais crítica da estrutura, é criado um modelo local desta parte, onde se pode observar a distribuição de tensões com mais detalhe, permitindo a identificação dos pontos de concentração de tensões. Desta forma, com a aplicação do método de análise global-local é possível a obtenção de resultados detalhados quanto aos esforços da estrutura com um menor custo computacional.

Essays on price dynamics

Esta tese tem como objetivo principal aproximar a evidencia empirica existente sobre os agregados macroeconomicos com as novas evidencias empiricas baseadas nos micro dados de precos ao consumidor, tendo como base os modelos padroes de rigidez de preco utilizados na literatura de politica monetaria. Para isso, esta tese utiliza a base de dados individuais de precos ao consumidor no Brasil fornecida pela Fundacao Getulio Vargas. Especificamente, esta tese foca em tres temas principais: a existencia de variac˜oes temporararias de precos, a heterogeneidade na rigidez de precos entre firmas de um mesmo setor e o formato das func˜oes hazard. Os resultados mostram que: existe de fato uma correlac˜ao entre as variaveis referentes as mudancas temporararias de precos e os agregados macroeconomicos; a heterogeneidade na rigidez de precos entre firmas de um mesmo setor apresenta efeitos significativos sobre a dinamica dos agregados macroeconomicos; e por fim, o formato mais geral da func˜ao hazard proposta nesta tese possibilita novas dinamicas dos agregados macroeconomicos.

Perspectivas da educação permanente

O mundo da educação sofre os estremecimentos e a perturbação de toda gestação nova e sua renovação tornou-se necessária pelas disfunções da prática educativa diante das exigências colocadas pela transformação das estruturas sócio-econômicas e pelo progresso científico-tecnológico. A aceleração das mudanças multiplica os problemas, gera crises, lança contínuos desafios ao homem e à sociedade como um todo e determina transformações substanciais na concepção da educação, da escola e de suas funções. A possibilidade antropológica de formação contínua ultrapassa os limites e as barreiras de um sistema educativo institucionalizado para se estender ao longo de toda a vida de um indivíduo. Contudo, as iniciativas de renovação pedagógica não podem contribuir para uma real transformação social, senão na medida em que elas estejam unidas às mudanças de mentalidade e de estruturas.

Influência de parâmetros moleculares em funções de correlação temporal na dinâmica de solvatação mecânica

No presente trabalho descrevemos nossos resultados relativos à investigação da dinâmica de solvatação mecânica por meio de simulações por dinâmica molecular, respeitando o regime da resposta linear, em sistemas-modelo de argônio líquido com um soluto monoatômico ou diatômico dissolvido. Estudamos sistematicamente a influência dos parâmetros moleculares dos solutos (tamanho, polarizabilidade) e da densidade frente a vários modelos de solvatação. Funções de Correlação Temporal da Energia de Solvatação foram calculadas com relação à correlações de n-corpos (n = 2; 3) distinguindo interações repulsivas e atrativas para ambos os sistemas líquidos. Também obtivemos segundas derivadas temporais dessas funções referindo-se à parcelas translacionais, rotacionais e roto-translacionais na solução do diatômico. Encontramos que funções de correlação temporal coletivas podem ser razoavelmente bem aproximadas por correlações binárias a densidades baixas e, a densidades altas, correlações ternárias tornam-se mais importantes produzindo um descorrelacionamento mais rápido das funções coletivas devido a efeitos de cancelamento parciais. As funções de correlação para interações repulsivas e atrativas exibem comportamentos dinâmicos independentes do modelo de solvatação devido a fatores de escalonamento linear que afetam apenas as amplitudes das dessas funções de correlação temporal. Em geral, os sistemas com grau de liberdade rotacional apresentam tempos de correlação mais curtos para a dinâmica coletiva e tempos de correlação mais longos para as funções binárias e ternárias. Finalmente, esse estudo mostra que os sistemas contendo o diatômico relaxam-se predominantemente por mecanismos translacionais binários em modelos de solvatação envolvendo alterações apenas na polarizabilidade do soluto, e por mecanismos rotacionais atrativos binários em modelos envolvendo alterações no comprimento de ligação.

O processo de humanização: as bases constitutivas da subjetividade

Iniciando-se por uma descrição longitudinal, seguida de interpretação "a posteriori", obteve-se categorias relativas a cada dimensão do psiquismo. Estas, em seu caráter de transversalidade, foram tomadas como esteios para as explicações acerca da passagem do sujeito da condição biológica à condição de sujeito falante. O norteamento dessa análise teve como ponto de partida a postulação acerca dessas dimensões, dos aspectos relativos as mesmas, e dos efetores referentes por essa passagem, cujo início tem a esses aspectos que respondem lugar quando a mãe se dispõe a apontar para o sujeito o limite entre o núcleo filogenético e as transformações sócio-histórica. A partir dessa delimitação, ocorre a assimilação de características, tipicamente humanas para, enfim, serem fixados pela lei paterna, os limites do “tornar-se humano”. Nesse ínterim, são observadas ocorrências relativas as dimensões mencionadas. Tratando-se de permanência de conteúdos mantidos fora do campo da simbolização, e excluídos da matriz imaginária, tem-se então a caracterização do real ~elo seu mecanismo específico: a foraclusão. Ainda se situa a pulsão no seu caráter impensável e inominável. Já as ações resultantes do recalcamento originário concorrem para a formação de marcas não-simbolizáveis que, como faceta do imaginário, se vinculam às informações do núcleo filogenético referentes aos protofantasmas, aqui considerado corno outra faceta desta dimensão. No que concerne ao simbólico, registra-se a significação como fundante da condição humana, resultante do recalcamento propriamente dito. Sendo assim, a humanização enquanto explicada pelo conceito de recalcamento, conforme evidencia a metapsicologia freudiana, não açambarça todas as nuances relativas ao humano; visto englobar apenas aquilo que é imaginado corno fantasma ou aquilo que e apresentado simbolicamente. Por isso, a utilização do conceito de foraclusão, da perspectiva lacaniana, justifica-se pelo fato de oferecer uma visão mais completa dessa dinâmica. Desse modo, se pode incluir no escopo do "tornar-se humano", aqueles elementos não passíveis de simbolização, bem como aqueles da realidade pulsional que não se apresentam no texto do fantasma; mas que, ao se constituírem como um "pano-de-fundo", possibilitam que a conteúdos simbólicos e imaginários viabilizem a socialização. De resto, a socialização é um processo que, originado do simbólico, se assenta no imaginário tomado como base. Para isto, o real, em seu efeito marginal, possibilita o retorno ao simbólico. A partir daí, se verifica a utilização dos elementos culturais, que são descobertos em função daquilo que a atividade fantasmática suscita.

Uma nova assinatura para glueballs

A representação de Fock-Tani é um formalismo de teoria de campos para tratar problemas envolvendo simultaneamente partículas compostas e seus constituintes. O formalismo foi originalmente desenvolvido para tratar problemas de física atômica e mais tarde estendido para problemas da física hadrônica. Nesta dissertação, inicialmente apresentamos uma revisão da Cromo dinâmica Quântica e dos modelos de quarks e de glúons constituintes. Revisamos também a representação de Fock-Tani para mésons e buscamos estendê-Ia para estados exóticos, mais precisamente para glueballs. Neste formalismo uma mudança de representação é implementada através de um operador unitário, tal que estados de glueballs no espaço de Fock, compostos por um par glúon-glúon, sejam descritos em termos de operadores de campo de glueballs elementares em um espaço de Hilbert estendido. A aplicação do operador unitário a um Hamiltoniano microscópico de gluons leva a um Hamiltoniano efetivo que descreve todos os possíveis processos envolvendo glúons e glueballs. Esse Hamiltoniano efetivo é utilizado para estudar as interações entre glueballs à baixa energia em um modelo de glúons constituintes, que interagem através da troca de um glúon virtual e são confinados por um potencial fenomenológico O méson J?C = 0++ pode ser descrito de duas formas completamente diferentes. Na primeira forma ele é descrito da maneira usual como sendo composto por um par qq. Na segunda abordagem este mesmo méson é descrito como sendo constituido por dois glúons (glueball). Também estudamos nesta dissertação o méson 2++ no mesmo contexto. Os resultados obtidos nesta dissertação indicam que o méson qq precisaria ter um raio quase igual ao dobro do respectivo raio do glueball para que suas seções de choque de espalhamento elástico fossem equivalentes. As diferenças acentuadas encontradas nas seções de choque de espalhamento elástico méson-méson e glueball-glueball, podem ser interpretadas como uma nova assinatura de glueballs.

Um estudo da otimização da geometria de um pára-quedas simplificado

Neste trabalho desenvolve-se um estudo numérico do fluxo de ar em torno da geometria de um pára-quedas tradicional simplificado, para alguns valores de Reynolds. O método baseia-se na solução das equações incompressíveis de Navier- Stokes discretizadas pelo método de diferenças finitas e integradas pelo método de Runge-Kutta. Utiliza-se o método dos contornos virtuais para representar a geometria numa malha cartesiana e o método de otimização não-linear dos poliedros flexíveis para otimização do coeficiente de arraste calculado através do código de dinâmica de fluidos computacional; esteé um método de busca multivariável, onde o pior vértice de um poliedro com n + 1 vérticesé substituído por um novo.