previous up next SYLLABUS  Previous: 1.3 Prototype problems  Up: 1.3 Prototype problems  Next: 1.3.2 Diffusion


1.3.1 Advection


Slide : [ advection - waves || VIDEO login]

Also called convection, advection models the streaming of infinitesimal elements in a fluid. It generally appears when a transport process is modeled in a Eulerian representation using the convective derivative

$\displaystyle \frac{d}{dt}f \equiv \frac{\partial f}{\partial t} + u\frac{\partial f}{\partial x} = 0$ (1.3.1#eq.1)

For a constant advection velocity $ u$ , the advection equation can be solved analytically $ f(x,t)=f_0(x-ut)$ $ \forall f_0\in\mathcal{C}(\Omega)$ , showing explicitly the underlying characteristic $ x=ut$ . Try the JBONE applet below to compute the advection of a Gaussian pulse using the Lagrangian CIP method from chapter 6.
JBONE applet:  press Start/Stop to simulate the advection of a Gaussian pulse. Verify that the pulse indeed propagates with an advection velocity u=1 as chosen in the input parameters.



Numerical experiments: advection
  1. Switch the initial condition from a Gaussian to Cosine and Box functions. Observe how the sharp edges from the box disappear in the numerical solution; this is in clear contradiction with the analytical solution.
  2. Reduce the Velocity input parameter and verify that the advection speed is reduced accordingly. Increase the velocity above unity and determine the threshold above which the numerical solution clearly becomes unstable.

Both experiments show that numerical simulations have to be used carefully, to work withing limits of applicability that will be discussed in the comming sections. Note that for a constant advection speed, the wave equation can be written in flux-conservative form that reminds an advection

$\displaystyle \frac{\partial^2 h}{\partial t^2} -u^2 \frac{\partial^2 h}{\parti...
...}\right)\cdot \left(\begin{array}{c}\!\!f\\ \!\!g \end{array}\right) \right] =0$ (1.3.1#eq.2)

This shows explicitly that the numerical methods for the advection equation can in principle be used also for wave problems.

SYLLABUS  Previous: 1.3 Prototype problems  Up: 1.3 Prototype problems  Next: 1.3.2 Diffusion

      
back up next contents bibliography Copyright © Lifelong-learners at 03:22:40, December 17th, 2017