A new extended matrix logarithm formulation for the simulation of viscoelastic fluids by spectral elements

This study presents a new algorithm developed in order to remove instabilities observed in the simulation of unsteady viscoelastic fluid flows in the framework of the spectral element method. In this study, we consider a particular model of the finite extensible nonlinear elastic family, FENE-P, but the method could be applied to other differential constitutive equations. Two distinct constraints for the FENE-P equation are imposed: (i) the square of the corresponding finite extensibility parameter of the polymer must be an upper limit for the trace of the conformation tensor and (ii) the eigenvalues of the conformation tensor should remain positive at all steps of the simulation. Negative eigenvalues cause the unbounded growth of instabilities in the flow. The proposed transformation is an extension of the matrix logarithm formulation originally presented by Fattal and Kupferman [1,2]. To evaluate the capability of this new algorithm with the classical conformation tensor, comprehensive studies have been done based on the linear stability analysis to show the influence of this method on the resulting eigenvalue spectra and explain its success to tackle high Weissenberg numbers. With this new method one can tackle high Weissenberg number flow at values of practical interest. A neat improvement of the computational algorithm with stable convergence has been demonstrated in this study. (C) 2010 Elsevier Ltd. All rights reserved.

Time-dependent algorithms for the simulation of viscoelastic flows with spectral element methods: applications and stability

Michel Deville, Nicolas Fietier

This paper presents the development of spectral element methods to simulate unsteady flows of viscoelastic fluids using a closed-form differential constitutive equation. The generation and decay Poiseuille planar flows are considered as benchmark problems to test the abilities of our computational method to deal with truly time-dependent flows. Satisfactory results converging toward steady-state regimes have been obtained for the flow through a four-to-one planar abrupt contraction with unsteady algorithms. Time-dependent simulations of viscoelastic flows are prone to numerical instabilities even for simple geometrical configurations. Possible methods to improve the numerical stability of the computational algorithms are discussed in view of the results carried out with numerical simulations for the flows through a straight channel and the four-to-one contraction.

2003

Spectral element methods for unsteady viscoelastic flows

Michel Deville, Nicolas Fietier

This paper presents the development of spectral element methods to simulate unsteady flows of viscoelastic fluids using a differential constitutive equation. Preliminary results for simple test problems are discussed as well as possible techniques to improve the numerical stability of the computational algorithms.

2000

Simulation of Time-Dependent Viscoelastic Fluid Flows by Spectral Elements

Azadeh Jafari

The research work reported in this dissertation is aimed to develop efficient and stable numerical schemes in order to obtain accurate numerical solution for viscoelastic fluid flows within the spectral element context. The present research consists in the transformation of a large class of differential constitutive models into an equation where the main variable is the logarithm of the conformation tensor or a quantity related to it in a simple way. Particular cases cover the Oldroyd-B fluid and the FENE-P model. Applying matrix logarithm formulation in the framework of the spectral element method is a new type of approach that according to our knowledge no one has implemented before. The reformulation of the classical constitutive equation using a new variable namely the logarithmic formulation, enforces the eigenvalues of the conformation tensor to remain positive for all steps of the simulation. However, satisfying the symmetric positive definiteness of the conformation tensor during the simulation is the necessary condition for stability; but definitely, it is not the sufficient condition to reach meaningful results. The main effort of this research is devoted to introduce a new algorithm in order to overcome the drawback of direct reformulating the classical constitutive equation to the logarithmic one. To evaluate the capability of the extended matrix logarithm formulation, comprehensive studies have been done based on the linear stability analysis to show the influence of this method on the resulting eigenvalue spectra and explain its success to tackle high Weissenberg numbers. With this new method one can treat high Weissenberg number flows at values of practical interest. One of the worst obstacles for numerical simulation of viscoelastic fluids is the presence of spurious modes during the simulation. At high Weissenberg number, many schemes suffer from instabilities and numerical convergence may not be attainable. This is often attributed to the presence of solution singularities due to the geometry, the dominant non-linear terms in the constitutive equations, or the change of type of the underlying mixed-form differential system. Refining the mesh proved to be not very helpful. In this study, to understand more deeply the mechanism of instability generation a comprehensive study about the growth of spurious modes with time evolution, mesh refinement, boundary conditions and Weissenberg number or any other affected parameters has been performed. Then to get rid of these spurious modes the filter based stabilization of spectral element methods proposed by Boyd was applied with success.

EPFL2011

Simulation of Time-Dependent Viscoelastic Fluid Flows by Spectral Elements

Azadeh Jafari

EPFL2011