Resolução numérica de escoamentos compressíveis empregando um método de partículas livre de malhas e o processamento em paralelo (CUDA)

Os métodos numéricos convencionais, baseados em malhas, têm sido amplamente aplicados na resolução de problemas da Dinâmica dos Fluidos Computacional. Entretanto, em problemas de escoamento de fluidos que envolvem superfícies livres, grandes explosões, grandes deformações, descontinuidades, ondas de choque etc., estes métodos podem apresentar algumas dificuldades práticas quando da resolução destes problemas. Como uma alternativa viável, existem os métodos de partículas livre de malhas. Neste trabalho é feita uma introdução ao método Lagrangeano de partículas, livre de malhas, Smoothed Particle Hydrodynamics (SPH) voltado para a simulação numérica de escoamentos de fluidos newtonianos compressíveis e quase-incompressíveis. Dois códigos numéricos foram desenvolvidos, uma versão serial e outra em paralelo, empregando a linguagem de programação C/C++ e a Compute Unified Device Architecture (CUDA), que possibilita o processamento em paralelo empregando os núcleos das Graphics Processing Units (GPUs) das placas de vídeo da NVIDIA Corporation. Os resultados numéricos foram validados e a eficiência computacional avaliada considerandose a resolução dos problemas unidimensionais Shock Tube e Blast Wave e bidimensional da Cavidade (Shear Driven Cavity Problem).

Investigations of nonhistone chromosomal proteins

The major nonhistone chromosomal proteins (NHC proteins) are a group of 14-20 acidic proteins associated with DNA in eukaryotic chromatin. In comparisons by SDS gel electrophoresis (molecular weight sieving) one observes a high degree of homology among the NHC protein fractions of different tissues from a given species. Tissue-specific protein bands are also observed. The appearance of a new NHC protein, A, in the NHC proteins of rat liver stimulated to divide by partial hepatectomy and of rat ascites cells suggests that this protein may play a role in preparing the cell for division. The NHC proteins of the same tissue from different species are also very similar. Quantitative but not qualitative changes in the NHC proteins of rat uterus are observed on stimulation (in vivo) with estrogen. These observations suggest that the major NHC proteins play a general role in chromatin structure and the regulation of genome expression; several may be enzymes of nucleic acid and histone metabolism and/or structural proteins analogous to histones. One such enzyme, a protease which readily and preferentially degrades histones, can be extracted from chromatin with 0.7 N NaCl.

Although the NHC proteins readily aggregate, they can be separated from histone and fractionated by ion exchange chromatography on Sephadex SE C-25 resin in 10 M urea-25% formic acid (pH 2.5). Following further purification, four fractions of NHC protein are obtained; two of these are single purified proteins, and the other two contain 4-6 and 4-7 different proteins. These NHC proteins show a ratio of acidic to basic amino acids from 2.7 to 1.2 and isoelectric points from apparently less than 3.7 to 8.0. These isolated fractions appear more soluble and easier to work with than any whole NHC protein preparation.

Identificação do padrão biomecânico de condutas motoras em escolares: o caso da marcha

