Multibody System Dynamics has been responsible for revolutionizing Mechanical Engineering Design by using mathematical models to simulate and optimize the dynamic behavior of a wide range of mechanical systems. These mathematical models not only can provide valuable informations about a system that could otherwise be obtained only by experiments with prototypes, but also have been responsible for the development of many model-based control systems. This work represents a contribution for dynamic modeling of multibody mechanical systems by developing a novel recursive modular methodology that unifies the main contributions of several Classical Mechanics formalisms. The reason for proposing such a methodology is to motivate the implementation of computational routines for modeling complex multibody mechanical systems without being dependent on closed source software and, consequently, to contribute for the teaching of Multibody System Dynamics in undergraduate and graduate levels. All the theoretical developments are based on and motivated by a critical literature review, leading to a general matrix form of the dynamic equations of motion of a multibody mechanical system (that can be expressed in terms of any set of variables adopted for the description of motions performed by the system, even if such a set includes redundant variables) and to a general recursive methodology for obtaining mathematical models of complex systems given a set of equations describing the dynamics of each of its uncoupled subsystems and another set describing the constraints among these subsystems in the assembled system. This work also includes some discussions on the description of motion (using any possible set of motion variables and admitting any kind of constraint that can be expressed by an invariant), and on the conditions for solving forward and inverse dynamics problems given a mathematical model of a multibody system. Finally, some examples of computational packages based on the novel methodology, along with some case studies, are presented, highlighting the contributions that can be achieved by using the proposed methodology. / A Dinâmica de Sistemas Multicorpos tem sido responsável por revolucionar projetos de Engenharia Mecânica pela utilização de modelos matemáticos para simulação e otimização do comportamento dinâmico de uma ampla gama de sistemas mecânicos. Estes modelos matemáticos não somente podem fornecer valiosas informações acerca de um sistema que caso contrário poderiam ser obtidas somente através de experimentos com protótipos, como também têm sido responsável pelo desenvolvimento de diversos sistemas de controle baseados em modelo. Este trabalho representa uma contribuição para a modelagem dinâmica de sistemas mecânicos multicorpos por meio do desenvolvimento de uma nova metodologia modular e recursiva que unifica as principais contribuições de diversos formalismos da Mecânica Clássica. A razão para propor tal metodologia é motivar a implementação de rotinas computacionais para a modelagem de sistemas mecânicos multicorpos complexos sem depender de pacotes de software de código fechado e, consequentemente, contribuir para o ensino de Dinâmica de Sistemas Multicorpos nos níveis de graduação e pós-graduação. Todos os desenvolvimentos teóricos são baseados em e motivados por uma revisão crítica da literatura, conduzindo a uma forma matricial geral das equações dinâmicas de movimento de um sistema mecânico multicorpos (que podem ser expressas em termos de qualquer conjunto de variáveis adotado para a descrição dos movimentos realizados pelo sistema, ainda que tal conjunto inclua variáveis redundantes) e a uma metodologia recursiva geral para a obtenção de modelos matemáticos de sistemas complexos, dado um conjunto de equações descrevendo a dinâmica de cada um de seus subsistemas desacoplados e outro descrevendo os vínculos entre estes subsistemas (no sistema) quando acoplado. Este trabalho também inclui algumas discussões acerca da descrição de movimentos (utilizando qualquer conjunto admissível de variáveis de movimento e admitindo qualquer tipo de vínculo que seja passível de descrição por invariantes), e das condições para a solução dos problemas de dinâmica direta e inversa dado um modelo matemático de um sistema multicorpos. Finalmente, alguns exemplos de pacotes computationais baseados na nova metodologia, juntamente com alguns estudos de caso, são apresentados, ressaltando as contribuições que podem ser alcançadas por meio do uso da metodologia proposta.

