Ce doctorat trouve son origine dans la compréhension des prises de décision et des comportements collectifs des animaux. Comment ces derniers parviennent-ils à effectuer des choix collectivement ? Comment les membres d’un groupe procèdent-ils pour synchroniser leurs comportements spatialement et temporellement ? Mon principal objectif a été de dégager, lors des déplacements collectifs et du fur rubbing chez le capucin moine, les évènements décisionnels dépendants de processus anonymes de ceux dépendants de processus liés à l’identité des individus et à leur réseau de relations sociales au sein du groupe. Dans les prises de décision collective relatives aux déplacements, les membres du groupe sont influencés dans leurs choix par leur identité sociale mais aussi par des mécanismes anonymes, de type mimétique. Le fur rubbing est également un comportement collectif dont les mécanismes sous-jacents incluent une dépendance interindividuelle de type mimétique. Des mécanismes similaires mettant en jeu des interactions entre individus basées sur des règles comportementales simples se retrouvent dans chacun des phénomènes collectifs étudiés. Ces résultats sont les premiers à démontrer l’émergence de prises de décision collective à partir de telles interactions anonymes dans un groupe de primates non humains. Ils permettent de faire le lien entre choix individuels et comportement collectif et de mieux concevoir comment un groupe de primates peut se coordonner, maintenir sa cohésion spatiale et synchroniser ses activités./How do animals reach collective consensus? How do group members spatially and temporally synchronise their behaviour? My main purpose was to demonstrate the respective roles of anonymous processes (contagion, mimetism) and individual-dependent processes (hierarchical rank, age, sex, kin, social relationships) in collective decision-making. During decision-making relating to collective movements, group members’ decisions depend on their social identity (individual-dependent mechanism) as well as anonymous processes. Fur rubbing is also a collective behaviour involving interindividual dependence with mimetic underlying mechanisms. We found similar mechanisms, involving interindividual interactions according to simple behavioural rules, in both collective phenomenon studied. These results are the first to demonstrate the emergence of collective decision-making based on anonymous interactions in a group of non human primates. They help to understand the link between individual choices and collective behaviour and to appreciate how a social group of primates maintain its spatial cohesion and synchronize its activities.

Dans une optique de maîtrise du procédé d'incinération des ordures ménagères et de ses possibles émissions polluantes, nous avons développé un modèle mathématique qui simule un lit d'ordures ménagères en combustion sur une grille mobile. Ce modèle décrit la plupart des phénomènes physicochimiques et thermiques intervenant lors de l'incinération : séchage et pyrolyse de la charge, combustion et gazéification du carbone résiduel, transferts thermiques, effondrement du lit, brassage... Il intègre également une description des mécanismes de volatilisation des métaux lourds et de formation des NOx. La cinétique de départ des métaux lourds est modélisée en tenant compte des différentes étapes de transport (transfert externe, diffusion intraparticulaire, volatilisation) au moyen de l'approche des temps caractéristiques additifs. Dans le cas simulé du cadmium, la prédiction d'une volatilisation quasi-complète est conforme aux résultats de la littérature. Le sous-modèle NOx prend en compte les mécanismes de formation thermique, prompt, combustible, par l'intermédiaire de N2O, ainsi que les mécanismes de réduction homogène par recombustion et hétérogène par le carbone résiduel. Les calculs révèlent que prédominent la formation par le mécanisme combustible et la destruction par la réduction hétérogène. Enfin, le modèle de lit a été utilisé pour étudier l'influence des divers paramètres opératoires : température, débit et distribution d'air primaire, taille des particules de déchets, brassage et schéma de brassage. Les résultats sont présentés et discutés en détail. L'influence des conditions opératoires sur l'efficacité du procédé et sur les émissions de Cd et NOx est analysée.

