|
[Sponsors] |
February 11, 2019, 10:32 |
Question on SIMPLEC in OpenFoam
|
#1 |
Senior Member
Join Date: Dec 2017
Posts: 153
Rep Power: 8 |
Dear all,
I have a question on SIMPLEC in openFOAM. Inspecting the source file I can find the following lines in pEqn.H: if (simple.consistent()) { rAtU = 1.0/(1.0/rAU - UEqn.H1()); // ok (rAtU=1/(aP-sum(a_nb))) phiHbyA += fvc::interpolate(rAtU() - rAU)*fvc::snGrad(p)*mesh.magSf(); // Ok, this defines pressure-correction coeffs later. HbyA -= (rAU - rAtU())*fvc::grad(p); // Why doing that ? } The line that creates me some troubles is the last one. It seems that before solving the pressure equation, it corrects the HbyA (which is sum(a_nb)/aP*U_nb, right? ) by a quantity depending on the previous iteration pressure field. Up to my Knowledge, this is not included in the original paper of Van Doormal and Raithby. I naturally not saying that it is wrong, but I am curious to know the theoretical motivation behind this correction. Any Hint? Thank you very much! |
|
February 11, 2019, 19:01 |
|
#2 | |
Senior Member
Andrew Somorjai
Join Date: May 2013
Posts: 175
Rep Power: 13 |
Quote:
http://openfoamwiki.net/index.php/SimpleFoam Best regards! |
||
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
OpenFOAM v3.0+ ?? | SBusch | OpenFOAM | 22 | December 26, 2016 14:24 |
OpenFOAM Training: Programming CFD Course 12-13 and 19-20 April 2016 | cfd.direct | OpenFOAM Announcements from Other Sources | 0 | January 14, 2016 10:19 |
[waves2Foam] Question about the coupling of OceanWave3D and openfoam | app1e | OpenFOAM Community Contributions | 2 | December 29, 2015 04:24 |
A general openfoam development question about boundary condition | kaifu | OpenFOAM | 10 | August 15, 2012 12:51 |
Cross-compiling OpenFOAM 1.7.0 on Linux for Windows 32 and 64bits with Mingw-w64 | wyldckat | OpenFOAM Announcements from Other Sources | 3 | September 8, 2010 06:25 |