A new method preserving the positive definiteness of a second order tensor variable in flow simulations with application to viscoelastic turbulence.

*(English)*Zbl 1242.76220Summary: We present a new numerical method for the simulation of the flow of complex fluids with internal microstructure that is described by a second order, positive definite, internal structural parameter, so that its positive definiteness is preserved in the simulations. An example of such an internal structural parameter is the conformation tensor, used in the modeling of viscoelastic flows, which characterizes the flow-induced molecular deformation. Numerical methods that guarantee the preservation of the positive definiteness of the conformation tensor help both the physical interpretation of the results and the numerical stability of the simulations. The new method presented here is based on the log-conformation representation of a second order tensor. It is implemented through the application of the Cayley-Hamilton theorem for second order tensors. When necessary, we also use an additional mapping that ensures the boundness of the conformation tensor. This approach has been applied in Direct Numerical Simulations (DNS) of viscoelastic turbulent channel flow, which has received great attention since the mid 1990s with ultimate goal to understand the phenomenon of maximum drag reduction. The algorithm uses a full 3D spectral representation of the spatial dependence for the flow and conformation variables and a second order accurate backward differentiation formula for the integration of the governing equations in time. The key issue for the successful implementation of the proposed scheme is a second order finite difference multigrid diffusion applied to the conformation field. Numerical diffusion has been always added to the hyperbolic evolution equations for the conformation tensor in spectral DNS of viscoelastic turbulent flows in order to remove a-physical high wavenumber instabilities induced due to the chaotic nature of the flow. The main advantage of the finite-difference implementation (as opposed to a spectral implementation as done before) is that it allows for the preservation of the positive definiteness of the conformation tensor.

##### MSC:

76M22 | Spectral methods applied to problems in fluid mechanics |

76A10 | Viscoelastic fluids |

76F10 | Shear flows and turbulence |

PDF
BibTeX
Cite

\textit{K. D. Housiadas} et al., Comput. Fluids 39, No. 2, 225--241 (2010; Zbl 1242.76220)

Full Text:
DOI

##### References:

[1] | Apostolakis, M.V.; Mavrantzas, V.G.; Beris, A.N., Stress gradient-induced migration effects in the taylor – couette flow of a dilute polymer solution, J non-newt fluid mech, 102, 409-445, (2002) · Zbl 1012.76007 |

[2] | Beris, A.N.; Edwards, B.J., Thermodynamics of flowing systems, (1994), Oxford University Press New York · Zbl 0800.76046 |

[3] | Beris, A.N.; Sureshkumar, R., Simulation of time-dependent viscoelastic channel Poiseuille flow at high Reynolds numbers, Chem eng sci, 51, 1451-1471, (1996) |

[4] | Bird, R.B.; Armstrong, R.C.; Hassager, O., Dynamics of polymeric fluids, vol. 1, (1987), Wiley New York |

[5] | Bird, R.B.; Curtiss, C.F.; Amstrong, R.C.; Hassanger, O., Dynamics of polymeric liquids, vol. 2, (1987), Wiley New York |

[6] | Canuto, C.; Hussaini, M.Y.; Quarteroni, A.; Zang, T.A., Spectral methods in fluid dynamics, (1988), Springer Berlin · Zbl 0658.76001 |

[7] | Crochet, M.J.; Legat, V., The consistent streamline-upwind/petrov – galerkin method for viscoelastic flow revisited, J non-newt fluid mech, 42, 283-299, (1992) · Zbl 0825.76055 |

[8] | De Angelis, E.; Casciola, C.M.; Piva, R., DNS of wall turbulence: dilute polymers and self-sustaining mechanisms, Comput fluids, 31, 495-507, (2002) · Zbl 1075.76556 |

[9] | Debae, F.; Legat, V.; Crochet, M.J., Practical evaluation of four mixed finite element methods for viscoelastic flow, J rheol, 38, 421-441, (1994) |

[10] | Deville, M.O.; Fischer, P.F.; Mund, E.H., High-order methods for incompressible fluid flow, (2002), Cambridge University Press Cambridge (MA) · Zbl 1007.76001 |

[11] | Dimitropoulos, C.D.; Sureshkumar, R.; Beris, A.N., Direct numerical simulation of viscoelastic turbulent channel flow exhibiting drag reduction: effect of the variation of rheological parameters, J non-newt fluid mech, 79, 433-468, (1998) · Zbl 0960.76057 |

[12] | Dupret, F.; Marchal, J.M., Loss of evolution in the flow of viscoelastic fluids, J non-newt fluid mech, 20, 143-171, (1986) · Zbl 0629.76008 |

[13] | El-Kareh, A.W.; Leal, L.G., Existence of solutions for all Deborah numbers for a non-Newtonian model modified to include diffusion, J non-newt fluid mech, 33, 257-287, (1989) · Zbl 0679.76004 |

[14] | Fattal, R.; Kupferman, R., Constitutive laws for the matrix-logarithm of the conformation tensor, J non-newt fluid mech, 123, 281-285, (2004) · Zbl 1084.76005 |

[15] | Fattal, R.; Kupferman, R., Time-dependent simulation of viscoelastic flows at high weissenberg number using the log-conformation representation, J non-newt fluid mech, 126, 23-37, (2005) · Zbl 1099.76044 |

[16] | Fox, R.O., Computational models for turbulent reacting flows, (2003), Cambridge University Press New York |

[17] | Housiadas, K.D.; Beris, A.N., Polymer-induced drag reduction: effects of the variations in elasticity and inertia in turbulent viscoelastic channel flow, Phys fluids, 15, 2369-2384, (2003) · Zbl 1186.76235 |

