\documentclass{article}
\usepackage{fullpage}
\usepackage{amsmath}
\usepackage{enumitem}
\usepackage{graphicx}

\newcommand{\bigo}{\ensuremath{\mathcal{O}}}

\begin{document}
\begin{center}
{\bf CME306 / CS205B Homework 7}
\end{center}
\section*{ENO-LLF}
\begin{equation} \label{eqn:burgers}
u_t + \left(\frac{u^2}{2}\right)_x = 0
\end{equation}
Consider Burgers' equation (above), discretized in a conservative way with forward-Euler:
\begin{equation}
\frac{u^{n+1}_i-u^n_i}{\Delta t} + \frac{\mathcal{F}^n_{i+1/2} - \mathcal{F}^n_{i-1/2}}{\Delta x}  = 0.
\end{equation}
Recall from lecture (and section) that we compute $\mathcal{F}^n_{i+1/2}$ the numerical flux via a variety of ENO schemes.

\begin{enumerate}
\item Recall that Lax-Friedrichs defines, for flux $\mathcal{F}_{i_0+1/2}$, $D_i^1H^\pm = f(\phi_i) \pm \alpha_{i_0+1/2}\phi_i$.  Show that this is equivalent to adding some viscosity to the solution and that this viscosity vanishes as $\Delta t, \Delta x \to 0$. \\ \\ \\ \\ \\ \\ \\ \\
\item In Local Lax-Friedrichs (or ENO-LLF), we take $\alpha_{i_0+1/2} = \max\{|\lambda_{i_0}|,|\lambda_{i_0+1}\}$.  A common mistake (one that shows up even in the literature) is to try to compute a global divided difference table as:
\begin{equation}
D_i^1H^\pm = f(u_i) \pm \alpha_{i+1/2}u_i
\end{equation}
Which is {\em very} different from the correct, {\bf local} divided difference table for $\mathcal{F}_{i_0+1/2}$,
\begin{equation}
D_i^1H^\pm = f(u_i) \pm \alpha_{i_0+1/2}u_i
\end{equation}
There are advantages to computing a global divided difference table, by saving the need to recompute quantities that've already been computed.  How can you do this for ENO-LLF while preserving the correct scheme?

\end{enumerate}

\pagebreak
\section*{Discrete Conservative Form}

\begin{enumerate}

\item Consider the following discretizations of Burger's equation~\eqref{eqn:burgers}.  Show that they are both consistent, but that only one of them is in conservative form.  Which one
  is in conservative form?
\begin{enumerate}
\item \[ \frac{u_j^{n+1} - u_j^n}{\Delta t} + \frac{(u_{j+1}^n)^2 - (u_j^n)^2}{2\Delta x} = 0 \]
\item \[ \frac{u_j^{n+1} - u_j^n}{\Delta t} + u_j^n\frac{u_{j+1}^n - u_j^n}{\Delta x} = 0 \]  \\ \\ \\ \\
\end{enumerate}
\item The equations given above are both consistent with Burger's equation, but only one of them gives shocks traveling at the right speed.  Why is this?
\\ \\ \\ \\
\item Consider the simple advection equation~\eqref{eqn:advection}.  Show that for this equation, ENO-Roe and ENO-LLF schemes will produce identical results.
\begin{equation} \label{eqn:advection}
u_t + au_x = 0
\end{equation}

\end{enumerate}

\end{document}
