Thu Sep 13 11:07:13 2018 NS2DE_TEST Python version: 3.6.5 Test the NS2DE library. R8VEC_AMAX_TEST Python version: 3.6.5 R8VEC_AMAX computes the maximum absolute value entry in an R8VEC. Input vector: 0: -5.63163 1: 9.12635 2: 6.59018 3: 1.23391 4: -1.69386 5: -8.67763 6: -4.84844 7: -7.80086 8: -9.12342 9: 2.67931 Max Abs = 9.12635 R8VEC_AMAX_TEST: Normal end of execution. R8VEC_AMIN_TEST Python version: 3.6.5 R8VEC_AMIN computes the minimum absolute entry in an R8VEC. Input vector: 0: -5.63163 1: 9.12635 2: 6.59018 3: 1.23391 4: -1.69386 5: -8.67763 6: -4.84844 7: -7.80086 8: -9.12342 9: 2.67931 Min Abs = 1.23391 R8VEC_AMIN_TEST: Normal end of execution. R8VEC_MAX_TEST Python version: 3.6.5 R8VEC_MAX computes the maximum entry in an R8VEC. Input vector: 0: -5.63163 1: 9.12635 2: 6.59018 3: 1.23391 4: -1.69386 5: -8.67763 6: -4.84844 7: -7.80086 8: -9.12342 9: 2.67931 Max = 9.12635 R8VEC_MAX_TEST: Normal end of execution. R8VEC_MIN_TEST Python version: 3.6.5 R8VEC_MIN computes the minimum entry in an R8VEC. Input vector: 0: -5.63163 1: 9.12635 2: 6.59018 3: 1.23391 4: -1.69386 5: -8.67763 6: -4.84844 7: -7.80086 8: -9.12342 9: 2.67931 Min = -9.12342 R8VEC_MIN_TEST: Normal end of execution. R8VEC_NORM_L2_TEST Python version: 3.6.5 R8VEC_NORM_L2 computes the L2 norm of an R8VEC. Input vector: 0: -5.63163 1: 9.12635 2: 6.59018 3: 1.23391 4: -1.69386 5: -8.67763 6: -4.84844 7: -7.80086 8: -9.12342 9: 2.67931 L2 norm = 20.3201 R8VEC_NORM_L2_TEST: Normal end of execution. R8VEC_PRINT_TEST Python version: 3.6.5 R8VEC_PRINT prints an R8VEC. Here is an R8VEC: 0: 123.456 1: 5e-06 2: -1e+06 3: 3.14159 R8VEC_PRINT_TEST: Normal end of execution. R8VEC_UNIFORM_AB_TEST Python version: 3.6.5 R8VEC_UNIFORM_AB computes a random R8VEC. -1 <= X <= 5 Initial seed is 123456789 Random R8VEC: 0: 0.31051 1: 4.73791 2: 3.97706 3: 2.37017 4: 1.49184 5: -0.603288 6: 0.545467 7: -0.340259 8: -0.737026 9: 2.80379 R8VEC_UNIFORM_AB_TEST: Normal end of execution. GRID_2D_TEST: Python version: 3.6.5 Generate a regular grid. 0 0 0 10.000000 5.000000 1 1 0 12.500000 5.000000 2 2 0 15.000000 5.000000 3 3 0 17.500000 5.000000 4 4 0 20.000000 5.000000 5 0 1 10.000000 5.500000 6 1 1 12.500000 5.500000 7 2 1 15.000000 5.500000 8 3 1 17.500000 5.500000 9 4 1 20.000000 5.500000 10 0 2 10.000000 6.000000 11 1 2 12.500000 6.000000 12 2 2 15.000000 6.000000 13 3 2 17.500000 6.000000 14 4 2 20.000000 6.000000 GRID_2D_TEST: Normal end of execution. UVP_LUCAS_TEST Python version: 3.6.5 Lucas Bystricky Flow: Estimate the range of velocity and pressure at the initial time T = 0, over the unit square. Kinematic viscosity NU = 1 Fluid density RHO = 1 Minimum Maximum U: -0.318305 0.318303 V: -0.991674 -0.000268753 P: 0 0 UVP_LUCAS_TEST: Normal end of execution. UVP_LUCAS_TEST2 Python version: 3.6.5 Lucas Bystricky Flow: Estimate the range of velocity and pressure on the boundary at the initial time T = 0, over the unit square. Kinematic viscosity NU = 1 Fluid density RHO = 1 Minimum Maximum U: -0.31831 0.31831 V: -0.999874 -0 P: 0 0 UVP_LUCAS_TEST2: Normal end of execution. RHS_LUCAS_TEST Python version: 3.6.5 Lucas Bystricky Flow Sample the Navier-Stokes right hand sides at the initial time T = 0, over the unit square. Kinematic viscosity NU = 1 Fluid density RHO = 1 Minimum Maximum Ur: -3.1575 3.15751 Vr: -8.79482 0.935149 Pr: 0 0 RHS_LUCAS_TEST: Normal end of execution. RESID_LUCAS_TEST Python version: 3.6.5 Lucas Bystricky Flow Sample the Navier-Stokes residuals at the initial time T = 0, over the unit square. Kinematic viscosity NU = 1 Fluid density RHO = 1 Minimum Maximum Ur: 0 0 Vr: 0 0 Pr: 0 0 RESID_LUCAS_TEST: Normal end of execution. NS2DE_GNUPLOT_LUCAS_TEST: Python version: 3.6.5 Generate the Lucas Bystricky velocity field on a regular grid. Store in GNUPLOT data and command files. Data written to "lucas_data.txt". Commands written to "lucas_commands.txt". NS2DE_GNUPLOT_LUCAS_TEST: Normal end of execution. NS2DE_MATPLOTLIB_LUCAS_TEST: Python version: 3.6.5 Generate a Lucas Bystricky Flow field on a regular grid. Display it using MATPLOTLIB NS2DE_MATPLOTLIB_LUCAS_TEST: Normal end of execution. UVP_POISEUILLE_TEST Python version: 3.6.5 Poiseuille Flow: Estimate the range of velocity and pressure at the initial time T = 0, over a channel region. Kinematic viscosity NU = 1 Fluid density RHO = 1 Minimum Maximum U: 0.010254 0.999996 V: 0 0 P: -11.9749 -0.0220605 UVP_POISEUILLE_TEST: Normal end of execution. UVP_POISEUILLE_TEST2 Python version: 3.6.5 Poiseuille Flow: Estimate the range of velocity and pressure on the boundary at the initial time T = 0, over a channel region. Kinematic viscosity NU = 1 Fluid density RHO = 1 Minimum Maximum U: 0 0.999898 V: 0 0 P: -12 -0 UVP_POISEUILLE_TEST2: Normal end of execution. RHS_POISEUILLE_TEST Python version: 3.6.5 Poiseuille Flow: Sample the Navier-Stokes right hand sides at the initial time T = 0, over a channel region. Kinematic viscosity NU = 1 Fluid density RHO = 1 Minimum Maximum Ur: 0 0 Vr: 0 0 Pr: 0 0 RHS_POISEUILLE_TEST: Normal end of execution. RESID_POISEUILLE_TEST Python version: 3.6.5 Poiseuille Flow: Sample the Navier-Stokes residuals at the initial time T = 0, over a channel region. Kinematic viscosity NU = 1 Fluid density RHO = 1 Minimum Maximum Ur: 0 0 Vr: 0 0 Pr: 0 0 RESID_POISEUILLE_TEST: Normal end of execution. NS2DE_GNUPLOT_POISEUILLE_TEST: Python version: 3.6.5 Generate a Poiseuille velocity field on a regular grid. Store in GNUPLOT data and command files. Data written to "poiseuille_data.txt". Commands written to "poiseuille_commands.txt". NS2DE_GNUPLOT_POISEUILLE_TEST: Normal end of execution. NS2DE_MATPLOTLIB_POISEUILLE_TEST: Python version: 3.6.5 Generate a Poiseuille velocity field on a regular grid. Display it using MATPLOTLIB NS2DE_MATPLOTLIB_POISEUILLE_TEST: Normal end of execution. PARAMETER_POISEUILLE_TEST Python version: 3.6.5 Poiseuille Flow Monitor solution norms for various values of NU, RHO. RHO affects the pressure scaling. RHO NU T ||U|| ||V|| ||P|| 1 1 0 0.02296 0 0.2197 1 1 0.2 0.02296 0 0.2197 1 1 0.4 0.02296 0 0.2197 1 1 0.6 0.02296 0 0.2197 1 1 0.8 0.02296 0 0.2197 1 1 1 0.02296 0 0.2197 0.01 1 0 0.02296 0 0.002197 0.01 1 0.2 0.02296 0 0.002197 0.01 1 0.4 0.02296 0 0.002197 0.01 1 0.6 0.02296 0 0.002197 0.01 1 0.8 0.02296 0 0.002197 0.01 1 1 0.02296 0 0.002197 0.0001 1 0 0.02296 0 2.197e-05 0.0001 1 0.2 0.02296 0 2.197e-05 0.0001 1 0.4 0.02296 0 2.197e-05 0.0001 1 0.6 0.02296 0 2.197e-05 0.0001 1 0.8 0.02296 0 2.197e-05 0.0001 1 1 0.02296 0 2.197e-05 NU affects the time scaling. RHO NU T ||U|| ||V|| ||P|| 1 1 0 0.02296 0 0.2197 1 1 0.2 0.02296 0 0.2197 1 1 0.4 0.02296 0 0.2197 1 1 0.6 0.02296 0 0.2197 1 1 0.8 0.02296 0 0.2197 1 1 1 0.02296 0 0.2197 1 0.1 0 0.02296 0 0.02197 1 0.1 0.2 0.02296 0 0.02197 1 0.1 0.4 0.02296 0 0.02197 1 0.1 0.6 0.02296 0 0.02197 1 0.1 0.8 0.02296 0 0.02197 1 0.1 1 0.02296 0 0.02197 1 0.01 0 0.02296 0 0.002197 1 0.01 0.2 0.02296 0 0.002197 1 0.01 0.4 0.02296 0 0.002197 1 0.01 0.6 0.02296 0 0.002197 1 0.01 0.8 0.02296 0 0.002197 1 0.01 1 0.02296 0 0.002197 1 0.001 0 0.02296 0 0.0002197 1 0.001 0.2 0.02296 0 0.0002197 1 0.001 0.4 0.02296 0 0.0002197 1 0.001 0.6 0.02296 0 0.0002197 1 0.001 0.8 0.02296 0 0.0002197 1 0.001 1 0.02296 0 0.0002197 PARAMETER_POISEUILLE_TEST: Normal end of execution. UVP_SPIRAL_TEST Python version: 3.6.5 Spiral Flow: Estimate the range of velocity and pressure at the initial time T = 0, over the unit square. Kinematic viscosity NU = 1 Fluid density RHO = 1 Minimum Maximum U: -0.0119441 0.0119905 V: -0.0120032 0.011991 P: 0.00296933 0.99743 UVP_SPIRAL_TEST: Normal end of execution. UVP_SPIRAL_TEST2 Python version: 3.6.5 Spiral Flow: Estimate the range of velocity and pressure on the boundary at the initial time T = 0, over the unit square. Kinematic viscosity NU = 1 Fluid density RHO = 1 Minimum Maximum U: 0 0 V: -0 -0 P: 0 1 UVP_SPIRAL_TEST2: Normal end of execution. RHS_SPIRAL_TEST Python version: 3.6.5 Spiral Flow: Sample the Navier-Stokes right hand sides at the initial time T = 0, over the unit square. Kinematic viscosity NU = 1 Fluid density RHO = 1 Minimum Maximum Ur: -0.761154 0.763694 Vr: 0.23586 1.76464 Pr: 0 0 RHS_SPIRAL_TEST: Normal end of execution. RESID_SPIRAL_TEST Python version: 3.6.5 Spiral Flow: Sample the Navier-Stokes residuals at the initial time T = 0, over the unit square. Kinematic viscosity NU = 1 Fluid density RHO = 1 Minimum Maximum Ur: 0 0 Vr: 0 0 Pr: 0 0 RESID_SPIRAL_TEST: Normal end of execution. NS2DE_GNUPLOT_SPIRAL_TEST: Python version: 3.6.5 Generate a Spiral Flow velocity field on a regular grid. Store in GNUPLOT data and command files. Data written to "spiral_data.txt". Commands written to "spiral_commands.txt". NS2DE_GNUPLOT_SPIRAL_TEST: Normal end of execution. NS2DE_MATPLOTLIB_SPIRAL_TEST: Python version: 3.6.5 Generate a Spiral Flow velocity field on a regular grid. Display it using MATPLOTLIB NS2DE_MATPLOTLIB_SPIRAL_TEST: Normal end of execution. PARAMETER_SPIRAL_TEST Python version: 3.6.5 Spiral Flow Monitor solution norms over time for various values of NU, RHO. RHO affects the pressure scaling. RHO NU T ||U|| ||V|| ||P|| 1 1 0 0.0001767 0.0001712 0.01798 1 1 0.2 0.0002121 0.0002055 0.01798 1 1 0.4 0.0002474 0.0002397 0.01798 1 1 0.6 0.0002828 0.000274 0.01798 1 1 0.8 0.0003181 0.0003082 0.01798 1 1 1 0.0003534 0.0003425 0.01798 0.01 1 0 0.0001767 0.0001712 0.0001798 0.01 1 0.2 0.0002121 0.0002055 0.0001798 0.01 1 0.4 0.0002474 0.0002397 0.0001798 0.01 1 0.6 0.0002828 0.000274 0.0001798 0.01 1 0.8 0.0003181 0.0003082 0.0001798 0.01 1 1 0.0003534 0.0003425 0.0001798 0.0001 1 0 0.0001767 0.0001712 1.798e-06 0.0001 1 0.2 0.0002121 0.0002055 1.798e-06 0.0001 1 0.4 0.0002474 0.0002397 1.798e-06 0.0001 1 0.6 0.0002828 0.000274 1.798e-06 0.0001 1 0.8 0.0003181 0.0003082 1.798e-06 0.0001 1 1 0.0003534 0.0003425 1.798e-06 NU affects the time scaling. RHO NU T ||U|| ||V|| ||P|| 1 1 0 0.0001767 0.0001712 0.01798 1 1 0.2 0.0002121 0.0002055 0.01798 1 1 0.4 0.0002474 0.0002397 0.01798 1 1 0.6 0.0002828 0.000274 0.01798 1 1 0.8 0.0003181 0.0003082 0.01798 1 1 1 0.0003534 0.0003425 0.01798 1 0.1 0 0.0001767 0.0001712 0.01798 1 0.1 0.2 0.0001803 0.0001747 0.01798 1 0.1 0.4 0.0001838 0.0001781 0.01798 1 0.1 0.6 0.0001873 0.0001815 0.01798 1 0.1 0.8 0.0001909 0.0001849 0.01798 1 0.1 1 0.0001944 0.0001884 0.01798 1 0.01 0 0.0001767 0.0001712 0.01798 1 0.01 0.2 0.0001771 0.0001716 0.01798 1 0.01 0.4 0.0001774 0.0001719 0.01798 1 0.01 0.6 0.0001778 0.0001723 0.01798 1 0.01 0.8 0.0001781 0.0001726 0.01798 1 0.01 1 0.0001785 0.000173 0.01798 1 0.001 0 0.0001767 0.0001712 0.01798 1 0.001 0.2 0.0001768 0.0001713 0.01798 1 0.001 0.4 0.0001768 0.0001713 0.01798 1 0.001 0.6 0.0001768 0.0001714 0.01798 1 0.001 0.8 0.0001769 0.0001714 0.01798 1 0.001 1 0.0001769 0.0001714 0.01798 PARAMETER_SPIRAL_TEST: Normal end of execution. UVP_TAYLOR_TEST Python version: 3.6.5 Estimate the range of velocity and pressure at the initial time T = 0, using a region that is the square centered at (1.5,1.5) with "radius" 1.0, Kinematic viscosity NU = 1 Fluid density RHO = 1 Minimum Maximum U: -0.997425 0.992956 V: -0.998769 0.998312 P: -0.498148 0.498566 UVP_TAYLOR_TEST: Normal end of execution. UVP_TAYLOR_TEST2 Python version: 3.6.5 Estimate the range of velocity and pressure on the boundary at the initial time T = 0, using a region that is the square centered at (1.5,1.5) with "radius" 1.0, Kinematic viscosity NU = 1 Fluid density RHO = 1 Minimum Maximum U: -0.999874 0.999874 V: -0.999874 0.999874 P: 0.000125864 0.5 UVP_TAYLOR_TEST2: Normal end of execution. RHS_TAYLOR_TEST Python version: 3.6.5 Taylor Flow: Sample the Navier-Stokes right hand sides at the initial time T = 0, using a region that is the square centered at (1.5,1.5) with "radius" 1.0, Kinematic viscosity NU = 1 Fluid density RHO = 1 Minimum Maximum Ur: 0 0 Vr: 0 0 Pr: 0 0 RHS_TAYLOR_TEST: Normal end of execution. RESID_TAYLOR_TEST Python version: 3.6.5 Sample the Taylor residuals at the initial time T = 0, using a region that is the square centered at (1.5,1.5) with "radius" 1.0, Kinematic viscosity NU = 1 Fluid density RHO = 1 Minimum Maximum Ur: 0 3.55271e-15 Vr: 0 3.55271e-15 Pr: 0 0 RESID_TAYLOR_TEST: Normal end of execution. NS2DE_GNUPLOT_TAYLOR_TEST: Python version: 3.6.5 Generate a Taylor velocity field on a regular grid. Store in GNUPLOT data and command files. Data written to "taylor_data.txt". Commands written to "taylor_commands.txt". NS2DE_GNUPLOT_TAYLOR_TEST: Normal end of execution. NS2DE_MATPLOTLIB_TAYLOR_TEST: Python version: 3.6.5 Generate a Taylor velocity field on a regular grid. Display it using MATPLOTLIB NS2DE_MATPLOTLIB_TAYLOR_TEST: Normal end of execution. PARAMETER_TAYLOR_TEST Python version: 3.6.5 Taylor Flow Monitor solution norms over time for various values of NU, RHO. RHO affects the pressure scaling. RHO NU T ||U|| ||V|| ||P|| 1 1 0 0.01571 0.01599 0.007905 1 1 0.2 0.0003032 0.0003086 2.943e-06 1 1 0.4 5.851e-06 5.955e-06 1.096e-09 1 1 0.6 1.129e-07 1.149e-07 4.081e-13 1 1 0.8 2.179e-09 2.217e-09 1.519e-16 1 1 1 4.204e-11 4.278e-11 5.658e-20 0.01 1 0 0.01571 0.01599 7.905e-05 0.01 1 0.2 0.0003032 0.0003086 2.943e-08 0.01 1 0.4 5.851e-06 5.955e-06 1.096e-11 0.01 1 0.6 1.129e-07 1.149e-07 4.081e-15 0.01 1 0.8 2.179e-09 2.217e-09 1.519e-18 0.01 1 1 4.204e-11 4.278e-11 5.658e-22 0.0001 1 0 0.01571 0.01599 7.905e-07 0.0001 1 0.2 0.0003032 0.0003086 2.943e-10 0.0001 1 0.4 5.851e-06 5.955e-06 1.096e-13 0.0001 1 0.6 1.129e-07 1.149e-07 4.081e-17 0.0001 1 0.8 2.179e-09 2.217e-09 1.519e-20 0.0001 1 1 4.204e-11 4.278e-11 5.658e-24 NU affects the time scaling. RHO NU T ||U|| ||V|| ||P|| 1 1 0 0.01571 0.01599 0.007905 1 1 0.2 0.0003032 0.0003086 2.943e-06 1 1 0.4 5.851e-06 5.955e-06 1.096e-09 1 1 0.6 1.129e-07 1.149e-07 4.081e-13 1 1 0.8 2.179e-09 2.217e-09 1.519e-16 1 1 1 4.204e-11 4.278e-11 5.658e-20 1 0.1 0 0.01571 0.01599 0.007905 1 0.1 0.2 0.01059 0.01078 0.003589 1 0.1 0.4 0.007135 0.007261 0.00163 1 0.1 0.6 0.004808 0.004893 0.0007399 1 0.1 0.8 0.00324 0.003297 0.000336 1 0.1 1 0.002183 0.002222 0.0001525 1 0.01 0 0.01571 0.01599 0.007905 1 0.01 0.2 0.01511 0.01537 0.007305 1 0.01 0.4 0.01452 0.01478 0.00675 1 0.01 0.6 0.01396 0.01421 0.006238 1 0.01 0.8 0.01342 0.01366 0.005764 1 0.01 1 0.0129 0.01313 0.005327 1 0.001 0 0.01571 0.01599 0.007905 1 0.001 0.2 0.01565 0.01593 0.007843 1 0.001 0.4 0.01559 0.01587 0.007781 1 0.001 0.6 0.01553 0.0158 0.00772 1 0.001 0.8 0.01547 0.01574 0.007659 1 0.001 1 0.01541 0.01568 0.007599 PARAMETER_TAYLOR_TEST: Normal end of execution. NS2DE_TEST: Normal end of execution. Thu Sep 13 11:07:18 2018