STEFAN NUMBER-INSENSITIVE NUMERICAL SIMULATION OF THE ENTHALPY METHOD FOR STEFAN PROBLEMS USING THE SPACE-TIME CONSERVATION ELEMENT AND SOLUTION ELEMENT METHOD

摘要

An improvement is introduced to the conservation element and solution element (CE/SE) phase change scheme presented previously. The improvement addresses a well known weakness in numerical simulations of the enthalpy method when the Stefan number, (the ratio of sensible to latent heat) is small (less than 0.1). Behavior of the improved scheme, at the limit of small Stefan numbers, is studied and compared with that of the original scheme. It is shown that high dissipative errors, associated with small Stefan numbers, do not occur using the new scheme.