Ce travail concerne la modélisation du mouvement des piétons via l’approche non régulière du contact dynamique des solides rigides et déformables. Une reformulation de cette approche est proposée en accord avec le formalisme de M.Frémond et celui de J.J.Moreau. L’approche proposée est basée sur la notion de percussion qui est l’intégrale de la force de contact au cours de la durée de la collision. Contrairement aux modèles classiques d’éléments discrets, il est supposé que les percussions ne peuvent être exprimées qu’en fonction de la vitesse avant l’impact. Cette hypothèse est vérifiée pour des lois de comportement classiques de la collision. Les équations de mouvement sont ensuite reformulées en tenant compte de multiples collisions simultanées. L’existence et l’unicité de la solution du nouveau modèle sont discutées en fonction de la régularité des forces (densité de Lebesgue apparaissant au cours de l’évolution régulière du système) et la régularité des percussions (Dirac-densité décrivant la collision). A la lumière des principes de la thermodynamique, une condition sur la percussion interne assurant que la collision est thermodynamiquement admissible, est établi. L’application aux collisions de disques rigides et à l’écoulement dans un sablier en forme d'entonnoir est présentée. L’approche est étendue au mouvement de la foule, en effet ; la circulation des piétons à travers les goulets d’étranglement est étudiée. Une analyse de sensibilité est effectuée pour étudier l’effet des paramètres d’un modèle de mouvement de foule discret 2D sur la nature des collisions et des temps d’évacuation des piétons. Nous avons identifié les paramètres qui régissent une collision de type piéton-piéton et étudié leurs effets sur le temps d'évacuation. Puis une expérience d’évacuation d’une salle avec une sortie de goulot d’étranglement est introduite et sa configuration est utilisée pour les simulations numériques. La question de l’estimation des forces de contact et de la pression générée dans une foule en mouvement est abordée à la fois d’un point de vue discret (un piéton est assimilé à un disque rigide) et continue (la foule est considérée comme un solide déformable). Une comparaison entre le modèle microscopique du second ordre (modèle discret 2D) et l’approche continue est présentée. Les forces de contact sont rigoureusement définies en tenant compte des contacts multiples et simultanés et le non chevauchement entre piétons. Nous montrons que pour une foule dense les percussions (saut de la quantité de mouvement correspondant au contact instantané) deviennent des forces de contact. Pour l’approche continue, la pression est calculée en fonction des contraintes volumiques et surfaciques. Et tenant compte les interactions non locales entre les piétons. Afin de rendre l’approche plus efficace, nous avons modélisé chaque piéton par un solide déformable, le cas unidimensionnel est étudié, une comparaison avec le cas discret est présentée pour un exemple d’écrasement d’une chaîne de piétons dans un obstacle fixe. La solution analytique des équations de contact est développée ce qui permet une calibration de paramètres du modèle et une étude asymptotique des solutions. La théorie non-régulière de la dynamique de solides déformables permet de calculer la vitesse réelle de la foule en tant qu’un milieu continu en tenant compte des interactions avec l’environnement et de la vitesse souhaitée. Une représentation macroscopique donnée par un problème couplé d’équations hyperbolique et elliptique. Une équation hyperbolique décrivant l’évolution de la densité de la foule dont la vitesse est calculée une équation elliptique, celle de l’évolution d’un solide déformable. Un résultat d’existence et unicité est développé concernant l’existence et l’unicité de la solution du problème couplé et la stabilité par rapport à la condition initiale et les conditions aux limites / This work concerns the modeling of pedestrian movement inspired by the non-smooth dynamics approach for the rigid and deformable solids. Firstly, a reformulation of the non-smooth approaches of M.Frémond and J.J.Moreau for rigid body dynamics is developed. The proposed theory relies on the notion of percussion which is the integral of the contact force during the duration of the collision. Contrary to classical discrete element models, it is here assumed that percussions can be only expressed as a function of the velocity before the impact. This assumption is checked for the usual mechanical constitutive laws for collisions. Motion equations are then reformulated taking into account simultaneous collisions of solids. The existence and uniqueness of the solution of the proposed model are discussed according to the regularity of both the forces (Lebesgue-density occurring during the regular evolution of the system) and the percussions (Dirac-density describing the collision). A condition on the internal percussion assuring that the collision is thermodynamically admissible is established. An application to the collision of rigid disks and the flow in a funnel-shaped hourglass is presented. The approach is extended to crowd motion, indeed; the circulation of pedestrians through the bottlenecks is studied and deals with to optimize evacuation and improve the design of pedestrian facilities. A sensitivity analysis is performed to study the effect of the parameters of a 2D discrete crowd movement model on the nature of pedestrian’s collision and on evacuation times. The question of estimation of contact forces and the pressure generated in a moving crowd is approached both from a discrete and continues point of view. A comparison between the second-order microscopic model (2D discrete model) and the continues approaches is presented. Contact forces are rigorously defined taking into account multiple, simultaneous contact and the non-overlapping condition between pedestrians. We show that for a dense crowd the percussions (moment umjump corresponding to instantaneous contact) become contact forces. For continuous approach, the pressure is calculated according to volume and surface constraints. This approach makes it possible to retain an admissible right-velocity (after impact), including both the non local interactions (at a distance interactions) between non neighbor pedestrians and the choice of displacement strategy of each pedestrian. Finally, two applications are presented : a one-dimensional simulation of an aligned pedestrian chain crashing into an obstacle, and a two-dimensional simulation corresponding to the evacuation of a room. In order to make the approach more efficient, we modeled each pedestrian with a deformable solid, the unidimensional case is studied a comparison with the discreet case is presented that corresponding to a crash of a pedestrian chain in a fixed obstacle is treated. The analytical solution of contact equations is developed for both approaches. This allows to calibrate the model parameters and offers an asymptotic study of the solutions. The non-smooth theory of deformable solids makes it possible to calculate the current velocity of the crowd as a continuous medium taking into account the interactions with the environment and their desired velocity. a macroscopic representation is developed through Hyperbolic – Elliptic Equations. indded;the crowd is described by its density whose evolution is given by a non local balance law. the current velocity involved in the equation is given by the collision equation of a deformable solid with a rigid plane. Firstly, we prove the well posedness of balance laws with a non smooth ux and function source in bounded domains, the existence of a weak entropic solution, it’s uniqueness and stability with respect to the initial datum and of the boundary datum. an application to crowdmodeling is presented

Orientadores: Yuri Dimitrov Bozhkov, Antonio Carlos Gilli Martins / Tese (doutorado) - Universidade Estadual de Campinas, Instituto de Matematica, Estatistica e Computação Cientifica / Made available in DSpace on 2018-09-24T19:39:04Z (GMT). No. of bitstreams: 1 Freire_IgorLeite_D.pdf: 977261 bytes, checksum: b8ba44493aeac3de0d37cdfff2fc581b (MD5) Previous issue date: 2008 / Resumo: Neste trabalho provamos um teorema que faz a classificacão completa dos grupos de simetrias de Lie da equação semilinear de Kohn - Laplace no grupo de Heisenberg tridimensional. Uma vez que tal equação possui estrutura variacional, determinamos quais são suas simetrias de Noether e a partir delas estabelecemos suas respectivas leis de conservação / Abstract: In this work, we carry out a complete group classification of Lie point symmetries of semilinear Kohn - Laplace equations on the three-dimensional Heisenberg group. Since this equation has variational structure, we determine which of its symmetries are Noether's symmetries. Then we establish their respectives conservation laws / Doutorado / Matematica Aplicada / Doutor em Matemática Aplicada

Nesta dissertação de Mestrado será apresentada uma poderosa técnica numérica, conhecida como método Level Set, capaz de simular e analisar movimentos de curvas em diferentes cenários físicos. Tal método - formulado por Osher e Sethian [1] - está sedimentado na seguinte idéia: representar uma determinada curva (ou superfície) Γ como a curva de nível zero (zero level set) de uma função Φ de maior dimensão (denominada função Level Set). A equação diferencial do tipo Hamilton-Jacobi que descreve a evolução da função Level Set é discretizada através da utilização de acurados esquemas hiperbólicos e, como resultado de tal acurácia, obtém-se uma formulação numérica capaz de tratar eficazmente mudanças topológicas e/ou descontinuidades que, eventualmente, podem surgir no decorrer da propagação da curva (ou superfície) de nível zero. Em virtude da eficácia e versatilidade do método Level Set, esta técnica numérica está sendo amplamente aplicada à diversas áreas científicas, incluindo mecânica dos fluidos, processamento de imagens e visão computacional, crescimento de cristais, geometria computacional e ciência dos materiais. Particularmente, o propósito deste trabalho equivale ao estudo dos fundamentos do método Level Set e, por fim, visa-se aplicar tal modelo numérico à problemas existentes na área de crescimento de cristais. [1] S. Osher and J. A. Sethian, Fronts propagating with curvature-dependent speed: Algorithms based on Hamilton-Jacobi formulations, J. Comp. Phys., 79:12, 1988. / In this dissertation, we present a powerful numerical technique known as Level Set Method for computing and analyzing moving fronts in different physical settings. The method -formulated by Osher and Sethian [1] - is based on the following idea: a curve (or surface) is embedded as the zero level set of a higher-dimensional function Φ (called level set function). Then, we can link the evolution of this function Φ to the propagation of the curve itself through a time-dependent initial value problem. At any time, the curve is given by the zero level set of the time-dependent level set function Φ. The evolution of the level set function Φ is described by a Hamilton-Jacobi type partial differential equation, which can be discretised by the use of accurate methods for hyperbolic equations. As a result, the Level Set Method is able to track complex curves that can develop large spikes, sharp corners or change its topology as they evolve. Because of its versatility and efficacy, this numerical technique has found applications in a large number of areas, including fluid mechanics, image processing and computer vision, crystal growth, computational geometry and materials science. Particularly, the aim of this dissertation has been to understand the fundamentals of Level Set Method and its final goal is compute the motion of bondaries in crystal growth using this numerical model. [1] S. Osher and J. A. Sethian, Fronts propagating with curvature-dependent speed: Algorithms based on Hamilton-Jacobi formulations, J. Comp. Phys., 79:12, 1988.

