07-Jan-2022 19:53:20
fd1d_advection_ftcs_test():
MATLAB/Octave version 9.8.0.1380330 (R2020a) Update 2
Test fd1d_advection_ftcs().
FD1D_ADVECTION_FTCS:
MATLAB/Octave version 9.8.0.1380330 (R2020a) Update 2
Solve the constant-velocity advection equation in 1D,
du/dt = - c du/dx
over the interval:
0.0 <= x <= 1.0
with periodic boundary conditions, and
with a given initial condition
u(0,x) = (10x-4)^2 (6-10x)^2 for 0.4 <= x <= 0.6
= 0 elsewhere.
We use a method known as FTCS:
FT: Forward Time : du/dt = (u(t+dt,x)-u(t,x))/dt
CS: Centered Space: du/dx = (u(t,x+dx)-u(t,x-dx))/2/dx
Number of nodes NX = 101
Number of time steps NT = 1000
Constant velocity C = 1
FD1D_ADVECTION_FTCS
Normal end of execution.
Graphics saved as "fd1d_advection_ftcs_test.png".
fd1d_advection_ftcs_test():
Normal end of execution.
07-Jan-2022 19:53:24