Este trabalho teve como objetivo adicionar ao programa computacional SAAM for Windows a capacidade de simular o transporte de poluentes em cursos d´água e substituiu com vantagens o modelo computacional QUAL2E para simulações em regime transiente para substâncias conservativas e não conservativas em casos de injeção instantânea. Desta maneira, os pesquisadores dispõem agora de uma ferramenta adicional para simulação de transporte de poluentes em cursos d´água. O SAAM foi desenvolvido originalmente para a simulação de conjuntos discretos, mas neste trabalho, o programa foi adaptado para a simulação de sistemas contínuos. O trabalho de implementação foi feito tomando-se por base a equação unidimensional de transporte de massa por advecção - dispersão adaptada ao processo de simulação do programa SAAM utilizando-se para isto o método de resolução de MacCormack. O SAAM foi implementado também com o recurso de avaliar as concentrações de DBO e OD de maneira simplificada. Foram analisados diversos cenários de transporte de massa através de simulação computacional utilizando o SAAM, o programa QUAL2E e quando possível, soluções analíticas. Os resultados obtidos nestas simulações foram bastante similares, com isso foi possível concluir que a implementação do programa SAAM foi um sucesso. / This work had as objective to implement the capacity of simulation of pollutant transport in watercourses to the program SAAM for Windows for time-variable simulations for the dispersion of conservative and non-conservative substances from instantaneous sources, as an advantageous alternative to QUAL2E. This way, researchers may have an additional tool for the simulation of pollutant transport in watercourses. SAAM was developed originally for the simulation of discrete systems; the program is now adapted for the simulation of continuous problems. This implementation was made through the use of the advection-dispersion mass transfer equation adapted to SAAM\'s simulation features, using, for numerical solution, the method of MacCormack. SAAM was also implemented for the evaluation of concentrations of BOD and DO in a simplified way. Several scenarios of mass transport were analyzed through computational simulation using SAAM, the program QUAL2E, and, when possible, analytical solutions. As the results obtained in these simulations were quite similar, it is possible to conclude that the proposed implementation to the program SAAM was successful.

Modelos matemáticos oferecem grandes benefícios para a compreensão dos mecanismos envolvidos nos processos de tratamento de águas residuárias uma vez que fornecem interpretações e possibilitam previsões de desempenho, comparações de alternativas de tratamento, otimização de futuras plantas ou o aprimoramento das existentes, podendo subsidiar a elaboração de projetos em escala real. Em virtude disto, nesta pesquisa visou-se o desenvolvimento de um modelo bioquímico-matemático para descrever o processo de degradação da glicose em um reator anaeróbio de leito fixo com fluxo ascendente, com a resultante produção biológica de hidrogênio por meio do processo de fermentação. O desenvolvimento do modelo foi baseado em estudos sobre a cinética bioquímica e as características hidrodinâmicas do sistema. Os parâmetros de ajuste do modelo aos dados experimentais foram as constantes de velocidade das reações bioquímicas envolvidas na produção de hidrogênio. A calibração foi realizada manualmente buscando minimizar o desvio global. Para a determinação dos parâmetros foi utilizada a técnica de geração de números aleatórios com distribuição de freqüência uniforme e em seguida, o método de inversão de matrizes. O modelo matemático se revelou bastante adequado para a previsão do perfil de concentrações ao longo do reator, e possibilitou a representação das rotas de utilização da matéria orgânica. A reação de oxidação do ácido propiônico pelas bactérias acidogênicas produtoras de hidrogênio constitui a principal via de produção de \'H IND.2\' no sistema. / Mathematical models bring benefits to the understanding of mechanisms involved on wastewater treatment processes because they provide interpretations and make possible performance predictions, evaluation of design alternatives, optimization of future plants or the improvement to existing systems. Therefore, in this work a mathematical model to describe the glucose degradation process, with hydrogen production through the fermentation, in an upflow anaerobic packed-bed reactor is developed. The model equations were based on studies of biochemical kinetics and hydrodynamics features of the system. The parameters considered were the rates of the biochemical reactions involved in the hydrogen production. The calibration was made through the minimization of the global deviation. The parameters determination was obtained with the use of a technique of generation of aleatory numbers, and after that, the method of matrices inversion for the solution of the system of linear equations. The mathematical model developed showed to be adequate for the concentrations prediction along the reactor, and it made possible the representation of the routes of organic material utilization. The oxidation reaction of propionic acid is the main hydrogen production route in the reactor.

