During this study two models were developed to predict growth of Pinus radiata D.Don plantations in Canterbury, New Zealand. The first, CanSPBL(1.2), is a model for whole rotations of stands owned by Selwyn Plantation Limited in Canterbury. The second model, CanSPBL(water) is a hybrid growth model for the Selwyn estate in Canterbury that incorporates an index of root zone water balance over the simulation period. An existing stand growth and yield model CanSPBL was examined using a validation dataset of PSP measurements that were not used in model fitting. Projection bias was shown for mean top height, basal area per hectare, and residual stand stocking particularly for stands at elevations exceeding 450 metres. The new model, CanSPBL(1.2) showed an increase in precision of 4 - 46% over CanSPBL(1.0) at a stand level. The components of the stand model include mean top height, basal area per hectare, stems per hectare, and diameter distribution. The mortality model was made in conjunction with managers at CanSPBL to exclude catastrophic mortality events from model projections. Data used for model fitting was filtered using a mortality index based on the -3/2 power law. An examination of this model with an independent dataset showed little apparent bias. The new model, CanSPBL(water) was developed to include an index of water balance over the simulation period. Water balance estimates were made using a sub model for root zone water balance included in the hybrid physiological model 3-PG (Landsberg and Waring, 1997). The new model showed an increase in precision of 1 - 4% over CanSPBL(1.2) at a stand level (with the exception of the model for maximum diameter which showed a decrease in precision of 0.78%) using climatic inputs that included yearly variation. However the model showed increases of precision from 0.5 to 8% (with the exception of maximum diameter again, showing a decrease in precision of 0.13%) using long term monthly average climatic inputs. The components of the stand model also include mean top height, basal area per hectare, stems per hectare, and diameter distribution. The mortality model was also fitted with a data set filtered using a mortality severity index based on the -3/2 power law to exclude catastrophic mortality events. An examination of this model with an independent dataset showed little apparent bias. Two models to predict a one sided canopy leaf area index (LAI) of radiata pine stands in the Canterbury Plains of New Zealand were also developed. The models were fitted using non-linear least squares regression of LAI estimates against stem measurements and stand characteristics. LAI estimates were derived from digital analysis of fisheye lens photography. The models were kept simple to avoid computational circularity for physiological modelling applications. This study included an objective comparison and validation of a range of model types. The models CANTY (Goulding, 1995), CanSPBL(1.2) (Pinjuv, 2005), CanSPBL-water (Pinjuv, 2005), and 3-PG (Landsberg and Waring, 1997) were compared and validated with the main criteria for comparison being each model s ability to match actual historical measurements of forest growth in an independent data set. Overall, the models CanSPBL(water), and CanSPBL(1.2) performed the best in terms of basal area and mean top height prediction. Both models CanSPBL(water), and CanSPBL(1.2) showed a slightly worse fit in predictions of stocking than did the model CANTY. The hybrid model 3PG showed a better fit for the prediction of basal area than the statistically based model CANTY, but showed a worse fit for the prediction of final stocking than all other models. In terms of distribution of residuals, CanSPBL(1.2) had overall the lowest skewness, kurtosis, and all model parameters tested significant for normality. 3PG performed the worst on average, in terms of the distribution of residuals, and all models tested positively for the normality of residual distribution.

