So far, I have been using dirichlet BC’s on no-slip walls to implicitly specify shear stress.
Is it possible to directly specify neumann BC’s on the walls for velocity without specifying the dirichlet?
So far, I have tried manually overriding the total stress term computed in the ghost.mako files with a viscous stress term corresponding to a given streamwise vel. gradient, but that has not yielded correct results in a model.
As an update, I believe I have identified the issue as stemming from insufficient grid resolution. But any verification of my procedure is appreciated.
Edit a flux.mako unique to the BC and specify the velocity gradient where arising in the viscous shear terms.
Edit a ghost.mako unique to the BC where for wall-tangent momentum terms, the inviscid flux is neglected to only enforce a neumann BC.
Specify the “bc_ldg_state” velocities on the ghost side to equal the interior.
An adjustment I have been looking to verify, which I believe is more proper and simple is
2.) Specify the “bc_rsolve_state” velocity components on the ghost-side to correspond to the wall-normal velocity vector.
3.) Specify the “bc_ldg_state” velocity components on the ghost-side to equal the wall-tangent velocity vector.