Applied Mathematics and Mathematical Physics Seminar 2010/2011 (sponsored by MITACS and PIMS)

This is the web page for the Applied Mathematics and Mathematical Physics Seminar at the Department of Mathematics and Statistics, University of Saskatchewan, Saskatoon, Canada. This seminar is organized and jointly run by W. Abou Salem, J. Brooke, A. Cheviakov, G. Patrick, A. Sowa, J. Szmigielski (Applied Mathematics/ Mathematical Physics) and R. Spiteri (Department of Computer Science).

For inquiries concerning the seminar send email to Jacek Szmigielski, szmigiel@math.usask.ca.



The seminar takes place in McLean Hall rm 242.1 on Wednesdays at 3:00 till 4:00 unless advertised differently.


APPLIED MATHEMATICS/ MATHEMATICAL PHYSICS SEMINAR 2007-2008 (sponsored by MITACS)


APPLIED MATHEMATICS/ MATHEMATICAL PHYSICS SEMINAR 2009-2010 (sponsored by MITACS)

Next meetings: August 11, 2011,

Previous meetings : July 7, 2011, June 30, 2011, May 20, 2011, April 20, 2011, March 25, 2011, March 23, 2011, March 14, 2011, March 3, 2011, February 15, 2011, February 4, 2011, February 2, 2011, November 24, 2010, November 17, 2010, November 10, 2010, November 4, 2010, November 3, 2010, November 2, 2010, September 23, 2010, September 9, 2010,


September 9, 2010

Speaker: Dr. Oluwaseun Sharomi Department of Mathematics University of Manitoba

Place: McLean Hall 242.1

Time: 3:00

Title: Mathematical Analysis of dynamics of Chlamydia trachomatis

Abstract

Chlamydia, caused by the bacterium Chlamydia trachomatis, is one of the most important sexually-transmitted infections globally. In addition to accounting for millions of cases every year, the disease causes numerous irreversible complications such as chronic pelvic pain, infertility in females and pelvic inflammatory disease. The talk will focus on the use of mathematical models, of the form of deterministic systems of non-linear differential equations, for gaining qualitative insight into the transmission dynamics and control of Chlamydia within an infected host in vivo) and in a population.

About The Speaker: Dr. Oluwaseun Sharomi received his Ph.D. in Mathematics from the University of Manitoba in 2010. In 2002/2003, he was awarded the United Bank Africa Prize, the LYNX Club Abeokuta Prize, the Professor Ishola Adamson Prize, and the University Prize for the Best Graduating Student in the Department of Mathematical Sciences, University of Agriculture, Abeokuta, Ogun State, Nigeria. His research interests include: mathematical modeling, analysis, and numerical methods for the spread of infectious diseases.


Back to top of page

September 23, 2010

Speaker: Professor George Corliss Department of Electrical and Computer Engineering Marquette University

Place: McLean Hall 242.1

Time: 3:00

Title: Propagating Uncertainties in Modeling Nonlinear Dynamic Systems

Abstract

In many practical applications of ordinary differential equations modeling physical phenomina, parameters and initial conditions are given with some uncertainties. How can we propagate these uncertainties rigorously through a solution approximation algorithm? We describe an approach that expands the solution in a Taylor model [Makino & Berz] in uncertain parameters and initial conditions. We evaluate the Taylor models using p-boxes [Ferson] and gradual numbers representing fuzzy numbers to represent the uncertainties in the state variables of the ODE. We give examples from reaction process dynamics to demonstrate the potential of this approach for studying the effect of uncertainties with imprecise probability distributions.


Back to top of page

November 2, 2010

Speaker: Dr Giulio Chiribella Perimeter Institute for Theoretical Physics, Waterloo, Ontario b>

Place: Arts 134

Time: 4:00

Title: Optimal Estimation of Quantum Signals in the Presence of Symmetry

Abstract

Quantum systems can be used as elementary gyroscopes that indicate directions in space or as elementary clocks that indicate moments in time. However, the laws of quantum mechanics impose fundamental precision limits to the corresponding measurements of orientation and time, limits that cannot be violated no matter how advanced our technology is. Assessing the exact value of these limits is important for many applications in interferometry, magnetometry, GPS systems, and the study of quantum communication protocols where the communicating parties try to establish a common reference frame. Finding the most precise estimate for the direction of an ensemble of atomic spins or for the phase of a quantum oscillator are instances of a very general problem: the optimal estimation of parameters pertaining to the action of symmetry groups. In this talk I will review the framework of quantum estimation theory in the presence of symmetries and highlight the basic group theoretic structures that underlie the optimal estimation strategies. In particular, I will highlight the role of quantum entanglement between the representation spaces and the multiplicity spaces associated to the group action of interest.

