ترغب بنشر مسار تعليمي؟ اضغط هنا

Modal stability analysis of viscoelastic channel and pipe flows using a well-conditioned spectral method

67   0   0.0 ( 0 )
 نشر من قبل Mihailo Jovanovic
 تاريخ النشر 2019
  مجال البحث فيزياء
والبحث باللغة English




اسأل ChatGPT حول البحث

Modal stability analysis provides information about the long-time growth or decay of small-amplitude perturbations around a steady-state solution of a dynamical system. In fluid flows, exponentially growing perturbations can initiate departure from laminar flow and trigger transition to turbulence. Although flow of a Newtonian fluid through a pipe is linearly stable for very large values of the Reynolds number ($Re sim 10^7$), a transition to turbulence often occurs for $Re$ as low as $1500$. When a dilute polymer solution is used in the place of a Newtonian fluid, the transitional value of the Reynolds number decreases even further. Using the spectral collocation method and Oldroyd-B constitutive equation, Garg et al. (Phys. Rev. Lett. 121:024502, 2018) claimed that such a transition in viscoelastic fluids is related to linear instability. Since differential matrices in the collocation method become ill-conditioned when a large number of basis functions is used, we revisit this problem using the well-conditioned spectral integration method. We show modal stability of viscoelastic pipe flow for a broad range of fluid elasticities and polymer concentrations, including cases considered by Garg et al. Similarly, we find that plane Poiseuille flow is linearly stable for cases where Garg et al. report instability. In both channel and pipe flows, we establish the existence of spurious modes that diverge slowly with finer discretization and demonstrate that these can be mistaken for grid-independent modes if the discretization is not fine enough.



قيم البحث

اقرأ أيضاً

Modal and nonmodal analyses of fluid flows provide fundamental insight into the early stages of transition to turbulence. Eigenvalues of the dynamical generator govern temporal growth or decay of individual modes, while singular values of the frequen cy response operator quantify the amplification of disturbances for linearly stable flows. In this paper, we develop well-conditioned ultraspherical and spectral integration methods for frequency response analysis of channel flows of Newtonian and viscoelastic fluids. Even if a discretization method is well-conditioned, we demonstrate that calculations can be erroneous if singular values are computed as the eigenvalues of a cascade connection of the frequency response operator and its adjoint. To address this issue, we utilize a feedback interconnection of the frequency response operator with its adjoint to avoid computation of inverses and facilitate robust singular value decomposition. Specifically, in contrast to conventional spectral collocation methods, the proposed method (i) produces reliable results in channel flows of viscoelastic fluids at high Weissenberg numbers ($sim 500$); and (ii) does not require a staggered grid for the equations in primitive variables.
Newtonian pipe flow is known to be linearly stable at all Reynolds numbers. We report, for the first time, a linear instability of pressure driven pipe flow of a viscoelastic fluid, obeying the Oldroyd-B constitutive equation commonly used to model d ilute polymer solutions. The instability is shown to exist at Reynolds numbers significantly lower than those at which transition to turbulence is typically observed for Newtonian pipe flow. Our results qualitatively explain experimental observations of transition to turbulence in pipe flow of dilute polymer solutions at flow rates where Newtonian turbulence is absent. The instability discussed here should form the first stage in a hitherto unexplored dynamical pathway to turbulence in polymer solutions. An analogous instability exists for plane Poiseuille flow.
The flow of viscoelastic fluids in channels and pipes remain poorly understood, particularly at low Reynolds numbers. Here, we investigate the flow of polymeric solutions in straight channels using pressure measurements and particle tracking. The law of flow resistance is established by measuring the flow friction factor $f_{eta}$ versus flow rate. Two regimes are found: a transitional regime marked by rapid increase in drag, and a turbulent-like regime characterized by a sudden decrease in drag and a weak dependence on flow rate. Lagrangian trajectories show finite transverse modulations not seen in Newtonian fluids. These curvature perturbations far downstream can generate sufficient hoop stresses to sustain the flow instabilities in the parallel shear flow.
Common modal decomposition techniques for flowfield analysis, data-driven modeling and flow control, such as proper orthogonal decomposition (POD) and dynamic mode decomposition (DMD) are usually performed in an Eulerian (fixed) frame of reference wi th snapshots from measurements or evolution equations. The Eulerian description poses some difficulties, however, when the domain or the mesh deforms with time as, for example, in fluid-structure interactions. For such cases, we first formulate a Lagrangian modal analysis (LMA) ansatz by a posteriori transforming the Eulerian flow fields into Lagrangian flow maps through an orientation and measure-preserving domain diffeomorphism. The development is then verified for Lagrangian variants of POD and DMD using direct numerical simulations (DNS) of two canonical flow configurations at Mach 0.5, the lid-driven cavity and flow past a cylinder, representing internal and external flows, respectively, at pre- and post-bifurcation Reynolds numbers. The LMA is demonstrated for several situations encompassing unsteady flow without and with boundary and mesh deformation as well as non-uniform base flows that are steady in Eulerian but not in Lagrangian frames. We show that LMA application to steady nonuniform base flow yields insights into flow stability and post-bifurcation dynamics. LMA naturally leads to Lagrangian coherent flow structures and connections with finite-time Lyapunov exponents (FTLE). We examine the mathematical link between FTLE and LMA by considering a double-gyre flow pattern. Dynamically important flow features in the Lagrangian sense are recovered by performing LMA with forward and backward (adjoint) time procedures.
In this numerical study, an original approach to simulate non-isothermal viscoelastic fluid flows at high Weissenberg numbers is presented. Stable computations over a wide range of Weissenberg numbers are assured by using the root conformation approa ch in a finite volume framework on general unstructured meshes. The numerical stabilization framework is extended to consider thermo-rheological properties in Oldroyd-B type viscoelastic fluids. The temperature dependence of the viscoelastic fluid is modeled with the time-temperature superposition principle. Both Arrhenius and WLF shift factors can be chosen, depending on the flow characteristics. The internal energy balance takes into account both energy and entropy elasticity. Partitioning is achieved by a constant split factor. An analytical solution of the balance equations in planar channel flow is derived to verify the results of the main field variables and to estimate the numerical error. The more complex entry flow of a polyisobutylene-based polymer solution in an axisymmetric 4:1 contraction is studied and compared to experimental data from the literature. We demonstrate the stability of the method in the experimentally relevant range of high Weissenberg numbers. The results at different imposed wall temperatures, as well as Weissenberg numbers, are found to be in good agreement with experimental data. Furthermore, the division between energy and entropy elasticity is investigated in detail with regard to the experimental setup.
التعليقات
جاري جلب التعليقات جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
mircosoft-partner

هل ترغب بارسال اشعارات عن اخر التحديثات في شمرا-اكاديميا