A quasi-static orthogonal method for solving spatial-time neutron kinetics equations in non-steady processes