Artificial compressibility method