To see the other types of publications on this topic, follow the link: Diffusion partial differential equations.
Dissertations / Theses on the topic 'Diffusion partial differential equations'
Create a spot-on reference in APA, MLA, Chicago, Harvard, and other styles
Consult the top 50 dissertations / theses for your research on the topic 'Diffusion partial differential equations.'
Next to every source in the list of references, there is an 'Add to bibliography' button. Press on it, and we will generate automatically the bibliographic reference to the chosen work in the citation style you need: APA, MLA, Harvard, Chicago, Vancouver, etc.
You can also download the full text of the academic publication as pdf and read online its abstract whenever available in the metadata.
Browse dissertations / theses on a wide variety of disciplines and organise your bibliography correctly.
1
Trojan, Alice von. "Finite difference methods for advection and diffusion." Title page, abstract and contents only, 2001. http://web4.library.adelaide.edu.au/theses/09PH/09phv948.pdf.
Full textAbstract:
Includes bibliographical references (leaves 158-163). Concerns the development of high-order finite-difference methods on a uniform rectangular grid for advection and diffuse problems with smooth variable coefficients. This technique has been successfully applied to variable-coefficient advection and diffusion problems. Demonstrates that the new schemes may readily be incorporated into multi-dimensional problems by using locally one-dimensional techniques, or that they may be used in process splitting algorithms to solve complicatef time-dependent partial differential equations.
2
Chow, Tanya L. M., of Western Sydney Macarthur University, and Faculty of Business and Technology. "Systems of partial differential equations and group methods." THESIS_FBT_XXX_Chow_T.xml, 1996. http://handle.uws.edu.au:8081/1959.7/43.
Full textAbstract:
This thesis is concerned with the derivation of similarity solutions for one-dimensional coupled systems of reaction - diffusion equations, a semi-linear system and a one-dimensional tripled system. The first area of research in this thesis involves a coupled system of diffusion equations for the existence of two distinct families of diffusion paths. Constructing one-parameter transformation groups preserving the invariance of this system of equations enables similarity solutions for this coupled system to be derived via the classical and non-classical procedures. This system of equation is the uncoupled in the hope of recovering further similarity solutions for the system. Once again, one-parameter groups leaving the uncoupled system invariant are obtained, enabling similarity solutions for the system to be elicited. A one-dimensional pattern formation in a model of burning forms the next component of this thesis. The primary focus of this area is the determination of similarity solutions for this reaction - diffusion system by means of one-parameter transformation group methods. Consequently, similarity solutions which are a generalisation of the solutions of the one-dimensional steady equations derived by Forbes are deduced. Attention in this thesis is then directed toward a semi-linear coupled system representing a predator - prey relationship. Two approaches to solving this system are made using the classical procedure, leading to one-parameter transformation groups which are instrumental in elicting the general similarity solution for this system. A triple system of equations representing a one-dimensional case of diffusion in the presence of three diffusion paths constitutes the next theme of this thesis. In association with the classical and non-classical procedures, the derivation of one-parameter transformation groups leaving this system invariant enables similarity solutions for this system to be deduced. The final strand of this thesis involves a one- dimensional case of the general linear system of coupled diffusion equations with cross-effects for which one-parameter transformation group methods are once more employed. The one-parameter groups constructed for this system prove instrumental in enabling the attainment of similarity solutions for this system to be accomplishedFaculty of Business and Technology
3
Garvie, Marcus Roland. "Analysis of a reaction-diffusion system of λ-w type." Thesis, Durham University, 2003. http://etheses.dur.ac.uk/4105/.
Full textAbstract:
The author studies two coupled reaction-diffusion equations of 'λ-w' type, on an open, bounded, convex domain Ω C R(^d) (d ≤ 3), with a boundary of class C², and homogeneous Neumann boundary conditions. The equations are close to a supercritical Hopf bifurcation in the reaction kinetics, and are model equations for oscillatory reaction-diffusion equations. Global existence, uniqueness and continuous dependence on initial data of strong and weak solutions are proved using the classical Faedo-Galerkin method of Lions and compactness arguments. The work provides a complete case study for the application of this method to systems of nonlinear reaction-diffusion equations. The author also undertook the numerical analysis of the reaction-diffusion system. Results are presented for a fully-practical piecewise linear finite element method by mimicking results in the continuous case. Semi-discrete and fully-discrete error estimates are proved after establishing a priori bounds for various norms of the approximate solutions. Finally, the theoretical results are illustrated and verified via the numerical simulation of periodic plane waves in one space dimension, and preliminary results representing target patterns and spiral solutions presented in two space dimensions.
4
Mbroh, Nana Adjoah. "On the method of lines for singularly perturbed partial differential equations." University of the Western Cape, 2017. http://hdl.handle.net/11394/5679.
Full textAbstract:
Magister Scientiae - MScMany chemical and physical problems are mathematically described by partial differential equations (PDEs). These PDEs are often highly nonlinear and therefore have no closed form solutions. Thus, it is necessary to recourse to numerical approaches to determine suitable approximations to the solution of such equations. For solutions possessing sharp spatial transitions (such as boundary or interior layers), standard numerical methods have shown limitations as they fail to capture large gradients. The method of lines (MOL) is one of the numerical methods used to solve PDEs. It proceeds by the discretization of all but one dimension leading to systems of ordinary di erential equations. In the case of time-dependent PDEs, the MOL consists of discretizing the spatial derivatives only leaving the time variable continuous. The process results in a system to which a numerical method for initial value problems can be applied. In this project we consider various types of singularly perturbed time-dependent PDEs. For each type, using the MOL, the spatial dimensions will be discretized in many different ways following fitted numerical approaches. Each discretisation will be analysed for stability and convergence. Extensive experiments will be conducted to confirm the analyses.
5
Palitta, Davide. "Preconditioning strategies for the numerical solution of convection-diffusion partial differential equations." Master's thesis, Alma Mater Studiorum - Università di Bologna, 2014. http://amslaurea.unibo.it/7464/.
Full textAbstract:
Il trattamento numerico dell'equazione di convezione-diffusione con le relative condizioni al bordo, comporta la risoluzione di sistemi lineari algebrici di grandi dimensioni in cui la matrice dei coefficienti è non simmetrica. Risolutori iterativi basati sul sottospazio di Krylov sono ampiamente utilizzati per questi sistemi lineari la cui risoluzione risulta particolarmente impegnativa nel caso di convezione dominante. In questa tesi vengono analizzate alcune strategie di precondizionamento, atte ad accelerare la convergenza di questi metodi iterativi. Vengono confrontati sperimentalmente precondizionatori molto noti come ILU e iterazioni di tipo inner-outer flessibile.
Nel caso in cui i coefficienti del termine di convezione siano a variabili separabili, proponiamo una nuova strategia di precondizionamento basata sull'approssimazione, mediante equazione matriciale, dell'operatore differenziale di convezione-diffusione. L'azione di questo nuovo precondizionatore sfrutta in modo opportuno recenti risolutori efficienti per equazioni matriciali lineari. Vengono riportati numerosi esperimenti numerici per studiare la dipendenza della performance dei diversi risolutori dalla scelta del termine di convezione, e dai parametri di discretizzazione.
6
Mateos, González Álvaro. "Asymptotic Analysis of Partial Differential Equations Arising in Biological Processes of Anomalous Diffusion." Thesis, Lyon, 2017. http://www.theses.fr/2017LYSEN069/document.
Full textAbstract:
Cette thèse est consacrée à l'analyse asymptotique d'équations aux dérivées partielles issues de modèles de déplacement sous-diffusif en biologie cellulaire. Notre motivation biologique est fondée sur les nombreuses observation récentes de protéinescytoplasmiques dont le déplacement aléatoire dévié de la diffusion Fickienne normale. Dans la première partie, nous étudions la décroissance auto-similaire de la solution d'une équation de renouvellement à queue lourde vers un état stationnaire. Les idéesmises en jeu sont inspirées de méthodes d'entropie relative. Nos principaux apports sont la preuve d'un taux de décroissance en norme L1 vers la loi de l'arc-sinus et l'introduction d'une fonction pivot spécifique dans une méthode d'entropie relative.La seconde partie porte sur la limite hyperbolique d'une équation de renouvellement structurée en âge et à sauts en espace. Nous y prouvons un résultat de « stabilité » : les solutions des problèmes rééchelonnés à ε > 0 convergent lorsque ε --> 0 vers la solution de viscosité de l'équation de Hamilton-Jacobi limite des problèmes à ε > 0. Les outilsmis en jeu proviennent de la théorie des équations de Hamilton-Jacobi.Ce travail présente trois idées intéressantes. La première est celle de prouver le résultat de convergence sur la condition de bord du problème plutôt que d'utiliser des fonctions test perturbées. La deuxième consiste en l'introduction de termes correcteurslogarithmiques en temps dans des estimations a priori ne découlant pas directementdu principe du maximum. Cela est dû à la non-existence d'un équilibre du problèmehomogène en espace. La troisième est une estimation précise de la décroissance de l'influence de la condition initiale sur le terme de renouvellement. Elle correspond à une estimation fine d'une version non-locale de la dérivée temporelle de la solution. Au cours de cette thèse, des simulations numériques de type Monte Carlo, schémas volumes finis, Lax-Friedrichs et Weighted Essentially Non Oscillating ont été réaliséesThis thesis is devoted to the asymptotic analysis of partial differential equations modelling subdiffusive random motion in cell biology. The biological motivation for this work is the numerous recent observations of cytoplasmic proteins whose random motion deviates from normal Fickian diffusion. In the first part, we study the self-similar decay towards a steady state of the solution of a heavy-tailed renewal equation. The ideas therein are inspired from relative entropy methods. Our main contributions are the proof of an L1 decay rate towards the arc-sine distribution and the introduction of a specific pivot function in a relative entropy method.The second part treats the hyperbolic limit of an age-structured space-jump renewal equation. We prove a "stability" result: the solutions of the rescaled problems at ε > 0 converge as ε --> 0 towards the viscosity solution of the limiting Hamilton-Jacobi equation of the ε > 0 problems. The main mathematical tools used come from the theory of Hamilton-Jacobi equations. This work presents three interesting ideas. The first is that of proving the convergence result on the boundary condition of the studied problem rather than using perturbed test functions. The second consists in the introduction of time-logarithmic correction termsin a priori estimates that do not follow directly from the maximum principle. That is due to the non-existence of a suitable equilibrium for the space-homogenous problem. The third is a precise estimate of the decay of the inuence of the initial condition on the renewal term. This is tantamount to a refined estimate of a non-local version of the time derivative of the solution. Throughout this thesis, we have performed numerical simulations of different types: Monte Carlo, finite volume schemes, Lax-Friedrichs schemes and Weighted Essentially Non Oscillating schemes
7
Manay, Siddharth. "Applications of anti-geometric diffusion of computer vision : thresholding, segmentation, and distance functions." Diss., Georgia Institute of Technology, 2003. http://hdl.handle.net/1853/33626.
Full text
8
Perella, Andrew James. "A class of Petrov-Galerkin finite element methods for the numerical solution of the stationary convection-diffusion equation." Thesis, Durham University, 1996. http://etheses.dur.ac.uk/5381/.
Full textAbstract:
A class of Petrov-Galerkin finite element methods is proposed for the numerical solution of the n dimensional stationary convection-diffusion equation. After an initial review of the literature we describe this class of methods and present both asymptotic and nonasymptotic error analyses. Links are made with the classical Galerkin finite element method and the cell vertex finite volume method. We then present numerical results obtained for a selection of these methods applied to some standard test problems. We also describe extensions of these methods which enable us to solve accurately for derivative values of the solution.
9
Lin, Xuelei. "Preconditioning techniques for all-at-once linear systems arising from advection diffusion equations." HKBU Institutional Repository, 2020. https://repository.hkbu.edu.hk/etd_oa/803.
Full textAbstract:
In this thesis, we mainly study preconditioning techniques for all-at-once linear systems arising from discretization of three types of time-dependent advection-diffusion equation: linear diffusion equation, constant-coefficients advection-diffusion equation, time-fractional sub-diffusion equation. The proposed preconditioners are used with Krylov subspace solvers. The preconditioner developed for linear diffusion equation is based on -circulant ap- proximation of temporal discretization. Diagonalizability, clustering of spectrum and identity-plus-low-rank decomposition are derived for the preconditioned matrix. We also show that generalized minimal residual (GMRES) solver for the preconditioned system has a linear convergence rate independent of matrix-size. The preconditioner for constant-coefficients advection-diffusion equation is based on approximating the discretization of advection term with a matrix diagonalizable by sine transform. Eigenvalues of the preconditioned matrix are proven to be lower and upper bounded by positive constants independent of discretization parameters. Moreover, as the preconditioner is based on spatial approximation, it is also applicable to steady-state problem. We show that GMRES for the preconditioned steady-state problem has a linear convergence rate independent of matrix size. The preconditioner for time-fractional sub-diffusion equation is based on approximat- ing the discretization of diffusion term with a matrix diagonalizable by sine transform. We show that the condition number of the preconditioned matrix is bounded by a constant independent of discretization parameters so that the normalized conjugate gradient (NCG) solver for the preconditioned system has a linear convergence rate independent of discretization parameters and matrix size. Fast implementations based on fast Fourier transform (FFT), fast sine transform (FST) or multigrid approximation are proposed for the developed preconditioners. Numerical results are reported to show the performance of the developed preconditioners
10
Shen, Wensheng. "Computer Simulation and Modeling of Physical and Biological Processes using Partial Differential Equations." UKnowledge, 2007. http://uknowledge.uky.edu/gradschool_diss/501.
Full textAbstract:
Scientific research in areas of physics, chemistry, and biology traditionally depends purely on experimental and theoretical methods. Recently numerical simulation is emerging as the third way of science discovery beyond the experimental and theoretical approaches. This work describes some general procedures in numerical computation, and presents several applications of numerical modeling in bioheat transfer and biomechanics, jet diffusion flame, and bio-molecular interactions of proteins in blood circulation.
A three-dimensional (3D) multilayer model based on the skin physical structure is developed to investigate the transient thermal response of human skin subject to external heating. The temperature distribution of the skin is modeled by a bioheat transfer equation. Different from existing models, the current model includes water evaporation and diffusion, where the rate of water evaporation is determined based on the theory of laminar boundary layer. The time-dependent equation is discretized using the Crank-Nicolson scheme. The large sparse linear system resulted from discretizing the governing partial differential equation is solved by GMRES solver.
The jet diffusion flame is simulated by fluid flow and chemical reaction. The second-order backward Euler scheme is applied for the time dependent Navier-Stokes equation. Central difference is used for diffusion terms to achieve better accuracy, and a monotonicity-preserving upwind difference is used for convective ones. The coupled nonlinear system is solved via the damped Newton's method. The Newton Jacobian matrix is formed numerically, and resulting linear system is ill-conditioned and is solved by Bi-CGSTAB with the Gauss-Seidel preconditioner.
A novel convection-diffusion-reaction model is introduced to simulate fibroblast growth factor (FGF-2) binding to cell surface molecules of receptor and heparan sulfate proteoglycan and MAP kinase signaling under flow condition. The model includes three parts: the flow of media using compressible Navier-Stokes equation, the transport of FGF-2 using convection-diffusion transport equation, and the local binding and signaling by chemical kinetics. The whole model consists of a set of coupled nonlinear partial differential equations (PDEs) and a set of coupled nonlinear ordinary differential equations (ODEs). To solve the time-dependent PDE system we use second order implicit Euler method by finite volume discretization. The ODE system is stiff and is solved by an ODE solver VODE using backward differencing formulation (BDF). Findings from this study have implications with regard to regulation of heparin-binding growth factors in circulation.
11
Marahrens, Daniel. "On some nonlinear partial differential equations for classical and quantum many body systems." Thesis, University of Cambridge, 2012. https://www.repository.cam.ac.uk/handle/1810/244203.
Full textAbstract:
This thesis deals with problems arising in the study of nonlinear partial differential equations arising from many-body problems. It is divided into two parts: The first part concerns the derivation of a nonlinear diffusion equation from a microscopic stochastic process. We give a new method to show that in the hydrodynamic limit, the particle densities of a one-dimensional zero range process on a periodic lattice converge to the solution of a nonlinear diffusion equation. This method allows for the first time an explicit uniform-in-time bound on the rate of convergence in the hydrodynamic limit. We also discuss how to extend this method to the multi-dimensional case. Furthermore we present an argument, which seems to be new in the context of hydrodynamic limits, how to deduce the convergence of the microscopic entropy and Fisher information towards the corresponding macroscopic quantities from the validity of the hydrodynamic limit and the initial convergence of the entropy. The second part deals with problems arising in the analysis of nonlinear Schrödinger equations of Gross-Pitaevskii type. First, we consider the Cauchy problem for (energy-subcritical) nonlinear Schrödinger equations with sub-quadratic external potentials and an additional angular momentum rotation term. This equation is a well-known model for superfluid quantum gases in rotating traps. We prove global existence (in the energy space) for defocusing nonlinearities without any restriction on the rotation frequency, generalizing earlier results given in the literature. Moreover, we find that the rotation term has a considerable influence in proving finite time blow-up in the focusing case. Finally, a mathematical framework for optimal bilinear control of nonlinear Schrödinger equations arising in the description of Bose-Einstein condensates is presented. The obtained results generalize earlier efforts found in the literature in several aspects. In particular, the cost induced by the physical work load over the control process is taken into account rather then often used L^2- or H^1-norms for the cost of the control action. We prove well-posedness of the problem and existence of an optimal control. In addition, the first order optimality system is rigorously derived. Also a numerical solution method is proposed, which is based on a Newton type iteration, and used to solve several coherent quantum control problems.
12
Qiao, Zhonghua. "Numerical solution for nonlinear Poisson-Boltzmann equations and numerical simulations for spike dynamics." HKBU Institutional Repository, 2006. http://repository.hkbu.edu.hk/etd_ra/727.
Full text
13
Lao, Kun Leng. "Multigrid algorithm based on cyclic reduction for convection diffusion equations." Thesis, University of Macau, 2010. http://umaclib3.umac.mo/record=b2148274.
Full text
14
Rolland, Guillaume. "Global existence and fast-reaction limit in reaction-diffusion systems with cross effects." Phd thesis, École normale supérieure de Cachan - ENS Cachan, 2012. http://tel.archives-ouvertes.fr/tel-00785757.
Full textAbstract:
This thesis is devoted to the study of parabolic systems of partial differential equations arising in mass action kinetics chemistry, population dynamics and electromigration theory. We are interested in the existence of global solutions, uniqueness of weak solutions, and in the fast-reaction limit in a reaction-diffusion system. In the first chapter, we study two cross-diffusion systems. We are first interested in a population dynamics model, where cross effects in the interactions between the different species are modeled by non-local operators. We prove the well-posedness of the corresponding system for any space dimension. We are then interested in a cross-diffusion system which arises as the fast-reaction limit system in a classical system for the chemical reaction C1+C2=C3. We prove the convergence when k goes to infinity of the solution of the system with finite reaction speed k to a global solution of the limit system. The second chapter contains new global existence results for some reaction-diffusion systems. For networks of elementary chemical reactions of the type Ci+Cj=Ck and under Mass Action Kinetics assumption, we prove the existence and uniqueness of global strong solutions, for space dimensions N<6 in the semi-linear case, and N<4 in the quasi-linear case. We also prove the existence of global weak solutions for a class of parabolic quasi-linear systems with at most quadratic non-linearities and with initial data that are only assumed to be nonnegative and integrable. In the last chapter, we generalize a global well-posedness result for reaction-diffusion systems whose nonlinearities have a "triangular" structure, for which we now take into account advection terms and time and space dependent diffusion coefficients. The latter result is then used in a Leray-Schauder fixed point argument to prove the existence of global solutions in a diffusion-electromigration system.
15
Yue, Wen. "Absolute continuity of the laws, existence and uniqueness of solutions of some SDEs and SPDEs." Thesis, University of Manchester, 2014. https://www.research.manchester.ac.uk/portal/en/theses/absolute-continuity-of-the-laws-existence-and-uniqueness-of-solutions-of-some-sdes-and-spdes(2bc80de8-7c36-453f-a7c2-69fa4ee0e705).html.
Full textAbstract:
This thesis consists of four parts. In the first part we recall some background theory that will be used throughout the thesis. In the second part, we studied the absolute continuity of the laws of the solutions of some perturbed stochastic differential equaitons(SDEs) and perturbed reflected SDEs using Malliavin calculus. Because the extra terms in the perturbed SDEs involve the maximum of the solution itself, the Malliavin differentiability of the solutions becomes very delicate. In the third part, we studied the absolute continuity of the laws of the solutions of the parabolic stochastic partial differential equations(SPDEs) with two reflecting walls using Malliavin calculus. Our study is based on Yang and Zhang \cite{YZ1}, in which the existence and uniqueness of the solutions of such SPDEs was established. In the fourth part, we gave the existence and uniqueness of the solutions of the elliptic SPDEs with two reflecting walls and general diffusion coefficients.
16
Ducasse, Romain. "Équations et systèmes de réaction-diffusion en milieux hétérogènes et applications." Thesis, Paris Sciences et Lettres (ComUE), 2018. http://www.theses.fr/2018PSLEE054/document.
Full textAbstract:
Cette thèse est consacrée à l'étude des équations et systèmes de réaction-diffusion dans des milieux hétérogènes. Elle est divisée en deux parties. La première est dédiée à l'étude des équations de réaction-diffusion dans des milieux périodiques. Nous nous intéressons en particulier aux équations posées dans des domaines qui ne sont pas l'espace entier $\mathbb{R}^{N}$, mais des domaines périodiques, avec des "obstacles". Dans un premier chapitre, nous étudions l'effet de la géométrie du domaine sur la vitesse d'invasion des solutions. Après avoir dérivé une formule de type Freidlin-Gartner, nous construisons des domaines où la vitesse d'invasion est strictement inférieure à la vitesse critique des fronts. Nous donnons également des critères géométriques qui garantissent l'existence de directions où l'invasion se produit à la vitesse critique. Dans le chapitre suivant, nous donnons des conditions nécessaires et suffisantes pour garantir que l'invasion ait lieu, après quoi nous construisons des domaines où des phénomènes intermédiaires (blocage, invasion orientée) se produisent. La deuxième partie de cette thèse est consacrée à l'étude de modèles décrivant l'influence de lignes à diffusion rapide (une route, par exemple) sur la propagation d'espèces invasives. Il a en effet été observé que certaines espèces, dont le moustique-tigre, envahissent plus rapidement que prévu certaines zones proches du réseau routier. Nous étudions deux modèles : le premier décrit l'influence d'une route courbe sur la propagation. Nous nous intéressons en particulier au cas de deux routes non-parallèles. Le second modèle décrit l'influence d'une route sur une niche écologique, en présence d'un changement climatique. Le résultat principal est que l'effet de la route est ambivalent : si la niche est stationnaire, alors l'effet de la route est délétère. Cependant, si la niche se déplace, suite à un changement climatique, nous montrons que la route peut permettre à une population de survivre. Pour étudier ce second modèle, nous développons une notion de valeur propre principale généralisée pour des systèmes de type KPP, et nous dérivons une inégalité de Harnack, qui est nouvelle pour ce type de systèmesThis thesis is dedicated to the study of reaction-diffusion equations and systems in heterogeneous media. It is divided into two parts. The first one is devoted to the study of reaction-diffusion equations in periodic media. We pay a particular attention to equations set on domains that are not the whole space $\mathbb{R}^{N}$, but periodic domains, with "obstacles". In a first chapter, we study how the geometry of the domain can influence the speed of invasion of solutions. After establishing a Freidlin-Gartner type formula, we construct domains where the speed of invasion is strictly less than the critical speed of fronts. We also give geometric criteria to ensure the existence of directions where the invasion occurs with the critical speed. In the second chapter, we give necessary and sufficient conditions to ensure that invasion occurs, and we construct domains where intermediate phenomena (blocking, oriented invasion) occur. The second part of this thesis is dedicated to the study of models describing the influence of lines with fast diffusion (a road, for instance) on the propagation of invasive species. Indeed, it was observed that some species, such as the tiger mosquito, invade faster than expected some areas along the road-network. We study two models : the first one describes the influence of a curved road on the propagation. We study in particular the case of two non-parallel roads. The second model describes the influence of a road on an ecological niche, in presence of climate change. The main result is that the effect of the road is ambivalent: if the niche is stationary, then effect of the road is deleterious. However, if the niche moves, because of a shifting climate, the road can actually help the population to persist. To study this model, we introduce a notion of generalized principal eigenvalue for KPP-type systems, and we derive a Harnack inequality, that is new for this type of systems
17
Le, Balc'h Kévin. "Contrôlabilité de systèmes de réaction-diffusion non linéaires." Thesis, Rennes, École normale supérieure, 2019. http://www.theses.fr/2019ENSR0016/document.
Full textAbstract:
Cette thèse est consacrée au contrôle de quelques équations aux dérivées partielles non linéaires. On s’intéresse notamment à des systèmes paraboliques de réaction-diffusion non linéaires issus de la cinétique chimique. L’objectif principal est de démontrer des résultats de contrôlabilité locale ou globale, en temps petit, ou en temps grand.Dans une première partie, on démontre un résultat de contrôlabilité locale à des états stationnaires positifs en temps petit, pour un système de réaction-diffusion non linéaire.Dans une deuxième partie, on résout une question de contrôlabilité globale à zéro en temps petit pour un système 2 × 2 de réaction-diffusion non linéaire avec un couplage impair.La troisième partie est consacrée au célèbre problème ouvert d’Enrique Fernández-Cara et d’Enrique Zuazua des années 2000 concernant la contrôlabilité globale à zéro de l’équation de la chaleur faiblement non linéaire. On démontre un résultat de contrôlabilité globale à états positifs en temps petit et un résultat de contrôlabilité globale à zéro en temps long.La dernière partie, rédigée en collaboration avec Karine Beauchard et Armand Koenig, est une incursion vers l’hyperbolique. On étudie des systèmes linéaires à coefficients constants, couplant une dynamique transport avec une dynamique parabolique. On identifie leur temps minimal de contrôle et l’influence de leur structure algébrique sur leurs propriétés de contrôleThis thesis is devoted to the control of nonlinear partial differential equations. We are mostly interested in nonlinear parabolic reaction-diffusion systems in reaction kinetics. Our main goal is to prove local or global controllability results in small time or in large time.In a first part, we prove a local controllability result to nonnegative stationary states in small time, for a nonlinear reaction-diffusion system.In a second part, we solve a question concerning the global null-controllability in small time for a 2 × 2 nonlinear reaction-diffusion system with an odd coupling term.The third part focuses on the famous open problem due to Enrique Fernndez-Cara and Enrique Zuazua in 2000, concerning the global null-controllability of the weak semi-linear heat equation. We show that the equation is globally nonnegative controllable in small time and globally null-controllable in large time.The last part, which is a joint work with Karine Beauchard and Armand Koenig, enters the hyperbolic world. We study linear parabolic-transport systems with constant coeffcients. We identify their minimal time of control and the influence of their algebraic structure on the controllability properties
18
Phan, Van Long Em. "Analyse asymptotique de réseaux complexes de systèmes de réaction-diffusion." Thesis, Le Havre, 2015. http://www.theses.fr/2015LEHA0012/document.
Full textAbstract:
Le fonctionnement d'un neurone, unité fondamentale du système nerveux, intéresse de nombreuses disciplines scientifiques. Il existe ainsi des modèles mathématiques qui décrivent leur comportement par des systèmes d'EDO ou d'EDP. Plusieurs de ces modèles peuvent ensuite être couplés afin de pouvoir étudier le comportement de réseaux, systèmes complexes au sein desquels émergent des propriétés. Ce travail présente, dans un premier temps, les principaux mécanismes régissant ce fonctionnement pour en comprendre la modélisation. Plusieurs modèles sont alors présentés, jusqu'à celui de FitzHugh-Nagumo (FHN), qui présente une dynamique très intéressante.C'est sur l'étude théorique mais également numérique de la dynamique asymptotique et transitoire du modèle de FHN en EDO, que se concentre la seconde partie de cette thèse. A partir de cette étude, des réseaux d'interactions d'EDO sont construits en couplant les systèmes dynamiques précédemment étudiés. L'étude du phénomène de synchronisation identique au sein de ces réseaux montre l'existence de propriétés émergentes pouvant être caractérisées par exemple par des lois de puissance. Dans une troisième partie, on se concentre sur l'étude du système de FHN dans sa version EDP. Comme la partie précédente, des réseaux d'interactions d'EDP sont étudiés. On entreprend dans cette partie une étude théorique et numérique. Dans la partie théorique, on montre l'existence de l'attracteur global dans l'espace L2(Ω)nd et on donne des conditions suffisantes de synchronisation. Dans la partie numérique, on illustre le phénomène de synchronisation ainsi que l'émergence de lois générales telles que les lois puissances ou encore la formation de patterns, et on étudie l'effet de l'ajout de la dimension spatiale sur la synchronisationThe neuron, a fundamental unit in the nervous system, is a point of interest in many scientific disciplines. Thus, there are some mathematical models that describe their behavior by ODE or PDE systems. Many of these models can then be coupled in order to study the behavior of networks, complex systems in which the properties emerge. Firstly, this work presents the main mechanisms governing the neuron behaviour in order to understand the different models. Several models are then presented, including the FitzHugh-Nagumo one, which has a interesting dynamic. The theoretical and numerical study of the asymptotic and transitory dynamics of the aforementioned model is then proposed in the second part of this thesis. From this study, the interaction networks of ODE are built by coupling previously dynamic systems. The study of identical synchronization phenomenon in these networks shows the existence of emergent properties that can be characterized by power laws. In the third part, we focus on the study of the PDE system of FHN. As the previous part, the interaction networks of PDE are studied. We have in this section a theoretical and numerical study. In the theoretical part, we show the existence of the global attractor on the space L2(Ω)nd and give the sufficient conditions for identical synchronization. In the numerical part, we illustrate the synchronization phenomenon, also the general laws of emergence such as the power laws or the patterns formation. The diffusion effect on the synchronization is studied
19
Bruna, Maria. "Excluded-volume effects in stochastic models of diffusion." Thesis, University of Oxford, 2012. http://ora.ox.ac.uk/objects/uuid:020c2d3e-5fef-478c-9861-553cd310daf5.
Full textAbstract:
Stochastic models describing how interacting individuals give rise to collective behaviour have become a widely used tool across disciplines—ranging from biology to physics to social sciences. Continuum population-level models based on partial differential equations for the population density can be a very useful tool (when, for large systems, particle-based models become computationally intractable), but the challenge is to predict the correct macroscopic description of the key attributes at the particle level (such as interactions between individuals and evolution rules). In this thesis we consider the simple class of models consisting of diffusive particles with short-range interactions. It is relevant to many applications, such as colloidal systems and granular gases, and also for more complex systems such as diffusion through ion channels, biological cell populations and animal swarms. To derive the macroscopic model of such systems, previous studies have used ad hoc closure approximations, often generating errors. Instead, we provide a new systematic method based on matched asymptotic expansions to establish the link between the individual- and the population-level models. We begin by deriving the population-level model of a system of identical Brownian hard spheres. The result is a nonlinear diffusion equation for the one-particle density function with excluded-volume effects enhancing the overall collective diffusion rate. We then expand this core problem in several directions. First, for a system with two types of particles (two species) we obtain a nonlinear cross-diffusion model. This model captures both alternative notions of diffusion, the collective diffusion and the self-diffusion, and can be used to study diffusion through obstacles. Second, we study the diffusion of finite-size particles through confined domains such as a narrow channel or a Hele–Shaw cell. In this case the macroscopic model depends on a confinement parameter and interpolates between severe confinement (e.g., a single- file diffusion in the narrow channel case) and an unconfined situation. Finally, the analysis for diffusive soft spheres, particles with soft-core repulsive potentials, yields an interaction-dependent non-linear term in the diffusion equation.
20
Patout, Florian. "Analyse asymptotique d'équations intégro-différentielles : modèles d'évolution et de dynamique des populations." Thesis, Lyon, 2019. http://www.theses.fr/2019LYSEN044/document.
Full textAbstract:
Cette thèse est consacrée à l’étude de phénomènes de propagation et de concentration dans des modèles d’équations intégro-différentielles venant de la écologie. On étudie certaines équations de réaction-diffusion non locales apparaissant en dynamique de populations, ainsi que des modèles représentant l’évolution Darwinienne avec un mode de reproduction sexué.Dans une première partie, nous étudions la propagation spatiale pour une équation de réaction-diffusion ou la dispersion opère via un noyau de convolution à queue lourde. Nous mesurons de manière précise l’accélération du front de propagation de la solution. Nous proposons également une échelle adaptée pour mesurer les «petites» mutations. Dans les deux cas nous utilisons le formalisme des équations de Hamilton-Jacobi.Dans un second temps nous étudions un modèle de génétique quantitative, avec un mode de reproduction sexuée. Un petit paramètre mesure la déviation entre le trait des descendants est la moyenne des traits des parents. Dans le régime où ce paramètre est petit nous étudions l’existence de solutions stationnaires, puis le problème de Cauchy lié à ce modèle. Les solutions se concentrent autour des optima de sélection, sous la forme de perturbations de distributions Gaussiennes avec petite variance fixée par le paramètre. Notre analyse généralise le cas linéaire de la reproduction asexuée en utilisant des outils d’analyse perturbative. Enfin dans une dernière partie nous fournissons des simulations numériques et des méthodes mathématiques pour étudier la dynamique interne des équilibres dans le régime de petite variance, pour les deux modes de reproduction : asexué et sexuéThis manuscript tackles propagation and concentration phenomena in different integro-differential equations with a background in ecology. We study non local reaction-diffusion equations from population dynamics, and models for Darwinian evolution with a sexual or asexual mode of reproduction, with a preference for the former.In a first part, we study spatial propagation for a reaction diffusion equation where dispersion acts through a fat tailed kernel. We measure accurately the acceleration of the propagation front of the population. We propose as well a scaling well adapted to “small mutations” when we consider the model in the context of adaptative dynamics. This scaling is very natural following the previous spatial investigation. In both cases we look at the long time behavior and we use the Hamilton-Jacobi framework. Then we turn our attention towards a quantitative genetics model, with a sexual mode of reproduction, imposed by the “infinitesimal operator”. In this non-linear setting, a small parameter tunes the deviation between the phenotypic trait of the offspring and the mean of the traits of the parents. In the regime where this parameter is small, we prove existence of stationary solutions, and their local uniqueness. We also provide an example of non-uniqueness in the case where the selection function admits several extrema. We prove that the solution concentrates around the points of minimum of the selection function. The analysis is carried by the small perturbations of special profiles : Gaussian distributions with small variance fixed by the parameter.We then study the stability of the Cauchy problem associated to the previous model. This time we prove that at all times, for a well prepared initial data, the solutions is arbitrary close to a Gaussian distribution with small variance. The proof follows the framework of the previous : we use perturbative analysis tools, but this time an even more precise description of the correctors is needed and we linearize the equation to obtain it. In a final part we show numerical simulations and different mathematical approaches to study inside dynamics of phenotypic lineages in the regime of small variance, with a moving environement
21
Bhikkaji, Bharath. "Model Reduction and Parameter Estimation for Diffusion Systems." Doctoral thesis, Uppsala : Acta Universitatis Upsaliensis : Univ.-bibl. [distributör], 2004. http://urn.kb.se/resolve?urn=urn:nbn:se:uu:diva-4252.
Full text
22
Ghazaryan, Anna R. "Nonlinear convective instability of fronts a case study /." Connect to resource, 2005. http://rave.ohiolink.edu/etdc/view?acc%5Fnum=osu1117552079.
Full textAbstract:
Thesis (Ph.D.)--Ohio State University, 2005.Title from first page of PDF file. Document formatted into pages; contains ix, 176 p.; also includes graphics. Includes bibliographical references (p. 172-176). Available online via OhioLINK's ETD Center
23
Duan, Xianglong. "Optimal transport and diffusion of currents." Thesis, Université Paris-Saclay (ComUE), 2017. http://www.theses.fr/2017SACLX054/document.
Full textAbstract:
Les travaux portent sur l'étude d'équations aux dérivées partielles à la charnière de la physique de la mécanique des milieux continus et de la géométrie différentielle, le point de départ étant le modèle d'électromagnétisme non-linéaire introduit par Max Born et Leopold Infeld en 1934 comme substitut aux traditionnelles équations linéaires de Maxwell. Ces équations sont remarquables par leurs liens avec la géométrie différentielle (surfaces extrémales dans l'espace de Minkowski) et ont connu un regain d'intérêt dans les années 90 en physique des hautes énergies (cordes et D-branes).Le travail se décompose en quatre chapitres.La théorie des systèmes paraboliques dégénérés d'EDP non-linéaires est fort peu développée, faute de pouvoir appliquer les principes de comparaison habituels (principe du maximum), malgré leur omniprésence dans de nombreuses applications (physique, mécanique, imagerie numérique, géométrie...). Dans le premier chapitre, on montre comment de tels systèmes peuvent être parfois dérivés, asymptotiquement, à partir de systèmes non-dissipatifs (typiquement des systèmes hyperboliques non-linéaires), par simple changement de variable en temps non-linéaire dégénéré à l'origine (où sont fixées les données initiales). L'avantage de ce point de vue est de pouvoir transférer certaines techniques hyperboliques vers les équations paraboliques, ce qui semble à première vue surprenant, puisque les équations paraboliques ont la réputation d'être plus facile à traiter (ce qui n'est pas vrai, en réalité, dans le cas de systèmes dégénérés). Le chapitre traite, comme prototype, du curve-shortening flow", qui est le plus simple des mouvements par courbure moyenne en co-dimension supérieure à un. Il est montré comment ce modèle peut être dérivé de la théorie des surfaces de dimension deux d'aire extrémale dans l'espace de Minkowski (correspondant aux cordes relativistes classiques) qui peut se ramener à un système hyperbolique. On obtient, presque automatiquement, l'équivalent parabolique des principes d'entropie relative et d'unicité fort-faible qu'il est, en fait, bien plus simple d'établir et de comprendre dans le cadre hyperbolique.Dans le second chapitre, la même méthode s'applique au système de Born-Infeld proprement dit, ce qui permet d'obtenir, à la limite, un modèle (non répertorié à notre connaissance) de Magnétohydrodynamique (MHD), où on retrouve à la fois une diffusivité non-linéaire dans l'équation d'induction magnétique et une loi de Darcy pour le champ de vitesse. Il est remarquable qu'un système d'apparence aussi lointaine des principes de base de la physique puisse être si directement déduit d'un modèle de physique aussi fondamental et géométrique que celui de Born-Infeld.Dans le troisième chapitre, un lien est établi entre des systèmes paraboliques et le concept de flot gradient de formes différentielles pour des métriques de transport. Dans le cas des formes volumes, ce concept a eu un succès extraordinaire dans le cadre de la théorie du transport optimal, en particulier après le travail fondateur de Felix Otto et de ses collaborateurs. Ce concept n'en est vraiment qu'à ses débuts: dans ce chapitre, on étudie une variante du «curve-shortening flow» étudié dans le premier chapitre, qui présente l'avantage d'être intégrable (en un certain sens) et de conduire à des résultats plus précis.Enfin, dans le quatrième chapitre, on retourne au domaine des EDP hyperboliques en considérant, dans le cas particulier des graphes, les surfaces extrémales de l'espace de Minkowski, de dimension et co-dimension quelconques. On parvient à montrer que les équations peuvent se reformuler sous forme d'un système élargi symétrique du premier ordre (ce qui assure automatiquement le caractère bien posé des équations) d'une structure remarquablement simple (très similaire à l'équation de Burgers) avec non linéarités quadratiques, dont le calcul n'a rien d'évidentOur work concerns about the study of partial differential equations at the hinge of the continuum physics and differential geometry. The starting point is the model of non-linear electromagnetism introduced by Max Born and Leopold Infeld in 1934 as a substitute for the traditional linear Maxwell's equations. These equations are remarkable for their links with differential geometry (extremal surfaces in the Minkowski space) and have regained interest in the 90s in high-energy physics (strings and D-branches).The thesis is composed of four chapters.The theory of nonlinear degenerate parabolic systems of PDEs is not very developed because they can not apply the usual comparison principles (maximum principle), despite their omnipresence in many applications (physics, mechanics, digital imaging, geometry, etc.). In the first chapter, we show how such systems can sometimes be derived, asymptotically, from non-dissipative systems (typically non-linear hyperbolic systems), by simple non-linear change of the time variable degenerate at the origin (where the initial data are set). The advantage of this point of view is that it is possible to transfer some hyperbolic techniques to parabolic equations, which seems at first sight surprising, since parabolic equations have the reputation of being easier to treat (which is not true , in reality, in the case of degenerate systems). The chapter deals with the curve-shortening flow as a prototype, which is the simplest exemple of the mean curvature flows in co-dimension higher than 1. It is shown how this model can be derived from the two-dimensional extremal surface in the Minkowski space (corresponding to the classical relativistic strings), which can be reduced to a hyperbolic system. We obtain, almost automatically, the parabolic version of the relative entropy method and weak-strong uniqueness, which, in fact, is much simpler to establish and understand in the hyperbolic framework.In the second chapter, the same method applies to the Born-Infeld system itself, which makes it possible to obtain, in the limit, a model (not listed to our knowledge) of Magnetohydrodynamics (MHD) where we have non-linear diffusions in the magnetic induction equation and the Darcy's law for the velocity field. It is remarkable that a system of such distant appearance of the basic principles of physics can be so directly derived from a model of physics as fundamental and geometrical as that of Born-Infeld.In the third chapter, a link is established between the parabolic systems and the concept of gradient flow of differential forms with suitable transport metrics. In the case of volume forms, this concept has had an extraordinary success in the field of optimal transport theory, especially after the founding work of Felix Otto and his collaborators. This concept is really only on its beginnings: in this chapter, we study a variant of the curve-shortening flow studied in the first chapter, which has the advantage of being integrable (in a certain sense) and lead to more precise results.Finally, in the fourth chapter, we return to the domain of hyperbolic EDPs considering, in the particular case of graphs, the extremal surfaces of the Minkowski space of any dimension and co-dimension. We can show that the equations can be reformulated in the form of a symmetric first-order enlarged system (which automatically ensures the well-posedness of the equations) of a remarkably simple structure (very similar to the Burgers equation) with quadratic nonlinearities, whose calculation is not obvious
24
Oliveira, Andrea Genovese de 1986. "Um sistema de equações parabólicas de reação-difusão modelando quimiotaxia." [s.n.], 2012. http://repositorio.unicamp.br/jspui/handle/REPOSIP/307414.
Full textAbstract:
Orientador: José Luiz BoldriniDissertação (mestrado) - Universidade Estadual de Campinas, Instituto de Matemática, Estatística e Computação Científica
Made available in DSpace on 2018-08-19T18:40:32Z (GMT). No. of bitstreams: 1 Oliveira_AndreaGenovesede_M.pdf: 1278255 bytes, checksum: f16ace92e18ff9cf5e8a4f8a66829f47 (MD5) Previous issue date: 2012
Resumo: Analisamos um sistema não linear parabólico de reação-difusão com duas equações definidas em ]0,T[x'ômega', (0 < T < 'infinito' e Q 'pertence' R³ limitado) e condições de fronteira do tipo Neumann. Tal sistema foi proposto para modelar o movimento de uma população de amebas unicelulares e tem como base o processo de locomoção chamado quimiotaxia positiva, na qual as amebas se movimentam em direção à região de alta concentração de uma certa substância química, que, neste caso, é produzida pelas próprias amebas. Embora adicionando os detalhes técnicos, este trabalho seguiu livremente o método de resolução proposto no artigo de A. Boy, Analysis for a System of Coupled Reaction-Diffusion Parabolic Equations Arising in Biology, Computers Math. Applic. Vol. 32, No. 4, páginas 15-21, 1996
Abstract: We will be analyzing a nonlinear parabolic reaction diffusion system with two equations, defined in ]0,T[x'omega', (0 < T < 'infinite' and Q 'belongs' R³) with Neumann boundary conditions. This system was proposed in order to model the movement of a population of single-cell amoebae and is based on the process of movement called chemotaxis, in which the amoebae move in the direction of the region of high concentration of a certain chemical substance, which, in this case, is produced by the amoebae themselves.While adding the technical details, this dissertation followed freely the solution method proposed in the paper: A. Boy, Analysis for a System of Coupled Reaction-Diffusion Parabolic Equations Arising in Biology, Computers Math. Applic. Vol. 32, No. 4, pages 15-21, 1996
Mestrado
Matematica
Mestre em Matemática
25
Houllier, Trescases Ariane. "Modélisation et Analyse Mathématique d'Equations aux Dérivées Partielles Issues de la Physique et de la Biologie." Thesis, Cachan, Ecole normale supérieure, 2015. http://www.theses.fr/2015DENS0037/document.
Full textAbstract:
Ce manuscrit présente des résultats d’analyse mathématique autour de deux exemples de problèmes singuliers d’équations aux dérivées partielles issus de la modélisation. I. Diffusion croisée en dynamique des populations. En dynamique des populations, les systèmes de réaction –diffusion croisée modélisent l’évolution de populations d’espèces en compétition avec un effet répulsif entre les individus. Pour ces systèmes fortement couplés, souvent non linéaires, une question aussi fondamentale que l’existence de solutions se révèle extrêmement complexe. Dans ce manuscrit, on introduit une approche basée sur des extensions récentes de lemmes de dualité et sur des méthodes d’entropie. On démontre l’existence de solutions faibles dans un cadre général de systèmes de réaction-diffusion croisée, ainsi que certaines propriétés qualitatives des solutions. II. Équation de Boltzmann en domaine borné. L’équation de Boltzmann, introduite en 1872, modélise la dynamique des gaz raréfiés hors équilibre. Malgré les nombreux résultats autour de la question de l’existence de solutions fortes proches de l’équilibre, très peu concernent le cas d’un domaine borné, situation pourtant fréquente dans les applications. Une raison de la difficulté du problème est l’irruption des singularités le long des trajectoires rasant le bord du domaine. Dans ce manuscrit, on présente une théorie de la régulation de l’équation de Boltzmann en domaine borné. Grâce à l’introduction d’une distance cinétique qui compense les singularités au bord, on montre des résultats de propagation de normes de Sobolev et de propagation C^1 en domaine convexe. En domaine non convexe, on montre un résultat de propagation de régularité BVThis manuscript presents results of mathematical analysis concerning two singular problems of partial differential equations coming from the modeling. I. Cross-diffusion in Population dynamics. In Population dynamics, reaction-cross diffusion systems model the evolution of the populations of competing species with a repulsive effect between individuals. For these strongly coupled, often non linear systems, a question as basic as the existence of solutions appears to be extremely complex. In this manuscript, we introduce an approach based on the most recent extensions of duality lemmas and on entropy methods. We prove the existence of weak solutions in a general setting of reaction-cross diffusion systems, as well as some qualitative properties of the solutions. II. Boltzmann equation in bounded domains The Boltzmann equation, introduced in 1872, model the evolution of a rarefied gas out of equilibrium. Despite the numerous results concerning the existence of strong solutions close to equilibrium, very few concern the case of bounded domain, though this situation is very useful in applications. A crucial reason of the difficulty of this problem is the formation of a singularity on the trajectories grazing the boundary. In this manuscript, we present a theory of the regularity of the Boltzmann equation in bounded domains. Thanks to the introduction of a kinetic distance which balances the singularity, we obtain results of propagation of Sobolev norms and C^1 propagation in convex domains. In non convex domains, we obtain the propagation of BV regularity
26
Leite, Jefferson Cruz dos Santos 1981. "Sistemas dinâmicos fuzzy aplicados a processos difusivos." [s.n.], 2011. http://repositorio.unicamp.br/jspui/handle/REPOSIP/306458.
Full textAbstract:
Orientador: Rodney Carlos BassaneziTese (doutorado) - Universidade Estadual de Campinas, Instituto de Matemática, Estatística e Computação Científica
Made available in DSpace on 2018-09-11T21:18:32Z (GMT). No. of bitstreams: 1 Leite_JeffersonCruzdosSantos_D.pdf: 39583926 bytes, checksum: ac69c5a564ed32a9d1eb58ac0e71c1fd (MD5) Previous issue date: 2011
Resumo: Neste trabalho definiremos solução fuzzy para problemas que envolvam difusão e exploraremos algumas propriedades importantes como unicidade e estabilidade dessas soluções. Basicamente estamos interessados em considerar algumas características importantes desses problemas difusivos como incertos, para isso, usaremos o conceito de numero fuzzy. Termos como coeficiente de difusão e condição inicial serão considerados como incertos e através da extensão de Zadeh aplicado a solução da equação determinística associada ao problema teremos a solução fuzzy. Serão obtidas também soluções via base de regras, utilizando sistemas dinâmicos pfuzzy, garantindo assim, uma maneira eficiente e prática de obtermos, boas respostas para os problemas, sem necessariamente termos as soluções explícitas. Aplicações desses resultados também serão apresentados
Abstract: This work will define fuzzy solution for problems involving di_usion and explore some important properties such as uniqueness and stability of these solutions. Basically we are interested in considering some important features of these diffusion problems as uncertain and, we use the concept of fuzzy numbers for this. Terms such as diffusion coefficient and initial condition are considered as uncertain and by the extension of Zadeh's solution applied to deterministic equation associated with the problem we have the fuzzy solution. Solutions for rule-base situations are also obtained, using p-fuzzy dynamic systems, thus guaranteeing an, efficient and practical way of obtaining adequate answers to the problems, not necessarily under the explicit solutions. Applications of these results will also be discussed
Doutorado
Matematica Aplicada
Doutor em Matemática Aplicada
27
Hepperger, Peter Thomas Verfasser], Claudia [Akademischer Betreuer] Klüppelberg, Rüdiger [Akademischer Betreuer] [Kiesel, and Fred Espen [Akademischer Betreuer] Benth. "Pricing and Hedging under High-Dimensional Jump-Diffusion Models using Partial Differential Equations / Peter Thomas Hepperger. Gutachter: Claudia Klüppelberg ; Rüdiger Kiesel ; Fred Espen Benth. Betreuer: Claudia Klüppelberg." München : Universitätsbibliothek der TU München, 2011. http://d-nb.info/1013436199/34.
Full text
28
Dimitry, Johan. "A mathematical study of convertible bonds." Thesis, KTH, Farkost och flyg, 2014. http://urn.kb.se/resolve?urn=urn:nbn:se:kth:diva-151312.
Full textAbstract:
A convertible bond (CB) is a financial derivative, a so called hybrid security. It is an issued contract from a company or a government, which is paid for up-front. The contract yields a known amount at the specified maturity date, unless the holder chooses to convert it into an amount of the underlying asset. This kind of financial products can have complex features affecting the contract price and the optimal exercising situation. The partial differential equation (PDE) approach used for pricing financial derivatives makes it possible to describe convertible bonds with a physical model, a reversed diffusion described by a parabolic PDE. One can sometimes find both analytical and numerical solutions for this type of PDEs and interpret the solutions from a financial point of view, as they suggest predictable behaviour of the contract price.
29
Chakrabarty, Nilaj. "Computational Study of Axonal Transport Mechanisms of Actin and Neurofilaments." Ohio University / OhioLINK, 2020. http://rave.ohiolink.edu/etdc/view?acc_num=ohiou1584441310326918.
Full text
30
Shabala, Alexander. "Mathematical modelling of oncolytic virotherapy." Thesis, University of Oxford, 2013. http://ora.ox.ac.uk/objects/uuid:cca2c9bc-cbd4-4651-9b59-8a4dea7245d1.
Full textAbstract:
This thesis is concerned with mathematical modelling of oncolytic virotherapy: the use of genetically modified viruses to selectively spread, replicate and destroy cancerous cells in solid tumours. Traditional spatially-dependent modelling approaches have previously assumed that virus spread is due to viral diffusion in solid tumours, and also neglect the time delay introduced by the lytic cycle for viral replication within host cells. A deterministic, age-structured reaction-diffusion model is developed for the spatially-dependent interactions of uninfected cells, infected cells and virus particles, with the spread of virus particles facilitated by infected cell motility and delay. Evidence of travelling wave behaviour is shown, and an asymptotic approximation for the wave speed is derived as a function of key parameters. Next, the same physical assumptions as in the continuum model are used to develop an equivalent discrete, probabilistic model for that is valid in the limit of low particle concentrations. This mesoscopic, compartment-based model is then validated against known test cases, and it is shown that the localised nature of infected cell bursts leads to inconsistencies between the discrete and continuum models. The qualitative behaviour of this stochastic model is then analysed for a range of key experimentally-controllable parameters. Two-dimensional simulations of in vivo and in vitro therapies are then analysed to determine the effects of virus burst size, length of lytic cycle, infected cell motility, and initial viral distribution on the wave speed, consistency of results and overall success of therapy. Finally, the experimental difficulty of measuring the effective motility of cells is addressed by considering effective medium approximations of diffusion through heterogeneous tumours. Considering an idealised tumour consisting of periodic obstacles in free space, a two-scale homogenisation technique is used to show the effects of obstacle shape on the effective diffusivity. A novel method for calculating the effective continuum behaviour of random walks on lattices is then developed for the limiting case where microscopic interactions are discrete.
31
Önskog, Thomas. "The Skorohod problem and weak approximation of stochastic differential equations in time-dependent domains." Doctoral thesis, Umeå universitet, Institutionen för matematik och matematisk statistik, 2009. http://urn.kb.se/resolve?urn=urn:nbn:se:umu:diva-25429.
Full textAbstract:
This thesis consists of a summary and four scientific articles. All four articles consider various aspects of stochastic differential equations and the purpose of the summary is to provide an introduction to this subject and to supply the notions required in order to fully understand the articles. In the first article we conduct a thorough study of the multi-dimensional Skorohod problem in time-dependent domains. In particular we prove the existence of cádlág solutions to the Skorohod problem with oblique reflection in time-independent domains with corners. We use this existence result to construct weak solutions to stochastic differential equations with oblique reflection in time-dependent domains. In the process of obtaining these results we also establish convergence results for sequences of solutions to the Skorohod problem and a number of estimates for solutions, with bounded jumps, to the Skorohod problem. The second article considers the problem of determining the sensitivities of a solution to a second order parabolic partial differential equation with respect to perturbations in the parameters of the equation. We derive an approximate representation of the sensitivities and an estimate of the discretization error arising in the sensitivity approximation. We apply these theoretical results to the problem of determining the sensitivities of the price of European swaptions in a LIBOR market model with respect to perturbations in the volatility structure (the so-called ‘Greeks’). The third article treats stopped diffusions in time-dependent graph domains with low regularity. We compare, numerically, the performance of one adaptive and three non-adaptive numerical methods with respect to order of convergence, efficiency and stability. In particular we investigate if the performance of the algorithms can be improved by a transformation which increases the regularity of the domain but, at the same time, reduces the regularity of the parameters of the diffusion. In the fourth article we use the existence results obtained in Article I to construct a projected Euler scheme for weak approximation of stochastic differential equations with oblique reflection in time-dependent domains. We prove theoretically that the order of convergence of the proposed algorithm is 1/2 and conduct numerical simulations which support this claim.
32
Schmidt, Daniel. "Kinetic Monte Carlo Methods for Computing First Capture Time Distributions in Models of Diffusive Absorption." Scholarship @ Claremont, 2017. https://scholarship.claremont.edu/hmc_theses/97.
Full textAbstract:
In this paper, we consider the capture dynamics of a particle undergoing a random walk above a sheet of absorbing traps. In particular, we seek to characterize the distribution in time from when the particle is released to when it is absorbed. This problem is motivated by the study of lymphocytes in the human blood stream; for a particle near the surface of a lymphocyte, how long will it take for the particle to be captured? We model this problem as a diffusive process with a mixture of reflecting and absorbing boundary conditions. The model is analyzed from two approaches. The first is a numerical simulation using a Kinetic Monte Carlo (KMC) method that exploits exact solutions to accelerate a particle-based simulation of the capture time. A notable advantage of KMC is that run time is independent of how far from the traps one begins. We compare our results to the second approach, which is asymptotic approximations of the FPT distribution for particles that start far from the traps. Our goal is to validate the efficacy of homogenizing the surface boundary conditions, replacing the reflecting (Neumann) and absorbing (Dirichlet) boundary conditions with a mixed (Robin) boundary condition.
33
Souza, Juliana Marta Rodrigues de 1985. "Estudo da dispersão de risco de epizootias em animais = o caso da influenza aviária." [s.n.], 2010. http://repositorio.unicamp.br/jspui/handle/REPOSIP/307276.
Full textAbstract:
Orientador: João Frederico da Costa Azevedo MeyerDissertação (mestrado) - Universidade Estadual de Campinas, Instituto de Matemática, Estatística e Computação Científica
Made available in DSpace on 2018-08-15T23:20:24Z (GMT). No. of bitstreams: 1 Souza_JulianaMartaRodriguesde_M.pdf: 3448446 bytes, checksum: c0a56c82b26926f022b1fbbb4b9e4fbe (MD5) Previous issue date: 2010
Resumo: Esta dissertação de mestrado do grupo de biomatemática do Instituto de Matemática Aplicada e Computacional da UNICAMP, com auxílio de Bolsa de mestrado da CNPq, é resultado de dois anos, 2008 e 2009, de estudo a respeito da dispersão do risco de contágio do H5N1. Após tratar brevemente da estrutura viral; do papel das aves que sofrem sua ação; dos problemas financeiros que o H5N1 traria ao Brasil e já inflingiu em outras nações; o trabalho concentra-se em modelar e simular um ambiente formado de duas populações de comportamento distinto. A primeira, de aves silvestre, livres, que podem migrar. A segunda população consiste de aves restritas ao controle de um criador; não voam, não se espalham além dos limites da pequena localidade onde são criadas para fins de subsistência. Cada uma das três subdivisões destas populações, de acordo com o status em relação à doença, é modelada por uma equação diferencial parcial, compondo um sistema cuja solução numérica, necessária por conta das descontinuidades das condições iniciais, prediz o comportamentos da infecção em função do tempo e do espaço. Dentre os resultados alcançados, destaca-se: o homem parece ter chance de conter o espalhamento do vírus. Para isso teria de sacrificar todos os animais de pequenas criações e, então indivíduos da população silvestre, mas a uma taxa menor do que eles são capazes de se reproduzir, ou seriam levados a extinção. Também estão contidos neste trabalho, o estudo dos estados estacionários do sistema e a estimativa de que o coeficiente de difusão do H5N1 assumiria valores entre 0,025 e 0,5 km²/dia
Abstract: This dissertation from the IMECC, UNICAMP, Biomathematical Group, with funds offered by CNPq, is the result of two years, 2008 and 2009, of study about the spreading of H5N1 risk of infection. After treating briefly the viral structure; the birds that suffer the virus; the financial problems that the disease would bring to Brazil and has already inflicted to other nations; this paper concentrates in modeling and simulating an environment composed by two distinct behaviour population. The first one is free wild birds, that migrate. The second population consists of birds restricted to a farmer control; they don't fly, don't spread beyond little farms limits where they are raised to subsistence purposes. After dividing each of these two populations in order three, acording to their status in relation to the H5N1 infection, they are modeled by means of Partial Differential Equation, composing a non-linear system which requires numerical solution because of descontinuous inicial conditions and predicts the infection behaviour in spatial and temporal terms. Among the results figure: Humans, by completely sacrifing small farms birds and, then, wild birds in smaller rate than they reproduce themselves, seems to have a chance of prevent the virus to spread even further. This paper also study stationary states and determine, through computational methods, the H5N1 coefficient range, among 0.025 and 0.5 km²/day
Mestrado
Biomatematica
Mestre em Matemática Aplicada
34
Roberts, Paul Allen. "Mathematical models of the retina in health and disease." Thesis, University of Oxford, 2015. http://ora.ox.ac.uk/objects/uuid:385f61c4-4ff1-45d3-bdb2-41338c174025.
Full textAbstract:
The retina is the ocular tissue responsible for the detection of light. Its extensive demand for oxygen, coupled with a concomitant elevated supply, renders this tissue prone to both hypoxia and hyperoxia. In this thesis, we construct mathematical models of the retina, formulated as systems of reaction-diffusion equations, investigating its oxygen-related dynamics in healthy and diseased states. In the healthy state, we model the oxygen distribution across the human retina, examining the efficacy of the protein neuroglobin in the prevention of hypoxia. It has been suggested that neuroglobin could prevent hypoxia, either by transporting oxygen from regions where it is rich to those where it is poor, or by storing oxygen during periods of diminished supply or increased uptake. Numerical solutions demonstrate that neuroglobin may be effective in preventing or alleviating hypoxia via oxygen transport, but that its capacity for oxygen storage is essentially negligible, whilst asymptotic analysis reveals that, contrary to the prevailing assumption, neuroglobin's oxygen affinity is near optimal for oxygen transport. A further asymptotic analysis justifies the common approximation of a piecewise constant oxygen uptake across the retina, placing existing models upon a stronger theoretical foundation. In the diseased state, we explore the effect of hyperoxia upon the progression of the inherited retinal diseases, known collectively as retinitis pigmentosa. Both numerical solutions and asymptotic analyses show that this mechanism may replicate many of the patterns of retinal degeneration seen in vivo, but that others are inaccessible to it, demonstrating both the strengths and weaknesses of the oxygen toxicity hypothesis. It is shown that the wave speed of hyperoxic degeneration is negatively correlated with the local photoreceptor density, high density regions acting as a barrier to the spread of photoreceptor loss. The effects of capillary degeneration and treatment with antioxidants or trophic factors are also investigated, demonstrating that each has the potential to delay, halt or partially reverse photoreceptor loss. In addition to answering questions that are not accessible to experimental investigation, these models generate a number of experimentally testable predictions, forming the first loop in what has the potential to be a fruitful experimental/modelling cycle.
35
Andreevska, Irena. "Mathematical modeling and analysis of options with jump-diffusion volatility." [Tampa, Fla.] : University of South Florida, 2008. http://purl.fcla.edu/usf/dc/et/SFE0002343.
Full text
36
Elias, Jan. "Modélisation mathématique du rôle et de la dynamique temporelle de la protéine p53 après dommages à l'ADN induits par les médicaments anticancéreux." Thesis, Paris 6, 2015. http://www.theses.fr/2015PA066253/document.
Full textAbstract:
Plusieurs modèles pharmacocinétiques-pharmacodynamiques moléculaires ont été proposés au cours des dernières décennies afin de représenter et de prédire les effets d'un médicament dans les chimiothérapies anticancéreuses. La plupart de ces modèles ont été développés au niveau de la population de cellules, puisque des effets mesurables peuvent y être observés beaucoup plus facilement que dans les cellules individuelles.Cependant, les véritables cibles moléculaires des médicaments se trouvent au niveau de la cellule isolée. Les médicaments utilisés soit perturbent l'intégrité du génome en provoquant des ruptures de brins de l'ADN et par conséquent initialisent la mort cellulaire programmée (apoptose), soit bloquent la prolifération cellulaire, par inhibition des protéines (cdks) qui permettent aux cellules de procéder d'une phase du cycle cellulaire à la suivante en passant par des points de contrôle (principalement en $G_1/S$ et $G_2/M$). Les dommages à l'ADN causés par les médicaments cytotoxiques ou la $\gamma$-irradiation activent, entre autres, les voies de signalisation contrôlées par la protéine p53 qui forcent directement ou indirectement la cellule à choisir entre la survie et la mort.Cette thèse vise à explorer en détail les voies intracellulaires impliquant la protéine p53, ``le gardien du génome", qui sont initiées par des lésions de l'ADN, et donc de fournir un rationnel aux cancérologues pour prédire et optimiser les effets des médicaments anticancéreux en clinique. Elle décrit l'activation et la régulation de la protéine p53 dans les cellules individuelles après leur exposition à des agents causant des dommages à l'ADN. On montre que les comportements dynamiques qui ont été observés dans les cellules individuelles peuvent être reconstruits et prédits par fragmentation des événements cellulaires survenant après lésion de l'ADN, soit dans le noyau, soit dans le cytoplasme. Ceci est mis en œuvre par la description du réseau des protéines à l'aide d'équations différentielles ordinaires (EDO) et partielles (EDP) impliquant plusieurs agents dont les protéines ATM, p53, Mdm2 et Wip1, dans le noyau aussi bien que dans le cytoplasme, et entre les deux compartiments. Un rôle positif de Mdm2 dans la synthèse de p53, qui a été récemment observé, est exploré et un nouveau mécanisme provoquant les oscillations de p53 est proposé. On pourra noter en particulier que le nouveau modèle rend compte d'observations expérimentales qui n'ont pas pu être entièrement expliquées par les modèles précédents, par exemple, l'excitabilité de p53.En utilisant des méthodes mathématiques, on observe de près la façon dont un stimulus (par exemple, une $\gamma$-irradiation ou des médicaments utilisés en chimiothérapie) est converti en un comportement dynamique spécifiques (spatio-temporel) de p53, en particulier que ces dynamiques spécifiques de p53, comme messager de l'information cellulaire, peuvent moduler le cycle de division cellulaire, par exemple provoquant l'arrêt du cycle ou l'apoptose. Des modèles mathématiques EDO et EDP de réaction-diffusion sont utilisés pour examiner comment le comportement (spatio-temporel) de p53 émerge, et nous discutons des conséquences de ce comportement sur les réseaux moléculaires, avec des applications possibles dans le traitement du cancer.Les interactions protéine-protéine sont considérées comme des réactions enzymatiques. On présente quelques résultats mathématiques pour les réactions enzymatiques, en particulier on étudie le comportement en temps grand du système de réaction-diffusion pour la réaction enzymatique réversible à l'aide d'une approche entropique. À notre connaissance, c'est la première fois qu'une telle étude est publiée sur ce sujetVarious molecular pharmacokinetic–pharmacodynamic models have been proposed in the last decades to represent and predict drug effects in anticancer therapies. Most of these models are cell population based models since clearly measurable effects of drugs can be seen on populations of (healthy and tumour) cells much more easily than in individual cells.The actual targets of drugs are, however, cells themselves. The drugs in use either disrupt genome integrity by causing DNA strand breaks and consequently initiate programmed cell death or block cell proliferation mainly by inhibiting proteins (cdks) that enable cells to proceed from one cell cycle phase to another. DNA damage caused by cytotoxic drugs or $\gamma$-irradiation activates, among others, the p53 protein-modulated signalling pathways that directly or indirectly force the cell to make a decision between survival and death.The thesis aims to explore closely intracellular pathways involving p53, ``the guardian of the genome", initiated by DNA damage and thus to provide oncologists with a rationale to predict and optimise the effects of anticancer drugs in the clinic. It describes p53 activation and regulation in single cells following their exposure to DNA damaging agents. We show that dynamical patterns that have been observed in individual cells can be reconstructed and predicted by compartmentalisation of cellular events occurring either in the nucleus or in the cytoplasm, and by describing protein interactions, using both ordinary and partial differential equations, among several key antagonists including ATM, p53, Mdm2 and Wip1, in each compartment and in between them. Recently observed positive role of Mdm2 in the synthesis of p53 is explored and a novel mechanism triggering oscillations is proposed. For example, new model can explain experimental observations that previous (not only our) models could not, e.g., excitability of p53.Using mathematical methods we look closely on how a stimulus (e.g., $\gamma$-radiation or drugs used in chemotherapy) is converted to a specific (spatio-temporal) pattern of p53 whereas such specific p53 dynamics as a transmitter of cellular information can modulate cellular outcomes, e.g., cell cycle arrest or apoptosis. Mathematical ODE and reaction-diffusion PDE models are thus used to see how the (spatio-temporal) behaviour of p53 is shaped and what possible applications in cancer treatment this behaviour might have. Protein-protein interactions are considered as enzyme reactions. We present some mathematical results for enzyme reactions, among them the large-time behaviour of the reaction-diffusion system for the reversible enzyme reaction treated by an entropy approach. To our best knowledge this is published for the first time
37
Jankowiak, Gaspard. "Etude asymptotique d'équations aux dérivées partielles de type diffusion non linéaire et inégalités fonctionnelles associées." Phd thesis, Université Paris Dauphine - Paris IX, 2014. http://tel.archives-ouvertes.fr/tel-01067226.
Full textAbstract:
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).
38
Seck, Ousmane. "Sur un modèle de diffusion non linéaire en dynamique des populations." Nancy 1, 1986. http://www.theses.fr/1986NAN10162.
Full textAbstract:
Les équations d'évolution. Étude dans L**(1) de solutions particulières. Les résultats d'existence connus. Estimations à priori. Relaxation des hypothèses de régularité et de compatibilité. Relaxation des hypothèses de stricte positivité
39
Rakotonasy, Solonjaka Hiarintsoa. "Modèle fractionnaire pour la sous-diffusion : version stochastique et edp." Phd thesis, Université d'Avignon, 2012. http://tel.archives-ouvertes.fr/tel-00839892.
Full textAbstract:
Ce travail a pour but de proposer des outils visant 'a comparer des résultats exp'erimentaux avec des modèles pour la dispersion de traceur en milieu poreux, dans le cadre de la dispersion anormale.Le "Mobile Immobile Model" (MIM) a été à l'origine d'importants progrès dans la description du transport en milieu poreux, surtout dans les milieux naturels. Ce modèle généralise l'quation d'advection-dispersion (ADE) e nsupposant que les particules de fluide, comme de solut'e, peuvent ˆetre immo-bilis'ees (en relation avec la matrice solide) puis relˆachées, le piégeage et le relargage suivant de plus une cin'etique d'ordre un. Récemment, une version stochastique de ce modèle a 'eté proposée. Malgré de nombreux succès pendant plus de trois décades, le MIM reste incapable de repr'esenter l''evolutionde la concentration d'un traceur dans certains milieux poreux insaturés. Eneffet, on observe souvent que la concentration peut d'ecroˆıtre comme unepuissance du temps, en particulier aux grands temps. Ceci est incompatible avec la version originale du MIM. En supposant une cinétique de piégeage-relargage diff'erente, certains auteurs ont propos'e une version fractionnaire,le "fractal MIM" (fMIM). C'est une classe d''equations aux d'eriv'ees par-tielles (e.d.p.) qui ont la particularit'e de contenir un op'erateur int'egral li'e'a la variable temps. Les solutions de cette classe d'e.d.p. se comportentasymptotiquement comme des puissances du temps, comme d'ailleurs cellesde l''equation de Fokker-Planck fractionnaire (FFPE). Notre travail fait partie d'un projet incluant des exp'eriences de tra¸cageet de vélocimétrie par R'esistance Magn'etique Nucl'eaire (RMN) en milieuporeux insatur'e. Comme le MIM, le fMIM fait partie des mod'eles ser-vant 'a interpréter de telles exp'eriences. Sa version "e.d.p." est adapt'eeaux grandeurs mesur'ees lors d'exp'eriences de tra¸cage, mais est peu utile pour la vélocimétrie RMN. En effet, cette technique mesure la statistiquedes d'eplacements des mol'ecules excit'ees, entre deux instants fixés. Plus précisément, elle mesure la fonction caractéristique (transform'ee de Fourier) de ces d'eplacements. Notre travail propose un outil d'analyse pour ces expériences: il s'agit d'une expression exacte de la fonction caract'eristiquedes d'eplacements de la version stochastique du mod'ele fMIM, sans oublier les MIM et FFPE. Ces processus sont obtenus 'a partir du mouvement Brown-ien (plus un terme convectif) par des changement de temps aléatoires. Ondit aussi que ces processus sont des mouvement Browniens, subordonnéspar des changements de temps qui sont eux-mˆeme les inverses de processusde L'evy non d'ecroissants (les subordinateurs). Les subordinateurs associés aux modèles fMIM et FFPE sont des processus stables, les subordinateursassoci'es au MIM sont des processus de Poisson composites. Des résultatsexp'erimenatux tr'es r'ecents on sugg'er'e d''elargir ceci 'a des vols de L'evy (plusg'en'eraux que le mouvement Brownien) subordonnés aussi.Le lien entre les e.d.p. fractionnaires et les mod'eles stochastiques pourla sous-diffusion a fait l'objet de nombreux travaux. Nous contribuons 'ad'etailler ce lien en faisant apparaˆıtre les flux de solut'e, en insistant sur une situation peu 'etudiée: nous examinons le cas o'u la cinétique de piégeage-relargage n'est pas la mˆeme dans tout le milieu. En supposant deux cinétiques diff'erentes dans deux sous-domaines, nous obtenons une version du fMIMavec un opérateur intégro-diff'erentiel li'e au temps, mais dépendant de la position.Ces r'esultats sont obtenus au moyen de raisonnements, et sont illustrés par des simulations utilisant la discrétisation d'intégrales fractionnaires etd'e.d.p. ainsi que la méthode de Monte Carlo. Ces simulations sont en quelque sorte des preuves numériques. Les outils sur lesquels elles s'appuient sont présentés aussi.
40
Landry, Richard S. Jr. "An Application of M-matrices to Preserve Bounded Positive Solutions to the Evolution Equations of Biofilm Models." ScholarWorks@UNO, 2017. https://scholarworks.uno.edu/td/2418.
Full textAbstract:
In this work, we design a linear, two step implicit finite difference method to approximate the solutions of a biological system that describes the interaction between a microbial colony and a surrounding substrate. Three separate models are analyzed, all of which can be described as systems of partial differential equations (PDE)s with nonlinear diffusion and reaction, where the biological colony grows and decays based on the substrate bioavailability. The systems under investigation are all complex models describing the dynamics of biological films. In view of the difficulties to calculate analytical solutions of the models, we design here a numerical technique to consistently approximate the system evolution dynamics, guaranteeing that nonnegative initial conditions will evolve uniquely into new, nonnegative approximations. This property of our technique is established using the theory of M-matrices, which are nonsingular matrices where all the entries of their inverses are positive numbers. We provide numerical simulations to evince the preservation of the nonnegative character of solutions under homogeneous Dirichlet and Neumann boundary conditions. The computational results suggest that the method proposed in this work is stable, and that it also preserves the bounded character of the discrete solutions.
41
Leguebe, Michael. "Modélisation de l'électroperméabilisation à l'échelle cellulaire." Thesis, Bordeaux, 2014. http://www.theses.fr/2014BORD0169/document.
Full textAbstract:
La perméabilisation des cellules à l’aide d’impulsions électriques intenses, appelée électroperméabilisation, est un phénomène biologique impliqué dans des thérapies anticancéreuses récentes. Elle permet, par exemple, d’améliorer l’efficacité d’une chimiothérapie en diminuant les effets secondaires, d’effectuer des transferts de gènes, ou encore de procéder à l’ablation de tumeurs. Les mécanismes de l’électroperméabilisation restent cependant encore méconnus, et l’hypothèse majoritairement admise par la communauté de formation de pores à la surface des membranes cellulaires est en contradiction avec certains résultats expérimentaux.Le travail de modélisation proposé dans cette thèse est basé sur une approche différente des modèles d’électroporation existants. Au lieu de proposer des lois sur les propriétés des membranes à partir d’hypothèses à l’échelle moléculaire, nous établissons des lois ad hoc pour les décrire, en se basant uniquement sur les informations expérimentales disponibles. Aussi, afin de rester au plus prèsde ces dernières et faciliter la phase de calibration à venir, nous avons ajouté un modèle de transport et de diffusion de molécules dans la cellule. Une autre spécificité de notre modèle est que nous faisons la distinction entre l’état conducteur et l’état perméable des membranes.Des méthodes numériques spécifiques ainsi qu’un code en 3D et parallèle en C++ ont été écrits et validés pour résoudre les équations aux dérivées partielles de ces différents modèles. Nous validons le travail de modélisation en montrant que les simulations reproduisent qualitativement les comportements observés in vitroCell permeabilization by intense electric pulses, called electropermeabilization, is a biological phenomenon involved in recent anticancer therapies. It allows, for example, to increase the efficacy of chemotherapies still reducing their side effects, to improve gene transfer, or to proceed tumor ablation. However, mechanisms of electropermeabilization are not clearly explained yet, and the mostly adopted hypothesis of the formation of pores at the membrane surface is in contradiction with several experimental results.This thesis modeling work is based on a different approach than existing electroporation models. Instead of deriving equations on membranes properties from hypothesis at the molecular scale, we prefer to write ad hoc laws to describe them, based on available experimental data only. Moreover, to be as close as possible to these data, and to ease the forthcoming work of parameter calibration, we added to our model equations of transport and diffusion of molecules in the cell. Another important feature of our model is that we differentiate the conductive state of membranes from their permeable state.Numerical methods, as well as a 3D parallel C++ code were written and validated in order to solve the partial differential equations of our models. The modeling work was validated by showing qualitative match between our simulations and the behaviours that are observed in vitro
42
Veruete, Mario. "Étude d'équations de réplication-mutation non locales en dynamique évolutive." Thesis, Montpellier, 2019. http://www.theses.fr/2019MONTS012/document.
Full textAbstract:
Nous analysons trois modèles non-locaux décrivant la dynamique évolutive d’un trait phénotypique continu soumis à l’action conjointe des mutations et de la sélection. Nous établissons l’existence et l’unicité des solutions du problème de Cauchy, et donnons la description du comportement en temps long de la solution. Dans le premier travail nous étudions l’équation du réplicateur-mutateur en domaine non borné et généralisons aux cas des valeurs sélectives confinantes les résultats connus dans le cas harmonique. À savoir, l’existence d’une unique solution globale, régulière, convergeant en temps long vers un profil universel ; pour cela, nous employons des techniques de décomposition spectrale d’opérateurs de Schrödinger. Le deuxième travail traite d’un modèle dont la valeur sélective est densité-dépendante. Afin de montrer le caractère bien posé de l’équation, nous combinons deux approches. La première est basée sur l’étude de la fonction génératrice des cumulants, satisfaisant une équation de transport non locale et permettant d’obtenir implicitement le trait moyen. La deuxième exploite un changement de variable (formule d’Avron-Herbst), permettant d’écrire la solution en termes du trait moyen et de la solution de l’équation de la chaleur avec même donnée initiale. Finalement, nous étudions un modèle dont le taux de mutation est proportionnel à la valeur moyenne du trait. Nous établissons un lien bijectif entre ce dernier modèle et le deuxième, permettant ainsi de décrire finement la dynamique de la solution. Nous montrons en particulier la croissance exponentielle du trait moyenWe analyze three non-local models describing the evolutionary dynamics of a continuous phenotypic trait undergoing the joint action of mutations and selection. We establish the existence and uniqueness of the solutions to the Cauchy problem, and give a description of the long-time behaviour of the solution. In the first work we study the replicator-mutator equation in the unbounded domain and generalize to cases of selective values confining the known results in the harmonic case. Namely, the existence of a unique global regular solution, converging towards a universal profile; for this, we use spectral decomposition techniques of Schrödinger operators. In the second work, we discuss a model whose fitness value is density-dependent. In order to show the well-posedness of the equation, we combine two approaches. The first is based on the study of the cumulant generating functions, satisfying a non-local transport equation and making it possible to implicitly obtain the average trait. The second uses a change of variable (Avron-Herbst formula), allowing the solution to be written in terms of the average trait and the solution of the heat equation with the same initial data. Finally, we study a model whose mutation rate is proportional to the average value of the trait. We establish a bijective link between this last model and the second, thus making it possible to describe the dynamics of the solution in detail. In particular, we show the exponential growth of the average trait
43
Laamri, El Haj. "Existence globale pour des systèmes de réaction-diffusion dans L**(1)." Nancy 1, 1988. http://www.theses.fr/1988NAN10164.
Full textAbstract:
On utilise des techniques L**(1) pour étudier l'existence globale en temps pour des systèmes de réaction-diffusion ou les deux propriétés suivantes sont satisfaites : - La positivité des solutions est préservée au cours du temps - la masse totate est préservée ou contrôlée au cours du temps. Ces systèmes proviennent de la modélisation des phénomènes chimiques, biologiques, écologiques, etc. . . Les idées de base relèvent de la théorie des opérateurs M-accrétifs dans L**(1) à résolvantes compactes
44
Stonkutė, Alina. "KINTAMO DIFUZIJOS KOEFICIENTO PARABOLINIŲ LYGČIŲ SPRENDIMAS SKAITINIAIS METODAIS." Master's thesis, Lithuanian Academic Libraries Network (LABT), 2010. http://vddb.laba.lt/obj/LT-eLABa-0001:E.02~2010~D_20100903_130346-27452.
Full textAbstract:
Magistro darbe sprendėme diferencialinę difuzijos lygtį naujais metodais. Išanalizavę standartinius kintamo difuzijos koeficiento parabolinių lygčių sprendimo metodus, mes šiame darbe pasiūlėme spręsti šias lygtis naudojant vadinamąsias „tilto“ funkcijas. Išbandėme dviejų rūšių „tilto“ funkcijas: hiperbolinio tangento ir trigonometrinio. Diferencialinės lygties sprendinio ieškojome per „tilto“ funkcijų ir polinomų sandaugų sumą: trigonometrinei „tilto“ funkcijai ir hiperbolinei tangento „tilto“ funkcijai. Gavome kompiuterinius sprendinius ir nustatėme tų sprendinių paklaidas. Palyginę trigonometrinio bei hiperbolinio tangento „tilto“ funkcijos paklaidų standartinius nuokrypius gavome, kad tikslesnis yra hiperbolinio tangento „tilto“ funkcijos metodas.Master thesis solved differential equation of diffusion of new techniques methods. Having analyzed the standard variable diffusion coefficient parabolic equation solution methods suggested in this work we solve these equations using the so-called "bridge" function. Tried two types of "bridge" functions: tangent hyperbolic and trigonometric. Differential equation, the solution we were looking for a "bridge" function and the amount of products of powers of polynomials: trigonometry "bridge" function and hyperbolic tangent of a "bridge" function. We have received computer-based solutions and the solutions found at the margins. A comparison of hyperbolic tangent trigonometric "bridge" function of the error standard deviations have received, the more accurate the hyperbolic tangent of a "bridge" function approach.
45
Dvořák, Radim. "Fyzikální modelování a simulace." Doctoral thesis, Vysoké učení technické v Brně. Fakulta informačních technologií, 2014. http://www.nusl.cz/ntk/nusl-261245.
Full textAbstract:
Disertační práce se zabývá modelováním znečištění ovzduší, jeho transportních a disperzních procesů ve spodní části atmosféry a zejména numerickými metodami, které slouží k řešení těchto modelů. Modelování znečištění ovzduší je velmi důležité pro předpověď kontaminace a pomáhá porozumět samotnému procesu a eliminaci následků. Hlavním tématem práce jsou metody pro řešení modelů popsaných parciálními diferenciálními rovnicemi, přesněji advekčně-difúzní rovnicí. Polovina práce je zaměřena na známou metodu přímek a je zde ukázáno, že tato metoda je vhodná k řešení určitých konkrétních problémů. Dále bylo navrženo a otestováno řešení paralelizace metody přímek, jež ukazuje, že metoda má velký potenciál pro akceleraci na současných grafických kartách a tím pádem i zvětšení přesnosti výpočtu. Druhá polovina práce se zabývá poměrně mladou metodou ELLAM a její aplikací pro řešení atmosférických advekčně-difúzních rovnic. Byla otestována konkrétní forma metody ELLAM společně s navrženými adaptacemi. Z výsledků je zřejmé, že v mnoha případech ELLAM překonává současné používané metody.
46
Mourad, Firas. "Estimation par méthodes inverses des paramètres de glissement et de diffusion des calottes glaciaires d'Antarctique." Thesis, Université Grenoble Alpes, 2020. http://www.theses.fr/2020GRALT021.
Full textAbstract:
Les modèles décrivant certains phénomènes naturels peuvent dépendre de paramètres non mesurables, d'où la nécessité de les estimer par méthodes inverses. Notre objectif est d'utiliser de telles techniques pour permettre une meilleure initialisation des modèles de simulations des calottes glaciaires en Antarctique. Cela permettra l'obtention de meilleures prévisions dans le cadre des études climatiques. Nous nous intéressons au paramètre de glissement basale qui caractérise le contact de la calotte glaciaire avec le socle rocheux. De même qu'au paramètre de diffusion qui dicte la dynamique au sein de l'équation différentielle partielle de continuité de masse décrivant son mouvement. Une approche basée sur la théorie de Lyapunov est proposée pour contrôler la convergence des modèles de transport inhomogènes 1D et 2D, vers un équilibre correspondant aux mesures de la topographie de surface de la calotte glaciaire de l'Antarctique. Notre travail propose une nouvelle loi pour l'inversion en 1D du coefficient de glissement basal. Nous utilisons également l'inversion adaptative de paramètres distribués pour récupérer le glissement basal depuis le paramètre de diffusion dans des modèles 1D et 2D. Ces deux méthodes sont testées sur des cas d'études et des données réelles. Nos résultats montrent que les méthodes proposées réussissent à inverser les paramètres de glissement et de diffusion tout en reproduisant les données disponiblesModels describing natural phenomena can depend on parameters that cannot be directly measured, hence the necessity to develop inverse techniques to determine them. Our goal is to utilize such techniques to enable better initialization of ice sheet models for Antarctica. This will help such models to produce better forecasts as part of climate studies. The parameters of interest are the basal sliding coefficient, which characterizes the contact of the ice sheet with the bed underneath, and the diffusion coefficient which dictates the dynamics within the mass-continuity partial differential equation describing the movement of ice sheets. A Lyapunov based approach is proposed to control the convergence of the 1D and 2D inhomogeneous transport models toward a feasible equilibrium matching the measurements of surface topography of the Antarctic ice sheet. Our work offers a new 1D update law for the basal sliding coefficient inversion. We also use adaptive distributed parameter inversion to retrieve basal sliding from diffusion in 1D and 2D models. These two methods are tested on study cases and real data. Our results show that the methods proposed are successful in inverting for sliding and diffusion while replicating the available data
47
Alshammari, Abdullah A. A. M. F. "Mathematical modelling of oxygen transport in skeletal and cardiac muscles." Thesis, University of Oxford, 2014. http://ora.ox.ac.uk/objects/uuid:65a34cb0-ef00-44c9-a04d-4147844c76ac.
Full textAbstract:
Understanding and characterising the diffusive transport of capillary oxygen and nutrients in striated muscles is key to assessing angiogenesis and investigating the efficacy of experimental and therapeutic interventions for numerous pathological conditions, such as chronic ischaemia. In articular, the influence of both muscle tissue and microvascular heterogeneities on capillary oxygen supply is poorly understood. The objective of this thesis is to develop mathematical and computational modelling frameworks for the purpose of extending and generalising the current use of histology in estimating the regions of tissue supplied by individual capillaries to facilitate the exploration of functional capillary oxygen supply in striated muscles. In particular, we aim to investigate the balance between local capillary supply of oxygen and oxygen demand in the presence of various anatomical and functional heterogeneities, by capturing tissue details from histological imaging and estimating or predicting regions of capillary supply. Our computational method throughout is based on a finite element framework that captures the anatomical details of tissue cross sections. In Chapter 1 we introduce the problem. In Chapter 2 we develop a theoretical model to describe oxygen transport from capillaries to uniform muscle tissues (e.g. cardiac muscle). Transport is then explored in terms of oxygen levels and capillary supply regions. In Chapter 3 we extend this modelling framework to explore the influence of the surrounding tissue by accounting for the spatial anisotropies of fibre oxygen demand and diffusivity and the heterogeneity in fibre size and shape, as exemplified by mixed muscle tissues (e.g. skeletal muscle). We additionally explore the effects of diffusion through the interstitium, facilitated--diffusion by myoglobin, and Michaelis--Menten kinetics of tissue oxygen consumption. In Chapter 4, a further extension is pursued to account for intracellular heterogeneities in mitochondrial distribution and diffusive parameters. As a demonstration of the potential of the models derived in Chapters 2--4, in Chapter 5 we simulate oxygen transport in myocardial tissue biopsies from rats with either impaired angiogenesis or impaired arteriolar perfusion. Quantitative predictions are made to help explain and support experimental measurements of cardiac performance and metabolism. In the final chapter we summarize the main results and indicate directions for further work.
48
Rodrigues, Daiana Aparecida. "Modelagem e solução numérica de equações reação-difusão em processos biológicos." Universidade Federal de Juiz de Fora, 2013. https://repositorio.ufjf.br/jspui/handle/ufjf/1153.
Full textAbstract:
Submitted by Renata Lopes (renatasil82@gmail.com) on 2016-04-11T19:27:27Z
No. of bitstreams: 1
daianaaparecidarodrigues.pdf: 8225936 bytes, checksum: 96ec323f343f92c319f4e261145f9c6a (MD5)Approved for entry into archive by Adriana Oliveira (adriana.oliveira@ufjf.edu.br) on 2016-04-24T03:34:16Z (GMT) No. of bitstreams: 1 daianaaparecidarodrigues.pdf: 8225936 bytes, checksum: 96ec323f343f92c319f4e261145f9c6a (MD5)
Made available in DSpace on 2016-04-24T03:34:16Z (GMT). No. of bitstreams: 1 daianaaparecidarodrigues.pdf: 8225936 bytes, checksum: 96ec323f343f92c319f4e261145f9c6a (MD5) Previous issue date: 2013-08-29
CAPES - Coordenação de Aperfeiçoamento de Pessoal de Nível Superior
Fenômenos biológicos são todo e qualquer evento que possa ser observado nos seres vivos. O estudo desses fenômenos permite propor explicações para o seu mecanismo, a m de entender as causas e efeitos. Pode-se citar como exemplos de fenômenos biológicos o comportamento das células como respiração, reprodução, metabolismo e morte celular. Equações de reação-difusão são frequentemente utilizadas para modelar fenômenos bioló- gicos. Sistemas de reação-difusão podem produzir padrões espaciais estáveis a partir de uma distribuição inicial uniforme esse fenômeno é conhecido como instabilidade de Turing. Este trabalho apresenta a análise da instabilidade de Turing bem como resultados numéricos para a solução de três modelos biológicos, modelo de Schnakenberg, modelo de glicólise e modelo da coagulação sanguínea. O modelo de Schnakenberg é utilizado para descrever uma reação química autocatalítica e o modelo de glicólise é relativo ao processo de degradação metabólica da molécula de glicose para proporcionar energia para o metabolismo celular, esses dois modelos são frequentemente relatados na literatura. O terceiro modelo é mais recente e descreve o fenômeno da coagulação sanguínea. Nas soluções numéricas se utiliza o método das linhas onde a discretização espacial é feita através de um esquema de diferenças nitas. O sistema de equações diferencias ordinárias resultante é resolvido por um esquema de integração adaptativo, com a utilização de pacote para computação cientí ca da linguagem Python, Scipy.
Biological phenomena are all and any event that can be observed in living beings. The study of these phenomena enables us to propose explanations for its mechanisms in order to understand causes and e ects. One can cite as examples of biological phenomena the behavior of cells as respiration, reproduction, metabolism and cell death. Reactiondi usion equations are often used to model biological phenomena. Reaction-di usion systems can produce stable spatial patterns from a uniform initial distribution, this phenomenon is known as Turing instability. This dissertation presents an analysis of the Turing instability as well as numerical results for the solution of three biological models, model Schnakenberg, model of glycolysis and model of blood coagulation. The Schnakenberg model is used to describe an autocatalytic chemical reaction and glycolysis model refers to the process of metabolic breakdown of the glucose molecule to provide energy for cellular metabolism, these two models are frequently reported in the literature. The third model is newer and describes the phenomenon of blood coagulation. The method of lines is used in the numerical solutions, where the spatial discretization is done through a nite di erence scheme. The resulting system of ordinary di erential equations is then solved by an adaptive integration scheme with the use of the package for scienti c computing of Python language, Scipy.
49
Dareiotis, Anastasios Constantinos. "Stochastic partial differential and integro-differential equations." Thesis, University of Edinburgh, 2015. http://hdl.handle.net/1842/14186.
Full textAbstract:
In this work we present some new results concerning stochastic partial differential and integro-differential equations (SPDEs and SPIDEs) that appear in non-linear filtering. We prove existence and uniqueness of solutions of SPIDEs, we give a comparison principle and we suggest an approximation scheme for the non-local integral operators. Regarding SPDEs, we use techniques motivated by the work of De Giorgi, Nash, and Moser, in order to derive global and local supremum estimates, and a weak Harnack inequality.
50
Ranner, Thomas. "Computational surface partial differential equations." Thesis, University of Warwick, 2013. http://wrap.warwick.ac.uk/57647/.
Full textAbstract:
Surface partial differential equations model several natural phenomena; for example in uid mechanics, cell biology and material science. The domain of the equations can often have complex and changing morphology. This implies analytic techniques are unavailable, hence numerical methods are required. The aim of this thesis is to design and analyse three methods for solving different problems with surface partial differential equations at their core. First, we define a new finite element method for numerically approximating solutions of partial differential equations in a bulk region coupled to surface partial differential equations posed on the boundary of this domain. The key idea is to take a polyhedral approximation of the bulk region consisting of a union of simplices, and to use piecewise polynomial boundary faces as an approximation of the surface and solve using isoparametric finite element spaces. We study this method in the context of a model elliptic problem. The main result in this chapter is an optimal order error estimate which is confirmed in numerical experiments. Second, we use the evolving surface finite element method to solve a Cahn- Hilliard equation on an evolving surface with prescribed velocity. We start by deriving the equation using a conservation law and appropriate transport formulae and provide the necessary functional analytic setting. The finite element method relies on evolving an initial triangulation by moving the nodes according to the prescribed velocity. We go on to show a rigorous well-posedness result for the continuous equations by showing convergence, along a subsequence, of the finite element scheme. We conclude the chapter by deriving error estimates and present various numerical examples. Finally, we stray from surface finite element method to consider new unfitted finite element methods for surface partial differential equations. The idea is to use a fixed bulk triangulation and approximate the surface using a discrete approximation of the distance function. We describe and analyse two methods using a sharp interface and narrow band approximation of the surface for a Poisson equation. Error estimates are described and numerical computations indicate very good convergence and stability properties.