Le démélange spectral est un des sujets majeurs de l’analyse d’images hyperspectrales. Ce problème consiste à identifier les composants macroscopiques présents dans une image hyperspectrale et à quantifier les proportions (ou abondances) de ces matériaux dans tous les pixels de l’image. La plupart des algorithmes de démélange suppose un modèle de mélange linéaire qui est souvent considéré comme une approximation au premier ordre du mélange réel. Cependant, le modèle linéaire peut ne pas être adapté pour certaines images associées par exemple à des scènes engendrant des trajets multiples (forêts, zones urbaines) et des modèles non-linéaires plus complexes doivent alors être utilisés pour analyser de telles images. Le but de cette thèse est d’étudier de nouveaux modèles de mélange non-linéaires et de proposer des algorithmes associés pour l’analyse d’images hyperspectrales. Dans un premier temps, un modèle paramétrique post-non-linéaire est étudié et des algorithmes d’estimation basés sur ce modèle sont proposés. Les connaissances a priori disponibles sur les signatures spectrales des composants purs, sur les abondances et les paramètres de la non-linéarité sont exploitées à l’aide d’une approche bayesienne. Le second modèle étudié dans cette thèse est basé sur l’approximation de la variété non-linéaire contenant les données observées à l’aide de processus gaussiens. L’algorithme de démélange associé permet d’estimer la relation non-linéaire entre les abondances des matériaux et les pixels observés sans introduire explicitement les signatures spectrales des composants dans le modèle de mélange. Ces signatures spectrales sont estimées dans un second temps par prédiction à base de processus gaussiens. La prise en compte d’effets non-linéaires dans les images hyperspectrales nécessite souvent des stratégies de démélange plus complexes que celles basées sur un modèle linéaire. Comme le modèle linéaire est souvent suffisant pour approcher la plupart des mélanges réels, il est intéressant de pouvoir détecter les pixels ou les régions de l’image où ce modèle linéaire est approprié. On pourra alors, après cette détection, appliquer les algorithmes de démélange non-linéaires aux pixels nécessitant réellement l’utilisation de modèles de mélange non-linéaires. La dernière partie de ce manuscrit se concentre sur l’étude de détecteurs de non-linéarités basés sur des modèles linéaires et non-linéaires pour l’analyse d’images hyperspectrales. Les méthodes de démélange non-linéaires proposées permettent d’améliorer la caractérisation des images hyperspectrales par rapport au méthodes basées sur un modèle linéaire. Cette amélioration se traduit en particulier par une meilleure erreur de reconstruction des données. De plus, ces méthodes permettent de meilleures estimations des signatures spectrales et des abondances quand les pixels résultent de mélanges non-linéaires. Les résultats de simulations effectuées sur des données synthétiques et réelles montrent l’intérêt d’utiliser des méthodes de détection de non-linéarités pour l’analyse d’images hyperspectrales. En particulier, ces détecteurs peuvent permettre d’identifier des composants très peu représentés et de localiser des régions où les effets non-linéaires sont non-négligeables (ombres, reliefs,...). Enfin, la considération de corrélations spatiales dans les images hyperspectrales peut améliorer les performances des algorithmes de démélange non-linéaires et des détecteurs de non-linéarités. / Spectral unmixing is one the major issues arising when analyzing hyperspectral images. It consists of identifying the macroscopic materials present in a hyperspectral image and quantifying the proportions of these materials in the image pixels. Most unmixing techniques rely on a linear mixing model which is often considered as a first approximation of the actual mixtures. However, the linear model can be inaccurate for some specific images (for instance images of scenes involving multiple reflections) and more complex nonlinear models must then be considered to analyze such images. The aim of this thesis is to study new nonlinear mixing models and to propose associated algorithms to analyze hyperspectral images. First, a ost-nonlinear model is investigated and efficient unmixing algorithms based on this model are proposed. The prior knowledge about the components present in the observed image, their proportions and the nonlinearity parameters is considered using Bayesian inference. The second model considered in this work is based on the approximation of the nonlinear manifold which contains the observed pixels using Gaussian processes. The proposed algorithm estimates the relation between the observations and the unknown material proportions without explicit dependency on the material spectral signatures, which are estimated subsequentially. Considering nonlinear effects in hyperspectral images usually requires more complex unmixing strategies than those assuming linear mixtures. Since the linear mixing model is often sufficient to approximate accurately most actual mixtures, it is interesting to detect pixels or regions where the linear model is accurate. This nonlinearity detection can be applied as a pre-processing step and nonlinear unmixing strategies can then be applied only to pixels requiring the use of nonlinear models. The last part of this thesis focuses on new nonlinearity detectors based on linear and nonlinear models to identify pixels or regions where nonlinear effects occur in hyperspectral images. The proposed nonlinear unmixing algorithms improve the characterization of hyperspectral images compared to methods based on a linear model. These methods allow the reconstruction errors to be reduced. Moreover, these methods provide better spectral signature and abundance estimates when the observed pixels result from nonlinear mixtures. The simulation results conducted on synthetic and real images illustrate the advantage of using nonlinearity detectors for hyperspectral image analysis. In particular, the proposed detectors can identify components which are present in few pixels (and hardly distinguishable) and locate areas where significant nonlinear effects occur (shadow, relief, ...). Moreover, it is shown that considering spatial correlation in hyperspectral images can improve the performance of nonlinear unmixing and nonlinearity detection algorithms.