La croissance du marché des pompes à chaleur contribue à l'augmentation de la thermosensibilité de la demande électrique. Il devient nécessaire de mieux comprendre l'impact des usages thermosensibles de l'électricité, notamment concernant ceux qui sont corrélés de manière non linéaire à la température extérieure. Dans cette optique, cette thèse vise à construire un cadre de modélisation qui permette i) d'analyser les facteurs d'influence de la thermosensibilité à partir d'une description physique des usages thermosensibles, et ii) de réaliser des diagnostics de ces paramètres d'influence tout en tenant compte des incertitudes associées. Une approche de modélisation hybride qui bénéficie des avantages de modèles statistiques et de modèles physiques est principalement employée pour répondre à ces questions.La première étape consiste à estimer la part thermosensible de la demande réelle par un modèle prédictif top-down. On développe ensuite un modèle d'analyse physique de la thermosensibilité à l'échelle régionale à partir de la thermique du bâtiment. On s'appuie notamment sur des modèles pseudo-physiques de performance de pompes à chaleur qui sont régressés sur des données constructeur ou des mesures de performances réelles. Un COP régional est déterminé pour l'ensemble des PAC installées. Enfin, les paramètres d'influence du modèle de thermosensibilité ainsi développé sont estimés à l'aide de l'approche bayésienne, qui offre un cadre pour le traitement de l'incertitude sous la forme de probabilités. Des coefficients équivalents de déperditions thermiques, une température intérieure équivalente ainsi que les parts du chauffage Joule et par PAC pour le parc de bâtiments régional ont été obtenus. / The growing heat pump market contributes to the increase in temperature sensitivity of electricity demand. It becomes necessary to understand the impact of temperature sensitive end-uses of electricity, including those which are correlated non-linearly to the outside temperature. In this context, this thesis aims to build a modeling framework to i) analyze the influencing factors of the temperature sensitivity of electricity demand from a physical description of temperature-sensitive equipment, and ii) to perform diagnoses of these parameters of influence by taking into account the associated uncertainties. A hybrid modeling approach that benefits the advantages of statistical models and physical models is used to answer these questions.Firstly, the temperature-sensitive part of electricity demand is estimated by a predictive top-down model. Then a physical model to analyze the temperature sensitivity at regional level is developed based on building thermal energy needs. A regional coefficient of performance (COP) is determined for the whole installed heat pumps by using pseudo-physical models which are regressed on manufacturer data or actual performance measures. Finally, the parameters of influence of the developed temperature sensitivity model are estimated using the Bayesian approach which provides a framework for the treatment of uncertainty in the form of probabilities. Equivalent coefficients of heat loss, an equivalent internal temperature, as well as the share of Joule heating and the share of heat pumps for the regional building stock are obtained.

Malgré l'augmentation constante des ressources informatiques dédiées au calcul scientifique, simuler des écoulements réactifs mettant en jeu une chimie complexe reste aujourd'hui encore un véritable challenge. L'objectif de cette thèse est le développement de la méthode Hybrid Transported-Tabulated Chemistry (HTTC), destinée aux simulations DNS/LES de flammes avec des mécanismes cinétiques détaillés, en offrant un temps de calcul acceptable. Cette nouvelle approche consiste à transporter les espèces majoritaires de l'écoulement, tandis que les espèces minoritaires sont extraites d'une table chimique. La méthode HTTC a été implémentée dans un code DNS/LES et validée sur des flammes 1D de méthane et de kérosène, mettant en évidence une réduction extrêmement importante du temps de calcul, comparé aux solveurs classiques de chimie détaillée. HTTC a ensuite été mis en œuvre avec succès sur des flammes triples de méthane en présence de forts gradients de fraction de mélange. L'impact des méthodes choisies pour prolonger la table chimique et pour calculer les variables de contrôle, utilisées pour paramétrer la table, a été étudiée avec une attention particulière. Un très bon accord a été trouvé avec les résultats de référence, obtenus avec un solveur de chimie détaillée. / Even if significant progress is being made to improve the power of high-performance computers, the numerical simulation of reactive flows involving complex chemistry is still a challenging task. The objective of this work is the development of the Hybrid Transported-Tabulated Chemistry method (HTTC), designed for the DNS/LES simulations of flames with detailed kinetic mechanisms, with an acceptable cost. This novel approach combines the transport of the main species in the flow with the tabulation of the radical species. It has been implemented in a DNS/LES code and validated on 1D methane and kerosene flames. The cost of the simulations has been considerably decreased, compared to classic detailed chemistry solvers. Then, simulations of methane edge flames, featuring large gradients of mixture fraction, have been performed with HTTC. In particular, the impact of the methods used to extend the chemical tables and to compute the control variables have been analyzed in details. A very good agreement has been found by comparison with detailed chemistry.