La majeure partie des travaux sur les problèmes d’ordonnancement se placent dans le contexte où les ressources sont disponibles en permanence. Ce qui en réalité n’est pas toujours le cas. Nous nous plaçons dans le contexte d’indisponibilités connues ; nous nous intéressons plus particulièrement aux problèmes de type job shop avec des périodes d’indisponibilité flexibles et des tâches pouvant éventuellement être interrompues par les périodes d’indisponibilité. L’intégration de ces contraintes rend les problèmes d’ordonnancement nettement plus difficiles à résoudre. La flexibilité que nous considérons peut être relative à au moins l’un des points suivants : déplacement de la période d’indisponibilité dans une fenêtre de temps, modification de la durée de la période d’indisponibilité, interruption d’une tâche par une période d’indisponibilité, ensuite reprise avec une éventuelle pénalité.Dans cette thèse, nous avons proposé des modèles mathématiques pour le problème. En plus de la résolution des problèmes considérés, le but de ces modélisations est de permettre d’analyser l'impact des différentes contraintes et d'évaluer la qualité des méthodes approchées que nous proposons. Ces dernières permettent de construire très rapidement un ordonnancement en se basant sur des règles de priorité. Les solutions sont aussi utilisées pour notre approche basée sur la génération de colonnes. Cette approche s’adapte bien à différents fonctions objectif et permet d'intégrer relativement facilement plusieurs contraintes. De nombreuses expérimentations ont été menées pour valider les méthodes proposées. / In most of the machine scheduling literature, resources are assumed to be continuously available which is not always true. We deal with the context of unavailability known a priori; we are particularly interested in job-shop scheduling problems with flexible unavailability periods and tasks that can eventually be interrupted by unavailability periods. Integrating these constraints increase the complexity of the scheduling problems. We deal with flexibility that is related to at least one of the following points: moving the unavailability period in a time window, modification of the duration of the unavailability period, interruption of a task by an unavailability period, then resumed with a possible penalty.In this thesis, we propose mathematical models for the problem. In addition to the resolution of the considered problems, the aim of this modeling is to allow for the analysis of the impact of different constraints and evaluation of the quality of the approximate methods and the column generation approach we develop. The approximate methods construct in very short time a schedule based on priority rules. The solutions are also used in our column generation approach. This approach adapts well to various objective functions and allows relatively easily for the integration of several constraints. Many experiments have been performed to validate the designed methods.

- Réalisé au centre de recherche de l'hospital du Sacré-Coeur de Montréal. - Programme conjoint entre Université de Montréal et École Polytechnique de Montréal. / La fibrillation auriculaire (FA) est la forme d’arythmie la plus fréquente chez les êtres humains. Les mécanismes qui gouvernent l’initiation et les manifestations de cette maladie sont complexes, de nature dynamique, incluant des interactions à travers multiples échelles temporelles et spatiales dans les oreillettes. Ceci conduit très souvent à des manifestations imprévisibles et à des phénomènes qui émergent à l’échelle de l’organe, et qui se reflètent à l’échelle de tout le torse. Pour remédier à ce problème, on peut effectuer une ablation par cathéter, qui consiste à créer sur le tissu auriculaire des lésions linéaires qui bloquent et contraignent la propagation électrique. Parfois, ces lignes se reconnectent quelque temps après l’intervention, ce qui mène à des récidives, nécessitant ainsi une nouvelle intervention. Le but de ce projet est de modéliser un suivi de l’onde P post-opératoire pour détecter de manière non-invasive la reconnexion des lignes d’ablation et ainsi prédire les récidives de fibrillation auriculaire. À l’aide d’un modèle mathématique des oreillettes et du thorax, les ondes P sont simulées avant et après ablation, ainsi qu’après reconnexion de certaines lignes d’ablation. Les résultats montrent que la morphologie et les caractéristiques de l’onde P, ainsi que la carte d’activation sont affectées significativement par l’ablation et les reconnexions subséquentes. Ces différences sont plus facilement détectables lorsque les reconnexions naissent sur la veine pulmonaire inférieure gauche. Les changements sont plus importants pour les électrodes placées sur certaines zones du torse, notamment dans le dos. Ces nouvelles données aident actuellement à la conception d’une étude clinique pour valider l’approche. / Atrial fibrillation (AF) is the most common form of arrhythmia in humans. The mechanisms governing the initiation and manifestations of that disease are complex, dynamic in nature, including interactions across multiple spatial and temporal scales in the atria. This often leads to unpredictable manifestations and phenomena that arise at the level of the organ, and are reflected across the entire torso. To remedy that problem, catheter ablation can be carried out, which consists in creating linear lesions which block and force the electrical propagation in the atrial tissue. Sometimes these lines reconnect after the procedure, which leads to atrial fibrillation recurrence, thus requiring a new intervention. The purpose of this work is to model the monitoring of the postoperative P wave to detect non-invasively the reconnection of ablation lines and to predict atrial fibrillation recurrences. Using a mathematical model of the atria and thorax, the P waves are simulated before and after ablation, as well as after reconnection of some ablation lines. The results show that the morphology and the characteristics of the P wave as well as the activation map are significantly affected by the ablation lines and the subsequent reconnections. These differences are more easily detected when reconnections arise on the left inferior pulmonary vein. The changes are most important in electrodes placed in certain areas of the torso, notably in the back. These new data are helping to plan a clinical study to validate the approach.

