Este trabalho descreve um modelo matemático para simulação da laringe humana durante a fonação. O objetivo foi produzir uma técnica computacional de grande escala de processamento para capturar os fenômenos fisiológicos que ocorrem na laringe durante a vocalização e servir de base para estudos mais aprofundados sobre esta importante estrutura do corpo humano. Usando o método dos elementos finitos como base para discretizar as equações dos tecidos musculares da laringe e das equações de Navier-Stokes, e um modelo de descrição da colisão entre as pregas vocais, o sinal glotal foi obtido a partir de diferentes geometrias de laringes com diferentes propriedades viscoelásticas. Os resultados confirmaram a teoria mioelástica-aerodinâmica que descreve a dinâmica da fonação, reproduzindo inclusive fenômenos fisiológicos que os modelos existentes são incapazes de simular. Estudos adicionais foram feitos para verificar a viabilidade do modelo para simular algumas doenças que danificam a laringe. / This work describes a mathematical model to simulate the human larynx during a phonation. The objective was to produce a large-scale computational technique to capture several physiological phenomena that take place on the larynx during the vocalization and to assist further studies about this important structure of the human body. Using the finite element methods as the way to discretize the muscle tissue equations of the larynx and the Navier-Stokes equations and a model to describe the collision between both vocal folds, the glottal signal for different larynx geometries with different viscoelastic properties was obtained. The results confirmed the myoelastic-aerodynamic theory which describes the dynamic of the phonation, also reproducing physiologic phenomena that current models are unable to simulate. Additional studies were conducted to confirm the feasibility of the model to simulate some diseases that affect the larynx.

O presente trabalho consiste em uma extensão do sistema FreeFlow-2D para simular escoamentos de fluidos não newtonianos bidimensionais com superfí cies livres, onde o fluido é descrito pelos modelos de Cross ou o modelo ``power-law\'\'. O método numérico empregado é o método GENSMAC. As equações governantes são aproximadas pelo método de diferenças finitas em uma malha deslocada e partículas marcadoras são utilizadas para a visualização do escoamento e localização da superfície livre. Resultados numéricos são apresentados. Em particular, a presente implementação é validada comparando-se a solução numérica com uma solução analítica / This work presents an extention of the Freeflow-2D system to non-Newtonian free surface flows. The governing equations are solved by the finite difference method on a staggered grid. Marker particles are used to describe the fluid providing the location and the visualization of the free surface. The methodology employed is based on the GENSMAC method. The fluid is modelled by the Cross and power-law models. Numerical examples are presented. The code is validated by making a comparison between analytical and numerical solutions

No presente trabalho estudam-se diferentes técnicas explícitas, em diferenças finitas, no emprego de algoritmos do tipo velocidade-pressão para simulação de escoamentos incompressíveis. O método de resolução da pressão de maneira direta e explícita, introduzido por Bravo e Claeyssen [BRA 97a], é analisado. Faz-se uma aproximação para o erro causado por esta técnica, e verifica-se como isto afeta a equação da continuidade. As simulações são realizadas na cavidade quadrada, comparando-se os diferentes métodos e validando as aproximações realizadas no estudo do método de resolução da pressão. Além disso, simula-se o escoamento em cavidades profundas e rasas, observandose a formação de vórtices e distribuição de energia cinética. Simulações do escoamento na cavidade cúbica também são apresentadas. / In this work different explicit technics in finite defferences in the application of velocity-pressure algorithm to simulate incompressible flows have been studied. The direct and explicit method of pressure resolution, introduced by Bravo anel Claeyssen [BRA 97a] is analyzed. An approximation to the error caused by this method is made, anel how this affects the continuity equation is verified. The simulations are maele in a square cavity, comparing the differents methoels anel valielating the approximations maele in the study of the pressure resolution method. Besieles this, flow in eleep and shallow cavities is simulateel, observing the formation of vortices and kinetic energy distribution. Simulations of the flow in the cubical cavity are also considered.

