6 October 2025 6:29:11.929 PM navier_stokes_2d_exact_test(): Fortran90 version Test navier_stokes_2d_exact(). uvp_cavity_test(): cavity flow Estimate the range of velocity and pressure at time T = 1.00000 over the interior of the [0,+1]x[0,+1] square. Kinematic viscosity NU = 1.00000 Fluid density RHO = 1.00000 Minimum Maximum U: -0.271337 0.992757 V: -0.383183 0.384406 P: -1745.53 2.71942 uvp_cavity_test2(): cavity flow Estimate the range of velocity and pressure at time T = 1.00000 over the boundary of the [0,+1]x[0,+1] square. Kinematic viscosity NU = 1.00000 Fluid density RHO = 1.00000 Minimum Maximum U: 0.00000 0.999796 V: -0.00000 -0.00000 P: -2045.33 2.88113 rhs_cavity_test(): cavity flow Sample the Navier-Stokes right hand sides at time T = 1.00000 over the interior of the [0,+1]x[0,+1] square. Kinematic viscosity NU = 1.00000 Fluid density RHO = 1.00000 Minimum Maximum Uf: -9526.55 3.86968 Vf: -7741.96 32.1482 Pf: 0.00000 0.00000 resid_cavity_test(): cavity flow Sample the Navier-Stokes residuals at time T = 1.00000 over the interior of the [0,+1]x[0,+1] square. Kinematic viscosity NU = 1.00000 Fluid density RHO = 1.00000 Minimum Maximum Ur: 0.00000 0.00000 Vr: 0.00000 0.00000 Pr: 0.00000 0.00000 velocity_cavity_test(): cavity flow Generate a velocity field on a regular grid. Store in GNUPLOT data and command files. Velocity scale factor S = 0.250000 Data written to "velocity_cavity_data.txt". Commands written to "velocity_cavity_commands.txt". uvp_exppoly_test(): exppoly flow Estimate the range of velocity and pressure at time T = 10.0000 over the interior of the [0,+1]x[0,+1] square. Kinematic viscosity NU = 1.00000 Fluid density RHO = 1.00000 Minimum Maximum U: -0.517205 0.515211 V: -0.516363 0.515634 P: -13885.9 8003.92 uvp_exppoly_test2(): exppoly flow Estimate the range of velocity and pressure at time T = 10.0000 over the boundary of the [0,+1]x[0,+1] square. Kinematic viscosity NU = 1.00000 Fluid density RHO = 1.00000 Minimum Maximum U: -0.00000 -0.00000 V: 0.00000 0.00000 P: -14022.5 8001.21 rhs_exppoly_test(): exppoly flow Sample the Navier-Stokes right hand sides at time T = 10.0000 over the interior of the [0,+1]x[0,+1] square. Kinematic viscosity NU = 1.00000 Fluid density RHO = 1.00000 Minimum Maximum Uf: -32.9100 32.8808 Vf: -69213.1 69206.5 Pf: 0.00000 0.00000 resid_exppoly_test(): exppoly flow Sample the Navier-Stokes residuals at time T = 10.0000 over the interior of the [0,+1]x[0,+1] square. Kinematic viscosity NU = 1.00000 Fluid density RHO = 1.00000 Minimum Maximum Ur: 0.00000 0.00000 Vr: 0.00000 0.00000 Pr: 0.00000 0.00000 velocity_exppoly_test(): exppoly flow Generate a velocity field on a regular grid. Store in GNUPLOT data and command files. Velocity scale factor S = 0.250000 Data written to "velocity_exppoly_data.txt". Commands written to "velocity_exppoly_commands.txt". uvp_exptrig_test(): exptrig flow Estimate the range of velocity and pressure at time T = 3.00000 over the interior of the [0,+1]x[0,+1] square. Kinematic viscosity NU = 1.00000 Fluid density RHO = 1.00000 Minimum Maximum U: -19.9834 19.9925 V: -19.8236 19.9999 P: -12.7747 7.29867 uvp_exptrig_test2(): exptrig flow Estimate the range of velocity and pressure at time T = 3.00000 over the boundary of the [0,+1]x[0,+1] square. Kinematic viscosity NU = 1.00000 Fluid density RHO = 1.00000 Minimum Maximum U: -0.491830E-14 0.301197E-30 V: -0.301197E-30 0.491830E-14 P: -12.7868 7.29616 rhs_exptrig_test(): exptrig flow Sample the Navier-Stokes right hand sides at time T = 3.00000 over the interior of the [0,+1]x[0,+1] square. Kinematic viscosity NU = 1.00000 Fluid density RHO = 1.00000 Minimum Maximum Uf: -2011.62 2014.45 Vf: -2042.94 2020.62 Pf: 0.00000 0.00000 resid_exptrig_test(): exptrig flow Sample the Navier-Stokes residuals at time T = 3.00000 over the interior of the [0,+1]x[0,+1] square. Kinematic viscosity NU = 1.00000 Fluid density RHO = 1.00000 Minimum Maximum Ur: 0.00000 0.00000 Vr: 0.00000 0.00000 Pr: 0.00000 0.00000 velocity_exptrig_test(): exptrig flow Generate a velocity field on a regular grid. Store in GNUPLOT data and command files. Velocity scale factor S = 0.474163E-02 Data written to "velocity_exptrig_data.txt". Commands written to "velocity_exptrig_commands.txt". uvp_gms_test(): gms flow Estimate the range of velocity and pressure at time T = 1.00000 over the interior of the [-1,+1]x[-1,+1] square. Kinematic viscosity NU = 1.00000 Fluid density RHO = 1.00000 Minimum Maximum U: -2.63917 2.63742 V: -2.64229 2.63874 P: -0.838318 0.840779 uvp_gms_test2(): gms flow Estimate the range of velocity and pressure at time T = 1.00000 over the boundary of the [-1,+1]x[-1,+1] square. Kinematic viscosity NU = 1.00000 Fluid density RHO = 1.00000 Minimum Maximum U: -0.647322E-15 0.647322E-15 V: -0.647322E-15 0.647322E-15 P: -0.841365 0.103050E-15 rhs_gms_test(): gms flow Sample the Navier-Stokes right hand sides at time T = 1.00000 over the interior of the [-1,+1]x[-1,+1] square. Kinematic viscosity NU = 1.00000 Fluid density RHO = 1.00000 Minimum Maximum Uf: -162.053 155.267 Vf: -157.207 160.660 Pf: 0.00000 0.00000 resid_gms_test(): gms flow Sample the Navier-Stokes residuals at time T = 1.00000 over the interior of the [-1,+1]x[-1,+1] square. Kinematic viscosity NU = 1.00000 Fluid density RHO = 1.00000 Minimum Maximum Ur: 0.00000 0.00000 Vr: 0.00000 0.00000 Pr: 0.00000 0.00000 velocity_gms_test(): gms flow Generate a velocity field on a regular grid. Store in GNUPLOT data and command files. Velocity scale factor S = 0.250000 Data written to "velocity_gms_data.txt". Commands written to "velocity_gms_commands.txt". uvp_lukas_test(): Lukas Bystricky flow Estimate the range of velocity and pressure at time T = 0.601347-153 using a region that is the unit square. Kinematic viscosity NU = 1.00000 Fluid density RHO = 1.00000 Minimum Maximum U: -0.318309 0.318309 V: -0.996972 -0.660024E-04 P: 0.00000 0.00000 uvp_lukas_test2(): Lukas Bystricky flow Estimate the range of velocity and pressure along the boundary, at time T = 0.00000 using a region that is the the unit square. Kinematic viscosity NU = 1.00000 Fluid density RHO = 1.00000 Minimum Maximum U: -0.318310 0.318310 V: -0.999874 -0.00000 P: 0.00000 0.00000 rhs_lukas_test(): Lukas Bystricky flow Sample the Navier-Stokes right hand sides at time T = 0.00000 using a region that is the unit square. Kinematic viscosity NU = 1.00000 Fluid density RHO = 1.00000 Minimum Maximum Uf: -3.15752 3.15752 Vf: -8.66418 0.739098 Pf: 0.00000 0.00000 resid_lukas_test() Lukas Bystricky flow Sample the Navier-Stokes residuals at time T = 0.395253-322 using a region that is the unit square. Kinematic viscosity NU = 1.00000 Fluid density RHO = 1.00000 Minimum Maximum Ur: 0.00000 0.00000 Vr: 0.00000 0.177636E-14 Pr: 0.00000 0.00000 velocity_lukas_test: Lukas Bystricky flow Generate a velocity field on a regular grid. Store in GNUPLOT data and command files. Velocity scale factor S = 0.250000 Data written to "velocity_lukas_data.txt". Commands written to "velocity_lukas_commands.txt". uvp_noflow_test(): noflow flow Estimate the range of velocity and pressure at time T = 1.00000 over the interior of the [-1,+1]x[-1,+1] square. Kinematic viscosity NU = 1.00000 Fluid density RHO = 1.00000 Minimum Maximum U: 0.00000 0.00000 V: 0.00000 0.00000 P: -0.774945 0.748334 uvp_noflow_test2(): noflow flow Estimate the range of velocity and pressure at time T = 1.00000 over the boundary of the [-1,+1]x[-1,+1] square. Kinematic viscosity NU = 1.00000 Fluid density RHO = 1.00000 Minimum Maximum U: 0.00000 0.00000 V: 0.00000 0.00000 P: -0.785398 0.785398 rhs_noflow_test(): noflow flow Sample the Navier-Stokes right hand sides at time T = 1.00000 over the interior of the [-1,+1]x[-1,+1] square. Kinematic viscosity NU = 1.00000 Fluid density RHO = 1.00000 Minimum Maximum Uf: -0.996304 0.989776 Vf: -0.971763 0.995023 Pf: 0.00000 0.00000 resid_noflow_test(): noflow flow Sample the Navier-Stokes residuals at time T = 1.00000 over the interior of the [-1,+1]x[-1,+1] square. Kinematic viscosity NU = 1.00000 Fluid density RHO = 1.00000 Minimum Maximum Ur: 0.00000 0.00000 Vr: 0.00000 0.00000 Pr: 0.00000 0.00000 19 September 2025_noflow_test: noflow flow Generate a pressure field on a regular grid. Store in GNUPLOT data and command files. Data written to "pressure_noflow_data.txt". Commands written to "pressure_noflow_commands.txt". uvp_poiseuille_test(): Poiseuille flow Estimate the range of velocity and pressure at time T = 0.601347-153 using a channel region. Kinematic viscosity NU = 1.00000 Fluid density RHO = 1.00000 Minimum Maximum U: 0.487872E-03 1.00000 V: 0.00000 0.00000 P: -11.9602 -0.263233E-02 uvp_poiseuille_test2(): Poiseuille flow Estimate the range of velocity and pressure along the boundary, at time T = 0.00000 using a channel region. Kinematic viscosity NU = 1.00000 Fluid density RHO = 1.00000 Minimum Maximum U: 0.00000 0.999898 V: 0.00000 0.00000 P: -12.0000 -0.00000 rhs_poiseuille_test(): Poiseuille flow Sample the Navier-Stokes right hand sides at time T = 0.00000 using a channel region. Kinematic viscosity NU = 1.00000 Fluid density RHO = 1.00000 Minimum Maximum Uf: 0.00000 0.00000 Vf: 0.00000 0.00000 Pf: 0.00000 0.00000 resid_poiseuille_test(): Poiseuille flow Sample the Navier-Stokes residuals at time T = 0.395253-322 using a channel region. Kinematic viscosity NU = 1.00000 Fluid density RHO = 1.00000 Minimum Maximum Ur: 0.00000 0.00000 Vr: 0.00000 0.00000 Pr: 0.00000 0.00000 velocity_poiseuille_test: Poiseuille flow Generate a velocity field on a regular grid. Store in GNUPLOT data and command files. Velocity scale factor S = 0.250000 Data written to "velocity_poiseuille_data.txt". Commands written to "velocity_poiseuille_commands.txt". PARAMETER_poiseuille_test Poiseuille Flow Monitor solution norms over time for various values of NU, RHO. Let RHO vary: RHO NU T ||U|| ||V|| ||P|| 1.000 1.000 0.0000 0.2287E-01 0.000 0.2145 1.000 1.000 0.2000 0.2287E-01 0.000 0.2145 1.000 1.000 0.4000 0.2287E-01 0.000 0.2145 1.000 1.000 0.6000 0.2287E-01 0.000 0.2145 1.000 1.000 0.8000 0.2287E-01 0.000 0.2145 1.000 1.000 1.0000 0.2287E-01 0.000 0.2145 0.1000E-01 1.000 0.0000 0.2287E-01 0.000 0.2145E-02 0.1000E-01 1.000 0.2000 0.2287E-01 0.000 0.2145E-02 0.1000E-01 1.000 0.4000 0.2287E-01 0.000 0.2145E-02 0.1000E-01 1.000 0.6000 0.2287E-01 0.000 0.2145E-02 0.1000E-01 1.000 0.8000 0.2287E-01 0.000 0.2145E-02 0.1000E-01 1.000 1.0000 0.2287E-01 0.000 0.2145E-02 0.1000E-03 1.000 0.0000 0.2287E-01 0.000 0.2145E-04 0.1000E-03 1.000 0.2000 0.2287E-01 0.000 0.2145E-04 0.1000E-03 1.000 0.4000 0.2287E-01 0.000 0.2145E-04 0.1000E-03 1.000 0.6000 0.2287E-01 0.000 0.2145E-04 0.1000E-03 1.000 0.8000 0.2287E-01 0.000 0.2145E-04 0.1000E-03 1.000 1.0000 0.2287E-01 0.000 0.2145E-04 Let NU vary: RHO NU T ||U|| ||V|| ||P|| 1.000 1.000 0.0000 0.2287E-01 0.000 0.2145 1.000 1.000 0.2000 0.2287E-01 0.000 0.2145 1.000 1.000 0.4000 0.2287E-01 0.000 0.2145 1.000 1.000 0.6000 0.2287E-01 0.000 0.2145 1.000 1.000 0.8000 0.2287E-01 0.000 0.2145 1.000 1.000 1.0000 0.2287E-01 0.000 0.2145 1.000 0.1000 0.0000 0.2287E-01 0.000 0.2145E-01 1.000 0.1000 0.2000 0.2287E-01 0.000 0.2145E-01 1.000 0.1000 0.4000 0.2287E-01 0.000 0.2145E-01 1.000 0.1000 0.6000 0.2287E-01 0.000 0.2145E-01 1.000 0.1000 0.8000 0.2287E-01 0.000 0.2145E-01 1.000 0.1000 1.0000 0.2287E-01 0.000 0.2145E-01 1.000 0.1000E-01 0.0000 0.2287E-01 0.000 0.2145E-02 1.000 0.1000E-01 0.2000 0.2287E-01 0.000 0.2145E-02 1.000 0.1000E-01 0.4000 0.2287E-01 0.000 0.2145E-02 1.000 0.1000E-01 0.6000 0.2287E-01 0.000 0.2145E-02 1.000 0.1000E-01 0.8000 0.2287E-01 0.000 0.2145E-02 1.000 0.1000E-01 1.0000 0.2287E-01 0.000 0.2145E-02 1.000 0.1000E-02 0.0000 0.2287E-01 0.000 0.2145E-03 1.000 0.1000E-02 0.2000 0.2287E-01 0.000 0.2145E-03 1.000 0.1000E-02 0.4000 0.2287E-01 0.000 0.2145E-03 1.000 0.1000E-02 0.6000 0.2287E-01 0.000 0.2145E-03 1.000 0.1000E-02 0.8000 0.2287E-01 0.000 0.2145E-03 1.000 0.1000E-02 1.0000 0.2287E-01 0.000 0.2145E-03 uvp_spiral_test(): Spiral flow Estimate the range of velocity and pressure at time T = 1.00000 over the unit square. Kinematic viscosity NU = 1.00000 Fluid density RHO = 1.00000 Minimum Maximum U: -0.120280E-01 0.120040E-01 V: -0.119688E-01 0.120020E-01 P: 0.800985E-03 0.999845 uvp_spiral_test2(): Spiral flow Estimate the range of velocity and pressure along the boundary, at time T = 0.00000 over the unit square. Kinematic viscosity NU = 1.00000 Fluid density RHO = 1.00000 Minimum Maximum U: 0.00000 0.00000 V: -0.00000 -0.00000 P: 0.00000 1.00000 rhs_spiral_test(): Spiral flow Sample the Navier-Stokes right hand sides at time T = 0.00000 over the unit square. Kinematic viscosity NU = 1.00000 Fluid density RHO = 1.00000 Minimum Maximum Uf: -0.765042 0.758558 Vf: 0.234496 1.76390 Pf: 0.00000 0.00000 resid_spiral_test(): Spiral flow Sample the Navier-Stokes residuals at time T = 0.395253-322 over the unit square. Kinematic viscosity NU = 1.00000 Fluid density RHO = 1.00000 Minimum Maximum Ur: 0.00000 0.111022E-15 Vr: 0.00000 0.444089E-15 Pr: 0.00000 0.00000 velocity_spiral_test: Spiral Flow Generate a velocity field on a regular grid. Store in GNUPLOT data and command files. Velocity scale factor S = 5.00000 Data written to "velocity_spiral_data.txt". Commands written to "velocity_spiral_commands.txt". PARAMETER_spiral_test Spiral Flow Monitor solution norms over time for various values of NU, RHO. Let RHO vary: RHO NU T ||U|| ||V|| ||P|| 1.000 1.000 0.0000 0.1666E-03 0.1784E-03 0.1809E-01 1.000 1.000 0.2000 0.1999E-03 0.2141E-03 0.1809E-01 1.000 1.000 0.4000 0.2332E-03 0.2498E-03 0.1809E-01 1.000 1.000 0.6000 0.2665E-03 0.2855E-03 0.1809E-01 1.000 1.000 0.8000 0.2998E-03 0.3211E-03 0.1809E-01 1.000 1.000 1.0000 0.3331E-03 0.3568E-03 0.1809E-01 0.1000E-01 1.000 0.0000 0.1666E-03 0.1784E-03 0.1809E-03 0.1000E-01 1.000 0.2000 0.1999E-03 0.2141E-03 0.1809E-03 0.1000E-01 1.000 0.4000 0.2332E-03 0.2498E-03 0.1809E-03 0.1000E-01 1.000 0.6000 0.2665E-03 0.2855E-03 0.1809E-03 0.1000E-01 1.000 0.8000 0.2998E-03 0.3211E-03 0.1809E-03 0.1000E-01 1.000 1.0000 0.3331E-03 0.3568E-03 0.1809E-03 0.1000E-03 1.000 0.0000 0.1666E-03 0.1784E-03 0.1809E-05 0.1000E-03 1.000 0.2000 0.1999E-03 0.2141E-03 0.1809E-05 0.1000E-03 1.000 0.4000 0.2332E-03 0.2498E-03 0.1809E-05 0.1000E-03 1.000 0.6000 0.2665E-03 0.2855E-03 0.1809E-05 0.1000E-03 1.000 0.8000 0.2998E-03 0.3211E-03 0.1809E-05 0.1000E-03 1.000 1.0000 0.3331E-03 0.3568E-03 0.1809E-05 Let NU vary: RHO NU T ||U|| ||V|| ||P|| 1.000 1.000 0.0000 0.1666E-03 0.1784E-03 0.1809E-01 1.000 1.000 0.2000 0.1999E-03 0.2141E-03 0.1809E-01 1.000 1.000 0.4000 0.2332E-03 0.2498E-03 0.1809E-01 1.000 1.000 0.6000 0.2665E-03 0.2855E-03 0.1809E-01 1.000 1.000 0.8000 0.2998E-03 0.3211E-03 0.1809E-01 1.000 1.000 1.0000 0.3331E-03 0.3568E-03 0.1809E-01 1.000 0.1000 0.0000 0.1666E-03 0.1784E-03 0.1809E-01 1.000 0.1000 0.2000 0.1699E-03 0.1820E-03 0.1809E-01 1.000 0.1000 0.4000 0.1732E-03 0.1855E-03 0.1809E-01 1.000 0.1000 0.6000 0.1766E-03 0.1891E-03 0.1809E-01 1.000 0.1000 0.8000 0.1799E-03 0.1927E-03 0.1809E-01 1.000 0.1000 1.0000 0.1832E-03 0.1963E-03 0.1809E-01 1.000 0.1000E-01 0.0000 0.1666E-03 0.1784E-03 0.1809E-01 1.000 0.1000E-01 0.2000 0.1669E-03 0.1788E-03 0.1809E-01 1.000 0.1000E-01 0.4000 0.1672E-03 0.1791E-03 0.1809E-01 1.000 0.1000E-01 0.6000 0.1676E-03 0.1795E-03 0.1809E-01 1.000 0.1000E-01 0.8000 0.1679E-03 0.1798E-03 0.1809E-01 1.000 0.1000E-01 1.0000 0.1682E-03 0.1802E-03 0.1809E-01 1.000 0.1000E-02 0.0000 0.1666E-03 0.1784E-03 0.1809E-01 1.000 0.1000E-02 0.2000 0.1666E-03 0.1784E-03 0.1809E-01 1.000 0.1000E-02 0.4000 0.1666E-03 0.1785E-03 0.1809E-01 1.000 0.1000E-02 0.6000 0.1667E-03 0.1785E-03 0.1809E-01 1.000 0.1000E-02 0.8000 0.1667E-03 0.1786E-03 0.1809E-01 1.000 0.1000E-02 1.0000 0.1667E-03 0.1786E-03 0.1809E-01 uvp_taylor_test(): Taylor flow Estimate the range of velocity and pressure at time T = 1.00000 using a region that is the square centered at (1.5,1.5) with "radius" 1.0, Kinematic viscosity NU = 1.00000 Fluid density RHO = 1.00000 Minimum Maximum U: -0.997599 0.997441 V: -0.998723 0.998932 P: -0.499144 0.498795 uvp_taylor_test2() Taylor flow Estimate the range of velocity and pressure along the boundary, at time T = 0.00000 using a region that is the square centered at (1.5,1.5) with "radius" 1.0, Kinematic viscosity NU = 1.00000 Fluid density RHO = 1.00000 Minimum Maximum U: -0.999874 0.999874 V: -0.999874 0.999874 P: 0.125864E-03 0.500000 rhs_taylor_test(): Taylor flow Sample the Navier-Stokes right hand sides at time T = 0.00000 using a region that is the square centered at (1.5,1.5) with "radius" 1.0, Kinematic viscosity NU = 1.00000 Fluid density RHO = 1.00000 Minimum Maximum Uf: 0.00000 0.00000 Vf: 0.00000 0.00000 Pf: 0.00000 0.00000 resid_taylor_test(): Taylor flow Sample the Navier-Stokes residuals at time T = 0.395253-322 using a region that is the square centered at (1.5,1.5) with "radius" 1.0, Kinematic viscosity NU = 1.00000 Fluid density RHO = 1.00000 Minimum Maximum Ur: 0.00000 0.355271E-14 Vr: 0.00000 0.355271E-14 Pr: 0.00000 0.00000 velocity_taylor_test: Taylor flow Generate a velocity field on a regular grid. Store in GNUPLOT data and command files. Velocity scale factor S = 0.100000 Data written to "velocity_taylor_data.txt". Commands written to "velocity_taylor_commands.txt". PARAMETER_taylor_test 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.000 1.000 0.0000 0.1642E-01 0.1578E-01 0.7546E-02 1.000 1.000 0.2000 0.3169E-03 0.3044E-03 0.2810E-05 1.000 1.000 0.4000 0.6115E-05 0.5874E-05 0.1046E-08 1.000 1.000 0.6000 0.1180E-06 0.1133E-06 0.3895E-12 1.000 1.000 0.8000 0.2277E-08 0.2187E-08 0.1450E-15 1.000 1.000 1.0000 0.4394E-10 0.4220E-10 0.5400E-19 0.1000E-01 1.000 0.0000 0.1642E-01 0.1578E-01 0.7546E-04 0.1000E-01 1.000 0.2000 0.3169E-03 0.3044E-03 0.2810E-07 0.1000E-01 1.000 0.4000 0.6115E-05 0.5874E-05 0.1046E-10 0.1000E-01 1.000 0.6000 0.1180E-06 0.1133E-06 0.3895E-14 0.1000E-01 1.000 0.8000 0.2277E-08 0.2187E-08 0.1450E-17 0.1000E-01 1.000 1.0000 0.4394E-10 0.4220E-10 0.5400E-21 0.1000E-03 1.000 0.0000 0.1642E-01 0.1578E-01 0.7546E-06 0.1000E-03 1.000 0.2000 0.3169E-03 0.3044E-03 0.2810E-09 0.1000E-03 1.000 0.4000 0.6115E-05 0.5874E-05 0.1046E-12 0.1000E-03 1.000 0.6000 0.1180E-06 0.1133E-06 0.3895E-16 0.1000E-03 1.000 0.8000 0.2277E-08 0.2187E-08 0.1450E-19 0.1000E-03 1.000 1.0000 0.4394E-10 0.4220E-10 0.5400E-23 NU affects the time scaling. RHO NU T ||U|| ||V|| ||P|| 1.000 1.000 0.0000 0.1642E-01 0.1578E-01 0.7546E-02 1.000 1.000 0.2000 0.3169E-03 0.3044E-03 0.2810E-05 1.000 1.000 0.4000 0.6115E-05 0.5874E-05 0.1046E-08 1.000 1.000 0.6000 0.1180E-06 0.1133E-06 0.3895E-12 1.000 1.000 0.8000 0.2277E-08 0.2187E-08 0.1450E-15 1.000 1.000 1.0000 0.4394E-10 0.4220E-10 0.5400E-19 1.000 0.1000 0.0000 0.1642E-01 0.1578E-01 0.7546E-02 1.000 0.1000 0.2000 0.1107E-01 0.1063E-01 0.3426E-02 1.000 0.1000 0.4000 0.7457E-02 0.7162E-02 0.1556E-02 1.000 0.1000 0.6000 0.5025E-02 0.4826E-02 0.7063E-03 1.000 0.1000 0.8000 0.3386E-02 0.3252E-02 0.3207E-03 1.000 0.1000 1.0000 0.2281E-02 0.2191E-02 0.1456E-03 1.000 0.1000E-01 0.0000 0.1642E-01 0.1578E-01 0.7546E-02 1.000 0.1000E-01 0.2000 0.1579E-01 0.1516E-01 0.6973E-02 1.000 0.1000E-01 0.4000 0.1518E-01 0.1458E-01 0.6443E-02 1.000 0.1000E-01 0.6000 0.1459E-01 0.1401E-01 0.5954E-02 1.000 0.1000E-01 0.8000 0.1402E-01 0.1347E-01 0.5502E-02 1.000 0.1000E-01 1.0000 0.1348E-01 0.1295E-01 0.5084E-02 1.000 0.1000E-02 0.0000 0.1642E-01 0.1578E-01 0.7546E-02 1.000 0.1000E-02 0.2000 0.1636E-01 0.1571E-01 0.7486E-02 1.000 0.1000E-02 0.4000 0.1629E-01 0.1565E-01 0.7427E-02 1.000 0.1000E-02 0.6000 0.1623E-01 0.1559E-01 0.7369E-02 1.000 0.1000E-02 0.8000 0.1617E-01 0.1553E-01 0.7311E-02 1.000 0.1000E-02 1.0000 0.1610E-01 0.1547E-01 0.7253E-02 uvp_vortex_test(): Vortex flow Estimate the range of velocity and pressure at time T = 2.50000 using a region that is the square centered at (1.5,1.5) with "radius" 1.0, Kinematic viscosity NU = 1.00000 Fluid density RHO = 1.00000 Minimum Maximum U: -0.999358 0.989479 V: -0.999817 0.999625 P: -0.499616 0.498737 uvp_vortex_test2(): Vortex flow Estimate the range of velocity and pressure along the boundary, at time T = 0.00000 using a region that is the square centered at (1.5,1.5) with "radius" 1.0, Kinematic viscosity NU = 1.00000 Fluid density RHO = 1.00000 Minimum Maximum U: -0.999874 0.999874 V: -0.999874 0.999874 P: 0.125864E-03 0.500000 rhs_vortex_test() Vortex flow Sample the Navier-Stokes right hand sides at time T = 0.00000 using a region that is the square centered at (1.5,1.5) with "radius" 1.0, Kinematic viscosity NU = 1.00000 Fluid density RHO = 1.00000 Minimum Maximum Uf: -19.6989 19.7180 Vf: -19.6791 19.7185 Pf: 0.00000 0.00000 resid_vortex_test(): Vortex flow Sample the Navier-Stokes residuals at time T = 0.395253-322 using a region that is the square centered at (1.5,1.5) with "radius" 1.0, Kinematic viscosity NU = 1.00000 Fluid density RHO = 1.00000 Minimum Maximum Ur: 0.555757E-02 19.7367 Vr: 0.407877E-04 19.6558 Pr: 0.00000 0.00000 velocity_vortex_test: Vortex flow Generate a velocity field on a regular grid. Store in GNUPLOT data and command files. Velocity scale factor S = 0.100000 Data written to "velocity_vortex_data.txt". Commands written to "velocity_vortex_commands.txt". PARAMETER_vortex_test Vortex flow Monitor solution norms over time for various values of NU, RHO. RHO affects the pressure scaling. RHO NU T ||U|| ||V|| ||P|| 1.000 1.000 0.0000 0.1545E-01 0.1613E-01 0.7915E-02 1.000 1.000 0.2000 0.1545E-01 0.1613E-01 0.7915E-02 1.000 1.000 0.4000 0.1545E-01 0.1613E-01 0.7915E-02 1.000 1.000 0.6000 0.1545E-01 0.1613E-01 0.7915E-02 1.000 1.000 0.8000 0.1545E-01 0.1613E-01 0.7915E-02 1.000 1.000 1.0000 0.1545E-01 0.1613E-01 0.7915E-02 0.1000E-01 1.000 0.0000 0.1545E-01 0.1613E-01 0.7915E-04 0.1000E-01 1.000 0.2000 0.1545E-01 0.1613E-01 0.7915E-04 0.1000E-01 1.000 0.4000 0.1545E-01 0.1613E-01 0.7915E-04 0.1000E-01 1.000 0.6000 0.1545E-01 0.1613E-01 0.7915E-04 0.1000E-01 1.000 0.8000 0.1545E-01 0.1613E-01 0.7915E-04 0.1000E-01 1.000 1.0000 0.1545E-01 0.1613E-01 0.7915E-04 0.1000E-03 1.000 0.0000 0.1545E-01 0.1613E-01 0.7915E-06 0.1000E-03 1.000 0.2000 0.1545E-01 0.1613E-01 0.7915E-06 0.1000E-03 1.000 0.4000 0.1545E-01 0.1613E-01 0.7915E-06 0.1000E-03 1.000 0.6000 0.1545E-01 0.1613E-01 0.7915E-06 0.1000E-03 1.000 0.8000 0.1545E-01 0.1613E-01 0.7915E-06 0.1000E-03 1.000 1.0000 0.1545E-01 0.1613E-01 0.7915E-06 NU affects the time scaling. RHO NU T ||U|| ||V|| ||P|| 1.000 1.000 0.0000 0.1545E-01 0.1613E-01 0.7915E-02 1.000 1.000 0.2000 0.1545E-01 0.1613E-01 0.7915E-02 1.000 1.000 0.4000 0.1545E-01 0.1613E-01 0.7915E-02 1.000 1.000 0.6000 0.1545E-01 0.1613E-01 0.7915E-02 1.000 1.000 0.8000 0.1545E-01 0.1613E-01 0.7915E-02 1.000 1.000 1.0000 0.1545E-01 0.1613E-01 0.7915E-02 1.000 0.1000 0.0000 0.1545E-01 0.1613E-01 0.7915E-02 1.000 0.1000 0.2000 0.1545E-01 0.1613E-01 0.7915E-02 1.000 0.1000 0.4000 0.1545E-01 0.1613E-01 0.7915E-02 1.000 0.1000 0.6000 0.1545E-01 0.1613E-01 0.7915E-02 1.000 0.1000 0.8000 0.1545E-01 0.1613E-01 0.7915E-02 1.000 0.1000 1.0000 0.1545E-01 0.1613E-01 0.7915E-02 1.000 0.1000E-01 0.0000 0.1545E-01 0.1613E-01 0.7915E-02 1.000 0.1000E-01 0.2000 0.1545E-01 0.1613E-01 0.7915E-02 1.000 0.1000E-01 0.4000 0.1545E-01 0.1613E-01 0.7915E-02 1.000 0.1000E-01 0.6000 0.1545E-01 0.1613E-01 0.7915E-02 1.000 0.1000E-01 0.8000 0.1545E-01 0.1613E-01 0.7915E-02 1.000 0.1000E-01 1.0000 0.1545E-01 0.1613E-01 0.7915E-02 1.000 0.1000E-02 0.0000 0.1545E-01 0.1613E-01 0.7915E-02 1.000 0.1000E-02 0.2000 0.1545E-01 0.1613E-01 0.7915E-02 1.000 0.1000E-02 0.4000 0.1545E-01 0.1613E-01 0.7915E-02 1.000 0.1000E-02 0.6000 0.1545E-01 0.1613E-01 0.7915E-02 1.000 0.1000E-02 0.8000 0.1545E-01 0.1613E-01 0.7915E-02 1.000 0.1000E-02 1.0000 0.1545E-01 0.1613E-01 0.7915E-02 navier_stokes_2d_exact_test(): Normal end of execution. 6 October 2025 6:29:11.971 PM