[18] | Housiadas, K.D.; Beris, A.N., An efficient fully implicit spectral scheme for DNS of turbulent viscoelastic channel flow, J non-newt fluid mech, 122, 243-262, (2004) · Zbl 1143.76330 |

[19] | Housiadas, K.D.; Beris, A.N., Characteristic scales and drag reduction evaluation in turbulent channel flow of non-constant viscosity viscoelastic fluids, Phys fluids, 16, 1581-1586, (2004) · Zbl 1186.76236 |

[20] | Housiadas, K.D.; Beris, A.N.; Handler, R.A., Viscoelastic effects on higher order statistics and on coherent structures in turbulent channel flow, Phys fluids, 17, 035106, (2005) · Zbl 1187.76219 |

[21] | Hulsen, M.A., A sufficient condition for the positive definite configuration tensor in differential models, J non-newt fluid mech, 38, 93-100, (1990) |

[22] | Hulsen, M.A.; Fattal, R.; Kupferman, R., Flow of viscoelastic fluids past a cylinder at high weissenberg number: stabilized simulations using matrix algorithms, J non-newt fluid mech, 127, 27-39, (2005) · Zbl 1187.76615 |

[23] | Jeanmart, H.; Winckelmans, G., Investigation of eddy-viscosity models modified using discrete filters: a simplified “regularized variational multiscale model“ and an “enhanced field model”, Phys fluids, 19, 055110, (2007) · Zbl 1146.76426 |

[24] | Joseph, D.D., Fluid dynamics of viscoelastic liquids, (1990), Springer New York · Zbl 0698.76002 |

[25] | Ma, X.; Symeonidis, V.; Karniadakis, G.E., A spectral vanishing viscosity method for stabilizing viscoelastic flows, J non-newt fluid mech, 115, 125-155, (2003) · Zbl 1063.76071 |

[26] | Marchal, J.M.; Crochet, M.J., A new mixed finite element method for calculating viscoelastic flow, J non-newt fluid mech, 26, 77-114, (1987) · Zbl 0637.76009 |

[27] | Min, T.; Choi, H.; Yoo, J.Y., Maximum drag reduction in a turbulent channel flow by polymer additives, J fluid mech, 490, 91-100, (2003) · Zbl 1063.76579 |

[28] | Min, T.; Yoo, J.Y.; Choi, H.; Joseph, D.D., Drag reduction by polymer additives in a turbulent channel flow, J fluid mech, 486, 213-238, (2003) · Zbl 1054.76041 |

[29] | Moin, P.; Kim, J., On the numerical solution of time-dependent viscous incompressible fluid flows involving solid boundaries, J comput phys, 35, 381-392, (1980) · Zbl 0425.76027 |

[30] | Orszag, S.A.; Kells, L.C., Transition to turbulence in plane Poiseuille and plane Couette flow, J fluid mech, 96, 159-205, (1980) · Zbl 0418.76036 |

[31] | Owens, R.G.; Phillips, T.N., Computational rheology, (2002), Imperial College Press London · Zbl 1015.76002 |

[32] | Peyret, R., Spectral methods for incompressible viscous flow, (2002), Springer New York · Zbl 1005.76001 |

[33] | Phillips, T.N.; Soliman, I.M., The influence matrix technique for the numerical spectral simulation of viscous incompressible flows, Numer meth PDE, 7, 9-24, (1991) · Zbl 0712.76071 |

[34] | Ptasinski, P.K.; Boersma, B.J.; Nieuwstadt, F.T.M.; Hulsen, M.A.; Van den Brule, B.H.A.A.; Hunt, J.C.R., Turbulent channel flow near maximum drag reduction: simulations, experiments and mechanisms, J fluid mech, 490, 251-291, (2003) · Zbl 1063.76580 |

[35] | Sagaut, P., Large eddy simulation for incompressible flows: an introduction, (2005), Springer Berlin |

[36] | Sureshkumar, R.; Beris, A.N., Effect of artificial stress diffusivity on the stability of numerical calculations and the flow dynamics of time-dependent viscoelastic flows, J non-newt fluid mech, 6, 53-80, (1995) |

[37] | Sureshkumar, R.; Beris, A.N.; Handler, R.A., Direct numerical simulation of turbulent channel flow of a polymer solution, Phys fluids, 9, 743-755, (1997) |

[38] | Tanner, R.I., Stresses in dilute solutions of bead-nonlinear-spring macromolecules. II. unsteady flows and approximate constitutive relations, Trans soc rheol, 19, 37-65, (1975) · Zbl 0393.76007 |

[39] | Tanner, R.I., Engineering rheology, (2000), Oxford University Press New York · Zbl 0533.76004 |

[40] | Toms BA. Observations on the flow of linear polymer solutions through straight tubes at large Reynolds numbers. In: Proc intl rheological congress, Holland, 1948, vol. 2. North-Holland; 1949, p. 35-141. |

[41] | Trottenberg, U.; Oosterlee, C.W.; Schuller, A., Multigrid, (2001), Academic Press London |

[42] | Tsubokura, M., On the outer large-scale motions of wall turbulence and their interaction with near-wall structures using large eddy simulation, Comput fluids, 38, 37-48, (2009) · Zbl 1237.76059 |

[43] | Vaithianathan, T.; Collins, L.R., Numerical approach to simulating turbulent flow of a viscoelastic polymer solution, J comput phys, 187, 1-21, (2003) · Zbl 1047.76524 |

This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. It attempts to reflect the references listed in the original paper as accurately as possible without claiming the completeness or perfect precision of the matching.