Abstract: The minisymposium seeks to bring together researchers to promote the exchange of ideas, and present recent developments on analysis and numerical methods for PDEs.
Organizers:
Saturday, October 19th, 2019 at 10:20  11:40 (205 Carver Hall)  

10:20  10:40  Cheng Wang, University of Massachusetts Dartmouth 
A Weakly Nonlinear, Energy Stable Scheme for the Strongly Anisotropic CahnHilliard System  
The strongly anisotropic CahnHilliard model is considered. In particular, a highly nonlinear anisotropic surface energy makes the PDE system very challenging at both the analytic and numerical levels. In this talk, a convexity analysis is performed to the surface energy potential, and a careful estimate reveals that all its second order functional derivatives stay uniformly bounded by a global constant. In turn, a linear approximation becomes available for the surface energy part, and a detailed estimate demonstrates the corresponding energy stability. Its combination with the implicit treatment of the nonlinear double well potential terms yields a weakly nonlinear, energy stable scheme for the whole system. Moreover, with a careful application of the global bound for the second order functional derivatives, an optimal rate convergence analysis becomes available, which is the first such result in this area. Some numerical simulation results are also presented in the talk. 

10:40  11:00  Yuan Liu, Wichita State University 
Conservative Total Variation Bounded Spectral Methods for Onedimensional Scalar Conservation Laws  
In this talk, we will introduce a class of total variation bounded (TVB) spectral methods for solving one dimensional scalar hyperbolic conservation laws. The main idea is to apply a TVB limiter to the existing spectral methods such that the total variation bounded. Numerical results are provided to demonstrate the capability of our proposed schemes. 

11:00  11:20  Zhuoran Wang, Colorado State University 
deal.II Implementation of a Twofield Finite Element Solver for Linear Poroelasticity on Quadrilateral Meshes  
This talk presents a finite element solver for linear poroelasticity on quadrilateral meshes based on the twofield model (solid displacement and fluid pressure) and its implementation on the deal.II platform. This solver combines the classical Lagrangian elements with reduced integration technique for the displacement in elasticity, the weak Galerkin elements for the pressure in Darcy flow through the implicit Euler discretization. The solver is penaltyfree and has fewer unknowns compared to other existing methods. Numerical experiments will be presented to show the solver is free of nonphysical pressure oscillations. This is joint work with Dr. James Liu and Dr. Simon Tavener. 

11:20  11:40  Hailiang Liu, Iowa State University 
Selection Dynamics for Deep Neural Networks  
Deep learning is machine learning using neural networks with many hidden layers, and it has become a primary tool in a wide variety of practical learning tasks, such as image classification, speech recognition, driverless cars, or game intelligence. As such, there is a pressing need to provide a solid mathematical framework to analyze various aspects of deep neural networks. This work introduces the mathematical formulation of deep residual neural networks as a PDE optimal control problem. We study the wellposedness, the large time solution behavior, and the characterization of the steady states for the forward. We state and prove optimality conditions for the inverse deep learning problem, using the HamiltonJacobiBellmann equation and the Pontryagin maximum principle. This serves to establish a mathematical foundation for investigating the algorithmic and theoretical connections between optimal control and deep learning. This is joint work with Peter Markowich (KAUST). 
Saturday, October 19th, 2019 at 14:40  16:00 (205 Carver Hall)  

14:40  15:00  Jay Mayfield, Iowa State University 
Asymptotic Methods for Wave Equations in Inhomogeneous Media  
We propose asymptotic methods that combine integral representations of the wavefield and approximations formulated by geometrical optics of the Green’s functions for solving the variable velocity wave equation. The wave is propagated as an integral with Green’s functions whose phase and amplitude terms are described as solutions of an Eikonal equation and a recurrent system of transport equations,respectively. The integral can be efficiently evaluated by the fast Fourier Transform after an appropriate lowrank approximation with Chebyshev polynomial interpolation. To truncate the infinite domain, we apply perfectly matched layers to modify the wave equation. Results and numerical experiments and presented and discussed. 

15:00  15:20  Lei Wang, University of Wisconsin, Milwaukee 
A KernelIndependent Treecode Based on Barycentric Lagrange Interpolation  
A kernelindependent treecode (KITC) is presented for fast summation of pairwise particle interactions. In general, treecodes replace particleparticle interactions by particlecluster interactions, and here we employ barycentric Lagrange interpolation at Chebyshev points to compute wellseparated particlecluster interactions. The scheme requires only kernel evaluations, is suitable for nonoscillatory kernels, and utilizes a scaleinvariance property of barycentric Lagrange interpolation. For a given level of accuracy, the treecode reduces the operation count for pairwise interactions from \(O(N^2)\) to \(O(N \log N)\), where \(N\) is the number of particles in the system. The algorithm is demonstrated in serial and parallel simulations for systems of regularized Stokeslets and rotlets in 3D, and numerical results show the treecode performance in terms of error, CPU time, and memory consumption. The KITC is a relatively simple algorithm with low memory consumption, and this enables a straightforward parallelization. 

15:20  15:40  Peimeng Yin, Wayne State University 
Unconditionally Energy Stable DG Schemes for the SwiftHohenberg Equation  
The SwiftHohenberg equation as a central nonlinear model in modern physics has a gradient flow structure. Here we introduce fully discrete discontinuous Galerkin (DG) schemes for a class of fourth order gradient flow problems, including the nonlinear SwiftHohenberg equation, to produce freeenergydecaying discrete solutions, irrespective of the time step and the mesh size. We exploit and extend the mixed DG method introduced in [H. Liu and P. Yin, J. Sci. Comput., 77: 467501, 2018] for the spatial discretization, and the ``Invariant Energy Quadratization" method for the time discretization. The resulting IEQDG algorithms are linear, thus they can be efficiently solved without resorting to any iteration method. We actually prove that these schemes are unconditionally energy stable. We present several numerical examples that support our theoretical results and illustrate the efficiency, accuracy and energy stability of our new algorithm. The numerical results on two dimensional pattern formation problems indicate that the method is able to deliver comparable patterns of high accuracy. 

15:40  16:00  Ruijun Zhao, Minnesota State University, Mankato 
An Efficient Numerical Method for an SIR Model with Directed Diffusion  
Population spatial movement is often modeled by random diffusion (to capture the intraspecies movement) and directed diffusion (to capture the interspecies movement). In this talk, we will first introduce an SIR model with directed spatial movement, particularly total population moving away from crowd and susceptible population moving away from infected group. The model is a degenerated system of secondorder partial differential equations, which poses challenges in designing efficient numerical methods to solve it. In the second part of this talk, we propose an efficient numerical scheme, the implicit integration factor (IIF) WENO scheme. We will discuss the order of convergence and its performance in solving this system. 
Saturday, October 19th, 2019 at 16:40  18:00 (205 Carver Hall)  

