The Owls package is a collection of python tools for data analysis and plotting of OpenFOAM cases. It provides the following basic functionality:
- A reader for OpenFOAM scalar and vector fields, set files and log files
- Converting OpenFOAM data to FoamFrames which are derived DataFrames from the pandas library
Owls can be installed via pip
pip install Owls # non development version probably out of date
For a more recent version install the development version
python setup.py install --user
If you have docker installed you can use the automaticly build docker container
docker pull greole/owls
It is recommended to have matplotlib installed, which itself depends on libpng and freetype. The simplest way is to just use your distros package manager to install matplotlib.
import Owls as ow
Owls provides convenient methods to import OpenFOAM data and experimental data i.e.
s1 = read_sets(
folder='/some/path/',
name='foo', # name of the case
)
where folder is the root folder of your OpenFOAM case. Additionally available read_eul
, read_lag
, and read_exp
, read_log
Owls is build on Pandas and you can use the standard features of Pandas like bracket notation to access the data
s1['T'] # access full Temperature row
s1.latest['T'] # access latest temperature row
# access latest temperature data
# at 'y0.1' location
# corresponds to sets/TIME/y0.1_T.xy
s1.latest.location('y0.1')
Data can be filtered using user defined function
s1.filter_locations(ow.isIn('axis')) # filter all location with axis in the name
s1.filter_fields("T", 1000, 2000) # return all fields where temperature is between 1000 and 2000
s1.latest.filter(name="Pos", func=lambda x: 0.1<x<0.2)
Data can be grouped to generate multiple plots at once
s1.by_index('Loc') # group into locations
s1.by_location() # group into locations
s1.by_time() # group by times
Data can be passed to any standard python plotting library i.e. matplotlib, ggplot, seaborn. Currently bokeh is choosen to generate on the fly data visualisations.
If .show is called the natural representation of the FoamFrame will be choosen, i.e. line plots for sets and scatter experimental data. To make standard calls concise the default value of x='Pos'.
s1.by_location().show(y=['v','u'], overlay="Field") # plot u and v for each location
s1.by_location().show(y=['v','u'], overlay="Group") # plot u and v for each location
you can pass bokeh arguments through the plot function for styling
s1.by_index('Loc').show(y='v', y_label='foo')
s1.by_index('Loc').scatter(x='u',y='v',legend='bar')