期刊
JOURNAL OF COMPUTATIONAL PHYSICS
卷 187, 期 1, 页码 1-21出版社
ACADEMIC PRESS INC ELSEVIER SCIENCE
DOI: 10.1016/S0021-9991(03)00028-7
关键词
turbulence; direct numerical simulations; non-newtonian; viscoelastic; FENE; FENE-P
In this paper, we present two new numerical algorithms for updating the equations of motion for a viscoelastic fluid that can be described by the finite extensible nonlinear elastic polymer model with the closure proposed by Peterlin (so called FENE-P model) in a transient calculation. In particular, our algorithms address two difficulties found in earlier formulations. First, the polymer extension, represented by the trace of the conformation tensor, can numerically exceed the finite extensible length causing the restoring spring force to change sign and the calculation to rapidly diverge. In our formulations, we have redefined the conformation tensor so that this possibility no longer exists. Secondly, the conformation tensor must remain symmetric and positive definite at all times for the calculation to remain stable. The accumulation of numerical errors can cause loss of this property, leading to the growth of Hadamard instabilities [J. Non-Newtonian Fluid Mech. 60 (1995) 53]. We present two matrix decompositions that enable us to construct the conformational tensor in a manner that ensures positive definiteness. Numerical tests of the new algorithms show significant departures from other approaches that rely on filtering to remove the instabilities. (C) 2003 Elsevier Science B.V. All rights reserved.
作者
我是这篇论文的作者
点击您的名字以认领此论文并将其添加到您的个人资料中。
推荐
暂无数据