La simulation numérique est devenue un outil central dans la modélisation de nombreux systèmes physiques tels que la dynamique des fluides, les plasmas, l'électromagnétisme, etc. L'existence de phénomènes multi-échelles rend l'intégration numérique de ces modèles très difficile du point de vue de la précision et du temps de calcul. En effet, dans les méthodes classiques d'intégration temporelle, le pas de temps est limité par la taille des plus petites mailles au travers d'une relation de type CFL. De plus, la forte disparité entre le pas de temps effectif et la condition CFL favorise les phénomènes de diffusion numérique. Dans la littérature, des nombreux algorithmes à pas de temps locaux (LTS) ont été développés. Pour la plupart des algorithmes LTS, les pas de temps locaux doivent être choisis parmi les fractions du pas de temps global. Nous présentons dans cette thèse une méthode asynchrone pour l'intégration explicite des équations différentielles multi-échelles. Cette méthode repose sur l'utilisation de critères de stabilité locaux, critères déterminés non pas globalement mais à partir de conditions CFL locales. De plus, contrairement aux schémas LTS, l'algorithme asynchrone permet la sélection de pas de temps indépendants pour chaque cellule de maillage. Cette thèse comporte plusieurs volets. Le premier concerne l'étude mathématique des propriétés du schéma asynchrone. Le deuxième a pour objectif d'étudier la montée en ordre, à la fois temporelle et spatiale, des méthodes asynchrones. De nombreux développements dans le cadre des méthodes de haute précision en temps ou en espace, telles que les méthodes de type Galerkin Discontinu, peuvent offrir un cadre naturel pour l'amélioration de la précision des méthodes asynchrones. Toutefois, les estimations garantissant l'ordre de précision de ces méthodes peuvent ne pas être directement compatibles avec l'aspect asynchrone. L'objectif de cette thèse est donc de développer un schéma numérique asynchrone d'ordre élevé mais qui permet également de limiter la quantité de calculs à effectuer. Le troisième volet de cette thèse se focalise sur l'application numérique puisqu'il concerne la mise en oeuvre de la méthode asynchrone dans la simulation des cas-tests représentatifs de problèmes multi-échelles. / Numerical simulation has become a central tool for the modeling of many physical systems (Fluid dynamics, plasmas, electromagnetism, etc). Multi-scale phenomena make the integration of these physical systems difficult in terms of accuracy and computational time. Numerical time-stepping integration techniques used for modeling such problems generally fall into two categories : explicit and implicit schemes. In the explicit schemes, all unknown variables are computed at the current time level from quantities already available. The time step used is then limited by the most restrictive CFL condition over the whole computation domain. In the implicit method the time step is no longer limited by the CFL conditions. However the scheme is generally not suitable for strongly coupled problems. To solve such problems, a number of local time-stepping (LTS) approaches have been developed. These methods are restricted by a local CFL condition rather than the traditional global CFL condition. For most of these LTS algorithms, local time steps are usually selected to be fractions of the global time step so that regular meeting points in time exist, and only little work is available on LTS methods with independent time steps. We present in this thesis an asynchronous method for the explicit integration of multi-scale partial differential equations. This method is restricted by a local CFL condition rather than the traditional global CFL condition. Moreover, contrary to other LTS methods, the asynchronous algorithm permits the selection of independent time steps in each mesh element. Our work consists of several components. The first one concerns the mathematical study of the properties of the asynchronous method. the objective of the second part is to study the improvement of the convergence rate for asynchronous methods. Many approaches in the context of high precision methods in time or in space, such as the Discontinuous Galerkin methods, may offer a natural setting to improve the precision of the asynchronous methods. However, the estimates ensuring the order of the accuracy of the method may not be directly compatible with the asynchronous aspect. Then, the objective is to develop a high order asynchronous numerical scheme which also preserves the computational time reduction. Finally, the third part is focused on the implementation of the asynchronous method and illustrate the advantages of the method on test-cases representative of multiscale problems.

