zhyang-dev/OpenFOAM-5.x
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

rhoCentralFoam: Correct handling of single-cell meshes

Browse Source 
Resolves bug-report http://openfoam.org/mantisbt/view.php?id=1548
This commit is contained in:
Henry
2015-03-13 09:51:33 +00:00
parent 43aec8425b
commit fd069939c1
1 changed files with 9 additions and 7 deletions
Download Patch File Download Diff File Expand all files Collapse all files

16
applications/solvers/compressible/rhoCentralFoam/directionInterpolate.H
View File

@ -27,15 +27,17 @@ tmp<GeometricField<Type, fvsPatchField, surfaceMesh> > interpolate
sf.rename(vf.name() + '_' + dir.name());
// Correct BCs of the positive (outgoing) fluxes
if (dir[0] > 0)
forAll(sf.boundaryField(), patchi)
{
forAll(sf.boundaryField(), patchi)
if
(
!sf.boundaryField()[patchi].coupled()
&& sf.boundaryField()[patchi].size()
&& dir.boundaryField()[patchi][0] > 0
)
{
if (!sf.boundaryField()[patchi].coupled())
{
sf.boundaryField()[patchi] =
vf.boundaryField()[patchi].patchInternalField();
}
sf.boundaryField()[patchi] =
vf.boundaryField()[patchi].patchInternalField();
}
}