Computational studies for the design of process equipment with complex geometries

This study combines several projects related to the flows in vessels with complex shapes representing different chemical apparata. Three major cases were studied. The first one is a two-phase plate reactor with a complex structure of intersecting micro channels engraved on one plate which is covered by another plain plate. The second case is a tubular microreactor, consisting of two subcases. The first subcase is a multi-channel two-component commercial micromixer (slit interdigital) used to mix two liquid reagents before they enter the reactor. The second subcase is a micro-tube, where the distribution of the heat generated by the reaction was studied. The third case is a conventionally packed column. However, flow, reactions or mass transfer were not modeled. Instead, the research focused on how to describe mathematically the realistic geometry of the column packing, which is rather random and can not be created using conventional computeraided design or engineering (CAD/CAE) methods. Several modeling approaches were used to describe the performance of the processes in the considered vessels. Computational fluid dynamics (CFD) was used to describe the details of the flow in the plate microreactor and micromixer. A space-averaged mass transfer model based on Fick’s law was used to describe the exchange of the species through the gas-liquid interface in the microreactor. This model utilized data, namely the values of the interfacial area, obtained by the corresponding CFD model. A common heat transfer model was used to find the heat distribution in the micro-tube. To generate the column packing, an additional multibody dynamic model was implemented. Auxiliary simulation was carried out to determine the position and orientation of every packing element in the column. This data was then exported into a CAD system to generate desirable geometry, which could further be used for CFD simulations. The results demonstrated that the CFD model of the microreactor could predict the flow pattern well enough and agreed with experiments. The mass transfer model allowed to estimate the mass transfer coefficient. Modeling for the second case showed that the flow in the micromixer and the heat transfer in the tube could be excluded from the larger model which describes the chemical kinetics in the reactor. Results of the third case demonstrated that the auxiliary simulation could successfully generate complex random packing not only for the column but also for other similar cases.

Multicomponent diffusion during Prato cheese ripening: mathematical modeling using the finite element method

The partial replacement of NaCl by KCl is a promising alternative to produce a cheese with lower sodium content since KCl does not change the final quality of the cheese product. In order to assure proper salt proportions, mathematical models are employed to control the product process and simulate the multicomponent diffusion during the reduced salt cheese ripening period. The generalized Fick's Second Law is widely accepted as the primary mass transfer model within solid foods. The Finite Element Method (FEM) was used to solve the system of differential equations formed. Therefore, a NaCl and KCl multicomponent diffusion was simulated using a 20% (w/w) static brine with 70% NaCl and 30% KCl during Prato cheese (a Brazilian semi-hard cheese) salting and ripening. The theoretical results were compared with experimental data, and indicated that the deviation was 4.43% for NaCl and 4.72% for KCl validating the proposed model for the production of good quality, reduced-sodium cheeses.

Intégration du concept de capacité de support d’un plan d’eau aux apports en phosphore à l’aménagement du territoire au Québec : Réalité ou utopie?

Ce mémoire a comme objectif de présenter une revue de la littérature et une analyse comparative des différents modèles existants pour le calcul de la capacité de support d’un plan d’eau aux apports en phosphore. Ce document a aussi pour but d’évaluer les différentes options pour l’utilisation de ce genre d’outil au Québec, ainsi que d’amorcer une réflexion sur les applications possibles de ce concept au niveau local, en lien avec l’aménagement du territoire des bassins versants. L’analyse comparative des modèles théoriques, combinée à la réflexion d’acteurs clés impliqués actuellement dans un processus qui consiste à tester leur utilisation dans les Laurentides au Québec, a permis de mettre en relief l’importance de calibrer les modèles régionalement. De plus, certains avantages semblent résider dans l’utilisation d’un modèle de nature empirique afin d’effectuer les prédictions sur les apports naturels et totaux en phosphore pour les lacs de cette région. Par contre, l’utilisation d’une approche de type « bilan de masse », s’avère tout de même indispensable afin de relativiser l’importance des différents apports anthropiques en phosphore du bassin versant. Dans l’avenir, l’utilisation de tels modèles permettra possiblement de justifier certaines mesures restrictives au développement des bassins versants, qui pourront s’insérer parmi les outils d’urbanisme actuels. Ce sont principalement les municipalités qui détiennent les pouvoirs nécessaires afin d’intégrer les prescriptions découlant de ce genre d’analyse à l’aménagement du territoire des bassins versants. Plusieurs outils d’urbanisme, tels que les plans d’aménagement d’ensemble (PAE) ou bien les règlements de zonage et de lotissement, semblent donner assez de latitude afin de permettre aux municipalités d’intervenir en ce sens. Toutefois, les modèles de capacité de support ne pourront pas être utilisés afin de calculer un nombre précis d’habitations qu’il est possible de construire dans le bassin versant d’un lac.

