Applied/ACMS: Difference between revisions

From UW-Math Wiki
Jump to navigation Jump to search
 
(72 intermediate revisions by 6 users not shown)
Line 10: Line 10:
<br>   
<br>   


== '''Fall 2024''' ==
== '''Spring 2025''' ==
{| class="wikitable"
{| cellpadding="8"
|+
! align="left" |Date
!Date
! align="left" |Speaker
!Speaker
! align="left" |Title
!Title
! align="left" |Host(s)
!Host(s)
|-
|-
|Sep 13*
|Jan 31
|[https://people.math.wisc.edu/~nchen29/ Nan Chen] (UW)
|[https://people.math.wisc.edu/~tgchandler/ Thomas Chandler] (UW)
|Intro. to Uncertainty Quantification (UQ) (tutorial)
|[[#Chandler|''Fluid–structure interactions in active complex fluids'']]
|Spagnolie
|Spagnolie
|-
|-
|Sep 20
|Feb 7
|[https://knewhall.web.unc.edu Katie Newhall] (UNC Chapel Hill)
|[https://afraser3.github.io/ Adrian Fraser] (Colorado)
|Energy landscapes, metastability, and transition paths
|[[#Fraser|''Destabilization of transverse waves by periodic shear flows'']]
|Rycroft
|Spagnolie
|-
|-
|Sep 27
|Feb 14
|[https://ptg.ukzn.ac.za Indresan Govender] (Mintek / Univ. of KwaZulu-Natal, South Africa)
|[https://jrluedtke.github.io/ Jim Luedtke] (UW)
|Granular flow modeling and visualization using nuclear imaging
|[[#Luedtke|Using integer programming for verification of binarized neural networks]]
|Rycroft
|Spagnolie
|-
|-
|Oct 4*
|Feb 21
|[https://sse.tulane.edu/math/people/hongfei-chen Hongfei Chen] (Tulane)
|[https://zhdankin.physics.wisc.edu/ Vladimir Zhdankin] (UW)
|
|[[#Zhdankin|Exploring astrophysical plasma turbulence with particle-in-cell methods]]
|Jean-Luc
|Spagnolie
|-
|-
|Oct 11 '''Colloquium in B239 at 4:00pm'''
|Feb 28
|[https://people.math.ethz.ch/~imikaela/ Mikaela Iacobelli] (ETH/IAS)
|[https://nmboffi.github.io/ Nick Boffi] (CMU)
|[[# TBA| TBA  ]]
|[[#Boffi|Generative modeling with stochastic interpolants]]
|Li
|Li, Rycroft
|-
|-
|Oct 18 '''Colloquium in B239 at 4:00pm'''
|Mar 7
|[https://galton.uchicago.edu/~guillaumebal/ Guillaume Bal] (U Chicago)
|[https://sites.lsa.umich.edu/shankar-lab/ Suraj Shankar] (Michigan)
|[[# TBA| TBA  ]]
|[[#Shankar|Designer active matter]]
| Li, Stechmann
|Spagnolie
|-
|-
|Oct 23 ('''Wednesday''')
|Mar 10
|[https://www.sandia.gov/ccr/staff/teresa-portone/ Teresa Portone] (Sandia)
|[https://www.math.kit.edu/csmm/~loevbak/en Emil Loevbak] (KIT)
|TBD
|[[#Loevbak|Discrete adjoint Monte Carlo for kinetic equations with reversible pseudorandom generators]]
|Stechmann
|Li
|-
|-
|Oct 25
|Mar 14
|[https://www.cs.cornell.edu/~damle/ Anil Damle] (Cornell)
|[https://lu.seas.harvard.edu/ Yue Lu] (Harvard) '''[Colloquium]'''
|Fine-grained Theory and Hybrid Algorithms for Randomized Numerical Linear Algebra
|[[#Lu|Nonlinear Random Matrices in Estimation and Learning: Equivalence Principles and Applications]]
|Li
|Li
|-
|-
| Nov 1
|Mar 21
|[https://research-hub.nrel.gov/en/persons/michael-sprague Michael Sprague] (NREL)
|[https://people.llnl.gov/vogman1 Genia Vogman] (LLNL)
|[[# TBA| TBA ]]
|[[#Vogman|TBA]]
|Spagnolie
|Li
|-
|-
| Nov 8
|Mar 28
|[https://personal.math.ubc.ca/~holmescerfon/ Miranda Holmes-Cerfon] (UBC)
|''Spring Break''
|
|
|
|Stechmann
|-
|-
| Nov 15*
|Apr 4
| [http://sun-yue.com Yue Sun] (UW–Madison)
|[https://mathsci.kaist.ac.kr/~donghwankim/ Donghwan Kim] (KAIST)
|
|TBA
| Rycroft
|Lyu
|-
|-
| Nov 22
|Apr 11
|[https://ibd.uchicago.edu/joinus/yenfellowship/ Ondrej Maxian] (U Chicago)
|[https://meche.mit.edu/people/faculty/pierrel@mit.edu Pierre Lermusiaux] (MIT)
|[[# TBA| TBA ]]
|[[#Lermusiaux|TBA]]
|Ohm & Spagnolie
|Chen
|-
|-
| Nov 29*
|Apr 18
|''Thanksgiving''
|[https://www.math.uci.edu/~jxin/ Jack Xin] (UC Irvine) '''[Colloquium]'''
|
|[[#Xin|TBA]]
|
|
|-
|-
| Dec 6
|Apr 25
|[https://www.simonsfoundation.org/people/ido-lavi/ Ido Lavi] (Flatiron)
|[https://www-users.cse.umn.edu/~bcockbur/ Bernardo Cockburn] (Minnesota)
|[[# TBA| TBA ]]
|[[#Cockburn|''Transforming stabilization into spaces'']]
|Spagnolie
| Stechmann, Fabien
|-
|May 2
|[https://sylviaherbert.com/ Sylvia Herbert] (UCSD)
|[[#Herbert|TBA]]
|Chen
|}
|}


Dates marked with an asterisk correspond to [https://uwbadgers.com/sports/football/schedule home football games of the UW–Madison Badgers]. On these dates it can be difficult to get a hotel room close to campus at short notice.
==Abstracts==
 
<div id="Chandler">
====Thomas G. J. Chandler (UW)====
Title: Fluid-structure interactions in active complex fluids
 
Fluid anisotropy is central to many biological systems, from rod-like bacteria that self-assemble into dense swarms that function as fluids, to the cell cytoskeleton where the active alignment of stiff biofilaments is crucial to cell division. Nematic liquid crystals provide a powerful model for studying these complex environments. However, large immersed bodies elastically frustrate these fluids, leading to intricate interactions. This frustration can be alleviated through body deformations, at the cost of introducing internal stresses. Additionally, active stresses, arising from particle motility or molecular activity, disrupt nematic order by driving flows. In this presentation, I will demonstrate how complex variables enable analytical solutions to a broad range of problems, offering key insights into the roles of body geometry, anchoring conditions, interaction dynamics, activity-induced flows, and body deformations in many biological settings.


== Abstracts ==
<div id="Fraser">
====Nan Chen (UW–Madison)====
====Adrian Fraser (Colorado)====
Title: Taming Uncertainty in a Complex World: The Rise of Uncertainty Quantification -- A Tutorial for Beginners
Title: Destabilization of transverse waves by periodic shear flows


I will provide a tutorial about uncertainty quantification (UQ) for those who have no background but are interested in learning more about this area. The talk will exploit many elementary examples, which are understandable to graduate students and senior undergraduates, to present the ideas of UQ. Topics include characterizing uncertainties using information theory, UQ in linear and nonlinear dynamical systems, UQ via data assimilation, the role of uncertainty in diagnostics, and UQ in advancing efficient modeling. The surprisingly simple examples in each topic explain why and how UQ is essential. Both Matlab and Python codes have been made available for these simple examples.
Periodic shear flows have the peculiar property that they are unstable to large-scale, transverse perturbations, and that this instability proceeds via a negative-eddy-viscosity mechanism (Dubrulle & Frisch, 1991). In this talk, I will show an example where this property causes transverse waves to become linearly unstable: a sinusoidal shear flow in the presence of a uniform, streamwise magnetic field in the framework of incompressible MHD. This flow is unstable to a KH-like instability for sufficiently weak magnetic fields, and uniform magnetic fields permit transverse waves known as Alfvén waves. Under the right conditions, these Alfvén waves become unstable, presenting a separate branch of instability that persists for arbitrarily strong magnetic fields which otherwise suppress the KH-like instability. After characterizing these waves with the help of a simple asymptotic expansion, I will show that they drive soliton-like waves in nonlinear simulations. With time permitting, I will discuss other fluid systems where similar dynamics are or may be found, including stratified flows and plasma drift waves.


'''Katie Newhall (UNC Chapel Hill)'''
<div id="Luedtke">
====Jim Luedtke (UW)====
Title: Using integer programming for verification of binarized neural networks


Title: Energy landscapes, metastability, and transition paths
Binarized neural networks (BNNs) are neural networks in which the weights are binary and the activation functions are the sign function. Verification of BNNs against input perturbation is one way to measure robustness of BNNs. BNN verification can be formulated as an integer linear optimization problem and hence can in theory be solved by state-of-the art methods for integer programming such as the branch-and-cut algorithm implemented in solvers like Gurobi. Unfortunately, the natural formulation is often difficult to solve in practice, even by the best such solvers, due to large integrality gap induced by its so-called "big-M" constraints. We present simple but effective techniques for improving the ability of the integer programming approach to solve the verification problem for BNNs. Along the way, we hope to illustrate more generally some of the strategies integer programmers use to attack difficult problems like this. We find that our techniques enable verifying BNNs against a higher range of input perturbation than using the natural formulation directly.


The concept of an energy landscape emerged in the 1930’s as a way to calculate chemical reaction rate constants via Henry Eyring’s transition state theory. Its use has expanded since then, remaining central to quantifying metastability (infrequent jumps between deterministically-stable, energy minimizing, states) that arises in noisy systems when the thermal energy is small relative to the energy barrier separating two states. In this talk, I will present extensions of this theory that I have developed and applied to physical and biological systems. The first is an infinite dimensional system for which I prove metastability is present in the absence of an energy barrier; I extend transition state theory to compute mean transition times. In the second, I derive a model for a spatially-extended magnetic system with spatially-correlated noise designed to sample the Gibbs distribution relative to a defined energy functional. In the third, I show a quasi-potential can be found and used to describe metastable transitions between stable clusters in a bead-spring polymer model of chromosome dynamics with additional stochastic binding pushing the system out of equilibrium.
This is joint work with Woojin Kim, Mathematics PhD student at UW-Madison.


<div id="Zhdankin">
====Vladimir Zhdankin (UW)====
Title: Exploring astrophysical plasma turbulence with particle-in-cell methods


'''Indresan Govender (Mintek / Univ. of KwaZulu Natal, South Africa)'''
Plasmas throughout the universe (as well as in the laboratory) tend to exist in turbulent, nonequilibrium states due to their "collisionless" nature. Described by the Vlasov-Maxwell equations in a six-dimensional phase space (of position and momentum), the basic physics of such plasmas is difficult to model from first principles. There remain open questions about entropy production, nonthermal particle acceleration, energy partition amongst different particle species, and more. Particle-in-cell simulations are a numerical tool that allow us to explore in depth the rich dynamics and statistical mechanics of collisionless plasmas, validating analytical speculation. I will describe some of the results from my group's work on this topic.


Title: Granular flow modeling and visualization using nuclear imaging
<div id="Boffi">
====Nick Boffi (CMU)====
Title: Generative modeling with stochastic interpolants


Despite its ubiquity, a complete theory to describe the underlying rheology of granular flows remains elusive. Central to this problem is the lack of detailed, in-situ measurements of the granular flow field. To this end, we present two non-invasive imaging techniques currently employed to measure the flow of individual grains within granular flow systems that span simple mono-sized flows of plastic beads to complex industrial mixture flows of rocks and slurry. The first technique employs diagnostic X-rays operated in biplanar mode to triangulate the motion of low-density granules in simplified flow systems to within a 3D spatial accuracy of 0.15 mm at tracking frequencies up to 100 Hz. The second—arguably the workhorse of our research operation—is the nuclear imaging technique of Positron Emission Particle Tracking (PEPT) which triangulates the back-to-back gamma rays emanating from radiolabeled particles to within a millimeter in 3D space at a millisecond timing resolution. PEPT can track the motion of any particle with a diameter greater than ∼20 microns. Both techniques are well suited to studying the flow of granular materials after the data is cast into volume and time averages consistent with the continuum framework. In this talk I will explore the many interesting analysis techniques employed to mapping out the complex flow regimes found in typical granular systems, and the insights they offer towards better understanding their rheological character. Examples explored will include rotating drum flows (wet and dry), shear cells and their industrial counterpart the IsaMill<sup>TM</sup>, hydrocyclone separator flows, and the motivation for tracking of multiple particles. The validation offered to numerical schemes like the Discrete Element Method will also be explored wherein we highlight the complimentary role that measurement and simulation play in unravelling the secrets of granular flows. Finally, and deviating somewhat from the imaging world, I will present our efforts towards utilizing granular flow modeling in real-time control of complex industrial flows encountered in mineral processing.
We introduce a class of generative models that unifies flows and diffusions. These models are built using a continuous-time stochastic process called a stochastic interpolant, which exactly connects two arbitrary probability densities in finite time. We show that the time-dependent density of the stochastic interpolant satisfies both a first-order transport equation and an infinite family of forward and backward Fokker-Planck equations with tunable diffusion coefficients. This viewpoint yields deterministic and stochastic generative models built dynamically from an ordinary or stochastic differential equation with an adjustable noise level. To formulate a practical algorithm, we discuss how the resulting drift functions can be characterized variationally and learned efficiently over flexible parametric classes such as neural networks. Empirically, we highlight the advantages of our formalism -- and the tradeoffs between deterministic and stochastic sampling -- through numerical examples in image generation, inverse imaging, probabilistic forecasting, and accelerated sampling.


==== Anil Damle (Cornell)====
<div id="Shankar">
Title: Fine-grained Theory and Hybrid Algorithms for Randomized Numerical Linear Algebra
====Suraj Shankar (Michigan)====
Title: Designer active matter


Randomized algorithms have gained increased prominence within numerical linear algebra and they play a key role in an ever-expanding range of problems driven by a breadth of scientific applications. In this talk we will explore two aspects of randomized algorithms by (1) providing experiments and accompanying theoretical analysis that demonstrate how their performance depends on matrix structures beyond singular values (such as coherence of singular subspaces), and (2) showing how to leverage those insights to build hybrid algorithms that blend favorable aspects of deterministic and randomized methods. A focus of this talk will be on methods that approximate matrices using subsets of columns. Relevant motivating applications will be discussed and numerical experiments will illuminate directions for further research.
Active matter, i.e., internally driven matter fueled by a sustained dissipation of free energy, is ubiquitous in the natural world. Examples range from bird flocks and human crowds to migrating cells and biopolymer gels, including synthetic systems like phoretic colloids and robots. While much is known about the emergent collective phenomena and complex dynamics that active matter exhibits, little is known about the inverse problem on how they can be controlled. I will discuss a few different vignettes on our recent efforts in controlling flows, forces and physical features of active materials, highlighting implications for the design of novel metamaterials and biomimetic constructs.


== Future semesters ==
<div id="Loevbak">
==== Emil Loevbak (KIT) ====
Title: Discrete adjoint Monte Carlo for kinetic equations with reversible pseudorandom generators


*[[Applied/ACMS/Fall2024|Fall 2024]]
Abstract: Kinetic equations, PDEs modeling particles in a position-velocity phase space, have many high-impact application areas, including nuclear fusion research and radiation therapy. In these applications, one often uses particle-based Monte Carlo methods to simulate the kinetic models. These methods solve the PDE by tracing sample particle trajectories through physical space in such a way that their ensemble distribution in phase-space corresponds with the solution of the PDE. One then uses these trajectories as samples to compute quantities such as the particles' mass-density, momentum, and energy as a function of space and time. These methods have the advantage of not constructing grids in the high-dimensional phase space but the drawback of producing computational results subject to a stochastic sampling error.
 
In this talk I consider PDE-constrained optimization, where a PDE is simulated with a Monte Carlo solver. Here, we compute gradients through a discrete adjoint approach. To ensure optimization convergence, it is imperative to ensure that the same particle trajectories are used when solving the original PDE when evaluating the objective functional and the adjoint PDE when computing gradients. I present an approach of using reversible random number generators to ensure path consistency, despite the adjoint PDE running backward in time. I first present this strategy using a didactic example using a 1D diffusion equation and then present some results from a fusion plasma-edge simulation case.
 
 
<div id="Lu">
==== Yue M. Lu (Harvard) ====
Title: Nonlinear Random Matrices in Estimation and Learning: Equivalence Principles and Applications
 
Abstract: In recent years, new classes of structured random matrices have emerged in statistical estimation and machine learning. Understanding their spectral properties has become increasingly important, as these matrices are closely linked to key quantities such as the training and generalization performance of large neural networks and the fundamental limits of high-dimensional signal recovery. Unlike classical random matrix ensembles, these new matrices often involve nonlinear transformations, introducing additional structural dependencies that pose challenges for traditional analysis techniques.
 
In this talk, I will present a set of equivalence principles that establish asymptotic connections between various nonlinear random matrix ensembles and simpler linear models that are more tractable for analysis. I will then demonstrate how these principles can be applied to characterize the performance of kernel methods and random feature models across different scaling regimes and to provide insights into the in-context learning capabilities of attention-based Transformer networks.
 
Bio: Yue M. Lu is a Harvard College Professor and Gordon McKay Professor of Electrical Engineering and Applied Mathematics at Harvard University. He has also held visiting appointments at Duke University (2016) and the École Normale Supérieure (ENS) in Paris (2019). His research focuses on the mathematical foundations of high-dimensional statistical estimation and learning. His contributions have been recognized with several best paper awards (IEEE ICIP, ICASSP, and GlobalSIP), the ECE Illinois Young Alumni Achievement Award (2015), and the IEEE Signal Processing Society Distinguished Lecturership (2022). He is a Fellow of the IEEE (Class of 2024).
 
 
<div id="Cockburn">
====Bernardo Cockburn (Minnesota)====
Title: Transforming stabilization into spaces


*[[Applied/ACMS/Spring2025|Spring 2025]]
In the framework of finite element methods for ordinary differential equations, we consider the continuous Galerkin method (introduced in 72) and the discontinuous Galerkin method (introduced in 73/74). We uncover the fact that both methods discretize the time derivative in exactly the same form, and discuss a few of its consequences. We end by briefly describing our ongoing work on the extension of this result to some Galerkin methods for partial differential equations.


== Archived semesters ==
== Archived semesters ==


*[[Applied/ACMS/Fall2024|Fall 2024]]
*[[Applied/ACMS/Spring2024|Spring 2024]]
*[[Applied/ACMS/Spring2024|Spring 2024]]
*[[Applied/ACMS/Fall2023|Fall 2023]]
*[[Applied/ACMS/Fall2023|Fall 2023]]

Latest revision as of 22:28, 7 March 2025


Applied and Computational Mathematics Seminar


Spring 2025

Date Speaker Title Host(s)
Jan 31 Thomas Chandler (UW) Fluid–structure interactions in active complex fluids Spagnolie
Feb 7 Adrian Fraser (Colorado) Destabilization of transverse waves by periodic shear flows Spagnolie
Feb 14 Jim Luedtke (UW) Using integer programming for verification of binarized neural networks Spagnolie
Feb 21 Vladimir Zhdankin (UW) Exploring astrophysical plasma turbulence with particle-in-cell methods Spagnolie
Feb 28 Nick Boffi (CMU) Generative modeling with stochastic interpolants Li, Rycroft
Mar 7 Suraj Shankar (Michigan) Designer active matter Spagnolie
Mar 10 Emil Loevbak (KIT) Discrete adjoint Monte Carlo for kinetic equations with reversible pseudorandom generators Li
Mar 14 Yue Lu (Harvard) [Colloquium] Nonlinear Random Matrices in Estimation and Learning: Equivalence Principles and Applications Li
Mar 21 Genia Vogman (LLNL) TBA Li
Mar 28 Spring Break
Apr 4 Donghwan Kim (KAIST) TBA Lyu
Apr 11 Pierre Lermusiaux (MIT) TBA Chen
Apr 18 Jack Xin (UC Irvine) [Colloquium] TBA
Apr 25 Bernardo Cockburn (Minnesota) Transforming stabilization into spaces Stechmann, Fabien
May 2 Sylvia Herbert (UCSD) TBA Chen

Abstracts

Thomas G. J. Chandler (UW)

Title: Fluid-structure interactions in active complex fluids

Fluid anisotropy is central to many biological systems, from rod-like bacteria that self-assemble into dense swarms that function as fluids, to the cell cytoskeleton where the active alignment of stiff biofilaments is crucial to cell division. Nematic liquid crystals provide a powerful model for studying these complex environments. However, large immersed bodies elastically frustrate these fluids, leading to intricate interactions. This frustration can be alleviated through body deformations, at the cost of introducing internal stresses. Additionally, active stresses, arising from particle motility or molecular activity, disrupt nematic order by driving flows. In this presentation, I will demonstrate how complex variables enable analytical solutions to a broad range of problems, offering key insights into the roles of body geometry, anchoring conditions, interaction dynamics, activity-induced flows, and body deformations in many biological settings.

Adrian Fraser (Colorado)

Title: Destabilization of transverse waves by periodic shear flows

Periodic shear flows have the peculiar property that they are unstable to large-scale, transverse perturbations, and that this instability proceeds via a negative-eddy-viscosity mechanism (Dubrulle & Frisch, 1991). In this talk, I will show an example where this property causes transverse waves to become linearly unstable: a sinusoidal shear flow in the presence of a uniform, streamwise magnetic field in the framework of incompressible MHD. This flow is unstable to a KH-like instability for sufficiently weak magnetic fields, and uniform magnetic fields permit transverse waves known as Alfvén waves. Under the right conditions, these Alfvén waves become unstable, presenting a separate branch of instability that persists for arbitrarily strong magnetic fields which otherwise suppress the KH-like instability. After characterizing these waves with the help of a simple asymptotic expansion, I will show that they drive soliton-like waves in nonlinear simulations. With time permitting, I will discuss other fluid systems where similar dynamics are or may be found, including stratified flows and plasma drift waves.

Jim Luedtke (UW)

Title: Using integer programming for verification of binarized neural networks

Binarized neural networks (BNNs) are neural networks in which the weights are binary and the activation functions are the sign function. Verification of BNNs against input perturbation is one way to measure robustness of BNNs. BNN verification can be formulated as an integer linear optimization problem and hence can in theory be solved by state-of-the art methods for integer programming such as the branch-and-cut algorithm implemented in solvers like Gurobi. Unfortunately, the natural formulation is often difficult to solve in practice, even by the best such solvers, due to large integrality gap induced by its so-called "big-M" constraints. We present simple but effective techniques for improving the ability of the integer programming approach to solve the verification problem for BNNs. Along the way, we hope to illustrate more generally some of the strategies integer programmers use to attack difficult problems like this. We find that our techniques enable verifying BNNs against a higher range of input perturbation than using the natural formulation directly.

This is joint work with Woojin Kim, Mathematics PhD student at UW-Madison.

Vladimir Zhdankin (UW)

Title: Exploring astrophysical plasma turbulence with particle-in-cell methods

Plasmas throughout the universe (as well as in the laboratory) tend to exist in turbulent, nonequilibrium states due to their "collisionless" nature. Described by the Vlasov-Maxwell equations in a six-dimensional phase space (of position and momentum), the basic physics of such plasmas is difficult to model from first principles. There remain open questions about entropy production, nonthermal particle acceleration, energy partition amongst different particle species, and more. Particle-in-cell simulations are a numerical tool that allow us to explore in depth the rich dynamics and statistical mechanics of collisionless plasmas, validating analytical speculation. I will describe some of the results from my group's work on this topic.

Nick Boffi (CMU)

Title: Generative modeling with stochastic interpolants

We introduce a class of generative models that unifies flows and diffusions. These models are built using a continuous-time stochastic process called a stochastic interpolant, which exactly connects two arbitrary probability densities in finite time. We show that the time-dependent density of the stochastic interpolant satisfies both a first-order transport equation and an infinite family of forward and backward Fokker-Planck equations with tunable diffusion coefficients. This viewpoint yields deterministic and stochastic generative models built dynamically from an ordinary or stochastic differential equation with an adjustable noise level. To formulate a practical algorithm, we discuss how the resulting drift functions can be characterized variationally and learned efficiently over flexible parametric classes such as neural networks. Empirically, we highlight the advantages of our formalism -- and the tradeoffs between deterministic and stochastic sampling -- through numerical examples in image generation, inverse imaging, probabilistic forecasting, and accelerated sampling.

Suraj Shankar (Michigan)

Title: Designer active matter

Active matter, i.e., internally driven matter fueled by a sustained dissipation of free energy, is ubiquitous in the natural world. Examples range from bird flocks and human crowds to migrating cells and biopolymer gels, including synthetic systems like phoretic colloids and robots. While much is known about the emergent collective phenomena and complex dynamics that active matter exhibits, little is known about the inverse problem on how they can be controlled. I will discuss a few different vignettes on our recent efforts in controlling flows, forces and physical features of active materials, highlighting implications for the design of novel metamaterials and biomimetic constructs.

Emil Loevbak (KIT)

Title: Discrete adjoint Monte Carlo for kinetic equations with reversible pseudorandom generators

Abstract: Kinetic equations, PDEs modeling particles in a position-velocity phase space, have many high-impact application areas, including nuclear fusion research and radiation therapy. In these applications, one often uses particle-based Monte Carlo methods to simulate the kinetic models. These methods solve the PDE by tracing sample particle trajectories through physical space in such a way that their ensemble distribution in phase-space corresponds with the solution of the PDE. One then uses these trajectories as samples to compute quantities such as the particles' mass-density, momentum, and energy as a function of space and time. These methods have the advantage of not constructing grids in the high-dimensional phase space but the drawback of producing computational results subject to a stochastic sampling error.

In this talk I consider PDE-constrained optimization, where a PDE is simulated with a Monte Carlo solver. Here, we compute gradients through a discrete adjoint approach. To ensure optimization convergence, it is imperative to ensure that the same particle trajectories are used when solving the original PDE when evaluating the objective functional and the adjoint PDE when computing gradients. I present an approach of using reversible random number generators to ensure path consistency, despite the adjoint PDE running backward in time. I first present this strategy using a didactic example using a 1D diffusion equation and then present some results from a fusion plasma-edge simulation case.


Yue M. Lu (Harvard)

Title: Nonlinear Random Matrices in Estimation and Learning: Equivalence Principles and Applications

Abstract: In recent years, new classes of structured random matrices have emerged in statistical estimation and machine learning. Understanding their spectral properties has become increasingly important, as these matrices are closely linked to key quantities such as the training and generalization performance of large neural networks and the fundamental limits of high-dimensional signal recovery. Unlike classical random matrix ensembles, these new matrices often involve nonlinear transformations, introducing additional structural dependencies that pose challenges for traditional analysis techniques.

In this talk, I will present a set of equivalence principles that establish asymptotic connections between various nonlinear random matrix ensembles and simpler linear models that are more tractable for analysis. I will then demonstrate how these principles can be applied to characterize the performance of kernel methods and random feature models across different scaling regimes and to provide insights into the in-context learning capabilities of attention-based Transformer networks.

Bio: Yue M. Lu is a Harvard College Professor and Gordon McKay Professor of Electrical Engineering and Applied Mathematics at Harvard University. He has also held visiting appointments at Duke University (2016) and the École Normale Supérieure (ENS) in Paris (2019). His research focuses on the mathematical foundations of high-dimensional statistical estimation and learning. His contributions have been recognized with several best paper awards (IEEE ICIP, ICASSP, and GlobalSIP), the ECE Illinois Young Alumni Achievement Award (2015), and the IEEE Signal Processing Society Distinguished Lecturership (2022). He is a Fellow of the IEEE (Class of 2024).


Bernardo Cockburn (Minnesota)

Title: Transforming stabilization into spaces

In the framework of finite element methods for ordinary differential equations, we consider the continuous Galerkin method (introduced in 72) and the discontinuous Galerkin method (introduced in 73/74). We uncover the fact that both methods discretize the time derivative in exactly the same form, and discuss a few of its consequences. We end by briefly describing our ongoing work on the extension of this result to some Galerkin methods for partial differential equations.

Archived semesters



Return to the Applied Mathematics Group Page