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)
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,
Place: McLean Hall 242.1
Time: 3:00
Title: Mathematical Analysis of dynamics of Chlamydia trachomatis
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.
Place: McLean Hall 242.1
Time: 3:00
Title: Propagating Uncertainties in Modeling Nonlinear Dynamic Systems
Place: Arts 134
Time: 4:00
Title: Optimal Estimation of Quantum Signals in the Presence of Symmetry
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.
Place: McLean Hall 242.1
Time: 3:30
Title:Quantum measurement theory: from quantum states to quantum networks
Place: Arts 134
Time: 4:00
Title: An Asympotic Analysis of Localized Solutions to Some Diffusive and Reaction-Diffusion Systems
Place: McLean Hall 242.1
Time: 3:00-4:00 pm
Title: Renormalization group approach to singular perturbation theory for nonlinear PDEs.
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.
Place: McLean Hall 242.1
Time: 3:00-4:00 pm
Title: Transport, freezing and melting in mushy layers
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.
Place: McLean Hall 242.1
Time: 3:00
Title: Domain decomposition for solving PDEs using RBF collocation methods
Place: McLean Hall 242.1
Time: 2:30-4:00 pm
Title: Stiffness analysis of cardiac cell models
Place: Arts 217
Time: 4:00-5:00 pm
Title: Waves and wave interactions: a paradigm in ocean waves
Place: Arts 217
Time: 3:30-4:30 pm
Title: Towards the Construction of Bosonic Many-body Models
(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.
Place: Arts 263
Time: 3:30-4:30 pm
Title: Branched polymers and Mayer expansions
Prof. Andrea Bertozzi,
University of California, Los Angeles, CA, USA
Place:Arts 217
Time: 3:30 pm
Title: Swarming by nature and by design
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.
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
Dr. Feride Tiglay
Fields Institute, Toronto, ON
Place:MCLH 242.1
Time: 2:30 pm
Title: The Periodic Cauchy Problem for Novikov's Equation
Speaker: Marco Merkli (Memorial University)
Place: McLean 242.2
Time: 3:00 pm
Title: Quantum Scattering Measurement
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
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
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
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
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.
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).
APPLIED MATHEMATICS/ MATHEMATICAL PHYSICS SEMINAR 2009-2010
(sponsored by MITACS)
September 9, 2010
Speaker: Dr. Oluwaseun Sharomi
Department of Mathematics
University of Manitoba
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.
Back to top of page
September 23, 2010
Speaker: Professor George Corliss
Department of Electrical and Computer Engineering
Marquette University
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>
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.
Back to top of page
November 3, 2010
Speaker: Dr Giulio Chiribella
Perimeter Institute for Theoretical Physics, Waterloo, Ontario
b>
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)
November 4, 2010
Speaker:
Michael Ward
University of British Columbia
Dept. of Mathematics b>
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
Abstract
Back to top of page
November 17, 2010
Samuel Butler, Department of Geological Sciences
Abstract
Back to top of page
November 24, 2010
Saeed Torabi Ziaratgahi, Department of Mathematics and Statistics
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
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
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)
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
Back to top of page
March 3, 2011
Professor David Brydges, Canada Research Chair (University of British Columbia,
Department of Mathematics)
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
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.
Back to top of page
March 23, 2011
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
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
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
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
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
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
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.
Back to top of page