Seminars and Colloquia by Series

Parareal methods for highly oscillatory ordinary differential equations

Series
Applied and Computational Mathematics Seminar
Time
Monday, April 13, 2015 - 14:00 for 1 hour (actually 50 minutes)
Location
Skiles 005
Speaker
Seong Jun KimGeorgia Tech
We introduce a new parallel in time (parareal) algorithm which couples multiscale integrators with fully resolved fine scale integration and computes highly oscillatory solutions for a class of ordinary differential equations in parallel. The algorithm computes a low-cost approximation of all slow variables in the system. Then, fast phase-like variables are obtained using the parareal iterative methodology and an alignment algorithm. The method may be used either to enhance the accuracy and range of applicability of the multiscale method in approximating only the slow variables, or to resolve all the state variables. The numerical scheme does not require that the system is split into slow and fast coordinates. Moreover, the dynamics may involve hidden slow variables, for example, due to resonances.

Control of Oscillators, Temporal Homogenization, and Energy Harvest by Super-Parametric Resonance

Series
Applied and Computational Mathematics Seminar
Time
Monday, April 6, 2015 - 14:00 for 1 hour (actually 50 minutes)
Location
Skiles 005
Speaker
Prof. Molei TaoGeorgia Tech School of Math.
We show how to control an oscillator by periodically perturbing its stiffness, such that its amplitude follows an arbitrary positive smooth function. This also motivates the design of circuits that harvest energies contained in infinitesimal oscillations of ambient electromagnetic fields. To overcome a key obstacle, which is to compensate the dissipative effects due to finite resistances, we propose a theory that quantifies how small/fast periodic perturbations affect multidimensional systems. This results in the discovery of a mechanism that reduces the resistance threshold needed for energy extraction, based on coupling a large number of RLC circuits.

A method of computation of 2D Fourier transforms and diffraction integrals with applications in vision science

Series
Applied and Computational Mathematics Seminar
Time
Monday, March 30, 2015 - 14:05 for 1 hour (actually 50 minutes)
Location
Skiles 005
Speaker
Professor Andrei Martinez-FinkelshteinUniversity of Almería
The importance of the 2D Fourier transform in mathematical imaging and vision is difficult to overestimate. For instance, the impulse response of an optical system can be defined in terms of diffraction integrals, that are in turn Fourier transforms of a function on a disk. There are several popular competing approaches used to calculate diffraction integrals, such as the extended Nijboer-Zernike (ENZ) theory. In this talk, an alternative efficient method of computation of two dimensional Fourier-type integrals based on approximation of the integrand by Gaussian radial basis functions is discussed. Its outcome is a rapidly converging series expansion for the integrals, allowing for their accurate calculation. The proposed method yields a reliable and fast scheme for simultaneous evaluation of such kind of integrals for several values of the defocus parameter, as required in the characterization of the through-focus optics.

A Non-convex Approach for Signal and Image Processing

Series
Applied and Computational Mathematics Seminar
Time
Tuesday, March 24, 2015 - 11:00 for 1 hour (actually 50 minutes)
Location
Skiles 005
Speaker
Prof. Yifei LouUT Dallas
A fundamental problem in compressed sensing (CS) is to reconstruct a sparsesignal under a few linear measurements far less than the physical dimensionof the signal. Currently, CS favors incoherent systems, in which any twomeasurements are as little correlated as possible. In reality, however, manyproblems are coherent, in which case conventional methods, such as L1minimization, do not work well. In this talk, I will present a novelnon-convex approach, which is to minimize the difference of L1 and L2 norms(L1-L2) in order to promote sparsity. Efficient minimization algorithms areconstructed and analyzed based on the difference of convex functionmethodology. The resulting DC algorithms (DCA) can be viewed as convergentand stable iterations on top of L1 minimization, hence improving L1 consistently. Through experiments, we discover that both L1 and L1-L2 obtain betterrecovery results from more coherent matrices, which appears unknown intheoretical analysis of exact sparse recovery. In addition, numericalstudies motivate us to consider a weighted difference model L1-aL2 (a>1) todeal with ill-conditioned matrices when L1-L2 fails to obtain a goodsolution. An extension of this model to image processing will be alsodiscussed, which turns out to be a weighted difference of anisotropic andisotropic total variation (TV), based on the well-known TV model and naturalimage statistics. Numerical experiments on image denoising, imagedeblurring, and magnetic resonance imaging (MRI) reconstruction demonstratethat our method improves on the classical TV model consistently, and is onpar with representative start-of-the-art methods.

Numerical schemes for stochastic backscatter in the inverse cascade of quasi-geostrophic turbulence

