Applied and Computational Mathematics Seminar

Department of Mathematics and Statistics



  

Spring 2022 Schedule
Parker Hall 358, Friday 3p.m. (CST)

For any questions or requests, please contact Selim Sukhtaiev (szs0266@auburn.edu)

 

Speaker Institution Date
Wenjing Liao Georgia Tech February 25
Rachidi Salako  University of Nevada, Las Vegas March 18 
Todd Arbogast The University of Texas at Austin March 25
Rihui Lan  University of South Carolina April 1 
 Chrysoula Tsogka University of California, Merced  April 15
 Rosemary Renault  Arizona State University  April 22 
 Jason Bramburger  George Mason University  May 6
Yuming Paul Zhang   UC San Diego  June 15
     
     
     
     
     
     

 

 

Jason Bramburger

Title:Deep learning of conjugate mappings

Abstract: Despite many of the most common chaotic dynamical systems being continuous in time, it is through discrete time mappings that much of the understanding of chaos is formed. Henri Poincaré first made this connection by tracking consecutive iterations of the continuous flow with a lower-dimensional, transverse subspace. The mapping that iterates the dynamics through consecutive intersections of the flow with the subspace is now referred to as a Poincaré map, and it is the primary method available for interpreting and classifying chaotic dynamics. Unfortunately, in all but the simplest systems, an explicit form for such a mapping remains outstanding. In this talk Jason will present a method of discovering explicit Poincaré mappings using deep learning to construct an invertible coordinate transformation into a conjugate representation where the dynamics are governed by a relatively simple chaotic mapping. The invertible change of variable is based on an autoencoder, which allows for dimensionality reduction, and has the advantage of classifying chaotic systems using the equivalence relation of topological conjugacies. We illustrate with low-dimensional systems such as the Rössler and Lorenz systems, while also demonstrating the utility of the method on the infinite-dimensional Kuramoto-Sivashinsky equation.
 

 

 

Rosemary Renault

Title:Examining the Feasibility of Computational Methods for Inverting Gravity and Magnetic Data: Numerical Methods and Joint Inversion

 

 

 

Chrysoula Tsogka

Title:Fast Signal Recovery from Quadratic Measurements.

Abstract:

We present a novel approach for recovering a sparse signal from quadratic measurements corresponding to a rank-one tensorization of the data vector. Such quadratic measurements, referred to as interferometric or cross-correlated data, naturally arise in many fields such as remote sensing, spectroscopy, holography and seismology. Compared to the sparse signal recovery problem that uses linear measurements, the unknown in this case is a matrix, $X=\rho \rho^*$, formed by the cross correlations of $\rho \in C^K$. This creates a bottleneck for the inversion since the number of unknowns grows quadratically in $K$. The main idea of the proposed approach is to reduce the dimensionality of the problem by recovering only the diagonal elements of the unknown matrix, $| \rho_i|^2$, $i=1,\ldots,K$.  The contribution of the off-diagonal terms $\rho_i \rho_j^*$ for $i \neq j$ to the data is treated as noise and is absorbed using the Noise Collector approach introduced in [Moscoso et al, The noise collector for sparse recovery in high dimensions, PNAS 117 (2020)]. With this strategy, we recover the unknown  by solving a convex linear problem whose cost is similar to the one that uses linear measurements. The proposed approach provides exact support recovery when the data is not too noisy, and there are no false positives for any level of noise.

 

Rihui Lan

Title:High-Order Multirate Explicit Time-Stepping Schemes for the Baroclinic-Barotropic Split Dynamics in Primitive Equations

Abstract:

