CFD Online Logo CFD Online URL
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Programming & Development

Solving stokes flow

Register Blogs Members List Search Today's Posts Mark Forums Read

LinkBack Thread Tools Search this Thread Display Modes
Old   March 13, 2018, 07:53
Default Solving stokes flow
Pablo Alarcón
Join Date: Mar 2018
Location: Liège
Posts: 59
Rep Power: 7
palarcon is on a distinguished road
Good afternoon

I'm a new programmer into OpenFOAM, and as a part of a larger development project, I want to start coding the Stokes flow (meaning steady state, without convection term, no body forces and laminar flow).
To do so, I started with icoFoam.C and then I implemented the SIMPLE algorithm from the description found on the OpenFOAM wiki

I'm not using simpleFoam because I don't need the turbulence either.

While compiling I've got the following two errors, which I think are related to the relaxation needed in the SIMPLE algorithm because icoFoam don't need that, but I have no idea how to fix that.

error: ‘class Foam::tmp<Foam::fvMatrix<Foam::Vector<double> > >’ has no member named ‘relax’

error: passing ‘const Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::Boundary’ as ‘this’ argument discards qualifiers [-fpermissive]

The beginning of my code is:

#include "fvCFD.H"
#include "singlePhaseTransportModel.H"
#include "simpleControl.H"

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

int main(int argc, char *argv[])
#include "setRootCase.H"
#include "createTime.H"
#include "createMesh.H"
#include "createFields.H"
#include "initContinuityErrs.H"

simpleControl simple(mesh);
palarcon is offline   Reply With Quote


icofoam, relax, simple algorithm, stokes flow

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On

Similar Threads
Thread Thread Starter Forum Replies Last Post
Problem with chtMultiregionFoam radiation boundary condition baran_foam OpenFOAM Running, Solving & CFD 10 December 17, 2019 17:36
pimpleDyMFoam computation randomly stops babapeti OpenFOAM Running, Solving & CFD 5 January 24, 2018 05:28
Foam::error::printStack after large bounding epsilon events StoneCFD727 OpenFOAM Running, Solving & CFD 1 January 1, 2018 17:04
High Courant Number @ icoFoam Artex85 OpenFOAM Running, Solving & CFD 11 February 16, 2017 13:40
Cannot run the code properly: very large time step continuity error crst15 OpenFOAM Running, Solving & CFD 9 December 14, 2014 18:17

All times are GMT -4. The time now is 14:16.