Commit 7e25c256 authored by Thomas Fetzer's avatar Thomas Fetzer
Browse files

[exercise-coupling-ff-pm]{turbulence] Improve calculation of velocity gradient for BJS walls

parent b1fa3e7b
......@@ -165,9 +165,6 @@ int main(int argc, char** argv) try
sol[stokesCellCenterIdx] = stokesSol[stokesCellCenterIdx];
sol[stokesFaceIdx] = stokesSol[stokesFaceIdx];
// TODO: update static wall properties
// TODO: update dynamic wall properties
darcyProblem->applyInitialSolution(sol[darcyIdx]);
auto solDarcyOld = sol[darcyIdx];
......@@ -175,6 +172,9 @@ int main(int argc, char** argv) try
couplingManager->init(stokesProblem, darcyProblem, sol);
// TODO: update static wall properties
// TODO: update dynamic wall properties
// the grid variables
using StokesGridVariables = typename GET_PROP_TYPE(StokesTypeTag, GridVariables);
auto stokesGridVariables = std::make_shared<StokesGridVariables>(stokesProblem, stokesFvGridGeometry);
......
......@@ -307,7 +307,7 @@ public:
*/
Scalar permeability(const SubControlVolumeFace& scvf) const
{
return couplingManager().couplingData().darcyPermeability(scvf);
return couplingManager().problem(CouplingManager::darcyIdx).spatialParams().permeabilityAtPos(scvf.center());
}
/*!
......
......@@ -164,12 +164,6 @@ int main(int argc, char** argv) try
auto solStokesOld = stokesSol;
sol[stokesCellCenterIdx] = stokesSol[stokesCellCenterIdx];
sol[stokesFaceIdx] = stokesSol[stokesFaceIdx];
// TODO: update static wall properties
// TODO: update dynamic wall properties
#if EXNUMBER >= 1
stokesProblem->updateStaticWallProperties();
stokesProblem->updateDynamicWallProperties(stokesSol);
#endif
darcyProblem->applyInitialSolution(sol[darcyIdx]);
auto solDarcyOld = sol[darcyIdx];
......@@ -177,6 +171,12 @@ int main(int argc, char** argv) try
auto solOld = sol;
couplingManager->init(stokesProblem, darcyProblem, sol);
// TODO: update static wall properties
// TODO: update dynamic wall properties
#if EXNUMBER >= 1
stokesProblem->updateStaticWallProperties();
stokesProblem->updateDynamicWallProperties(stokesSol);
#endif
// the grid variables
using StokesGridVariables = typename GET_PROP_TYPE(StokesTypeTag, GridVariables);
......
......@@ -338,7 +338,7 @@ public:
*/
Scalar permeability(const SubControlVolumeFace& scvf) const
{
return couplingManager().couplingData().darcyPermeability(scvf);
return couplingManager().problem(CouplingManager::darcyIdx).spatialParams().permeabilityAtPos(scvf.center());
}
/*!
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment