flame_ode_test
21-Jan-2019 10:06:10

FLAME_ODE_TEST:
  MATLAB version
  Test FLAME_ODE.

BASE_RUN
  Solve and plot the flame ODE for delta = 0.01.


  A copy of the plot was saved as "base_run.png"

BASE_RUN:
  Normal end of execution.

UNIFORM_RUN
  Vary the value of delta using a log-uniformly
  distributed factor between 1/2 and 2.
  Our quantity of interest Q is the time at which the
  solution reaches 0.99.
  Plot the solution curves, and plot Delta versus Q

  U = 0.629447, factor = 1.54697, DELTA = 0.0154697
  Y(T) = 0.99 at T = 72.3672
  U = 0.811584, factor = 1.75514, DELTA = 0.0175514
  Y(T) = 0.99 at T = 64.5763
  U = -0.746026, factor = 0.596244, DELTA = 0.00596244
  Y(T) = 0.99 at T = 176.404
  U = 0.826752, factor = 1.77369, DELTA = 0.0177369
  Y(T) = 0.99 at T = 63.9594
  U = 0.264718, factor = 1.2014, DELTA = 0.012014
  Y(T) = 0.99 at T = 91.2236
  U = -0.804919, factor = 0.572394, DELTA = 0.00572394
  Y(T) = 0.99 at T = 183.438
  U = -0.443004, factor = 0.735602, DELTA = 0.00735602
  Y(T) = 0.99 at T = 144.412
  U = 0.093763, factor = 1.06715, DELTA = 0.0106715
  Y(T) = 0.99 at T = 101.815
  U = 0.915014, factor = 1.88559, DELTA = 0.0188559
  Y(T) = 0.99 at T = 60.5635
  U = 0.929777, factor = 1.90498, DELTA = 0.0190498
  Y(T) = 0.99 at T = 60.0132
  Multiple solutions plotted in file "uniform_run.png".
  Quantity of Interest plotted in file "uniform_qoi.png".

UNIFORM_RUN
  Normal end of execution

QOI_QUAD
  Use quadrature to estimate the expected value of our
  quantity of interest Q, the time at which the
  combustion solution reaches 0.99.

  Using Clenshaw-Curtis quadrature of orders 1 through 12

   N      Estimated Q

   1       108.176
   2       133.169
   3       116.507
   4       116.326
   5       116.375
   6       116.367
   7       116.363
   8       116.367
   9       116.372
  10        116.36
  11       116.368
  12        116.37

QOI_QUAD
  Normal end of execution

FLAME_ODE_TEST:
  Normal end of execution.

21-Jan-2019 10:06:15
diary off