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

Ce travail de thèse porte sur le développement de la méthode des vitesses discrètes pour la résolution numérique de équations cinétiques modèles, BGK, S modèle et ES modèle, qui représentent les différentes approximations de l’équation de Boltzmann. / This work is devoted to the development of the numerical resolution of the kinetic model equations such as BGK, S-model, ES-model by the discrete velocity method. The different approximations of the Boltzmann equation are presented

La contribution du marketing management à la stratégie est une problématique qui a suscité un grand intérêt, aussi bien dans la communauté des chercheurs que celle des praticiens. L’orientation marché a été un sujet d’investigation fertile dans ce courant de recherche et des débats ont été lancés sur sa capacité à influencer le design stratégique de l’entreprise. Toutefois, peu d’études empiriques ont analysé la relation entre l’orientation stratégique et les comportements orientés marché. Le présent travail tente de contribuer à la théorie de l’interface en examinant cette relation et son impact sur la performance de l’entreprise. Il intègre par ailleurs l’analyse du poids de certains modérateurs qui influencent cette association. Cette recherche doctorale présente une étude empirique qui répond à ces questionnements. Nous avons menée une étude qui a testé notre modèle auprès de 254 senior managers et nous avons opté pour une modélisation en équations structurelles. Nos résultats débouchent sur de nombreuses implications managériales. / Recognized as a topic with important implications for executive practice, the contribution of marketing to business strategy management is a fertile topic of investigation. Within this research stream the role of market orientation has deserved special attention. Prior research raised a substantial debate regarding the market orientation capability to influence strategy design. However, few have investigated this interaction effect and little is understood of the business performance implications of market orientation‘s contribution to business strategy formation. This research addresses theseshortfalls in the literature by examining relationships between market orientation and strategic orientation, their interaction effect on performance, and some factors that moderate their effects on performance. In an attempt to address these key issues, this doctoral research presents an empirical investigation. We conducted a survey of 254 senior executives about their practices adopted in their respective firms. We next performed a structural equation modeling analysis. The results highlight many managerial implications for the practitioners.

Ce travail est consacré à l'étude du comportement en temps grand d'équations aux dérivées partielles de type parabolique. Plus particulièrement, on s'intéresse à des équations non linéaires de type diffusion, qui interviennent dans de nombreux modèles issus de la physique (par exemple l'équation des milieux poreux) ou de la biologie (par exemple le modèle de Patlak-Keller-Segel pour la chimiotaxie). Dans les chapitres I et II on s'intéresse à une amélioration de l'inégalité de Sobolev à travers son inégalité duale, l'inégalité de Hardy-Littlewood-Sobolev, dans le cadre du laplacien ordinaire et du laplacien fractionnaire, respectivement. Le chapitre III est un passage en revue de l'inégalité d'Onofri, qui joue le rôle de l'inégalité de Sobolev pour la dimension deux. De nouveaux résultats sont apportés, dont certains sont étendus aux variétés riemanniennes au chapitre IV. Enfin, le chapitre V traite des états stationnaires de deux modèles paraboliques, utilisés pour l'étude du déplacement de foules et la modélisation en biologie (chimiotaxie). / This work is dedicated to the study of the large time behaviour of some parabolic type partial differential equations. More specifically, we look into non linear diffusion equations that appear in a number of models arising in physics (e.g. the porous medium equation) or biology (e.g. the Patlak-Keller-Segel model for chemotaxis)Chapters I and II deal with an improved Sobolev inequality by means of its dual, the Hardy-Littlewood-Sobolev inequality, in the framework of the standard and fractional Laplacian, respectively. Chapter III is a review of the Onofri inequality,which acts as the Sobolev inequality for dimension two. New results are provided, and some of them are extended to Riemannian manifolds in Chapter IV. Finally, Chapter V deals with the stationary states of two parabolic models, used for thestudy of crowd motion and modeling in biologie (chemotaxis).