A presente pesquisa aborda um problema de roteirização e programação de veículos incorporando uma nova restrição operacional: a requisição simultânea de múltiplos veículos para atendimento da demanda. Trata-se de uma característica encontrada em operações de apoio à exploração de petróleo \"offshore\", em que mais de uma embarcação é requerida para executar tarefas de reboque e lançamento de linhas de ancoragem. Esta imposição, somada às restrições de janela de tempo, precedência entre tarefas, autonomia das embarcações e atendimento integral da demanda, configuram este problema. A programação é orientada pela minimização dos custos variáveis da operação e dos custos associados ao nível de serviço no atendimento. Este problema é uma variação do problema clássico de roteirização e programação de veículos com janela de tempo, de classe NP-Difícil. Nesta pesquisa, propõe-se modelar e resolver o problema em escala real por meio do algoritmo \"branch and cut\" acoplado às heurísticas de busca em vizinhança \"local branching\" e \"variable neighborhood search\". Para gerar as soluções iniciais será empregado o método \"feasibility pump\" e uma heurística construtiva. / This research focuses a fleet scheduling problem with new operational constraints: each task requiring multiple types of vehicles simultaneously. This kind of operation occurs in offshore exploitation and production sites, when more than one vessel is needed to accomplish the tugging and mooring of oil platforms. Other constraints are maintained such as time windows, precedence between tasks, route duration and the demand attendance. The solution schedules are cost oriented, which encompasses the routing variable costs and the customer service costs. This is a variation of the classical fleet routing and scheduling, which is an NP-Hard problem. This research aims to solve the real scale problem through a combined use of branch and cut strategy with local search algorithms such as local branching and variable neighborhood search. An efficient heuristic rule will be used in order to generate initial solutions using the feasibility pump method.

O processamento térmico de alimentos visa garantir qualidade e inocuidade pela inativação de micro-organismos e enzimas. Geralmente ocorre sobreprocessamento, levando à perda de características nutricionais e sensoriais e à elevação de custos operacionais. Foi desenvolvida a modelagem matemática (modelos térmico, mássico, hidráulico e reativo) de um equipamento bitubular para o processamento térmico contínuo de um alimento líquido homogêneo em regime laminar difusivo, para determinação das distribuições de temperatura e letalidade. No modelo foram incorporados parâmetros para quantificar os efeitos das difusividades térmica e mássica efetivas na direção radial ao escoamento. O modelo contempla três seções: aquecimento; tubo de retenção e resfriamento, em que as seções de aquecimento e resfriamento são trocadores de calor modulares. Foram considerados os seguintes volumes de controle: alimento, tubo interno, fluido de utilidade, tubo externo, isolamento térmico e ar ambiente. Balanços diferenciais de energia e massa foram aplicados para modelar as trocas térmicas e a distribuição de micro-organismos, enzimas ou nutrientes. Condições de contorno foram aplicadas para manter a continuidade entre as seções e módulos. O modelo foi resolvido através do método de diferenças finitas, discretizando as componentes axiais e radiais. No estudo de caso foram feitas simulações para o processamento térmico de suco de amora, quantificando os efeitos dos parâmetros difusivos e das variáveis de processo sobre os atributos de qualidade nutricional (antocianina) e qualidade microbiológica (levedura). Na validação dos modelos térmico e mássico os efeitos de mistura e difusão na direção radial foram quantificados através do ajuste da condutividade térmica efetiva e número de Peclet modificado. Utilizou-se um pasteurizador em escala laboratorial e os fluidos: mistura glicerina/água 80% (newtoniano) e solução de CMC 1% (pseudoplástico) nas vazões de 10 a 50 L/h. Os resultados indicam que o modelo proposto representou satisfatoriamente as distribuições de temperaturas, letalidade e comportamento experimental com reduzida complexidade computacional. / The thermal processing ensures quality and safety by inactivation of micro-organisms and enzymes. Generally, over-processing occurs, resulting in sensorial and nutritional losses and operating costs increase. A mathematical model (thermal, mass, hydraulic and reactive) of a double-pipe pasteurizer has developed for the continuous thermal processing of a liquid food in laminar diffusive regime, to obtain the temperature and lethality distributions. Parameters to quantify the effects of thermal and mass effective diffusivities in the radial direction were incorporated in the model. The model comprises three zones: heating; holding tube and cooling, where the heating and cooling zones are modular double-pipe heat exchangers. The following control volumes were considered: food, inner tube, utility fluid, outer tube, thermal insulation and ambient air. Energy and mass differential balances were used for modeling the heat transfer and the destruction of micro-organisms, enzymes or nutrients. Boundary conditions were applied to maintain the continuity between zones and modules. The finite difference method was used to discretize the model in the axial and radial directions. A simulation study case for blueberry juice thermal processing was made, analyzing the effects of the diffusive parameters and processing variables of the thermal treatment in the nutritional quality (anthocyanin) and microbiological quality (yeasts) attributes. In the validations of the heat and mass models, the mixture and diffusive effects in the radial direction were quantified by adjusting the heat conductivity and the modified Peclet number. A laboratory scale pasteurizer was used with the following fluids: glycerin/water mixture 80% (newtonian) and CMC 1% solution (pseudoplastic) in flows rates between 10 and 50 L/h. The obtained results indicate that the proposed model successfully represented the temperature and lethality distributions and experimental behavior with reduced computational complexity.

