01-Jul-2023 19:20:06

burgers_steady_viscous_test():
  MATLAB/Octave version 5.2.0
  Test burgers_steady_viscous().

bsv_test01():
  Solution of steady viscous Burgers equation.

  Step  |F(U)|

     0  0.9
     1  1.3177
     2  0.140345
     3  0.00158582
     4  1.58645e-07
  Graphics saved as "bsv_test01.png".

  U(X0) = 0 estimated at X0 = 2.13852e-14

bsv_test02():
  Solution of steady viscous Burgers equation.
  Consider a variety of values of viscosity nu.

  Using NU = 0.8
  Using NU = 0.4
  Using NU = 0.2
  Using NU = 0.1
  Using NU = 0.05
  Using NU = 0.025

BURGERS_STEADY_VISCOUS - Warning!
  The Newton iteration did not converge.

  Graphics saved as "bsv_test02.png"

bsv_test03():
  Solution of steady viscous Burgers equation.
  Vary the left boundary condition ALPHA around the value +1.

  Using ALPHA = 0.96
  Using ALPHA = 0.98
  Using ALPHA = 0.99
  Using ALPHA = 0.995
  Using ALPHA = 1
  Using ALPHA = 1.005
  Using ALPHA = 1.01
  Using ALPHA = 1.02
  Using ALPHA = 1.04

  Graphics saved as "bsv_test03.png".

bsv_test04():
  Solution of steady viscous Burgers equation.
  Vary the left boundary location A around the value -1.

  Using A = -1.04
  Using A = -1.02
  Using A = -1.01
  Using A = -1.005
  Using A = -1
  Using A = -0.995
  Using A = -0.99
  Using A = -0.98
  Using A = -0.96

  Saved plot to file "bsv_test04.png".

bsv_test05():
  For the Burgers equation on [A,B] with viscosity NU and 
  boundary conditions U(A)=ALPHA, U(B) = BETA,
  with ALPHA and BETA of opposite sign,
  let X0 be the point where the solution U changes sign.
  Sample and plot the functional relationship X0(ALPHA).

  Saved plot to file "bsv_test05.png".

bsv_test06():
  For the Burgers equation on [A,B] with viscosity NU and 
  boundary conditions U(A)=ALPHA, U(B) = BETA,
  with ALPHA and BETA of opposite sign,
  let X0 be the point where the solution U changes sign.
  Assume ALPHA is Gaussian with mean 0 and standard deviation 0.05.
  Estimate E(X0(ALPHA)) using M Gaussian samples.

     M    E(X0(ALPHA)) estimate

    16       0.0618514
    32      -0.0160235
    64      -0.0240451
   128       0.0642659
   256       0.0292172
   512      -0.0186115
  1024        0.011543

bsv_test07():
  For the Burgers equation on [A,B] with viscosity NU and 
  boundary conditions U(A)=ALPHA, U(B) = BETA,
  with ALPHA and BETA of opposite sign,
  let X0 be the point where the solution U changes sign.
  Assume ALPHA is Gaussian with mean 0 and standard deviation 0.05.
  Estimate Var(X0(ALPHA)) using M Gaussian samples.

     M    Var(X0(ALPHA)) estimate

    16        0.314728
    32        0.324146
    64          0.3545
   128        0.341258
   256        0.352361
   512        0.335533
  1024        0.343424

bsv_test08():
  Compare bsv() and bsv_upwind().
  Upwinding is a scheme which reduces the numerical oscillations
  that can occur as the viscosity in the Burgers equation is decreased.

  The distortion caused by upwinding is visible for N = 21, NU = 0.1.
  The oscillations caused by NOT upwinding are visible for N = 21, NU = 0.01.

  Graphics saved as "bsv_test08.png".
  Graphics saved as "tanh_plot.png"

burgers_steady_viscous_test():
  Normal end of execution.

01-Jul-2023 19:22:11