In order to treat the multiple time scales of ocean dynamics in an efficient manner, the baroclinic-barotropic splitting technique has been widely used for solving the primitive equations for ocean modeling. Based on the framework of strong stability-preserving Runge-Kutta  approach, we propose two high-order  multirate  explicit time-stepping schemes (SSPRK2-SE and SSPRK3-SE) for the resulting split system in this paper. The  proposed schemes allow for a large time step  to be used for  the three-dimensional  baroclinic (slow) mode and a small time step for the two-dimensional barotropic (fast)  mode,  in which each of the two mode solves just need to satisfy their respective CFL conditions for numerical stability. Specifically, at each time step, the  baroclinic velocity is first computed  by advancing the baroclinic mode and fluid thickness of the system with the large time-step  and the assistance of  some intermediate approximations of the baroctropic mode obtained by substepping with the small time step; then the barotropic velocity is corrected by using  the small time step to re-advance the barotropic mode under an improved barotropic forcing produced by interpolation of the forcing terms from the preceding baroclinic mode solves; lastly, the fluid thickness is updated by coupling the baroclinic and barotropic velocities. Additionally, numerical inconsistencies on the discretized sea surface height caused by the mode splitting  are  relieved  via a reconciliation process with carefully calculated flux deficits. Two benchmark tests from the ``MPAS-Ocean" platform are  carried out to numerically demonstrate  the performance and parallel scalability of the proposed SSPRK-SE schemes.

 

 

Todd Arbogast (with Chieh-Sen Huang, Ming-Hsien Kuo, and Xikai Zhao)

Title:Polynomial and RBF Finite Volume WENO Reconstructions for Solving Conservation Laws

Abstract: We analyze weighted essentially nonoscillatory (WENO) reconstructions of functions from finite volume average values, basing the approximations on either polynomials or radial basis functions (RBFs). In the polynomial case, we present new multi-level WENO reconstructions with adaptive order (WENO-AO). We give conditions under which the polynomial reconstructions achieve optimal order accuracy for both smooth solutions and solutions with discontinuities. The theory of existence of finite volume RBF approximations is developed, and numerical evidence is given to show that the RBF approximation converges to the same order as a polynomial approximation when the RBF is infinitely differentiable. Specific multiquadric RBFs on stencils of 2 and 3 mesh cells are proven to have this convergence property. These WENO-AO and RBF-WENO-AO reconstructions are applied to develop finite volume schemes for hyperbolic conservation laws on nonuniform meshes over multiple space dimensions. Numerical examples show that the schemes maintain proper accuracy and achieve the essentially non-oscillatory property when solving hyperbolic conservation laws.

 

 

Rachidi B. Salako

Title:Coexistence of pathogens of a multi-strain epidemic model.

Abstract: We study some multi-strain reaction-diffusion epidemic model dynamics with spatially heterogeneous infection and recovery rates. The nonexistence of coexistence endemic equilibrium (EE) solutions and the final extinction of one or multiple strains are proved under some proper assumptions on the model’s parameters. However, when some of these assumptions fail to hold for the case of the two strains model, we show that both strains of the disease may coexist, and solutions of the initial value problem stabilize at the unique coexistence endemic equilibrium. Finally, the asymptotic behaviors of the coexistence endemic equilibrium solutions, as the diffusion rate of the populations approaches zero, are also investigated, where the spatial segregation of multiple strains is found and determined.

 

 

Wenjing Liao

Title: Statistical learning theory of deep neural networks for data with low-dimensional structures

Abstract: In the past decade, deep learning has made astonishing breakthroughs in various real-world applications. It is a common   belief that deep neural networks are good at learning various geometric structures hidden in data sets, such as rich local regularities,   global symmetries, or repetitive patterns. One of the central interests in deep learning theory is to understand why deep neural   networks are successful, and how they utilize low-dimensional data structures. In this talk, I will present some statistical learning   theory for deep ReLU networks where data exhibit low-dimensional structures, such as lying on a low-dimensional manifold. The   learning tasks include regression, classification and learning operators between Hilbert spaces. When data are sampled on a low-   dimensional manifold, the sample complexity crucially depends on the intrinsic dimension of the manifold instead of the ambient   dimension of the data. These results demonstrate that deep neural networks are adaptive to low-dimensional geometric structures of   data sets.