Resumo não disponível.

Para a pesquisa, utilizou-se um concreto asfáltico com ligante modificado com borracha ensaiado nas dependências do LAPAV-UFRGS, sob compressão diametral de cargas cíclicas, cujo procedimento é brevemente descrito, associando-se os seus procedimentos com as hipóteses assumidas na modelagem, o que é importante na concepção do modelo empregado. Apresenta-se e discute-se o manuseio dos dados experimentais para o modelo, bem como a escolha do modelo constitutivo adotado, no qual se optou por um modelo viscoelástico linear anisotrópico simplificado, cujas equações base e métodos para calibração são apresentados, discutidos e comparados, bem como os resultados das calibrações dos três primeiros ciclos de carga-descanso, no qual se obteve uma excelente descrição do comportamento do material, o que foi comprovado nas previsões dos dois ciclos seguintes. As previsões de resultados a ensaios sob diferentes condições aos das calibrações mostraram resultados qualitativamente bons às freqüências de ciclo de 1 Hz, comprovados na comparação das curvas de deslocamentos normalizadas. Para a pesquisa, utilizou-se um concreto asfáltico com ligante modificado com borracha ensaiado nas dependências do LAPAV-UFRGS, sob compressão diametral de cargas cíclicas, cujo procedimento é brevemente descrito, associando-se os seus procedimentos com as hipóteses assumidas na modelagem, o que é importante na concepção do modelo empregado. Apresenta-se e discute-se o manuseio dos dados experimentais para o modelo, bem como a escolha do modelo constitutivo adotado, no qual se optou por um modelo viscoelástico linear anisotrópico simplificado, cujas equações base e métodos para calibração são apresentados, discutidos e comparados, bem como os resultados das calibrações dos três primeiros ciclos de carga-descanso, no qual se obteve uma excelente descrição do comportamento do material, o que foi comprovado nas previsões dos dois ciclos seguintes. Para a pesquisa, utilizou-se um concreto asfáltico com ligante modificado com borracha ensaiado nas dependências do LAPAV-UFRGS, sob compressão diametral de cargas cíclicas, cujo procedimento é brevemente descrito, associando-se os seus procedimentos com as hipóteses assumidas na modelagem, o que é importante na concepção do modelo empregado. Apresenta-se e discute-se o manuseio dos dados experimentais para o modelo, bem como a escolha do modelo constitutivo adotado, no qual se optou por um modelo viscoelástico linear anisotrópico simplificado, cujas equações base e métodos para calibração são apresentados, discutidos e comparados, bem como os resultados das calibrações dos três primeiros ciclos de carga-descanso, no qual se obteve uma excelente descrição do comportamento do material, o que foi comprovado nas previsões dos dois ciclos seguintes.

Esse trabalho tem por objetivo geral o desenvolvimento de uma metodologia de modelagem numérica que represente o escoamento e o fenômeno de mistura em um modelo físico de panela de aciaria de base elíptica. Os objetivos específicos do trabalho são: o estudo dos coeficientes das forças de interação líquido-bolha, dos modelos de turbulência e da mudança do formato da base da panela de circular para elíptica. O escoamento e o fenômeno de mistura foram calculados através do método de Volume Finitos baseado em Elementos por meio do software CFX5.7 da Ansys. Dados da literatura e ensaios em modelo físico, realizados em laboratório, auxiliaram na validação dos modelos numéricos. O estudo dos coeficientes das forças de não-arrasto mostrou que os resultados da distribuição de ar ao longo da altura do banho mudam com a variação dos coeficientes. No final, coeficientes para 3 configurações de panelas em diferentes vazões de ar foram encontrados. Com relação ao estudo dos modelos de turbulência, observou-se que para a solução do escoamento e do fenômeno de mistura em uma panela de base circular, o k-ε é o modelo de turbulência mais indicado. Por outro lado, para uma panela de base elíptica, o modelo RSM mostrou-se o mais adequado. Finalmente, com relação ao estudo da mudança do formato da base da panela, observou-se que os tempos de mistura de uma panela de base elíptica são maiores que uma de base circular e aumentam à medida que a vazão de ar diminui.