Estudamos a dinâmica de formação de bolhas em fluidos viscosos. Desenvolvemos novas técnicas experimentais para aprofundar o estudo do sistema de um bico borbulhador no fundo de um cilindro preenchido por uma solução de água e glicerina. Neste sistema, controlamos o fluxo de ar por meio de um fluxímetro de massa e uma válvula solenoide. Medimos simultaneamente: o tempo entre duas bolhas sucessivas usando um laser focalizado em um foto- diodo, cuja luz é espalhada com a passagem das bolhas; as ondas de pressão que o borbulhamento causa no sistema de injeção usando microfones colocados entre a válvula solenoide e o bico; e a geometria das bolhas em formação, através da análise de imagens coletadas por uma câmera de alta velocidade. Motivados pelos resultados deste sistema, construímos um modelo para a dinâmica de formação de bolhas, baseado em primeiros princípios, que descreve a rota de adição de períodos observada no experimento. Investigamos também um novo sistema experimental, que tem dois bicos borbulhadores dentro do mesmo tubo. Usamos as mesmas técnicas experimentais desenvolvidas para o experimento de um bico. Identificamos que os borbulhamentos sincronizam apesar de seus sistemas de injeção de ar serem independentes, concluímos, portanto, que o acoplamento deve-se à circulação que o líquido introduz. Estudamos a sincronização por meio da construção de diagramas de bifurcação e espaços de parâmetros mostrando as regiões de sincronização. Adaptamos também o modelo de um bico para descrever este sistema. Descobrimos que a interação entre as duas dinâmicas de formação de bolhas é maior conforme aumentamos a altura da coluna líquida. / We studied the bubble formation dynamics in viscous fluids. We developed new experimental techniques to provide a better understanding of the system with one bubbling nozzle on the bottom of a cylinder filled with a water and glycerol solution. In this system, we control the air flow with a mass flowme- ter and a solenoid valve. We have measured simultaneously: the time between successive bubbles, using a laser focused in a photo-diode that has his light scattered when one bubble passes between them; the pressure waves that the bubbling causes in the air injecting system using microphones placed between the solenoid valve and the nozzle; and the geometry of the growing bubbles, through the analysis of the images collected by a high speed camera. The re- sults from this system motivated us to make a model for the bubbles formation dynamics that describes the period-adding route observed in the experiment and is based on first principles. We also investigated a new experimental system, that has two bubbling nozzles inside the same cylinder. We used the same experimetal techniques developed to the one nozzle experiment. We identified that the bubbling on the nozzles synchronizes due the coupling introduced by the liquid circulation, as the air injecting systems are independent. We studied the synchronization by constructing bifurcation diagrams and parameter spaces that show the syn- chronized regions. We also adapted the one nozzle bubbling model to decribe this system. We discovered that the interaction between the two bubble forma- tion dynamics is stronger when we increase the liquid height.