La drépanocytose, ou anémie falciforme, présente une variabilité interindividuelle considérable, conditionnée par de multiples facteurs, dynamiques et interactifs, depuis le niveau moléculaire jusqu’au niveau du patient. L’hémoglobine drépanocytaire, ou hémoglobine S (HbS, tétramère a2bS 2), est un mutant de l’hémoglobine A (a2b2) : elle possède à sa surface une valine (hydrophobe) substituant un acide glutamique natif (négativement chargé). Cette mutation entraîne l’agrégation de l’HbS désoxygénée en polymères, ainsi que l’altération des propriétés de l’érythrocyte -dont sa rhéologie et ses interactions avec les différentes cellules vasculaires. C’est pourquoi la polymérisation de l’HbS constitue un facteur étiologique clef, sinon le primum movens, de la drépanocytose, et une hypothèse thérapeutique (étayée par l’observation) postule que la réduction des fibres intra-érythrocytaires de HbS pourrait améliorer le statut clinique des patients en abaissant la fréquence et la sévérité des crises vasoocclusives. Dans l’optique de mieux comprendre et de mieux gérer la variabilité individuelle drépanocytaire, il apparaît donc indispensable de disposer, en premier lieu, d’une description réaliste de la polymérisation de l’HbS. L’objectif de ce travail de thèse est la mise en place et la validation d’un modèle mathématique de la polymérisation de l’HbS désoxygénée, en tant que processus cinétiquethermodynamique, sous l’influence de la concentration et de la température –les deux facteurs modulateurs les plus importants. A partir d’un modèle existant, mais linéaire et incomplet (Ferrone et al., 1985), nous avons procédé à son implémentation, à sa correction et à sa mise à jour, ainsi qu’à l’évaluation quantitative de ses performances dynamiques, par intégration complète et simulation numérique (Simulink©). Ceci nous a permis de réaliser un diagnostic et d’effectuer un certain nombre de raffinements, concernant en particulier (i) la voie de nucléation hétérogène (formation de néo-fibres sur les fibres préexistantes), (ii) la non-idéalité de la solution protéique de HbS, induite par le volume exclus des fibres polymères (coefficients d’activité calculé à partir de la « théorie des particules convexes »), ainsi que (iii) la structuration spatiale des polymères en domaines. Le modèle développé dans ce travail servira de base pour une description (i) de l’influence dynamique de l’oxygénation et des hémoglobines non-polymérisantes sur la polymérisation de HbS, puis (ii) des polymères de HbS sur les propriétés membranaires et rhéologiques de l’érythrocyte drépanocytaire. / Sickle cell disease pathology exhibits a strong interindividual variability, which depends upon multiple, dynamic and interacting factors, from the molecular to the patient level. Sickle hemoglobin, hemoglobin S (HbS, a2bS 2 tetramer), is a mutant of HbA (a2b2), with a surface valine (hydrophobic) substituting a native glutamic acid (negatively charged). Such a mutation endows deoxygenated HbS with the propensity to agregate into polymers, altering erythrocyte properties –including its rheology and its interactions with vascular and circulatory cells. Thus HbS polymerization is a key etiological factor of sickle cell disease, if not the primum movens. Indeed, one therapeutical hypothesis (supported by observation) postulates that the reduction of intra-erythrocytic HbS fibers could improve patients clinical status by lowering the frequency and the severity of vasooclusive crisis. In order to better understand and manage sickle cell disease variability, it is essential to have a realistic description of HbS polymerization. This work aims at developing and validating a mathematical model of deoxygenated HbS polymerization, as a kinetic and thermodynamic process under the influence of concentration and temperature –the two most important modulators. Building on an existing, but linearized and uncomplete (Ferrone et al., 1985) model, we have implemented, corrected and updated, and quantitatively evaluated its dynamical performances: this was done by full numerical integration using Simulink©. This allowed us to make several improvements, related in particular to : (i) the heterogeneous nucleation pathway (seeding and formation of new fibers from pre-existing ones), (ii) the non-ideality of the HbS protein solution, caused by polymer fibers excluded volume (activity coefficients were calculated with the CPT, Convex Particle Theory), and (iii) the spatial organization of polymers into domains. The model developped in this work will ground the description of the dynamic influence (i) oxygenation and non-polymerizing hemoglobins, (ii) HbS polymers interactions with membrane and consequences upon rheological properties of sickle cell erythrocyte.

