20110301-SIAMCSE.tex

%\documentclass[handout]{beamer}
\documentclass{beamer}

\mode<presentation>
{
\usetheme{default}
\usefonttheme[onlymath]{serif}
%\usetheme{Singapore}
%\usetheme{Warsaw}
%\usetheme{Malmoe}
% \useinnertheme{circles}
% \useoutertheme{infolines}
% \useinnertheme{rounded}

\setbeamercovered{transparent=5}
}

\usepackage[english]{babel}
\usepackage[latin1]{inputenc}
\usepackage{textpos,alltt,listings,multirow,ulem,siunitx}
\newcommand\hmmax{0}
\newcommand\bmmax{0}
\usepackage{bm}

% font definitions, try \usepackage{ae} instead of the following
% three lines if you don't like this look
\usepackage{mathptmx}
\usepackage[scaled=.90]{helvet}
%\usepackage{courier}
\usepackage[T1]{fontenc}
\usepackage{tikz}
\usetikzlibrary[shapes,shapes.arrows,arrows,shapes.misc,fit,positioning]

% \usepackage{pgfpages}
% \pgfpagesuselayout{4 on 1}[a4paper,landscape,border shrink=5mm]

\input{JedMacros.tex}

\title{Implicit solution of free surface flows in glaciology}

\author{Jed Brown}


% - Use the \inst command only if there are several affiliations.
% - Keep it simple, no one is interested in your street address.
\institute[ETH Z\"urich]
{
  Laboratory of Hydrology, Hydraulics, and Glaciology \\
  ETH Zürich
}

\date[2011-03-01]{2011-03-01}

% This is only inserted into the PDF information catalog. Can be left
% out.
\subject{Talks}


% If you have a file called "university-logo-filename.xxx", where xxx
% is a graphic format that can be processed by latex or pdflatex,
% resp., then you can add a logo as follows:

% \pgfdeclareimage[height=0.5cm]{university-logo}{university-logo-filename}
% \logo{\pgfuseimage{university-logo}}



% Delete this, if you do not want the table of contents to pop up at
% the beginning of each subsection:
% \AtBeginSubsection[]
% {
% \begin{frame}<beamer>
% \frametitle{Outline}
% \tableofcontents[currentsection,currentsubsection]
% \end{frame}
% }

% If you wish to uncover everything in a step-wise fashion, uncomment
% the following command:

%\beamerdefaultoverlayspecification{<+->}

\begin{document}
\lstset{language=C}
\normalem

\begin{frame}
\titlepage
\end{frame}

% \begin{frame}
%   \frametitle{Outline}
%   \tableofcontents
%   % You might wish to add the option [pausesections]
% \end{frame}

% \section{Practicalities}

\input{slides/GroundingLine/BindschadlerNaturalHistory.tex}

\section[Hydrostatic]{A robust multigrid solver for the hydrostatic equations}
\input{slides/THI/Equations.tex}
\begin{frame}
  \includegraphics[width=\textwidth]{figures/THI/x-shear}
\end{frame}
\begin{frame}
  \includegraphics[width=\textwidth]{figures/THI/x-80km-m16p2l6-ew} \\
  Grid-sequenced Newton-Krylov solution of test $X$.  The solid lines denote nonlinear iterations, and the dotted lines with $\times$ denote linear residuals.
\end{frame}
\input{slides/THI/Y5kmClip.tex}
\begin{frame}
  \begin{figure}
    \includegraphics[width=\textwidth]{figures/THI/y-10km-m10p6l5-ew}
    \centering\caption{Grid sequenced Newton-Krylov convergence for test $Y$.
    The ``cliff'' has \SI{58}{\degree} angle in the red line ($12\times 125$ meter elements), \SI{73}{\degree} for the cyan line ($6\times 62$ meter elements).}\label{fig:testy}
  \end{figure}
\end{frame}
\begin{frame}
  \begin{figure}
    \includegraphics[width=\textwidth]{figures/THI/linear4}
    \centering\caption{Average number of Krylov iterations per nonlinear iteration.  Each nonlinear system was solved to a relative tolerance of $10^{-2}$.}\label{fig:linear}
  \end{figure}
\end{frame}

\input{slides/THI/Shaheen.tex}
\begin{frame}
  \includegraphics[width=0.8\textwidth]{figures/THI/c-steady-crop}
\end{frame}

% \section{Other models}
% \input{slides/GroundingLine/SSA.tex}
% \begin{frame}{Non-Newtonian Stokes system: velocity $\bm u$, pressure $p$}
% \begin{columns}
% \begin{column}{0.5\textwidth}
%   \alert{\begin{align*}
%     -\nabla \cdot(\eta D\uu) + \nabla p - \ff &= 0 \\
%     \nabla \cdot \uu &= 0
%   \end{align*}}
% \end{column}
% \begin{column}{0.5\textwidth}
%     \begin{align*}
%       D\uu &= \tfrac 1 2 \left(\nabla \uu + (\nabla \uu)^T \right) \\
%       \gamma(D\uu) &= \tfrac 1 2 D\uu \tcolon D\uu \\
%       \eta(\gamma) &= B(\theta,\dotsc)\big(\gamma_0 + \gamma \big)^{\frac{\mathfrak{p}-2}{2}} \\
%       \mathfrak{p} &= 1 + \tfrac{1}{\mathfrak{n}} \approx \tfrac 4 3 \\
%       T &= \bm 1 - \bm n \otimes \bm n \\
%     \end{align*}
% \end{column}
% \end{columns}
% \vspace{-1.5em}
%     with boundary conditions
%     \begin{align*}
%       (\eta D\bm u - p\bm 1)\cdot\bm n =
%       \begin{cases}\bm 0 & \text{free surface} \\
%         -\rho_w z \bm n & \text{ice-ocean interface}\end{cases} \\
%       \bm u = \bm 0\qquad\qquad \text{frozen bed}, \theta < \theta_0 \\
%       \left. \begin{aligned}
%           \bm u \cdot \bm n &= \bm g_{\text{melt}}(T\uu,\dotsc) \\
%           T (\eta D\bm u - p\bm 1)\cdot\bm n &= \bm g_{\text{slip}}(T \bm u,\dotsc) \end{aligned}\right\}
%       \text{nonlinear slip}, \theta \ge \theta_0 \\
%     \end{align*}
%     \vspace{-3em}
%     \[ \bm g_{\text{slip}}(T\uu) = \beta_{\mathfrak{m}}(\dotsc) \lvert T\bm u \rvert^{\mathfrak{m}-1} T \bm u \]
%     Navier $\mathfrak{m}=1$, \quad Weertman $\mathfrak{m}\approx \frac 1 3$, \quad Coulomb $\mathfrak{m}=0$.
% \end{frame}

