Time domain simulation of a piano. Part 2: Numerical aspects

This article is the second of a series of two papers devoted to the numerical simulation of the piano. It concerns the numerical aspects of the work, the implementation of a piano code and the presentation of corresponding simulations. The main difficulty is time discretization and stability is achieved via energy methods. Numerical illustrations are provided for a realistic piano and compared to experimental recordings.