Mitochondria are cell organelles which play an essential role in the cell energy supply providing the universal high energetic molecule ATP which is used in numerous energy consuming processes. The core of the ATP production, oxidative phosphorylation (OXPHOS) consists of four enzyme complexes (respiratory chain) which establish, driven by redox reactions, a proton gradient over the inner mitochondrial membrane. The ATP-synthase uses this electrochemical gradient to phosphorylate ADP to ATP. Dysfunctioning of an OXPHOS complex can have severe consequences for the energy metabolism and cause rare but incurable dysfunctions in particular tissues with a high energy demand such as brain, heart, kidney and skeleton muscle. Moreover mitochondria are linked to widespread diseases like diabetes, cancer, Alzheimer and Parkinson. Further, reactive oxygen species which are a by-product of the respiratory chain, are supposed to play a crucial role in aging. The aim of this work is to provide a realistic model of OXPHOS which shall help understanding and predicting the interactions within the OXPHOS and how a local defect (enzyme deficiency or modification) is expressed globally in mitochondrial oxygen consumption and ATP synthesis. Therefore we chose a bottom-up approach. In a first step different types of rate equations were analyzed regarding their ability to describe the steady state kinetics of the isolated respiratory chain complexes in the absence of the proton gradient. Here Michaelis-Menten like rate equations were revealed to be appropriate for describing their behavior over a wide range of substrate and product concentrations. For the validation of the equations and the parameter estimation we have performed kinetic measurements on bovine heart submitochondrial particles. The next step consisted in the incorporation of the proton gradient into the rate equations, distributing its influence among the kinetic parameters such that reasonable rates were obtained in the range of physiological electrochemical potential differences. In the third step, these new individual kinetic rate expressions for the OXPHOS complexes were integrated in a global model of oxidative phosphorylation. The new model could fit interrelated data of oxygen consumption, the transmembrane potential and the redox state of electron carriers. Furthermore, flux inhibitor titration curves can be well reproduced, which validates its global responses to local effects. This model may be of great help to understand the increasingly recognized role of mitochondria in many cell processes and diseases as illustrated by some simulations proposed in this work. / Les mitochondries sont l’usine à énergie de la cellule. Elles synthétisent l’ATP à partir d’une succession de réactions d’oxydo-réduction catalysées par quatre complexes respiratoires qui forment la chaîne respiratoire. Avec la machinerie de synthèse d’ATP l’ensemble constitue les oxydations phosphorylantes (OXPHOS). Le but de ce travail est de bâtir un modèle des OXPHOS basé sur des équations de vitesse simples mais thermodynamiquement correctes, représentant l’activité des complexes de la chaîne respiratoire (équations de type Michaelis- Menten). Les paramètres cinétiques de ces équations sont identifiés en utilisant les cinétiques expérimentales de ces complexes respiratoires réalisées en absence de gradient de proton. La phase la plus délicate de ce travail a résidé dans l’introduction du gradient de protons dans ces équations. Nous avons trouvé que la meilleure manière était de distribuer l’effet du gradient de proton sous forme d’une loi exponentielle sur l’ensemble des paramètres, Vmax et Km pour les substrats et les produits. De cette manière, j’ai montré qu’il était possible de représenter les variations d’oxygène, de ΔΨ et de ΔpH trouvés dans la littérature. De plus, contrairement aux autres modèles, il fut possible de simuler les courbes de seuil observées expérimentalement lors de la titration du flux de respiration par l’inhibiteur d’un complexe respiratoire donné.Ce modèle pourra présenter un très grand intérêt pour comprendre le rôle de mieux en mieux reconnu des mitochondries dans de nombreux processus cellulaires, tels que la production d’espèces réactives de l’oxygène, le vieillissement, le diabète, le cancer, les pathologies mitochondriales etc. comme l’illustrent un certain nombre de prédictions présentées dans ce travail.

Les traitements existants des métastases hépatiques du cancer colorectal montrent une efficacité insuffisante. Le projet GR5FU visait à améliorer cette efficacité et consistait à délivrer le 5-fluorouracile (5FU) dans le foie via son encapsulation dans des globules rouges (GR). Dans ce contexte, la modélisation visait à prédire la quantité de 5FU à encapsuler dans les GR pour atteindre une efficacité équivalente à celle du 5FU standard. Dans cette thèse, nous avons construit et implémenté un modèle mathématique multi-échelle qui relie l’injection du 5FU à son efficacité sur la croissance tumorale en intégrant sa pharmacocinétique et son mécanisme d’action intracellulaire. Des simulations de population de ce modèle, s’appuyant sur des paramètres de la littérature, nous ont permis (i) de reproduire des résultats cliniques montrant le pouvoir prédictif de l’enzyme Thymidylate Synthase (TS) et (ii) d’identifier deux prédicteurs potentiels de la réponse au 5FU à l’échelle d’une population virtuelle, en complément du niveau de TS : la vitesse de croissance tumorale et le métabolisme intracellulaire des pyrimidines. Nous avons également analysé, à l’aide de modèles à effets mixtes, (i) la croissance in vivo de la tumeur intra-hépatique VX2 sans traitement, tenant lieu de modèle animal de métastase hépatique, et (ii) la distribution plasmatique et hépatique du 5FU chez l’animal. Cette modélisation statistique nous a permis d’identifier les modèles décrivant des données expérimentales, d’estimer les paramètres de ces modèles et leur variabilité, et de générer une meilleure connaissance de la croissance de la tumeur VX2 et de la pharmacocinétique animale du 5FU, en particulier hépatique. Dans cette thèse, nous avons illustré comment l’intégration du métabolisme d’un médicament et de son mécanisme d’action dans un modèle global et la simulation de ce modèle à l’échelle d’une population virtuelle, constituent une approche prometteuse pour optimiser le développement d’hypothèses thérapeutiques innovantes en collaboration avec des expérimentateurs. / Existing treatments for liver metastases of colorectal cancer show a lack of efficacy. In order to improve the prognosis of patients, the GR5FU project has been implemented. It consisted in delivering the drug 5-fluorouracil (5FU) in the liver via its encapsulation in red blood cells (RBC) to increase its efficacy / toxicity ratio. In this context, the modeling aimed at predicting the amount of 5FU to encapsulate in RBC to achieve an efficacy equivalent to standard 5FU. In this thesis, we have created and implemented a multiscale mathematical model that links the injection of 5FU to its efficacy on tumor growth by integrating its pharmacokinetics and mechanism of intracellular action. Population simulations of this model, using parameters from the literature, allowed us (i) to reproduce clinical results showing the predictive power of TS enzyme level and (ii) to identify two potential predictors of response to 5FU at the level of a population of virtual patients, in addition to TS level. We also analyzed, using mixed effects models, (i) the in vivo growth of intrahepatic VX2 tumor without treatment, serving as an animal model of liver metastasis, and (ii) the distribution of 5FU in the animal’s organism. This statistical modelization enabled us to identify the models describing experimental data, to estimate the parameters of these models and their variability, and generate a better knowledge of VX2 tumor growth and animal 5FU pharmacokinetics. In this thesis, we illustrated how the integration of drug metabolism and its mechanism of action in a global model and the simulation of this model at the scale of a virtual population, form a promising approach to optimize the development of innovative therapeutic hypotheses in collaboration with experimentalists.

