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')