O objetivo deste estudo foi avaliar o efeito de formulações de rações para frangos de corte utilizando programação linear e não linear sobre desempenho, características de carcaça e aspectos econômicos. Um mil e oitenta machos da linhagem Cobb 500™, com um dia de idade, foram distribuídos em delineamento em blocos casualizados, com seis tratamentos e seis repetições de 30 aves cada. Os tratamentos foram: 1) ração formulada por programação linear com exigência nutricional das Tabelas Brasileiras para Aves e Suínos; 2) ração formulada por programação linear com exigência nutricional do Manual de Manejo Cobb™; 3) ração formulada por programação linear com exigência nutricional de modelos matemáticos; 4) ração formulada por programação não linear em cenário normal de mercado; 5) ração formulada por programação não linear em cenário favorável de mercado; e 6) ração formulada por programação não linear em cenário desfavorável de mercado. Considerou-se cenário de mercado normal os valores do milho, do farelo de soja e do preço pago pelo kg de frango vivo em 26/03/2015, para o interior do Estado de São Paulo. Para o cenário favorável, os valores dos ingredientes baixaram 10% e o preço pago pelo kg de frango vivo subiu 10% e, para o cenário desfavorável, ocorreu o inverso. Foram avaliadas as características de desempenho e econômicas aos 21 e 42 dias de idade e também as características de carcaça aos 42 dias. Os dados foram submetidos à análise de variância e, em caso de significância (P<0,05), aplicou-se o teste de Scott-Knott (5%). Os melhores resultados de desempenho e de características de carcaça foram associados aos tratamentos oriundos da programação linear, cujas concentrações de nutrientes foram mais elevadas, dos quais proteína bruta, aminoácidos e fósforo disponível são determinantes do desempenho superior de frangos de corte. Os tratamentos oriundos da programação não linear geraram melhores resultados econômicos, principalmente no cenário de mercado favorável, que apontou o menor custo para produção de 1 kg de peso corporal e maior lucro. Conclui-se que a formulação de ração com o uso da programação não linear, embora não maximize o desempenho, possibilita a obtenção de rações que otimizam o lucro. / The aim of this study was to evaluate the effect of feed formulations for broilers using linear and nonlinear programming on performance, carcass characteristics and economic aspects. One thousand and eighty, one-day-old male broilers, Cobb 500 ™ strain, were distributed in a randomized block design, with six treatments and six replicates of 30 birds each. The treatments were: 1) feed formulated by linear programming with nutritional requirement of the Brazilian Tables for Poultry and Swine; 2) feed formulated by linear programming with nutritional requirement of Management Cobb™ Manual; 3) feed formulated by linear programming with nutritional requirement of mathematical models; 4) feed formulated by nonlinear programming on normal market scenario; 5) feed formulated by nonlinear programming on favorable market scenario; and 6) feed formulated by nonlinear programming on unfavorable market scenario. In the normal market scenario values of corn, soybean meal and the price paid for kg of live broiler were that practiced on March 26, 2015, for the state of São Paulo. To the favorable scenario, the values of the ingredients decreased 10% and the price paid for kg of live broiler increased by 10% and, for the unfavorable scenario, the opposite happened. The performance and economic characteristics at 21 and 42 days, and carcass characteristics at 42 days were evaluated. Data were submitted to analysis of variance and, in case of significance (P<0.05), the Scott-Knott test (5%) were applied. Best performance and carcass characteristics results were associated with treatments derived from linear programming, whose nutrient concentrations were higher, of which crude protein, amino acids and available phosphorus are determining the superior performance of broilers. Treatments arising from the non-linear programming generated better economic results, especially in the favorable market scenario, which showed the lowest cost to produce 1 kg of body weight and higher profit. It was concluded that feed formulation using non-linear programming, even though does not maximize performance, allows to feed formulation that optimize profit.