Résumé indisponible. / In this thesis we have concerned ourselves with university timetabling problems both course timetabling and examination timetabling problems. Most of the timetabling problems are computationally NP-complete problems, which means that the amount of computation required to find solutions increases exponentially with problem size. These are idiosyncratic nature problems, for example different universities have their own set of constraints, their own definition of good timetable, feasible timetable and their own choice about the use of constraint type (as a soft or hard constraint). Unfortunately, it is often the case that a problem solving approach which is successfully applied for one specific problem may not become suitable for others. This is a motivation, we propose a generalized problem which covers many constraints used in different universities or never used in literature. Many university timetabling problems are sub problems of this generalized problem. Our proposed algorithms can solve these sub problems easily, moreover constraints can be used according to the desire of user easily because these constraints can be used as reference to penalty attached with them as well. It means that give more penalty value to hard constraints than soft constraint. Thus more penalty value constraints are dealt as a hard constraint by algorithm. Our algorithms can also solve a problem in two phases with little modification, where in first phase hard constraints are solved. In this work we have preferred and used two phase technique to solve timetabling problems because by using this approach algorithms have broader search space in first phase to satisfy hard constraints while not considering soft constraints at all. Two types of algorithms are used in literature to solve university timetabling problem, exact algorithms and approximation algorithms. Exact algorithms are able to find optimal solution, however in university timetabling problems exact algorithms constitute brute-force style procedures. And because these problems have the exponential growth rates of the search spaces, thus these kinds of algorithms can be applied for small size problems. On the other side, approximation algorithms may construct optimal solution or not but they can produce good practically useable solutions. Thus due to these factors we have proposed approximation algorithms to solve university timetabling problem. We have proposed metaheuristic based techniques to solve timetabling problem, thus we have mostly discussed metaheuristic based algorithms such as evolutionary algorithms, simulated annealing, tabu search, ant colony optimization and honey bee algorithms. These algorithms have been used to solve many other combinatorial optimization problems other than timetabling problem by modifying a general purpose algorithmic framework. We also have presented a bibliography of linear integer programming techniques used to solve timetabling problem because we have formulated linear integer programming formulations for our course and examination timetabling problems. We have proposed two stage algorithms where hard constraints are satisfied in first phase and soft constraints in second phase. The main purpose to use this two stage technique is that in first phase hard constraints satisfaction can use more relax search space because in first phase it does not consider soft constraints. In second phase it tries to satisfy soft constraints when maintaining hard constraints satisfaction which are already done in first phase. (...)

Le procédé de refusion sous laitier électroconducteur (Electro Slag Remelting ou ESR) est aujourd’hui largement utilisé pour la production d’alliages métalliques à haute valeur ajoutée, comme les aciers spéciaux ou les superalliages base nickel. La modélisation mathématique et la simulation numérique du procédé ESR présentent un grand intérêt puisque les études expérimentales sur installations industrielles sont coûteuses et souvent difficiles à mettre en oeuvre. Ainsi, afin d’améliorer la compréhension et la maîtrise de la conduite d’une refusion, un modèle prédictif a été développé dans le cadre de cette étude. Il décrit les transferts couplés de chaleur et de quantité de mouvement lors de la croissance et de la solidification d’un lingot, en géométrie axisymétrique. La résolution des équations est basée sur une approche de type volumes finis. Le modèle tient compte de l’effet Joule dans le laitier résistif, des forces électromagnétiques et de la turbulence éventuelle de l’écoulement des phases liquides. La zone pâteuse est traitée comme un milieu poreux. Le modèle permet notamment de prédire la formation de la peau de laitier solide qui entoure le laitier et le lingot. Par ailleurs, il offre l’avantage de simuler le comportement du lingot et du laitier après la coupure finale du courant.Le développement s’est accompagné d’une importante étape de validation. Quatre refusions à l’échelle industrielle ont ainsi été réalisées à l’aciérie des Ancizes (Aubert&Duval). Les observations expérimentales ont ensuite été confrontées aux résultats du calcul. La comparaison a montré que le modèle peut être utilisé afin de prédire le comportement du procédé, à condition d’accorder une attention particulière à l’estimation des propriétés thermophysiques du métal, et surtout du laitier. Enfin, afin d’illustrer l’utilisation du modèle comme support à la compréhension du procédé, nous avons étudié l’influence de la variation de paramètres opératoires tels que la profondeur d’immersion de l’électrode, le taux de remplissage ou la pression de l’eau de refroidissement. / Electro Slag Remelting (ESR) is widely used for the production of high-value-added alloys such as special steels or nickel-based superalloys. Because of high trial costs and complexity of the process, trial-and-error based approaches are not well suitable for fundamental studies and optimization of the process.Consequently, a transient-state numerical model which accounts for electromagnetic phenomena and coupled heat and momentum transfers in an axisymmetrical geometry has been developed. The model simulates the continuous growth of the electroslag remelted ingot through a mesh-splitting method. In addition, solidification of the metal and slag is modelled by an enthalpy-based technique. A turbulence model is implemented to compute the motion of liquid phases (slag and metal), while the mushy zone is described as a porous medium whose permeability varies with the liquid fraction, thus enabling an accurate calculation of solid/liquid interaction. The coupled partial differential equations are solved using a finite-volume technique.Computed results are compared to experimental observation of 4 industrial remelted ingots fully dedicated to the model validation step. Pool depth and shape are particularly investigated in order to validate the model. Comparison shows that the model can be used as a predictive tool to analyse the process behavior. Nevertheless, it is necessary to pay a particular attention to the estimation of the thermophysical properties of metal and especially slag.These results provide valuable information about the process performance and influence of operating parameters. In this way, we present some examples of model use as a support to analyse the influence of operating parameters. We have studied the variation of electrode immersion depth, fill ratio and water pressure in the cooling circuit.

