12-Jun-2023 17:31:58

burgers_time_viscous_test():
  MATLAB/Octave version 5.2.0
  Test burgers_time_viscous().

btv_test01():
  Test burgers_time_viscous() with the gaussian initial condition.

  Initial condition: gaussian
  Number of space nodes = 81
  Number of time steps = 200
  Final time T_MAX = 2
  Viscosity = 0.01
  Boundary condition = 4
  Graphics saved as "btv_test01.png"

btv_test02():
  Test burgers_time_viscous() with gaussian initial condition.
  Now we use the clamped boundary condition.

  Initial condition: gaussian
  Number of space nodes = 81
  Number of time steps = 200
  Final time T_MAX = 2
  Viscosity = 0.01
  Boundary condition = 0
  Graphics saved as "btv_test02.png"

btv_test03():
  Test burgers_time_viscous() with gaussian initial condition.
  Use a Neumann condition on the right.

  Initial condition: gaussian
  Number of space nodes = 81
  Number of time steps = 200
  Final time T_MAX = 2
  Viscosity = 0.01
  Boundary condition = 3
  Graphics saved as "btv_test03.png"

btv_test04():
  Test burgers_time_viscous() with the shock initial condition.
  Use periodic boundaries.

  Initial condition: shock
  Number of space nodes = 81
  Number of time steps = 300
  Final time T_MAX = 3
  Viscosity = 0.01
  Boundary condition = 4
  Graphics saved as "btv_test04.png"

btv_test05():
  Test burgers_time_viscous() with the expansion initial condition.
  Use periodic boundaries.

  Initial condition: expansion
  Number of space nodes = 81
  Number of time steps = 200
  Final time T_MAX = 2
  Viscosity = 0.01
  Boundary condition = 4
  Graphics saved as "btv_test05.png"

btv_test06():
  Test burgers_time_viscous() with the spline initial condition.
  Use periodic boundaries.

  Initial condition: spline
  Number of space nodes = 81
  Number of time steps = 200
  Final time T_MAX = 2
  Viscosity = 0.01
  Boundary condition = 4
  Graphics saved as "btv_test06.png"

btv_test07():
  Test burgers_time_viscous() with the gaussian initial condition.
  Plot the solutions U(X,T) as a surface.

  Initial condition: gaussian
  Number of space nodes = 81
  Number of time steps = 200
  Final time T_MAX = 2
  Viscosity = 0.01
  Boundary condition = 4
  Graphics saved as "btv_test07.png"

btv_test08():
  Test BURGERS_TIME_VISCOUS with the spline initial condition.
  Use periodic boundaries.
  Plot the solution profile at the final time.

  Initial condition: spline
  Number of space nodes = 81
  Number of time steps = 200
  Final time T_MAX = 2
  Viscosity = 0.01
  Boundary condition = 4
  Graphics saved as "btv_test08.png"

btv_test09():
  Test burgers_time_viscous() with the spike initial condition.

  Initial condition: spike
  Number of space nodes = 81
  Number of time steps = 400
  Final time T_MAX = 4
  Viscosity = 0.01
  Boundary condition = 4
  Graphics saved as "btv_test09.png"

btv_test10():
  Test burgers_time_viscous() with the shock initial condition.
  BC: U(left)=0, U'(right)=0 boundary.

  Initial condition: shock
  Number of space nodes = 800
  Number of time steps = 10000
  Final time T_MAX = 2
  Viscosity = 0.01
  Boundary condition = 1
  Graphics saved as "btv_test10.png"

burgers_time_viscous_test():
  Normal end of execution.

12-Jun-2023 17:32:54