No presente trabalho são apresentadas análises de amostras de proteínas totais e aminoácidos empregando-se os reagentes p-benzoquinona e tetraamin cobre (II). No caso do reagente p-benzoquinona, as análises foram feitas através de um sistema inédito em fluxo, no qual a mistura reacional passava através de um reator composto de um tubo capilar de aço inoxidável aquecido num banho de glicerina. Os derivados formados passavam através de uma cela de fluxo acoplada a um espectrofotômetro, onde eram feitas as leituras das absorbâncias. Este método se mostrou mais rápido em relação àquele feito em banho-maria. A utilização de um sistema de análise por injeção em fluxo (FIA) com o reagente tetraamin cobre (II), permitiu que as análises fossem feitas de forma rápida, a baixo custo e com a vantagem de não necessitarem de aquecimento, como no caso da p-benzoquinona. Os resultados das análises de amostras de albumina humana obtidos em ambos os métodos acima apresentaram boa concordância quando comparados àqueles obtidos pelo método Kjeldahl para as mesmas amostras. / This work describes two methods for analysis of total proteins and aminoacids using both p-benzoquinone (PBQ) and tetraamin copper (II) reagents. With the p-benzoquinone reagent the method was performed by an original flow system made with stainless steel capillary tubing, conveniently heated into a glycerin bath, where the reactional mixture travels in its way to the detector. The derivate products are carried out to a flow cell adjusted to a spectrophotometer where the absorbances are measured. This method was efficient, with a fair cost, and providing results very faster than in the batch mode of operation. The second method, using the tetraamin copper (II) reagent, was performed by a common flow injection system (FIA). When compared with the PBQ method, the FIA tetramin copper (II) method was also rapid, efficient, with a fair cost and with the great advantage of simplicity, once it is performed at room temperature. The results obtained with human albumine samples using both methods are in agreement with those ones obtained by the Kjeldahl method for the same samples.

Gas-solid fluidized bed reactors are widely applied in chemical and energy industries, and their design and scale-up are virtually empirical, extremely expensive and time consuming. This scenario has motivated the development of alternative theoretical tools, and two-fluid modeling, where gas and particulate are both treated as interpenetrating continuum phases, has appeared as a most promising approach. Owing to the large domains to be resolved in real-scale fluidized bed reactors, only filtered modeling approaches are feasible, and closure models become necessary to recover sub-grid effects that are filtered by the very coarse numerical grids that are imposed owing to computational limitations. Those closure models, which in hydrodynamic formulations account for filtered interphase momentum exchanges and filtered and residual stresses in the phases, can be derived from results of highly resolved simulations (HRS) performed over small size domains under refined numerical grids. One widely practiced approach consists of applying two-fluid modeling under micro-scale defined closures, generally known as microscopic two-fluid modeling. This approach includes microscopic closures for solid phase stresses derived from the kinetic theory of granular flows (KTGF), which accounts for kinetic-collisional effects only, and is adequate to dilute flows. Otherwise, the conventional KTGF does not account for interparticle friction effects, and its application to dense flow conditions is quite questionable. In this work a literature available modified version of KTGF is applied which also accounts for interparticle friction, and highly resolved simulations are performed for dense flow conditions in order to evaluate the effects of friction over relevant filtered parameters (namely effective drag coefficient, filtered and residual stresses). Ranges of domain average solid volume fractions and gas Reynolds numbers are considered (macro-scale conditions) embracing dense gas-solid fluidized flows from suspensions up to pneumatic transport. The MFIX open source code is used in all the simulations, which are performed over 2D periodical domains for a unique monodisperse particulate. The HRS results (i.e. meso-scale flow fields) are filtered over regions compatible with grid sizes in large scale simulations, and the relevant filtered parameters of concern are derived and classified by ranges of other filtered parameters taken as independent variables (filtered solid volume fraction, filtered slip velocity, and filtered kinetic energy of solid velocity fluctuations, which are referred to as markers). Results show that the relevant filtered parameters of concern are well correlated to all of those filtered markers, and also to all of the imposed macro-scale conditions. Otherwise, interparticle friction showed no significant effects over any filtered parameter. It is recognized that this issue clearly requires further investigation notably regarding the suitability of the markers that were assumed for classifying the filtered results. The current work is intended as a contribution for future developments of more accurate closure models for large scale simulations of gas-solid fluidized flows. / Reatores de leito fluidizado de escoamento gás-sólido são largamente utilizados nas indústrias química e de energia, e o seu projeto e escalonamento são virtualmente empíricos, extremamente caros e demorados. Este cenário tem motivado o desenvolvimento de ferramentas teóricas alternativas, e a modelagem de dois fluidos, onde gás e particulado são ambos tratados com fases contínuas interpenetrantes, tem surgido como uma aproximação muito promissora. Devido aos grandes domínios a serem resolvidos em reatores de leito fluidizado de escala real, apenas aproximações de modelagem filtradas são viáveis, e modelos de fechamento tornam-se necessários para recuperar efeitos sub-malha que são filtrados pelas malhas numéricas grosseiras que são impostas devido as limitações computacionais. Estes modelos de fechamento, que em formulações hidrodinâmicas respondem principalmente por trocas de momentum filtradas entre fases e tensões filtradas e residuais nas fases, podem ser obtidos de resultados de simulações altamente resolvidas (SAR) realizadas em domínios de dimensões reduzidas sob malhas numéricas refinadas. Uma aproximação largamente praticada consiste na aplicação de modelagem de dois fluidos sob fechamentos definidos na micro-escala, genericamente conhecida como modelagem microscópica de dois fluidos. Esta aproximação inclui fechamentos microscópicos para tensões da fase sólida obtidos da teoria cinética dos escoamentos granulares (TCEG), que considera apenas efeitos cinéticos-colisionais, e é adequada para escoamentos diluídos. Por outro lado, a TCEG convencional não leva em conta efeitos de fricção interpartículas, e sua aplicação para condições densas de escoamento é bastante questionável. Neste trabalho aplica-se uma versão modificada da TCEG disponível na literatura que também leva em conta fricção interpartículas, e simulações altamente resolvidas são realizadas para condições de escoamentos densos visando avaliar os efeitos da fricção sobre os parâmetros filtrados relevantes (coeficiente de arrasto efetivo, tensões filtradas e residuais). Considera-se faixas de frações volumétricas de sólido e números de Reynolds do gás médios no domínio (condições de macro-escala) abrangendo escoamentos gás-sólido fluidizados densos desde suspensões até transporte pneumático. O código aberto MFIX é utilizado em todas as simulações, que foram executadas sobre domínios periódicos 2D para um único particulado monodisperso. Os resultados das SAR (i.e., campos de escoamento de meso-escala) foram filtrados sobre regiões compatíveis com tamanhos de malha praticados em simulações de grandes escalas, e os parâmetros filtrados relevantes de interesse são calculados e classificados por faixas de outros parâmetros filtrados tomados como variáveis independentes (fração volumétrica de sólido filtrada, velocidade de deslizamento filtrada, e energia cinética das flutuações de velocidade da fase sólida filtrada, que são referidos como marcadores). Os resultados mostram que os parâmetros filtrados relevantes de interesse são bem correlacionados com todos os marcadores, e também com todas as condições de macro-escala impostas. Por outro lado, a fricção interpartículas não mostrou efeitos significativos sobre qualquer parâmetro filtrado. Reconhece-se que este aspecto claramente requer investigações adicionais, notadamente com respeito à adequação dos marcadores que foram considerados para classificação dos resultados filtrados. O trabalho corrente é posto como uma contribuição para o desenvolvimento futuro de modelos de fechamento mais acurados para simulações de grandes escalas de escoamentos gás-sólido fluidizados.