No estudo da propagação de uma doença infecciosa, diz-se que sua transmissão ocorre horizontalmente, quando um indivíduo suscetível tem um contato direto ou indireto com um indivíduo infeccioso. Algumas doenças, entretanto, também podem ser transmitidas verticalmente, entendendo-se que, neste caso, a doença é transmitida a um indivíduo, ao ser gerado por uma mãe infecciosa. Fazendo uso de modelos epidemiológicos determinísticos básicos, envolvendo sistemas de equações diferenciais ordinárias, nosso principal objetivo, neste trabalho, consiste em investigar qual o papel da transmissão vertical na propagação de doenças causadas por microparasitas. Diversas formas de inclusão de transmissão vertical são apresentadas e, em cada modelo estudado, investigamos a existência e a estabilidade local dos estados de equilíbrio da população hospedeira, identificamos os parâmetros e limiares que caracterizam a dinâmica do sistema, e completamos as informações decorrentes dos resultados analíticos com a apresentação de soluções numéricas do mesmo. Por fim, comparamos os efeitos da transmissão horizontal com aqueles decorrentes da transmissão vertical.

Neste trabalho nos concentraremos na análise de.um aspecto particular dos processos econômicos: modelar a competição entre diferentes agentes em um ambiente onde as trocas de recursos entre eles acontecem em um modo conservativo, ou seja, um modelo de mercado de trocas conservativas (MMTC). Nós apresentamos um modelo simplificado para a exploração de recursos por agentes econômicos interagentes,onde cada agente recebe uma fração randômica dos recursos disponíveis. Em cada passo de tempo, o agente mais pobre, isto é, aquele com a menor riqueza do mercado, sofre uma mudança randômica na sua riqueza.. Essa dinâmica de mínimo garante que o agente mais pobre tenha chance de mudar seu estado econômico. Qualquer riqueza que é ganha (ou perdida) por este agente é debitada (ou creditada) entre seus dois primeiros vizinhos na rede, de modo que a riqueza total do mercado permanece constante. Depois de um longo transiente o sistema se auto-organiza em um estado crítico auto-organizado (SOC) que maximiza a performance média de cada participante. No estado SOC quase todos os agentes têm riqueza acima de um certo limiar, um valor de riqueza crítico nt = 0,4, que é chamado de linha de pobreza. Acima deste limiar, nosso modelo exibe um novo tipo de condensação de riqueza; uma distriibuição de riqueza exponencial onde muito poucos agentes extremamente ricos são estáveis no tempo, enquanto o restante permanece no que chamamos de classe média. A redistribuição de riqueza é então evidente. O progresso econômico na sociedade é uniforme, embora lento. A evolução temporal dos eventos de redistribuição de riqueza entre os agentes segue o padrão de equilíbrio pontuado, com intervalos relativamente longos de tranquilidade intercalados por períodos de explosões de atividade. Além disso, a distribuição da duração das avalanches segueuma lei de potência com expoente o - 1 que corresponde a um espectro de freqüência 1/f que se refere a correlações temporais extremamente longas. Ambos os fenômenos são característicos de sistemas no estado SOC. Estudamos também uma solução não geográfica de campo médio (interações globais), que gera umu distribuição de riqueza completamente diferente que é quase linear acima de um limiar crítico, nt = 0,2, (bem menor do que no caso de interações locais) A globalização aumentou as diferenças econômicas, pois houve um aumento do número de agentes mais ricos e também um crecismento da miséria dos mais pobres, que se tornaram em maior número. Neste caso se pode identificar claramente 3 classes sociais: a rica, a classe média e a pobre. De certo modo, o modelo local (não globalizado) corresponde a um tipo de mundo feudal, onde os barões locais mantêm o seu domínio por longos períodos de tempo, sobre uma popução de agentes semipobres. Com isso, o MMTC apresentado aqui fornece uma descrição simples da redistribuição de riqueza nos estágios iniciais da história econômica da humanidade, e indica algumas das possíveis forças motrizes que influenciaram esse processo de redistribuição. Nós acreditamos que essas conclusões podem ser de interesse; em vista do presente debate sobre os aspectos positivos e negativos da globalização. / In this work we concentrate our analysis in a particular aspect of the economic processes: modelling the competition among different agents in an environment where all exchanges of resources between agents take place in a conservative manner, that is, a conservative exchange market model (CEMM). We present a simplified model for the exploitation of resources by interacting economic agents, where eache agent receives a random fraction of the available resources. At each time step, the poorest agent, i.e., the one with the minimum wealth in the market, suffers a random change in its wealth. This minimum dynamic ensures that the poorest agent has a chance to change its economic welfare. Whatever wealth is gained (or lost) by the poorest agent is equally debited (or credited) between its two nearest neighbors in the lattice, making the total wealth constant. After a long transient, the system self-organizes into a self-organized critical state (SOC) that maximizes the average performance of each participant. In the SOC state almost all agents have wealth beyond a certain thershold, a critical value of wealth n=0,4, that is called poverty line. Above this threshold our model exhibits a new kind of wealth condensation; an exponential wealth distribution where the very few extremelly rich agents are stable in time, while the rest remains in what we call a middle class. Wealth redistribution is then evident. The economic progress in the society is steady, even if slow. The temporal evolution of the wealth redistribution events among the agents follows a punctuated equilibrium pattern, with relatively long periods of stasis separated by bursts of activity. Moreover, the distribution of avalanches lifetime follows a power law with exponent a - 1, which corresponds to a 1/f frequency spectrum, that is related to extremelly long temporal correlations. Both phenomena are characteristic of systems in the SOC state. We also study the a-geographic mean field solution (global interactions), wich generates a completely different wealth distribution that is almost linear beyond a critical threshold, nr=0,2, (much lower than in the local interaction case). Globalization increases the economic inequalities, because there is a relatively great fraction of wealth agents but there is also a higher number of poorer agents with lower income. In this case it is possible to clearly identify three social classes: the rich, the middle, and the poor one. In a sense, the local model (non-globalized) corresponds to a kind of feudal world, where local barons mantain their dominance for long periods of time, over a population of semi-poor agents. Therefore, the CEMM presented here provides a simple description of wealth redistribution in the early stages of human economic history, and indicates some of the possible driving forces that influenced this redistribution process. We believe that these conclusions may be of interest in view of the present debate over the goods and evils of globalization.