Cette thèse traite de questions mathématiques posées par des problèmes issus du micromagnétisme ; un thème central en est les champs de vecteur de rotationnel nul et de norme 1, qu'on voit naturellement apparaître comme configurations minimisant des énergies micromagnétiques.Le premier chapitre est motivé par la question suivante : peut-on, en dimension plus grande que deux, caractériser les champs de vecteur de rotationnel nul et de norme 1 par une formulation cinétique ?Une telle formulation a d'abord été introduite en dimension 2 dans l'article \cite{Jabin_Otto_Perthame_Line_energy_2002} de Jabin, Otto et Perthame où elle apparaît naturellement dans le cadre de la minimisation d'une énergie de type Ginzburg-Landau. Ignat et De Lellis ont ensuite montré dans \cite{DeLellis_Ignat_Regularizing_2014} qu'une telle formulation cinétique caractérise les champs de rotationnel nul et de norme 1 possédant une certaine régularité en dimension 2. Le premier chapitre de cette thèse est consacré à l'étude d'une formulation cinétique similaire en dimension quelconque ; le résultat principal en est qu'en dimension strictement plus grande que 2, cette fomulation cinétique ne caractérise non plus tous les champs de rotationnel nul et de norme 1, mais seulement les champs constants ou les vortex.La caractérsation cinétique des champs de vecteur de rotationnel nul et de norme 1 en dimension 2,prouvée par De Lellis et Ignat et que nous venons de mentionner reposait sur la notion d'entropie.Ayant obtenu une formulation cinétique en dimension quelconque, il était naturel de vouloir l'exploiter un tentant d'étendre également la notion d'entropie aux dimensions supérieures à 2. C'est ce à quoi est consacré le deuxième chapitre de cette thèse ; nous y définissons en particulier une notion d'entropie en dimension quelconque. Le point central en est la caractérisation de ces entropies par un système d'\équations aux dérivées partielles, et leur description complète en dimension 3, ainsi que la preuve pour ces entropies de propriétés tout à fait semblables à celles des entropies deux dimensionnelles.Le troisième chapitre de cette thèse, qui expose les résultats d'un travail en collaboration avec Antonin Monteil, s'intéresse à la minimisation d'\'energies de type Aviles-Giga de la forme $\mathcal_f(m)=\int_f(|m^+-m^-|)$ o\`u $m$ est un champ de rotationnel nul et de norme 1 et où $J(m)$ désigne les lignes de saut de $m$. Deux questions classiques se posent pour ce type d'énergie : la solution de viscosité de l'équation eikonale est-elle un minimiseur et l'énergie est-elle semi-continue inférieurement pour une certaine topologie. Le résutat principal de cette partie est un construction, qui nous permet en particulier de répondre par la négative à ces deux questions dans les cas où $f(t)= t^p$ avec $p \in ]0,1[$ en donnant une condition nécessaire sur $f$ pour que $\mathcal_f$ soit semi-continue inférieurement.Enfin, le dernier chapitre de cette thèse est consacré à l'étude d'une variante de l'énergie de Ginzburg-Landau introduite par Béthuel, Brezis et Helein où on a remplacé la condition de bord par une pénalisation dépendant d'un paramètre. Nous y décrivons le comportement asymptotique de l'énergie minimale qui, suivant la valeur de ce paramètre, soit se comporte comme l'énergie de Ginzburg-Landau classique en privilégiant une configuration vortex, soit privilégie au contraire une configuration singulière suivant une ligne. / This thesis is motivated by mathematical questions arising from micromagnetism. One would say that a central topic of this thesis is curl-free vector fields taking value into the sphere. Such fields naturally arise as minimizers of micromagnetic-type energies. The first part of this thesis is motivated by the following question : can we find a kinetic formulation caracterizing curl-free vector fields taking value into the sphere in dimension greater than 2 ? Such a formulation has been found in two dimension by Jabin, Otto and Perthame in \cite. De Lellis and Ignat used this formulation in \cite{DeLellis_Ignat_Regularizing_2014} to caracterize curl-free vector fields taking value into the sphere with a given regularity. The main result of this part is the generalization of their kinetic formulation in any dimension and the proof that if $d>2$, this formulation caracterizes only constant vector fields and vorteces, i. e. vector fields of the form $\pm \frac$. The second part of this thesis is devoted to a generalization of the notion of \textit, which plays a key role in the article of De Lellis and Ignat we talked about above. We give a definition of entropy in any dimension, and prove properties quite similar to those enjoyed by the classical two-dimensional entropy. The third part of this thesis, which is the result of a joint work with Antonin Monteil, is about the study of an Aviles-Giga type energy. The main point of this part is a necessary condition for such an energy to be lower semi continuous. We give in particular an example of energy of this type for which the viscosity solution of the eikonal equation is \textit a minimizer. The last part, finally is devoted to the study of a Ginzburg-Landau type energy where we replace the boundary condition of the classical Ginzburg-Landau energy introduced by Béthuel, Brezis and Helein by a penalization within the energy at the critical scaling depending on a parameter. The core result of this part is the description of the asymptotic of the minimal energy, which, depending on the parameter, favorizes vortices-like configuration like in the classical Ginzburg-Landau case, or configurations singular along a line.

Cette thèse est consacrée à l'étude de problèmes de stabilisation exponentielle par retour d'état ou "feedback" des équations de Navier-Stokes dans un domaine borné Ω ⊂ Rd, d = 2 ou 3. Le cas d'un contrôle localisé sur la frontière du domaine est considéré. Le contrôle s'exprime en fonction du champ de vitesse à l'aide d'une loi de feedback non-linéaire. Celle-ci est fournie grâce aux techniques d'estimation a priori via la procédure de Faedo-Galerkin laquelle consiste à construire une suite de solutions approchées en utilisant une base de Galerkin adéquate. Cette loi de feedback assure la décroissance exponentielle de l'énergie du problème discret correspondant et grâce au résultat de compacité, nous passons à la limite dans le système satisfait par les solutions approchées. Le chapitre 1 étudie le problème de stabilisation des équations de Navier- Stokes autour d'un état stationnaire donné, tandis que le chapitre 2 examine le problème de stabilisation autour d'un état non-stationnaire prescrit. Le chapitre 3 est consacré à l'étude de la stabilisation du problème de Navier-Stokes avec des conditions aux bords mixtes (Dirichlet- Neumann) autour d'un état d'équilibre donné. Enfin, nous présentons dans le chapitre 4, des résultats numériques dans le cas d'un écoulement autour d'un obstacle circulaire / In this thesis we study the exponential stabilization of the two and three-dimensional Navier- Stokes equations in a bounded domain Ω, by means of a boundary control. The Control is expressed in terms of the velocity field by using a non-linear feedback law. In order to determine a feedback law, we consider an extended system coupling the Navier-Stokes equations with an equation satisfied by the control on the domain boundary. While most traditional approaches apply a feedback controller via an algebraic Riccati equation, the Stokes-Oseen operator or extension operators, a Galerkin method is proposed instead in this study. The Galerkin method permits to construct a stabilizing boundary control and by using energy a priori estimation technics, the exponential decay is obtained. A compactness result then allows us to pass to the limit in the nonlinear system satisfied by the approximated solutions. Chapter 1 deals with the stabilization problem of the Navier-Stokes equations around a given steady state, while Chapter 2 examines the stabilization problem around a prescribed non-stationary state. Chapter 3 is devoted to the stabilization of the Navier-Stokes problem with mixed-boundary conditions (Dirichlet-Neumann), around to a given steady-state. Finally, we present in Chapter 4, numerical results in the case of a flow around a circular obstacle

Dans cette thèse nous étudions trois modèles décrivant la dynamique de l’écoulement d’un fluide à densité variable, dans des échelles spatio-temporelles grandes. Dans ce cadre, le mouvement relatif induit par des forces extérieures,comme la force de Coriolis ou la poussée hydrostatique, s’avère être beaucoup plus important que le mouvement intrinsèque du fluide induit par le transport des particules. Une tel déséquilibre contraint ainsi le mouvement, induisant des structures persistantes dans l’écoulement du fluide.D’un point de vue mathématique, l’une des difficultés consiste en l’étude des perturbations induites par les forces extérieures, qui se propagent à grande vitesse.Ce type d’analyse peut être effectué au moyen de plusieurs outils mathématiques ;on choisit ici d’employer des techniques caractéristiques de l’analyse de Fourier,comme l’analyse des propriétés dispersives des intégrales oscillantes.Tout au long de cette thèse, on se restreint à considérer des domaines spatiaux sans frontière : c’est le cas de l’espace entier, ou encore de l’espace périodique. Les modèles considérés sont donc les suivants: équations primitives dont les nombres de Froude et de Rossby sont comparables,et pour lesquelles la diffusion verticale est nulle, fluides stratifiés dans un régime à faible nombre de Froude, fluides faiblement compressibles et tournants dans un régime où les nombres de Mach et de Rossby sont comparables.On prouve que ces systèmes propagent globalement dans le temps des donnés peu régulières. Nous n’imposons jamais de condition de petitesse sur les données initiales. Toutefois, on prendra en compte certaines hypothèses spécifiques de régularité, lorsque des raisons techniques l’imposent. / In this thesis we discuss three models describing the dynamics of density-dependent fluids in long lifes pans and on a planetary scale. In such setting the relative displacement induced by various external physical forces, such as the Coriolis force and the stratification buoyancy, is far more relevant than the intrinsic motion generated by the collision of particles of the fluid itself. Such disproportion of balance limits hence the motion, inducing persistent structures in the velocity flow.On a mathematical level one of the main difficulties relies in giving a full description of the perturbations induced by the external forces, which propagate at high speed. This analysis can be performed by the aid of several tools, we chose here to adopt techniques characteristic of harmonic analysis, such as the analysis of the dispersive properties of highly oscillating integrals.All along the thesis we consider boundary-free, three-dimensional domains, and inspecific we study only the case in which the domain in either the whole space or the periodic space . The models we consider are the following ones : primitive equations with comparable Froude and Rossby number and zero vertical diffusivity, density-dependent stratified fluids in low Froude number regime, weakly compressible and fast rotating fluid in a regime in which Mach and Rossbynumber are comparable. We prove that these systems propagate globally-in-time data with low-regularity. Nosmallness assumption is ever made, specific constructive hypothesis are assumed on the initial data when required.

Le handicap de l’enfant peut avoir des conséquences sur des dimensions de la vie familiale : santé physique et psychologique, relations sociales et familiales, temps de prise en charge, emploi et poids économique. Mieux connaître les facteurs liés à cet impact permettrait d'en limiter les conséquences négatives sur la vie familiale.L’impact familial de la paralysie cérébrale (PC) d’adolescents européens a été étudié en 2009 et des facteurs personnels, sociaux et environnementaux de leurs familles également recueillis. Un questionnaire sur l’impact familial du handicap dans 5 des 7 dimensions identifiées dans la littérature a été utilisé. Etendu aux dimensions manquantes par l’ajout de 4 questions, l’outil amélioré a été validé dans le cadre de cette enquête européenne.Un modèle théorique fondé sur le modèle d’adaptation Double ABCX a été confronté aux données. Sont mis en relation des éléments de perturbations de l’équilibre familial, liées à l’enfant ou à l’environnement (A), des ressources disponibles pour faire face (B), la perception de conséquences négatives (C) et le degré d’adaptation familiale, mesurée par la détresse psychologique (X). La modélisation par équations structurelles a permis de valider un modèle estimant des coefficients de corrélation significatifs entre les variables selon une méthode itérative par moindres carrés pondérés (WLSMV). 57,4% de la variance de la détresse psychologique a été expliquée par des perturbations (A) : âge avancé de l’adolescent (bêta=0.19), absence d’emploi des parents (bêta=0.25), sévérité de la déficience motrice (bêta=0.23), sévérité des troubles du comportement de l’adolescent (bêta=0.45) ; par des ressources (B) : accès à des modes de garde extérieurs (bêta=-0.15), degré de dysfonctionnement familial (bêta=0.40), sentiments positifs du parent par rapport au handicap (bêta=-0.16) ; et enfin par l’impact négatif perçu (C) sur la santé des membres de la famille (bêta=0.45). Des éléments de (A) et de (B) influençaient l’impact négatif perçu dans d’autres dimensions de vie familiale. D’après ce modèle, l’adaptation de la famille face au handicap de l’adolescent semble se jouer plus autour de la qualité du fonctionnement familial que sur la mise en place de compensations physiques ou financières. Il est donc important de mettre en place des actions permettant aux familles d’améliorer les différentes dimensions de ce fonctionnement. L’utilisation d’outils de mesure validés peut aider à évaluer et sélectionner des interventions les plus efficaces pour gérer les troubles du comportement chez l’enfant en situation de handicap, augmenter l’autonomie de l’enfant, favoriser le bon fonctionnement familial, ou promouvoir les attitudes parentales positives face au handicap.L’absence d’emploi influençant le bien-être psychologique, cette situation a été étudiée pour les mères d’enfants en situation de handicap. A partir de données disponibles, nous avons constaté qu’en tenant compte de nombreuses caractéristiques de l’enfant et de la famille, la présence d’une déficience intellectuelle chez l’enfant et un bas niveau d’études de la mère ou le non-emploi du père influençaient significativement le risque de non-emploi de la mère.Par ailleurs, les mères d’enfants en situation de handicap sévère se sont avérées proportionnellement plus nombreuses que la population générale à être sans emploi (47,2% contre 31,4%). En revanche, les mères d’enfants avec PC interviewées lors d’une enquête n’étaient pas plus sans emploi que les mères de la population générale, remettant en question l’usage de la PC comme modèle du handicap, pour l’étude des conséquences familiales.L’étude de l’impact du handicap est encore trop rare en Europe et les effets de possibles interventions mal connues. Il est important de développer la recherche dans ce domaine pour améliorer le quotidien de ces familles et offrir à l’enfant en situation de handicap le meilleur environnement possible pour son développement. / The child's disability could impact family life in several dimensions: physical and psychological health, social and family relationships, time management, employment and economic weight. Knowledge about the factors related to this effect will limit the negative impact on family life.The family impact of the cerebral palsy (CP) in European adolescents was studied in 2009 and the personal, social and environmental factors of these families were also collected. A questionnaire measuring 5 dimensions of the 7 identified in literature was used. Extended to missing dimensions by adding four questions, this improved tool has been validated in the framework of this European survey.A theoretical model based on the Double ABCX model of family adaptation has been compared to the situation observed by these families, linking the stressors related to the child or the family environment (A), the resources available to cope with (B), the negative perceived impact in different dimensions of the family life (C) and the level of family adaptation, measured by the parental psychological distress (X). Structural equation modelling was used to validate a model estimating significant correlation coefficients between variables by an iterative method using weighted least squares (WLSMV). 57.4% of the variance in psychological distress was explained by stressors (A): age of the adolescent (beta=0.19), parental unemployment (beta=0.25), the severity of motor impairment (beta=0.23), the severity of behavioural disorders in adolescents (beta=0.45); by resources (B), the access to support in the childcare (beta=-0.15), the level of family dysfunction (beta=0.40), the parental positive attitude towards disability (beta=-0.16); and finally by the negative perceived impact (C) on the health of family members (beta=0.45). Some elements from (A) and (B) explained the perceived family impact in various dimensions of family life. Based on this model, in adolescence the family adaptation to disability seems to be more played out on quality of family functioning than on physical adaptation or financial compensation. Thus it is important to implement process to enable families to improve different aspects of this functioning. The use of validated measurement tools can help to assess and select the most effective interventions to manage behavioural problems in children with disabilities, to increase the autonomy of the child, to promote family functioning or parental positive attitudes towards disability.A lack of employment affecting the psychological well-being, this has been studied for mothers of children with disabilities. From an analysis of available data, we found that taking into account many characteristics of the child and family, the presence of an intellectual impairment in children and a low level of education of the mother or the father’s nonemployment appeared to significantly influence the risk of non-employment of the mother. On the other hand, the mothers of children with severe disabilities were more likely than the general population to be not employed (47.2% against 31.4%). But in contrast, mothers of children with CP interviewed during a survey were not more unemployed than mothers in the general population; calling into question the use of CP as a model of disability to study family consequences.The study of the impact of disability is still too rare in Europe and the effects of possible interventions little known. It is important to develop research in this area to improve the lives of these families and provide children with disabilities the best possible environment for its development.The study of the impact of disability is too rare in Europe and little is known on the effects of possible interventions. It is important to develop research in this area to improve the life for these families and to provide for disabled children the best possible environment for its development.

La modélisation structurale consiste à approximer les structures géologiques du sous-sol en un modèle numérique afin d'en visualiser la géométrie et d'y effectuer des calculs d'estimation et de prédiction. L'approche implicite de la modélisation structurale utilise des données de terrain interprétées pour construire une fonction volumétrique sur le domaine d'étude qui représente la géologie. Cette fonction doit honorer les observations, interpoler entre ces dernières, et extrapoler dans les zones sous-échantillonnées tout en respectant les concepts géologiques. Les méthodes actuelles portent cette interpolation soit sur les données, soit sur un maillage. Ensuite, le problème de modélisation est posé selon la discrétisation choisie : par krigeage dual sur les points de donnée ou en définissant un critère de rugosité sur les éléments du maillage. Dans cette thèse, nous proposons une formulation continue de la modélisation structurale par méthodes implicites. Cette dernière consiste à minimiser une somme de fonctionnelles arbitraires. Les contraintes de donnée sont imposées avec des fonctionnelles discrètes, et l'interpolation est contrôlée par des fonctionnelles continues. Cette approche permet de (i) développer des liens entre les méthodes existantes, (ii) suggérer de nouvelles discrétisations d'un même problème de modélisation, et (iii) modifier le problème de modélisation pour mieux honorer certains cas géologiques sans dépendre de la discrétisation. Nous portons également une attention particulière à la gestion des discontinuités telles que les failles et les discordances. Les méthodes existantes nécessitent soit la création de zones volumétriques avec des géométries complexes, soit la génération d'un maillage volumétrique dont les éléments sont conformes aux surfaces de discontinuité. Nous montrons, en explorant des méthodes sans maillage locales et des concepts de réduction de maillage, qu'il est possible d'assurer l'interpolation des structures tout en réduisant les contraintes liées à la gestion des discontinuités. Deux discrétisations de notre problème de minimisation sont suggérées : l'une utilise les moindres carrés glissants avec des critères optiques pour la gestion des discontinuités, et l'autre utilise des fonctions issues de la méthode des éléments finis avec le concept de nœuds fantômes pour les discontinuités. Une étude de sensibilité et une comparaison des deux méthodes sont proposées en 2D, ainsi que quelques exemples en 3D. Les méthodes développées dans cette thèse ont un grand impact en termes d'efficacité numérique et de gestion de cas géologiques complexes. Par exemple, il est montré que notre problème de minimisation au sens large apporte plusieurs solutions pour la gestion de cas de plis sous-échantillonnés et de variations d'épaisseur dans les couches stratigraphiques. D'autres applications sont également présentées tels que la modélisation d'enveloppe de sel et la restauration mécanique. / Implicit structural modeling consists in approximating geological structures into a numerical model for visualization, estimations, and predictions. It uses numerical data interpreted from the field to construct a volumetric function on the domain of study that represents the geology. The function must fit the observations, interpolate in between, and extrapolate where data are missing while honoring the geological concepts. Current methods support this interpolation either with the data themselves or using a mesh. Then, the modeling problem is posed depending on these discretizations: performing a dual kriging between data points or defining a roughness criterion on the mesh elements. In this thesis, we propose a continuous formulation of implicit structural modeling as a minimization of a sum of generic functionals. The data constraints are enforced by discrete functionals, and the interpolation is controlled by continuous functionals. This approach enables to (i) develop links between the existing methods, (ii) suggest new discretizations of the same modeling problem, and (iii) modify the minimization problem to fit specific geological issues without any dependency on the discretization. Another focus of this thesis is the efficient handling of discontinuities, such as faults and unconformities. Existing methods require either to define volumetric zones with complex geometries, or to mesh volumes with conformal elements to the discontinuity surfaces. We show, by investigating local meshless functions and mesh reduction concepts, that it is possible to reduce the constraints related to the discontinuities while performing the interpolation. Two discretizations of the minimization problem are then suggested: one using the moving least squares functions with optic criteria to handle discontinuities, and the other using the finite element method functions with the concept of ghost nodes for the discontinuities. A sensitivity analysis and a comparison study of both methods are performed in 2D, with some examples in 3D. The developed methods in this thesis prove to have a great impact on computational efficiency and on handling complex geological settings. For instance, it is shown that the minimization problem provides the means to manage under-sampled fold structures and thickness variations in the layers. Other applications are also presented such as salt envelope surface modeling and mechanical restoration.

Nous étudions dans ce travail la modélisation de la cristallisation d'une chambre magmatique fermée. Le modèle physique et mathématique couple trois phénomènes élémentaires : solidification, sédimentation et réactions chimiques entre solide et liquide magmatique. L'adimensionnement des équations aux dérivées partielles permet de faire apparaître deux paramètres exprimant les ratios respectifs de la vitesse de solidification sur la vitesse de déplacement et de la cinétique d'échange chimique sur la vitesse de déplacement. La vitesse de déplacement relatif entre le solide et le liquide, la vitesse de solidification et la loi de partage à l'équilibre sont supposées connues ; la loi d'équilibre chimique peut être non linéaire et s'appliquer aux éléments majeurs. Le modèle est écrit pour un constituant chimique. Il est scindé en deux sous modèles, celui de la cristallisation/sédimentation et celui du transport réactif. Le premier est exprimé par une équation de nature hyperbolique et est résolu par un schéma à trois points, le second est résolu par des schémas décentrés. Le code de calcul est écrit en langage Fortran 90, il est ensuite validé par des méthodes théoriques telles que la méthode des courbes caractéristiques, par des calculs analytiques ou par des considérations qualitatives. Les simulations numériques montrent que, pour certaines valeurs des paramètres adimensionnés et pour certaines formes de la loi d'équilibre chimique, la composition chimique de la chambre peut présenter une répartition non uniforme des compositions, en particulier bimodale (deux valeurs des concentrations sont privilégiées), à partir de conditions initiales homogènes. Le degré de cette bimodalité dépend notamment de la forme de la loi de partage à l'équilibre. Ce modèle fournit un cadre intellectuel pour discuter les phénomènes responsables de la variété de compositions des roches magmatiques, notamment dans une même province. Il montre en particulier que le couplage entre trois phénomènes élémentaires suffit à rendre compte de la bimodalité ou plus généralement de l'apparition de discontinuités de composition, sans faire intervenir de phénomène additionnel.