16:40  17:00  Shuwang Li, Illinois Institute of Technology 
Weakly Nonlinear Investigation of the SaffmanTaylor Instability in a Radial Spreading Annulus  
We study viscous fingering formation in an immiscible threelayer HeleShaw problem, where two coupled interfaces spread radially. More specifically, we examine how the initial distance between the inner and outer interfaces (i.e., annulus's thickness) influence the shape of the emerging fingering patterns, mainly on its impact on finger tipbroadening and finger tipsplitting events for both the interfaces. The problem is tackled by employing a perturbative modecoupling approach and we focus our attention on weakly nonlinear stages of the dynamics. By considering different values of viscosity ratios, we perform our analysis in three different scenarios: (a) both interfaces are unstable, (b) inner interface is unstable while the outer is stable, and (c) inner interface is stable and the outer one is unstable. Our weakly nonlinear findings indicate dramatically differences in the dynamics of the coupled interfaces as we explore these three cases of interest. Full nonlinear simulations confirm these findings. 

17:00  17:20  Songting Luo, Iowa State University 
An Effective Approach for High Frequency Helmholtz Equation through the Steadystate of Timedependent Schrodinger Equation  
For computing the solution of the high frequency Helmholtz equation, geometrical optics (GO) approximation provides an efficient way to achieve the goal. However, GO approximation is only valid locally near the primary source and fails in the region where caustics occur. To capture the caustics faithfully, we will compute the solution as the steadystate solution of an appropriate timedependent Schrodinger equation. Numerical experiments will be presented to demonstrate the methods. 

17:20  17:40  Stefan Schnake, University of Oklahoma 
\(H^1\)norm Error Estimate for a Nonstandard Finite Element Approximation of Secondorder Linear Elliptic PDEs in Nondivergence Form  
This talk establishes the optimal \(H^1\)norm error estimate for a nonstandard finite element method for approximating \(H^2\) strong solutions of second order linear elliptic PDEs in nondivergence form with continuous coefficients. To circumvent the difficulty of lacking an effective duality argument for this class of PDEs, a new analysis technique is introduced; the crux of it is to establish an \(H^1\)norm stability estimate for the finite element approximation operator which mimics a similar estimate for the underlying PDE operator recently established by the authors and its proof is based on a freezing coefficient technique and a topological argument. Moreover, both the \(H^1\)norm stability and error estimate also hold for the linear finite element method. 

17:40  18:00  Wei Zhu, University of Alabama, Tuscaloosa 
Fast algorithm for \(L^p\)norm of mean curvature based image denoising models  
In this talk, we will discuss a new class of imaging denoising models by using the \(L^p\)norm of mean curvature of image graphs as regularizers with \( 1 < p \le 2 \). The motivation of introducing such models is to add stronger regularizations than that of the original mean curvature based image denoising model in order to remove noise more eﬃciently. To minimize these variational models, we develop a novel augmented Lagrangian method, and one thus just needs to solve two linear elliptic equations to ﬁnd saddle points of the associated augmented Lagrangian functionals. Speciﬁcally, we linearize the nonlinear term in one of the two subproblems and minimize a proximallike functional that can be easily treated. We prove that the minimizer of the substitute functional does reduce the value of the original subproblem under certain conditions. Numerical results are presented to illustrate the features of the proposed models and also the eﬃciency of the designed algorithm. 
Sunday, October 20th, 2019 at 10:20  11:40 (205 Carver Hall)  

10:20  10:40  Wumaier Maimaitiyiming, Iowa State University 
Positivitypreserving Schemes for a Reduced PoissonNernstPlanck System  
The PoissonNernstPlanck (PNP) system is a widely accepted model for simulation of ionic channels. We design, analyze, and numerically validate a second order unconditional positivitypreserving scheme for solving a reduced PNP system, which can well approximate the three dimensional ion channel problem. Positivity of numerical solutions is proven to hold true independent of the size of time steps and the choice of the Poisson solver. The scheme is easy to implement without resorting to any iteration method. Several numerical examples further confirm the positivitypreserving property, and demonstrate the accuracy, efficiency, and robustness of the proposed scheme. 

10:40  11:00  Yijin Gao, Iowa State University 
Solving Timedependent Schrodinger Equation Numerically with Fast Huygens Sweeping Methods  
Solving timedependent Schrodinger equation numerically is always challenging in that the wave function is oscillatory and the domain is infinite. To deal with the problem, we present Strang operator splitting based fast huygens sweeping method. Also, the perfectly matched layer approach is used to limit the computation onto a bounded subdomain. The method is asymptotic and combines Strang operator splitting technique, the Huygen's principle and the WKBJ approximation. Implementation details will be explained with formulations, and numerical examples will be used to demonstrate the method. 

11:00  11:20  Colin Jensen, Colorado State University 
Two Approaches for Numerical Implementation of Helmholtz Decomposition  
Helmholtz decomposition is a method of expressing a smooth vector field as the sum of a divergencefree vector field and a curlfree vector field. This decomposition has important applications in physics, engineering, and data science. We will present two methods of numerically implementing the Helmholtz decomposition. The crux of both approaches is to numerically solve the Poisson equation with pure Neumann boundary conditions. Our first approach accomplishes this with a straightforward finite difference method. Our second approach uses the novel weak Galerkin finite element method. The two methods will be illustrated by numerical examples and compared briefly. We will also discuss Matlab and Python implementations of these methods. 

Florida International University 

