Translational frames for OpenFOAM

An equivalent of MFR for translational frames for OpenFOAM solvers.

• At the moment instrumentation of solvers is necessary

Copyright information

Copyright (C) 2013-2015
Robert Sawko
Cranfield University, UK

License

This program is free software: you can redistribute it and/or modify it under the terms of the GNU General Public License as published by the Free Software Foundation, either version 3 of the License, or (at your option) any later version.

This program is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details.

You should have received a copy of the GNU General Public License along with this program. If not, see http://www.gnu.org/licenses/.

Compilation

Using GNU g++ should work fine by

# Load OpenFOAM environement
cd <source>
mkdir build #for out of source build
cd build
cmake ../
make #or make -j 4 for a four thread build

On Astral we need to use Intel compilers and old libraries. I created a special CMake variable which will switch between old and new compilation rules and files for different stuff.

cd <source>
mkdir build #for out of source build
cd build
cmake -DCMAKE_CXX_COMPILER=icpc ../
make

Solver instrumentation

Create the translational frame e.g.

autoPtr<translationalFrame> frameOfReference (
    translationalFrame::New(mesh)
);

// This is in order to force the BC change on fixed boundaries
frameOfReference.registerVelocity(U); 

Update the translational frame somewhere in the time or steady state convergence loop. Pressure field and the deviator of the stress tensor have to be passed into the function.

frameOfReference.update(p, turbulence->devReff());