implement dirichlet boundary conditions

This commit is contained in:
Tor Harald Sandve 2022-10-19 09:16:35 +02:00
parent 6f54bd9310
commit 2bd0575523

View File

@ -352,7 +352,7 @@ private:
}
const bool nonzero_massrate = massrate != VectorBlock(0.0);
if (free || nonzero_massrate) {
const auto& exFluidState = problem_().initialFluidState(myIdx);
const auto& exFluidState = problem_().boundaryFluidState(myIdx, dir_id);
BoundaryConditionData bcdata{free ? BCType::FREE : BCType::RATE,
massrate,
exFluidState.pvtRegionIndex(),