# Applied/ACMS/absS18

# ACMS Abstracts: Spring 2018

### Thomas Fai (Harvard)

*The Lubricated Immersed Boundary Method*

Many real-world examples of fluid-structure interaction, including the transit of red blood cells through the narrow slits in the spleen, involve the near-contact of elastic structures separated by thin layers of fluid. The separation of length scales between these fine lubrication layers and the larger elastic objects poses significant computational challenges. Motivated by the challenge of resolving such multiscale problems, we introduce an immersed boundary method that uses elements of lubrication theory to resolve thin fluid layers between immersed boundaries. We apply this method to two-dimensional flows of increasing complexity, including eccentric rotating cylinders and elastic vesicles near walls in shear flow, to show its increased accuracy compared to the classical immersed boundary method. We present preliminary simulation results of cell suspensions, a problem in which near-contact occurs at multiple levels, such as cell-wall, cell-cell, and intracellular interactions, to highlight the importance of resolving thin fluid layers in order to obtain the correct overall dynamics.

### Michael Herty (RWTH-Aachen)

*Opinion Formation Models and Mean field Games Techniques*

Mean-Field Games are games with a continuum of players that incorporate the time dimension through a control-theoretic approach. Recently, simpler approaches relying on reply strategies have been proposed. Based on an example in opinion formation modeling we explore the link between differentiability notions and mean-field game approaches. For numerical purposes a model predictive control framework is introduced consistent with the mean-field game setting that allows for efficient simulation. Numerical examples are also presented as well as stability results on the derived control.

### Lee Panetta (Texas A&M)

*Traveling waves and pulsed energy emissions seen in numerical simulations of electromagnetic wave scattering by ice crystals*

The numerical simulation of single particle scattering of electromagnetic energy plays a fundamental role in remote sensing studies of the atmosphere and oceans, and in efforts to model aerosol "radiative forcing" processes in a wide variety of models of atmospheric and climate dynamics, I will briefly explain the main challenges in the numerical simulation of single particle scattering and describe how work with 3-d simulations of scattering of an incident Gaussian pulse, using a Pseudo-Spectral Time Domain method to numerically solve Maxwell’s Equations, led to an investigation of episodic bursts of energy that were observed at various points in the near field during the decay phase of the simulations. The main focus of the talk will be on simulations in dimensions 1 and 2, simple geometries, and a single refractive index (ice at 550 nanometers). The periodic emission of pulses is easy to understand and predict on the basis of Snell’s laws in the 1-d case considered. In much more interesting 2-d cases, simulations show traveling waves within the crystal that give rise to pulsed emissions of energy when they interact with each other or when they enter regions of high surface curvature. The time-dependent simulations give a more dynamical view of "photonic nanojets" reported earlier in steady-state simulations in other contexts, and of energy release in "morphology-dependent resonances."

### Haizhao Yang (National University of Singapore)

*A Unified Framework for Oscillatory Integral Transform: When to use NUFFT or Butterfly Factorization?*

This talk introduces fast algorithms of the matvec $g=Kf$ for $K\in \mathbb{C}^{N\times N}$, which is the discretization of the oscillatory integral transform $g(x) = \int K(x,\xi) f(\xi)d\xi$ with a kernel function $K(x,\xi)=\alpha(x,\xi)e^{2\pi i\Phi(x,\xi)}$, where $\alpha(x,\xi)$ is a smooth amplitude function , and $\Phi(x,\xi)$ is a piecewise smooth phase function with $O(1)$ discontinuous points in $x$ and $\xi$. A unified framework is proposed to compute $Kf$ with $O(N\log N)$ time and memory complexity via the non-uniform fast Fourier transform (NUFFT) or the butterfly factorization (BF), together with an $O(N)$ fast algorithm to determine whether NUFFT or BF is more suitable. This framework works for two cases: 1) explicite formulas for the amplitude and phase functions are known; 2) only indirect access of the amplitude and phase functions are available. Especially in the case of indirect access, our main contributions are: 1) an $O(N\log N)$ algorithm for recovering the amplitude and phase functions is proposed based on a new low-rank matrix recovery algorithm; 2) a new stable and nearly optimal BF with amplitude and phase functions in form of a low-rank factorization (IBF-MAT) is proposed to evaluate the matvec $Kf$. Numerical results are provided to demonstrate the effectiveness of the proposed framework.

### Eric Keaveny (Imperial College London)

*Linking the micro- and macro-scales in populations of swimming cells*

Swimming cells and microorganisms are as diverse in their collective dynamics as they are in their individual shapes and swimming mechanisms. They are able to propel themselves through simple viscous fluids, as well as through more complex environments where they must interact with other microscopic structures. In this talk, I will describe recent simulations that explore the connection between dynamics at the scale of the cell with that of the population in the case where the cells are sperm. In particular, I will discuss how the motion of the sperm’s flagella can greatly impact the overall dynamics of their suspensions. Additionally, I will discuss how in complex environments, the density and stiffness of structures with which the cells interact impact the effective diffusion of the population.

### Molei Tao (Georgia Tech)

*Explicit high-order symplectic integration of nonseparable Hamiltonians: algorithms and long time performance*

Symplectic integrators preserve the phase-space volume and have favorable performances in long time simulations. Methods for an explicit symplectic integration have been extensively studied for separable Hamiltonians (i.e., H(q,p)=K(p)+V(q)), and they lead to both accurate and efficient simulations. However, nonseparable Hamiltonians also model important problems, such as non-Newtonian mechanics and nearly integrable systems in action-angle coordinates. Unfortunately, implicit methods had been the only available symplectic approach for general nonseparable systems.

This talk will describe a recent result that constructs explicit and arbitrary high-order symplectic integrators for arbitrary Hamiltonians. Based on a mechanical restraint that binds two copies of phase space together, these integrators have good long time performance. More precisely, based on backward error analysis, KAM theory, and some additional multiscale analysis, a pleasant error bound is established for integrable systems. This bound is then demonstrated on a conceptual example and the Schwarzschild geodesics problem. For nonintegrable systems, some numerical experiments with the nonlinear Schrodinger equation will be discussed.

### Boualem Khouider (UVic)

*Title TBA*

Abstract TBA