A constituição de uma base de dados normativa da marcha é essencial para o diagnóstico e o tratamento de padrões atípicos da locomoção (SUTHERLAND et al., 1997). Não obstante, são escassas as informações relativas aos padrões normais da marcha de crianças (GANLEY e POWERS, 2005), carência ainda mais evidente no tocante à produção acadêmica sobre o padrão biomecânico da locomoção da população de crianças brasileiras. Nesse sentido, especialistas em análise de marcha alertam para o fato que crianças de diferentes populações podem exibir diferentes padrões de marcha de acordo com os grupos étnicos das quais foram extraídas (MORENO-HERNÁNDEZ et al., 2010). Assim sendo, o objetivo do presente estudo foi descrever o padrão biomecânico da marcha de crianças hígidas brasileiras entre 6 e 11 anos de idade. Cento e vinte e duas crianças hígidas, entre seis e 11 anos de idade foram aleatoriamente recrutadas de um universo de 328 alunos. Os sujeitos foram alocados em três grupos etários: Grupo 1 (6-7 anos), Grupo 2 (8-9 anos) e Grupo 3 (10-11 anos). Para o registro das imagens das marchas das crianças foi utilizado um sistema de captura bidimensional de movimento a uma freqüência de aquisição de 30 Hz, composto por uma câmera Sony modelo HC 46 posicionada ortogonalmente a 6 metros da pista. Marcadores esféricos reflexivos de 20mm de diâmetro foram fixados em ambos os lados do corpo dos participantes. Os valores em bruto das coordenadas dos marcadores foram transformadas em coordenadas globais 2D (CALDWELL et al., 2004) e processadas no software SkillSpector (Versão 1.0). A estratégia de Hof (1996) foi utilizada para a normalização dos dados da marcha. Os comprimentos de passo e passada apresentaram uma tendência de aumento com o avanço da idade até 8-9 anos de idade, ao passo que a cadência dos passos apresentou uma tendência de diminuição até o mesmo período. Os números não-dimensionais não apresentaram qualquer tendência de alteração com o avanço da idade. Os três grupos etários apresentaram trajetórias angulares articulares semelhantes. O presente estudo constitui ação pioneira no que tange à descrição do padrão cinemático da marcha de crianças hígidas brasileiras entre 6 e 11 anos de idade. Assim sendo, consideramos que um primeiro passo foi dado no sentido da constituição de uma base de dados normativa da locomoção desses indivíduos.

Spin-related phenomena in magnetic and spin-orbit couppled bose-einstein condensates

129 p.

Habilidades sociais de pais e mães de crianças do Rio de Janeiro: estudo preliminar sobre sua relação com as práticas e estilos parentais

O estilo parental tem influência no desenvolvimento humano e pode ser definido como um padrão de comportamentos identificáveis nos pais quando se relacionam com seus filhos. A literatura tem utilizado alguns parâmetros para compreender os estilos parentais. Um deles se caracteriza pelo binômio: limites-afeto ou responsividade e exigência, delimitando três estilos parentais: Autoritário, Permissivo e Autoritativo. O estilo autoritário caracteriza-se por pais que dão muitos limites e pouco afeto. Estes pais procuram manter a ordem e o controle da família, devendo a obediência ser alcançada ainda que sob padrões punitivos. O estilo permissivo caracteriza-se por pais que dão pouco limite e muito afeto. Nesse caso, há uma crença naturalista segundo a qual a criança deve expressar livremente suas necessidades e aprender por si só. Estes pais não se vêem como responsáveis por ensinar ou modificar comportamentos infantis. O estilo autoritativo caracteriza-se por pais que dão muito limite e muito afeto. O cuidador apresenta condutas intermediárias entre as permissivas e autoritárias, estimulando verbalizações emocionais, permitindo que os filhos participem das decisões, estimulando argumentações, autonomia e disciplina. A empatia é considerada uma habilidade essencial para a manutenção dos laços sociais, que propicia a percepção da necessidade do outro e ao mesmo tempo favorece a sensação de ser compreendido. A assertividade é a habilidade de expressar e defender sua opinião com firmeza e segurança, sem, contudo desrespeitar o outro. Esta pesquisa tem como hipótese a idéia de que pais autoritativos seriam mais empáticos e assertivos do que pais autoritários. O objetivo deste estudo foi correlacionar os Estilos Parentais s Habilidades Sociais: Assertividade e Empatia. Participaram da pesquisa 64 pais, 21 homens e 43 mulheres. Os instrumentos utilizados foram: o Inventário de Habilidades Sociais (IHS-Del Prette), ao Inventario de Empatia (IE FALCONE) e ao Parental Autority Questionary (PAQ). A Análise dos dados foi feita através do SPSS. Foram efetuadas correlações r de Person entre os fatores do PAQ e os fatores do IHS e do IE. Os resultados apontaram para uma correlação positiva (r = 0,30 ; p < 0,05) entre o fator Enfrentamento e Auto-afirmação do IHS e o estilo Autoritativo. Uma relação negativa (r = -0,24 ; p < 0,05) foi observada entre este mesmo fator e o estilo Autoritário. Em relação a empatia, os fatores Tomada de Perspectiva (r = 0,33 ; p < 0,01), Sensibilidade Afetiva (r = 0,25 ; p < 0,05) e Autruismo (r = 0,25 ; p < 0,05) correlacionaram-se positivamente ao estilos Autoritativo, Em relação ao fator Autoritário, foi verificada uma correlação negativa ( r = -0,24 ; p < 0,05) com o fator Altruísmo do IE. A partir deste resultado, foram realizadas três análises de regressão linear múltipla (stepwise) considerando cada um dos fatores do PAQ como variáveis dependentes dos fatores do IHS e do IE. Os Fatores Enfrentamento e Auto-afirmação do IHS e o fator Autruísmo do IE junto explicam 15% da variância do Fator autoritativo. Os resultados apresentados são discutidos a luz da psicologia do desenvolvimento, psicologia social e neuropsicologia.

