A lack of reliable data treatment method has been for several decades the bottleneck of viscosity measurement by disturbance amplitude damping method of shock waves. In this work the finite difference method is firstly applied to obtain the numerical solutions for disturbance amplitude damping behavior of sinusoidal shock front in inviscid and viscous flow. When water shocked to 15 GPa is taken as an example, the main results are as follows: (1) For inviscid and lower viscous flows the numerical method gives results in good agreement with the analytic solutions under the condition of small disturbance (a 0 /λ=0.02); (2) For the flow of viscosity beyond 200 Pa s (η = κ) the analytic solution is found to overestimate obviously the effects of viscosity. It is attributed to the unreal pre-conditions of analytic solution by Miller and Ahrens; (3) The present numerical method provides an effective tool with more confidence to overcome the bottleneck of data treatment when the effects of higher viscosity in experiments of Sakharov and flyer impact are expected to be analyzed, because it can in principle simulate the development of shock waves in flows with larger disturbance amplitude, higher viscosity, and complicated initial flow. viscosity coefficient, two-dimensional Eulerian flow, shock front, sinusoidal disturbance PACS: 66.20.Cy, 62.50.+p Viscosity of substance at high pressures is important for the description of dynamic evolution in the deeper earth [1,2] and of flow development of shock waves [3,4], which is described by shear and bulk coefficients in Newtonian viscous fluid theory [5]. However, data of these quantities are very lacking for metals and silicates because direct measurements become very difficult at high pressures of Mbar and high temperatures of several thousand Kelvin.The so-called oscillatory damping method was put forward for such an extreme condition, which relies on the correlation between viscosity of shock compressed state and amplitude damping and oscillation of an initial sinusoidal disturbance on shock front in concerned substance. Sakharov and Mineev et al. [4,[6][7][8][9] especially designed experiments for generating shock waves with sinusoidal geometrical disturbance on front and measuring the amplitude damping and oscillation by a combination of explosive loading and high-speed imaging techniques. In their analysis, Zaidel's [10] analytic solution for the disturbance development of shock front were adopted to constrain the shear viscosity coefficients of several materials. Two decades later, Miller and Ahrens [5] re-examined Zaidel's approximations in details. They considered the effect of bulk viscosity, and got an analytic solution for a more general non-uniform initial condition. It should be pointed out that some preconditions were still put on their solution: (a) weak viscosity; (b) small-amplitude disturbance; (c) ideal discontinuity of the shock front; (d) a simple initial flow. These harsh requirements can not be satisfied in actual experiments [4,[6][7][8][9][11][12][1...