This talk is motivated by work on the existence theory for viscoelasticity of Kelvin-Voigt type with non-convex stored energies (joint with K. Koumatos (U. of Sussex), C. Lattanzio and S. Spirito (U. of LAquila)), which shows propagation of H1-regularity for the deformation gradient of weak solutions for semiconvex stored energies. It turns out that weak solutions with deformation gradient in H1 are in fact unique in two-space dimensions, providing a striking analogy to corresponding results in the theory of 2D Euler equations with bounded vorticity.

While weak solutions still exist for initial data in L2, oscillations on the deformation gradi- ent can now persist and propagate in time. This can be seen via a counterexample indicating that for non-monotone stress-strain relations in 1-d oscillations of the strain lead to solutions with sustained oscillations. The existence of sustained oscillations in hyperbolic-parabolic system is then studied in several examples motivated by viscoelasticity and thermoviscoelas- ticity. Sufficient conditions for persistent oscillations are developed for linear problems, and examples in some nonlinear systems of interest. In several space dimensions oscillatory exam- ples are associated with lack of rank-one convexity of the stored energy. Nonlinear examples in models with thermal effects are also developed.

The formulation of Mean Field Games (MFG) via partial differential equations typically requires continuous differentiability of the Hamiltonian in order to determine the advective term in the Kolmogorov--Fokker--Planck equation for the density of players. However, in many cases of practical interest, the underlying optimal control problem may exhibit bang-bang controls, which typically lead to nondifferentiable Hamiltonians. In this talk we will present results on the analysis and numerical approximation of stationary second-order MFG systems for the general case of convex, Lipschitz, but possibly nondifferentiable Hamiltonians. In particular, we will propose a generalization of the MFG system as a Partial Differential Inclusion (PDI) based on interpreting the derivative of the Hamiltonian in terms of subdifferentials of convex functions. We present results that guarantee the existence of unique weak solutions to the stationary MFG PDI under a monotonicity condition similar to one that has been considered previously by Lasry and Lions. Moreover, we will propose a monotone finite element discretization of the weak formulation of the MFG PDI, and present results that confirm the strong H^1-norm convergence of the approximations to the value function and strong L^q-norm convergence of the approximations to the density function. The performance of the numerical method will be illustrated in experiments featuring nonsmooth solutions. This talk is based on joint work with my supervisor Iain Smears.

I will start with two very brief surveys.  First is a class of problems, namely variational inequalities (VIs), which generalize PDE problems, and second is a class of solver algorithms, namely full approximation storage (FAS) nonlinear multigrid for PDEs.  Motivation for applying FAS to VIs is demonstrated in the standard mathematical model for glacier surface evolution, a very general VI problem relevant to climate modeling.  (Residuals for this nonlinear and non-local VI problem are computed by solving a Stokes model.)  Some existing nonlinear multilevel VI schemes, based on global (Newton) linearization would seem to be less suited to such general VI problems.  From this context I will sketch some work-in-progress toward the scalable solutions of nonlinear and nonlocal VIs by an FAS-type multilevel method.

Arguably some of the most interesting phenomena in fluid dynamics, both from a mathematical and a physical perspective, stem from the interplay between a fluid and its boundaries. This talk will present some examples of how boundary effects lead to remarkable outcomes.  Singularities can form in finite time as a consequence of the continuum assumption when material surfaces are in smooth contact with horizontal boundaries of a fluid under gravity. For fluids with chemical solutes, the presence of boundaries impermeable to diffusion adds further dynamics which can give rise to self-induced flows and the formation of coherent structures out of scattered assemblies of immersed bodies. These effects can be analytically and numerically predicted by simple mathematical models and observed in “simple” experimental setups.