Charge and spin transport in graphene devices

170 p.

Condições de contorno albedo para cálculos globais de reatores nucleares térmicos com o modelo de ordenadas discretas a dois grupos de energia

Como eventos de fissão induzida por nêutrons não ocorrem nas regiões nãomultiplicativas de reatores nucleares, e.g., moderador, refletor, e meios estruturais, essas regiões não geram potência e a eficiência computacional dos cálculos globais de reatores nucleares pode portanto ser aumentada eliminando os cálculos numéricos explícitos no interior das regiões não-multiplicativas em torno do núcleo ativo. É discutida nesta dissertação a eficiência computacional de condições de contorno aproximadas tipo albedo na formulação de ordenadas discretas (SN) para problemas de autovalor a dois grupos de energia em geometria bidimensional cartesiana. Albedo, palavra de origem latina para alvura, foi originalmente definido como a fração da luz incidente que é refletida difusamente por uma superfície. Esta palavra latina permaneceu como o termo científico usual em astronomia e nesta dissertação este conceito é estendido para reflexão de nêutrons. Este albedo SN nãoconvencional substitui aproximadamente a região refletora em torno do núcleo ativo do reator, pois os termos de fuga transversal são desprezados no interior do refletor. Se o problema, em particular, não possui termos de fuga transversal, i.e., trata-se de um problema unidimensional, então as condições de contorno albedo, como propostas nesta dissertação, são exatas. Por eficiência computacional entende-se analisar a precisão dos resultados numéricos em comparação com o tempo de execução computacional de cada simulação de um dado problema-modelo. Resultados numéricos para dois problemas-modelo com de simetria são considerados para ilustrar esta análise de eficiência.

Mugilid Fish Are Sentinels of Exposure to Endocrine Disrupting Compounds in Coastal and Estuarine Environments

Effects on fish reproduction can result from a variety of toxicity mechanisms first operating at the molecular level. Notably, the presence in the environment of some compounds termed endocrine disrupting chemicals (EDCs) can cause adverse effects on reproduction by interfering with the endocrine system. In some cases, exposure to EDCs leads to the animal feminization and male fish may develop oocytes in testis (intersex condition). Mugilid fish are well suited sentinel organisms to study the effects of reproductive EDCs in the monitoring of estuarine/marine environments. Up-regulation of aromatases and vitellogenins in males and juveniles and the presence of intersex individuals have been described in a wide array of mullet species worldwide. There is a need to develop new molecular markers to identify early feminization responses and intersex condition in fish populations, studying mechanisms that regulate gonad differentiation under exposure to xenoestrogens. Interestingly, an electrophoresis of gonad RNA, shows a strong expression of 5S rRNA in oocytes, indicating the potential of 5S rRNA and its regulating proteins to become useful molecular makers of oocyte presence in testis. Therefore, the use of these oocyte markers to sex and identify intersex mullets could constitute powerful molecular biomarkers to assess xenoestrogenicity in field conditions.