Nesta dissertação de Mestrado será apresentada uma poderosa técnica numérica, conhecida como método Level Set, capaz de simular e analisar movimentos de curvas em diferentes cenários físicos. Tal método - formulado por Osher e Sethian [1] - está sedimentado na seguinte idéia: representar uma determinada curva (ou superfície) Γ como a curva de nível zero (zero level set) de uma função Φ de maior dimensão (denominada função Level Set). A equação diferencial do tipo Hamilton-Jacobi que descreve a evolução da função Level Set é discretizada através da utilização de acurados esquemas hiperbólicos e, como resultado de tal acurácia, obtém-se uma formulação numérica capaz de tratar eficazmente mudanças topológicas e/ou descontinuidades que, eventualmente, podem surgir no decorrer da propagação da curva (ou superfície) de nível zero. Em virtude da eficácia e versatilidade do método Level Set, esta técnica numérica está sendo amplamente aplicada à diversas áreas científicas, incluindo mecânica dos fluidos, processamento de imagens e visão computacional, crescimento de cristais, geometria computacional e ciência dos materiais. Particularmente, o propósito deste trabalho equivale ao estudo dos fundamentos do método Level Set e, por fim, visa-se aplicar tal modelo numérico à problemas existentes na área de crescimento de cristais. [1] S. Osher and J. A. Sethian, Fronts propagating with curvature-dependent speed: Algorithms based on Hamilton-Jacobi formulations, J. Comp. Phys., 79:12, 1988. / In this dissertation, we present a powerful numerical technique known as Level Set Method for computing and analyzing moving fronts in different physical settings. The method -formulated by Osher and Sethian [1] - is based on the following idea: a curve (or surface) is embedded as the zero level set of a higher-dimensional function Φ (called level set function). Then, we can link the evolution of this function Φ to the propagation of the curve itself through a time-dependent initial value problem. At any time, the curve is given by the zero level set of the time-dependent level set function Φ. The evolution of the level set function Φ is described by a Hamilton-Jacobi type partial differential equation, which can be discretised by the use of accurate methods for hyperbolic equations. As a result, the Level Set Method is able to track complex curves that can develop large spikes, sharp corners or change its topology as they evolve. Because of its versatility and efficacy, this numerical technique has found applications in a large number of areas, including fluid mechanics, image processing and computer vision, crystal growth, computational geometry and materials science. Particularly, the aim of this dissertation has been to understand the fundamentals of Level Set Method and its final goal is compute the motion of bondaries in crystal growth using this numerical model. [1] S. Osher and J. A. Sethian, Fronts propagating with curvature-dependent speed: Algorithms based on Hamilton-Jacobi formulations, J. Comp. Phys., 79:12, 1988.

