Cette thèse est consacrée à l'étude de plusieurs problèmes d'équations aux dérivées partielles non-linéaires.<p><p>La première partie (chapitres 1-2-3) traite d'un problème trouvant son origine en biologie mathématique, à savoir l'étude de la survie à long terme d'une population dont l'évolution est gouvernée par une équation parabolique non-linéaire. Dans le modèle considéré, le mécanisme de diffusion est contrôlé par le p-Laplacien, la non-linéarité est de type logistique et fait intervenir un poids m pouvant changer de signe, et les conditions aux limites sont de flux nul. Le poids m correspond à une répartition des ressources devant permettre la survie de la population. Dans le chapitre 1, nous déterminons entre autres un critère de survie à long terme faisant intervenir la valeur propre principale du p-Laplacien avec poids m. Cette valeur propre apparait, plus précisément, comme la valeur limite d'un paramètre en-dessous de laquelle toute solution positive de l'équation converge vers zéro lorsque t tend vers l'infini. Ceci nous conduit naturellement au problème de minimiser la valeur propre en question lorsque m varie dans une classe adéquate de poids. Dans le chapitre 2, nous prouvons l'existence de minimiseurs et montrons que ces derniers satisfont une propriété de type “bang-bang”. Plusieurs propriétés de montonie sont aussi étudiées dans des situations géométriques particulières, et une caractérisation complète est donnée en dimension 1. Le chapitre 3 est consacré à l'élaboration de simulations numériques, où l'algorithme utilisé combine un méthode de plus grande pente avec une représentation de certains ensembles comme ensembles de niveaux.<p><p>La deuxième sujet de cette thèse (chapitre 4) est un problème elliptique faisant intervenir l'opérateur de courbure moyenne. Nous nous intéressons à l'existence et à la multiplicité de solutions nodales de ce problème. Nous montrons que, si un certain paramètre de l'équation est suffisamment grand, il existe une solution nodale qui change de signe exactement deux fois. Nous établissons également l'existence d'un nombre arbitrairement grand de solutions nodales. Enfin, dans le cas particulier où le domaine est une boule, un résultat de brisure de symétrie est obtenu, résultat qui induit l'existence d'au moins deux solutions à deux domaines nodaux. / Doctorat en Sciences / info:eu-repo/semantics/nonPublished

Wound healing is a complex process, in which a sequence of interrelated events at both the cell and tissue levels interact and contribute to the reduction in wound size. For diabetic patients, many of these processes are compromised, so that wound healing slows down and in some cases halts. In this thesis we develop a series of increasingly detailed mathematical models to describe and investigate healing of full thickness skin wounds. We begin by developing a time-dependent ordinary differential equation model. This phenomenological model focusses on the main processes contributing to closure of a full thickness wound: proliferation in the epidermis and growth and contraction in the dermis. Model simulations suggest that the relative contributions of growth and contraction to healing of the dermis are altered in diabetic wounds. We investigate further the balance between growth and contraction by developing a more detailed, spatially-resolved model using continuum mechanics. Due to the initial large retraction of the wound edge upon injury, we adopt a non-linear elastic framework. Morphoelasticity theory is applied, with the total deformation of the material decomposed into an addition of mass and an elastic response. We use the model to investigate how interactions between growth and stress influence dermal wound healing. The model reveals that contraction alone generates unrealistically high tension in the dermal tissue and, hence, volumetric growth must contribute to healing. We show that, in the simplified case of homogeneous growth, the tissue must grow anisotropically in order to reduce the size of the wound and we postulate mechanosensitive growth laws consistent with this result. After closure the surrounding tissue remodels, returning to its residually stressed state. We identify the steady state growth profile associated with this remodelled state. The model is used to predict the outcome of rewounding experiments as a method of quantifying the amount of stress in the tissue and the application of pressure treatments to control tissue synthesis. The thesis concludes with an extension to the spatially-resolved mechanical model to account for the effects of the biochemical environment. Partial differential equations describing the dynamics of fibroblasts and a regulating growth factor are coupled to equations for the tissue mechanics, described in the morphoelastic framework. By accounting for biomechanical and biochemical stimuli the model allows us to formulate mechanistic laws for growth and contraction. We explore how disruption of mechanical and chemical feedback can lead to abnormal wound healing and use the model to identify specific treatments for normalising healing in these cases.

