Fromm's is C.F.D. scheme

\begin{displaymath}u_j^{n+1}=u_j^n+\frac{\Delta t} {\Delta
x}[F_{j-1/2}^{n+1/2}-F_{j+1/2}^{n+1/2}] \end{displaymath}

where

\begin{displaymath}F_{j+1/2}=a(u_j^n-1/2*(1-\sigma)\Delta u_j)
\end{displaymath}

where $\Delta u_j/\Delta x$ is the slope of the line segment in jth cell and $\Delta u_j=(u_{j+1}^n-u_{j-1}^n)/2$

Slope limiter is a slope limited method, we choose $\Delta u_j$ depending on the nearby values of u.

Idea: Don't allow the interpolation function

\begin{displaymath}u_I(x)=u_j^n+(x-x_j)/\Delta x *\Delta \widetilde{u_j}
\end{displaymath}

to introduce new MAX or MIN's in the discrete data

\begin{displaymath}U^n=(u_0,\cdots,u_{M-1}) \end{displaymath}

we can do this by evaluationg uI(x) at the edges of the jth cell xj-1/2 and xj+1/2 and compare with uj-1,uj,uj+1:

for eg.

\begin{displaymath}if u_I(x_{j+1/2})=u_j^n+1/2\Delta u_j >u_j,u_{j+1}
\end{displaymath}

(where $\Delta u_j=(u_{j+1}-u_{j-1})/2$,then we set $\Delta
\widetilde{u_j}=0 $consider all possible slopes you can get from uj-1,uj,uj+1 :

\begin{displaymath}\Delta u_j^L=u_j-u_{j-1}\qquad
\Delta u_j^R=u_{j+1}-u_j\qquad
\Delta u^C=(u_{j+1}-u_{j-1})/2
\end{displaymath}

(when you use $\Delta u_j=0$ in Fromm's method you get first order upwind.)

The Van Leer solpes $\Delta u_j^{VL}$ is given by:

\begin{displaymath}\Delta^{VL}u_j= \left \{ \begin {array} {r c}
S_j \cdot min(2...
..., if \quad lopes\quad change\quad sign
\end {array}
\right .
\end{displaymath}

where $\phi =(u^n_{j+1}-u_j^n)(u_j^n-u^n_{j-1})$

Sj=sign(uj+1-uj-1)

About this document ...

This document was generated using the LaTeX2HTML translator Version 98.1p1 release (March 2nd, 1998)

Copyright © 1993, 1994, 1995, 1996, 1997, Nikos Drakos, Computer Based Learning Unit, University of Leeds.

The command line arguments were:
latex2html -split 0 notes2-9.tex.

The translation was initiated by Wenlong Jin on 1999-02-10


Wenlong Jin
1999-02-10