Approche cartésienne pour le calcul du vent en terrain complexe avec application à la propagation des feux de forêt

La méthode de projection et l'approche variationnelle de Sasaki sont deux techniques permettant d'obtenir un champ vectoriel à divergence nulle à partir d'un champ initial quelconque. Pour une vitesse d'un vent en haute altitude, un champ de vitesse sur une grille décalée est généré au-dessus d'une topographie donnée par une fonction analytique. L'approche cartésienne nommée Embedded Boundary Method est utilisée pour résoudre une équation de Poisson découlant de la projection sur un domaine irrégulier avec des conditions aux limites mixtes. La solution obtenue permet de corriger le champ initial afin d'obtenir un champ respectant la loi de conservation de la masse et prenant également en compte les effets dûs à la géométrie du terrain. Le champ de vitesse ainsi généré permettra de propager un feu de forêt sur la topographie à l'aide de la méthode iso-niveaux. L'algorithme est décrit pour le cas en deux et trois dimensions et des tests de convergence sont effectués.

Numerical Simulation of Two-Part Underwater Towing System

The motion instability is an important issue that occurs during the operation of towed underwater vehicles (TUV), which considerably affects the accuracy of high precision acoustic instrumentations housed inside the same. Out of the various parameters responsible for this, the disturbances from the tow-ship are the most significant one. The present study focus on the motion dynamics of an underwater towing system with ship induced disturbances as the input. The study focus on an innovative system called two-part towing. The methodology involves numerical modeling of the tow system, which consists of modeling of the tow-cables and vehicles formulation. Previous study in this direction used a segmental approach for the modeling of the cable. Even though, the model was successful in predicting the heave response of the tow-body, instabilities were observed in the numerical solution. The present study devises a simple approach called lumped mass spring model (LMSM) for the cable formulation. In this work, the traditional LMSM has been modified in two ways. First, by implementing advanced time integration procedures and secondly, use of a modified beam model which uses only translational degrees of freedoms for solving beam equation. A number of time integration procedures, such as Euler, Houbolt, Newmark and HHT-α were implemented in the traditional LMSM and the strength and weakness of each scheme were numerically estimated. In most of the previous studies, hydrodynamic forces acting on the tow-system such as drag and lift etc. are approximated as analytical expression of velocities. This approach restricts these models to use simple cylindrical shaped towed bodies and may not be applicable modern tow systems which are diversed in shape and complexity. Hence, this particular study, hydrodynamic parameters such as drag and lift of the tow-system are estimated using CFD techniques. To achieve this, a RANS based CFD code has been developed. Further, a new convection interpolation scheme for CFD simulation, called BNCUS, which is blend of cell based and node based formulation, was proposed in the study and numerically tested. To account for the fact that simulation takes considerable time in solving fluid dynamic equations, a dedicated parallel computing setup has been developed. Two types of computational parallelisms are explored in the current study, viz; the model for shared memory processors and distributed memory processors. In the present study, shared memory model was used for structural dynamic analysis of towing system, distributed memory one was devised in solving fluid dynamic equations.

A review of dissolved oxygen modelling techniques for lowland rivers

This review introduces the methods used to simulate the processes affecting dissolved oxygen (DO) in lowland rivers. The important processes are described and this provides a modelling framework to describe those processes in the context of a mass-balance model. The process equations that are introduced all require (reaction) rate parameters and a variety of common procedures for identifying those parameters are reviewed. This is important because there is a wide range of estimation techniques for many of the parameters. These different techniques elicit different estimates of the parameter value and so there is the potential for a significant uncertainty in the model's inputs and therefore in the output too. Finally, the data requirements for modelling DO in lowland rivers are summarised on the basis of modelling the processes described in this review using a mass-balance model. This is reviewed with regard to what data are available and from where they might be obtained. (C) 2003 Elsevier Science B.V. All rights reserved.

Ice-sheet contributions to future sea-level change

