# FASTMath Seminar Series

### The purpose of the seminar series is to invite domain scientists to talk about their math challenges in order to motivate and encourage collaborations between science teams and FASTMath.

### May 20, 2022, 4:00pm-5:00pm Eastern

Moderators: Yang Liu, LBNL; Ben Whitney, ORNL

Speakers:

Michael Zingale, Stony Brook University: "Algorithmic Improvements for Coupling Hydrodynamics and Reactions in Astrophysical Flows"

*Energy production in stellar environments is dominated by thermonuclear energy release, which can take place in quiet convective flows or explosive environments. Multidimensional models of astrophysical reactive flows have many challenges: capturing the relevant length and timescales, and keeping the different physics inputs coupled together over the course of the simulation. Reaction networks are stiff, requiring implicit time integrators, while the hydrodynamics is treated explicitly. The standard approach in the field is operator splitting of the reactions, but this can break down when the reactions are vigorous. I will discuss some new time integration methods we have been developing in the AMReX-Astrophysics Castro code (https://github.com/amrex-astro/Castro). Two alternate time-integration methods are being developed. The first follows the traditional spectral deferred corrections (SDC) method of using low order approximations to the time updates to build a high-order solution through iteration. With this approach we can do fully fourth-order reactive hydrodynamics simulations. We are also exploring a simpler approach that uses some of the ideas of SDC together with our existing hydrodynamics and reaction solvers to eliminate the operator splitting error. Both methods will be shown and comparisons to the traditional operator splitting approach will be discussed.*Yuxi Chen, Princeton University: "Magnetohydrodynamics with Embedded Particle-in-Cell Model and Its Application to Magnetic Reconnection"

*It is challenging to capture kinetic phenomena in global simulations due to significant differences between the kinetic scales and global scales. The magnetohydrodynamics with embedded particle-in-cell model (MHD-EPIC) is developed to incorporate local kinetic physics into global simulations. MHD-EPIC combines the physics capability of a particle-in-cell (PIC) code and the efficiency of an MHD model by coupling a semi-implicit PIC code with an MHD model. The PIC code is used to cover regions where kinetic effects are important, while the MHD model handles the rest part of the simulation domain. So far, the application of the MHD-EPIC model has been focused on magnetic reconnection, which is a ubiquitous phenomenon in many laboratory and astrophysical systems, including fusion devices, stellar atmospheres, and (exo-)planetary and pulsar magnetospheres. In this presentation, I will first discuss the model coupling approach and then show an application of the MHD-EPIC model to the magnetic island coalescence problem and its comparison with the full PIC simulation.*

### March 8, 2022, 11:00am-12:00pm Eastern

Moderators: Ahmed Attia, ANL; Juliane Mueller, LBNL

Speaker:

James Amundson, FNAL: "Computational Topics in Particle Accelerator Simulation"

*Particle accelerators are enabling technology for scientific research high energy physics, nuclear physics, and many other fields. The most computationally difficult aspect of simulation the dynamics of particle accelerator beams is the inclusion of collective effects. I will give an overview of the challenges presented in the simulation of collective effects in beam dynamics. A major application of beam dynamic simulations is in the design optimization of new and/or upgraded accelerators. I will also describe topics in accelerator design optimization and recent progress in addressing them.*

### February 16, 2022, 3:00pm-4:00pm Eastern

Moderators: Cody Balos, LLNL; Cameron Smith, RPI

Speakers:

David Green, ORNL: "ASCR-Relevant Challenges of the Fusion Energy RF-SciDAC Partnership"

*The RF-SciDAC Center has the goal of predicting how magnetically confined fusion plasmas respond to externally applied Radio Frequency (RF) power in the context of heating and current drive for fusion reactors. Amongst the present and near-term ASCR-relevant challenges the Center is facing are the following topics: (i) solving time-harmonic (indefinite) Maxwell's equations with non-linear boundary conditions in 3D domains with high geometric fidelity boundary shapes - this necessitates degrees-of-freedom reduction via high order (elements and meshing), adaptivity, preconditioning to enable robust scaling of iterative solvers, and simplification of the interaction with CAD representations of domain boundaries; (ii) solving fluid transport in highly anisotropic (magnetized plasma) mediums, again with high geometric fidelity domain boundaries which are not aligned with the anisotropy - our approach here, to enable robust meshing of variable simulation geometries, is a combination of an unstructured mesh with high order elements to resolve the anisotropy without introducing unwanted transport. This approach requires the development of preconditioners which reduce the condition number of the resulting system for high order / high anisotropy cases; (iii) high-dimensional PDEs describing kinetic transport - presently we are investigating adaptive sparse-grids; (iv) non-local kernels for the plasma conductivity in the wave equation, which have previously been handled via Fourier spectral methods but which are insufficient when incorporating high geometric fidelity in the domain boundary; and (v) robust frameworks for the acceleration of code-to-code Picard type iteration. In this presentation, we will give an overview of these challenges, the status of our progress on solving them and identify opportunities for collaboration.*Jeff Candy, General Atomics: "Eulerian Gyrokinetics: Spectral and pseudospectral discretization in CGYRO"

### January 31, 2022, 4:00pm-5:00pm Eastern

Moderators: Yang Liu, LBNL; Ben Whitney, ORNL

Speakers:

Thomas Maier, ORNL: "Math Needs for Materials Quantum Monte Carlo Applications"

*Quantum Monte Carlo (QMC) methods are widely used to study the finite temperature behavior of strongly interacting electron systems. Here I will use the DCA++ application that implements a dynamic cluster QMC algorithm as an example to discuss several challenges that could benefit from collaborations with applied mathematicians. These include acceleration of the simulations (AI/ML), extraction of real time dynamics from the imaginary time QMC data (inverse problems), and managing large tensors that describe the effective interaction between electrons and provide the deepest insight into the physics of the system (linear eigensolvers and data compression).*Martin Head-Gordon, UC Berkeley/LBNL: "Can applied math help solve the electron correlation problem in computational quantum chemistry?"

*This talk will provide a discussion of some issues that the applied mathematics community may be able to address as part of the fundamental challenge of efficiently evaluating the correlation energy of electrons in molecules. Computational quantum chemistry seeks to approximate the many-electron Schrodinger equation by introducing (i) the algebraic approximation (a one-particle basis set to convert a partial differential equation into algebraic equations), and (ii) the correlation model (which truncates those equations in a rational and efficient way). These approximations together reduce the formal exponential compute complexity of the Schrodinger equation to polynomial. However, the development of fast algorithms for electron correlation offers the prospect of further lowering the non-linear polynomial scaling of compute costs, and is an area where applied math is poised to play an important role in our SciDAC scientific partnership. I will present some of the challenges and the possibilities using one of the simplest correlation models, second order Moller-Plesset (MP2) theory, which has compute costs scaling formally with the 5th power of molecule size.*

### October 26, 2021, 4:00pm-5:00pm Eastern

Moderators: Cody Balos, LLNL; Roel Van Beeumen, LBNL

Speakers:

Eirik Endeve, ORNL: "Considerations for time integration methods in nuclear astrophysics applications"

*Core-collapse supernovae and binary neutron star mergers are astrophysical events of considerable interest in nuclear physics. They are dominant sources of heavy elements in the Universe, and emit photons, neutrinos, and gravitational waves. Gaining insights into the physical processes driving these multi-messenger events relies heavily on modeling and simulation. The models solve a coupled system of equations for self-gravity, magneto-hydrodynamics for nuclear matter, nuclear reaction networks, and neutrino transport, and involve a wide range of spatial and temporal scales. In this talk, we will give a brief overview of the physical models and some of the numerical methods targeted by the Exascale Computing Project’s ExaStar team. To initiate discussions, we will place particular emphasis on considerations relevant to the design of time integration methods for deployment in nuclear astrophysics applications.*Robert Edwards, Jefferson Lab: "Computing the Properties of Matter"

*I will review some of the numerical challenges faced in using lattice field theory techniques for computations in Nuclear Physics.*

### September 14, 2021, 4:00pm-5:00pm Eastern

Moderators: Cameron Smith, RPI; Ben Whitney, ORNL

Speakers:

Stephen Price, LANL: "Future numerical and computational challenges in DOE ice sheet modeling"

*In this talk, I will briefly summarize DOE progress on ice sheet model development and applications during the past 10 years, under joint BER and ASCR funding. I will then discuss remaining challenges and improvements that future FASTMath research efforts might contribute to, with a focus on four areas: 1) numerical and computational methods, 2) meshing, 3) optimization and uncertainty quantification approaches, and 4) model physics. Specific context will be provided through examples of past and ongoing efforts in the area of ice sheet modeling.*Peter Bosler, SNL: "Compact, Performance-Portable Semi-Lagrangian Methods for E3SM"

