Return to search

An Application of M-matrices to Preserve Bounded Positive Solutions to the Evolution Equations of Biofilm Models

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.

Identiferoai:union.ndltd.org:uno.edu/oai:scholarworks.uno.edu:td-3532
Date20 December 2017
CreatorsLandry, Richard S., Jr.
PublisherScholarWorks@UNO
Source SetsUniversity of New Orleans
Detected LanguageEnglish
Typetext
Formatapplication/pdf
SourceUniversity of New Orleans Theses and Dissertations

Page generated in 0.0016 seconds