This study seeks to develop a method that generalises the use of Adams-Bashforth to solve or treat partial differential equations with local and non-local differentiation by deriving a two-step Adams-Bashforth numerical scheme in Laplace space. The resulting solution is then transformed back into the real space by using the inverse Laplace transform. This is a powerful numerical algorithm for fractional order derivative. The error analysis for the method is studied and presented. The numerical simulations of the method as applied to the four-dimensional model, Caputo-Lu-Chen model and the wave equation are presented. In the analysis, the bifurcation dynamics are discussed and the periodic doubling processes that eventually caused chaotic behaviour (butterfly attractor) are shown. The related graphical simulations that show the existence of fractal structure that is characterised by chaos and usually called strange attractors are provided. For the Caputo-Lu-Chen model, graphical simulations have been realised in both integer and fractional derivative orders. / Mathematical Sciences / M. Sc. (Applied Mathematics)

Cell migration, the active translocation of cells is involved in various biological processes, e.g. development of tissues and organs, tumor invasion and wound healing. Cell migration behavior can be divided into two distinct classes: single cell migration and collective cell migration. Single cell migration describes the migration of cells without interaction with other cells in their environment. Collective cell migration is the joint, active movement of multiple cells, e.g. in the form of strands, cohorts or sheets which emerge as the result of individual cell-cell interactions. Collective cell migration can be observed during branching morphogenesis, vascular sprouting and embryogenesis. Experimental studies of single cell migration have been extensive. Collective cell migration is less well investigated due to more difficult experimental conditions than for single cell migration. Especially, experimentally identifying the impact of individual differences in cell phenotypes on individual cell migration behavior inside cell populations is challenging because the tracking of individual cell trajectories is required. In this thesis, a novel mathematical modeling approach, individual-based lattice-gas cellular automata (IB-LGCA), that allows to investigate the migratory behavior of individual cells inside migrating cell populations by enabling the tracking of individual cells is introduced. Additionally, stochastic differential equation (SDE) approximations of individual cell trajectories for IB-LGCA models are constructed. Such SDE approximations allow the analytical description of the trajectories of individual cells during single cell migration. For a complete analytical description of the trajectories of individual cell during collective cell migration the aforementioned SDE approximations alone are not sufficient. Analytical approximations of the time development of selected observables for the cell population have to be added. What observables have to be considered depends on the specific cell migration mechanisms that is to be modeled. Here, partial integro-differential equations (PIDE) that approximate the time evolution of the expected cell density distribution in IB-LGCA are constructed and coupled to SDE approximations of individual cell trajectories. Such coupled PIDE and SDE approximations provide an analytical description of the trajectories of individual cells in IB-LGCA with density-dependent cell-cell interactions. Finally, an IB-LGCA model and corresponding analytical approximations were applied to investigate the impact of changes in cell-cell and cell-ECM forces on the migration behavior of an individual, labeled cell inside a population of epithelial cells. Specifically, individual cell migration during the epithelial-mesenchymal transition (EMT) was considered. EMT is a change from epithelial to mesenchymal cell phenotype which is characterized by cells breaking adhesive bonds with surrounding epithelial cells and initiating individual migration along the extracellular matrix (ECM). During the EMT, a transition from collective to single cell migration occurs. EMT plays an important role during cancer progression, where it is believed to be linked to metastasis development. In the IB-LGCA model epithelial cells are characterized by balanced cell-cell and cell-ECM forces. The IB-LGCA model predicts that the balance between cell-cell and cell-ECM forces can be disturbed to some degree without being accompanied by a change in individual cell migration behavior. Only after the cell force balance has been strongly interrupted mesenchymal migration behavior is possible. The force threshold which separates epithelial and mesenchymal migration behavior in the IB-LGCA has been identified from the corresponding analytical approximation. The IB-LGCA model allows to obtain quantitative predictions about the role of cell forces during EMT which in the context of mathematical modeling of EMT is a novel approach.

