摘要

A modified least square finite difference (LSFD) method is proposed to solve the gas film lubrication equation in the head/disk interface in hard disk drives (HDDs). In the modified LSFD scheme, around a reference node there are 24 regular distribution supporting points, and a new scaling matrix is used to replace the one used in the general LSFD method. A least square technique is used to smooth the first-order derivative of flying height where discontinuities occur. Numerical examples show that the modified LSFD method can overcome divergence problem of the general LSFD method in complicated cases and has higher computational efficiency.