Pour les accélérateurs de particules à hautes énergies, le système de verrouillage est une partie clé de la protection de la machine. Le verrouillage de la machine est l’inhibition du faisceau dès lors qu’un équipement critique tombe en panne et/ou qu’un faisceau est de faible qualité. Pour un système de verrouillage, sa sûreté de fonctionnement est la caractéristique la plus importante. Cette thèse présente le développement d’un système de verrouillage pour les collisionneurs linéaires avec une application au projet CLIC (Compact Linear Collider). Son élaboration s’appuie sur la norme d’ingénierie IEEE 1220 et se décline en quatre parties. Tout d’abord, les spécifications sont établies. Une attention particulière est portée sur la sûreté de fonctionnement, plus précisément, la fiabilité et la disponibilité du système. La deuxième étape est la proposition d’un design. Celui-ci est basé sur une analyse fonctionnelle, les interfaces du système et l’architecture du CLIC. Troisièmement, une étude de faisabilité est effectuée en appliquant les concepts dans un environnement opérationnel. Finalement, la dernière étape est la vérification matérielle. Le but est de prouver que le design proposé est capable de remplir le cahier des charges établi. / For high energy accelerators, the interlock system is a key part of the machine protection. The interlock principle is to inhibit the beam either on failure of critical equipment and/or on low beam quality evaluation. The dependability of such a system is the most critical parameter. This thesis presents the design of an dependable interlock system for linear collider with an application to the CLIC (Compact Linear Collider) project. This design process is based on the IEEE 1220 standard and is is divided in four steps. First,the specifications are established, with a focus on the dependability, more precisely the reliability and the availability of the system. The second step is the design proposal based on a functional analysis, the CLIC and interfaced systems architecture. Third, the feasibility study is performed, applying the concepts in an accelerator facility. Finally, the last step is the hardware verification. Its aim is to prove that the proposed design is able to reach the requirements.

Parmi les techniques usuelles de fouille de données, peu sont celles capables de tirer avantage de la complémentarité des dimensions pour des données sous forme de tableaux à plusieurs dimensions. A l'inverse les techniques de décomposition tensorielle recherchent spécifiquement les processus sous-jacents aux données, qui permettent d'expliquer les données dans toutes les dimensions. Les travaux rapportés dans ce manuscrit traitent de l'amélioration de l'interprétation des résultats de la décomposition tensorielle canonique polyadique par l'ajout de connaissances externes au modèle de décomposition, qui est par définition un modèle aveugle n'utilisant pas la connaissance du problème physique sous-jacent aux données. Les deux premiers chapitres de ce manuscrit présentent respectivement les aspects mathématiques et appliqués des méthodes de décomposition tensorielle. Dans le troisième chapitre, les multiples facettes des décompositions sous contraintes sont explorées à travers un formalisme unifié. Les thématiques abordées comprennent les algorithmes de décomposition, la compression de tenseurs et la décomposition tensorielle basée sur les dictionnaires. Le quatrième et dernier chapitre présente le problème de la modélisation d'une variabilité intra-sujet et inter-sujet au sein d'un modèle de décomposition contraint. L'état de l'art en la matière est tout d'abord présenté comme un cas particulier d'un modèle flexible de couplage de décomposition développé par la suite. Le chapitre se termine par une discussion sur la réduction de dimension et quelques problèmes ouverts dans le contexte de modélisation de variabilité sujet. / Among commonly used data mining techniques, few are those which are able to take advantage of the multiway structure of data in the form of a multiway array. In contrast, tensor decomposition techniques specifically look intricate processes underlying the data, where each of these processes can be used to describe all ways of the data array. The work reported in the following pages aims at incorporating various external knowledge into the tensor canonical polyadic decomposition, which is usually understood as a blind model. The first two chapters of this manuscript introduce tensor decomposition techniques making use respectively of a mathematical and application framework. In the third chapter, the many faces of constrained decompositions are explored, including a unifying framework for constrained decomposition, some decomposition algorithms, compression and dictionary-based tensor decomposition. The fourth chapter discusses the inclusion of subject variability modeling when multiple arrays of data are available stemming from one or multiple subjects sharing similarities. State of the art techniques are studied and expressed as particular cases of a more general flexible coupling model later introduced. The chapter ends on a discussion on dimensionality reduction when subject variability is involved, as well a some open problems.