Estudos aplicados ao processo de filtração tangencial têm sido o objetivo de muitos trabalhos devido a sua ampla capacidade de adequação aos mais diversos processos industriais. Os mecanismos de transferência associados a esse processo envolvem basicamente o escoamento em tubos permeáveis (ou membranas) com fenômenos de transferência de massa presentes junto à superfície. Esta tese de doutorado apresenta uma modelagem numérica capaz de representar o escoamento de fluidos newtonianos e não-newtonianos em tubos permeáveis aplicada ao processo de filtração tangencial, para os regimes de escoamento laminar e turbulento. No regime turbulento, utilizou-se o modelo de comprimento de mistura de Prandtl. O modelo numérico envolveu as equações de conservação da massa, da quantidade de movimento e das espécies químicas, acopladas a adequadas condições de fronteira. Em particular, por se tratar de um modelo numérico, analisou-se a discretização dos termos convectivos por meio da implementação de três esquemas: WACEB, CUBISTA e QUICKEST adaptativo. Os resultados produzidos pela implementação dos esquemas convectivos foram analisados por meio de comparações com resultados analíticos e experimentais da literatura, para regimes de escoamento laminar e turbulento. De acordo com as comparações realizadas, o esquema QUICKEST adaptativo apresentou melhor desempenho na modelagem desse escoamento. Diversas simulações numéricas geraram resultados, os quais foram comparados com expressões analíticas e dados experimentais da literatura e com dados produzidos pelo laboratório do Núcleo de Engenharia Térmica e Fluidos (NETeF) da USP São Carlos. Verificou-se, assim, que o modelo matemático produziu resultados compatíveis com o fenômeno estudado, tendo-se, portanto, uma ferramenta para descrição do problema convectivo mássico do escoamento em tubos permeáveis. / Studies applied to the crossflow filtration process have been the goal of several researches due to its large capacity for adaptation to the various industrial processes. The mechanism of transfer associated with this process basically involves fiows in permeable tubes (or membranes) with the phenomenon of mass transfer on the surface. This doctoral thesis presents a numerical modeling able to represent the fiow of Newtonian and non-Newtonian fiuids in permeable tubes applied to the process of crossflow filtration, in laminar and turbulent flow regimes. In the turbulent regime the Prandtl mixing length model was used. The numerical model involved the mass conservation, momentum conservation and mass transport equations coupled to the appropriate boundary conditions. In particular, as it is a numerical model, the discretization of convective terms was analyzed through the implementation of three schemes: WACEB, CUBISTA and adaptative QUICKEST. The results produced by the implementation of the convective schemes were analyzed through comparisons with experimental and analytical results from the literature, for laminar and turbulent flow regImes. According to the comparisons made, adaptative QUICKEST scheme showed the best performance in the modeling in this flow. Several numerical simulations generated results, which were compared with experimental data and analytical expressions of the literature and with data produced by the laboratory of the Center of Thermal Thermal and Fluids Engineering (NETeF) of the EESC- USP/São Cados. It was possible to verify, that the mathematical model produced results consistent with the phenomenon studied, and can be considered a tool for the description of the mass convective problem flow in permeavel tubes.

