A bifurcation study of laminar flow in a curved tube of rectangular cross-section

Abstract

The spiralling flow within a curved tube of rectangular cross-section normally forms two distinct cells in the plane of the cross-section, but a transition to a different, four-cell flow is known to occur at a critical value of the axial pressure gradient. This paper shows that, for a square cross-section, the transition is a result of a complex structure of multiple, symmetric and asymmetric solutions. The structure is revealed by solving extended systems of equations for steady-state, fully developed, laminar flow, in a finite-element approximation, to locate exactly the positions of singular points. Continuation methods are used to trace the paths of these points as the aspect ratio and radius of curvature vary. For a square cross-section, multiple, symmetric flows of both two and four cells are found in two distinct ranges of axial pressure gradient q. In addition, a pair of asymmetric solutions arise from a symmetry-breaking bifurcation point on the primary flow branch. The paths of limit points and symmetry-breaking bifurcation points are obtained as the aspect ratio γ varies, and a transcritical bifurcation point is found at γ = 1.43. For larger aspect ratios the secondary four-cell branch is disconnected from that of the primary two-cell flow. The bifurcation set in the two parameters q and γ has a complex structure, with a number of higher-order singularities; in particular, the path of limit points is found to cross in the manner of an unfolded swallowtail catastrophe. A stability analysis shows that all multiple solutions except the two-cell type are unstable with respect to either symmetric or antisymmetric perturbations. For values of the aspect ratio less than 1.43 there is a range of axial pressure gradients for which there is no stable flow. The critical Dean numbers of the singular points are found to vary strongly at small values of radius of curvature β, but the bifurcation set remains qualitatively the same. Previous work is interpreted in the light of the present results.

On the finite Dean problem: linear theory

“…We note that the critical Taylor numbers T c found here are consistent with the results of Winters [38], who solved the full equations numerically. Winters also used bifurcation detection and continuation methods to follow the solution branches as the parameters were changed, and the results of Kachoyan [28] agree qualitatively with the numerical work of Winters.…”

“…This type of flow is intended to model flow in curved rectangular pipes, and flow in helically coiled ducts with small pitch. Recent numerical and experimental work on this problem has been carried out by, for example, Winters [38], Hille et al [25], Joseph et al [27], Ghia and Sokhey [21], De Vriend [18] and Cheng et al [13].…”

“…Hence, any of the flow patterns described here are not necessarlty expected to eventuate, but are demonstrations of the variety of solutions made possible by a restriction of the domain. A weakly-nonlinear analysis will be given in a future paper (Kachoyan [28]) which will investigate the growth of steady solutions and which will include a comparison with the above previous work, especially that of Winters [38].…”

Traveling-wave solutions of the flow in a curved-square duct

“…Although all the solution branches in Fig. 2 have been already found by, for example, Winters, 4 Wang and Yang, 9 and Mondal et al, 5 they are newly calculated for the present study.…”

“…Numerical and experimental investigations 4,5 have shown that many types of steady and unsteady solutions appear as the Dean number, Dn, increases for the curved-rectangular-duct flow, where Dn is proportional to the Reynolds number multiplied by the square root of the duct curvature. If Dn is very small, a steady solution, which has two secondary vortices ͑two-vortex solution͒ in a cross section, is realized in not only the numerical simulations but the experimental observations.…”