About the speaker: Dr Chiribella was, just recently, awarded the Hermann Weyl Prize at the 28th International Colloquium on Group Theoretical Methods in Physics, Northumbria University, Newcastle upon Tyne, UK, July 2010.


Back to top of page

November 3, 2010

Speaker: Dr Giulio Chiribella Perimeter Institute for Theoretical Physics, Waterloo, Ontario b>

Place: McLean Hall 242.1

Time: 3:30

Title:Quantum measurement theory: from quantum states to quantum networks

Abstract

The main aspect in which quantum theory is richer than classical probability theory is that the quantum framework contains an explicit description of the measurement process, through the concepts of POVM (positive operator-valued measure) and quantum instrument. POVMs and quantum instruments provide the probability of outcomes and the input-output evolution of quantum states in a measurement, respectively. However, measuring quantum systems is not the only type of measurement-like operation one can perform: for example, one can try to measure properties of quantum devices, like the gain of an amplifier or the loss in an optical fibre. In this talk I will review the theoretical description and the physical implementation of these general measurement processes, which has been provided very recently with the theory of quantum combs and testers [1]. As examples of application of the theory, I will present the optimal architectures for the estimation of an unknown group transformation, and the optimal estimation-disturbance trade-off in the estimation of an unknown unitary dynamics. [1] G. Chiribella, G. M. D'Ariano, and P. Perinotti, Theoretical framework for quantum networks, Phys. Rev. A 80, 022339 (2009)

Back to top of page

November 4, 2010

Speaker: Michael Ward University of British Columbia Dept. of Mathematics b>

Place: Arts 134

Time: 4:00

Title: An Asympotic Analysis of Localized Solutions to Some Diffusive and Reaction-Diffusion Systems

Abstract

A survey of the development of a unified singular perturbation methodology to analyze some linear and nonlinear PDE models of diffusion and reaction-diffusion type with localized solutions is presented. Specific results from this theory are given for three diverse applications. The first problem is to determine the mean first passage time (MFPT) for free diffusion from within a sphere to small localized traps on its boundary. In the context of cellular signal transduction, the results predict the time-scale needed for a diffusing molecule to arrive at localized signalling compartments on the boundary of a biological cell. From a mathematical viewpoint, the problem of optimizing this MFPT is shown to be closely related to the well-known Fekete point problem of finding the minimum energy configuration of repelling point charges on the surface of a sphere. Secondly, in the context of spatial ecology, a long-standing problem is to determine the persistence threshold for extinction of a species in a heterogeneous spatial landscape consisting of either favorable or unfavorable local habitats. For a 2-D spatial landscape consisting of such localized patches, and in the context of the diffusive logistic model, this extinction threshold is calculated asymptotically and the effects of both habitat fragmentation and habitat location on the persistence threshold are obtained. From a mathematical viewpoint, the persistence threshold represents the principal eigenvalue of an indefinite weight singularly perturbed eigenvalue problem. Finally, the dynamics, stability, and self-replication behavior of localized spot-type solutions to the well-known Gray-Scott reaction-diffusion model of chemical physics in a two-dimensional domain are discussed. Reduced ODE systems for the dynamics of spots are given together with phase diagrams in parameter space classifying the different types of spot instabilities. In this lecture I will emphasize the common mathematical features in the analysis of these three problems, most notably the role of the Neumann Green's function for the Laplacian. Applications of the results to more theoretical questions in PDE and spectral theory will also be emphasized.


Back to top of page

November 10, 2010

Walid Abou Salem, Department of Mathematics and Statistics

Place: McLean Hall 242.1

Time: 3:00-4:00 pm

Title: Renormalization group approach to singular perturbation theory for nonlinear PDEs.

Abstract

I discuss the rigorous application of the renormalization group method to (singular) perturbation theory for nonlinear partial differential equations. As a paradigm, I consider the concrete example of the nonlinear Schrodinger equation with quadratic nonlinearity in three spatial dimensions. I show how to obtain an approximate solution using the RG method together with an estimate of the difference between the true and approximate solutions. The analysis applies to differential equations where (space-time) resonances are present.


Back to top of page

November 17, 2010

Samuel Butler, Department of Geological Sciences