Indiana University-Purdue University Indianapolis (IUPUI) / Neural signal processing to decode neural activity has been an active research area in the last few decades. The next generation of advanced multi-electrode neuroprosthetic devices aim to detect a multiplicity of channels from multiple electrodes, making the relatively time-critical processing problem massively parallel and pushing the computational demands beyond the limits of current embedded digital signal processing (DSP) techniques. To overcome these limitations, a new hybrid computational technique was explored, the Extended Analog Computer (EAC). The EAC is a digitally confgurable analog computer that takes advantage of the intrinsic ability of manifolds to solve partial diferential equations (PDEs). They are extremely fast, require little power, and have great potential for mobile computing applications. In this thesis, the EAC architecture and the mechanism of the formation of potential/current manifolds was derived and analyzed to capture its theoretical mode of operation. A new mode of operation, resistance mode, was developed and a method was devised to sample temporal data and allow their use on the EAC. The method was validated by demonstration of the device solving linear diferential equations and linear functions, and implementing arbitrary finite impulse response (FIR) and infinite impulse response (IIR) linear flters. These results were compared to conventional DSP results. A practical application to the neural computing task was further demonstrated by implementing a matched filter with the EAC simulator and the physical prototype to detect single fiber action potential from multiunit data streams derived from recorded raw electroneurograms. Exclusion error (type 1 error) and inclusion error (type 2 error) were calculated to evaluate the detection rate of the matched filter implemented on the EAC. The detection rates were found to be statistically equivalent to that from DSP simulations with exclusion and inclusion errors at 0% and 1%, respectively.

This thesis is concerned with the investigation of numerical methods for the solution of the Hamilton-Jacobi-Bellman (HJB) equations arising in European and American option pricing with proportional transaction costs. We first consider the problem of computing reservation purchase and write prices of a European option in the model proposed by Davis, Panas and Zariphopoulou [19]. It has been shown [19] that computing the reservation purchase and write prices of a European option involves solving three different fully nonlinear HJB equations. In this thesis, we propose a penalty approach combined with a finite difference scheme to solve the HJB equations. We first approximate each of the HJB equations by a quasi-linear second order partial differential equation containing two linear penalty terms with penalty parameters. We then develop a numerical scheme based on the finite differencing in both space and time for solving the penalized equation. We prove that there exists a unique viscosity solution to the penalized equation and the viscosity solution to the penalized equation converges to that of the original HJB equation as the penalty parameters tend to infinity. We also prove that the solution of the finite difference scheme converges to the viscosity solution of the penalized equation. Numerical results are given to demonstrate the effectiveness of the proposed method. We extend the penalty approach combined with a finite difference scheme to the HJB equations in the American option pricing model proposed by Davis and Zarphopoulou [20]. Numerical experiments are presented to illustrate the theoretical findings.

A equação de advecção-dispersão possui grande importância na engenharia e nas ciências aplicadas. No entanto, como é bem conhecido, a obtenção de uma solução numérica apropriada para essa equação é um problema desafiador tanto para engenheiros como para matemáticos, físicos e outros profissionais que trabalham com a modelagem de fenômenos associados a ela. Muitos métodos numéricos desenvolvidos podem apresentar uma série de inconvenientes, tais como oscilações, dispersão e/ou dissipação numérica e instabilidade, além de serem inapropriados para determinadas condições de contorno. O presente trabalho apresenta e analisa a metodologia DFLD-exp, uma nova abordagem para a obtenção de soluções semi-analíticas da equação de advecção-dispersão, a qual utiliza um tipo particular de diferenças finitas para a discretização espacial juntamente com técnicas de exponencial de matrizes para a resolução temporal. Uma cuidadosa análise numérica mostra que a metodologia resultante é não-oscilatória, essencialmente não-dispersiva e não-dissipativa, e incondicionalmente estável. Resoluções de vários exemplos numéricos, através de um código desenvolvido em linguagem MATLAB, confirmam os resultados teóricos. / The advection-dispersion equation has been very important in engineering and the applied sciences. However, the obtainment of an appropriate numerical solution to that equation has been challenging problem to engineers, mathematicians, physicians and others that work in the modeling of phenomena associate to advection-dispersion equation. Many developed numerical methods may produce a succession of mistakes, just as oscillations, numerical dispersion and/or dissipation, instability and those methods also may be inappropriate to determined boundary conditions. The present work shows and analyses the DFLD-exp methodology, a new way to obtain semi-analytic solutions to advection-dispersion equation, that make use of a particular form of finite differencing to the spatial discretization with techniques of matrix exponential to the time solving. A detailed numerical analysis shows the methodology is non-oscillatory, essentially non-dispersive and non-dissipative, and unconditionally stable. Resolutions of any numerical examples, by a computational code developed in MATLAB language, confirm the theoretical results.