Análise proteômica comparativa dos componentes da superfície celular e do secretoma de Aspergillus fumigatus

Aspergillus fumigatus é o principal agente etiológico da aspergilose invasiva, infecção fúngica oportunista com altas taxas de mortalidade afetando, principalmente, pacientes com neutropenia profunda e prolongada. Durante o processo de invasão e disseminação características desta infecção sistêmica, os conídios do fungo inalados e não eliminados pelas células do sistema imune inato diferenciam-se em hifas que, por sua vez, são angioinvasivas. Pouco se conhece sobre as moléculas da parede celular envolvidas na patogênese do A. fumigatus e/ou secretadas por este patógeno. Neste contexto, este trabalho procura ampliar o entendimento desta doença através do estudo de proteínas diferencialmente expressas na superfície de A. fumigatus durante a morfogênese. Foi utilizada uma abordagem proteômica e foram estudados extratos de superfície de células de A. fumigatus em diferentes estágios durante o processo de filamentação. Estas células foram denominadas, de acordo com o tempo de cultivo e a morfologia, como: TG6h (tubo germinativo), H12h ou H72h (hifas). As proteínas de superfície celular foram extraídas, a partir de células intactas, por tratamento brando com o agente redutor DTT (ditiotreitol). Observou-se que o perfil funcional das proteínas expressas por H12h e H72h foi similar, com exceção de proteínas relacionadas à resposta ao estresse, enquanto o perfil para TG6h apresentou diferenças significativas para vários grupos funcionais de proteínas quando comparado às hifas. Desta forma, foram realizados experimentos de proteômica diferencial entre tubo germinativo (TG6h) e a hifa madura (H72h), pela técnica de DIGE (differential gel electrophoresis). Os resultados revelaram que entre as proteínas diferencialmente expressas, aquelas relacionadas às vias de biossíntese e outras denominadas multifuncionais encontraram-se superexpressas em TG6h. Em relação às proteínas de resposta a estresse, observou-se que algumas HSPs eram mais expressas neste morfotipo, enquanto a MnSOD, relativa à resposta ao estresse oxidativo, era mais abundante na hifa. Com exceção da PhiA, integrante da parede celular, as proteínas identificadas como diferencialmente expressas na superfície do A. fumigatus não possuem sinal para secreção identificável, enquadrando-se nas proteínas atípicas de superfície. Foi verificada a integridade da membrana celular após tratamento com DTT, bem como a marcação por biotina das proteínas extraídas, o que comprovou sua localização superficial na célula fúngica. Hipóteses de que estas proteínas sejam endereçadas à parede celular por via secretória alternativa sustentam estes dados. Estas evidências foram confirmadas pelo fato de não terem sido encontradas as mesmas proteínas da superfície na análise do secretoma do A. fumigatus. Além disso, todas as proteínas caracterizadas no secretoma apresentavam sinal de secreção determinado pelo FunSecKB (www.proteomics.ysu.edu/secretomes/fungi.php). A análise do secretoma foi realizada utilizando-se a cepa selvagem AF293 e a mutante ∆prtT, mutante para um fator de transcrição que atua na regulação da secreção de proteases. Os resultados revelam a ALP1 como expressa na cepa selvagem, assim como outras proteases importantes para virulência e desenvolvimento da célula fúngica, estando suprimidas quando o gene prtT foi deletado.

Simulação de escoamentos incompressíveis empregando o método Smoothed Particle Hydrodynamics utilizando algoritmos iterativos na determinação do campo de pressões