O presente trabalho tem por objetivo obter uma metodologia para a simulação do escoamento no interior da bifurcação da artéria carótida. Como a análise completa deste tipo de fluxo não é uma tarefa óbvia. simplificações quanto ao fluido e método numérico são necessárias. Desta forma, o escoamento analisado no presente trabalho, no interior da artéria carótida, é considerado como sendo bidimensional, incompressível, viscoso, em regime permanente, com paredes rígidas e newtoniana com viscosidade e massa específ'ica um pouco superiores a da água. Estas hipóteses são baseadas no fato de que o sangue para altas taxas de cisalhamento, como as encontradas na artéria carótida, comporta-se como um fluido newtoniana incompressível semelhante a água; o que não ocorre para baixas taxas de cisalhamento como em algumas outras artérias ou veias do corpo humano. Adota-se, neste caso, o método de diferenças finitas para a formulação incompressível, onde as equações de conservação da quantidade de movimento são utilizadas para obter os componentes do vetor velocidade e a equação da continuidade juntamente com a quantidade do movimento para obter a pressão, baseado no processo de integração temporal de Runge-Kutta de três estágios e com aproximação espacial de segunda ordem. Emprega-se ainda coordenadas generalizadas devido à forma geométrica da carótida ser complexa. Cuidados adicionais são tomados na aplicação das condições de contorno. Resultados numéricos são apresentados para três geometrias: um duto simples, uma artéria carótida e um duto com obstrução (ressalto). O escoamento no interior do duto visa a verificação do código confeccionado pois este problema possui solução analítica. Com o mesmo código analisa-se regiões que apresentam distúrbios do fluxo, com a presença de recirculação e aumento da tensão de cisalhamento, no interior da artéria carótida, e por fim analisa-se a influência do número de Reynolds, variando entre 100 e 500, para um modelo de estenose (estreitamento) no interior de um duto. Observações importantes como redução da velocidade na região de bifurcação e na região do sinus da artéria carótida são verificadas podendo auxiliar a prática da medicina na detecção e/ou prevenção de doenças como aterosclerose que são causas de morte principalmente em adultos. / The goal o f the present work is to obtain a numerical methodology for the tlow simulatíon through a carotid artery bifurcation. As the complete numerical flow analysis is not easy and obvious, simplifications related to the tluid and the numerical method are necessary (or indispensable). In this way, the tlow situation ts considered to be two-dimensionaL incompressible, viscous. steady state, with rigid walls and for a Newtonian fluid with density and viscosity slightly greater than for the water. These hypothesis are based on the fact that the blood. submitted to high strain deformations. as in the carotid artery, behaves as a Newtonian incompressible fluid (it seems like water); what does not occur at many other parts ofthe human body. ln order to solve this flow problem the finite difference method based on the incompressible formulation is employed. In this method the momentum equations are used to obtain the velocity vector components and the continuity equation, together with the momentum equations, to obtain the pressure. This numerical technique is based on the three stages Runge-Kutta time stepping scheme for a second order spatial discretization. Body fitted coordinates are used dueto the geometry complexity ofthe carotid artery. Additional care has to be taken with the implementation ofboundary conditions. Numerical results are obtained for three geometries: a duct, a carotid artery and a blocked duct. The simple duct flow was analyzed because of its geo"\etrical simplicity and in order to verify the code written in cartesian and generalized coordinates. With the same code the carotid artery flow with recirculation zones and high strain rates was analysed and finally the flow into a blocked duct for Reynolds numbers between 100 and 500 (a stenosis model) was analysed. Important observations as the decrease o f velocity at the bifurcation region and in the sinus of the internai carotid artery are obtained which can help medicai practice in the detection and/or prevention of diseases like the atherosclerosis which kills mainly adults.

