Transient supercritical droplet evaporation with emphasis on the effects of equation of state

This paper reports a numerical investigation of droplet evaporation in a su percritical environment. A comprehensive physical-numerical model is develo ped to simulate the transcritical and supercritical droplet vaporization ph enomena. The model is based on the time-dependent conservation equations fo r both liquid and gas phases, pressure-dependent variable thermophysical pr operties, and a detailed treatment of liquid-vapor phase equilibrium at the droplet surface. The numerical solution of the two-phase equations employs an arbitrary Eulerian-Lagrangian, explicit-implicit method with a dynamica lly adaptive mesh. The first part of the study examines the capability of d ifferent equations of state (EOS) for predicting the phase equilibrium and transcritical droplet vaporization behavior. Predictions using the Redlich- Kwong (RK) EOS are shown to be markedly different from those using the Peng -Robinson (PR) and Soave-Redlich-Kwong (SRK) EOS. Results for the phase-equ ilibrium of a n-heptane-nitrogen system indicate that compared to PR- and S RK-EOS, the RK-EOS predicts higher fuel vapor concentration, higher liquid- phase solubility of nitrogen, lower critical-mixing-state temperature, and lower enthalpy of vaporization. As a consequence, it significantly overpred icts droplet vaporization rates and, thus, underpredicts droplet lifetimes compared to those predicted by PR- and SRK-EOS, as well as compared to expe rimental data. Furthermore, using RK-EOS, attainment of the critical mixing state at the droplet surface is predicted earlier in droplet lifetime comp ared with that using the other two EOS. In contrast, predictions using the PR-EOS show excellent agreement with experimental data over a wide range of ambient conditions. The PR-EOS is then used for a detailed investigation o f the transcritical droplet vaporization phenomena. Results indicate that a t low to moderate ambient temperatures, the droplet lifetime first increase s and then decreases as the ambient pressure is increased, while at high am bient temperatures, the droplet lifetime decreases monotonically with press ure. This behavior is in accord with the published experimental results. Th e numerical model is also used to obtain the minimum pressure required for the attainment of critical mixing state at the droplet surface for a given ambient temperature. (C) 1999 Elsevier Science Ltd. All rights reserved.