In this work, the distributed-order time fractional 2D Sobolev equation is introduced. The orthonormal Bernoulli polynomials, as a renowned family of basis functions, are employed to solve this problem. To effectively use ...