A demanda mundial de minério de ferro tem aumentado consideravelmente nos últimos anos e mais recentemente tem sofrido queda do valor de venda, provocando a necessidade de melhoria de desempenho das instalações de beneficiamento mineral. A automação de processos é uma das áreas que promove esta melhoria de desempenho, uma vez que as estratégias de controle automatizadas reduzem a variabilidade dos processos a partir da minimização da interferência do operador. A automação mantém o sistema sobre vigília todo o tempo reproduzindo as ações do operador, padronizando-as em uma frequência que não seria possível para seres humanos. Estas estratégias, normalmente, são definidas empiricamente a partir de experiência e do conhecimento dos técnicos de operação. O objetivo geral do trabalho é propor um método científico para definição de estratégia de controle automático para um circuito de moagem, desenvolvido através da utilização de tecnologia de simulação de processos e análise multivariada. Objetivos específicos visam a obtenção de balanço de massas, granulométrico e de água do circuito de moagem (i) e de modelos matemáticos que descrevam o comportamento do circuito de moagem considerando as variáveis críticas (ii) e determinação estruturada de estratégia de controle automatizada de processo (iii). A otimização permitiu atingir a produtividade esperada, acima de 680t/h, a meta para a especificação granulométrica do produto, 50,0m, e atingir a meta de recuperação em massa, próxima de 90%. A sequência utilizada neste trabalho pode ser considerada como uma metodologia científica para definição de estratégias de controle automatizadas para processamento mineral. Este método poderá ser empregado em sistemas de grande importância como, neste caso, a moagem de minério de ferro. / The world demand for iron ore has increased considerably in the past years and most recently the sale value has decreased, causing the need for performance improvement of mineral processing facilities. Process automation is one of the areas that promote this performance improvement since the automated control strategies reduce process variability by minimizing the operator interference keeping the system under constant surveillance all the while reproducing the operators’ actions, standardizing them in a frequency that would not be possible for humans. These strategies are usually defined empirically from experience and knowledge of operating technicians. The overall objective of this study is to propose a scientific method to define an automatic control strategy for a grinding circuit, developed through the use of processes simulation technology and multivariate analysis. Specific objectives aim to obtaining mass balance, particle size distribution and water in the grinding circuit (i), and mathematical models describing the grinding circuit behavior considering the critical variables (ii) and the structured determination of process's automated control strategy (iii). The optimization was able to achieve the expected productivity, above 680th, the target for particle size specification of the product, 50,0m, and achieve mass recovery target, close to 90%.The sequence used in this work can be considered as a scientific methodology to define automated control strategies for mineral processing. This method can be used in systems of great importance as, in this case, the iron ore grinding.

