In an optimal visiting problem, we want to control a trajectory that has to pass as close as possible to a collection of target points or regions. We introduce a hybrid control-based approach for the classic problem where the trajectory can switch between a group of discrete states related to the targets of the problem. The model is subsequently adapted to a mean-field game framework, that is when a huge population of agents plays the optimal visiting problem with a controlled dynamics and with costs also depending on the distribution of the population. In particular, we investigate a single continuity equation with possible sinks and sources and the field possibly depending on the mass of the agents. The same problem is also studied on a network framework. More precisely, we study a mean-field game model by proving the existence of a suitable definition of an approximated mean-field equilibrium and then we address the passage to the limit.

The present work seeks to fully investigate, describe and characterize the distinct flow regimes existing within a mixing vessel at various rotational speeds. This investigation is computational in nature and simulates the flow within a baffled tank containing a Rushton turbine of the standard configuration. For a Re based on impeller diameter and blade rotational speed (Re â ¡ Ï ND2/μ) the following flow regimes were identified and investigated in detail: Reverse/reciprocating flows at very low Re (<10); stalled flows at low Re (â 10); laminar pumping flow for higher Re and transitional pumping flow (10 squared < Re <10 to the 4th). For the three Re numbers 1, 10 and 28, it was found that for the higher Re number (28), the flow exhibited the familiar outward pumping action associated with radial impellers under turbulent flow conditions. However, as the Re number decreases, the net radial flow during one impeller revolution was reduced and for the lowest Re number a reciprocating motion with negligible net pumping was observed. In order to elucidate the physical mechanism responsible for the observed flow pattern at low Re, the forces acting on a fluid element in the radial direction were analyzed. Based on this analysis, a simplified quasi-analytic model of the flow was developed that gives a satisfactory qualitative, as well as quantitative representation of the flow at very low Re. Investigation of the transitional flow regime (Re â 3000) includes a compilation and characterization of ensemble and turbulent quantities such as the Reynolds stress components, dissipation length η and time scales Ï , as well a detailed investigation of the near-impeller flow and trailing vortex. Calculation and compilation of all terms in the turbulent kinetic energy transport equation was performed (including generation and the illusive turbulent pressure work). Specifically, the most important transport mechanism was turbulent convection/diffusion from the impeller disk-plane/trailing vortex region. Mean flow transport of turbulent kinetic energy was primarily towards the impeller disk-plane and radially outward from the trailing vortex region. The turbulent pressure work was found to partially counteract turbulent convection. Turbulent dissipation followed by turbulent viscous work were found to be the least important mechanism responsible for turbulent transport with both terms being maximized within the vortex region and at the disk-plane down-stream from the vortices. / Ph. D.

