zhyang-dev/openfoam
1
0
Fork 0
mirror of https://develop.openfoam.com/Development/openfoam.git synced 2025-11-28 03:28:01 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

Removed boundary setting in DsmcCloud field calculation, calling correctBoundaryConditions in solver instead.

Browse Source
This commit is contained in:
graham
2009-03-11 11:35:03 +00:00
parent dc0d41f816
commit 3a73f3d1aa
2 changed files with 10 additions and 34 deletions
Download Patch File Download Diff File Expand all files Collapse all files

11
applications/solvers/dsmc/dsmcFoam/dsmcFoam.C
View File

@ -56,21 +56,30 @@ int main(int argc, char *argv[])
dsmc.evolve();
// Retrieve field data from dsmcCloud
// Retrieve flow field data from dsmcCloud
rhoN = dsmc.rhoN();
rhoN.correctBoundaryConditions();
rhoM = dsmc.rhoM();
rhoM.correctBoundaryConditions();
dsmcRhoN = dsmc.dsmcRhoN();
dsmcRhoN.correctBoundaryConditions();
momentum = dsmc.momentum();
momentum.correctBoundaryConditions();
linearKE = dsmc.linearKE();
linearKE.correctBoundaryConditions();
internalE = dsmc.internalE();
internalE.correctBoundaryConditions();
iDof = dsmc.iDof();
iDof.correctBoundaryConditions();
// Retrieve surface field data from dsmcCloud
q = dsmc.q();

33
src/lagrangian/dsmc/clouds/Templates/DsmcCloud/DsmcCloudI.H
View File

@ -358,10 +358,6 @@ Foam::DsmcCloud<ParcelType>::rhoN() const
rhoN *= nParticle_/mesh().cellVolumes();
trhoN().boundaryField() = trhoN().boundaryField().boundaryInternalField();
trhoN().correctBoundaryConditions();
return trhoN;
}
@ -399,10 +395,6 @@ Foam::DsmcCloud<ParcelType>::rhoM() const
rhoM *= nParticle_/mesh().cellVolumes();
trhoM().boundaryField() = trhoM().boundaryField().boundaryInternalField();
trhoM().correctBoundaryConditions();
return trhoM;
}
@ -438,11 +430,6 @@ Foam::DsmcCloud<ParcelType>::dsmcRhoN() const
dsmcRhoN[cellI]++;
}
tdsmcRhoN().boundaryField() =
tdsmcRhoN().boundaryField().boundaryInternalField();
tdsmcRhoN().correctBoundaryConditions();
return tdsmcRhoN;
}
@ -485,11 +472,6 @@ Foam::DsmcCloud<ParcelType>::momentum() const
momentum *= nParticle_/mesh().cellVolumes();
tmomentum().boundaryField() =
tmomentum().boundaryField().boundaryInternalField();
tmomentum().correctBoundaryConditions();
return tmomentum;
}
@ -527,11 +509,6 @@ Foam::DsmcCloud<ParcelType>::linearKE() const
linearKE *= nParticle_/mesh().cellVolumes();
tlinearKE().boundaryField() =
tlinearKE().boundaryField().boundaryInternalField();
tlinearKE().correctBoundaryConditions();
return tlinearKE;
}
@ -569,11 +546,6 @@ Foam::DsmcCloud<ParcelType>::internalE() const
internalE *= nParticle_/mesh().cellVolumes();
tinternalE().boundaryField() =
tinternalE().boundaryField().boundaryInternalField();
tinternalE().correctBoundaryConditions();
return tinternalE;
}
@ -612,11 +584,6 @@ Foam::DsmcCloud<ParcelType>::iDof() const
iDof *= nParticle_/mesh().cellVolumes();
tiDof().boundaryField() =
tiDof().boundaryField().boundaryInternalField();
tiDof().correctBoundaryConditions();
return tiDof;
}