Extensive use of the Amudarya river waters for irrigation has severely impacted semi-natural ecosystems along its course and in its delta region. Currently, new strategies are searched for multi-user and multi-objective water management to mitigate ecological and socio-economic deterioration. The GIS-based modeling framework, TUGAI, has been developed to support exploration of alternative water management strategies for the Amudarya river delta and to analyze their ecological implications. Available information of heterogeneous type and quality on resource availability and habitat demands of deltaic ecosystems has been integrated into a comprehensive tool by a hybrid approach. A multi-objective water allocation model, AmuEPIC, has been combined with simple, spatially-explicit statistical and rule-based models of landscape dynamics, AmuGIS, and an ecological assessment procedure based on a fuzzy habitat suitability index model for riverine Tugai forests, TugaiHSI. Users can develop scenarios of alternative water management strategies for a time period of up to 30 years and compare their ecological effects. The tool facilitates a first quick assessment of the response of the delta environment to water management measures in a problem-oriented way. It assists in structuring the problem of water allocation to the environment, facilitates analysis of tradeoffs and uncertainties, fosters discussion between stakeholders and supports a goal finding process. Results of scenario analysis demonstrate solutions to given management tasks, which can serve as goals for implementation of measures in reality. First testing results indicate that there is a potential for increase of water discharge for environmental needs, while, at the same time, providing irrigation and other water users with sufficient water.

Résumé : Les véhicules aériens et terrestres sont constitués de systèmes bâtis complexes. La structure principale est généralement composée de panneaux légers renforcés par des éléments rigides. Cette solution de conception est répandue parce qu’elle allie la force et un faible poids. Cependant, on sait qu’elle offre des résultats vibroacoustiques médiocres, c’est à dire que l’effet des perturbations externes qui touchent le système peut générer un niveau de bruit excessif à l’intérieur de la cabine des passagers. C’est une préoccupation majeure chez les fabricants, parce que ce niveau de bruit nuit sensiblement au confort ressenti par les clients et peut causer de la fatigue chez les conducteurs et les pilotes. Pour cette raison, les composants passifs constitués de matériaux dissipatifs assemblés en mode multicouche sont généralement intégrés à la structure. Ces assemblées bordées intègrent surtout des matériaux poroélastiques, qui sont plutôt répandus, grâce à l’agencement intéressant de bonnes propriétés d’isolation sonore et de faible poids. L’intégration en amont des traitements de contrôle du bruit au processus de conception est la clé de succès d’un produit. Pour ce faire, des outils pratiques numériques en mesure de capter le comportement dynamique des systèmes vibroacoustiques impliquant les structures, les cavités et les matériaux d’insonorisation sont requis. D’une part, la modélisation de ces systèmes couplés en utilisant des procédés à base d’éléments finis peut être, bien que précis, irréalisable pour des applications pratiques. D’autre part, les approches analytiques telles que la méthode de matrice de transfert sont souvent préférées grâce à leur facilité d’utilisation, même si elles manquent de précision en raison des hypothèses rigoureuses inhérentes au cadre analytique. Dans ce contexte, les procédures de structuration hybrides sont récemment devenues très populaires. En effet, les différentes techniques de modélisation sont généralement recherchées pour décrire les systèmes vibroacoustiques complexes arbitraires sur la plus large gamme de fréquences possible. L’objectif du projet proposé est de mettre au point un cadre hybride offrant une mé- thodologie simple pour tenir compte des traitements de contrôle du bruit dans l’analyse vibroacoustique par éléments finis. A savoir, le modèle de calcul qui en découle conserve la souplesse et la précision de la méthode des éléments finis en bénéficiant de la simplicité et de l’efficacité de la méthode de matrice de transfert pour obtenir une réduction de la charge de calcul pour la modélisation de composants acoustiques passifs. La performance de la méthode pour prédire la réponse vibroacoustique de structures planes homogènes avec des traitements acoustiques attachées est évaluée. Les résultats démontrent que la méthode hybride proposée est très prometteuse, parce qu’elle permet une réduction de l’effort de calcul tout en conservant suffisamment de précision par rapport à l’analyse complète par éléments finis. En outre, la méthode de matrice de transfert proposée de modélisation des traitements de contrôle des bruits est générale, comme on peut l’appliquer dans d’autres cadres outre l’application de l’élément fini considéré dans ce travail. // Abstract : Aerial and terrestrial vehicles consist of complex built-up systems. The main structure is typically made of light panels strengthened by stiffer components. Such design solution is widely used as it combines strength and low weight. However, it is known to give poor vibroacoustic performances, i.e. the effect of external disturbances acting on the system may generate an excessive noise level inside the passengers cabin. This is a main concern for the manufacturers, as it significantly affects the comfort experienced by the costumers and may fatigue drivers and pilots. For this reason, passive components consisting of dissipative materials assembled in a multilayer fashion are typically integrated within the structure. These lined assemblies mainly involve poroelastic materials, which are commonly used thanks to the appealing combination of good noise insulation properties and low weight. The early integration of noise control treatments in the design process is the key to a successful product. For this purpose, practical numerical tools able to capture the dynamic behavior of vibroacoustic systems involving structures, cavities and noise proofing materials are demanded. On the one hand, modeling such coupled systems using finite element based methods can be, albeit accurate, time consuming for practical applications. On the other hand, analytical approaches such as the transfer matrix method are often preferred thanks to their ease of use, although they suffer from a lack of accuracy due to the stringent assumptions inherent within the analytical framework. In this context, hybrid substructuring procedures have recently become quite popular. Indeed, different modeling techniques are typically sought to describe arbitrarily complex vibroacoustic systems over the widest possible frequency range. The aim of this thesis is to devise a hybrid framework providing a simple methodology to account for noise control treatments in vibroacoustic finite element analysis. Namely, the resulting computational model retains the flexibility and accuracy of the finite element method while taking advantage from the simplicity and efficiency of the transfer matrix method to obtain a reduction of the computational burden in the modeling of passive acoustic components. The performance of the method in predicting the vibroacoustic response of flat structures with attached homogeneous acoustic treatments is assessed. The results prove that the proposed hybrid methodology is very promising, as it allows for a reduction of the computational effort while preserving enough accuracy with respect to full finite element analysis. Furthermore, the proposed transfer matrix based methodology for noise control treatments modeling is general, as it can be used in alternative frameworks besides the finite element application considered in this work.

Les SPA sont des moteurs piézoélectriques inertiels amplifiés. Ils fonctionnent en convertissant des vibrations asymétriques en μ-déplacements par frottement, grâce au phénomène d’adhérence-glissement. Ces vibrations et ce frottement soulèvent des problématiques liées à l’usure et au bruit. Pour y répondre, ces travaux étudient les comportements tribologique et vibratoire de ces moteurs. Dans un premier temps, un tribomètre pion-plan, actionné par un SPA, a été développé. Il a permis d’observer, in-situ et en fonctionnement, les débits solides interfaciaux: le 3ème corps. Ces observations directes, couplées à des analyses post-mortem ont permis de décrypter les sollicitations tribologiques et de multiplier par dix la durée de vie des moteurs. Dans un second temps, des analyses numérique et expérimentale ont été proposées afin de mieux appréhender le comportement vibratoire des SPA. Le modèle existant, à constantes localisées, a été modifié en une version hybride, constantes localisées-FEM, plus versatile. Les tests expérimentaux ont permis de déterminer le principal contributeur acoustique et de réduire le niveau de bruit de seize dBA. Finalement, la double approche tribologie/mécanique vibratoire a permis de mieux saisir les subtilités des SPA. Elle a notamment mené au développement de structures innovantes tels qu’un moteur à trois degrés de liberté et deux moteurs rotatifs. / SPA are amplified piezoelectric inertia motors. They operate by converting asymmetrical vibrations into μ-displacements by means of friction through stick-slip. These vibrations and this friction raise issues related to wear and noise. So, the present work investigates both the tribological and vibratory behaviours of these motors. First, a pin-on-pad tribometer, actuated by a SPA, was developed. It led us to observe, in-situ and in operation, interfacial solid flows: the 3rd body. These direct observations carried out in conjunction with post-mortem analyses resulted in decrypting the tribological solicitations and in increasing by ten the lifetime of the motors. Second, numerical and experimental considerations have been proposed to better grasp the vibratory behaviour of SPA. The existing lumped model was modified to a more versatile hybrid lumped-FEM model. The experimental tests made it possible to determine the main acoustic contributor and to reduce the noise level by sixteen dBA. Finally, this dual - tribology/vibratory - approach provided a better understanding of the SPA intricacies. In particular, it led to develop innovative structures such as a motor with three degrees of freedom and two rotary motors.