Place: McLean Hall 242.1

Time: 3:00-4:00 pm

Title: Transport, freezing and melting in mushy layers

Abstract

Mushy layers consist of a porous medium in which the solid matrix and interstitial liquid are close to thermodynamic equilibrium. Mass can be transferred between the solid and liquid by melting and freezing. Mushy layers are found in nature in magma chambers and sea ice and have been postulated to occur at the Earth's inner-core outer core, and outer-core mantle boundaries. They are also found in industrial settings including metal castings. In this talk, I will compare transport processes in mushy layers with those in non-reactive porous media and give expressions for effective transport velocities and diffusion coefficients and look at the degree of phase change that occurs for given system parameters. Some numerical simulations for a mushy system in which natural convection occurs will also shown and discussed.


Back to top of page

November 24, 2010

Saeed Torabi Ziaratgahi, Department of Mathematics and Statistics

Place: McLean Hall 242.1

Time: 3:00

Title: Domain decomposition for solving PDEs using RBF collocation methods

Abstract

The well-known finite difference, finite element, and finite volume methods for solving partial differential equations are based on a mesh discretization that may be a complicated and time-consuming process, particularly for complex, higher-dimensional geometries. The meshfree or meshless methods try to circumvent the cumbersome issues of mesh generation. One of the most common basic meshfree approximation methods is the Radial Basis Function (RBF) collocation method. Originally, the motivation for this method came from applications in geodesy, geophysics, mapping, and meteorology. Later, applications were found in areas such as the numerical solution of ODEs and PDEs, artificial intelligence, learning theory, neural networks, signal processing, sampling theory, statistics, finance, and optimization. RBF collocation methods are simple to implement because the collocation points need not have any connectivity requirement. They scale well with spatial dimension making them attractive for modeling high-dimensional problems. They also possess a high rate of convergence. For small to moderate-sized problems, RBF collocation methods outperform traditional methods, but for large problems, the resultant coefficient matrix is highly ill-conditioned, hindering the applicability of the RBF collocation methods. One of the best remedies to ill-conditioning problem is Domain Decomposition (DD) method, which splits the original domain into smaller sub-domains and solves the sub-problems in parallel. In this talk, we first describe the interpolation by RBFs, which is used to construct RBF collocation methods. Then we describe the RBF collocation methods for solving ODEs and PDEs. After that we present different type of DD methods and combine them with RBF collocation methods. Finally we show the efficiency of the approach by numerical examples.



Back to top of page

February 2, 2011

Ray Spiteri, Department of Computer Science

Place: McLean Hall 242.1

Time: 2:30-4:00 pm

Title: Stiffness analysis of cardiac cell models

Abstract

The electrophysiology in a cardiac cell can be modelled as a system of ordinary differential equations. The efficient solution of these systems is important because they must be solved many times as sub-problems of tissue- or organ-level simulations of cardiac electrophysiology. The wide variety of existing cardiac cell models encompasses many different properties, including the complexity of the model and the degree of stiffness. Accordingly, no single numerical method can be expected to be the most efficient for every model. In this talk, I discuss the stiffness properties of a range of cardiac cell models and discuss the implications for their numerical solution. This analysis allows us to select or design numerical methods that are highly effective for a given model and hence outperform commonly used methods.


Back to top of page

February 4, 2011

Professor Walter Craig (McMaster University, Department of Mathematics and Statistics ) PIMS distinguished Applied Mathematics Colloquium

Place: Arts 217

Time: 4:00-5:00 pm

Title: Waves and wave interactions: a paradigm in ocean waves

Abstract

Wave phenomena occur on an enormous range of scales, from the sub-quantum mechanical to the astrophysical. This talk will discuss the problem of free surface waves in water, which is a classical problem in mathematical hydrodynamics. We will discuss some of the common scale-independent features of wave prop- agation and wave interaction in this context, including detailed descriptions of nonlinear wave collisions, and the beginnings of a rigorous kinetic theory for a regime of wave turbulence.


Back to top of page

February 15, 2011

Professor Joel Feldman (University of British Columbia, Department of Mathematics)

Place: Arts 217

Time: 3:30-4:30 pm

Title: Towards the Construction of Bosonic Many-body Models

Abstract

Bosonic many-body models are an important class of mathematical models that are used to study gases of bosonic particles at very low temperatures. I will introduce

(i) the physical systems,

(ii) one way Physicists formulate the models in question,