Nesse trabalho, foi desenvolvido um simulador numérico (C/C++) para a resolução de escoamentos de fluidos newtonianos incompressíveis, baseado no método de partículas Lagrangiano, livre de malhas, Smoothed Particle Hydrodynamics (SPH). Tradicionalmente, duas estratégias são utilizadas na determinação do campo de pressões de forma a garantir-se a condição de incompressibilidade do fluido. A primeira delas é a formulação chamada Weak Compressible Smoothed Particle Hydrodynamics (WCSPH), onde uma equação de estado para um fluido quase-incompressível é utilizada na determinação do campo de pressões. A segunda, emprega o Método da Projeção e o campo de pressões é obtido mediante a resolução de uma equação de Poisson. No estudo aqui desenvolvido, propõe-se três métodos iterativos, baseados noMétodo da Projeção, para o cálculo do campo de pressões, Incompressible Smoothed Particle Hydrodynamics (ISPH). A fim de validar os métodos iterativos e o código computacional, foram simulados dois problemas unidimensionais: os escoamentos de Couette entre duas placas planas paralelas infinitas e de Poiseuille em um duto infinito e foram usadas condições de contorno do tipo periódicas e partículas fantasmas. Um problema bidimensional, o escoamento no interior de uma cavidade com a parede superior posta em movimento, também foi considerado. Na resolução deste problema foi utilizado o reposicionamento periódico de partículas e partículas fantasmas.

Foto-assemblage

FOTO-ASSEMBLAGE consiste em nomenclatura sugestionada para definir os trabalhos que tenho produzido a partir da junção de fotografias digitais. As elaborações e fundamentações desses trabalhos representam também o cerne das pesquisas que resultaram na presente dissertação. Em princípio, o termo foto-assemblage haveria de referir-se a questões técnicas ou formais dessa prática. Contudo, ao desenvolver as pesquisas alguns procedimentos acabaram por determinar certas nuances que revelaram aspectos comuns também em seus e conteúdos. Como resultado de construções artísticas juntando fotografias desde 2009, cheguei às composições sintéticas aqui apresentadas, construídas a partir de duas fotografias. Aventei o nome foto-assemblage por observar nas imagens resultantes ressalvas que as distinguiriam de certas convenções atribuídas à ideia de fotografia. Ao mesmo tempo, as referidas imagens proporiam um possível desdobramento ao entendimento de assemblage enquanto técnica artística. Ainda que não seja uma regra, fotografias revelam imagens de momentos. Em sua relação com a compreensão humana de tempo ou espaço, fotografias quase sempre contêm instâncias mínimas. Fotografias, contudo, podem ser também compreendidas como uma contração de um percurso de tempo. Toda imagem fotográfica pode ser assimilada como resultante de determinados acontecimentos anteriores e mesmo tida como elemento gerador de conseqüências futuras. Seguindo esse entendimento, o que proponho com a foto-assemblage é que essa lide com segmentos de tempo ou de espaço contidos numa mesma imagem. Essas fotografias originárias ganhariam uma nova atribuição, sendo retiradas de seu contexto original, serviriam de balizas do percurso de tempo ou espaço suprimido e subjetivado entre elas. Poeticamente, eventos ocorridos entre as fotografias originárias estariam contidos nas imagens produzidas. O termo assemblage foi incorporado às artes a partir de 1953, por Jean Dubuffet, para descrever trabalhos que seriam algo mais do que simples colagem. A ideia de assemblage se baseia no princípio de que todo e qualquer material ou objeto colhido de nosso mundo cotidiano pode ser incorporado a uma obra de arte, criando um novo conjunto, sem que perca seu sentido original. Esse objeto é arrancado de seu uso habitual e inserido num novo contexto, tecendo laços de relação com os demais elementos, construindo narrativas num novo ambiente, o da obra. Na ideia da foto-assemblage, entretanto, é sugerido uso das imagens fotográficas originárias não como objetos que estariam em um mundo cotidiano, mas sim como imagem na concepção do que seria uma entidade mental. Adoto como que uma visão mágica onde as imagens originárias e básicas estariam numa outra dimensão, num plano bidimensional, não manipulável por nós habitantes da tridimensionalidade. Nesse ambiente imaginário ou não, as fotografias são assentadas consolidando a foto-assemblage. Quando a foto-assemblage se concretiza, se corporifica numa mídia, sendo impressa para uma contemplação, ai então, passaria a integrar nosso mundo tridimensional. O resultado poderia ser admitido como um híbrido, uma terceira coisa, a partir de duas que já não se dissociam mais no ensejo de uma compreensão estética. Ao final da dissertação, apresento experiências práticas que resultaram em quatro séries de imagens em foto-assemblage. Cada série enfatiza aspectos peculiares do que denomino paisagem expandida, representando percursos de tempo, espaço ou trajetos entre o mundo concreto e mundos do inconsciente.