Conservation laws play an important role in science. The aim of this thesis is to provide an overview and develop new methods for constructing conservation laws using Lie group theory. The derivation of conservation laws for invariant variational problems is based on Noether’s theorem. It is shown that the use of Lie-Bäcklund transformation groups allows one to reduce the number of basic conserved quantities for differential equations obtained by Noether’s theorem and construct a basis of conservation laws. Several examples on constructing a basis for some well-known equations are provided. Moreover, this approach allows one to obtain new conservation laws even for equations without Lagrangians. A formal Lagrangian can be introduced and used for computing nonlocal conservation laws. For self-adjoint or quasi-self-adjoint equations nonlocal conservation laws can be transformed into local conservation laws. One of the fields of applications of this approach is electromagnetic theory, namely, nonlocal conservation laws are obtained for the generalized Maxwell-Dirac equations. The theory is also applied to the nonlinear magma equation and its nonlocal conservation laws are computed.

Conservation laws play an important role in science. The aim of this thesis is to provide an overview and develop new methods for constructing conservation laws using Lie group theory. The derivation of conservation laws for invariant variational problems is based on Noether’s theorem. It is shown that the use of Lie-Bäcklund transformation groups allows one to reduce the number of basic conserved quantities for differential equations obtained by Noether’s theorem and construct a basis of conservation laws. Several examples on constructing a basis for some well-known equations are provided. Moreover, this approach allows one to obtain new conservation laws even for equations without Lagrangians. A formal Lagrangian can be introduced and used for computing nonlocal conservation laws. For self-adjoint or quasi-self-adjoint equations nonlocal conservation laws can be transformed into local conservation laws. One of the fields of applications of this approach is electromagnetic theory, namely, nonlocal conservation laws are obtained for the generalized Maxwell-Dirac equations. The theory is also applied to the nonlinear magma equation and its nonlocal conservation laws are computed. / <p>Thesis for the degree of Doctor of Philosophy</p>

The Vlasov-Poisson system is most commonly used to model the movement of charged particles in a plasma or of stars in a galaxy. It consists of a kinetic equation known as the Vlasov equation coupled with a force determined by the Poisson equation. The system in Euclidean space is well-known and has been extensively studied under various assumptions. In this paper, we derive the Vlasov-Poisson equations assuming the particles exist only on the 2-sphere, then take an in-depth look at particles which initially lie along a great circle of the sphere. We show that any great circle is an invariant set of the equations of motion and prove that the total energy, number of particles, and entropy of the system are conserved for circular initial distributions. / Graduate

Submitted by Johnny Rodrigues (johnnyrodrigues@ufcg.edu.br) on 2018-08-10T18:01:20Z No. of bitstreams: 1 KEYTT AMARAL DA SILVA - DISSERTAÇÃO PPGMAT 2015..pdf: 1966719 bytes, checksum: d55ff8700252c9540c54209c808e22a3 (MD5) / Made available in DSpace on 2018-08-10T18:01:20Z (GMT). No. of bitstreams: 1 KEYTT AMARAL DA SILVA - DISSERTAÇÃO PPGMAT 2015..pdf: 1966719 bytes, checksum: d55ff8700252c9540c54209c808e22a3 (MD5) Previous issue date: 2015-08 / Neste trabalho apresentamos a construção detalhada da solução do Problema de Riemann associado à um sistema de leis de conservação de um problema não estritamente hiperbólico, proveniente da modelagem matemática de um escoamento unidimensional bifásico num meio poroso em que as fases são óleo e água com polímero dissolvido, para dados iniciais arbitrários no espaço de estados. A construção da solução do sistema é baseada na solução da equação de Buckley−Leverett para cada nível de concentração constante de polímero e nas curvas integrais de uma campo característico linearmente degenerado que dá origem as chamadas ondas de contato. / We present the detailed construction of the Riemann problem solution associate to a system of conservation laws of a non−strictly hyperbolic problem, from mathematical modeling of a one-dimensional two-flow in a porous medium filled by oil and water with dissolved polymer, for arbitrary initial data in the state space. The construction of the system solution is based on the solution Buckley−Leverett equation for each level constant polymer concentration and on the integral curves of a linearly degenerated field characteristic that gives rise to the so-called contact waves.

