This paper describes the solvability of Dirichlet problems for Laplace's equation when the boundary data is not smooth enough for the existence of a weak solution in H-1(Omega). Scales of spaces of harmonic functions and of boundary traces are defined and the solutions are characterized as limits of classical harmonic functions in special norms. The generalized harmonic functions, and their norms, are defined using series expansions involving harmonic Steklov eigenfunctions on the domain. It is shown that the usual trace operator has a continuous extension to an isometric isomorphism of specific spaces. This provides a characterization of the generalized solutions of harmonic Dirichlet problems. Numerical simulations of a model problem are described. This problem is related to the dewetting of thin films and the associated phenomenology is described.