A influência da espessura óssea na medição da altura da borda alveolar em imagens de tomografia computadorizada de feixe cônico e espiral

Imagens de tomografia computadorizada (TC) permitem a visualização, sem distorções ou sobreposições, do complexo maxilo-facial, principalmente do osso alveolar. Estudos demonstraram boa reprodutibilidade e precisão da mensuração da altura da borda alveolar, todavia a influência da espessura óssea ainda é pouco descrita. Através da comparação com a mensuração direta, o objetivo deste estudo foi avaliar a precisão, reprodutibilidade e a influência da espessura óssea, na mensuração da altura da borda alveolar em imagens volumétricas e imagens bidimensionais multiplanares em TC de feixe cônico (TCFC) e em TC espiral (TCE). Utilizando 10 mandíbulas secas de humanos, 57 dentes anteriores foram tomografados em equipamentos iCAT (Imaging Science International, Hatfield, PA, EUA) e Brilliance 64 canais (Philips Eletronics, Eindhoven, Holanda), ambos utilizando voxels de 0,25 mm. Através de imagens volumétricas (3D) e imagens bidimensionais (2D) de cortes multiplanares, foi comparada a mensuração da altura da borda alveolar dessas imagens com a mensuração direta nas mandíbulas, feita por vestibular e lingual, por três avaliadores, com o auxílio de um paquímetro, totalizando 114 bordas alveolares medidas. Alta reprodutibilidade intra-avaliador (0,999 a 0,902) e interavaliador (0,998 e 0,868) foi observada através do índice de correlação intraclasse (ICC). Observou-se alta correlação entre a mensuração direta e indireta da altura da borda alveolar em imagens 2D, sendo r=0,923** e 0,916**, e em imagens 3D, com r=0,929** e 0,954*, em TCFC e TCE, respectivamente. Imagens 2D superestimam a altura da borda alveolar em 0,32 e 0,49 mm e imagens 3D em 0,34 e 0,30 mm, em TCFC e TCE respectivamente. Quando o osso alveolar apresenta espessura de no mínimo 0,6 mm a média da diferença entre medidas diretas e indiretas é de 0,16 e 0,28 mm em imagens 2D e de 0,12 e 0,03 mm em imagens 3D para TCFC e TCE respectivamente, sendo que 95% do limite de concordância varia de -0,46 a 0,79 mm e -0,32 a 0,88 mm em imagens 2D, e de -0,64 a 0,67 mm e -0,57 a 0,62 mm em imagens 3D, para TCFC e TCE respectivamente. Quando o osso alveolar é mais fino do que 0,6 mm a TC é imprecisa, pois 95% do limite de concordância variou de -1,74 a 5,42 mm e -1,64 a 5,42 mm em imagens 2D, e de -3,70 a 4,28 mm e -3,49 a 4,25 mm em imagens 3D, para TCFC e TCE respectivamente. Conclui-se que a mensuração da altura da borda alveolar através de imagens tomográficas apresenta alta reprodutibilidade, sendo que quando a borda alveolar apresenta pelo menos 0,6 mm, a precisão da mensuração é alta, todavia quando esta espessura é menor do que 0,6 mm a técnica é imprecisa.

Teoria de Landau-Ginzburg para o estado supercondutor nemático