\input{slides/GroundingLine/ALEFormulation.tex}

\begin{frame}{Stokes challenges}
  \begin{block}{Mass conservation is critical}
    \begin{itemize}
    \item Staggered grid finite difference (hard to deal with geometry)
    \item Stabilized methods (conservation artifacts when non-smooth)
    \item Inf-sup stable mixed finite element method
      \begin{itemize}
      \item Use discontinuous pressure to enforce local mass conservation
      \item Inf-sup constant decays like $\sqrt{\epsilon}$ for $Q_k-P_{k-1}^{\text{disc}}$
      \item Sub-optimal order of accuracy for $Q_k-Q_{k-2}^{\text{disc}}$
      \end{itemize}
    \end{itemize}      
  \end{block}
  \begin{block}{Solving saddle-point problems}
    \begin{itemize}
    \item Not uniformly elliptic: solvers are much less robust
    \item Standard preconditioners do not work
    \item Coupled multigrid with Vanka smoothers offer best performance,\\
      \quad not robust for stretched grids or anisotropic viscosity
    \item Block preconditioners require approximate commutators, \\
      \quad fragile for strong anisotropy and non-smooth viscosity
    \end{itemize}
  \end{block}
\end{frame}

\section[Slip]{Slip boundary conditions on bumpy surfaces}
\input{slides/GroundingLine/Normals.tex}

\input{slides/GroundingLine/ALEBlockForm.tex}

\section{Coupling}
\input{slides/PETSc/Coupling.tex}

\begin{frame}{Outlook}
  \begin{block}{}
    \begin{itemize}
    \item We have textbook multigrid efficiency for hydrostatic equations
    \item Finally a good algebraic interface for tightly-coupled multiphysics
    \item Technical challenges for Stokes
    \item Local conservation is critical, well-balanced slip
    \item Singularities: reentrant corners, transition from frozen \\ \quad 
      to slip bounadry conditions, grounded margins, grounding lines
    \item Stiff geometric coupling terms
    \item IMEX time integration: additive Runge-Kutta
    \end{itemize}    
  \end{block}
  \begin{block}{Tools}
    \begin{itemize}
    \item PETSc\ \url{http://mcs.anl.gov/petsc}
      \begin{itemize}\item ML, Hypre, MUMPS
      \end{itemize}
    \item ITAPS \url{http://itaps.org}
      \begin{itemize}\item MOAB, CGM, Lasso
      \end{itemize}
    \end{itemize}
  \end{block}
\end{frame}

\end{document}