(iii) why the models are expected to exhibit some very interesting behaviour,

(iv) why the models can be expected to be very difficult to deal with mathematically rigorously, and

(v) the first steps in a programme to construct the models mathematically rigorously.

This is joint work with Tadeusz Balaban of Rutgers University and Horst Knoerrer and Eugene Trubowitz of the ETH-Zurich.


Back to top of page

March 3, 2011

Professor David Brydges, Canada Research Chair (University of British Columbia, Department of Mathematics)

Place: Arts 263

Time: 3:30-4:30 pm

Title: Branched polymers and Mayer expansions

Abstract

The Mayer expansion is a power series expansion that has a central place in statistical mechanics. It is also full of combinatorial miracles that relate it to graphs, forests and branched polymers. I will discuss the background, the results, and some open problems.


Back to top of page

March 14, 2011

PIMS Applied Mathematics Series

Prof. Andrea Bertozzi, University of California, Los Angeles, CA, USA

Place:Arts 217

Time: 3:30 pm

Title: Swarming by nature and by design

Abstract

The cohesive movement of a biological population is a commonly observed natural phenomenon. With the advent of platforms of unmanned vehicles, such phenomena have attracted a renewed interest from the engineering community. This talk will cover a survey of the speaker's research and related work in this area ranging from aggregation models in nonlinear partial differential equations to control algorithms and robotic testbed experiments. We conclude with a discussion of some interesting problems for the applied mathematics community.

About The Speaker

Andrea Bertozzi is a mathematician who is known for her interdisciplinary work with computer scientists, physicists, and engineers. Much of her work has, in one way or another, examined the behavior of thin liquid films on hard surfaces. In tandem with physicists and engineers, she has worked at the Argonne National Laboratory and at Duke University constructing mathematical models that explain this and other physical phenomena. Born in 1965, in Boston, Massachusetts, to William and Norma Bertozzi, Andrea was encouraged by both of her parents to study and attend university. Her father, a professor of physics at the Massachusetts Institute of Technology, encouraged her to pursue her interest in the sciences. In 1991, she married Bradley Koetje, a management consultant. Bertozzi knew from an early age that she was interested in mathematics. Even in the first grade, she was captivated by the rudimentary math that was being taught and pushed to learn more. By high school, she had begun to learn advanced math and was concentrating on theory and abstract concepts, which she found to be the most interesting part of mathematics. After graduating from high school in Lexington, Massachusetts, in 1983, Bertozzi enrolled in Princeton University to study mathematics. She also studied a considerable amount of physics, although she took no degree in that subject. She earned her B.A. in math in 1987 and remained at Princeton to complete an M.S. in 1988 and a Ph.D. in 1991. After completing her Ph.D., Bertozzi took a position as L. E. Dickson Instructor of Mathematics at the University of Chicago. At Chicago, Bertozzi first became interested in the mathematics of thin films. She began working with a group of physicists who were studying mathematical models that described the behavior of phenomena that were similar to thin films. Gradually, the problem centered specifically on a mathematical description of liquids flowing on a solid surface. This was an area of mathematics that had not received much attention but had been researched by physicists since the 1960s. Bertozzi remained at the University of Chicago until 1995 when she was offered the position of associate professor at Duke University in Durham, North Carolina. Then during 1995-96, she worked at the Argonne National Laboratory, located outside of Chicago in Argonne, Illinois. Here, as a MARIA GOEPPERT MAYER Distinguished Scholar, she continued her work in the field of scientific computing, which she had begun at the University of Chicago. The purpose of scientific computing is to create computer models that simulate physical processes on the computer. In this way, virtual experiments that can mimic actual physical conditions are created. At Argonne, Bertozzi continued her study of the mathematical-physical properties of thin liquids on dry surfaces. This problem, which seems relatively simple, is actually complicated. A liquid applied to a dry surface will not spread evenly but will pool and spread onto the surface in fingerlike rivulets. Bertozzi worked on a set of partial differential equations, also called evolution equations because this kind of math describes an event occurring over time, that fit a model for film-coating behavior into mathematical terms. This work, although basic research, may someday be helpful for industries such as the microchip-manufacturing sector, which needs to understand this coating process in making their complicated and delicate product. After her year at the Argonne Lab, Bertozzi returned to her job as associate professor of mathematics at Duke University in 1996. In 1998, she became associate professor of mathematics and physics, and in 1999, she became a full professor in both disciplines. Currently, she is director of Duke’s Center for Nonlinear and Complex Systems, an interdisciplinary research center that includes scientists from the disciplines of math, biology, engineering, medical sciences, and environmental studies. In addition to her studies of thin films on hard surfaces, Bertozzi works in more general problems of fluid dynamics. Bertozzi was recognized for her work by the Sloan Foundation, which awarded her a research fellowship in 1995. In 1996, she was presented the Presidential Early Career Award for Scientists and Engineers by the U.S. Office of Naval Research. Cambridge University Press published her book, coauthored with Andrew Majda, Vorticity and Incompressible Flow, in 2000.