On s'intéresse dans cette thèse à la couverture des produits dérivés dans des marchés incomplets. L'approche choisie peut se voir comme une extension des travaux de M. Schweizer sur la minimisation locale du risque quadratique. En effet, tout en restant dans le cadre de la modélisation des actifs par des semimartingales, notre méthode consiste à remplacer le critère de risque quadratique par un critère de risque plus général, sous la forme d'une fonctionnelle convexe du coût local. Nous obtenons d'abord des résultats d'existence, d'unicité et de caractérisation des stratégies optimales dans un marché sans friction, en temps discret et en temps continu. Puis nous explicitons ces stratégies dans le cadre de modèles de diffusion avec et sans sauts. Nous étendons également notre méthode au cas où la liquidité n'est plus infinie. Enfin nous montrons par le biais de simulations numériques les effets du choix de la fonctionnelle de risque sur la constitution du portefeuille optimal. / This thesis deals with the issue of hedging contingent claims in incomplete markets. The way we tackle this issue may be seen as an extension of M. Schweizer's work on quadratic local risk-minimization. Indeed, while still modelling assets as semimartingales, our method relies on the introduction of a convex function of the local costs to assess risk, thus relaxing the quadratic assumption. The results we obtain are existence and uniqueness results first and characterizations of optimal strategies in a frictionless market, both in discrete and continuous time settings. We then make those strategies explicit by using diffusion models with and without jumps. We further extend our approach in the case when liquidity is given through a stochastic supply curve. Finally we show the effect of the choice of different risk functions on the optimal portfolio by numerically solving the optimality equations.

Cette thèse traite de la réduction d’ordre de modèle du phénomène communément rencontré dans la modélisation de microsystèmes, à savoir, dans la littérature anglaise, le « squeeze-film damping ». Dans un premier chapitre sont présentées les différentes méthodes de réduction d’ordre de modèle. Dans le cas des systèmes linéaires, elles ont un cadre théorique bien établi. Ces méthodes peuvent être adaptées pour les systèmes non-linéaires. La validité des modèles réduits résultants sera alors réduite à un certain espace des phases, leur établissement faisant intervenir certaines trajectoires particulières servant d’apprentissage. On présente finalement la méthode des modes normaux non-linéaires dont les modèles résultants ne dépendent pas d’une trajectoire d’apprentissage. Au chapitre 2, on s’intéresse plus particulièrement au phénomène de « squeeze-film damping » régi par l’équation de Reynolds. Après avoir détaillé son établissement à partir de certaines hypothèses, on décrit les différentes méthodes de résolution de l’équation linéaire puis non-linéaire de la littérature. On compare ensuite les résultats d’un modèle de l’équation de Reynolds à des simulations éléments finis de l’équation de Navier-Stokes afin de valider les hypothèses faites pour la dérivation de l’équation de Reynolds. On propose ensuite une résolution originale par changement de variable. On étudie aussi plusieurs autres résolutions possibles ainsi que plusieurs bases de projection parmi celles décrites dans le premier chapitre. Le chapitre 3 est consacré à la modélisation du problème couplé que constitue le micro-interrupteur MEMS qui est un candidat au remplacement des interrupteurs à base de transistors dans les communications RF. Sa modélisation fait intervenir trois domaines, la mécanique, l’électrostatique, et la fluidique à travers l’équation de Reynolds. Après voir décrit les différents modèles de la littérature, on propose un modèle réduit couplé dont le modèle fluidique est basé sur le modèle établi au chapitre 2. Ce modèle est validé par rapport à des modèles différences finies et à des résultats expérimentaux de la littérature.Enfin le quatrième chapitre traite de la réduction du coût d’évaluation du modèle réduit couplé de micro-interrupteur du chapitre 3. La première méthode proposée consiste à trouver une fonction d’approximation de la projection de la force fluidique sur le premier mode mécanique, fonction des coordonnées modales mécaniques position et vitesse. Cette méthode ne se révèle valable que dans le cas incompressible. Dans le cas compressible, la résolution de l’équation de Reynolds restant obligatoire, on utilise la méthode de Rewienski et al. qui consiste à linéariser par morceaux les fonctions régissant la dynamique. Une autre méthode de linéarisation par morceaux, tirant parti d’une particularité du modèle du chapitre 2 et permettant de s’affranchir d’une trajectoire d’apprentissage, est également proposée. / This thesis deals with reduced-order modelling of squeeze-film damping, a fluidic phenomenon that is commonly encountered in MEMS. In the first chapter, reduced-order modelling methods are presented. For linear systems, well-established theories exist. They can be adapted to nonlinear systems. However, the resulting reduced-order models are valid in a certain region of the state-space only, depending on the training trajectory. The method of nonlinear normal modes, which does not depend on a training trajectory is also introduced. Chapter two is focused on the squeeze-film damping phenomenon governed by the Reynolds equation. We first establish the equation from appropriate hypotheses, and then present the different resolutions of its linear and nonlinear form found in literature. The results from a model based on the Reynolds equation are then compared to results from a finite element Navier-Stokes model, in order to validate the various hypotheses made. An original method of resolution based on a change of variable is then proposed. Several other method of resolution are studied as well as different projection bases amongst those presented in chapter one.Chapter three is dedicated to the modelling a micro-switch, a candidate to the replacement of switches based on transistors in RF communications. Its modelling implies the coupling of three domains: mechanics, electrostatics, and fluidics with Reynolds equation. Following a description of the models from literature, a coupled model is proposed, the fluidic model being the one presented in chapter two. This model is validated compared to finite difference models as well as experimental data from the literature.Finally, the fourth chapter aims at reducing the evaluation cost of the coupled micro-switch model established in chapter three. The first method consists in finding an approximation function of the projection of the fluidic force on the first linear mechanical mode as a function of the mechanical modal coordinates, position and speed. This method is applicable in the incompressible case only. In the compressible case, the Reynolds equation has to be solved. The method of Rewienski and al., which consists in piecewise-linearizing the functions governing the dynamics, is used. Another method based on a piecewise-linear approach, taking advantage of the particular structure of the model presented in chapter two, thus not depending on a training trajectory, is proposed.

Thèse numérisée par la Direction des bibliothèques de l'Université de Montréal.

Ce travail de thèse contribue à l'étude acoustique expérimentale des propriétés élastiques des milieux granulaires et de leurs évolutions lentes avec le temps faisant suite à des sollicitations mécaniques. Grâce à un montage expérimental et des protocoles expérimentaux précisément décrits, une étude des caractéristiques acoustiques linéaires et non linéaires d'une tranche de milieu granulaire est réalisée. Des zones fréquentielles caractéristiques de la fonction de transfert acoustique sont identifiées. Dans une zone basse fréquence, les propriétés élastiques du squelette solide (les billes et leurs contacts) régissent le comportement de la fonction de transfert acoustique. A plus hautes fréquences, seuls les paramètres acoustiques du fluide équivalent jouent un rôle. Il est montré que le couplage des modes de propagation au niveau de la détection du signal se manifeste à la transition entre ces zones caractéristiques. Les effets non linéaires d'auto-action et génération d'harmonique 2 dans cette tranche sont analysés pour différentes compacités. L'étude des effets de mémoire d'un échantillon granulaires est réalisée lors de la modification du protocole de compaction (l'amplitude de la sollicitation mécanique change de valeur). Cette modification entraîne des variations soudaines des paramètres élastiques et dissipatifs, linéaires et non linéaires, qui sont interprétées en termes de forces de contacts et distribution de forces de contacts. Pour la première fois, les effets de mémoire sont analysés via l'élasticité du milieu granulaire et non simplement via sa géométrie. Enfin, la relaxation lente des propriétés élastiques, linéaires et non linéaires, est observée suite à une unique sollicitation mécanique, un "tap". L'augmentation du paramètre élastique linéaire, la diminution de l'atténuation linéaire et les diminutions des paramètres non linéaires sont recueillies avec une résolution en temps de l'ordre de la seconde. Les rôles de la température et de l'hygrométrie de l'air ambiant dans les temps caractéristiques de relaxation qui sont de l'ordre de plusieurs minutes sont analysés. / This work is a contribution to the experimental acoustic study of the elastic properties of granular media and their slow evolutions during time after mechanical solicitations. With an experimental setup and experimental methods specially adapted, a study of the linear and nonlinear acoustic characteristics of a granular slab is done. Characteristic frequency regions of the acoustic transfer function are identified. In a low frequency band (<10 kHz) the elastic properties of the solid skeleton (beads and theirs contacts) govern the behavior of the acoustic transfer function. In a higher frequency range (15 - 30 kHz), only the acoustic parameters of the equivalent fluid play a role. It is shown that the coupling of propagation modes in the detected signal is visible at the transition between these characteristic frequency regions. The nonlinear effects of self-action and harmonic generation in this slab are analyzed for several compacities. Then the study of the memory effects in a granular sample is performed during the modification of the compaction process (the amplitude of the mechanical solicitation is changed). This modification makes sudden variations of the linear and nonlinear elastic and dissipative parameters. These modifications are interpreted in terms of contact forces and contact force distribution. For the first time, the memory effects are not simply analyzed via packing geometry effects but via elasticity of the granular packing. Finally, the slow relaxation of the linear and nonlinear elastic properties is observed after one single mechanical solicitation called "tap". The increase in the linear elastic parameter, the decrease in the linear attenuation and the decrease in the nonlinear parameters are observed with a less than one second time resolution. Roles of the temperature and the hygrometry of the ambient air in the characteristic relaxation times which are lasting several minutes are analyzed.

Cette thèse est consacrée à l'étude de la théorie des barrières pour les systèmes non linéaires sous contraintes d'entrées et d'état. La principale contribution concerne la généralisation au cas de contraintes mixtes, c'est-à-dire dépendant des entrées et de l'état de façon couplée. Ce type de contraintes apparaît souvent dans les applications et dans les systèmes différentiellement plats sous contraintes. On prouve un théorème du type principe du minimum qui permet de construire la barrière et l'ensemble admissible associé. De plus, dans le cas d'intersection de plusieurs trajectoires ainsi construites, on démontre que les points intersections transversaux sont des points d'arrêt de la barrière. Ces résultats sont utilisés pour calculer l'ensemble admissible d'un pendule inversé avec un câble non-rigide monté sur un chariot, la contrainte correspondant au fait que le câble reste tendu. Ce problème correspond en fait à la détermination de l'ensemble potentiellement sûr dans le cadre des systèmes hybrides. / This thesis deals with the theory of barriers in input and state constrained nonlinear systems. Our main contribution is a generalisation to the case where the constraints are mixed, that is they depend on both the input and the state in a coupled way. Constraints of this type often appear in applications, as well as in constrained flat systems. We prove a minimum-like principle that allows the construction of the barrier and the associated admissible set. Moreover, in case of intersection of some of the trajectories involved in this principle, we prove that such transversal intersection points are stopping points of the barrier.We demonstrate the utility of all the theoretical contributions by finding the admissible set for the inverted pendulum on a cart with a non-rigid cable, the constraint being that the cable remains taut. Note that this problem corresponds to the determination of potentially safe sets in hybrid systems.

This work deals with the different manifestations of cooperativity in pedopsychiatry in connection with severe disorders in adolescence such as suicide attempts and runaways, and developmental processes in early childhood. Our principal hypothesis is that cooperativity, nonlinearity and the resulting feedback loops underlie the appearance of suicidal peaks, the transitions between different mental states in a population of susceptible, suicidal and runaway individuals, or finally the selfish to egalitarian transition in young children. The principal methodology put forward is mathematical modeling, which allows one to identify the principal mechanisms present, to assess the role of the environmental constraints and to make predictions. This provides in turn insights on possible prevention, managemant and therapy strategies. Using this approach we have shown that suicide contagion can take an explosive form that is dramatically accelerated by the presence of virtual contacts as it happens, in particular, in communication via the Internet. Different scenarios have been explored, inspired by typical cases encountered in clinical practice when taking in charge adolescents in crisis. Regarding the population of susceptible, suicidal and runaways, a model has been developed whose analysis shows that depending on the parameter values and the size of the initial population the system can evolve toward states where the populations of the suicidal and of the runaways are negligible or, on the contrary, reach appreciable values. Finally, we have shown that in absence of interactions between individuals the selfish to egalitarian transition in children occurs gradually, whereas in the presence of cooperativity the transition is more abrupt. Our results suggest that when taking in charge adolescents in a state of risk contagion and transition between mental states must be recognized as additional iatrogenic risks, which can be minimized by privileging small non-mixed groups. / Doctorat en Sciences médicales (Médecine) / info:eu-repo/semantics/nonPublished