Les dynamiques des systèmes modélisés par des équations aux dérivées partielles (EDPs) en dimension infinie sont largement liées aux réseaux physiques. La synthèse de la commande et l'analyse de la stabilité de ces systèmes sont étudiées dans cette thèse. Les systèmes singulièrement perturbés, contenant des échelles de temps multiples sont naturels dans les systèmes physiques avec des petits paramètres parasitaires, généralement de petites constantes de temps, les masses, les inductances, les moments d'inertie. La théorie des perturbations singulières a été introduite pour le contrôle à la fin des années $1960$, son assimilation dans la théorie du contrôle s'est rapidement développée et est devenue un outil majeur pour l'analyse et la synthèse de la commande des systèmes. Les perturbations singulières sont une façon de négliger la transition rapide, en la considérant dans une échelle de temps rapide séparée. Ce travail de thèse se concentre sur les systèmes hyperboliques linéaires avec des échelles de temps multiples modélisées par un petit paramètre de perturbation. Tout d'abord, nous étudions une classe de systèmes hyperboliques linéaires singulièrement perturbés. Comme le système contient deux échelles de temps, en mettant le paramètre de la perturbation à zéro, deux sous-systèmes, le système réduit et la couche limite, sont formellement calculés. La stabilité du système complet de lois de conservation implique la stabilité des deux sous-systèmes. En revanche un contre-exemple est utilisé pour illustrer que la stabilité des deux sous-systèmes ne suffit pas à garantir la stabilité du système complet. Cela montre une grande différence avec ce qui est bien connu pour les systèmes linéaires en dimension finie modélisés par des équations aux dérivées ordinaires (EDO). De plus, sous certaines conditions, l'approximation de Tikhonov est obtenue pour tels systèmes par la méthode de Lyapunov. Plus précisément, la solution de la dynamique lente du système complet est approchée par la solution du système réduit lorsque le paramètre de la perturbation est suffisamment petit. Deuxièmement, le théorème de Tikhonov est établi pour les systèmes hyperboliques linéaires singulièrement perturbés de lois d'équilibre où les vitesses de transport et les termes sources sont à la fois dépendant du paramètre de la perturbation ainsi que les conditions aux bords. Sous des hypothèses sur la continuité de ces termes et sous la condition de la stabilité, l'estimation de l'erreur entre la dynamique lente du système complet et le système réduit est obtenue en fonction de l'ordre du paramètre de la perturbation. Troisièmement, nous considérons des systèmes EDO-EDP couplés singulièrement perturbés. La stabilité des deux sous-systèmes implique la stabilité du système complet où le paramètre de la perturbation est introduit dans la dynamique de l'EDP. D'autre part, cela n'est pas valable pour le système où le paramètre de la perturbation est présent dans l'EDO. Le théorème Tikhonov pour ces systèmes EDO-EDP couplés est prouvé par la technique de Lyapunov. Enfin, la synthèse de la commande aux bords est abordée en exploitant la méthode des perturbations singulières. Le système réduit converge en temps fini. La synthèse du contrôle aux bords est mise en œuvre pour deux applications différentes afin d'illustrer les résultats principaux de ce travail. / Systems modeled by partial differential equations (PDEs) with infinite dimensional dynamics are relevant for a wide range of physical networks. The control and stability analysis of such systems become a challenge area. Singularly perturbed systems, containing multiple time scales, often occur naturally in physical systems due to the presence of small parasitic parameters, typically small time constants, masses, inductances, moments of inertia. Singular perturbation was introduced in control engineering in late $1960$s, its assimilation in control theory has rapidly developed and has become a tool for analysis and design of control systems. Singular perturbation is a way of neglecting the fast transition and considering them in a separate fast time scale. The present thesis is concerned with a class of linear hyperbolic systems with multiple time scales modeled by a small perturbation parameter. Firstly we study a class of singularly perturbed linear hyperbolic systems of conservation laws. Since the system contains two time scales, by setting the perturbation parameter to zero, the two subsystems, namely the reduced subsystem and the boundary-layer subsystem, are formally computed. The stability of the full system implies the stability of both subsystems. However a counterexample is used to illustrate that the stability of the two subsystems is not enough to guarantee the full system's stability. This shows a major difference with what is well known for linear finite dimensional systems. Moreover, under certain conditions, the Tikhonov approximation for such system is achieved by Lyapunov method. Precisely, the solution of the slow dynamics of the full system is approximated by the solution of the reduced subsystem for sufficiently small perturbation parameter. Secondly the Tikhonov theorem is established for singularly perturbed linear hyperbolic systems of balance laws where the transport velocities and source terms are both dependent on the perturbation parameter as well as the boundary conditions. Under the assumptions on the continuity for such terms and under the stability condition, the estimate of the error between the slow dynamics of the full system and the reduced subsystem is the order of the perturbation parameter. Thirdly, we consider singularly perturbed coupled ordinary differential equation ODE-PDE systems. The stability of both subsystems implies that of the full system where the perturbation parameter is introduced into the dynamics of the PDE system. On the other hand, this is not true for system where the perturbation parameter is presented to the ODE. The Tikhonov theorem for such coupled ODE-PDE systems is proved by Lyapunov technique. Finally, the boundary control synthesis is achieved based on singular perturbation method. The reduced subsystem is convergent in finite time. Boundary control design to different applications are used to illustrate the main results of this work.

