From 0912c6c7904a886dfbaf7c651a0e23199523f2a5 Mon Sep 17 00:00:00 2001 From: henry Date: Fri, 27 Jun 2008 17:30:28 +0100 Subject: [PATCH] Made consistent with respect to pd BC and ddtPhiCorr --- .../solvers/heatTransfer/buoyantFoam/pEqn.H | 3 --- .../solvers/heatTransfer/buoyantSimpleFoam/pEqn.H | 3 --- .../heatTransfer/buoyantSimpleRadiationFoam/pEqn.H | 3 --- .../heatTransfer/chtMultiRegionFoam/fluid/pEqn.H | 14 ++++---------- 4 files changed, 4 insertions(+), 19 deletions(-) diff --git a/applications/solvers/heatTransfer/buoyantFoam/pEqn.H b/applications/solvers/heatTransfer/buoyantFoam/pEqn.H index 2421bd3bf8..2a1929eb72 100644 --- a/applications/solvers/heatTransfer/buoyantFoam/pEqn.H +++ b/applications/solvers/heatTransfer/buoyantFoam/pEqn.H @@ -1,8 +1,5 @@ bool closedVolume = pd.needReference(); -//pd.boundaryField() == -// p.boundaryField() - rho.boundaryField()*gh.boundaryField() - pRef.value(); - rho = thermo->rho(); volScalarField rUA = 1.0/UEqn.A(); diff --git a/applications/solvers/heatTransfer/buoyantSimpleFoam/pEqn.H b/applications/solvers/heatTransfer/buoyantSimpleFoam/pEqn.H index 8accc61b7f..6f7aa582f0 100644 --- a/applications/solvers/heatTransfer/buoyantSimpleFoam/pEqn.H +++ b/applications/solvers/heatTransfer/buoyantSimpleFoam/pEqn.H @@ -1,6 +1,3 @@ -pd.boundaryField() == - p.boundaryField() - rho.boundaryField()*gh.boundaryField() - pRef.value(); - volScalarField rUA = 1.0/UEqn().A(); U = rUA*UEqn().H(); UEqn.clear(); diff --git a/applications/solvers/heatTransfer/buoyantSimpleRadiationFoam/pEqn.H b/applications/solvers/heatTransfer/buoyantSimpleRadiationFoam/pEqn.H index 8accc61b7f..6f7aa582f0 100644 --- a/applications/solvers/heatTransfer/buoyantSimpleRadiationFoam/pEqn.H +++ b/applications/solvers/heatTransfer/buoyantSimpleRadiationFoam/pEqn.H @@ -1,6 +1,3 @@ -pd.boundaryField() == - p.boundaryField() - rho.boundaryField()*gh.boundaryField() - pRef.value(); - volScalarField rUA = 1.0/UEqn().A(); U = rUA*UEqn().H(); UEqn.clear(); diff --git a/applications/solvers/heatTransfer/chtMultiRegionFoam/fluid/pEqn.H b/applications/solvers/heatTransfer/chtMultiRegionFoam/fluid/pEqn.H index 9e07ab170c..3e071ce8ba 100644 --- a/applications/solvers/heatTransfer/chtMultiRegionFoam/fluid/pEqn.H +++ b/applications/solvers/heatTransfer/chtMultiRegionFoam/fluid/pEqn.H @@ -1,11 +1,6 @@ { bool closedVolume = false; - //pdf[i].boundaryField() == - // thermof[i].p().boundaryField() - // - rhof[i].boundaryField()*ghf[i].boundaryField() - // - pRef.value(); - rhof[i] = thermof[i].rho(); volScalarField rUA = 1.0/UEqn().A(); @@ -15,11 +10,11 @@ fvc::interpolate(rhof[i]) *( (fvc::interpolate(Uf[i]) & fluidRegions[i].Sf()) -// + fvc::ddtPhiCorr(rUA, rhof[i], Uf[i], phif[i]) + + fvc::ddtPhiCorr(rUA, rhof[i], Uf[i], phif[i]) ) - fvc::interpolate(rhof[i]*rUA*ghf[i]) - *fvc::snGrad(rhof[i]) - *fluidRegions[i].magSf(); + *fvc::snGrad(rhof[i]) + *fluidRegions[i].magSf(); // Solve pressure difference # include "pdEqn.H" @@ -55,8 +50,7 @@ initialMassf[i] ) - fvc::domainIntegrate(thermof[i].psi()*thermof[i].p()) - ) - /fvc::domainIntegrate(thermof[i].psi()); + )/fvc::domainIntegrate(thermof[i].psi()); rhof[i] = thermof[i].rho(); }