A presente pesquisa tem como objetivo realizar um estudo teórico e desenvolver simulações computacionais envolvendo a dinâmica de sistemas gás-sólido. O foco principal do trabalho é a modelagem do particulado através da análise das forças de contato entre partículas de materiais granulares utilizando modelos contínuos baseados na mecânica dos solos e na teoria cinética dos escoamentos granulares (sistemas grandes com muitas partículas, formulação Euleriana - Volumes Finitos) e modelos discretos baseados nas características físicas dos materiais (sistemas intermediários e número limitado de partículas, formulação Lagrangeana - Método dos Elementos Discretos). Investigam-se os modelos existentes na literatura com intuito de melhorar os modelos contínuos e discretos baseados na interação entre as partículas que caracterizam a dinâmica do particulado em sistemas gás-sólido. Propõe-se uma nova abordagem para a determinação do coeficiente de rigidez da mola baseada em uma equivalência entre os modelos lineares e não-lineares. Utiliza-se o código fonte MFIX para realizar simulações computacionais da dinâmica de sistemas gás-sólido, analisando o processo de fluidização, mistura e segregação de partículas, influência das correlações de arrasto, e análise das forças de contato entre as partículas através do novo método para a determinação do coeficiente de rigidez da mola . Os resultados obtidos são comparados com dados numéricos e experimentais da literatura. / The purpose of the present study is to perform a theoretical study and develop numerical simulations involving dynamic in gas-solid systems. The focus of the work is the modeling of particulate matter using continuous models based on soil mechanics and the kinetic theory of granular flows (large systems with many particles, Eulerian formulation - Finite Volume) and discrete models based on physical characteristics of the particles (intermediate systems and limited number of particles, Lagrangian formulation - Discrete Element Method). It is proposed a new approach to determine the normal spring stiffness coefficient of the linear model through the numerical solution for the overlap between particles in non-linear models. The linear spring stiffness is determined using an equivalence between the linear and the non-linear models. It is used the MFIX computational code to perform numerical simulations of the dynamics of gas-solid systems. It is analyzed the processes of fluidization, mixing and particle segregation and the influence of drag correlations. The proposed approach for normal spring stiffness coefficient is applied in the numerical simulations of two problems: single freely falling particle and bubbling fluidized bed. The results were compared with numerical and experimental data from literature.