Flow and transport phenomena in the subsurface often span a wide range of length (nanometers to kilometers) and time (nanoseconds to years) scales, and frequently arise in applications of CO₂ sequestration, pollutant transport, and near-well acid stimulation. Reliable field-scale predictions depend on our predictive capacity at each individual scale as well as our ability to accurately propagate information across scales. Pore-scale modeling (coupled with experiments) has assumed an important role in improving our fundamental understanding at the small scale, and is frequently used to inform/guide modeling efforts at larger scales. Among the various methods, there often exists a trade-off between computational efficiency/simplicity and accuracy. While high-resolution methods are very accurate, they are computationally limited to relatively small domains. Since macroscopic properties of a porous medium are statistically representative only when sample sizes are sufficiently large, simple and efficient pore-scale methods are more attractive. In this work, two Eulerian pore-network models for simulating single-phase flow and solute transport are developed. The models focus on capturing two key pore-level mechanisms: a) partial mixing within pores (large void volumes), and b) shear dispersion within throats (narrow constrictions connecting the pores), which are shown to have a substantial impact on transverse and longitudinal dispersion coefficients at the macro scale. The models are verified with high-resolution pore-scale methods and validated against micromodel experiments as well as experimental data from the literature. Studies regarding the significance of different pore-level mixing assumptions (perfect mixing vs. partial mixing) in disordered media, as well as the predictive capacity of network modeling as a whole for ordered media are conducted. A mortar domain decomposition framework is additionally developed, under which efficient and accurate simulations on even larger and highly heterogeneous pore-scale domains are feasible. The mortar methods are verified and parallel scalability is demonstrated. It is shown that they can be used as “hybrid” methods for coupling localized pore-scale inclusions to a surrounding continuum (when insufficient scale separation exists). The framework further permits multi-model simulations within the same computational domain. An application of the methods studying “emergent” behavior during calcite precipitation in the context of geologic CO₂ sequestration is provided. / text

Ce mémoire de thèse présente deux nouvelles approches pour l’analyse et la commande des systèmes non-linéaires complexes, comme les systèmes dynamiques à commutation de la classe des convertisseurs d’énergie électrique. Ces systèmes ont plusieurs modes de fonctionnement et ont un point de fonctionnement désiré qui, en général, n’est le point d’équilibre d’aucun des modes. Dans cette classe de systèmes, la commutation d’un mode de fonctionnement à un autre est commandée selon une loi qui doit être synthétisée. Par conséquent, la synthèse de commande implique l’étude des conditions qui permettent à un cycle limite stable de s’établir au voisinage du point de fonctionnement désiré, puis de la trajectoire de commande qui permet de l’atteindre en respectant les contraintes physiques de comportement (courant maximum supporté par les composants,. . .) ou les contraintes de temps (durée minimum entre deux commutations,. . .). Le cycle limite sera qualifié d’hybride car il est composé de plusieurs dynamiques(deux dans ces travaux).La première méthode développée s’appuie sur les propriétés géométriques des champs de vecteurs et est une extension d’une partie des travaux de thèse de Manon au LAGEP. Une condition nécessaire et suffisante d’existence et de stabilité d’un cycle limite hybride composé d’une séquence de deux modes de fonctionnement dans IR2 est présentée. Ce cycle définit la région finale à atteindre par le système depuis son état initial, par une trajectoire déterminée de manière optimale selon un critère donné (durée totale, énergie dépensée, . . .). La méthode proposée est appliquée aux convertisseurs d’énergie Buck et Buck-Boost alimentant une charge résistive. Une extension à IRn a été proposée et démontrée. Elle est illustrée sur un système non-linéaire dans IR3.La deuxième méthode est développée dans IR2 et basée sur la théorie de Lyapunov, bien connue en automatique pour étudier la stabilité des systèmes non-linéaires et concevoir des commandes stabilisantes.Il s’agit de déterminer par une approche géométrique, une fonction de Lyapunov quadratique commune aux deux modes de fonctionnement du système, qui permette d’obtenir un cycle limite hybride stable le plus proche possible du point de fonctionnement désiré et une commande stabilisante directe des interrupteurs / This PhD-thesis presents two new approaches for the analysis and control of complex nonlinearsystems, such as switching dynamic systems of the class of power converters. These systems have severalmodes of operation and a desired operating point, which, in general, is not the equilibrium point of anymode. In this class of systems, switching from one mode to another is controlled by a switching law tobe designed. Therefore, the synthesis of a control law involves the study of the conditions that allow astable limit cycle to settle near the desired operating point, and of the control trajectory to reach thislimit cycle and stabilize on it, meeting the constraints dues to the physical behavior (maximum currentsupported by the components, . . .) or time constraints (minimum duration between two switchings, . . .).The limit cycle is called hybrid because it is composed of several dynamics (two in this work).The first method is based on the geometric properties of vector fields and is an extension of part ofthe PhD-thesis of Manon at LAGEP. A necessary and sufficient condition of existence and stability of ahybrid limit cycle consisting of a sequence of two operating modes, is presented in IR2. This limit cycledefines the final region to be reached by the system from its initial state, along a trajectory determinedoptimally according to a given criterion (total duration, energy expended, . . .). This method is applied tothe Buck and Buck-Boost power converters with a resistive load. An extension to IRn has been proposedand demonstrated. It is illustrated on a nonlinear system in IR3.The second method is developed in IR2, based on the Lyapunov theory, well known in automatic controlfor studying the stability of nonlinear systems and designing stabilizing control methods. The objective isto design, with a geometric approach, a quadratic Lyapunov function common to both modes of operation,which defines a stable hybrid limit cycle closest to the desired operating point and a direct stabilizingcontrol of the switches.

