Applied/ACMS: Difference between revisions
Line 5: | Line 5: | ||
*'''When:''' Fridays at 2:25pm (except as otherwise indicated) | *'''When:''' Fridays at 2:25pm (except as otherwise indicated) | ||
*'''Where:''' 901 Van Vleck Hall | *'''Where:''' 901 Van Vleck Hall | ||
*'''Organizers:''' [ | *'''Organizers:''' [https://math.wisc.edu/staff/fabien-maurice/ Maurice Fabien], [https://people.math.wisc.edu/~rycroft/ Chris Rycroft], and [https://www.math.wisc.edu/~spagnolie/ Saverio Spagnolie], | ||
*'''To join the ACMS mailing list:''' Send mail to [mailto:acms+join@g-groups.wisc.edu acms+subscribe@g-groups.wisc.edu]. | |||
<br> | <br> | ||
== Fall | == '''Fall 2024''' == | ||
{| class="wikitable" | |||
{| | |+ | ||
! | !Date | ||
! | !Speaker | ||
! | !Title | ||
! | !Host(s) | ||
|- | |- | ||
| | |Sep 13* | ||
|[ | |[https://people.math.wisc.edu/~nchen29/ Nan Chen] (UW) | ||
| | |Intro. to Uncertainty Quantification (UQ) (tutorial) | ||
|Spagnolie | |||
| | |||
|- | |- | ||
| | |Sep 20 | ||
|[ | |[https://knewhall.web.unc.edu Katie Newhall] (UNC Chapel Hill) | ||
| | |Energy landscapes, metastability, and transition paths | ||
| | |Rycroft | ||
|- | |- | ||
| | |Sep 27 | ||
|[ | |[https://ptg.ukzn.ac.za Indresan Govender] (Mintek / Univ. of KwaZulu-Natal, South Africa) | ||
| | |Granular flow modeling and visualization using nuclear imaging | ||
| | |Rycroft | ||
|- | |- | ||
| | |Oct 4* | ||
|[ | |[https://sse.tulane.edu/math/people/hongfei-chen Hongfei Chen] (Tulane) | ||
| | |Investigating Hydrodynamics of Choanoflagellate Colonies: A Reduced Model Approach | ||
|Jean-Luc | |Jean-Luc | ||
|- | |- | ||
|Oct | |Oct 11 '''Colloquium in B239 at 4:00pm''' | ||
|[ | |[https://people.math.ethz.ch/~imikaela/ Mikaela Iacobelli] (ETH/IAS) | ||
| | |[[# TBA| TBA ]] | ||
| | |Li | ||
|- | |- | ||
|Oct | |Oct 18 '''Colloquium in B239 at 4:00pm''' | ||
|[ | |[https://galton.uchicago.edu/~guillaumebal/ Guillaume Bal] (U Chicago) | ||
| | |[[#Bal| Speckle formation of laser light in random media: The Gaussian conjecture ]] | ||
| | | Li, Stechmann | ||
|- | |- | ||
|Oct | |Oct 23 ('''Wednesday''') | ||
|[ | |[https://www.sandia.gov/ccr/staff/teresa-portone/ Teresa Portone] (Sandia) | ||
| | |[[#Portone | Beyond parametric uncertainty: quantifying model-form uncertainty in model predictions ]] | ||
| | |Stechmann | ||
|- | |- | ||
|Oct | |Oct 25 | ||
|[ | |[https://www.cs.cornell.edu/~damle/ Anil Damle] (Cornell) | ||
| | |[[#Damle | Fine-grained Theory and Hybrid Algorithms for Randomized Numerical Linear Algebra ]] | ||
| | |Li | ||
|- | |- | ||
| | | Nov 1 | ||
|[ | |[https://research-hub.nrel.gov/en/persons/michael-sprague Michael Sprague] (NREL) | ||
| | |[[#Sprague| Exascale supercomputing and predictive wind energy simulations ]] | ||
| | |Spagnolie | ||
|- | |- | ||
| | | Nov 8 | ||
|[ | |[https://personal.math.ubc.ca/~holmescerfon/ Miranda Holmes-Cerfon] (UBC) | ||
| | |[[#Holmes-Cerfon | The dynamics of particles with ligand-receptor contacts ]] | ||
| | |Stechmann | ||
|- | |- | ||
| | | Nov 15* | ||
|[http:// | | [http://sun-yue.com Yue Sun] (UW–Madison) | ||
| | |[[#Holmes-Cerfon | Simulating fluid–structure interaction: A tale of two methods ]] | ||
| | | Rycroft | ||
|- | |- | ||
|Nov | | Nov 22 | ||
|[ | |[https://ibd.uchicago.edu/joinus/yenfellowship/ Ondrej Maxian] (U Chicago) | ||
| | |[[#Maxian | From slender body numerics to patterning the cell cortex: two stories of actin filament dynamics ]] | ||
| | |Ohm & Spagnolie | ||
|- | |- | ||
|Nov | | Nov 29* | ||
|''Thanksgiving'' | |||
|'' | |||
| | | | ||
| | |||
|- | |||
| Dec 6 | |||
|[https://www.simonsfoundation.org/people/ido-lavi/ Ido Lavi] (Flatiron) | |||
|[[#Lavi| Emergence of large-scale patterns in active matter: from nematic fluids to multicellular systems ]] | |||
|Spagnolie | |||
|} | |} | ||
< | 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 == | |||
===Nan Chen (UW–Madison)=== | |||
Title: Taming Uncertainty in a Complex World: The Rise of Uncertainty Quantification -- A Tutorial for Beginners | |||
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. | |||
===Katie Newhall (UNC Chapel Hill)=== | |||
Title: Energy landscapes, metastability, and transition paths | |||
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. | |||
===Indresan Govender (Mintek / Univ. of KwaZulu Natal, South Africa)=== | |||
Title: Granular flow modeling and visualization using nuclear imaging | |||
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. | |||
===Hongfei Chen (Tulane)=== | |||
Title: Investigating Hydrodynamics of Choanoflagellate Colonies: A Reduced Model Approach | |||
Abstract: Choanoflagellates, eukaryotes with a distinctive cellular structure consisting of a cell body, a flagellum, and a collar of microvilli, exhibit fascinating biological behavior. While many species exist as single cells, some form colonies, with the species ''C. Flexa'' standing out for its ability to dynamically transition its flagella between positions inside and outside the colony. | |||
Modeling the hydrodynamics of these colonies ideally requires detailed representations of each cell’s flagellum, microvilli, and body. However, the computational cost of simulating colonies with hundreds of cells makes this approach very expensive. To address this, we propose a reduced modeling framework that simplifies each cell to a force dipole while retaining key hydrodynamic features. | |||
Our force dipole model is calibrated against detailed computational simulations that account for the complete cellular structure. We show that this reduced model closely matches experimental data for non-deforming, free-swimming colonies. We further investigate how colony swimming and feeding performance depend on the flagellar position relative the colony, cell density, and overall colony shape. Finally, we explore the impact of the wall for flagella-in colonies, which are frequently observed in laboratory settings. | |||
<div id="Bal"> | |||
===Guillaume Bal (Chicago)=== | |||
Title: Speckle formation of laser light in random media: The Gaussian conjecture | |||
A widely accepted conjecture in the physical literature states that classical wave-fields propagating in random media over large distances eventually follow a complex circular Gaussian distribution. In this limit, the wave intensity becomes exponentially distributed, which corroborates the speckle patterns of, e.g., laser light observed in experiments. This talk reports on recent results settling the conjecture in the weak-coupling, paraxial regime of wave propagation. The limiting macroscopic Gaussian wave-field is fully characterized by a correlation function that satisfies an unusual diffusion equation. | |||
The paraxial model of wave propagation is an approximation of the Helmholtz model where backscattering has been neglected. It is mathematically simpler to analyze but quite accurate practically for wave-fields that maintain a beam-like structure as in the application of laser light propagating in turbulent atmospheres. | |||
The derivation of the limiting model is first obtained in the Itô-Schrödinger regime, where the random medium is replaced by its white noise limit. The resulting stochastic PDE has the main advantage that finite dimensional statistical moments of the wave-field satisfy closed form equations. The proof of the derivation of the macroscopic model is based on showing that these moment solutions are asymptotically those of the Gaussian limit, on obtaining a stochastic continuity (and tightness) result, and on establishing that moments in the paraxial and the Itô-Schrödinger regimes are asymptotically close. | |||
This is joint work with Anjali Nair. | |||
<div id="Portone"> | |||
===Teresa Portone (Sandia)=== | |||
Title: Beyond parametric uncertainty: quantifying model-form uncertainty in model predictions | |||
Uncertainty quantification (UQ) is the science of characterizing, quantifying, and reducing | |||
uncertainties in mathematical models. It is critical for informing decisions, because it provides a measure | |||
of confidence in model predictions, given the uncertainties present in the model. While approaches to | |||
characterize uncertainties in model parameters, boundary and initial conditions are well established, it is | |||
less clear how to address uncertainties arising when the equations of a mathematical model are | |||
themselves uncertain—that is, when there is model-form uncertainty. Model-form uncertainty often | |||
arises in models of complex physical phenomena where (1) simplifications for computational tractability | |||
or (2) lack of knowledge lead to unknowns in the governing equations for which appropriate | |||
mathematical forms are unknown or may not exist. In this talk, I briefly introduce major concepts in UQ, | |||
then I discuss approaches to characterize model-form uncertainty and its impact on model predictions. | |||
<div id="Damle"> | |||
=== Anil Damle (Cornell) === | |||
Title: Fine-grained Theory and Hybrid Algorithms for Randomized Numerical Linear Algebra | |||
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. | |||
<div id="Sprague"> | |||
=== Michael Sprague (NREL) === | |||
Title: Exascale supercomputing and predictive wind energy simulations | |||
The predictive simulation modern wind turbines and wind farms is a high-performance-computing (HPC) grand challenge. Wind turbines are the largest rotating machines in the world, with rotor diameters exceeding 200 meters, and with heights reaching well into the atmospheric boundary layer. To address this grand challenge, the U.S. Department of Energy (DOE) Wind Energy Technologies Office and the DOE Exascale Computing Project have been supporting the creation of the ExaWind modeling and simulation environment since 2016. ExaWind is composed of the incompressible-flow computational-fluid-dynamics (CFD) solvers AMR-Wind and Nalu-Wind and the wind-turbine-dynamics solver OpenFAST. ExaWind codes have been developed with performance portability as a priority, with the first U.S. exascale computer, Frontier, being our target platform. Frontier relies on graphical processing units (GPUs) for acceleration, which presents a major challenge to codes designed for CPUs. In this talk I will give a historical overview of the Exascale Computing Project, an eight-year $1.7 billion project. I will show results from our ExaWind challenge problem on Frontier and describe the strong-scaling challenges, and I will describe the challenges of modeling and simulating floating offshore wind turbines. I will also give my perspectives on life as a Research Scientist in Applied Mathematics at a DOE national laboratory. | |||
<div id="Holmes-Cerfon"> | |||
=== Miranda Holmes-Cerfon (UBC) === | |||
Title: The dynamics of particles with ligand-receptor contacts | |||
One way to glue objects together at the nanoscale or microscale is by ligand-receptor interactions, where short sticky hair-like ligands stick to receptors on another surface, much like velcro on the nanoscale. Such interactions are common in biological systems, such as white blood cells, virus particles, cargo in the nuclear pore complex, etc, and they are also useful in materials science, where coating colloids with single-stranded DNA creates particles with programmable interactions. In these systems, the ligand-receptor interactions not only hold particles together, but also influence their dynamics. How do such particles move? Do they “roll” on each others’ surfaces, as is commonly thought? Or could they slide? And does it matter? In this talk I will introduce our modelling and experimental efforts aimed at understanding the coarse-grained dynamics of particles with ligand-receptor interactions. Our models predict these interactions can change the particles' effective diffusion by orders of magnitude. Our experiments, using DNA-coated colloids, verify this dramatic dynamical slowdown, but also show other dynamical features not yet captured by our models, which suggest new avenues for exploration. | |||
<div id="Sun"> | |||
=== Yue Sun (UW–Madison) === | |||
Title: Simulating fluid–structure interaction: A tale of two methods | |||
Computational approaches have become essential for complementing experimental and theoretical methods in the study of fluid–structure interaction (FSI)—from matching specific experimental conditions to creating digital twins for exploring otherwise unattainable data, to developing adaptable domain-specific methods. In this talk, I will discuss our collaborative work on developing two FSI methods for experimental digital twin applications and general method development. | |||
The first part will highlight our collaboration with the Prigozhin Group at Harvard to create a 3D digital twin of the cryo-plunging process. Using ''cryoflo'', a massively parallelized code with adaptive mesh refinement (AMR) built on AMReX, we model fluid–structure interactions and heat transfer between biological samples and cryogen. This simulation captures critical cooling dynamics, providing insights that inform experimental protocols for cryo-electron microscopy (cryo-EM). | |||
The second part will focus on general method development. Over the past decade, Rycroft ''et al.'' introduced the reference map technique (RMT), a fully Eulerian method for modeling finite-strain deformation in FSI. Here, we integrate the RMT with the lattice Boltzmann (LB) method, introducing a new approach (LBRMT) to simulate finite-strain solids on the LB’s fixed Eulerian grid. We demonstrate LBRMT’s capabilities by modeling interactions among multiple solid structures in fluids, showcasing its adaptability for various FSI scenarios such as collective behavior in active and soft matter. | |||
<div id="Maxian"> | |||
=== Ondrej Maxian (U Chicago) === | |||
Title: From slender body numerics to patterning the cell cortex: two stories of actin filament dynamics | |||
Actin filaments are the main ingredient in the cell cytoskeleton, which controls cell division, motility, and structure. In this talk, I will present two projects whose shared goal is to determine how microscopic dynamics of actin shape larger-scale behaviors of the cell cortex. In the first part, I will detail a new general purpose simulation package for fiber dynamics which accounts for filament inextensibility, Brownian motion, and nonlocal hydrodynamics. I will focus in particular on how to formulate a mobility matrix (force-velocity relationship) which is positive definite (necessary for Brownian motion) and has cost independent of the fiber slenderness (necessary for efficient simulation), then demonstrate how the package can be used to simulate cross-linked actin networks and sedimenting fiber arrays. In the second part, I will present a model for how actin filaments shape their own homeostasis through biochemical coupling with the protein RhoA. I will introduce an activator-inhibitor model for RhoA/actin coupling, then use a Bayesian inverse framework to infer the distribution of actin dynamics parameters associated with experimental data in ''C. elegans'' and starfish embryos. The inferred parameter values demonstrate how varying actin kinetics can explain changing patterns of RhoA excitability observed across multiple experimental systems. | |||
<div id='Lavi'> | |||
=== Ido Lavi (Flatiron) === | |||
Title: Emergence of large-scale patterns in active matter: from nematic fluids to multicellular systems | |||
Active systems exhibit a fascinating interplay of chaos, order, and collective dynamics. If time permits, I will explore two stories of emergent behavior: the dynamical arrest in defect-free active nematics and the role of intercellular adhesions in coordinating multicellular systems. | |||
In the realm of active nematics, our large-scale simulations reveal that defect-free active turbulence can transition into a dynamically arrested state. We find that the flow alignment coupling, which determines how nematics reorient under shear, acts as a control parameter that tunes the contrast between contractile and extensile systems. In contractile systems, it amplifies chaotic jets, while in extensile systems, it promotes a tree-like network of persistent streams, aligned with a maze of nematic domain walls. These findings highlight a novel role for topological defects and invite experimental investigation of defect-free systems. They also suggest an intriguing mechanism by which chaos is harnessed to generate patterns. | |||
Turning to multicellular systems, adhesion molecules like E-cadherins play a key role in tissue mechanics by dynamically linking the cytoskeletons of neighboring cells. Far from being passive glue, these molecules couple actin filaments across cells, enabling coordinated collective behavior. To explore these dynamics, we developed a theoretical model describing cytoskeletons as contractile gels, with boundary conditions controlled by E-cadherin linkers attached to actin on both sides. Our finite element simulations reveal emergent patterns such as global polarization, anti-polarization, ring-like arrangements, and transient supracellular networks of actin cables. Beyond specific predictions, this study provides a mathematical framework for understanding how intracellular activity and intercellular adhesion feedback drive emergent multicellular dynamics. | |||
== Future semesters == | |||
*[[Applied/ACMS/Fall2024|Fall 2024]] | |||
*[[Applied/ACMS/Spring2025|Spring 2025]] | |||
== Archived semesters == | == Archived semesters == | ||
*[[Applied/ACMS/Spring2024|Spring 2024]] | |||
*[[Applied/ACMS/Fall2023|Fall 2023]] | |||
*[[Applied/ACMS/Spring2023|Spring 2023]] | |||
*[[Applied/ACMS/Fall2022|Fall 2022]] | |||
*[[Applied/ACMS/Spring2022|Spring 2022]] | |||
*[[Applied/ACMS/Fall2021|Fall 2021]] | |||
*[[Applied/ACMS/Spring2021|Spring 2021]] | |||
*[[Applied/ACMS/Fall2020|Fall 2020]] | |||
*[[Applied/ACMS/Spring2020|Spring 2020]] | |||
*[[Applied/ACMS/Fall2019|Fall 2019]] | |||
*[[Applied/ACMS/Spring2019|Spring 2019]] | |||
*[[Applied/ACMS/Fall2018|Fall 2018]] | |||
*[[Applied/ACMS/Spring2018|Spring 2018]] | |||
*[[Applied/ACMS/Fall2017|Fall 2017]] | |||
*[[Applied/ACMS/Spring2017|Spring 2017]] | |||
*[[Applied/ACMS/Fall2016|Fall 2016]] | |||
*[[Applied/ACMS/Spring2016|Spring 2016]] | |||
*[[Applied/ACMS/Fall2015|Fall 2015]] | |||
*[[Applied/ACMS/Spring2015|Spring 2015]] | |||
*[[Applied/ACMS/Fall2014|Fall 2014]] | |||
*[[Applied/ACMS/Spring2014|Spring 2014]] | |||
*[[Applied/ACMS/Fall2013|Fall 2013]] | |||
*[[Applied/ACMS/Spring2013|Spring 2013]] | |||
*[[Applied/ACMS/Fall2012|Fall 2012]] | |||
*[[Applied/ACMS/Spring2012|Spring 2012]] | *[[Applied/ACMS/Spring2012|Spring 2012]] | ||
*[[Applied/ACMS/Fall2011|Fall 2011]] | *[[Applied/ACMS/Fall2011|Fall 2011]] |
Latest revision as of 22:55, 26 November 2024
Applied and Computational Mathematics Seminar
- When: Fridays at 2:25pm (except as otherwise indicated)
- Where: 901 Van Vleck Hall
- Organizers: Maurice Fabien, Chris Rycroft, and Saverio Spagnolie,
- To join the ACMS mailing list: Send mail to acms+subscribe@g-groups.wisc.edu.
Fall 2024
Dates marked with an asterisk correspond to 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
Nan Chen (UW–Madison)
Title: Taming Uncertainty in a Complex World: The Rise of Uncertainty Quantification -- A Tutorial for Beginners
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.
Katie Newhall (UNC Chapel Hill)
Title: Energy landscapes, metastability, and transition paths
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.
Indresan Govender (Mintek / Univ. of KwaZulu Natal, South Africa)
Title: Granular flow modeling and visualization using nuclear imaging
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 IsaMillTM, 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.
Hongfei Chen (Tulane)
Title: Investigating Hydrodynamics of Choanoflagellate Colonies: A Reduced Model Approach
Abstract: Choanoflagellates, eukaryotes with a distinctive cellular structure consisting of a cell body, a flagellum, and a collar of microvilli, exhibit fascinating biological behavior. While many species exist as single cells, some form colonies, with the species C. Flexa standing out for its ability to dynamically transition its flagella between positions inside and outside the colony.
Modeling the hydrodynamics of these colonies ideally requires detailed representations of each cell’s flagellum, microvilli, and body. However, the computational cost of simulating colonies with hundreds of cells makes this approach very expensive. To address this, we propose a reduced modeling framework that simplifies each cell to a force dipole while retaining key hydrodynamic features.
Our force dipole model is calibrated against detailed computational simulations that account for the complete cellular structure. We show that this reduced model closely matches experimental data for non-deforming, free-swimming colonies. We further investigate how colony swimming and feeding performance depend on the flagellar position relative the colony, cell density, and overall colony shape. Finally, we explore the impact of the wall for flagella-in colonies, which are frequently observed in laboratory settings.
Guillaume Bal (Chicago)
Title: Speckle formation of laser light in random media: The Gaussian conjecture
A widely accepted conjecture in the physical literature states that classical wave-fields propagating in random media over large distances eventually follow a complex circular Gaussian distribution. In this limit, the wave intensity becomes exponentially distributed, which corroborates the speckle patterns of, e.g., laser light observed in experiments. This talk reports on recent results settling the conjecture in the weak-coupling, paraxial regime of wave propagation. The limiting macroscopic Gaussian wave-field is fully characterized by a correlation function that satisfies an unusual diffusion equation.
The paraxial model of wave propagation is an approximation of the Helmholtz model where backscattering has been neglected. It is mathematically simpler to analyze but quite accurate practically for wave-fields that maintain a beam-like structure as in the application of laser light propagating in turbulent atmospheres.
The derivation of the limiting model is first obtained in the Itô-Schrödinger regime, where the random medium is replaced by its white noise limit. The resulting stochastic PDE has the main advantage that finite dimensional statistical moments of the wave-field satisfy closed form equations. The proof of the derivation of the macroscopic model is based on showing that these moment solutions are asymptotically those of the Gaussian limit, on obtaining a stochastic continuity (and tightness) result, and on establishing that moments in the paraxial and the Itô-Schrödinger regimes are asymptotically close.
This is joint work with Anjali Nair.
Teresa Portone (Sandia)
Title: Beyond parametric uncertainty: quantifying model-form uncertainty in model predictions
Uncertainty quantification (UQ) is the science of characterizing, quantifying, and reducing uncertainties in mathematical models. It is critical for informing decisions, because it provides a measure of confidence in model predictions, given the uncertainties present in the model. While approaches to characterize uncertainties in model parameters, boundary and initial conditions are well established, it is less clear how to address uncertainties arising when the equations of a mathematical model are themselves uncertain—that is, when there is model-form uncertainty. Model-form uncertainty often arises in models of complex physical phenomena where (1) simplifications for computational tractability or (2) lack of knowledge lead to unknowns in the governing equations for which appropriate mathematical forms are unknown or may not exist. In this talk, I briefly introduce major concepts in UQ, then I discuss approaches to characterize model-form uncertainty and its impact on model predictions.
Anil Damle (Cornell)
Title: Fine-grained Theory and Hybrid Algorithms for Randomized Numerical Linear Algebra
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.
Michael Sprague (NREL)
Title: Exascale supercomputing and predictive wind energy simulations
The predictive simulation modern wind turbines and wind farms is a high-performance-computing (HPC) grand challenge. Wind turbines are the largest rotating machines in the world, with rotor diameters exceeding 200 meters, and with heights reaching well into the atmospheric boundary layer. To address this grand challenge, the U.S. Department of Energy (DOE) Wind Energy Technologies Office and the DOE Exascale Computing Project have been supporting the creation of the ExaWind modeling and simulation environment since 2016. ExaWind is composed of the incompressible-flow computational-fluid-dynamics (CFD) solvers AMR-Wind and Nalu-Wind and the wind-turbine-dynamics solver OpenFAST. ExaWind codes have been developed with performance portability as a priority, with the first U.S. exascale computer, Frontier, being our target platform. Frontier relies on graphical processing units (GPUs) for acceleration, which presents a major challenge to codes designed for CPUs. In this talk I will give a historical overview of the Exascale Computing Project, an eight-year $1.7 billion project. I will show results from our ExaWind challenge problem on Frontier and describe the strong-scaling challenges, and I will describe the challenges of modeling and simulating floating offshore wind turbines. I will also give my perspectives on life as a Research Scientist in Applied Mathematics at a DOE national laboratory.
Miranda Holmes-Cerfon (UBC)
Title: The dynamics of particles with ligand-receptor contacts
One way to glue objects together at the nanoscale or microscale is by ligand-receptor interactions, where short sticky hair-like ligands stick to receptors on another surface, much like velcro on the nanoscale. Such interactions are common in biological systems, such as white blood cells, virus particles, cargo in the nuclear pore complex, etc, and they are also useful in materials science, where coating colloids with single-stranded DNA creates particles with programmable interactions. In these systems, the ligand-receptor interactions not only hold particles together, but also influence their dynamics. How do such particles move? Do they “roll” on each others’ surfaces, as is commonly thought? Or could they slide? And does it matter? In this talk I will introduce our modelling and experimental efforts aimed at understanding the coarse-grained dynamics of particles with ligand-receptor interactions. Our models predict these interactions can change the particles' effective diffusion by orders of magnitude. Our experiments, using DNA-coated colloids, verify this dramatic dynamical slowdown, but also show other dynamical features not yet captured by our models, which suggest new avenues for exploration.
Yue Sun (UW–Madison)
Title: Simulating fluid–structure interaction: A tale of two methods
Computational approaches have become essential for complementing experimental and theoretical methods in the study of fluid–structure interaction (FSI)—from matching specific experimental conditions to creating digital twins for exploring otherwise unattainable data, to developing adaptable domain-specific methods. In this talk, I will discuss our collaborative work on developing two FSI methods for experimental digital twin applications and general method development.
The first part will highlight our collaboration with the Prigozhin Group at Harvard to create a 3D digital twin of the cryo-plunging process. Using cryoflo, a massively parallelized code with adaptive mesh refinement (AMR) built on AMReX, we model fluid–structure interactions and heat transfer between biological samples and cryogen. This simulation captures critical cooling dynamics, providing insights that inform experimental protocols for cryo-electron microscopy (cryo-EM).
The second part will focus on general method development. Over the past decade, Rycroft et al. introduced the reference map technique (RMT), a fully Eulerian method for modeling finite-strain deformation in FSI. Here, we integrate the RMT with the lattice Boltzmann (LB) method, introducing a new approach (LBRMT) to simulate finite-strain solids on the LB’s fixed Eulerian grid. We demonstrate LBRMT’s capabilities by modeling interactions among multiple solid structures in fluids, showcasing its adaptability for various FSI scenarios such as collective behavior in active and soft matter.
Ondrej Maxian (U Chicago)
Title: From slender body numerics to patterning the cell cortex: two stories of actin filament dynamics
Actin filaments are the main ingredient in the cell cytoskeleton, which controls cell division, motility, and structure. In this talk, I will present two projects whose shared goal is to determine how microscopic dynamics of actin shape larger-scale behaviors of the cell cortex. In the first part, I will detail a new general purpose simulation package for fiber dynamics which accounts for filament inextensibility, Brownian motion, and nonlocal hydrodynamics. I will focus in particular on how to formulate a mobility matrix (force-velocity relationship) which is positive definite (necessary for Brownian motion) and has cost independent of the fiber slenderness (necessary for efficient simulation), then demonstrate how the package can be used to simulate cross-linked actin networks and sedimenting fiber arrays. In the second part, I will present a model for how actin filaments shape their own homeostasis through biochemical coupling with the protein RhoA. I will introduce an activator-inhibitor model for RhoA/actin coupling, then use a Bayesian inverse framework to infer the distribution of actin dynamics parameters associated with experimental data in C. elegans and starfish embryos. The inferred parameter values demonstrate how varying actin kinetics can explain changing patterns of RhoA excitability observed across multiple experimental systems.
Ido Lavi (Flatiron)
Title: Emergence of large-scale patterns in active matter: from nematic fluids to multicellular systems
Active systems exhibit a fascinating interplay of chaos, order, and collective dynamics. If time permits, I will explore two stories of emergent behavior: the dynamical arrest in defect-free active nematics and the role of intercellular adhesions in coordinating multicellular systems.
In the realm of active nematics, our large-scale simulations reveal that defect-free active turbulence can transition into a dynamically arrested state. We find that the flow alignment coupling, which determines how nematics reorient under shear, acts as a control parameter that tunes the contrast between contractile and extensile systems. In contractile systems, it amplifies chaotic jets, while in extensile systems, it promotes a tree-like network of persistent streams, aligned with a maze of nematic domain walls. These findings highlight a novel role for topological defects and invite experimental investigation of defect-free systems. They also suggest an intriguing mechanism by which chaos is harnessed to generate patterns.
Turning to multicellular systems, adhesion molecules like E-cadherins play a key role in tissue mechanics by dynamically linking the cytoskeletons of neighboring cells. Far from being passive glue, these molecules couple actin filaments across cells, enabling coordinated collective behavior. To explore these dynamics, we developed a theoretical model describing cytoskeletons as contractile gels, with boundary conditions controlled by E-cadherin linkers attached to actin on both sides. Our finite element simulations reveal emergent patterns such as global polarization, anti-polarization, ring-like arrangements, and transient supracellular networks of actin cables. Beyond specific predictions, this study provides a mathematical framework for understanding how intracellular activity and intercellular adhesion feedback drive emergent multicellular dynamics.
Future semesters
Archived semesters
- Spring 2024
- Fall 2023
- Spring 2023
- Fall 2022
- Spring 2022
- Fall 2021
- Spring 2021
- Fall 2020
- Spring 2020
- Fall 2019
- Spring 2019
- Fall 2018
- Spring 2018
- Fall 2017
- Spring 2017
- Fall 2016
- Spring 2016
- Fall 2015
- Spring 2015
- Fall 2014
- Spring 2014
- Fall 2013
- Spring 2013
- Fall 2012
- Spring 2012
- Fall 2011
- Spring 2011
- Fall 2010
Return to the Applied Mathematics Group Page