Creating solver without pressure corrections

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

 October 1, 2015, 16:38 Creating solver without pressure corrections #1 Member   Lisandro Maders Join Date: Feb 2013 Posts: 50 Rep Power: 4 Hi, I am trying to create a solver that solves only the following equations: 1) Continuity 2) Momentum 3) Generic transport equation of variable Z The fact is that my pressure field is constant, it almost does not vary. And also I take the density from a table based on Z values. So I do not need to make that pressure corrections that SIMPLE algorithm and others do. So how to proper write the .C solver file? Or if you can suggest me to take a look in some specific code that has this simplification it would be good as well. Regards, Lisandro

 October 1, 2015, 20:15 What I have done so far.. #2 Member   Lisandro Maders Join Date: Feb 2013 Posts: 50 Rep Power: 4 Here is what I have done so far, but it is dividing by zero (at least I interpreted the error like this) somewhere in rhoEqn.. Code: ```int main(int argc, char *argv[]) { #include "setRootCase.H" #include "createTime.H" #include "createMesh.H" #include "initFGM.H" #include "createFields.H" simpleControl simple(mesh); // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // #include "CourantNo.H" // #include "compressibleCourantNo.H" while (simple.loop()) { Info<< "Time = " << runTime.timeName() << nl << endl; while (simple.correctNonOrthogonal()) { Info << " Antes de rhoEqn.H " << endl; Info << " Escrevendo rho " << rho << endl; Info << " Escrevendo phi " << phi << endl; #include "rhoEqn.H" #include "UEqn.H" #include "pvEqn.H" #include "ztEqn.H" } if (runTime.write()) { #include "writeFGMfields.H" #include "writeThermoPropertyFields.H" } Info<< "ExecutionTime = " << runTime.elapsedCpuTime() << " s" << " ClockTime = " << runTime.elapsedClockTime() << " s" << nl << endl; } Info<< "End\n" << endl; return 0; }``` Anything incredibly wrong at a first look? Jamamoto takasi likes this.

 October 7, 2015, 09:31 Still weird #3 Member   Lisandro Maders Join Date: Feb 2013 Posts: 50 Rep Power: 4 Jamamoto, I interpreted your "like" as: that looks right! However, the results I am getting for velocity field are very weird, which tells me I missed something in this simplification. Could you or someone else give me some help on this? Is there any solver in OpenFOAM that does not make use of such simplification (no pressure/momentum corrections).. Thanks!

 Tags correction, momentum, pressure, simple, simplification

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post Akshay OpenFOAM Programming & Development 8 October 3, 2015 09:22 user0314 FLUENT 3 September 2, 2011 16:12 aamer Main CFD Forum 2 June 15, 2011 12:45 sankarv OpenFOAM 0 April 4, 2010 18:06 DS & HB Main CFD Forum 0 January 8, 2000 16:00

All times are GMT -4. The time now is 00:24.