*The ASCR/BER partnership SciDAC Project, \emph{Non-hydrostatic dynamics with characteristic discontinuous Galerkin methods}, develops semi-Lagrangian (SL) algorithms and associated software for passive tracer transport in E3SM that are both (a) tailored for the advanced architectures of current and anticipated DOE LCF computing platforms and (b) provably successful at achieving the fundamental required traits of a climate model's transport scheme: conservation, accuracy, tracer consistency, shape preservation, and computational efficiency. Integrating SL transport into E3SM frequently requires algorithmic improvements in other parts of the code (time stepping, in particular) and exposes opportunities for investigating non-hydrostatic effects and radiative-convective equilibrium at high resolution. In this talk, we examine the impacts and applications of this work on the E3SM Atmosphere Model's (EAM) version 2 and discuss ongoing work with MPAS-Ocean targeted at the version 3 biogeochemistry science campaign. We conclude with discussion questions associated with this work that may be relevant for future FastMath collaborations.*

### July 29, 2021, 4:00pm-5:00pm Eastern

Moderators: Ahmed Attia, ANL; Yang Liu, LBNL

Speakers:

S.C. Jardin, PPPL: "The M3D-C1 code"

*We describe the M3D-C1 code, emphasizing the linear solvers and opportunities for uncertainty quantification. M3D-C1 is an initial-value magnetohydrodynamic code based on 3D finite elements with C1 continuity and semi-implicit time-stepping. The finite element mesh is unstructured in the (R,Z) plane, but structured in the toroidal angle φ. We take advantage of the structure by utilizing a block-Jacobi preconditioner based on multiple instances of SuperLU_dist for the full 3D GMRES iterative sparse matrix solves. The solution depends in a complex way on the initial conditions, boundary conditions, and the values of the dissipation coefficients (particle diffusivity, thermal conductivity, viscosity, electrical resistivity) which also determine resolution requirements.*CS Chang, PPPL: "Needs for Advanced Linear Systems and UQ in the Nonequilibrium Fusion Kinetic Code XGC"

*XGC started as a partnership code between fusion application SciDAC Center and the ASCR SciDAC Institutes a decade and a half ago to solve the nonequilibrium kinetic physics problems in the boundary region of magnetic fusion reactors utilizing particle-in-cell method on leadership class computers. XGC has been making difficult scientific discoveries that have not been possible by other methods. XGC is now entering into new multiscale interaction regimes that solve together the small-scale/small-amplitude kinetic turbulence, the large-scale/large-amplitude fluid instabilities, and the background plasma profile dynamics. The soon-to-arrive exascale computers will be powerful enough to enable such a multiscale simulation and prediction of ITER plasmas in realistic geometry. New physics capabilities may require new mathematical tools that can optimize XGC’s performance on the new exascale hardware/software architectures. Uncertainty quantification of the exascale computing results and reduction to digital-twins/surrogate-models are another serious class of topics that XGC is facing. In this discussion, we will emphasize the needs for advanced linear and UQ systems. One advantageous UQ property of XGC is that, even though the simulations are expensive, the input dimensionality is low (only on the order of a dozen) and a single simulation can provide many data points since the background plasma itself is time-evolving. Successful digital-twins and surrogate-models can be highly valuable in experimental design.*