A neutralizacao de pH e utilizada nas industrias, para garantir o descarte seguro de euentes. As plantas de neutralizacao de pH sao um problema complexo de controle, visto que a planta segue um modelo nao-linear e apresenta caractersticas variantes no tempo, o que demanda sua correta modelagem para o projeto de sistemas de controle ecientes. No entanto, a teoria referente a modelagem de pH nao e facilmente aplicada na pratica, resultando frequentemente em modelos que nao predizem corretamente o comportamento dinamico da planta. O primeiro objetivo deste trabalho foi modelar matematicamente a Planta Piloto de Neutralizacao de pH do Laboratorio de Controle de Processos Industriais (LCPI), utilizando uma metodologia que possa ser aplicada para obter o modelo matematico de outras plantas de neutralizacao de pH. Inicialmente a Planta Piloto de Neutralizacao de pH do LCPI foi modelada de acordo com a abordagem fenomenologica, utilizando-se os princpios de conservacao de massa, da eletroneutralidade e os conceitos de equilbrio qumico. Em seguida, o modelo foi ajustado aos dados experimentais do processo (abordagem emprica), utilizando-se curvas de titulacao dos inuentes e distribuicoes de tempos de residencia do reator. Atraves de experimentos, vericou-se que o modelo representou, de forma satisfatoria, a dinamica real da Planta Piloto de Neutralizacao de pH do LCPI. Ademais, este modelo foi utilizado para alcancar o segundo objetivo deste trabalho: projetar um sistema de controle de pH, o qual foi composto por um observador nao-linear e um controlador baseado em modelo. Esta estrutura de controle foi testada experimentalmente, onde certicou-se que os requisitos de controle foram satisfeitos. / The pH neutralization is used in industry to discard properly the wastewater, ensuring the environment preservation. The pH neutralization is a complex control problem, as the model of the plant presents a strong nonlinearity and time varying characteristics, which demands a proper modeling in order to design ecient control systems. However, the application of the theory related to pH modeling is not a trivial task and may result in models that can not predict the plant dynamics. The rst objective of this research was to model the pH Neutralization Pilot Plant, of the Laboratory of Industrial Processes Control (LCPI), using a methodology that could be replicated to model other pH neutralization plants. Initially, the pH Neutralization Pilot Plant was modeled with the phenomenological approach, utilizing rst principles, such as the mass conservation, electroneutrality and chemical equilibrium. Moreover, the model was adjusted to represent the process observed data (empirical approach), as its titration curves of the inuent streams and its reactor residence time distribution. Through experiments, it was veried that the model could represent adequately the real process dynamics. Furthermore, this model was used to achieve the second objective of this research: to design a pH control system, which was composed of a nonlinear observer and a modelbased control. This control structure was tested experimentally, ensuring that the control requirements were satised.

A hemólise induzida mecanicamente em sistemas de bombeamento centrífugo é um processo complexo para o qual ainda se almeja um modelo matemático que possa auxiliar eficazmente o desenvolvimento e otimização destes dispositivos médicos. Este trabalho propõe um equacionamento matemático do escoamento do fluido sanguíneo no interior de bombas centrífugas, levando-se em consideração as variações do campo de tensões através de uma abordagem reológica não-newtoniana. A destruição da fase globular é modelada matematicamente por uma distribuição estatística. Tal modelo bifásico simplificado é utilizado em simulações numéricas de experimentos virtuais propostos e padronizados pela comunidade científica internacional especializada. A validação do modelo é realizada através de comparação com resultados referenciais. Os resultados obtidos demonstram um nova possibilidade de solução numérico-computacional para o problema da hemólise por trauma mecânico. / Mechanically induced hemolysis in centrifugal pumping systems is a complex process for which a mathematical model is still desired. It can effectively aid the development and optimization of these medical devices. This work proposes a mathematical modeling for the blood flow into centrifugal pumps, taking into account the variations of the stress field through a non-Newtonian rheological approach. The destruction of the globular phase is mathematically modeled by a statistical distribution. Such a simplified two-phase model is used in numerical simulations of virtual experiments proposed and standardized by the international scientific community. The validation of the model is performed through comparison with reference results. The results obtained show a new mathematical solution for the problem of hemolysis induced by mechanical trauma.