A equação de advecção-dispersão possui grande importância na engenharia e nas ciências aplicadas. No entanto, como é bem conhecido, a obtenção de uma solução numérica apropriada para essa equação é um problema desafiador tanto para engenheiros como para matemáticos, físicos e outros profissionais que trabalham com a modelagem de fenômenos associados a ela. Muitos métodos numéricos desenvolvidos podem apresentar uma série de inconvenientes, tais como oscilações, dispersão e/ou dissipação numérica e instabilidade, além de serem inapropriados para determinadas condições de contorno. O presente trabalho apresenta e analisa a metodologia DFLD-exp, uma nova abordagem para a obtenção de soluções semi-analíticas da equação de advecção-dispersão, a qual utiliza um tipo particular de diferenças finitas para a discretização espacial juntamente com técnicas de exponencial de matrizes para a resolução temporal. Uma cuidadosa análise numérica mostra que a metodologia resultante é não-oscilatória, essencialmente não-dispersiva e não-dissipativa, e incondicionalmente estável. Resoluções de vários exemplos numéricos, através de um código desenvolvido em linguagem MATLAB, confirmam os resultados teóricos. / The advection-dispersion equation has been very important in engineering and the applied sciences. However, the obtainment of an appropriate numerical solution to that equation has been challenging problem to engineers, mathematicians, physicians and others that work in the modeling of phenomena associate to advection-dispersion equation. Many developed numerical methods may produce a succession of mistakes, just as oscillations, numerical dispersion and/or dissipation, instability and those methods also may be inappropriate to determined boundary conditions. The present work shows and analyses the DFLD-exp methodology, a new way to obtain semi-analytic solutions to advection-dispersion equation, that make use of a particular form of finite differencing to the spatial discretization with techniques of matrix exponential to the time solving. A detailed numerical analysis shows the methodology is non-oscillatory, essentially non-dispersive and non-dissipative, and unconditionally stable. Resolutions of any numerical examples, by a computational code developed in MATLAB language, confirm the theoretical results.

<p align="justify">Une approche multidimensionelle pour la résolution numérique des équations d'Euler et de Navier-Stokes sur maillages non-structurés est proposée. Dans une première partie, un exposé complet des schémas de distribution, dits de "fluctuation-splitting" ,est décrit, comprenant une étude comparative des schémas décentrés, positifs et de 2ème ordre, pour résoudre l'équation de convection à coefficients constants, ainsi qu'une étude théorique et numérique de la précision des schémas sur maillages réguliers et distordus. L'extension à des lois de conservation non-linéaires est aussi abordée, et une attention particulière est portée au problème de la linéarisation conservative. Dans une deuxième partie, diverses discrétisations des termes visqueux pour l'équation de convection-diffusion sont développées, avec pour but de déterminer l'approche qui offre le meilleur compromis entre précision et coût. L'extension de la méthode aux systèmes des lois de conservation, et en particulier à celui des équations d'Euler de la dynamique des gaz, représente le noyau principal de la thèse, et est abordée dans la troisième partie. Contrairement aux schémas de distribution classiques, qui reposent sur une extension formelle du cas scalaire, l'approche développée ici repose sur une décomposition du résidu par élément en équations scalaires, modélisant le transport de variables caracteristiques. La difficulté vient du fait que les équations d'Euler instationnaires ne se diagonalisent pas, et admettent une infinité de solutions élémentaires (ondes simples) se propageant dans toutes les directions d'espace. En régime stationnaire, en revanche, les équations se diagonalisent complètement dans le cas des écoulements supersoniques, et partiellement dans le cas des écoulements subsoniques. Ainsi, les équations sous forme conservative peuvent être remplacées par un système équivalent comprenant deux équations totalement découplées, exprimant l'invariance de l'entropie et de l'enthalpie totale le long des lignes de courant, et deux autres équations, modélisant les effets purement acoustiques. En régime supersonique, celles-ci se découplent aussi, et expriment la convection le long des lignes de Mach d'invariants de Riemann généralisés. La discrétisation de ces équations par des schémas scalaires décentrés permet de simuler des écoulements continus et discontinus avec une grande précision et sans oscillations. Finalement, dans une dernière partie, l'extension aux équations de Navier-Stokes est abordée, et la discrétisation des termes visqueux par une approche éléments finis est proposée. Les résultats numériques confirment la précision et la robustesse de la méthode.</p> / Doctorat en sciences appliquées / info:eu-repo/semantics/nonPublished

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.