Pete Casazza : Applications of Hilbert space frames
- Applied Math and Analysis ( 116 Views )Hilbert space frames have traditionally been used in signal/image processing. Recently, there have arisen a variety of new applications to speeding up the internet, producing cell phones which won't fade, quantum information theory, distributed processing and more. We will review the fundamentals of frame theory and then look at the myriad of applications of frames.
Tom Beale : Uniform error estimates for fluid flow with moving boundaries using finite difference methods
- Applied Math and Analysis ( 98 Views )Recently there has been extensive development of numerical methods for fluid flow interacting with moving boundaries or interfaces, using regular finite difference grids which do not conform to the boundaries. Simulations at low Reynolds number have demonstrated that, with certain choices in the design of the method, the velocity can be accurate to about O(h^2) while discretizing near the interface with truncation error as large as O(h). We will describe error estimates which verify that such accuracy can be achieved in a simple prototype problem, even near the interface, using corrections to difference operators as in the immersed interface method. We neglect errors in the interface location and derive uniform estimates for the fluid velocity and pressure. We will first discuss maximum norm estimates for finite difference versions of the Poisson equation and diffusion equation with a gain of regularity. We will then describe the application to the Navier-Stokes equations.
Paul Tupper : The Relation Between Shadowing and Approximation in Distribution
- Applied Math and Analysis ( 151 Views )In computational physics, molecular dynamics refers to the computer simulation of a material at the atomic level. I will consider classical deterministic molecular dynamics in which large Hamiltonian systems of ordinary differential equations are used, though many of the same issues arise with other models. Given its scientific importance there is very little rigorous justification of molecular dynamics. From the viewpoint of numerical analysis it is surprising that it works at all. The problem is that individual trajectories computed by molecular dynamics are accurate for only small time intervals, whereas researchers trust the results over very long time intervals. It has been conjectured that molecular dynamics trajectories are accurate over long time intervals in some weak statistical sense. Another conjecture is that numerical trajectories satisfy the shadowing property: that they are close over long time intervals to exact trajectories with different initial conditions. I will explain how these two views are actually equivalent to each other, after we suitably modify the concept of shadowing.
Casey Rodriguez : The Radiative Uniqueness Conjecture for Bubbling Wave Maps
- Applied Math and Analysis ( 179 Views )One of the most fundamental questions in partial differential equations is that of regularity and the possible breakdown of solutions. We will discuss this question for solutions to a canonical example of a geometric wave equation; energy critical wave maps. Break-through works of Krieger-Schlag-Tataru, Rodnianski-Sterbenz and Rapha ̈el-Rodnianski produced examples of wave maps that develop singularities in finite time. These solutions break down by concentrating energy at a point in space (via bubbling a harmonic map) but have a regular limit, away from the singular point, as time approaches the final time of existence. The regular limit is referred to as the radiation. This mechanism of breakdown occurs in many other PDE including energy critical wave equations, Schro ̈dinger maps and Yang-Mills equations. A basic question is the following: • Can we give a precise description of all bubbling singularities for wave maps with the goal of finding the natural unique continuation of such solutions past the singularity? In this talk, we will discuss recent work (joint with J. Jendrej and A. Lawrie) which is the first to directly and explicitly connect the radiative component to the bubbling dynamics by constructing and classifying bubbling solutions with a simple form of prescribed radiation. Our results serve as an important first step in formulating and proving the following Radiative Uniqueness Conjecture for a large class of wave maps: every bubbling solution is uniquely characterized by it’s radiation, and thus, every bubbling solution can be uniquely continued past blow-up time while conserving energy.
Blair Sullivan : Can we Quantify & Exploit Tree-like Intermediate Structure in Complex Networks?
- Applied Math and Analysis ( 108 Views )Large complex networks naturally represent relationships in a variety of settings, e.g. social interactions, computer/communication networks, and genomic sequences. A significant challenge in analyzing these networks has been understanding the intermediate structure those properties not captured by metrics which are local (e.g. clustering coefficient) or global (e.g. degree distribution). It is often this structure which governs the dynamic evolution of the network and behavior of diffusion-like processes on it. Although there is a large body of empirical evidence suggesting that complex networks are often tree-like at intermediate to large size-scales (e.g. work of Boguna et al in physics, Kleinberg on internet routing, and Chung & Lu on power-law graphs), it remains a challenge to take algorithmic advantage of this structure in data analysis. We discuss several approaches and heuristics for quantifying and elucidating tree-like structure in networks, including various tree-decompositions and Gromov's delta hyperbolicity. These approaches were developed with very different "tree-like" applications in mind, and thus we discuss the strengths and short-comings of each in the context of complex networks and how each might aid in identifying intermediate-scale structure in these graphs.
Jean-Philippe Thiran : Multimodal signal analysis for audio-visual speech recognition
- Applied Math and Analysis ( 111 Views )After a short introduction presenting our group and our main research topics, I will address the problem of audio-visual speech recognition, i.e. a typical example of multimodal signal analysis, when we want to extract and exploit information coming from two different but complementary signals: an audio and a video channel. We will discuss two important aspects of this analysis. We will first present a new feature extraction algorithm based in information theoretical principles, and show its performances, compared to other classical approaches, in our multimodal context. Then we will discuss multimodal information fusion, i.e. how to combine information from those two channels for optimal classification.
Peter Diao : Model-Free Consistency of Graph Partitioning using Dense Graph Limits
- Applied Math and Analysis ( 141 Views )The beautiful work of Borgs, Chayes, Lovasz, Sos, Szegedy, Vesztergombi, and many others on dense graph limits has received quite a bit of attention in pure math as well as statistics and machine learning. In this talk we will review some of the previous work on dense graph limits and then present recent work on providing a more robust mathematical framework for proving the statistical consistency of graph partitioning algorithms such as spectral clustering. A striking feature of our approach is that it is model-free, compared to the popular iid paradigm. Our results are thus broadly applicable in real-world settings, where it is notoriously difficult to obtain relevant models for network data, and observations are not independent. At the end, I will discuss implications for how mathematical foundations can be developed for other modern data analysis techniques. This is joint work with Dominique Guillot, Apoorva Khare, and Bala Rajaratnam. Preprint available at https://arxiv.org/abs/1608.03860.
Sung Ha Kang : Efficient methods for curvature based variational imaging models
- Applied Math and Analysis ( 97 Views )Starting with an introduction to multiphase image segmentation, this talk will focus on inpainting and illusory contour using variational models with curvature terms. Recent developments of fast algorithms, based on operator splitting, augmented Lagrangian, and alternating minimization, enabled us to efficiently solve functional with higher order terms. Main ideas of the models and algorithms, some analysis and numerical results will be presented.
Jonathan Hermon : Mixing and hitting times - theory and applications
- Applied Math and Analysis ( 108 Views )We present a collection of results, based on a novel operator maximal inequality approach, providing precise relations between the time it takes a Markov chain to converge to equilibrium and the time required for it to exit from small sets. These refine results of Aldous and Lovasz & Winkler. Among the applications are: (1) A general characterization of an abrupt convergence to equilibrium phenomenon known as cutoff. Specializing this to Ramanujan graphs and trees. (2) Proving that the return probability decay is not geometrically robust (resolving a problem of Aldous, Diaconis - Saloff-Coste and Kozma). (3) Random walk in evolving environment
John Neu : Resonances in Geometric Optics
- Applied Math and Analysis ( 98 Views )Consider wavefront propagation in the plane, in a medium whose propagation speed is doubly periodic. Think of a "wavefront" as the moving boundary between "light" and "darkness." There are "macroscopic plane wavefronts", for which the wavefront is everywhere and always bounded close to a moving line with constant normal velocity. The normal velocity depends on direction. Some nice function of angle, and we ought to compute it, no? To contemplate this calculation, visualize an infinite planar vineyard, with a square lattice of grape stakes. Gaze outward from the the grape stake at the origin. In some directions, the line of sight is blocked by a stake, and we'll call these directions "rational." The rational directions are dense, but of zero measure on the unit circle. The simplest formal asymptotics of the direction dependent speed produces a strange result: A formal series, one term for each rational direction. The graph of an individual term as a function of angle looks like the amplitude response of a forced oscillator, with a plus infinity vertical asymptote as you scan through the rational direction. Nevertheless, the series converges absolutely for almost all directions on the unit circle. (In the vineyard, almost all lines of sight escape to infinity.) At the outset, the prognosis seems: "Difficult and obscure." The problem of direction dependent speed does not LOOK like a venue for formal classical asymptotics, but that's what this talk proposes. The ingredients are standard matched asymptotic expansions, and baby number theory concerning period cells of the vineyard.
Edward Waymire : Dispersion in the Presence of Interfacial Discontinuities
- Applied Math and Analysis ( 97 Views )This talk will focus on probability questions arising in the geophysical and biological sciences concerning dispersion in highly heterogeneous environments, as characterized by abrupt changes (discontinuities) in the diffusion coefficient. Some specific phenomena observed in laboratory and field experiments involving breakthrough curves (first passage times), occupation times, and local times will be addressed. This is based on joint work with Thilanka Appuhamillage, Vrushali Bokil, Enrique Thomann, and Brian Wood at Oregon State University.
Hien Tran : HIV Model Analysis under Optimal Control Based Treatment Strategies
- Applied Math and Analysis ( 142 Views )In this talk, we will introduce a dynamic mathematical model that describes the interaction of the immune system with the human immunodeficiency virus (HIV). Using optimal control theory, we will illustrate that optimal dynamic multidrug therapies can produce a drug dosing strategy that exhibits structured treatment interruption, a regimen in which patients are cycled on and off therapy. In addition, sensitivity analysis of the model including both classical sensitivity functions and generalized sensitivity functions will be presented. Finally, we will describe how stochastic estimation can be used to filter and estimate states and parameters from noisy data. In the course of this analysis it will be shown that automatic differentiation can be a powerful tool for this type of study.
Yifei Lou : Nonconvex Approaches in Data Science
- Applied Math and Analysis ( 95 Views )Although big data is ubiquitous in data science, one often faces challenges of small data, as the amount of data that can be taken or transmitted is limited by technical or economic constraints. To retrieve useful information from the insufficient amount of data, additional assumptions on the signal of interest are required, e.g. sparsity (having only a few non-zero elements). Conventional methods favor incoherent systems, in which any two measurements are as little correlated as possible. In reality, however, many problems are coherent. I will present a nonconvex approach that works particularly well in the coherent regime. I will also address computational aspects in the nonconvex optimization. Various numerical experiments have demonstrated advantages of the proposed method over the state-of-the-art. Applications, ranging from super-resolution to low-rank approximation, will be discussed.
Mary Lou Zeeman : Modeling the Menstrual Cycle:How does estradiol initiate the LH surge?
- Applied Math and Analysis ( 142 Views )In vertebrates, ovulation is triggered by a surge of luteinizing hormone (LH) from the pituitary. The precise mechanism by which rising estradiol (E2) from the ovaries initiates the LH surge in the human menstrual cycle remains a mystery. The mystery is due in part to the bimodal nature of estradiol feedback action on LH secretion, and in part to disagreement over the site of the feedback action.
We will describe a differential equations model in which the mysterious bimodality of estradiol action arises from the electrical connectivity of a network of folliculo-stellate cells in the pituitary. The mathematical model is based as closely as possible on current experimental data, and is being used to design and conduct new experiments. No biological background will be assumed.
Alina Chertock : Numerical Methods for Chemotaxis and Related Models
- Applied Math and Analysis ( 83 Views )Chemotaxis is a movement of micro-organisms or cells towards the areas of high concentration of a certain chemical, which attracts the cells and may be either produced or consumed by them. In its simplest form, the chemotaxis model is described by a system of nonlinear PDEs: a convection-diffusion equation for the cell density coupled with a reaction-diffusion equation for the chemoattractant concentration. It is well-known that solutions of such systems may develop spiky structures or even blow up in finite time provided the total number of cells exceeds a certain threshold. This makes development of numerical methods for chemotaxissystems extremely delicate and challenging task. In this talk, I will present a new family of high-order finite-volume finite-difference methods for the Keller-Segel chemotaxis system and several related models. Applications of the proposed methods to the classical Patlak-Keller-Segel model, its extensions to the two-species case as well as to the coupled chemotaxisfluid system will also be discussed.
Christoph Ortner : Multi-scale simulation of crystal defects
- Applied Math and Analysis ( 110 Views )PART 1: I will construct a mathematical model of a defect embedded in an infinite homogeneous crystal. I will then establish a regularity result for minimisers, which given the crucial information on which approximation schemes are based. As an elementary application of this framework I will prove convergence rates for two computational schemes: (1) clamped far-field and (2) coupling to harmonic far-field model.
PART 2: The conditions under which the theory of Part 1 holds are separability and locality of the total energy. In Part 2 I will show how for a tight-binding model (a minimalistic electronic structure model) these two condition arise. This analysis raises some interesting (open) questions.
PART 3: Finally, I will use the theory developed in PART 1 and PART 2 to construct and analyse a new family of QM/MM embedding schemes with rigorous error estimates.
Zhennan Zhou : On the classical limit of a time-dependent self-consistent field system: analysis and computation
- Applied Math and Analysis ( 101 Views )We consider a coupled system of Schroedinger equations, arising in quantum mechanics via the so-called time-dependent self-consistent field method. Using Wigner transformation techniques we study the corresponding classical limit dynamics in two cases. In the first case, the classical limit is only taken in one of the two equations, leading to a mixed quantum-classical model which is closely connected to the well-known Ehrenfest method in molecular dynamics. In the second case, the classical limit of the full system is rigorously established, resulting in a system of coupled Vlasov-type equations. In the second part of our work, we provide a numerical study of the coupled semiclassically scaled Schroedinger equations and of the mixed quantum-classical model obtained via Ehrenfest's method. A second order (in time) method is introduced for each case. We show that the proposed methods allow time steps independent of the semi-classical parameter(s) while still capturing the correct behavior of physical observables. It also becomes clear that the order of accuracy of our methods can be improved in a straightforward way.
Yifeng Yu : Random Homogenization of Non-Convex Hamilton-Jacobi Equations in 1d
- Applied Math and Analysis ( 99 Views )I will present the proof of the random homogenization of general coercive Hamiltonian in 1d with the form as H(p,x,\omega)=H(p)+V(x, \omega). Some interesting and complex phenomena associated with non-convex Hamiltonian will also be discussed. This is a joint work with Scott Armstrong and Hung Tran.
Johann Guilleminot : Stochastic Modeling and Simulations of Random Fields in Computational Nonlinear Mechanics
- Applied Math and Analysis ( 87 Views )Accounting for system-parameter and model uncertainties in computational models is a highly topical issue at the interface of computational mechanics, materials science and probability theory. In addition to the construction of efficient (e.g. Galerkin-type) stochastic solvers, the construction, calibration and validation of probabilistic representations are now widely recognized as key ingredients for performing accurate and robust simulations. This talk is specifically focused on the modeling and simulation of spatially-dependent properties in both linear and nonlinear frameworks. Information-theoretic models for matrix-valued random fields are first introduced. These representations are typically used, in solid mechanics, to define tensor-valued coefficients in elliptic stochastic partial differential operators. The main concepts and tools are illustrated, throughout this part, by considering the modeling of elasticity tensors fluctuating over nonpolyhedral geometries, as well as the modeling and identification of random interfaces in polymer nanocomposites. The latter application relies, in particular, on a statistical inverse problem coupling large-scale Molecular Dynamics simulations and a homogenization procedure. We then address the probabilistic modeling of strain energy functions in nonlinear elasticity. Here, constraints related to the polyconvexity of the potential are notably taken into account in order to ensure the existence of a stochastic solution. The proposed framework is finally exemplified by considering the modeling of various soft biological tissues, such as human brain and liver tissues.
Anna Mazzucato : Explicit parametrices for time-dependent Fokker-Planck equations
- Applied Math and Analysis ( 88 Views )We construct explicit approximate Green's functions of time-dependent, linear Fokker-Planck equations in terms of Dyson series, Taylor expansions, and exact commutator formulas. Our method gives an approximate solution that is accurate to arbitrary order in time in the short-time limit, and it can be extended to large time by bootstrapping. I will also present some numerical results showing that our algorithm works well also for degenerate equations such as those arising in pricing of contingent claims. This is joint work with Victor Nistor and Wen Cheng.
Matthew Jacobs : A fast approach to optimal transport: the back-and-forth method
- Applied Math and Analysis ( 195 Views )Given two probability measures and a transportation cost, the optimal transport problem asks to find the most cost efficient way to transport one measure to the other. Since its introduction in 1781 by Gaspard Monge, the optimal transport problem has found applications in logistics, economics, physics, PDEs, and more recently data science. However, despite sustained attention from the numerics community, solving optimal transport problems has been a notoriously difficult task. In this talk I will introduce the back-and-forth method, a new algorithm to efficiently solve the optimal transportation problem for a general class of strictly convex transportation costs. Given two probability measures supported on a discrete grid with n points, the method computes the optimal map in O(n log(n)) operations using O(n) storage space. As a result, the method can compute highly accurate solutions to optimal transportation problems on spatial grids as large as 4096 x 4096 and 384 x 384 x 384 in a matter of minutes. If time permits, I will demonstrate an extension of the algorithm to the simulation of a class of gradient flows. This talk is joint work with Flavien Leger.
Shi Jin : Asymptotic-preseving schemes for the Boltzmann equation and relative problems with multiple scales
- Applied Math and Analysis ( 100 Views )We propose a general framework to design asymptotic preserving schemes for the Boltzmann kinetic and related equations. Numerically solving these equations are challenging due to the nonlinear stiff collision (source) terms induced by small mean free or relaxation time. We propose to penalize the nonlinear collision term by a BGK-type relaxation term, which can be solved explicitly even if discretized implicitly in time. Moreover, the BGK-type relaxation operator helps to drive the density distribution toward the local Maxwellian, thus naturally imposes an asymptotic-preserving scheme in the Euler limit. The scheme so designed does not need any nonlinear iterative solver or the use of Wild Sum. It is uniformly stable in terms of the (possibly small) Knudsen number, and can capture the macroscopic fluid dynamic (Euler) limit even if the small scale determined by the Knudsen number is not numerically resolved. We will show how this idea can be applied to other collision operators, such as the Landau-Fokker-Planck operator, Ullenbeck-Urling model, and in the kinetic-fluid model of disperse multiphase flows, and can be implemented in the Monte-Carlo framework in which is time step is not limited by the possibly small mean free time.
Marina Iliopoulou : Algebraic aspects of harmonic analysis
- Applied Math and Analysis ( 96 Views )When we want to understand a geometric picture, finding the zero set of a polynomial hiding in it can be very helpful: it can reveal structure and allow computations. Polynomial partitioning, developed by Guth and Katz, is a technique to find such a nice algebraic hypersurface. Polynomial partitioning has revolutionised discrete incidence geometry in the recent years, thanks to the fact that interaction of lines with algebraic hypersurfaces is well-understood. Recently, however, Guth discovered agreeable interaction between tubes and algebraic hypersurfaces, and thus used polynomial partitioning to improve on the 3-dim restriction problem. In this talk, we will present polynomial partitioning via a discrete analogue of the Kakeya problem, and discuss its potential to be extensively used in harmonic analysis.
Jun Kitagawa : Free discontinuity regularity and stability in optimal transport
- Applied Math and Analysis ( 91 Views )Regularity of solutions in the optimal transport problem requires very rigid hypotheses (e.g., convexity of certain sets). When such conditions are not available, one can consider the question of partial regularity, in other words, the in-depth analysis of the structure of singular sets. In this talk, I will discuss the regularity of the set of ``free singularities`` which arise in an optimal transport problem with inner product cost, from a connected set to a disconnected set, along with the stability of such sets under suitable perturbations of the data involved. Some of these results are proven via a non-smooth implicit function theorem for convex functions, which is of independent interest. This talk is based on joint work with Robert McCann.
Mark Embree : CUR Matrix Factorizations: Algorithms, Analysis, Applications
- Applied Math and Analysis ( 98 Views )Interpolatory matrix factorizations provide alternatives to the singular value decomposition (SVD) for obtaining low-rank approximations; this class includes the CUR factorization, where the C and R matrices are subsets of columns and rows of the target matrix. While interpolatory approximations lack the SVD's optimality, their ingredients are easier to interpret than singular vectors: since they are copied from the matrix itself, they inherit the data's key properties (e.g., nonnegative/integer values, sparsity, etc.). We shall provide an overview of these approximate factorizations, describe how they can be analyzed using interpolatory projectors, and introduce a new method for their construction based on the Discrete Empirical Interpolation Method (DEIM). (This talk describes joint work with Dan Sorensen (Rice).)
Xiaoming Wang : Large Prandtl Number Behavior of the Boussinesq System
- Applied Math and Analysis ( 134 Views )We consider large Prandtl number behavior of the Boussinesq system for Rayleigh-B\'enard convection at large time. We first show that the global attractors of the Boussinesq system converge to that of the infinite Prandtl number model. This is accomplished via a generalization of upper semi-continuity property with respect to parameters of dissipative dynamical systems to the case of singular limit of two time scale problems of relaxation type. We then show that stationary statistical properties (in terms of invariant measures) of the Boussinesq system converge to that of the infinite Prandtl number model. In particular, we derive a new upper bound on heat transport in the vertical direction (the Nusselt number) for the Boussinesq system. The new upper bound agrees with the recent physically optimal upper bound on the infinite Prandtl number model at large Prandtl number. We will also comment on possible noise induced stability and its relation to the E-Mattingly-Sinai theory.
Yossi Farjoun : Solving Conservation Law and Balance Equations by Particle Management
- Applied Math and Analysis ( 101 Views )Conservation equations are at the heart of many interesting and important problems. Examples come from physics, chemistry, biology, traffic and many more. Analytically, hyperbolic equations have a beautiful structure due to the existence of characteristics. These provide the possibility of transforming a conservation PDE into a system of ODE and thus greatly reducing the computational effort required to solve such problems. However, even in one dimension, one encounters problems after a short time.
The most obvious difficulty that needs to be dealt with has to do with the creation of shocks, or in other words, the crossing of characteristics. With a particle based method one would like to avoid a situation when one particle overtakes a neighboring one. However, since shocks are inherent to many hyperbolic equations and relevant to the problems that one would like to solve, it would be good not to ``smooth away'' the shock but rather find a good representation of it and a good solution for the offending particles.
In this talk I will present a new particle based method for solving (one dimensional, scalar) conservation law equations. The guiding principle of the method is the conservative property of the underlying equation. The basic method is conservative, entropy decreasing, variation diminishing and exact away from shocks. A recent extension allows solving equations with a source term, and also provides ``exact'' solutions to the PDE. The method compares favorably to other benchmark solvers, for example CLAWPACK, and requires less computation power to reach the same resolution. A few examples will be shown to illustrate the method, with its various extensions. Due to the current limitation to 1D scalar, the main application we are looking at is traffic flow on a large network. Though we still hope to manage to extend the method to either systems or higher dimensions (each of these extensions has its own set of difficulties), I would be happy to discuss further possible applications or suggestions for extensions.
Yoshiaki Teramoto : Benard-Marangoni problem of heat convection with free surface
- Applied Math and Analysis ( 95 Views )When a fluid layer is heated from below with temperature larger than a certain critical value, the convective motion appears in the fluid. The convection caused by the thermocapillary effect is called Benard-Marangoni heat convection. The thermocapillary effect is the dependence of the surface tension on the temperature. Near a hot spot on a free surface of fluid a thermocapillary tangential stress generates a fluid motion. In this talk the mathematical model system for this convection is explained. The Oberbeck-Boussinesq approximation is used for the system and the upper boundary is a free surface with surface tension which depends on the temperature. After formulating the linearized problem around the conductive state, we derive the resolvent estimates which guarantee the sectorial property. Stationary and Hopf bifurcations (periodic solutions) are proved to exist depending on the parameters (Raylegh and Marangoni numbers).