O fenômeno da dispersão de poluentes precisa ser compreendido e analisado pro- fundamente para que estudos de impactos ambientais possam ser realizados para projetar e prever situações. Então, obter um modelo analítico de dispersão de poluente realista torna- se interessante, pois permite avaliar com maior precisão o impacto ambiental da liberação de poluentes na atmosfera, além de abrir novos horizontes para a pesquisa. A Camada Limite Planetária (CLP) é o domínio de interesse da grande maioria dos modelos, já que nela os fenômenos turbulentos estão presentes. Uma característica destes fenômenos são os vórtices e turbilhões chamados de Estruturas Coerentes (ECs) e que são dominados por uma fase. O modelo analítico mais estudado é o que tem por base a equação advecção-difusão onde diversas simplificações como médias de Reynolds e o fechamento Fickiano são aplicados. Essas simplificações tornam o modelo determinístico e linear, mesmo o fenômeno sendo estocástico e não linear. Para resgatar algumas características do fenômeno turbulento neste trabalho sugeriu-se a inclusão de uma fase na equação advecção-difusão, através de um coeficiente difusivo complexo. Isso porque estruturas coerentes são características turbulentas dominadas pela fase. Comparou-se aqui os modelos com coeficiente difusivo real e complexo para as mesmas condições podendo assim observar qualitativamente a inclusão da fase no modelo que reproduz uma característica do escoamento turbulento, apresentando um comportamento mais realista. Apesar de não podemos garantir que esta é a maneira mais acertada de incluir a fase na equação advecção-difusão, com certeza ela nos trás um grande benefício que é a garantia de sempre ter soluções semi-positivas definidas, compatíveis com distribuições, para a representação da concentração. / The phenomenon of dispersion of pollutants needs to be analyzed and deeply un- derstood so that environmental impact studies can be performed to design and predict situations. Then, to obtain an analytical model of realistic pollutant dispersion becomes interesting because it allows to accurately assess the environmental impact of the release of pollutants into the atmosphere and to open new horizons for research. The planetary bound- ary layer (PBL) is the domain of interest of the majority of models, since this is domain where turbulent phenomena are present. A feature of these phenomena are the vortexes and eddies that are Coherent Structures (CSs) and are dominated by a phase. The most studied analytical model is based on the advection-diffusion equation where several simpli cations such as Reynolds averages and Fickian closure are applied. These simpli cations render the model deterministic and linear, although the phenomenon is stochastic and nonlinear. To recover some characteristics of turbulent phenomena in this work a phase was included in the advection-diffusion equation by a complex diffusion coefficient. This is because coherents structures are turbulent characteristics and dominated by phase. We compared the models with real and complex diffusion coefficient for the same conditions, thus being able to observe qualitatively the inclusion of phase in the model that reproduces a characteristic of turbulent ow, presenting a more realistic behavior. Although we can not guarantee that this is the adequate way to include the phase in the advection-diffusion equation, it certainly brings us bene t that is the guarantee to always have at least semi-positive solutions, compatible with distributions that represent the concentration.

A identificação de modelos consiste na atividade de construir modelos matemáticos representativos de sistemas reais. Diversas técnicas estão disponíveis para a sua realização, as quais dependem do nível de informação preliminares disponíveis referentes ao sistema em estudo. Neste trabalho técnicas de identificação caixa-cinza, as quais utilizam estruturas de modelos fenomenológicos e dados experimentais foram estudadas, dando enfoque a duas etapas principais: análise de identificabilidade e estimação de parâmetros. A análise de identificabilidade é uma ferramenta fundamental na identificação e deve preceder a etapa de estimação dos parâmetros. Seu objetivo é verificar, a partir de uma estrutura e de dados experimentais quais parâmetros de um modelo podem ser identificados. Neste trabalho as principais metodologias disponíveis na literatura foram analisadas. Posteriormente, foi desenvolvida uma nova metodologia, fundamentada em conceitos de teoria de controle, a qual visa suprir as limitações inerentes aos demais métodos estudados. A estimação dos parâmetros foi realizada considerando dois métodos de otimização dinâmica distintos: o método single-shooting e método simultâneo. Além disso, nesta mesma etapa diferentes formulações de função objetivo foram testadas, verificando o efeito da minimização do quadrado da derivada do erro e comparando com a abordagem tradicional, a qual não faz uso da derivada. Para ilustrar e testar a aplicação das metodologias, dois estudos de caso foram analisados: Reator de Fermentação Alcóolica Contínua e Reator Semi-Batelada de Williams-Otto. Em ambos os casos, os dados experimentais foram gerados a partir da simulação do modelo sob parâmetros, condições operacionais e inferências conhecidas, a fim de se verificar o efeito de cada uma destas peculiaridades. Para consolidar os estudos desenvolvidos, uma aplicação real também foi testada através da identificação do modelo de um reator de polimerização de PEAD utilizando dados de uma planta industrial. / Model identification is the activity of constructing representative mathematical models of real systems. Several techniques are available for its completion, which depend on the level of information available about the studied system. In this work techniques for gray-box identification, which use both phenomenological models and experimental data, were studied by focusing on two main steps: identifiability analysis and parameter estimation. Identifiability analysis is a fundamental tool in the identification process and must precede the parameter estimation. Your goal is to verify, from a model structure and experimental data which parameters of a model can be identified. In this work the most important methodologies available in the literature were analyzed. Subsequently, a new methodology was developed, based on concepts of control theory, which aims to overcome the limitations inherent in other methods. The parameters estimation was made considering two distinct dynamic optimization methods: the single-shooting method and simultaneous method. Furthermore, in this same step, several formulations of objective function were tested, checking the effect of minimizing the square of the derivative of the error when compared to the traditional approach, which does not make use of the derivative. To elucidate the methodology, two case studies were proposed: Continuous Alcoholic Fermentation Reactor and Williams-Otto Semi-Batch Reactor. In both cases, the experimental data were generated from the simulation of the model under know parameters, operational conditions and inference in order to verify the effect of each of these peculiarities. To consolidate the studies developed a real application was also tested, and this identifies the model of a polymerization reactor HDPE using data from an industrial plant.

