from snake.openfoam.simulation import OpenFOAMSimulation
from snake.cuibm.simulation import CuIBMSimulation
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(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='forceCoefficientsCompareKrishnanEtAl2014')
dataframe = simulation.create_dataframe_forces(display_strouhal=True,
display_coefficients=True)
dataframe2 = krishnan.create_dataframe_forces(display_strouhal=True,
display_coefficients=True,
coefficient=2.0)
dataframe = dataframe.append(dataframe2)
print(dataframe)
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)
# author: Olivier Mesnard ([email protected]) # description: Plots the instantaneous force coefficients # and compare to results from previous simulation. # Run this script from the simulation directory. from snake.openfoam.simulation import OpenFOAMSimulation simulation = OpenFOAMSimulation(description='present') 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) other = OpenFOAMSimulation(description='other', directory='') other.read_forces(display_coefficients=True) other.get_mean_forces(limits=[32.0, 64.0]) other.get_strouhal(limits=[32.0, 64.0], order=200) simulation.plot_forces(display_coefficients=True, display_extrema=True, order=200, limits=(0.0, 80.0, 0.0, 3.0), other_simulations=other, other_coefficients=1.0, save_name='forceCoefficientsCompareOther') dataframe = simulation.create_dataframe_forces(display_strouhal=True, display_coefficients=True) dataframe2 = other.create_dataframe_forces(display_strouhal=True, display_coefficients=True) dataframe = dataframe.append(dataframe2) print(dataframe)
# file: plotForceCoefficients.py # author: Olivier Mesnard ([email protected]) # description: Plots the instantaneous force coefficients. # Run this script from the simulation directory. from snake.openfoam.simulation import OpenFOAMSimulation simulation = OpenFOAMSimulation() 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) simulation.plot_forces(display_coefficients=True, display_extrema=True, order=200, limits=(0.0, 80.0, 0.0, 3.0), save_name='forceCoefficients') simulation.plot_forces(display_coefficients=True, display_extrema=True, order=200, limits=(0.0, 100.0, 0.0, 3.0), save_name='forceCoefficientsExtended') dataframe = simulation.create_dataframe_forces(display_strouhal=True, display_coefficients=True) print(dataframe)
"""
Post-processes the force coefficients from a OpenFOAM simulation.
This script reads the forces, computes the mean forces within a given range,
computes the Strouhal number within a range, plots the force coefficients,
saves the figure, and prints a data-frame that contains the mean values.
"""
from snake.openfoam.simulation import OpenFOAMSimulation
simulation = OpenFOAMSimulation()
simulation.read_forces(display_coefficients=True)
time_limits = (32.0, 64.0)
simulation.get_mean_forces(limits=time_limits)
simulation.get_strouhal(limits=time_limits, order=200)
simulation.plot_forces(display_coefficients=True,
display_extrema=True, order=200,
limits=(0.0, 80.0, 0.0, 3.0),
save_name='forceCoefficients',
style='mesnardo',
show=True)
dataframe = simulation.create_dataframe_forces(display_strouhal=True,
display_coefficients=True)
print(dataframe)
"""
Post-processes the force coefficients from a OpenFOAM simulation.
This script reads the forces, computes the mean forces within a given range,
computes the Strouhal number within a range, plots the force coefficients,
saves the figure, and prints a data-frame that contains the mean values.
"""
from snake.openfoam.simulation import OpenFOAMSimulation
simulation = OpenFOAMSimulation()
simulation.read_forces(display_coefficients=True)
time_limits = (32.0, 64.0)
simulation.get_mean_forces(limits=time_limits)
simulation.get_strouhal(limits=time_limits, order=200)
simulation.plot_forces(display_coefficients=True,
display_extrema=True,
order=200,
limits=(0.0, 80.0, 0.0, 3.0),
save_name='forceCoefficients',
style='mesnardo',
show=True)
dataframe = simulation.create_dataframe_forces(display_strouhal=True,
display_coefficients=True)
print(dataframe)
# and compare to results from previous simulation.
# Run this script from the simulation directory.
from snake.openfoam.simulation import OpenFOAMSimulation
simulation = OpenFOAMSimulation(description="present")
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)
other = OpenFOAMSimulation(description="other", directory="")
other.read_forces(display_coefficients=True)
other.get_mean_forces(limits=[32.0, 64.0])
other.get_strouhal(limits=[32.0, 64.0], order=200)
simulation.plot_forces(
display_coefficients=True,
display_extrema=True,
order=200,
limits=(0.0, 80.0, 0.0, 3.0),
other_simulations=other,
other_coefficients=1.0,
save_name="forceCoefficientsCompareOther",
)
dataframe = simulation.create_dataframe_forces(display_strouhal=True, display_coefficients=True)
dataframe2 = other.create_dataframe_forces(display_strouhal=True, display_coefficients=True)
dataframe = dataframe.append(dataframe2)
print(dataframe)
from snake.openfoam.simulation import OpenFOAMSimulation
simulation = OpenFOAMSimulation(description='present')
simulation.read_forces(display_coefficients=True)
time_limits = (32.0, 64.0)
simulation.get_mean_forces(limits=time_limits)
simulation.get_strouhal(limits=time_limits, order=200)
other = Simulation(description='',
directory='',
software='')
other.read_forces()
other.get_mean_forces(limits=time_limits)
other.get_strouhal(limits=time_limits, order=200)
simulation.plot_forces(display_coefficients=True,
display_extrema=True, order=200,
limits=(0.0, 80.0, 0.0, 3.0),
other_simulations=other,
other_coefficients=1.0,
save_name='forceCoefficientsCompareOther',
style='mesnardo',
show=True)
dataframe = simulation.create_dataframe_forces(display_strouhal=True,
display_coefficients=True)
dataframe2 = other.create_dataframe_forces(display_strouhal=True,
display_coefficients=True)
print(dataframe.append(dataframe2))
# file: plotForceCoefficients.py # author: Olivier Mesnard ([email protected]) # description: Plots the instantaneous force coefficients. # Run this script from the simulation directory. from snake.openfoam.simulation import OpenFOAMSimulation simulation = OpenFOAMSimulation() 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) simulation.plot_forces(display_coefficients=True, display_extrema=True, order=200, limits=(0.0, 80.0, 0.0, 3.0), save_name='forceCoefficients') simulation.plot_forces(display_coefficients=True, display_extrema=True, order=200, limits=(0.0, 100.0, 0.0, 3.0), save_name='forceCoefficientsExtended') dataframe = simulation.create_dataframe_forces(display_strouhal=True, display_coefficients=True) print(dataframe)