def ns2de_matplotlib_spiral_test ( ): #*****************************************************************************80 # ## NS2DE_MATPLOTLIB_SPIRAL_TEST plots a Spiral Flow velocity field. # # Licensing: # # This code is distributed under the GNU LGPL license. # # Modified: # # 31 January 2015 # # Author: # # John Burkardt # from grid_2d import grid_2d from uvp_spiral import uvp_spiral print '' print 'NS2DE_MATPLOTLIB_SPIRAL_TEST:' print ' Generate a Spiral Flow velocity field on a regular grid.' print ' Display it using MATPLOTLIB' x_lo = 0.0 x_hi = 1.0 x_num = 21 y_lo = 0.0 y_hi = 1.0 y_num = 21 [ x, y ] = grid_2d ( x_num, x_lo, x_hi, y_num, y_lo, y_hi ) nu = 1.0 rho = 1.0 n = x_num * y_num t = 0.0 [ u, v, p ] = uvp_spiral ( nu, rho, n, x, y, t ) header = 'spiral' s = 5.0 ns2de_matplotlib ( header, n, x, y, u, v, s ) print '' print 'NS2DE_MATPLOTLIB_SPIRAL_TEST:' print ' Normal end of execution.' return
def ns2de_gnuplot_spiral_test ( ): #*****************************************************************************80 # ## NS2DE_GNUPLOT_SPIRAL_TEST plots a Spiral Flow velocity field. # # Licensing: # # This code is distributed under the GNU LGPL license. # # Modified: # # 31 January 2015 # # Author: # # John Burkardt # from grid_2d import grid_2d from uvp_spiral import uvp_spiral print '' print 'NS2DE_GNUPLOT_SPIRAL_TEST:' print ' Generate a Spiral Flow velocity field on a regular grid.' print ' Store in GNUPLOT data and command files.' x_lo = 0.0 x_hi = 1.0 x_num = 21 y_lo = 0.0 y_hi = 1.0 y_num = 21 [ x, y ] = grid_2d ( x_num, x_lo, x_hi, y_num, y_lo, y_hi ) nu = 1.0 rho = 1.0 n = x_num * y_num t = 0.0 [ u, v, p ] = uvp_spiral ( nu, rho, n, x, y, t ) header = 'spiral' s = 5.0 ns2de_gnuplot ( header, n, x, y, u, v, s ) print 'NS2DE_GNUPLOT_SPIRAL_TEST:' print ' Normal end of execution.' return
def parameter_spiral_test(): #*****************************************************************************80 # ## PARAMETER_SPIRAL_TEST monitors solution norms over time for various values of NU, RHO. # # Location: # # http://people.sc.fsu.edu/~jburkardt/py_src/navier_stokes_2d_exact/parameter_spiral_test.py # # Licensing: # # This code is distributed under the GNU LGPL license. # # Modified: # # 31 January 2015 # # Author: # # John Burkardt # import numpy as np from r8vec_norm_l2 import r8vec_norm_l2 from r8vec_uniform_ab import r8vec_uniform_ab from uvp_spiral import uvp_spiral print '' print 'PARAMETER_SPIRAL_TEST' print ' Spiral Flow' print ' Monitor solution norms over time for various' print ' values of NU, RHO.' n = 1000 xy_lo = 0.0 xy_hi = 1.0 seed = 123456789 x, seed = r8vec_uniform_ab(n, xy_lo, xy_hi, seed) y, seed = r8vec_uniform_ab(n, xy_lo, xy_hi, seed) # # Vary RHO. # print '' print ' RHO affects the pressure scaling.' print '' print ' RHO NU T ||U|| ||V|| ||P||' print '' nu = 1.0 rho = 1.0 for j in range(0, 3): for k in range(0, 6): t = k / 5.0 u, v, p = uvp_spiral(nu, rho, n, x, y, t) u_norm = r8vec_norm_l2(n, u) / n v_norm = r8vec_norm_l2(n, v) / n p_norm = r8vec_norm_l2(n, p) / n print ' %10.4g %10.4g %8.4g %10.4g %10.4g %10.4g' \ % ( rho, nu, t, u_norm, v_norm, p_norm ) print '' rho = rho / 100.0 # # Vary NU. # print '' print ' NU affects the time scaling.' print '' print ' RHO NU T ||U|| ||V|| ||P||' print '' nu = 1.0 rho = 1.0 for i in range(0, 4): for k in range(0, 6): t = k / 5.0 u, v, p = uvp_spiral(nu, rho, n, x, y, t) u_norm = r8vec_norm_l2(n, u) / n v_norm = r8vec_norm_l2(n, v) / n p_norm = r8vec_norm_l2(n, p) / n print ' %10.4g %10.4g %8.4g %10.4g %10.4g %10.4g' \ % ( rho, nu, t, u_norm, v_norm, p_norm ) print '' nu = nu / 10.0 print '' print 'PARAMETER_TAYLOR_TEST:' print ' Normal end of execution.' return
def parameter_spiral_test ( ): #*****************************************************************************80 # ## PARAMETER_SPIRAL_TEST monitors solution norms over time for various values of NU, RHO. # # Location: # # http://people.sc.fsu.edu/~jburkardt/py_src/navier_stokes_2d_exact/parameter_spiral_test.py # # Licensing: # # This code is distributed under the GNU LGPL license. # # Modified: # # 31 January 2015 # # Author: # # John Burkardt # import numpy as np from r8vec_norm_l2 import r8vec_norm_l2 from r8vec_uniform_ab import r8vec_uniform_ab from uvp_spiral import uvp_spiral print '' print 'PARAMETER_SPIRAL_TEST' print ' Spiral Flow' print ' Monitor solution norms over time for various' print ' values of NU, RHO.' n = 1000 xy_lo = 0.0 xy_hi = 1.0 seed = 123456789 x, seed = r8vec_uniform_ab ( n, xy_lo, xy_hi, seed ) y, seed = r8vec_uniform_ab ( n, xy_lo, xy_hi, seed ) # # Vary RHO. # print '' print ' RHO affects the pressure scaling.' print '' print ' RHO NU T ||U|| ||V|| ||P||' print '' nu = 1.0 rho = 1.0 for j in range ( 0, 3 ): for k in range ( 0, 6 ): t = k / 5.0 u, v, p = uvp_spiral ( nu, rho, n, x, y, t ) u_norm = r8vec_norm_l2 ( n, u ) / n v_norm = r8vec_norm_l2 ( n, v ) / n p_norm = r8vec_norm_l2 ( n, p ) / n print ' %10.4g %10.4g %8.4g %10.4g %10.4g %10.4g' \ % ( rho, nu, t, u_norm, v_norm, p_norm ) print '' rho = rho / 100.0 # # Vary NU. # print '' print ' NU affects the time scaling.' print '' print ' RHO NU T ||U|| ||V|| ||P||' print '' nu = 1.0; rho = 1.0; for i in range ( 0, 4 ): for k in range ( 0, 6 ): t = k / 5.0 u, v, p = uvp_spiral ( nu, rho, n, x, y, t ) u_norm = r8vec_norm_l2 ( n, u ) / n v_norm = r8vec_norm_l2 ( n, v ) / n p_norm = r8vec_norm_l2 ( n, p ) / n print ' %10.4g %10.4g %8.4g %10.4g %10.4g %10.4g' \ % ( rho, nu, t, u_norm, v_norm, p_norm ) print '' nu = nu / 10.0 print '' print 'PARAMETER_TAYLOR_TEST:' print ' Normal end of execution.' return