diff --git a/modules/navier_stokes/include/fvkernels/INSFVTKESourceSink.h b/modules/navier_stokes/include/fvkernels/INSFVTKESourceSink.h
index f61b494a3a54..c7dabd1d8714 100644
--- a/modules/navier_stokes/include/fvkernels/INSFVTKESourceSink.h
+++ b/modules/navier_stokes/include/fvkernels/INSFVTKESourceSink.h
@@ -72,81 +72,3 @@ class INSFVTKESourceSink : public FVElementalKernel
   std::map<const Elem *, std::vector<const FaceInfo *>> _face_infos;
   ///@}
 };
-
-// //* This file is part of the MOOSE framework
-// //* https://www.mooseframework.org
-// //*
-// //* All rights reserved, see COPYRIGHT for full restrictions
-// //* https://github.com/idaholab/moose/blob/master/COPYRIGHT
-// //*
-// //* Licensed under LGPL 2.1, please see LICENSE for details
-// //* https://www.gnu.org/licenses/lgpl-2.1.html
-
-// #pragma once
-
-// #include "FVElementalKernel.h"
-// #include "INSFVVelocityVariable.h"
-
-// /**
-//  * Computes source the sink terms for the turbulent kinetic energy.
-//  */
-// class INSFVTKESourceSink : public FVElementalKernel
-// {
-// public:
-//   static InputParameters validParams();
-
-//   virtual void initialSetup() override;
-
-//   INSFVTKESourceSink(const InputParameters & parameters);
-
-// protected:
-//   ADReal computeQpResidual() override;
-
-// protected:
-//   /// The dimension of the domain
-//   const unsigned int _dim;
-
-//   /// x-velocity
-//   const Moose::Functor<ADReal> & _u_var;
-//   /// y-velocity
-//   const Moose::Functor<ADReal> * _v_var;
-//   /// z-velocity
-//   const Moose::Functor<ADReal> * _w_var;
-
-//   /// epsilon - dissipation rate of TKE
-//   const Moose::Functor<ADReal> & _epsilon;
-
-//   /// Density
-//   const Moose::Functor<ADReal> & _rho;
-
-//   /// Dynamic viscosity
-//   const Moose::Functor<ADReal> & _mu;
-
-//   /// Turbulent dynamic viscosity
-//   const Moose::Functor<ADReal> & _mu_t;
-
-//   /// Wall boundaries
-//   const std::vector<BoundaryName> & _wall_boundary_names;
-
-//   /// Maximum mixing length allowed for the domain
-//   const Real _max_mixing_length;
-
-//   /// Linearized model?
-//   const bool _linearized_model;
-
-//   /// No equilibrium treatement
-//   const bool _non_equilibrium_treatment;
-
-//   /// C_mu constant
-//   const Real _C_mu;
-
-//   // Production Limiter Constant
-//   const Real _C_pl;
-
-//   ///@{
-//   /// Maps for wall treatement
-//   std::map<const Elem *, bool> _wall_bounded;
-//   std::map<const Elem *, std::vector<Real>> _dist;
-//   std::map<const Elem *, std::vector<const FaceInfo *>> _face_infos;
-//   ///@}
-// };