Abstract
We develop an effective numerical scheme to capture hydrodynamic modes in general classical anharmonic chains. This scheme is based on the hydrodynamic theory suggested by Ernst-Hauge-van Leeuwen, which takes full role of pressure fluctuations into account. With this scheme we show that the traditional pictures given by the current nonlinear fluctuating hydrodynamic theory are valid only when the system's pressure is zero and the pressure fluctuations are weak. For nonvanishing pressure, the hydrodynamic modes can, however, respond to small and large pressure fluctuations and relax in some distinct manners. Our results shed new light on understanding thermal transport from the perspective of hydrodynamic theory.