Accurate simulation of ice-sheet surface mass balance requires higher spatial resolution than is afforded by typical atmosphere-ocean general circulation models (AOGCMs), owing, in particular, to the need to resolve the narrow and steep margins where the majority of precipitation and ablation occurs. We have developed a method for calculating mass-balance changes by combining ice-sheet average time-series from AOGCM projections for future centuries, both with information from high-resolution climate models run for short periods and with a 20 km ice-sheet mass-balance model. Antarctica contributes negatively to sea level on account of increased accumulation, while Greenland contributes positively because ablation increases more rapidly. The uncertainty in the results is about 20% for Antarctica and 35% for Greenland. Changes in ice-sheet topography and dynamics are not included, but we discuss their possible effects. For an annual- and area-average warming exceeding 4.5 +/- 0.9 K in Greenland and 3.1 +/- 0.8 K in the global average, the net surface mass balance of the Greenland ice sheet becomes negative, in which case it is likely that the ice sheet would eventually be eliminated, raising global-average sea level by 7 m.

Recent Glacier Retreat Exceeds Internal Variability

Glacier fluctuations exclusively due to internal variations in the climate system are simulated using downscaled integrations of the ECHAM4/OPYC coupled general circulation model (GCM). A process-based modeling approach using a mass balance model of intermediate complexity and a dynamic ice flow model considering simple shearing flow and sliding are applied. Multimillennia records of glacier length fluctuations for Nigardsbreen (Norway) and Rhonegletscher (Switzerland) are simulated using autoregressive processes determined by statistically downscaled GCM experiments. Return periods and probabilities of specific glacier length changes using GCM integrations excluding external forcings such as solar irradiation changes, volcanic, or anthropogenic effects are analyzed and compared to historical glacier length records. Preindustrial fluctuations of the glaciers as far as observed or reconstructed, including their advance during the “Little Ice Age,” can be explained by internal variability in the climate system as represented by a GCM. However, fluctuations comparable to the present-day glacier retreat exceed any variation simulated by the GCM control experiments and must be caused by external forcing, with anthropogenic forcing being a likely candidate.

Os efeitos sistêmicos de uma reestruturação produtiva

A presente dissertação tem como objetivo estudar o caso da empresa Lupatech CSL no que se refere, fundamentalmente, ao seu processo produtivo. Busca-se também discutir a possibilidade de se promoverem melhorias e mudanças que possam oferecer saídas que oportunizem, através de uma reestruturação produtiva, diferenciações competitivas ao negócio da empresa. Nesse sentido, considerando a análise sistêmica do ambiente, propõe-se a adoção de alguns princípios de gestão do M odelo Toyota de Produção, cujo um dos princípios preconiza uma relação de efetiva parceria com os fornecedores. Para tanto, considerando a permanente necessidade de as empresas buscarem seu equilíbrio e distinção no mercado, procurou-se fazer uma análise do pensamento sistêmico e outras fundamentações de convergência teórica com este trabalho. Além disso, a fim de melhor entender os prós e contras do pensamento toyotista de administração, buscou-se traçar uma investigação desde o antecessor modelo de produção em massa taylorista-fordista até os últimos acontecimentos que abalaram mundialmente a imagem da montadora Toyota. Acredita-se, ao final, que a análise dissertativa do caso da Lupatech poderá fortalecer e auxiliar a empresa em sua tomada de decisão, tanto no campo operacional como na definição de uma proposta estratégica de gestão.

Resistência e performance na arte de João Cota

Researching about the art of tell and the tales told by João Cota is somehow revisiting the oral tradition and social practice in the story-telling art. It takes into consideration the resistance this art still exerts, mainly by using the performance of oral transmission and receptiveness to tales. The study of this practice contributes to the vivacious and dynamcis permanence of this authentic and traditional storyteller of his repertory and of his form to tellin our culture. Story-telling is part of the humankind living heritage communicated by means of popular wisdom. Despite the risk of vanishing into thin air, along with their narrators, the tales still manage to resist the contemporary mass culture model. How long further will stories like the ones narrated by João Cota be able to resist to strong and stronger structures dictated by writing and other communication means? João Cota s practice in story-telling will be studied not only as a proposal to identify the presence of this practice and the oral cultural resistance but also, through the performance prospective, to identify the oral transmission and receptiveness to the tales that are part of this storyteller s repertoire. In other words: what he tells, how he tells it, and why he tells it. The advent of new technologies such as internet, through which people can easily communicate with others in different parts of the world, and the greater and greater expansion in the writing skill concept interfere the maintenance of the oral tradition elements present in João Cota s narratives and inserted in the Brazilian culture. This has become more visible in the latest decades although we still notice the living tradition and permanence of the story-telling practice in several parts of the country through their wise storytellers. Our research target will require - in each of its study stages reference to works by several theoreticians namely Paul Zumthor, Mikhail Bakhtin, Câmara Cascudo, theoreticians from the receptiveness aesthetic, from the written culture history, from oral cultures and reading practices, from tradition and the Brazilian Culture of oral story-telling. In order to get to know and draw a profile of this storyteller, we ve chosen to use the comprehensive interview method by French Sociologist Jean Claude Kaufmann. The originality in the method consists of the qualitative data put together in situ , concentrated on the storyteller s narratives/speeches recorded on tape, which will be the focal point of this study. Our analysis method is based on tireless sessions of listening to interviews out of which we gathered information related to the storyteller, his practice in telling the tales, and his repertoire