Orientador: Jun Takahashi / Tese (doutorado) - Universidade Estadual de Campinas, Instituto de Fisica Gleb Wataghin / Made available in DSpace on 2018-08-23T16:49:52Z (GMT). No. of bitstreams: 1 Souza_RafaelDerradide_D.pdf: 5902855 bytes, checksum: ef026256d77f665d1d9d346c08ad7ec8 (MD5) Previous issue date: 2013 / Resumo: Nas colisões nucleares relativísticas, devido ao comportamento coletivo do sistema, é esperado que a anisotropia espacial inicial da região de uperposição entre os núcleos incidentes seja convertida em uma anisotropia final na distribuição de momento das partículas produzidas. Em particular, o coeficiente de fluxo elíptico, dado pelo segundo harmônico (v2) da decomposição de Fourier da distribuição azimutal das partículas finais produzidas, é um observável sensível à anisotropia geométrica inicial. Neste trabalho, estudamos o comportamento do coeficiente v2 obtido a partir de eventos gerados com o código de evolução hidrodinâmica NeXSPheRIO simulando colisões Au+Au a (SNN)1/2 = 200 GeV. O NeXSPheRIO permite o controle da condição inicial, da evolução hidrodinâmica e gera as partículas finais como em um evento real. Portanto, foi possível aplicar as mesmas técnicas de análise de dados utilizadas experimentalmente. Estudamos os efeitos de flutuações na condição inicial no observável v2 e testamos alguns observáveis experimentais propostos que seriam sensíveis a estas flutuações. Como os observáveis também são sensíveis a outros tipos de correlações, geralmente chamadas de non-flow, e considerando que a maior contribuição pode vir de jatos de partículas gerados nessas colisões, estudamos os efeitos de jatos nesses observáveis. Construímos um modelo simplificado de duas componentes a partir da superposição de eventos NeXSPheRIO, com eventos de colisões p+p produzidos pelo gerador de eventos Pythia. Resultados deste modelo descrevem melhor a dependência de v2 com o momento transversal quando comparados com dados experimentais. Determinamos também que, ao contrário do que se esperava, os observáveis sensíveis à flutuação na condição inicial são pouco afetados pela inclusão dos jatos. Este resultado sugere que os efeitos de non-flow, como considerados aqui, são desprezíveis nas medidas de flutuações de v2, indicando que os resultados obtidos pelos experimentos sejam dominados pelas flutuações da condição inicial / Abstract: Due to the collective behavior of the system created in relativistic nuclear collisions, it is expected that the initial spatial anisotropy of the overlap region between the two incoming nuclei converts into an anisotropy of the momentum distribution of the outgoing particles. In particular, the elliptic flow coefficient, which is given by the second harmonic (v2) of the Fourier decomposition of the final particle azimuthal distribution, is an observable sensitive to the initial geometric anisotropy. In this work, we studied the behavior of the elliptic flow coefficient obtained from events generated using the hydrodynamic-based code NeXSPheRIO for Au+Au collisions at (SNN)1/2= 200 GeV. The NeXSPheRIO code allows the control of the initial condition, the hydrodynamic evolution and produces the outgoing particles just like in a real event. Therefore, it was possible to apply the same data analysis techniques used experimentally. We studied the effects of initial condition fluctuations on flow observables and tested some experimental observables thought to be sensitive to the fluctuations. Since these observables are also sensitive to other types of correlations, usually referred to as non-flow, and considering that the major contribution to them may come from jets of particles generated in the collision, we also studied the effects of jets on these observables. A simplified two-component model was built by superimposing events produced with the NeXSPheRIO code and p+p collision events produced with the Pythia event generator. Results from our model present a better description of the ?2 dependence with transverse momentum when compared to experimental data. Also, opposite to what was expected, the observables sensitive to the fluctuations are not much sensitive to the inclusion of jets. This result may suggest that the contribution of non-flow, as considered here, is negligible in the observed flow fluctuations measurement, which indicates that the values observed in the experiments are mainly dominated by fluctuations from the initial condition / Doutorado / Física / Doutor em Ciências