Biopolímeros polihidroxialcanoatos (PHAs) são poliésteres sintetizados por inúmeras bactérias e armazenados em corpos de inclusão nas células como reserva de energia, são produzidos a partir de fontes renováveis, e sua principal característica é a biodegradabilidade. Neste trabalho realizou-se a produção do PHA poli(3- hidroxibutirato) (P3HB) em cultivos submersos da bactéria Bacillus megaterium em meio mineral, utilizando-se sacarose como fonte de carbono e sulfato de amônio como fonte de nitrogênio. Esta bactéria apresenta comportamento diferente das bactérias comumente utilizadas para produção de biopolímero, sendo capaz de produzir grandes quantidades de P(3HB) sem limitação de nitrogênio. Experimentos em estufa incubadora rotatória foram realizados, segundo planejamento fatorial 22 composto central rotacional e três repetições no ponto central, a fim de se avaliar a melhor concentração inicial de sacarose e a razão inicial carbono/nitrogênio. Nestes experimentos, obteve-se um acúmulo de P(3HB) nas células de cerca de 75%. Não houve limitação de nitrogênio e sacarose para a maioria dos casos e observou-se uma forte correlação entre o acúmulo de P(3HB) e o pH do meio. Em biorreator foram realizados experimentos com controle de oxigênio dissolvido e pH, e com a condição de substrato inicial escolhida em função do planejamento de experimentos realizado. Obteve-se acúmulo de P(3HB) de 30%, inferior ao obtido em estufa incubadora rotatória. Modelos fenomenológicos de cinética microbiana, produção de polímero e consumo de substratos foram estudados, propondo-se modificações. Os parâmetros dos modelos foram estimados com os dados obtidos experimentalmente no biorreator. O modelo proposto obteve um desempenho superior aos demais estudados e com um número menor de parâmetros. Simulações foram realizadas a fim de verificar o comportamento do modelo frente a variações nas condições iniciais. / Polyhydroxyalkanoates (PHAs) biopolymers are polyesters synthesized by numerous bacteria and are accumulated on inclusion bodies in the cytoplasm of the cells as energy storage, they are produced from renewable carbon sources and its manly characteristic is the biodegradability. In this work it was carried the production of PHA poly(3-hydroxybutyrate) (P3HB) out in submerse cultures of bacteria Bacillus megaterium in the mineral medium, using sucrose as carbon source and ammonium sulphate as nitrogen source. This bacteria shows different behavior from bacteria commonly used on biopolymer production, and it is able to produce large amounts of P(3HB) without nitrogen limitation. Shaker experiments were realized, according to a central composite design and three repetitions in the central point, in order to evaluate the best initial sucrose concentration and carbon-to-nitrogen ratio. In these experiments, P(3HB) was accumulated in the cell about 75%. In most of the cases, there was no limitation in the sucrose or nitrogen, and a strong correlation between P(3HB) accumulated and the medium pH was observed. Bioreactor experiments were carried out with controlled dissolved oxygen and pH. Substrate initial concentration was chosen according to the experiment design. P(3HB) accumulation in bioreactor was lower than in shaker, about 30%. Phenomenological models of microbial kinetics, polymer production, and substrate consumption were studied and improvements were proposed. The models parameters were estimated using the experiment data from the bioreactor. The proposed model showed a better performance than the available models, with a smaller number of parameters. Simulations were performed with the purpose of verifying the model behavior for initial conditions changes.