Submitted by Renata Lopes (renatasil82@gmail.com) on 2016-02-22T15:53:26Z No. of bitstreams: 1 williammassayukisakaguchiyamashita.pdf: 9277047 bytes, checksum: 3a0de46103c4b3001459c13047dfdb1a (MD5) / Approved for entry into archive by Adriana Oliveira (adriana.oliveira@ufjf.edu.br) on 2016-02-26T14:09:40Z (GMT) No. of bitstreams: 1 williammassayukisakaguchiyamashita.pdf: 9277047 bytes, checksum: 3a0de46103c4b3001459c13047dfdb1a (MD5) / Made available in DSpace on 2016-02-26T14:09:40Z (GMT). No. of bitstreams: 1 williammassayukisakaguchiyamashita.pdf: 9277047 bytes, checksum: 3a0de46103c4b3001459c13047dfdb1a (MD5) Previous issue date: 2014-03-25 / CAPES - Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / FAPEMIG - Fundação de Amparo à Pesquisa do Estado de Minas Gerais / Nas últimas décadas, a incidência global da dengue tem crescido dramaticamente favorecida pelo aumento da mobilidade humana e da urbanização. O estudo da população do mosquito é de grande importância para a saúde pública em países como o Brasil, onde as condições climáticas e ambientais são favoráveis para a propagação desta doença. Este trabalho baseia-se no estudo de modelos matemáticos que tratam do ciclo de vida do mosquito da dengue usando equações diferencias parciais. Nós investigamos a existência de solução na forma de onda viajante para ambos os modelos. Nós usamos um método semi-analítico combinando técnicas de Sistemas Dinâmicos (como a seção de Poincaré e análise local com base no Teorema de Hartman-Grobman) e integração numérica usando Matlab. / In recent decades the global incidence of dengue has grown dramatically by increased human mobility and urbanization. The study of the mosquito population is of great importance for public health in countries like Brazil, where climatic and environmental conditions are favorable for the propagation of this disease. This work is based on the study of mathematical models dealing with the life cycle of the dengue mosquito using partial differential equations. We investigate the existence of a solution in the form of travelling wave for both models. We use a semi-analytical method combining dynamical systems techniques (e.g. Poincaré section and local analysis based on Hartman-Grobman theorem) and numerical integration using Matlab.