Les nouveaux services de télécommunication accessibles en mobilité entraînent une demande croissante de débit et amènent aux limites de capacités des réseaux actuels. Pour répondre à cette demande, les opérateurs cellulaires déploient de nouvelles technologies et densifient leurs réseaux en ajoutant des petites cellules, où les stations de base sont plus proches des utilisateurs. Compte tenu de la visibilité directe et de la faible hauteur des antennes, le canal de propagation des petites cellules est différent de celui des macro-cellules et les modèles stochastiques de canal existants ne sont plus adaptés. C'est dans ce contexte que se situe cette thèse, qui explore les problématiques de dynamique dans le canal de propagation et dans l'utilisation des ressources radio. Nous proposons dans une première partie de caractériser les dynamiques dans le canal de propagation de petite cellule en milieu extérieur. Cette caractérisation se base sur des prédictions déterministes, où une partie de l'environnement est généré de manière aléatoire afin de procéder à des études statistiques. Pour cela, nous mettons en place une plateforme de simulation, capable de générer des objets statiques, de les déplacer et de les prendre en compte dans des outils de prédiction de propagation adaptés aux petites cellules. Une phase de mesures de l'impact du trafic de véhicules sur la puissance du signal permet de valider les outils de propagation développés dans la plateforme. L'objectif de la deuxième partie de la thèse est d'évaluer l'impact de ces nouvelles topologies sur les performances des réseaux hétérogènes. Nous présentons un modèle hybride de propagation qui prend en compte les variations de signal dues au trafic de véhicules en combinant les approches déterministes et statistiques, pour avoir les avantages de résultats précis et d'un temps de calcul rapide. Ce modèle hybride est utilisé par la suite dans un processus de simulation de réseau cellulaire 3G LTE au niveau système, où nous évaluons l'impact du trafic de véhicules sur le débit maximal, la consommation d'énergie et le respect de la qualité de service.