Dans une optique de maîtrise du procédé d’incinération des ordures ménagères et de ses possibles émissions polluantes, nous avons développé un modèle mathématique qui simule un lit d’ordures ménagères en combustion sur une grille mobile. Ce modèle décrit la plupart des phénomènes physicochimiques et thermiques intervenant lors de l’incinération : séchage et pyrolyse de la charge, combustion et gazéification du carbone résiduel, transferts thermiques, effondrement du lit, brassage… Il intègre également une description des mécanismes de volatilisation des métaux lourds et de formation des NOx. La cinétique de départ des métaux lourds est modélisée en tenant compte des différentes étapes de transport (transfert externe, diffusion intraparticulaire, volatilisation) au moyen de l’approche des temps caractéristiques additifs. Dans le cas simulé du cadmium, la prédiction d’une volatilisation quasi-complète est conforme aux résultats de la littérature. Le sous-modèle NOx prend en compte les mécanismes de formation thermique, prompt, combustible, par l’intermédiaire de N2O, ainsi que les mécanismes de réduction homogène par recombustion et hétérogène par le carbone résiduel. Les calculs révèlent que prédominent la formation par le mécanisme combustible et la destruction par la réduction hétérogène. Enfin, le modèle de lit a été utilisé pour étudier l’influence des divers paramètres opératoires : température, débit et distribution d’air primaire, taille des particules de déchets, brassage et schéma de brassage. Les résultats sont présentés et discutés en détail. L’influence des conditions opératoires sur l’efficacité du procédé et sur les émissions de Cd et NOx est analysée / As a tool for controlling the Municipal Solid Waste (MSW) incineration process and its possible pollutant emissions, a mathematical model of the MSW bed burning on travelling grate of an incinerator was developed. The model describes most of the physico-chemical and thermal phenomena taking place in incineration like the drying and pyrolysis of the feed, combustion and gasification of char, oxidation of pyrolysis gases, heat transfer, bed shrinking, feed stirring, etc. Also described in the model are the mechanisms of Heavy Metals (HM) volatilization and NOx formation. Kinetics of HM release was modelled using the approach of additive reaction times accounting for the various transport mechanisms involved: external transfer, intra-particle diffusion and actual volatilization. In the case simulated, i.e. of Cd, almost total volatilization is predicted, which is confirmed by literature findings. The NOx sub-model takes into account most of the common mechanisms of formation like thermal, prompt, fuel, N2O intermediate and also NOx reduction by homogeneous reburning and heterogeneous reduction by char. Calculations show that NOx formation is predominated by the fuel mechanism and destruction by the heterogeneous reduction. Finally, the bed model was applied to study the influence of various operating parameters like flow rate, temperature and distribution of air under grates, waste particle size, feed stirring and the stirring scheme. The results are presented and discussed in detail and the influence of operating conditions on process efficiency and on emissions of Cd and NOx is analyzed