|
[Sponsors] |
thermo.rho() -= psi*p_rgh in buoyantPimpleFoam |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
November 11, 2015, 00:18 |
thermo.rho() -= psi*p_rgh in buoyantPimpleFoam
|
#1 |
Member
Yan Wang
Join Date: May 2015
Location: Beijing
Posts: 41
Rep Power: 10 |
Hi everyone,
I already have a good understanding of buoyantSimpleFoam, but in the pEqn.H of buoyantPimpleFoam, the following code makes it difficult to understand the logic of the solver: Code:
// Thermodynamic density needs to be updated by psi*d(p) after the // pressure solution - done in 2 parts. Part 1: thermo.rho() -= psi*p_rgh; volScalarField rAU(1.0/UEqn.A()); surfaceScalarField rAUf("rAUf", fvc::interpolate(rho*rAU)); volVectorField HbyA("HbyA", U); HbyA = rAU*UEqn.H(); surfaceScalarField phig(-rAUf*ghf*fvc::snGrad(rho)*mesh.magSf()); Many thanks, Yan Wang
__________________
Blog: http://blog.sina.com.cn/multiphyzks RG:https://www.researchgate.net/profile/Yan_Wang154 |
|
November 21, 2015, 01:50 |
|
#2 | |
Member
Yan Wang
Join Date: May 2015
Location: Beijing
Posts: 41
Rep Power: 10 |
Quote:
It seems that I have solved the question.
__________________
Blog: http://blog.sina.com.cn/multiphyzks RG:https://www.researchgate.net/profile/Yan_Wang154 |
||
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Fan in buoyantPimpleFoam | tommo97 | OpenFOAM Running, Solving & CFD | 2 | January 12, 2016 16:23 |
buoyantPimpleFoam high Courant numbers | lcbuijs | OpenFOAM Running, Solving & CFD | 5 | March 23, 2015 04:35 |
buoyantPimpleFoam with thermalBaffle1D: strange temperature decrease | donQi | OpenFOAM Running, Solving & CFD | 0 | July 14, 2014 20:30 |
buoyantPimpleFoam Convergence Issues | joel.lehikoinen | OpenFOAM | 1 | December 5, 2013 14:58 |
Questions about buoyantPimpleFoam and rhoPimpleFoam | Mojtaba.a | OpenFOAM Running, Solving & CFD | 6 | August 1, 2012 04:50 |