diff --git a/src/finiteVolume/Make/files b/src/finiteVolume/Make/files index 971740c7da..12a617fc38 100644 --- a/src/finiteVolume/Make/files +++ b/src/finiteVolume/Make/files @@ -388,6 +388,7 @@ $(basicSource)/basicSource/IObasicSourceList.C $(basicSource)/actuationDiskSource/actuationDiskSource.C $(basicSource)/radialActuationDiskSource/radialActuationDiskSource.C $(basicSource)/explicitSource/explicitSource.C +$(basicSource)/explicitSetValue/explicitSetValue.C LIB = $(FOAM_LIBBIN)/libfiniteVolume diff --git a/src/finiteVolume/cfdTools/general/fieldSources/basicSource/actuationDiskSource/actuationDiskSourceTemplates.C b/src/finiteVolume/cfdTools/general/fieldSources/basicSource/actuationDiskSource/actuationDiskSourceTemplates.C index a5c8208d7b..1b1436767e 100644 --- a/src/finiteVolume/cfdTools/general/fieldSources/basicSource/actuationDiskSource/actuationDiskSourceTemplates.C +++ b/src/finiteVolume/cfdTools/general/fieldSources/basicSource/actuationDiskSource/actuationDiskSourceTemplates.C @@ -48,14 +48,14 @@ void Foam::actuationDiskSource::addActuationDiskAxialInertialResistance E.xx() = uniDiskDir.x(); E.yy() = uniDiskDir.y(); E.zz() = uniDiskDir.z(); - const vectorField U1((1.0 - a)*U); + forAll(cells, i) { - T[i] = 2.0*rho[cells[i]]*diskArea_*mag(U1[cells[i]])*a/(1.0 - a); + T[i] = 2.0*rho[cells[i]]*diskArea_*mag(U[cells[i]])*a/(1.0 - a); } forAll(cells, i) { - Usource[cells[i]] += ((Vcells[cells[i]]/V())*T[i]*E) & U1[cells[i]]; + Usource[cells[i]] += ((Vcells[cells[i]]/V())*T[i]*E) & U[cells[i]]; } } diff --git a/src/finiteVolume/cfdTools/general/fieldSources/basicSource/basicSource/basicSource.C b/src/finiteVolume/cfdTools/general/fieldSources/basicSource/basicSource/basicSource.C index 6bca9b3e30..114e6e58cb 100644 --- a/src/finiteVolume/cfdTools/general/fieldSources/basicSource/basicSource/basicSource.C +++ b/src/finiteVolume/cfdTools/general/fieldSources/basicSource/basicSource/basicSource.C @@ -116,10 +116,15 @@ void Foam::basicSource::setCellSet() label globalCellI = returnReduce(cellI, maxOp