The thermal convection of a Jeffreys fluid subjected to a plane Poiseuille flow in a fluid-porous system composed of a fluid layer and a porous layer is studied in the paper.A linear stability analysis and a Chebyshev...The thermal convection of a Jeffreys fluid subjected to a plane Poiseuille flow in a fluid-porous system composed of a fluid layer and a porous layer is studied in the paper.A linear stability analysis and a Chebyshevτ-QZ algorithm are employed to solve the thermal mixed convection.Unlike the case in a single layer,the neutral curves of the two-layer system may be bi-modal in the proper depth ratio of the two layers.We find that the longitudinal rolls(LRs)only depend on the depth ratio.With the existence of the shear flow,the effects of the depth ratio,the Reynolds number,the Prandtl number,the stress relaxation,and strain retardation times on the transverse rolls(TRs)are also studied.Additionally,the thermal instability of the viscoelastic fluid is found to be more unstable than that of the Newtonian fluid in a two-layer system.In contrast to the case for Newtonian fluids,the TRs rather than the LRs may be the preferred mode for the viscoelastic fluids in some cases.展开更多
This study develops a direct optimal growth algorithm for three-dimensional transient growth analysis of perturbations in channel flows which are globally stable but locally unstable. Different from traditional non-mo...This study develops a direct optimal growth algorithm for three-dimensional transient growth analysis of perturbations in channel flows which are globally stable but locally unstable. Different from traditional non-modal methods based on the Orr- Somrnerfeld and Squire (OSS) equations that assume simple base flows, this algorithm can be applied to arbitrarily complex base flows. In the proposed algorithm, a re- orthogonalization Arnoldi method is used to improve orthogonality of the orthogonal basis of the Krylov subspace generated by solving the linearized forward and adjoint Navier-Stokes (N-S) equations. The linearized adjoint N-S equations with the specific boundary conditions for the channel are derived, and a new convergence criterion is pro- posed. The algorithm is then applied to a one-dimensional base flow (the plane Poiseuille flow) and a two-dimensional base flow (the plane Poiseuille flow with a low-speed streak) in a channel. For one-dimensional cases, the effects of the spanwise width of the chan- nel and the Reynolds number on the transient growth of perturbations are studied. For two-dimensional cases, the effect of strength of initial low-speed streak is discussed. The presence of the streak in the plane Poiseuille flow leads to a larger and quicker growth of the perturbations than that in the one-dimensional case. For both cases, the results show that an optimal flow field leading to the largest growth of perturbations is character- ized by high- and low-speed streaks and the corresponding streamwise vortical structures. The lift-up mechanism that induces the transient growth of perturbations is discussed. The performance of the re-orthogonalization Arnoldi technique in the algorithm for both one- and two-dimensional base flows is demonstrated, and the algorithm is validated by comparing the results with those obtained from the OSS equations method and the cross- check method.展开更多
基金Project supported by the National Natural Science Foundation of China(Nos.11702135,11271188,and 11672164)the Natural Science Foundation of Jiangsu Province of China(No.BK20170775)+1 种基金the China Postdoctoral Science Foundation(No.2016M601798)the Jiangsu Planned Project for Postdoctoral Research Funds of China(No.1601169B)。
文摘The thermal convection of a Jeffreys fluid subjected to a plane Poiseuille flow in a fluid-porous system composed of a fluid layer and a porous layer is studied in the paper.A linear stability analysis and a Chebyshevτ-QZ algorithm are employed to solve the thermal mixed convection.Unlike the case in a single layer,the neutral curves of the two-layer system may be bi-modal in the proper depth ratio of the two layers.We find that the longitudinal rolls(LRs)only depend on the depth ratio.With the existence of the shear flow,the effects of the depth ratio,the Reynolds number,the Prandtl number,the stress relaxation,and strain retardation times on the transverse rolls(TRs)are also studied.Additionally,the thermal instability of the viscoelastic fluid is found to be more unstable than that of the Newtonian fluid in a two-layer system.In contrast to the case for Newtonian fluids,the TRs rather than the LRs may be the preferred mode for the viscoelastic fluids in some cases.
基金supported by the National Natural Science Foundation of China(No.11372140)
文摘This study develops a direct optimal growth algorithm for three-dimensional transient growth analysis of perturbations in channel flows which are globally stable but locally unstable. Different from traditional non-modal methods based on the Orr- Somrnerfeld and Squire (OSS) equations that assume simple base flows, this algorithm can be applied to arbitrarily complex base flows. In the proposed algorithm, a re- orthogonalization Arnoldi method is used to improve orthogonality of the orthogonal basis of the Krylov subspace generated by solving the linearized forward and adjoint Navier-Stokes (N-S) equations. The linearized adjoint N-S equations with the specific boundary conditions for the channel are derived, and a new convergence criterion is pro- posed. The algorithm is then applied to a one-dimensional base flow (the plane Poiseuille flow) and a two-dimensional base flow (the plane Poiseuille flow with a low-speed streak) in a channel. For one-dimensional cases, the effects of the spanwise width of the chan- nel and the Reynolds number on the transient growth of perturbations are studied. For two-dimensional cases, the effect of strength of initial low-speed streak is discussed. The presence of the streak in the plane Poiseuille flow leads to a larger and quicker growth of the perturbations than that in the one-dimensional case. For both cases, the results show that an optimal flow field leading to the largest growth of perturbations is character- ized by high- and low-speed streaks and the corresponding streamwise vortical structures. The lift-up mechanism that induces the transient growth of perturbations is discussed. The performance of the re-orthogonalization Arnoldi technique in the algorithm for both one- and two-dimensional base flows is demonstrated, and the algorithm is validated by comparing the results with those obtained from the OSS equations method and the cross- check method.