MCTP Home
Research Experience
     for Undergraduates
Return to:
Department Home
UNM
|
     
     
     
PROPOSED REU PROJECTS
The following are sample REU projects proposed by the faculty of the
Department of Mathematics and Statistics at UNM. They serve as guidelines for
choosing a project.
Interpolation using polynomial of 2 variables,
mentored by Michael Nakamaye.
Lagrange interpolation gives the equation of a function of
the form y = f(x), where f(x) is a
polynomial, passing through specified points in the
Euclidean plane. If one looks instead for a more general polynomial of
the form p(x,y) = 0, and in addition looks for higher multiplicities at
the specified points, much less is known. This problem provides a nice
mix of algebra and geometry and plent of room for creative effort!
Self--accelerating flows for the Navier--Stokes equations,
mentored by Jens Lorenz.
Prof. Lorenz and Dr. Randy Ott recently succeeded in
constructing an incompressible velocity field $u(x)$ which accelerates
under its own pressure. (The phenomenon of self--acceleration is necessary,
but not sufficient, for blow--up to occur.
The blow--up problem for Navier--Stokes is one of the Clay prize problems.)
One can introduce parameters in the self---accelerating flow and
ask: How large can one make the acceleration when certain norms of
$u(x)$ (e.g., the energy) remain fixed?
A good visualization of the self--accelerating flow is also of interest.
Simulation of a 1D polymer system, mentored by Daniel Appelo.
The aim of this project is to simulate a basic polymer problem such as
the self-assembly of a spherical micelle or the phase transition in a
oil-water mixture. Such systems are be described by the so called
self-consistent-field equations which can be solved numerically by
a combination of solution of non-linear equations, solution of
linear parabolic partial differential equations and computation of
integrals using quadrature. As an extension the detailed description
of the pathway between states in the system can be modeled by the
so called string method. This requires the solution of ordinary
differential equations and interpolation between system states.
Discrete and dyadic harmonic analysis, mentored by Cristina Pereyra.
Advanced undergraduate students who have
been exposed to basic ideas of Fourier analysis and are acquainted
with the Haar basis, can start exploring
a number of important concepts in harmonic analysis in the discrete
(purely linear algebra) and in the dyadic setting,
where the theory is much cleaner and transparent,
see Pereyra & Ward (to appear). One can analyze
the finite dimensional cases and study, analytically or numerically,
what happens as the dimension
goes to infinity, and justify the appearance of certain
functional analysis tools.
One can also explore and justify very powerful algorithms, based
on subtle time/frequency
arguments (for example the non-standard form
introduced by Beylkin, Coifman and Rokhlin (1991), which requires
good understanding of the multiresolution analysis, and almost orthogonality
arguments), in fact students can be lead
to discover by themselves these algorithms in simple case studies.
G. Beylkin, R. Coifman, V. Rokhlin (1991).
{Fast wavelet transforms and numerical
algorithms.I} Comm. Pure Appl. Math., 44(2), p.141--183.
M.C. Pereyra, L. Ward (2012),
{\em Harmonic Analysis: from Fourier to Wavelets}.
Undergraduate STML series of the AMS (420+pp).
3d vortex filament motion, mentored by M. Nitsche.
Analisa Calini used local induction approximation to obtain knotted steady
solutions of vortex filaments.
We will compute 3d simulations of regularized filament evolution and
compare the results to analytical solutions.
Evolution of elliptical vortex sheets, mentored by M. Nitsche.
An initially flat circular vortex sheet inducing potential flow
past a disk evolves under its self-induced motion, and rolls
up to form one vortex ring (Krasny & Nitsche 2002).
An initially spherical vortex sheet inducing potential flow
past a sphere also rolls up into a vortex ring. However, unlike the flat
sheet case, this ring sheds an infinite set of self-similar vortices
(Nitsche 2001).
We will investigate the evolution of a family of initially elliptical
sheets spanning the range between the disk and the sphere and
determine how the vortex ring shedding depends on the initial ellipticity.
R. Krasny, M. Nitsche (2002). {The onset of chaos in vortex
sheet flow}, J. Fluid Mech.\ {\bf 454}, pp 47--69 (2002).
M. Nitsche (2001). {Self-similar shedding of vortex rings},
J. Fluid Mech.\ {\bf 435}, pp 397--407.
Discrete-time modeling of mosquito-borne
disease and control, mentored by Helen Wearing.
Extend the simple discrete-time model of infectious disease to include human and mosquito populations and sources of stochasticity, using the mosquito-borne disease, dengue, as a case study. Conduct sensitivity analysis of the model and investigate the impacts of mosquito control. Potentially estimate parameters by fitting the model to epidemic data.
Lift and drag of airfoils, mentored by Peter Vorobieff.
The undergraduate research project will involve conducting experiments
with a wind and water tunnel, analyzing their results to determine the
dependence of lift and drag on the flow speed and angle of attack,
where applicable, for several bodies: a symmetric airfoil, an
asymmetric airfoil, a cylinder, and a sphere. In addition, the
students will be required to produce flow-visualization images with
the digital camera and to comment on the transient patterns.
Calculus in the ring of integers, mentored by Alexandru Buium and Santiago Simanca.
There is an arithmetic analogue of differential calculus
in which functions of a real variable are replaced by
integer numbers
and in which the derivative operator on functions is replaced by a
Fermat quotient operatorĀ $\delta: Z \rightarrow Z$,
$\delta x=(x-x^p)/p$. Cf. A. Buium, Inventiones Math., 122, 2 (1995),
309-340. This can be generalised
to integers in numberĀ fields
with bounded ramification.
Accordingly there is a notion of (non-linear)
differential operator acting on such integers.
One possible project would be to extend this theory to
integers in fields with unbounded ramification. This might
have applications to diophantine problems (such as
Mazur's question on effective bounds for the torsion
of curves in Jacobians; such applications are known in
the unramified case, cf. A.Buium, Duke Math. J. 82, 2 (1996), 349-367.)
Another possible project would be to express various well-known
arithmetically flavored functions (such as power residue symbols,
Hecke characters, e.t.c.)
as (non-linear) ``differential operators'' in this arithmetic theory.
This would be part of a more general conjectural picture explained in:
A. Buium,
Trans. AMS 357, 3 (2005), 901-964. Finally one could investigate
arithmetic analogues of various classical ordinary or partial
differential equations
and/or techniques related to such equations (along the lines of
the papers by A. Buium and S. Simanca,
Arithmetic partial differential equations I and II, Advances in Math, 225 (2010).
Investigation of fluid cavitation, mentored by Deborah Sulsky.
Cavitation was first studied by Lord Rayleigh in the late 19th century.
The collapse of cavities is a relatively low energy event, but it is highly
localized and therefore can erode metals and shorten the useful life of
pumps or propellors. However, cavitation can also be harnessed for
productive purposes such as ultrasonic cleaning, or the breakdown of
suspended particles in paints or milk. Cavitation also plays a role in the
destruction of kidney stones in shock wave lithotripsy.
Despite the wide range of scenarios in which cavitation can occur, both
constructively and destructively, little is understood about the onset of
cavitation.
We propose to use datasets from high resolution numerical simulations
(already computed and available). A student will be responsible for taking
the computed stresses from numerical simulations and transforming them to
principal stresses and evaluating various cavitation criteria proposed in
the literature (eg. Joseph, 1998) and extensions of these criteria
suggested by our own research in failure modeling (eg. Schreyer, et al.,
2006).
D. D. Joseph (1998). Cavitation and the state of stress in a flowing
fluid, J. Fluid Mech., 366, pp 367-378.
H. L. Schreyer, D. L. Sulsky, L. B. Munday, M. D. Coon, and R. Kwok (2006).
Elastic-decohesive constitutive model for sea ice, J. Geophys. Res.,
C11S26, doi:10.1029/205JC003334.
Spectral shift functions, mentored by A. Skripka.
In problems of quantum mechanics (for instance, in computation of the energy of a
system or time evolution of a system), one encounters linear operators [that is,
infinite-dimensional generalizations of matrices] that arise as functions of other
operators. It can be important to know how a small perturbation of an operator effects
the change of a function of the operator. Some aspects of this change are reflected in
spectral shift functions, which originate from I.M. Lifshits's work on theoretical
physics of 1947 - 1952. Existence of higher order spectral shift functions, which
allow to encompass more general perturbations, has recently been proved by D. Potapov,
A. Skripka, and F. Sukochev, but properties of these functions are yet to be
understood. Testing/discovering (and, later, proving) properties of these functions
should start with finite matrices (of low dimensions). To be ready to start working on
a project, a student should possess knowledge of basic linear algebra. Familiarity
with mathematics software (Matlab, Maple, or Mathematica) would also be useful.
|