O objetivo geral deste projeto é propor um modelo bidimensional que descreva o novo estado supercondutor, que apresenta simetria de cristal líquido, chamado de supercondutor nemático. O estudo começa com uma introdução sobre a teoria de Landau-Ginzburg das transições de fase, onde são discutidos conceitos como parâmetro de ordem e as ordens das transições de fase, que são essenciais para o desenvolvimento deste projeto. Em seguida, é feita uma discussão sobre as principais características dos supercondutores como a resistência zero, o efeito Meissner-Ochsenfeld, os tipos de supercondutores, o surgimento de vórtices e uma análise sobre a teoria de Landau-Ginzburg para transição de fase metal-supercondutor. Após isto, é feita uma abordagem sobre os principais tipos de cristais líquidos, com destaque ao cristal líquido nemático, onde é desenvolvida a teoria de Landau-Ginzburg para transição de fase isotrópica-nemática e um estudo sobre o surgimento de disclinações no cristal líquido nemático em duas dimensões. Por fim, é apresentado o modelo proposto para descrever o estado supercondutor nemático, com a construção da teoria de Landau-Ginzburg, o estudo do acoplamento entre as fases e os defeitos topológicos presentes nesse estado.

Importance de certaines caractéristiques biologiques dans la structuration génétique des espèces de poissons: le cas de Ethmalosa fimbriata et Sarotherodon melanotheron

Six populations of Ethmalosa fimbriata and six of Sarotherodon melanotheron have been analysed using enzymatic electrophoresis. The study of gene flow intensity in these two species indicate that: - In Ethmalosa fimbriata, a migratory species with high fecundity and pelagic eggs, there is a high gene flow between populations (3 Nm 83). - In Sarotherodon melanotheron, a sedentary and mouthbrooder species with low fecundity, there is a low gene flow between populations (1 Nm 4).

Modelagem e simulação do escoamento imiscível em meios porosos fractais descritos pela equação de Kozeny-Carman Generalizada

O presente trabalho trata do escoamento bifásico em meios porosos heterogêneos de natureza fractal, onde os fluidos são considerados imiscíveis. Os meios porosos são modelados pela equação de Kozeny-Carman Generalizada (KCG), a qual relaciona a porosidade com a permeabilidade do meio através de uma nova lei de potência. Esta equação proposta por nós é capaz de generalizar diferentes modelos existentes na literatura e, portanto, é de uso mais geral. O simulador numérico desenvolvido aqui emprega métodos de diferenças finitas. A evolução temporal é baseada em um esquema de separação de operadores que segue a estratégia clássica chamada de IMPES. Assim, o campo de pressão é calculado implicitamente, enquanto que a equação da saturação da fase molhante é resolvida explicitamente em cada nível de tempo. O método de otimização denominado de DFSANE é utilizado para resolver a equação da pressão. Enfatizamos que o DFSANE nunca foi usado antes no contexto de simulação de reservatórios. Portanto, o seu uso aqui é sem precedentes. Para minimizar difusões numéricas, a equação da saturação é discretizada por um esquema do tipo "upwind", comumente empregado em simuladores numéricos para a recuperação de petróleo, o qual é resolvido explicitamente pelo método Runge-Kutta de quarta ordem. Os resultados das simulações são bastante satisfatórios. De fato, tais resultados mostram que o modelo KCG é capaz de gerar meios porosos heterogêneos, cujas características permitem a captura de fenômenos físicos que, geralmente, são de difícil acesso para muitos simuladores em diferenças finitas clássicas, como o chamado fenômeno de dedilhamento, que ocorre quando a razão de mobilidade (entre as fases fluidas) assume valores adversos. Em todas as simulações apresentadas aqui, consideramos que o problema imiscível é bidimensional, sendo, portanto, o meio poroso caracterizado por campos de permeabilidade e de porosidade definidos em regiões Euclideanas. No entanto, a teoria abordada neste trabalho não impõe restrições para sua aplicação aos problemas tridimensionais.