Les pathologies cardiovasculaires sont la première cause des décès dans le monde. Il est donc vital de les étudier afin d’en comprendre les mécanismes et pouvoir prévenir et traiter plus efficacement ces maladies. Cela passe donc par la compréhension de l’anatomie, de la structure et du mouvement du coeur. Dans le cadre de cette thèse, nous nous sommes intéressés dans un premier temps au modèle de Gabarit Déformable Élastique qui permet d’extraire l’anatomie et le mouvement cardiaques. Le Gabarit Déformable Élastique consiste à représenter le myocarde par un modèle de forme a priori donné que l’on déforme élastiquement pour l’adapter à la forme spécifique du coeur du patient Dans le premier chapitre de cette thèse, nous utilisons une méthode de perturbation singulière permettant la segmentation avec précision de l’image. Nous avons démontré que si l’on faisait tendre vers 0 les coefficients de l’élasticité, le modèle mathématique convergeait vers une solution permettant la segmentation. Dans le cadre d’une formulation au sens des moindres carrés il est nécessaire de disposer d’une méthode numérique performante pour résoudre l’équation du transport au sens des moindres carrés. La méthode des éléments finis pour traiter les phénomènes de transport ne permet pas d’avoir un principe du maximum faible, sauf si l’opérateur aux dérivées partielles en temps est séparé de l’opérateur aux dérivées partielles en espace. Dans le chapitre 2 de la thèse nous considérons une formulation au sens des moindres carrés espace-temps et nous proposons de résoudre un problème sous contraintes afin de récupérer un principe du maximum discret. Le dernier objectif de la thèse est le suivi dynamique d’images cardiaques ou la reconstruction anatomique du coeur à partir de coupes 2D dans le plan orthogonal au grand axe du cœur. La méthode mathématique que nous utilisons pour cela est le transport optimal. Dans le chapitre 3 nous analysons les performances de l’algorithme proposé par Peyré pour calculer le transport optimal de nos images. La résolution numérique du transport optimal est un problème difficile et couteux en temps de calcul. C’est pourquoi nous proposons une méthode adaptative pour le calcul de l’opérateur proximal de la fonction à minimiser permettant de diviser par quatre le nombre nécessaire des itérations pour que l’algorithme converge. / Cardiovascular disease are the leading cause of death worldwide. It is therefore vital to study them in order to understand the mechanisms and to prevent and treat these diseases more effectively. Therefore it requires an understanding of anatomy, structure and motion of the heart. In this thesis, we are interested in a first time at the Deformable Elastic Template model which can extract the cardiac anatomy and movement. The Elastic Deformable template is to represent the myocardium by a shape model a priori given that it elastically deforms to fit the specific shape of the patient’s heart. In the first chapter of this thesis, we use a singular perturbation method for accurately segmenting the image. We have demonstrated that if we did tend to 0 the coefficients of elasticity, the mathematical model converge to a solution to the problem of segmentation. As part of a formulation to the least squares sense it is necessary to have an efficient numerical method for solving the transport equation in the least squares sense. The finite element method to treat transport phenomena can not have a weak maximum principle, unless the operator of partial time is separated from the operator of partial space.In Chapter 2 of the thesis, we consider a least squares formulation of space-time and we propose to solve the problem constraints to recover a discrete maximum principle. The final objective of this thesis is the dynamic monitoring of cardiac images or anatomical reconstruction of the heart from 2D slices orthogonal to the long axis of the heart level. The mathematical method we use for this is the optimal transport. In Chapter 3 we analyze the performance of the algorithm proposed by Peyré to calculate the optimal transport of our images. The numerical resolution of optimal transport is a difficult and costly in computation time problem. That is why we propose an adaptive method for determining the proximity operator of the function to be minimized to divide by four the number of iterations required for the algorithm converges.

Les méthodes particulaires sont des méthodes numériques adaptées à la résolution d'équations de conservation. Leur principe consiste à introduire des particules ``numériques'' conservant localement l'inconnue sur un petit volume, puis à les transporter le long de leur trajectoire. Lorsqu'un terme source est présent dans les équations, l'évolution de la solution le long des caractéristiques est prise en compte par une intéraction entre les particules. Ces méthodes possèdent de bonnes propriétés de conservation et ne sont pas soumises aux conditions habituelles de CFL qui peuvent être contraignantes pour les méthodes Eulériennes. Cependant, une contrainte de recouvrement entre les particules doit être satisfaite pour vérifier des propriétés de convergence de la méthode. Pour satisfaire cette condition de recouvrement, un remaillage périodique des particules est souvent utilisé. Elle consiste à recréer régulièrement de nouvelles particules uniformément réparties, à partir de celles ayant été advectées à l'itération précédente. Quand cette étape de remaillage est effectuée à chaque pas de temps, l'analyse numérique de ces méthodes particulaires remaillées nécessite d'être reconsidérée, ce qui représente l'objectif de ces travaux de thèse. Pour mener à bien cette analyse, nous nous basons sur une analogie entre méthodes particulaires avec remaillage et schémas de grille. Nous montrons que pour des grands pas de temps les schémas numériques obtenus souffrent d'une perte de précision. Nous proposons des méthodes de correction, assurant la consistance des schémas en tout point de grille, le pas de temps étant contraint par une condition sur le gradient du champ de vitesse. Cette méthode est construite en dimension un. Des techniques de limitation sont aussi introduites de manière à remailler les particules sans créer d'oscillations en présence de fortes variations de la solution. Enfin, ces méthodes sont généralisées aux dimensions plus grandes que un en s'inspirant du principe de splitting d'opérateurs. Les applications numériques présentées dans cette thèse concernent la résolution de l'équation de transport sous forme conservative en dimension un à trois, dans des régimes linéaires ou non-linéaires. / Particle methods are numerical methods designed to solve advection dominated conservation equations. Their principle is to introduce ``numerical'' particles that concentrate the unknown locally on a small volume, and to transport them along their trajectories. These methods have good conservation properties and are not subject to the usual CFL conditions that can be binding for the Eulerian methods. However, an overlap condition must be satisfied between the particles to ensure convergence properties of the method. To satisfy this condition, a periodic remeshing of the particles is often used. New particles uniformly distributed are created on a regular mesh. When this remeshing step is performed at every time step, numerical analysis of particle methods needs to be revisited. This is the purpose of this thesis. To carry out this analysis, we rely on an analogy between remeshed particle methods and grid schemes. We show that for large time step the numerical schemes have a loss of accuracy. We propose correction methods wich ensure consistency at any grid point, provided the time step satisfies a condition based on the gradient of the velocity field. Limitation techniques are also introduced to remesh particles without creating any oscillations in the presence of strong variations of the solution. Finally, these methods are generalized to dimensions greater than one. Numerical example on various transport equations are given to illustrate the benefit of the proposed algorithms.

Le développement de la technologie des composants CMOS ultimes à grille ultra-courte (L < 20 nm) se heurte à de nombreuses difficultés technologiques, mais également à des limites thermiques qui perturbent notablement les règles de mise à l'échelle communément employées jusqu'à présent. Les fortes densités de courant obtenues dans des zones actives aussi réduites génèrent un important échauffement local (par effet Joule), lié à l'émission de phonons par les porteurs chauds, qui peut conduire à des réductions très sensibles des performances, voire à des défaillances. Ce phénomène est identifié comme un des plus critiques pour la poursuite de l'augmentation de la densité d'intégration des circuits. Cela est particulièrement crucial dans les technologies SOI (silicium sur isolant), où la présence de l'isolant enterré constitue un frein à l'évacuation de la chaleur. À l'échelle nanométrique, l'étude théorique de ces phénomènes d'échauffement n'est plus possible par des modèles macroscopiques (coefficient de diffusion de la chaleur) mais nécessite une description microscopique détaillée des transferts de chaleur qui sont localement hors d’équilibre. Il s'agit donc de modéliser de façon appropriée, non seulement le transport électronique et la génération de phonons, mais aussi le transport de phonons hors équilibre et les interactions phonons-phonons et électrons-phonons.Le formalisme de l’équation de transport de Boltzmann (BTE) est très bien adapté à l'étude de ce problème. En effet, il est largement utilisé depuis des années pour l'étude du transport des particules chargées dans les composants semi-conducteurs. Ce formalisme est beaucoup moins standard pour étudier le transport des phonons. Une des problématiques de ce travail concerne le couplage de la résolution de la BTE des phonons avec celle des électrons.Ce travail de thèse a développé un algorithme de calcul du transport de phonons par résolution directe de la BTE des phonons. Cet algorithme de transport de phonon a été couplé au transport électronique simulé grâce au logiciel "MONACO" basé sur une résolution statistique (ou Monte Carlo) de la BTE. Finalement, ce nouveau simulateur électrothermique a été utilisé pour étudier les effets d’auto échauffement dans des nano-transistors. L’intérêt principal de ces travaux est de permettre une analyse du transport electro-thermique au-delà d’une approche macroscopique (respectivement formalisme de Fourier pour la thermique et dérive-diffusion pour le courant). En effet, il donne accès aux distributions de phonons dans le dispositif et pour chaque mode de phonon. En particulier, ce simulateur apporte une meilleure compréhension des effets des électrons chauds sur les points chauds et leur relaxation dans les accès. / The ultra-short gate (LG < 20 nm) CMOS components (Complementary Metal-Oxide-Semiconductor) face thermal limitations due to significant local heating induced by phonon emission by hot carriers in active regions of reduced size. This phenomenon, called self-heating effect, is identified as one of the most critical for the continuous increase in the integration density of circuits. This is especially crucial in SOI technology (silicon on insulator), where the presence of the buried insulator hinders the dissipation of heat.At the nanoscale, the theoretical study of these heating phenomena, which cannot be led using the macroscopic models (heat diffusion coefficient), requires a detailed microscopic description of heat transfers that are locally non-equilibrium. It is therefore appropriate to model, not only the electron transport and the phonon generation, but also the phonon transport and the phonon-phonon and electron-phonon interactions. The formalism of the Boltzmann transport equation (BTE) is very suitable to study this problem. In fact, it is widely used for years to study the transport of charged particles in semiconductor components. This formalism is much less standard to study the transport of phonons. One of the problems of this work concerns the coupling of the phonon BTE with the electron transport.In this context, wse have developed an algorithm to calculate the transport of phonons by the direct solution of the phonon BTE. This algorithm of phonon transport was coupled with the electron transport simulated by the simulator "MONACO" based on a statistical (Monte Carlo) solution of the BTE. Finally, this new electro-thermal simulator was used to study the self-heating effects in nano-transistors. The main interest of this work is to provide an analysis of electro-thermal transport beyond a macroscopic approach (Fourier formalism for thermal transport and the drift-diffusion approach for electric current, respectively). Indeed, it provides access to the distributions of phonons in the device for each phonon mode. In particular, the simulator provides a better understanding of the hot electron effects at the hot spots and of the electron relaxation in the access.

Nesta dissertação, são apresentados os seguintes modelos matemáticos de transporte de nêutrons: a equação linearizada de Boltzmann e a equação da difusão de nêutrons monoenergéticos em meios não-multiplicativos. Com o objetivo de determinar o período fluxo escalar de nêutrons, é descrito um método espectronodal que gera soluções numéricas para o problema de difusão em geometria planar de fonte fixa, que são livres de erros de truncamento espacial, e que conjugado com uma técnica de reconstrução espacial intranodal gera o perfil detalhado da solução. A fim de obter o valor aproximado do fluxo angular de nêutrons em um determinado ponto do domínio e em uma determinada direção de migração, descreve-se também um método de reconstrução angular baseado na solução analítica da equação unidimensional de transporte de nêutrons monoenergéticos com espalhamento linearmente anisotrópico com aproximação sintética de difusão nos termos de fonte por espalhamento. O código computacional desenvolvido nesta dissertação foi implementado na plataforma livre Scilab, e para ilustrar a eficiência do código criado,resultados numéricos obtidos para três problemas-modelos são apresentados / We describe a method to determine the neutron scalar flux in a slab using monoenergetic diffusion model. To achieve this goal we used three ingredients in the computational code that we developed on the Scilab platform: (i) a spectral nodal method that generates numerical solution for the one-speed slab-geometry fixed-source difusion problem with no spatial truncation errors; (ii) a spatial reconstruction scheme to yield detailed proile of the coarse-mesh solution; and (iii) an angular reconstruction scheme to yield approximately the neutron angular flux profile within the slab. The angular reconstruction scheme is based on the analytical solution of the neutron transport equation in slab geometry with linearly anisotropic scattering and diffusion approximation for the scattering source terms. Numerical results are given to illustrate the efficiency of the offered code

Thickened tailings production and disposal continue to grow in importance in the mining industry. In particular, the transport of oil sands tailings is of interest in this study. These tailings must be in a homogeneous state (non-segregating) during pipeline flow and subsequent discharge. Tailings are often transported in an open channel or flume. Slurries containing both clay and coarse sand particles typically exhibit non-Newtonian rheological behaviour. The prediction of the flow behaviour of these slurries is complicated by the limited research activity in this area. As a result, the underlying mechanisms of solids transport in these slurries are not well understood. To address this deficiency, experimental studies were conducted with kaolin clay slurries containing coarse sand in an open circular channel.<p> A numerical model has been developed to predict the behaviour of coarse solid particles in laminar, open channel, non-Newtonian flows. The model involves the simultaneous solution of the Navier-Stokes equations and a scalar concentration equation describing the behaviour of coarse particles within the flow. The model uses the theory of shear-induced particle diffusion (Phillips et al., 1992) to provide a number of relationships to describe the diffusive flux of coarse particles within laminar flows. A sedimentation flux has been developed and incorporated into the Phillips et al. (1992) model to account for gravitational flux of particles within the flow. Previous researchers (Gillies et al., 1999) have shown that this is a significant mechanism of particle migration.<p> The momentum and concentration partial differential equations have been solved numerically by applying the finite volume method. The differential equations are non-linear, stiff and tightly coupled which requires a novel means of analysis. Specific no-flux, no-slip and no-shear boundary conditions have been applied to the channel walls and free surface to produce simulated velocity and concentration distributions. The results show that the model is capable of predicting coarse particle settling in laminar, non-Newtonian, open channel flows. The results of the numerical simulations have been compared to the experimental results obtained in this study, as well as the experimental results of previous studies in the literature.

