Topology optimization based on the high-order numerical manifold method by implementing a four-node quadrilateral element