Made available in DSpace on 2017-07-10T18:08:16Z (GMT). No. of bitstreams: 1 Thiago DArisbo.pdf: 2108504 bytes, checksum: 7b8aad29ec7d75a6fd370e54a95cd849 (MD5) Previous issue date: 2011-02-24 / Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / Ion exchange is a process that is used in the treatment of aqueous industrial effluents containing organic compounds and heavy metals. The fixed bed columns are longer applied by allowing the process to occur continuously (cycles of regeneration). The design and process optimization of the ion exchange column requires the use of mathematical models. Phenomenological models of these systems involve the solution of partial differential and algebraic equations. The equilibrium data for ion exchange processes are usually described by the Mass Action Law (MAL), which can be considered non-ideality of aqueous and solid phases. Artificial Neural Networks (ANN) are being used successfully for the study of equilibrium data because they are empirical models and don t demand a mathematical rigor. This work aimed to evaluate the applicability of the hybrid model to describe the dynamics of ion exchange in fixed beds of binary systems. This system consists of partial differential equations obtained from mass balance in fluid phases in the ion exchanger and ANN to describe the balance. LAM was adjusted to experimental data of ion exchange equilibrium and then were generated 4200 data sets for each binary pair studied, which served as training for RNA. We tested networks with different structures, with one and two input layers. The 3-3-2 structure was used in the simulations of the hybrid model because it was the best represented the systems during the training phase. The differential equations were solved by the lines method. A computer program in FORTRAN language was developed for solving the model equations. DASSL subroutine was used to solve the equations. The performance of the hybrid model was evaluated from the results obtained with the phenomenological model, in which case the equilibrium description was made with the use of MAL. It also was the analysis of results from the comparison of experimental data. To evaluate the model we used data from the literature of ion exchange in Amberlite IR 120 resin on the systems Cu-Na and Zn-Na and in NaY zeolite on Fe-Na and Zn-Na. Both models were efficient to describe the dynamics of ion-exchange fixed bed columns, and the hybrid model had the advantage of the reduced computational time (82% reduction on average) as a result of not needing to solve a nonlinear equation. / A troca iônica é um processo muito utilizado no tratamento de efluentes industriais aquosos contendo compostos orgânicos e metais pesados. As colunas de leito fixo são mais aplicadas por permitir que o processo ocorra de maneira contínua (ciclos de regeneração). O projeto e a otimização de processos de troca iônica em coluna requer o uso de modelos matemáticos. Os modelos fenomenológicos destes sistemas envolvem a resolução de equações diferenciais parciais e algébricas. Os dados de equilíbrio de processos de troca iônica geralmente são descritos pela Lei da Ação das Massas (LAM), na qual podem ser consideradas as não idealidades das fases aquosa e sólida. As Redes Neurais Artificiais (RNA) estão sendo utilizadas com sucesso para o estudo destes dados de equilíbrio por serem modelos empíricos e não demandarem tal rigor matemático. Esta dissertação teve por objetivo avaliar a aplicabilidade do modelo híbrido para descrever a dinâmica do processo de troca iônica em leito fixo de sistemas binários. Este sistema é constituído de equações diferenciais parciais obtidas por meio de balanço de massa nas fases fluida e no trocador iônico e de RNA para descrever o equilíbrio. A LAM foi ajustada a dados experimentais de equilíbrio de troca iônica e, então, foram gerados conjuntos de 4200 dados para cada par binário estudado, os quais serviram como treinamento para a RNA. Foram testadas redes com diferentes estruturas, com uma e com duas camadas de entrada. A estrutura 3-3-2 foi utilizada nas simulações do modelo híbrido, pois foi a que melhor representou os sistemas na etapa de treinamento. As equações diferenciais foram resolvidas pelo método das linhas. Um programa computacional em linguagem FORTRAN foi desenvolvido para a resolução das equações do modelo. Foi utilizada a sub-rotina DASSL para resolver as equações. O desempenho do modelo híbrido foi avaliada a partir dos resultados obtidos com o modelo fenomenológico, sendo que neste caso a descrição do equilíbrio foi feita pelo uso da LAM. Também foi feita a análise dos resultados a partir da comparação dos dados experimentais. Para avaliar o modelo foram utilizados dados da literatura de troca iônica em resina Amberlite IR 120 dos sistemas Cu-Na e Zn-Na e na zeólita NaY dos sistemas Fe-Na e Zn-Na. Ambos os modelos foram eficientes para descrever a dinâmica de troca iônica de colunas de leito fixo, sendo que o modelo híbrido apresentou como vantagem o menor tempo computacional (82% de redução em média) em decorrência de não necessitar resolver a equação não-linear.