zhulianhua / dugksFoam

An OpenFOAM solver for Boltzmann model equation using discrete unified gas kinetic scheme

Geek Repo:Geek Repo

Github PK Tool:Github PK Tool

Repeating the calculation?

XcodeRole opened this issue · comments

In the function Foam::discreteVelocity::updateGHsurf() , gSurf_ = (1.0 - relaxFactor)gSurf_ + relaxFactorgEq; has already calculated all the boundary fields (operator overloading calculates both the interior and boundary fields). Why is there a need to explicitly iterate through forAll(gSurfPatch, facei) and calculate the boundary fields again?