if snake.__version__ != '0.1.2': warnings.warn('The figures were originally created with snake-0.1.2, ' + 'you are using snake-{}'.format(snake.__version__)) directory = os.path.join(os.path.dirname(__file__), 'h0.002') simulation = IBAMRSimulation(description='h=0.002', directory=directory) simulation.read_forces() simulation.get_mean_forces(limits=[32.0, 64.0]) directory = os.path.join(os.path.dirname(__file__), 'h0.004') other = IBAMRSimulation(description='h=0.004', directory=directory) other.read_forces() other.get_mean_forces(limits=[32.0, 64.0]) dataframe = simulation.create_dataframe_forces(display_coefficients=True, coefficient=-2.0) dataframe2 = other.create_dataframe_forces(display_coefficients=True, coefficient=-2.0) dataframe = dataframe.append(dataframe2) print(dataframe) # Plots the instantaneous force coefficients. 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$ (fine - $h=0.002$)', color='#377eb8', linestyle='-',
'Re2000AoA30')) 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', 'flyingSnake2dRe2000AoA30', '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')
directory = os.path.join(os.path.dirname(__file__), 'cfl0.1') simulation = IBAMRSimulation(description='cfl=0.1', directory=directory) simulation.read_forces() simulation.get_mean_forces(limits=[32.0, 64.0]) directory = os.path.join(os.path.dirname(__file__), 'cfl0.3') loose = IBAMRSimulation(description='cfl=0.3', directory=directory) loose.read_forces() loose.get_mean_forces(limits=[32.0, 64.0]) dataframe = simulation.create_dataframe_forces(display_coefficients=True, coefficient=-2.0) dataframe2 = loose.create_dataframe_forces(display_coefficients=True, coefficient=-2.0) dataframe = dataframe.append(dataframe2) print(dataframe) # Plots the instantaneous force coefficients 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$ ($CFL=0.1$)', color='#377eb8', linestyle='-', linewidth=2,
simulation_directory = dirs['ibamr_forceCoefficientsRe2000']['AoA30'] simulation = IBAMRSimulation(description='IBAMR', 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', 'flyingSnake2dRe2000AoA30', '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') ax.plot(simulation.forces[0].times,
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) simulation.plot_forces(indices=[1], 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, save_name='liftCoefficientCompareKrishnanEtAl2014') dataframe = simulation.create_dataframe_forces(indices=[1], display_strouhal=True, display_coefficients=True, coefficient=-2.0) dataframe2 = krishnan.create_dataframe_forces(indices=[1], display_strouhal=True, display_coefficients=True, coefficient=2.0) dataframe = dataframe.append(dataframe2) print(dataframe)