diff --git a/amr-wind/projection/incflo_apply_nodal_projection.cpp b/amr-wind/projection/incflo_apply_nodal_projection.cpp index e33efd3bd2..ad12059c8c 100644 --- a/amr-wind/projection/incflo_apply_nodal_projection.cpp +++ b/amr-wind/projection/incflo_apply_nodal_projection.cpp @@ -355,12 +355,11 @@ void incflo::ApplyProjection( } // Need to apply custom Neumann funcs for inflow-outflow BC - // after setting the inflow vels above. + // after setting the inflow vels above + // and then enforce solvability by matching outflow to inflow. if (!proj_for_small_dt and !incremental and velocity.has_inout_bndry()) { velocity.apply_bc_funcs(amr_wind::FieldState::New); - } - if (velocity.has_inout_bndry()) { amr_wind::nodal_projection::enforce_inout_solvability( velocity, m_repo.mesh().Geom(), m_repo.num_active_levels()); }