Nos últimos anos, o ensino de fisiologia tem sido fortemente beneficiado pelo desenvolvimento de modelos matemáticos e simuladores de paciente capazes de reproduzir com segurança partes específicas ou sistemas fisiológicos completos. Estudos mostram que aulas teóricas aliadas a simulações conseguem potencializar o nível de compreensão dos conceitos envolvidos. O principal objetivo deste trabalho consiste no desenvolvimento de um simulador do sistema cardiovascular, capaz de representar de forma didática a dinâmica do ciclo cardíaco. A disponibilidade de simuladores deste tipo é reduzida, e limitada a modelos simplificados ou interfaces pouco amigáveis. De modo a conseguir maior flexibilidade nas situações simuladas e um adequado grau de proximidade com o sistema real, decidiu-se pela utilização de um modelo matemático da fisiologia cardiovascular para o cálculo das variáveis e parâmetros atuantes no plano de fundo do simulador. O desenvolvimento deste modelo constitui o segundo objetivo deste trabalho. Devido ao grau de complexidade desejado, optou-se pela utilização de um modelo existente na literatura, acrescido dos parâmetros que fossem considerados ausentes. O resultado foi um modelo com oito compartimentos vasculares e quatro compartimentos representando as câmaras do coração, todas com atividade pulsátil, controlados pelo mecanismo reflexo baroceptor, controle de curto prazo que atua na regulação da pressão arterial. O simulador (desenvolvido em Visual C# com interface em Macromedia Flash) permite a alteração individual dos parâmetros vasculares e a simulação de hemorragias, bem como a visualização de gráficos de pressão, fluxo e volume em qualquer compartimento, e a construção de alças pressão-volume para os átrios e ventrículos. Além disso, o aplicativo resultante possibilita a inclusão futura de novos recursos e ferramentas, como tutoriais e simulação de patologias. / During the last years, physiology learning has been strongly favored by the development of mathematical models and patient simulators capable of safely reproduce specific parts or complete physiological systems. Studies show that traditional classes together with simulations are able to in-crease the comprehension of concepts involved. The main objective of this work is to develop a simu-lator of the cardiovascular system capable of representing the cardiac cycle dynamics in a didactic fashion. There are just a few available simulators of this kind, and the existing ones are based on simplified models or unfriendly interfaces. In order to achieve high flexibility in simulated scenarios and an adequate level of realism, we decided to use a mathematical model of the cardiovascular physiology to calculate variables and parameters acting in the background of the application. The development of this model constitutes the second objective of the present work. Due to the desired level of complexity, we decided to use an existing model found in the literature, improved with addi-tional parameters we found necessary. The result is a model with eight vascular compartments and all four cardiac chambers, with pulsatile behavior, controlled by the baroreflex mechanism, a short term control that regulates arterial pressure. The simulator (developed in Visual C# with an interface built in Macromedia Flash) allows the user to change individual parameters and simulate blood losses, as well as visualize press, flow and volume graphs from any compartment and also pressure-volume loops from the cardiac chambers. Moreover, the resulting application is open to future inclusion of new resources and tools, such as tutorials and pathology simulation.

Nas indústrias de celulose, o processo dominante é a polpação Kraft. A polpação ou o cozimento de madeira em digestores constitui a operação central da produção de celulose, onde se processam as reações de deslignificação dos cavacos de madeira. Em geral, aí os modelos são empíricos e baseados em experimentos em escala de laboratório, e o \"scale-up\" para nível de digestores industriais torna-se tedioso e às vezes impraticável, pelo número de experimentos necessários e pela dificuldade de reprodução de todas as condições operacionais. Propõe-se neste trabalho um modelo matemático de polpação Kraft de eucalipto. O modelo consiste num sistema de equações diferenciais parciais que descrevem os efeitos combinados da cinética de polpação e da difusão de reagentes, como um reator de mistura perfeita com fases heterogêneas. Apresentam-se as previsões das simulações, comparam-se seus valores com dados experimentais. Este modelo é heterogêneo e difere dos modelos homogêneos usualmente propostos para a polpação. O primeiros permite uma abordagem das influências da espessura, massa específica e porosidade de cavacos, e também da difusão de reagentes nos cavacos durante a polpação. O presente trabalho dá ainda alguns exemplos de aplicação e tece comentários quanto ao modelo. / Kraft pulping is the dominant process f the pulp industry. The pulping of wood in digesters is the central operation of the pulp production, which consists of the delignification of the wood chips. Usually, the models of the delignification are empiricals and derived from laboratory scale experiments . The scale-up for industrial applications is tedious and sometimes impracticable, requiring a large number of experiments and presenting difficulty in reproducing all the operational condictions. The present work suggests a mathematical model for the Kraft pulping of the eucalyptus. This model consists of a system of partial differential equations describing the combined effects of the pulping kinetics and the diffusion of reactants, in form of continuous stirred tank reactor with heterogeneous phases. The prediction of the values given by simulation and their comparison with experimental data are presented. Differing from the homogeneous model, the heterogeneous allows an approach to estimate the influence of the chip thickness, density, porosity and that of the diffusion of the reactants in the chips, during pulping. The present work gives further some examples of application and comments on the model.