time_limits = (32.0, 64.0) simulation.get_mean_forces(limits=time_limits) simulation.get_strouhal(limits=time_limits, order=200) krishnan = CuIBMSimulation(description='Krishnan et al. (2014)') filepath = os.path.join(os.environ['SNAKE'], 'resources', 'flyingSnake2d_cuibm_anush', 'flyingSnake2dRe2000AoA35', 'forces') krishnan.read_forces(file_path=filepath) krishnan.get_mean_forces(limits=time_limits) krishnan.get_strouhal(limits=time_limits, order=200) simulation.plot_forces(display_coefficients=True, coefficient=2.0, display_extrema=True, order=200, limits=(0.0, 80.0, 0.0, 3.0), other_simulations=krishnan, other_coefficients=2.0, style='mesnardo', save_name='forceCoefficientsCompareKrishnanEtAl2014') dataframe = simulation.create_dataframe_forces(display_strouhal=True, display_coefficients=True, coefficient=2.0) dataframe2 = krishnan.create_dataframe_forces(display_strouhal=True, display_coefficients=True, coefficient=2.0) print(dataframe.append(dataframe2))
simulation = OpenFOAMSimulation(description='IcoFOAM') simulation.read_forces(display_coefficients=True) simulation.get_mean_forces(limits=[32.0, 64.0]) simulation.get_strouhal(limits=[32.0, 64.0], order=200) krishnan = CuIBMSimulation(description='Krishnan et al. (2014)') krishnan.read_forces(file_path='{}/resources/flyingSnake2d_cuibm_anush/' 'flyingSnake2dRe2000AoA35/forces' ''.format(os.environ['SNAKE'])) krishnan.get_mean_forces(limits=[32.0, 64.0]) krishnan.get_strouhal(limits=[32.0, 64.0], order=200) simulation.plot_forces(indices=[1], display_coefficients=True, display_extrema=True, order=200, limits=(0.0, 80.0, 0.0, 3.0), other_simulations=krishnan, other_coefficients=2.0, save_name='liftCoefficientCompareKrishnanEtAl2014') dataframe = simulation.create_dataframe_forces(indices=[1], display_strouhal=True, display_coefficients=True) dataframe2 = krishnan.create_dataframe_forces(indices=[1], display_strouhal=True, display_coefficients=True, coefficient=2.0) dataframe = dataframe.append(dataframe2) print(dataframe)
'revision86-cusp-0.4.0', 'flyingSnake2dRe2000AoA35_20160502', 'numericalSolution')) simulation.read_forces() simulation.get_mean_forces(limits=[32.0, 64.0]) simulation.get_strouhal(limits=[32.0, 64.0], order=200) krishnan = CuIBMSimulation(description='Krishnan et al. (2014)') krishnan.read_forces(file_path='{}/resources/flyingSnake2d_cuibm_anush/' 'flyingSnake2dRe2000AoA35/forces' ''.format(os.environ['SNAKE'])) krishnan.get_mean_forces(limits=[32.0, 64.0]) krishnan.get_strouhal(limits=[32.0, 64.0], order=200) dataframe = simulation.create_dataframe_forces(display_strouhal=True, display_coefficients=True, coefficient=2.0) dataframe2 = krishnan.create_dataframe_forces(display_strouhal=True, display_coefficients=True, coefficient=2.0) dataframe = dataframe.append(dataframe2) print(dataframe) pyplot.style.use( os.path.join(os.environ['SNAKE'], 'snake', 'styles', 'snakeReproducibility.mplstyle')) fig, ax = pyplot.subplots(figsize=(6, 4)) ax.grid(True, zorder=0) ax.set_xlabel('non-dimensional time-unit') ax.set_ylabel('force coefficients') ax.plot(simulation.forces[0].times,
directory=simulation_directory) simulation.read_forces() simulation.get_mean_forces(limits=[32.0, 64.0]) # Computes the mean force coefficients from the cuIBM simulation # with grid-spacing h=0.006 in the uniform region. # The force coefficients are averaged between 32 and 64 time-units. simulation_directory = os.path.join(os.path.dirname(__file__), 'h0.006_vatol16_patol8_dt0.0002') simulation2 = CuIBMSimulation(description='h=0.006', directory=simulation_directory) simulation2.read_forces() simulation2.get_mean_forces(limits=[32.0, 64.0]) # Creates a table with the time-averaged force coefficients. dataframe = simulation.create_dataframe_forces(display_coefficients=True, coefficient=2.0) dataframe2 = simulation2.create_dataframe_forces(display_coefficients=True, coefficient=2.0) dataframe = dataframe.append([dataframe2]) print(dataframe) # Plots the instantaneous force coefficients from the simulation using shifted # markers and from the simulation using exact markers. fig, ax = pyplot.subplots(figsize=(6, 4)) ax.grid(True, zorder=0) ax.set_xlabel('non-dimensional time-unit') ax.set_ylabel('force coefficients') ax.plot(simulation2.forces[0].times, 2.0 * simulation2.forces[0].values, label='$C_d$ - $h=0.006$', color='#377eb8',
simulation.read_forces() simulation.get_mean_forces(limits=[32.0, 64.0]) # Computes the mean force coefficients from the cuIBM simulation # with grid-spacing h=0.006 in the uniform region and atol=1.0E-08 for the # Poisson solver. # The force coefficients are averaged between 32 and 64 time-units. simulation_directory = os.path.join(os.path.dirname(__file__), 'h0.006_vatol16_patol8_dt0.0002') simulation2 = CuIBMSimulation(description='atol=1.0E-08', directory=simulation_directory) simulation2.read_forces() simulation2.get_mean_forces(limits=[32.0, 64.0]) # Creates a table with the time-averaged force coefficients. dataframe = simulation.create_dataframe_forces(display_coefficients=True, coefficient=2.0) dataframe2 = simulation2.create_dataframe_forces(display_coefficients=True, coefficient=2.0) dataframe = dataframe.append([dataframe2]) print(dataframe) # Plots the instantaneous force coefficients from the simulation using shifted # markers and from the simulation using exact markers. fig, ax = pyplot.subplots(figsize=(6, 4)) ax.grid(True, zorder=0) ax.set_xlabel('non-dimensional time-unit') ax.set_ylabel('force coefficients') ax.plot(simulation.forces[0].times, 2.0*simulation.forces[0].values, label='$C_d$ - $atol = 10^{-6}$', color='#377eb8', linestyle='-',
directory=simulation_directory) simulation.read_forces() simulation.get_mean_forces(limits=[32.0, 64.0]) simulation.get_strouhal(limits=[32.0, 64.0], order=200) krishnan = CuIBMSimulation(description='Krishnan et al. (2014)') krishnan.read_forces(file_path=os.path.join(os.environ['SNAKE'], 'resources', 'flyingSnake2d_cuibm_anush', 'flyingSnake2dRe1000AoA35', 'forces')) krishnan.get_mean_forces(limits=[32.0, 64.0]) krishnan.get_strouhal(limits=[32.0, 64.0], order=200) dataframe = simulation.create_dataframe_forces(display_strouhal=True, display_coefficients=True, coefficient=2.0) dataframe2 = krishnan.create_dataframe_forces(display_strouhal=True, display_coefficients=True, coefficient=2.0) dataframe = dataframe.append(dataframe2) print(dataframe) pyplot.style.use(os.path.join(os.environ['SNAKE'], 'snake', 'styles', 'snakeReproducibility.mplstyle')) fig, ax = pyplot.subplots(figsize=(6, 4)) ax.grid(True, zorder=0) ax.set_xlabel('non-dimensional time-unit') ax.set_ylabel('force coefficients')
other.get_strouhal(limits=[32.0, 64.0], order=200) revision86 = CuIBMSimulation(description='cuIBM (old) - CUSP-0.4.0', directory=os.path.join(os.environ['HOME'], 'snakeReproducibilityPackages', 'cuibm', 'revision86', 'Re2000AoA35')) revision86.read_forces(file_path=os.path.join(revision86.directory, 'numericalSolution', 'forces')) revision86.get_mean_forces(limits=[32.0, 64.0]) revision86.get_strouhal(limits=[32.0, 64.0], order=200) dataframe = simulation.create_dataframe_forces(display_strouhal=True, display_coefficients=True, coefficient=2.0) dataframe2 = other.create_dataframe_forces(display_strouhal=True, display_coefficients=True, coefficient=2.0) dataframe3 = revision86.create_dataframe_forces(display_strouhal=True, display_coefficients=True, coefficient=2.0) dataframe = dataframe.append([dataframe2, dataframe3]) print(dataframe) pyplot.style.use(os.path.join(os.environ['SNAKE'], 'snake', 'styles', 'snakeReproducibility.mplstyle')) fig, ax = pyplot.subplots(figsize=(6, 4))
simulation = OpenFOAMSimulation(description='IcoFOAM') simulation.read_forces(display_coefficients=True) simulation.get_mean_forces(limits=[32.0, 64.0]) simulation.get_strouhal(limits=[32.0, 64.0], order=200) krishnan = CuIBMSimulation(description='Krishnan et al. (2014)') krishnan.read_forces(file_path='{}/resources/flyingSnake2d_cuibm_anush/' 'flyingSnake2dRe2000AoA35/forces' ''.format(os.environ['SNAKE'])) krishnan.get_mean_forces(limits=[32.0, 64.0]) krishnan.get_strouhal(limits=[32.0, 64.0], order=200) simulation.plot_forces(indices=[1], display_coefficients=True, display_extrema=True, order=200, limits=(0.0, 80.0, 0.0, 3.0), other_simulations=krishnan, other_coefficients=2.0, save_name='liftCoefficientCompareKrishnanEtAl2014') dataframe = simulation.create_dataframe_forces(indices=[1], display_strouhal=True, display_coefficients=True) dataframe2 = krishnan.create_dataframe_forces(indices=[1], display_strouhal=True, display_coefficients=True, coefficient=2.0) dataframe = dataframe.append(dataframe2) print(dataframe)
'Re2000AoA35')) other.read_forces() other.get_mean_forces(limits=[32.0, 64.0]) other.get_strouhal(limits=[32.0, 64.0], order=200) revision86 = CuIBMSimulation( description='cuIBM (old) - CUSP-0.4.0', directory=os.path.join(os.environ['HOME'], 'snakeReproducibilityPackages', 'cuibm', 'revision86', 'Re2000AoA35')) revision86.read_forces(file_path=os.path.join(revision86.directory, 'numericalSolution', 'forces')) revision86.get_mean_forces(limits=[32.0, 64.0]) revision86.get_strouhal(limits=[32.0, 64.0], order=200) dataframe = simulation.create_dataframe_forces(display_strouhal=True, display_coefficients=True, coefficient=2.0) dataframe2 = other.create_dataframe_forces(display_strouhal=True, display_coefficients=True, coefficient=2.0) dataframe3 = revision86.create_dataframe_forces(display_strouhal=True, display_coefficients=True, coefficient=2.0) dataframe = dataframe.append([dataframe2, dataframe3]) print(dataframe) pyplot.style.use( os.path.join(os.environ['SNAKE'], 'snake', 'styles', 'snakeReproducibility.mplstyle')) fig, ax = pyplot.subplots(figsize=(6, 4)) ax.grid(True, zorder=0)
directory=simulation_directory) simulation.read_forces() simulation.get_mean_forces(limits=[32.0, 64.0]) simulation.get_strouhal(limits=[32.0, 64.0], order=200) krishnan = CuIBMSimulation(description='Krishnan et al. (2014)') krishnan.read_forces(file_path=os.path.join(os.environ['SNAKE'], 'resources', 'flyingSnake2d_cuibm_anush', 'flyingSnake2dRe1000AoA35', 'forces')) krishnan.get_mean_forces(limits=[32.0, 64.0]) krishnan.get_strouhal(limits=[32.0, 64.0], order=200) dataframe = simulation.create_dataframe_forces(display_strouhal=True, display_coefficients=True, coefficient=2.0) dataframe2 = krishnan.create_dataframe_forces(display_strouhal=True, display_coefficients=True, coefficient=2.0) dataframe = dataframe.append(dataframe2) print(dataframe) pyplot.style.use(os.path.join(os.environ['SNAKE'], 'snake', 'styles', 'snakeReproducibility.mplstyle')) fig, ax = pyplot.subplots(figsize=(6, 4)) ax.grid(True, zorder=0) ax.set_xlabel('non-dimensional time-unit') ax.set_ylabel('force coefficients')