Orientador: Denis José Schiozer / Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Mecânica / Made available in DSpace on 2018-08-20T08:15:37Z (GMT). No. of bitstreams: 1 Paiva_HernaniPetroni_M.pdf: 4091307 bytes, checksum: 9c0f72cc029a22e1f5a96e7587ad9ab5 (MD5) Previous issue date: 2012 / Resumo: A recuperação de petróleo em reservatórios naturalmente fraturados apresenta-se como um risco de projeto, sobretudo em sistemas molháveis ao óleo ou com molhabilidade intermediária, já que a simulação utilizando deslocamento por injeção de água indica significativa redução da recuperação. As fraturas representam descontinuidades do meio poroso e possuem efeito capilar e condutividade hidráulica distintos, o que altera sensivelmente o comportamento do escoamento no reservatório, e também os mecanismos físicos envolvidos no processo de recuperação. A simulação de reservatórios fraturados é geralmente realizada com o modelo de dupla porosidade, que está implementado nos principais simuladores comerciais. Neste modelo os processos físicos envolvidos na recuperação são representados pela função de transferência entre matriz e fratura. No entanto, os simuladores comerciais utilizam diferentes funções de transferências com diferentes modelos para representar o processo de recuperação. Neste trabalho, foi construído um simulador de dupla porosidade no qual foram implementadas as funções de transferência de Kazemi et al. (1976), Sabathier et al. (1998) e Lu et al. (2008) para comparação dos resultados de recuperação utilizando deslocamento por injeção de água em sistemas molháveis ao óleo ou com molhabilidade intermediária. A comparação entre as funções de transferência foi realizada para diferentes combinações de processos físicos, mostrando que há significativo aumento de recuperação em sistemas molháveis ao óleo ou com molhabilidade intermediária, especialmente em reservatórios totalmente descontínuos quando o deslocamento ocorre por embebição concorrente, coerentemente com o resultado experimental de Firoozabadi (2000). As funções de transferência implementadas, associadas ao modelo de dupla porosidade, foram também comparadas a simuladores comerciais e a um modelo de fraturas discretas refinado, obtendo-se, entretanto, resultados distintos, mostrando que os diferentes modelos de função de transferência fornecem diferentes resultados. A injeção de água mostrou-se um método de recuperação efetivo mesmo em reservatórios totalmente descontínuos em sistemas molháveis ao óleo ou com molhabilidade intermediária quando há deslocamento por embebição concorrente com gradientes de pressão nos blocos de matriz suficientemente elevados. Os resultados são sensíveis aos parâmetros de caracterização e variam de acordo com o processo físico utilizado. Portanto, a caracterização de reservatórios naturalmente fraturados deve ser realizada levando-se em consideração os fenômenos físicos e os modelos utilizados para representá-los / Abstract: The naturally fractured reservoir recovery is a project risk specially in oil-wet or intermediate-wet systems because of the simulations results under waterflood displacement. Fractures are porous medium discontinuities with distinct capillarity and hydraulic conductivity properties that change the reservoir flow behaviour as well the physical mechanisms acting in petroleum recovery. Double-porosity models are generally used in fractured reservoir simulation and have been implemented in the major commercial reservoir simulators. The physical processes acting in petroleum recovery are represented in double-porosity models by matrix-fracture transfer functions. Commercial simulators have their own transfer function implementations, and as a result different kinetics and final recoveries are attained. In this work, a double porosity simulator was built with Kazemi et al. (1976), Sabathier et al. (1998) and Lu et al. (2008) transfer function implementations and their recovery results compared using waterflood displacement in oil-wet or intermediate-wet systems. The results of transfer function comparisons show recovery improvements in oil-wet or intermediate-wet systems under different physical processes combination, particularly in fully discontinuous porous medium when concurrent imbibition takes place, coherent with Firoozabadi (2000) experimental results. Furthermore, the implemented transfer functions, related to a double-porosity model, were compared to double-porosity commercial simulators models, as well a discrete fracture model with refined grid, showing differences between them. Waterflood can be an effective recovery method even in fully discontinuous media for oil-wet or intermediate-wet systems where concurrent imbibition takes place with high enough pressure gradients across the matrix blocks. These results are sensitive to reservoir characterization parameters whose sensitiveness depends on the physical process employed. Naturally fractured reservoir characterization must consider the physical phenomena occurring during recovery and the models used to represent them / Mestrado / Reservatórios e Gestão / Mestre em Ciências e Engenharia de Petróleo