Michael Gratton : Transient and self-similar dynamics in thin film coarsening
- Applied Math and Analysis ( 144 Views )Coarsening is the phenomenon where many objects (water drops, molecular islands, particles in a freezing liquid) becoming smaller in number but larger in size in an orderly way. This talk will examine modeling one such system, nanoscopic liquid drops, through three models: a PDE for the fluid, a coarsening dynamical system for the drops, and an LSW-type ensemble model for the distribution of drops. We will find self-similar solutions for the drop population valid for intermediate times and discuss transient effects that can delay the self-similar scaling.
Pierre Degond : Asymptotic-Preserving numerical methods for variable-scale problems. Examples from fluids and plasma dynamics
- Applied Math and Analysis ( 100 Views )Multiscale problems are often treated via asymptotic of homogenization techniques: one first determines the asymptotic limit and then finds an appropriate numerical methods to solve it. Variable scale problems which exhibit a continuous variation of the perturbation parameter from a finite to an infinitesimal value cannot be solved by this method alone. They require the coupling of the asymptotic problem to the original one across the region of scale variation. This coupling is often quite complex and lacks robustness. Asymptotic-Preserving methods represent an alternative to the coupling strategy and provide a way to resolve the original problem without resorting to its asymptotic limit. They provide a systematic methodology to resolve multiscale problems even in situations where the asymptotic limit is quite complex. We will provide examples of this methodology for the treatment of the low-Mach number regime, of quasineutrality in plasmas, large magnetic fields or strong anisotropy in diffusion equations.
Alexei Novikov : Diffusion in fluid flows
- Applied Math and Analysis ( 108 Views )Enhacement of diffusion and reaction processes by advection is a classical subject that has been extensively studied by both physicists and mathematicians. In this talk I will consider reaction-diffusion on a bounded domain in the presence of a fast incompressible flow. I will also consider convective systems in a bounded domain, in which viscous fluids described by the Stokes system are coupled using the Boussinesq approximation to a reaction-advection-diffusion equation for the temperature. For such systems enhacement of diffusion implies, in particular, that the explosion threshold for power-like nonlinearities tends to infinity in the large Rayleigh number limit.
Karin Leiderman : A Spatial-Temporal Model of Platelet Deposition and Blood Coagulation Under Flow
- Applied Math and Analysis ( 141 Views )In the event of a vascular injury, a blood clot will form to prevent bleeding. This response involves two intertwined processes: platelet aggregation and coagulation. Activated platelets are critical to coagulation in that they provide localized reactive surfaces on which many of the coagulation reactions occur. The final product from the coagulation cascade directly couples the coagulation system to platelet aggregation by acting as a strong activator of platelets and cleaving blood-borne fibrinogen into fibrin which then forms a mesh to help stabilize platelet aggregates. Together, the fibrin mesh and the platelet aggregates comprise a blood clot, which in some cases, can grow to occlusive diameters. Transport of coagulation proteins to and from the vicinity of the injury is controlled largely by the dynamics of the blood flow. It is crucial to learn how blood flow affects the growth of clots, and how the growing masses, in turn, feed back and affect the fluid motion. We have developed the first spatial-temporal model of platelet deposition and blood coagulation under flow that includes detailed descriptions of the coagulation biochemistry, chemical activation and deposition of blood platelets, as well as the two-way interaction between the fluid dynamics and the growing platelet mass.
Ingrid Daubechies : Surface Comparison With Mass Transportation
- Applied Math and Analysis ( 144 Views )In many applications, ranging from computer animation to biology, one wants to quantify how similar two surfaces are to each other. In the last few years, the Gromov-Haussdorff distance has been applied to this problem; this gives good results, but turns out to be very heavy computationally. This talk proposes a different approach, in which (disk-like) 2-dimensional surfaces (typically embedded in 3-dimensional Euclidean space) are first mapped conformally to the unit disk, and the corresponding conformal densities are then compared via optimal mass transportation,. This mass transportation problem differs from the standard case in that we require the solution to be invariant under global Moebius transformations. The metric we construct also defines meaningful intrinsic distances between pairs of "patches" in the two surfaces, which allows automatic alignment of the surfaces. Numerical experiments on "real-life" surfaces to demonstrate possible applications in natural sciences will be shown as well. This is joint work with Yaron Lipman.
Shahed Sharif : Who wants to be a millionaire?
- Applied Math and Analysis ( 93 Views )The Mordell-Weil theorem shows that the rational points on an elliptic curve defined over the field of rational numbers is a finitely generated abelian group. The Birch and Swinnerton-Dyer conjecture relates the rank of this group to a number of analytic and algebraic invariants of the curve. (More generally it considers an elliptic curve defined over a number field.) The conjecture is one of the Millennium Prize problems and the Clay Institute is offering a reward of 1 million dollars for a solution. This talk will be an introduction to the conjecture. In following weeks we will have lectures explaining each of the terms in the formula.
Andrea Bertozzi : Geometry based image processing - a survey of recent results
- Applied Math and Analysis ( 101 Views )I will present a survey of recent results on geometry-based image processing. The topics will include wavelet-based diffuse interface methods, pan sharpening and hyperspectral sharpening, and sparse image representation.
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.
Alexandr Labovschii : High accuracy numerical methods for fluid flow problems and turbulence modeling
- Applied Math and Analysis ( 98 Views )We present several high accuracy numerical methods for fluid flow problems and turbulence modeling.
First we consider a stabilized finite element method for the Navier-Stokes equations which has second order temporal accuracy. The method requires only the solution of one linear system (arising from an Oseen problem) per time step.
We proceed by introducing a family of defect correction methods for the time dependent Navier-Stokes equations, aiming at higher Reynolds' number. The method presented is unconditionally stable, computationally cheap and gives an accurate approximation to the quantities sought.
Next, we present a defect correction method with increased time accuracy. The method is applied to the evolutionary transport problem, it is proven to be unconditionally stable, and the desired time accuracy is attained with no extra computational cost.
We then turn to the turbulence modeling in coupled Navier-Stokes systems - namely, MagnetoHydroDynamics. We consider the mathematical properties of a model for the simulation of the large eddies in turbulent viscous, incompressible, electrically conducting flows. We prove existence, uniqueness and convergence of solutions for the simplest closed MHD model. Furthermore, we show that the model preserves the properties of the 3D MHD equations.
Lastly, we consider the family of approximate deconvolution models (ADM) for turbulent MHD flows. We prove existence, uniqueness and convergence of solutions, and derive a bound on the modeling error. We verify the physical properties of the models and provide the results of the computational tests.
John Stockie : Porous immersed boundaries
- Applied Math and Analysis ( 116 Views )Porous, deformable membranes are encountered in a wide range of applications including red blood cells, vesicles, porous wave makers, and parachutes. The "immersed boundary method" has already proven to be a versatile and robust approach for simulating the interaction of impermeable, elastic structures with an incompressible fluid flow. We demonstrate how to extend the method to handle porous boundaries by incorporating an explicit porous slip velocity that is determined by Darcy's law. We derive a simple, radially-symmetric exact solution, which is then used to validate numerical simulations of porous membranes in two dimensions.
Guillaume Bal : Some convergence results in equations with random coefficients
- Applied Math and Analysis ( 102 Views )The theory of homogenization for equations with random coefficients is now quite well-developed. What is less studied is the theory for the correctors to homogenization, which asymptotically characterize the randomness in the solution of the equation and as such are important to quantify in many areas of applied sciences. I will present recent results in the theory of correctors for elliptic and parabolic problems and briefly mention how such correctors may be used to improve reconstructions in inverse problems. Homogenized (deterministic effective medium) solutions are not the only possible limits for solutions of equations with highly oscillatory random coefficients as the correlation length in the medium converges to zero. When fluctuations are sufficiently large, the limit may take the form of a stochastic equation and stochastic partial differential equations (SPDE) are routinely used to model small scale random forcing. In the very specific setting of a parabolic equation with large, Gaussian, random potential, I will show the following result: in low spatial dimensions, the solution to the parabolic equation indeed converges to the solution of a SPDE, which however needs to be written in a (somewhat unconventional) Stratonovich form; in high spatial dimension, the solution to the parabolic equation converges to a homogenized (hence deterministic) equation and randomness appears as a central limit-type corrector. One of the possible corollaries for this result is that SPDE models may indeed be appropriate in low spatial dimensions but not necessarily in higher spatial dimensions.
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.
Michael Minion : Parallel in time integration: parareal and deferred corrections
- Applied Math and Analysis ( 101 Views )The efficient parallelization of numerical methods for ordinary or partial differential equations in the temporal direction is an intriguing possibility that has of yet not been fully realized despite decades of investigation. For partial differential equations, virtually all large scale computations now employ parallelization across space, and there are freely available computational tools and libraries to aid in the development of spatially parallelized codes. Conversely, parallelization in the temporal direction is rarely even considered. I will discuss a relatively recent parallel strategy called the parareal algorithm that has generated a renewed wave of interest in time parallelization. I will show how the iterative structure of the parareal algorithm can be interpreted as a particular form of deferred corrections and then present a modified parareal strategy based on spectral deferred corrections that can significantly reduce the computational cost of the method. Finally I will make some observations as to why parallel in time methods may be attractive in the future.
Margaret Beck : Nonlinear stability of time-periodic viscous shocks
- Applied Math and Analysis ( 126 Views )In order to understand the nonlinear stability of many types of time-periodic traveling waves on unbounded domains, one must overcome two main difficulties: the presence of zero eigenvalues that are embedded in the continuous spectrum and the time-periodicity of the associated linear operator. I will outline these issues and show how they can be overcome in the context of time-periodic Lax shocks in systems of viscous conservation laws. The method involves the development of a contour integral representation of the linear evolution, similar to that of a strongly continuous semigroup, and detailed pointwise estimates on the resultant Greens function, which are sufficient for proving nonlinear stability under the necessary assumption of spectral stability.
Andrew J. Bernoff : Domain Relaxation in Langmuir Films
- Applied Math and Analysis ( 132 Views )We report on an experimental and theoretical study of a molecularly thin polymer Langmuir layers on the surface of a Stokesian subfluid. Langmuir layers can have multiple phases (fluid, gas, liquid crystal, isotropic or anisotropic solid); at phase boundaries a line tension force is observed. By comparing theory and experiment we can estimate this line tension. We first consider two co-existing fluid phases; specifically a localized phase embedded in an infinite secondary phase. When the localized phase is stretched (by a transient stagnation flow), it takes the form of a bola consisting of two roughly circular reservoirs connected by a thin tether. This shape relaxes to the minimum energy configuration of a circular domain. The tether is never observed to rupture, even when it is more than a hundred times as long as it is thin. We model these experiments by taking previous descriptions of the full hydrodynamics (primarily those of Stone & McConnell and Lubensky & Goldstein), identifying the dominant effects via dimensional analysis, and reducing the system to a more tractable form. The result is a free boundary problem where motion is driven by the line tension of the domain and damped by the viscosity of the subfluid. The problem has a boundary integral formulation which allows us to numerically simulate the tether relaxation; comparison with the experiments allows us to estimate the line tension in the system. We also report on incorporating dipolar repulsion into the force balance and simulating the formation of "labyrinth" patterns.
Michael Siegel : Modeling, analysis, and computations of the influence of surfactant on the breakup of bubbles and drops in a viscous fluid
- Applied Math and Analysis ( 114 Views )We present an overview of experiments, numerical simulations, and mathematical analysis of the breakup of a low viscosity drop in a viscous fluid, and consider the role of surface contaminants, or surfactants, on the dynamics near breakup. As part of our study, we address a significant difficulty in the numerical computation of fluid interfaces with soluble surfactant that occurs in the important limit of very large values of bulk Peclet number Pe. At the high values of Pe in typical fluid-surfactant systems, there is a narrow transition layer near the drop surface or interface in which the surfactant concentration varies rapidly, and its gradient at the interface must be determined accurately to find the drops dynamics. Accurately resolving the layer is a challenge for traditional numerical methods. We present recent work that uses the narrowness of the layer to develop fast and accurate `hybrid numerical methods that incorporate a separate analytical reduction of the dynamics within the transition layer into a full numerical solution of the interfacial free boundary problem.
Lucy Zhang : Modeling and Simulations of Fluid and Deformable-Structure Interactions in Bio-Mechanical Systems
- Applied Math and Analysis ( 154 Views )Fluid-structure interactions exist in many aspects of our daily lives. Some biomedical engineering examples are blood flowing through a blood vessel and blood pumping in the heart. Fluid interacting with moving or deformable structures poses more numerical challenges for its complexity in dealing with transient and simultaneous interactions between the fluid and solid domains. To obtain stable, effective, and accurate solutions is not trivial. Traditional methods that are available in commercial software often generate numerical instabilities.
In this talk, a novel numerical solution technique, Immersed Finite Element Method (IFEM), is introduced for solving complex fluid-structure interaction problems in various engineering fields. The fluid and solid domains are fully coupled, thus yield accurate and stable solutions. The variables in the two domains are interpolated via a delta function that enables the use of non-uniform grids in the fluid domain, which allows the use of arbitrary geometry shapes and boundary conditions. This method extends the capabilities and flexibilities in solving various biomedical, traditional mechanical, and aerospace engineering problems with detailed and realistic mechanics analysis. Verification problems will be shown to validate the accuracy and effectiveness of this numerical approach. Several biomechanical problems will be presented: 1) blood flow in the left atrium and left atrial appendage which is the main source of blood in patients with atrial fibrillation. The function of the appendage is determined through fluid-structure interaction analysis, 2) examine blood cell and cell interactions under different flow shear rates. The formation of the cell aggregates can be predicted when given a physiologic shear rate.
Elizabeth L. Bouzarth : Modeling Biologically Inspired Fluid Flow Using RegularizedSingularities and Spectral Deferred Correction Methods
- Applied Math and Analysis ( 146 Views )The motion of primary nodal cilia present in embryonic development resembles that of a precessing rod. Implementing regularized singularities to model this fluid flow numerically simulates a situation for which colleagues have exact mathematical solutions and experimentalists have corresponding laboratory studies on both the micro- and macro-scales. Stokeslets are fundamental solutions to the Stokes equations, which act as external point forces when placed in a fluid. By strategically distributing regularized Stokeslets in a fluid domain to mimic an immersed boundary (e.g., cilium), one can compute the velocity and trajectory of the fluid at any point of interest. The simulation can be adapted to a variety of situations including passive tracers, rigid bodies and numerous rod structures in a fluid flow generated by a rod, either rotating around its center or its tip, near a plane. The exact solution allows for careful error analysis and the experimental studies provide new applications for the numerical model. Spectral deferred correction methods are used to alleviate time stepping restrictions in trajectory calculations. Quantitative and qualitative comparisons to theory and experiment have shown that a numerical simulation of this nature can generate insight into fluid systems that are too complicated to fully understand via experiment or exact numerical solution independently.
Jason Metcalf : Strichartz estimates on Schwarzschild black hole backgrounds
- Applied Math and Analysis ( 131 Views )In this talk, we will present some recent work on dispersive estimates for wave equations on Schwarzschild black hole backgrounds. We in particular will discuss Strichartz estimates and localized energy estimate. This is from a joint work with Jeremy Marzuola, Daniel Tataru, and Mihai Tohaneanu.
Ralph Smith : Model Development and Control Design for High Performance Nonlinear Smart Material Systems
- Applied Math and Analysis ( 141 Views )High performance transducers utilizing piezoceramic, electrostrictive, magnetostrictive or shape memory elements offer novel control capabilities in applications ranging from flow control to precision placement for nanoconstruction. To achieve the full potential of these materials, however, models, numerical methods and control designs which accommodate the constitutive nonlinearities and hysteresis inherent to the compounds must be employed. Furthermore, it is advantageous to consider material characterization, model development, numerical approximation, and control design in concert to fully exploit the novel sensor and actuator capabilities of these materials in coupled systems.
In this presentation, the speaker will discuss recent advances in the development of model-based control strategies for high performance smart material systems. The presentation will focus on the development of unified nonlinear hysteresis models, inverse compensators, reduced-order approximation techniques, and nonlinear control strategies for high precision or high drive regimes. The range for which linear models and control methods are applicable will also be outlined. Examples will be drawn from problems arising in structural acoustics, high speed milling, deformable mirror design, artificial muscle development, tendon design to minimize earthquake damage, and atomic force microscopy.
Seung-Yeal Ha : Uniform L^p-stability problem for the Boltzmann equation
- Applied Math and Analysis ( 146 Views )The Boltzmann equation governs the dynamics of a dilute gas. In this talk, I will address the L^p-stability problem of the Boltzmann equation near vacuum and a global Maxwellian. In a close-to-vacuum regime, I will explain the nonlinear functional approach motivated by Glimm's theory in hyperbolic conservation laws. This functional approach yields the uniform L^1-stability estimate. In contrast, in a close-to-global maxwellian regime, I will present the L^2-stability theory which establishes the uniform L^2-stability of several classical solutions.
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.
Paolo Aluffi : Chern class identities from string theory
- Applied Math and Analysis ( 143 Views )(joint with Mboyo Esole) String theory considerations lead to a non-trivial identity relating the Euler characteristics of an elliptically fibered Calabi-Yau fourfold and of certain related surfaces. After giving a very sketchy idea of the physics arguments leading to this identity, I will present a Chern class identity which confirms it, generalizing it to arbitrary dimension and to varieties that are not necessarily Calabi-Yaus. The relevant loci are singular, and this plays a key role in the identity.
Svetlana Tlupova : Numerical Solutions of Coupled Stokes and Darcy Flows Based on Boundary Integrals
- Applied Math and Analysis ( 151 Views )Coupling between free fluid flow and flow through porous media is important in many industrial applications, such as filtration, underground water flow in hydrology, oil recovery in petroleum engineering, fluid flow through body tissues in biology, to name a few.
Stokes flows appear in many applications where the fluid viscosity is high and/or the velocity and length scales are small. The flow through a porous medium can be described by Darcy's law. A region that contains both requires a careful coupling of these different systems at the interface through appropriate boundary conditions.
Our objective is to develop a method based on the boundary integral formulation for computing the fluid/porous medium problem with higher accuracy using fundamental solutions of Stokes and Darcy's equations. We regularize the kernels to remove the singularity for stability of numerical calculations and eliminate the largest error for higher accuracy.