CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Programming & Development

Solving stokes flow

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   March 13, 2018, 08:53
Default Solving stokes flow
  #1
Member
 
Pablo Alarcón
Join Date: Mar 2018
Location: Liège
Posts: 53
Rep Power: 4
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

https://openfoamwiki.net/index.php/O...hm_in_OpenFOAM

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’
UEqn.relax();

error: passing ‘const Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::Boundary’ as ‘this’ argument discards qualifiers [-fpermissive]
p.boundaryField().updateCoeffs();

The beginning of my code is:

Code:
#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

Reply

Tags
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 18:36
pimpleDyMFoam computation randomly stops babapeti OpenFOAM Running, Solving & CFD 5 January 24, 2018 06:28
Foam::error::printStack after large bounding epsilon events StoneCFD727 OpenFOAM Running, Solving & CFD 1 January 1, 2018 18:04
High Courant Number @ icoFoam Artex85 OpenFOAM Running, Solving & CFD 11 February 16, 2017 14:40
Cannot run the code properly: very large time step continuity error crst15 OpenFOAM Running, Solving & CFD 9 December 14, 2014 19:17


All times are GMT -4. The time now is 10:04.