Crucial properties of the moment closure model FENE-QE

We investigate four crucial properties for testing and evaluating a moment closure approximation of the FENE dumbbell model for dilute polymer solutions: non-negative configuration distribution function, energy dissipation, accuracy of approximation and computational expense. Through mathematical analysis, numerical experiments and comparisons with closure model FENE-P and FENE-YDL, we prove that the FENE-QE approximation has non-negative configuration distribution function, approximates the energy dissipation behavior of original kinetic theory and provides good accuracy. To improve the efficiency of this closure approximation, we introduce a piecewise linear approximation technique that greatly reduces the computational cost. This extension of FENE-QE, FENE-QE-PLA, is the closure model we recommend for simulating dilute polymer solutions.