0
Research Papers: Techniques and Procedures

Numerical Simulation of Heat Transfer Process of Viscoelastic Fluid Flow at High Weissenberg Number by Log-Conformation Reformulation

[+] Author and Article Information
Hong-Na Zhang

School of Energy Science and Engineering,
Harbin Institute of Technology,
Harbin 150001, China
e-mail: zhanghn@hit.edu.cn

Dong-Yang Li

School of Energy Science and Engineering,
Harbin Institute of Technology,
Harbin 150001, China
e-mail: lidongyang@hit.edu.cn

Xiao-Bin Li

Mem. ASME
School of Energy Science and Engineering,
Harbin Institute of Technology,
Harbin 150001, China
e-mail: lixb@hit.edu.cn

Wei-Hua Cai

Mem. ASME
School of Energy Science and Engineering,
Harbin Institute of Technology,
Harbin 150001, China
e-mail: caiwh@hit.edu.cn

Feng-Chen Li

Mem. ASME
School of Energy Science and Engineering,
Harbin Institute of Technology,
Harbin 150001, China
e-mail: lifch@hit.edu.cn

1Both authors contributed equally to the paper.

2Corresponding authors.

Contributed by the Fluids Engineering Division of ASME for publication in the JOURNAL OF FLUIDS ENGINEERING. Manuscript received November 28, 2016; final manuscript received March 9, 2017; published online June 28, 2017. Assoc. Editor: Mark F. Tachie.

J. Fluids Eng 139(9), 091402 (Jun 28, 2017) (12 pages) Paper No: FE-16-1777; doi: 10.1115/1.4036592 History: Received November 28, 2016; Revised March 09, 2017

Viscoelastic fluids are now becoming promising candidates of microheat exchangers’ working medium due to the occurrence of elastic instability and turbulence at microscale. This paper developed a sound solver for the heat transfer process of viscoelastic fluid flow at high Wi, and this solver can be used to design the multiple heat exchangers with viscoelastic fluids as working medium. The solver validation was conducted by simulating four fundamental benchmarks to assure the reliability of the established solver. After that, the solver was adopted to study the heat transfer process of viscoelastic fluid flow in a curvilinear channel, where apparent heat transfer enhancement (HTE) by viscoelastic fluid was achieved. The observed heat transfer enhancement was attributed to the occurrence of elastic turbulence which continuously mix the hot and cold fluids by the twisting and wiggling flow motions.

FIGURES IN THIS ARTICLE
<>
Copyright © 2017 by ASME
Your Session has timed out. Please sign back in to continue.

References

Figures

Grahic Jump Location
Fig. 7

Time-averaged streamwise velocity component U in the y–z midplane at x = 0 for different Wi: (a) Wi = 0, (b) Wi = 1, (c)Wi = 2, (d) Wi = 5, and (e) Wi = 10

Grahic Jump Location
Fig. 3

(a) Schematic of cross slot geometry, (b) contour of velocity magnitude, and (c) contour of the trace of the conformation tensor together with streamline

Grahic Jump Location
Fig. 4

(a) Schematic of startup of Poiseuille flow; and comparison between the numerical results and analytical solution at (b) E = 10; (c) E = 50 for Oldroyd-B fluid with β = 0.01, U0 is centerline velocity, and t is the dimensionless time

Grahic Jump Location
Fig. 5

Schematic of (a) the computational domain and (b) the mesh of 3D serpentine microchannel. The origin of the coordinate system for this channel is its center point. (Reprinted with permission from Li et al. [51]. Copyright 2017 by Springer.)

Grahic Jump Location
Fig. 6

Temporal evolutions of (a) the streamwise velocity component U, (b) the radial velocity component −V, (c) the dimensionless temperature, at the monitoring point (the center point of the channel) at various Wi, and (d) spectrum of fluctuations of radial velocity component and dimensionless temperature at the monitoring points at Wi = 5

Grahic Jump Location
Fig. 1

Flowchart of the numerical procedures

Grahic Jump Location
Fig. 2

Schematic of sphere sedimentation in a tube

Grahic Jump Location
Fig. 8

Time-averaged heat transfer performance of viscoelastic fluid flow at various Wi

Grahic Jump Location
Fig. 9

Distributions of dimensionless temperature and tr(C) in the horizontal plane (x–y midplane at z = 0) and vertical plane (y–z midplane at x = 0) for Newtonian fluid flow (Wi = 0) and viscoelastic fluid flow (Wi = 10)

Grahic Jump Location
Fig. 10

Numerical results in the horizontal centerline of center cross section for five sets of mesh resolution: (a) time-averaged dimensionless velocity profiles and (b) time-averaged dimensionless temperature profiles, at Re = 1, Pr = 1000, Wi = 20, and β = 0.2. (Reprinted with permission from Li et al. [51]. Copyright 2017 by Springer.)

Tables

Errata

Discussions

Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In