A general analytical solution for the multidimensional transient linear hydraulic diffusivity equation is presented and its application is discussed in this work. The analytical solution is valid for heterogeneous and anisotropic porous media in any orthogonal coordinate system within a regular boundary domain. A source term, dependent on position and time, was included in the partial differential equation to represent production and/or injection wells. A general boundary condition is considered to account for a constant pressure boundary, a no-flow boundary (symmetry or impermeable boundary) or any prescribed flow boundary. The Generalized Integral Transform Technique (GITT) was used to obtain the analytical solution of the proposed problem. This technique transforms a complex partial differential equation into a system of simpler ordinary differential equations that can be solved by analytical, numerical or hybrid methods. Here, we consider a fully analytical solution to the problem. The general solution is then applied to a one-phase tank model in a homogeneous and isotropic porous medium for validation and analysis. A wide variety of applications can be envisioned and discussed for the proposed general solution: well testing in regular and irregular domains, analysis of local reservoir heterogeneities, analysis of reservoir anisotropy degree, analysis of stream lines in naturally fractured and hydraulic fractured reservoirs, formation damage quantification, reservoir shape factor calculations, oil recovery in different waterflood patterns, amongst many others. These applications are discussed and analyzed in the light of the presented general analytical solution.