The mechanisms of thermal transport in defect-free silicon nanostructures are examined using a combination of lattice dynamics (LD) calculations and the Boltzmann transport equation (BTE). To begin, the thermal conductivity reduction in thin films is examined using a hierarchical method that first predicts phonon transport properties using LD calculations, and then solves the phonon BTE using the lattice Boltzmann method. This approach, which considers all of the phonons in the first Brillouin-zone, is used to assess the suitability of common assumptions used to reduce the computational effort. Specifically, we assess the validity of: (i) neglecting the contributions of optical modes, (ii) the isotropic approximation, (iii) assuming an averaged bulk mean-free path (i.e., the Gray approximation), and (iv) using the Matthiessen rule to combine the effect of different scattering mechanisms. Because the frequency-dependent contributions to thermal conductivity change as the film thickness is reduced, assumptions that are valid for bulk are not necessarily valid for thin films. Using knowledge gained from this study, an analytical model for the length-dependence of thin film thermal conductivity is presented and compared to the predictions of the LD-based calculations. The model contains no fitting parameters and only requires the bulk lattice constant, bulk thermal conductivity, and an acoustic phonon speed as inputs. By including the mode-dependence of the phonon lifetimes resulting from phonon-phonon and phonon-boundary scattering, the model predictions capture the approach to the bulk thermal conductivity better than predictions made using Gray models based on a single lifetime. Both the model and the LD-based method are used to assess a procedure commonly used to extract bulk thermal conductivities from length-dependent molecular dynamics simulation data. Because the mode-dependence of thermal conductivity is not included in the derivation of this extrapolation procedure, using it can result in significant error. Finally, phonon transport across a silicon/vacuum-gap/silicon structure is modelled using lattice dynamics and Landauer theory. The phonons transmit thermal energy across the vacuum gap via atomic interactions between the leads. Because the incident phonons do not encounter a classically impenetrable potential barrier, this mechanism is not a tunneling phenomenon. The heat flux due to phonon transport can be 4 orders of magnitude larger than that due to photon transport predicted from near-field radiation theory.

Thickened tailings production and disposal continue to grow in importance in the mining industry. In particular, the transport of oil sands tailings is of interest in this study. These tailings must be in a homogeneous state (non-segregating) during pipeline flow and subsequent discharge. Tailings are often transported in an open channel or flume. Slurries containing both clay and coarse sand particles typically exhibit non-Newtonian rheological behaviour. The prediction of the flow behaviour of these slurries is complicated by the limited research activity in this area. As a result, the underlying mechanisms of solids transport in these slurries are not well understood. To address this deficiency, experimental studies were conducted with kaolin clay slurries containing coarse sand in an open circular channel.<p> A numerical model has been developed to predict the behaviour of coarse solid particles in laminar, open channel, non-Newtonian flows. The model involves the simultaneous solution of the Navier-Stokes equations and a scalar concentration equation describing the behaviour of coarse particles within the flow. The model uses the theory of shear-induced particle diffusion (Phillips et al., 1992) to provide a number of relationships to describe the diffusive flux of coarse particles within laminar flows. A sedimentation flux has been developed and incorporated into the Phillips et al. (1992) model to account for gravitational flux of particles within the flow. Previous researchers (Gillies et al., 1999) have shown that this is a significant mechanism of particle migration.<p> The momentum and concentration partial differential equations have been solved numerically by applying the finite volume method. The differential equations are non-linear, stiff and tightly coupled which requires a novel means of analysis. Specific no-flux, no-slip and no-shear boundary conditions have been applied to the channel walls and free surface to produce simulated velocity and concentration distributions. The results show that the model is capable of predicting coarse particle settling in laminar, non-Newtonian, open channel flows. The results of the numerical simulations have been compared to the experimental results obtained in this study, as well as the experimental results of previous studies in the literature.

Gallium nitride (GaN) has been targeted for use in high power (>30 W/mm) and high frequency (>160 GHz) application due to its wide band gap and its large break down field. One of the most significant advances in GaN devices has evolved from the AlGaN/GaN high electron mobility transistor (HEMT). As a result of the large power densities being applied to these devices there can develop intense hot spots near areas of highest electric field. The hot spot phenomenon has been linked to a decrease in device reliability through a range of degradation mechanisms. In order to minimize the effect that hot spot temperatures have on device reliability a detailed understanding of relevant transport mechanisms must be developed. This study focuses on two main aspects of phonon transport within GaN devices. The first area of focus was to establish an understanding of phonon relaxation times within bulk GaN. These relaxation times were calculated from an application of Fermi's Golden Rule and explicitly conserve energy and crystal momentum. This analysis gives insight into the details behind the macroscopic thermal conductivity parameter. Once relaxation times for GaN were established a multiscale phonon transport modeling methodology was developed that allowed the Boltzmann Transport Equation to be coupled to the energy equation. This coupling overcomes some computational limits and allows for nanoscale phenomena to be resolved within a macroscopic domain. Results of the transport modeling were focused on benchmarking the coupling method as well as calculating the temperature distribution within an operating 6 finger HEMT.