Back to top of page

March 23, 2011

Dr. Feride Tiglay Fields Institute, Toronto, ON

Place:MCLH 242.1

Time:4:00 pm

Title: Integrable evolution equations on spaces of tensor densities and their peakon solutions

Abstract

In a pioneering paper V. Arnold presented a general framework within which it is possible to employ geometric and Lie theoretic techniques to study a variety of equations of interest in mathematical physics. I will describe how to extend his formalism using tensor densities and introduce two integrable PDE. One of the equations turns out to be closely related to the inviscid Burgers equation while the other has not been identified in any form before. These two PDE possess all the hallmarks of integrability: the existence of a Lax pair formulation, a bihamiltonian structure, the presence of an infinite family of conserved quantities and the ability to write down explicitly some of its solutions. I will also talk about local well-posedness of the corresponding Cauchy problem and some results on blow-up as well as global existence of solutions. Time permitting, I will describe the peakon solutions for these equations.


Back to top of page

March 25, 2011

Dr. Feride Tiglay Fields Institute, Toronto, ON

Place:MCLH 242.1

Time: 2:30 pm

Title: The Periodic Cauchy Problem for Novikov's Equation

Abstract

We study the periodic Cauchy problem for an integrable equation with cubic nonlinearities introduced by Novikov. Like the Camassa–Holm and Degasperis–Procesi equations, Novikov’s equation has Lax pair representations and admits peakon solutions, but it has nonlinear terms that are cubic, rather than quadratic. We show the local well-posedness of the problem in Sobolev spaces and existence and uniqueness of solutions for all time using orbit invariants. Furthermore, we prove a Cauchy–Kowalevski type theorem for this equation, which establishes the existence and uniqueness of real analytic solutions.


Back to top of page

April 20, 2011

Speaker: Marco Merkli (Memorial University)

Place: McLean 242.2

Time: 3:00 pm

Title: Quantum Scattering Measurement

Abstract

We consider a quantum system (scatterer) which interacts with a sequence of identical, independent quantum systems (scattering probes). The interaction is sequential, one by one. After leaving the scatterer, a quantum measurement is perfomed on each probe. The measurement outcomes form a random process. We analyze the asymptotic properties of this process, such as the probability of convergence. If the process converges, then the scatterer is driven to a final state determined by the measurement outcomes. We also examine large deviations for the average of the measurements. We illustrate the concepts and results on the (truncated) Jaynes-Cummings model, where both the scatterer and the probes are spins 1/2, representing two degrees of freedom active in the scattering process.


Back to top of page

May 20, 2011

Speaker: Calin Atanasiu (EURATOM MEdC Association, Bucharest and Max-Planck Institut für Plasmaphysik, Garching bei Munich, Germany)

Place: McLean 242.1

Time: 2:00 pm

Title: Special aspects of MHD calculations in tokamaks

Abstract

Magneto-hydrodynamic equilibrium and stability calculations (analytical and numerical) for real diverted tokamak configurations are presented. (a) Two families of exact analytical solutions of the Grad-Shafranov equation are presented by specifying the highest polynomial dependence of the plasma current density on the solution in such a way that the Grad-Shafranov equation becomes a linear inhomogeneous differential equation. (b) By introducing a "cast function" in a classical flux coordinate system, in the presence of a separatrix, the solution of the equilibrium equation - the unknown moments - is determined by the difference between the real flux surface contours and those described by the cast functions only. Thus, the necessary number of moments is small enough to make computations time-efficient. (c) For instability calculation of tearing and external kink type, the expression of the potential energy has been written in terms of the perturbation of the flux function, and performing an Euler minimization, a system of ordinary differential equations in that perturbation has been obtained. For a diverted configuration, the usual vanishing boundary conditions for the perturbed flux function at the magnetic axis and at infinity can no longer be used. An approach to fix "natural" boundary conditions for the perturbed flux function just at the plasma boundary has been developed; this replaces the vanishing boundary conditions at infinity. Special attention is given to the stabilization of external kink modes in the presence of a conducting wall - the resistive wall modes, the most dangerous instability for the future tokamak reactors.