Series
Applied and Computational Mathematics Seminar
Time
Monday, March 23, 2015 - 14:05 for 1 hour (actually 50 minutes)
Location
Skiles 005
Speaker
Yoonsang LeeCourant Institute of Mathematical Sciences
Backscatter is the process of energy transfer from small to large scales in turbulence; it is crucially important in the inverse energy cascades of two-dimensional and quasi-geostrophic turbulence, where the net transfer of energy is from small to large scales. A numerical scheme for stochastic backscatter in the two-dimensional and quasi-geostrophic inverse kinetic energy cascades is developed and analyzed. Its essential properties include a local formulation amenable to implementation in finite difference codes and non-periodic domains, smooth behavior at the coarse grid scale, and realistic temporal correlations, which allows detailed numerical analysis, focusing on the spatial and temporal correlation structure of the modeled backscatter. The method is demonstrated in an idealized setting of quasi-geostrophic turbulence using a low-order finite difference code, where it produces a good approximation to the results of a spectral code with more than 5 times higher nominal resolution. This is joint work with I. Grooms and A. J. Majda

Implicit interface boundary integral methods

Series
Applied and Computational Mathematics Seminar
Time
Friday, March 13, 2015 - 14:00 for 1 hour (actually 50 minutes)
Location
Skiles 168
Speaker
Richard TsaiUniversity of Texas at Austin
I will present a new approach for computing boundary integrals that are defined on implicit interfaces, without the need of explicit parameterization. A key component of this approach is a volume integral which is identical to the integral over the interface. I will show results applying this approach to simulate interfaces that evolve according to Mullins-Sekerka dynamics used in certain phase transition problems. I will also discuss our latest results in generalization of this approach to summation of unstructured point clouds and regularization of hyper-singular integrals.

Transition path processes and coarse-graining of stochastic system

Series
Applied and Computational Mathematics Seminar
Time
Monday, March 9, 2015 - 14:00 for 1 hour (actually 50 minutes)
Location
Skiles 005
Speaker
Prof. Jianfeng LuDuke University
Understanding rare events like transitions of chemical system from reactant to product states is a challenging problem due to the time scale separation. In this talk, we will discuss some recent progress in mathematical theory of transition paths. In particular, we identify and characterize the stochastic process corresponds to transition paths. The study of transition path process helps to understand the transition mechanism and provides a framework to design and analyze numerical approaches for rare event sampling and simulation.

Existence and Stability of Radially Symmetric Solutions to the Swift--Hohenberg Equation

Series
Applied and Computational Mathematics Seminar
Time
Monday, March 2, 2015 - 14:00 for 1 hour (actually 50 minutes)
Location
Skiles 005
Speaker
Professor Scott McCallaMontana State University
The existence, stability, and bifurcation structure of localized radially symmetric solutions to the Swift--Hohenberg equation is explored both numerically through continuation and analytically through the use of geometric blow-up techniques. The bifurcation structure for these solutions is elucidated by formally treating the dimension as a continuous parameter in the equations. This reveals a family of solutions with an anomalous amplitude scaling that is far larger than expected from a formal scaling in the far field. One key advantage of the geometric blow-up techniques is that a priori knowledge of this scaling is unnecessary as it naturally emerges from the construction. The stability of these patterned states will also be discussed.

Nonnegative Inverse Eigenvalue and Singular Value Problems

Series
Applied and Computational Mathematics Seminar
Time
Monday, February 16, 2015 - 14:00 for 1 hour (actually 50 minutes)
Location
Skiles 005
Speaker
Prof. Matthew LinNational Chung Cheng University, Georgia Tech

Please Note: Reference[1] Moody T. Chu , Nonnegative Inverse Eigenvalue and Singular Value Problems, SIAM J. Numer. Anal (1992).[2] Wei Ma and Zheng-J. Bai, A regularized directional derivative-based Newton method for inverse singular value problems, Inverse Problems (2012).

Nonnegative inverse eigenvalue and singular value problems have been a research focus for decades. It is true that an inverse problem is trivial if the desired matrix is not restricted to any structure. This talk is to present two numerical procedures, based on a conquering procedure and an alternating projection process, to solve inverse eigenvalue and singular value problems for nonnegative matrices, respectively. In theory, we also discuss the existence of nonnegative matrices subject to prescribed eigenvalues and singular values. Though the focus of this talk is on inverse eigenvalue and singular value problems with nonnegative entries, the entire procedure can be straightforwardly applied to other types of structure with no difficulty.

Methods to solve Kohn-Sham equations for electron density

Series
Applied and Computational Mathematics Seminar
Time
Monday, February 9, 2015 - 14:00 for 1 hour (actually 50 minutes)
Location
Skiles 005
Speaker
Timo EirolaAalto University, Helsinki, Finland
We consider three different approaches to solve the equations for electron density around nuclei particles. First we study a nonlinear eigenvalue problem and apply Quasi-Newton methods to this. In many cases they turn to behave better than the Pulay mixer, which widely used in physics community. Second we reformulate the problem as a minimization problem on a Stiefel manifold. One that formed from mxn matrices with orthonormal columns. Then for Quasi-Newton techniques one needs to transfer the secant conditions to the new tangent space, when moving on the manifold. We also consider nonlinear conjugate gradients in this setting. This minimization approach seems to work well especially for metals, which are known to be hard. Third (if time permits) we add temperature (the first two are for ground state). This means that we need to include entropy in the energy and optimize also with respect to occupation numbers. Joint work with Kurt Baarman and Ville Havu.

Pages