Energy and angular momentum of the gravitational field in the teleparallel geometry

The Hamiltonian formulation of the teleparallel equivalent of general relativity is considered. Definitions of energy, momentum and angular momentum of the gravitational field arise from the integral form of the constraint equations of the theory. In particular, the gravitational energy-momentum is given by the integral of scalar densities over a three-dimensional spacelike hypersurface. The definition for the gravitational energy is investigated in the context of the Kerr black hole. In the evaluation of the energy contained within the external event horizon of the Kerr black hole, we obtain a value strikingly close to the irreducible mass of the latter. The gravitational angular momentum is evaluated for the gravitational field of a thin, slowly rotating mass shell. © 2002 The American Physical Society.

Estimativa da radiação solar ultravioleta em Botucatu/SP Brasil utilizando técnicas de aprendizado de máquina

Pós-graduação em Agronomia (Energia na Agricultura) - FCA

A influência do Arco do Desmatamento sobre o ciclo hidrológico da Amazônia

Este trabalho buscou quantificar algumas alterações no ciclo hidrológico da Amazônia propiciadas pelo desmatamento da região, principalmente devido à faixa conhecida como "Arco do Desmatamento". Neste sentido, foram realizados experimentos numéricos utilizando o modelo BRAMS (Brazilian Regional Atmospheric Modeling System) tendo o submodelo de vegetação dinâmica GEMTM (General Energy and Mass Transport Model) a ele acoplado. Foram investigados os impactos causados pelo Arco do Desmatamento atual em relação à floresta intacta, bem como as futuras modificações, causados pelo avanço do desmatamento até o ano de 2050. Como condições de contorno na superfície para o modelo BRAMS, foram usados cenários oriundos de modelos empíricos de desmatamento para os anos de 2002 e 2050. Os resultados mostraram que o avanço do Arco do Desmatamento até 2050 tem uma complexa relação com as variáveis analisadas. Por exemplo, a precipitação apresentou distribuição espacial heterogênea, com áreas de anomalias positivas e negativas que se mostraram coerentes com as anomalias de outras variáveis, como a evapotranspiração e a divergência de umidade. Também foram encontradas algumas áreas que evidenciaram as possíveis influências dos grandes rios e topografia da região sobre essa precipitação. Os balanços de radiação e energia também foram afetados pelo desmatamento, sendo que grande parte das alterações é devido à mudança do albedo da superfície, a qual ocorre com a substituição da floresta pela pastagem. Quanto às alterações propiciadas pelo Arco do Desmatamento atual em relação à floresta intacta, foi observado que todas as variáveis foram afetadas, entretanto a maioria dos impactos ainda é localizada sobre a área mais afetada pelo desmatamento, diferentemente dos resultados encontrados para o cenário de 2050, onde as modificações já se dão a nível regional.

Avaliação da dinâmica temporal da evaporação no reservatório de Itumbiara, GO, utilizando dados obtidos por sensoriamento remoto

The object of this work was to study the dynamics of evaporation in the Itumbiara reservoir, located in Central Brazil, using MODIS-derived water surface temperature (product MOD11A1) and meteorological data acquired over the water surface. The evaporation rates were derived from latent heat flux, estimated through a mass transfer model. The estimates were carried out for the period between 1/1/2010 and 31/12/2010. The results showed that evaporation rate tends to increase from January to September and then decrease from September to December. The evaporation rate reached values near 20 mm day-1 in Itumbiara reservoir during the dry season in 2010. The mean evaporation rate for the wet season was 3.66 mm day-1 and 8.25 mm day-1 for the dry season. The total water volume evaporated from Itumbiara reservoir during 2010 was estimated at about 1.7 billion m³ (2,300 mm) which represents 10% of total reservoir volume. The results suggest that advection is the main transport mechanism which drives the evaporation in Itumbiara. The convective processes contribute secondarily to evaporation in Itumbiara reservoir.

Um estudo da recente transformação na estrutura social do Brasil: nova classe média?

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)