Back to top of page

June 30, 2011

Speaker: Prof. Jean-Francois Ganghoffer, LEMTA, Nancy University, France

Place: McLean 242.1

Time: 11:00 am

Title: Mechanics and thermodynamics of surface growth. Application to bone remodeling

Abstract

The surface growth of biological tissues is presently analyzed at the continuum scale of tissue elements, adopting the framework of the thermodynamics of surfaces and in line with Eshelbian mechanics. From a kinematic viewpoint, growth is assumed to occur in a moving referential configuration, considered as an open evolving domain exchanging mass, work, and energy with its environment. The growing surface is endowed with a superficial excess concentration of moles, which is ruled by an appropriate kinetic equation. The material surface forces for growth are evaluated versus a surface Eshelby stress, the curvature tensor of the growing surface, the gradient of the chemical energy of nutrients and the applied superficial force field. A system of coupled field equations is written for the superficial density of minerals, their concentration and the surface velocity of the growing surface. Application of the developed formalism to bone external remodeling highlights the interplay between transport phenomena and generation of surface mechanical forces. The model is able to describe both bone growth and resorption, according to the respective magnitude of the chemical and mechanical contributions to the material surface driving force for growth.


Back to top of page

July 7, 2011

Speaker: Ahmed Kaffel, Postdoctoral Fellow, Department of Computer Science, University of Saskatchewan

Place: McLean 242.1

Time: 2:00 pm

Title: On the Stability of plane viscoelastic shear flows in the limit of infinite Weissenberg and Reynolds numbers

Abstract

Elastic effects on the hydrodynamic instability of inviscid parallel shear flows are investigated through a linear stability analysis. We focus on the upper convected Maxwell model in the limit of infinite Weissenberg and Reynolds numbers. Specifically, we study the effects of elasticity on the instability of a few classes of simple parallel flows, specifically plane Poiseuille and Couette flows, the hyperbolic-tangent shear layer and the Bickley jet. The equation for stability is derived and solved numerically using the Chebyshev collocation spectral method. This algorithm is computationally efficient and accurate in reproducing the eigenvalues. We consider flows bounded by walls as well as flows bounded by free surfaces. In the inviscid, nonelastic case all the flows we study are unstable for free surfaces. In the case of wall bounded flow, there are instabilities in the shear layer and Bickley jet flows. In all cases, the effect of elasticity is to reduce and ultimately suppress the inviscid instability. The numerical solutions are compared with the analysis of the long wave limit and excellent agreement is shown between the analytical and the numerical solutions. We found flows which are long wave stable, but nevertheless unstable to wave numbers in a certain finite range. While elasticity is ultimately stabilizing, this effect is not monotone; there are instances where a small amount of elasticity actually destabilizes the flow. The linear stability in the short wave limit of shear flows bounded by two parallel free surfaces is investigated. Unlike the plane Couette flow which has no short wave instability, we show that plane Poiseuille flow has two unstable eigenmodes localized near the free surfaces which can be combined into an even and an odd eigenfunctions. The derivation of the asymptotics of these modes shows that our numerical eigenvalues are in agreement with the analytic formula and that the difference between the two eigenvalues tends to zero exponentially with the wavenumber α.


Back to top of page

August 11, 2011

Speaker: Professor Abba Gumel, Department of Mathematics, University of Manitoba.

Place: McLean 242.1

Time: 2:00 pm

Title: Dynamically consistent finite-difference methods for differential equations

Abstract

Standard numerical integrators, such as the Runge-Kutta family of explicit finite-difference methods, are known to exhibit numerous scheme-dependent instabilities and generally fail to preserve some of the main essential qualitative features (such as positivity, boundedness, asymptotic stability, and bifurcation properties) of the governing continuous system they approximate. The talk will address the problem of designing appropriate discrete-time models that are dynamically consistent with the corresponding continuous-time model they approximate. Some models arising from modeling real-life phenomena in the natural and engineering sciences will be discussed.

Biodata: Abba Gumel is a Professor of Mathematics at the University of Manitoba. His research work is based on the design and analysis (qualitative and quantitative) of models for the spread and control of emerging and re-emerging diseases of public health significance. His homepage is http://home.cc.umanitoba.ca/~gumelab.


Back to top of page