Taylor-Couette Flow Research Articles

MATRICS: The implicit matrix-free Eulerian hydrodynamics solver

Context. There exists a zoo of different astrophysical fluid dynamics solvers, most of which are based on an explicit formulation and hence stability-limited to small time steps dictated by the Courant number expressing the local speed of sound. With this limitation, the modeling of low-Mach-number flows requires small time steps that introduce significant numerical diffusion, and a large amount of computational resources are needed. On the other hand, implicit methods are often developed to exclusively model the fully incompressible or 1D case since they require the construction and solution of one or more large (non)linear systems per time step, for which direct matrix inversion procedures become unacceptably slow in two or more dimensions. Aims. In this work, we present a globally implicit 3D axisymmetric Eulerian solver for the compressible Navier–Stokes equations including the energy equation using conservative formulation and a fully simultaneous approach. We use the second-order-in-time backward differentiation formula for temporal discretization as well as the κ scheme for spatial discretization. We implement different limiter functions to prohibit the occurrence of spurious oscillations in the vicinity of discontinuities. Our method resembles the well-known monotone upwind scheme for conservation laws (MUSCL). We briefly present efficient solution methods for the arising sparse and nonlinear system of equations. Methods. To deal with the nonlinearity of the Navier–Stokes equations we used a Newton iteration procedure in which the required Jacobian matrix-vector product was reconstructed with a first-order finite difference approximation to machine precision in a matrix-free way. The resulting linear system was solved either completely matrix-free with a combination of a sufficient Krylov solver and an approximate Jacobian preconditioner or semi-matrix-free with the incomplete lower upper factorization technique as a preconditioner. The latter was dependent on a standalone approximation of the Jacobian matrix, which was optionally calculated and needed solely for the purpose of preconditioning. Results. We show our method to be capable of damping sound waves and resolving shocks even at Courant numbers larger than one. Furthermore, we prove the method’s ability to solve boundary value problems like the cylindrical Taylor-Couette flow (TC), including viscosity, and to model transition flows. To show the latter, we recover predicted growth rates for the vertical shear instability, while choosing a time step orders of magnitude larger than the explicit one. Finally, we verify that our method is second order in space by simulating a simplistic, stationary solar wind.

Nonlinear Asymptotic Stability and Transition Threshold for 2D Taylor–Couette Flows in Sobolev Spaces

In this paper, we investigate the stability of the 2-dimensional (2D) Taylor–Couette (TC) flow for the incompressible Navier–Stokes equations. The explicit form of velocity for 2D TC flow is given by u=(Ar+Br)(-sinθ,cosθ)T\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$u=(Ar+\\frac{B}{r})(-\\sin \ heta , \\cos \ heta )^T$$\\end{document} with (r,θ)∈[1,R]×S1\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$(r, \ heta )\\in [1, R]\ imes \\mathbb {S}^1$$\\end{document} being an annulus and A, B being constants. Here, A, B encode the rotational effect and R is the ratio of the outer and inner radii of the annular region. Our focus is the long-term behavior of solutions around the steady 2D TC flow. While the laminar solution is known to be a global attractor for 2D channel flows and plane flows, it is unclear whether this is still true for rotating flows with curved geometries. In this article, we prove that the 2D Taylor–Couette flow is asymptotically stable, even at high Reynolds number (Re∼ν-1\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$Re\\sim \ u ^{-1}$$\\end{document}), with a sharp exponential decay rate of exp(-ν13|B|23R-2t)\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\exp (-\ u ^{\\frac{1}{3}}|B|^{\\frac{2}{3}}R^{-2}t)$$\\end{document} as long as the initial perturbation is less than or equal to ν12|B|12R-2\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\ u ^\\frac{1}{2} |B|^{\\frac{1}{2}}R^{-2}$$\\end{document} in Sobolev space. The powers of ν\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\ u $$\\end{document} and B in this decay estimate are optimal. It is derived using the method of resolvent estimates and is commonly recognized as the enhanced dissipative effect. Compared to the Couette flow, the enhanced dissipation of the rotating Taylor–Couette flow not only depends on the Reynolds number but also reflects the rotational aspect via the rotational coefficient B. The larger the |B|, the faster the long-time dissipation takes effect. We also conduct space-time estimates describing inviscid-damping mechanism in our proof. To obtain these inviscid-damping estimates, we find and construct a new set of explicit orthonormal basis of the weighted eigenfunctions for the Laplace operators corresponding to the circular flows. These provide new insights into the mathematical understanding of the 2D Taylor–Couette flows.

Hydrodynamics-engineered polyaniline nanofibers on graphene nanosheets for high-performance pH sensors

Real-time sweat monitoring in the human body using wearable sensors is promising for noninvasive and continuous healthcare applications. For reliable and accurate sweat analysis, further improvement in the response time and long-term stability of sensors is required. This paper reports an electrochemical pH sensor based on polyaniline/graphene (PANi-Gr) with rapid electron and ion transfer. The implementation of the Taylor–Couette flow approach allowed the efficient exfoliation of graphite into Gr nanosheets and the rapid polymerization of aniline monomers on the Gr surface, resulting in nanofiber PANi-deposited Gr nanocomposites. PANi-Gr had a large surface area and porous structure, which are favorable for rapid and stable electrochemical reactions. The screen-printed PANi-Gr-based pH sensor exhibited a Nernstian sensitivity of 61.91 mV/pH, rapid response time of <1 s, and long-term stability upon monitoring with minimal potential drift (0.04 mV h−1) under strong acidic and basic conditions. Furthermore, the on-body performance of the wearable pH sensor integrated with a wireless electronic device was analyzed via the noninvasive real-time monitoring of human sweat pH during stationary cycling. The developed sensor can be used in wearable real-time sweat monitoring devices for medical and healthcare applications.

Process Intensification of the Continuous Synthesis of Bio-Derived Monomers for Sustainable Coatings Using a Taylor Vortex Flow Reactor.

We describe the optimization and scale-up of two consecutive reaction steps in the synthesis of bio-derived alkoxybutenolide monomers that have been reported as potential replacements for acrylate-based coatings (Sci. Adv.2020, 6, eabe0026). These monomers are synthesized by (i) oxidation of furfural with photogenerated singlet oxygen followed by (ii) thermal condensation of the desired 5-hydroxyfuranone intermediate product with an alcohol, a step which until now has involved a lengthy batch reaction. The two steps have been successfully telescoped into a single kilogram-scale process without any need to isolate the 5-hydroxyfuranone between the steps. Our process development involved FTIR reaction monitoring, FTIR data analysis via 2D visualization, and two different photoreactors: (i) a semicontinuous photoreactor based on a modified rotary evaporator, where FTIR and 2D correlation spectroscopy (2D-COS) revealed the loss of the methyl formate coproduct, and (ii) our fully continuous Taylor Vortex photoreactor, which enhanced the mass transfer and permitted the use of near-stoichiometric equivalents of O2. The use of in-line FTIR monitoring and modeling greatly accelerated process optimization in the Vortex reactor. This led to scale-up of the photo-oxidation in 85% yield with a projected productivity of 1.3 kg day-1 and a space-time yield of 0.06 mol day-1 mL-1. Higher productivities could be achieved while sacrificing yield (e.g., 4 kg day-1 at 40% yield). The use of superheated methanol at 200 °C in a pressurized thermal flow reactor accelerated the second step, the thermal condensation of 5-hydroxyfuranone, from a 20 h batch reflux reaction (0.5 L, 85 g) to a space time of <1 min in a reactor only 3 mL in volume operating with projected productivities of >700 g day-1. Proof of concept for telescoping the two steps was established with an overall two-step yield of 67%, producing a process with a projected productivity of 1.1 kg day-1 for the methoxybutenolide monomer without any purification of the 5-hydroxyfuranone intermediate.

Coherent structures of elastoinertial instabilities in Taylor–Couette flows

We combine flow visualisation techniques and particle image velocimetry to experimentally investigate the higher-order transition to elastoinertial turbulence of Boger fluids ( $El = 0.11\unicode{x2013}0.34$ ) in Taylor–Couette flows. The observed route to turbulence is associated with the appearance of chaotic inflow jets, termed flame patterns, for increasing inertia, and stable structures of solitons, known as diwhirls, for decreasing inertia. We also report for the first time spatially and temporally resolved flow fields in the meridional plane for the three characteristic viscoelastic flow regimes (diwhirls, flame patterns and elastoinertial turbulence). We observe in all cases coherent structures of dynamically independent solitary vortex pairs. The stability of these coherent structures is jet-dominated and can be mainly ascribed to the high extension of the polymer chains in the inflow boundaries in the $r$ – $z$ plane. Solitary pairs are self-sustained when created through random events and do not split; instead, they merge when moving sufficiently close and annihilate when hoop stresses are not sufficient to sustain them. The highly localised and random events result in highly fluctuating, chaotic flow states. We estimate the decay exponent of spatial power spectral density, illustrating a universal scaling of $-2.5$ for elastoinertial turbulence. Based on our observations and in an effort to unify and combine precedent theories with our results, we suggest a mechanism for the origins of elastoinertial instabilities, accounting for both the effect of elasticity on the vortex formation and the effect of increasing/decreasing inertia on the flow dynamics.

Controlling Couette flow by alternating axial mass flux

This paper presents numerical simulations of the Taylor vortex flow under the influence of an externally applied alternating axial mass flux (through-flow) in a Taylor–Couette system with axial periodic boundary conditions. Such an axially modulating flow can lead to a significant variation in the onset of primary instabilities. Depending on the system parameters, the effect can be both stabilizing and destabilizing, i.e., shifting the bifurcation threshold to larger or smaller control parameters, respectively. It is found that the system response around the primary instability is sensitive to and critically influenced by an alternating mass flux, particularly the modulation frequency. We show that such an alternating axial flow represents an easily and, more importantly, precisely controllable key parameter to change the nonlinear system response from subcritical to supercritical behavior and vice versa. Furthermore, we observe different parameter regimes with regular and irregular intermittent flow dynamics.

Energy Budget Characterisation of the Optimal Disturbance in Stratified Shear Flow

Stratified Taylor–Couette flow (STCF) undergoes transient growth. Recent studies have shown that there exists transient amplification in the linear regime of counter-rotating STCF. The kinetic budget of the optimal transient perturbation is analysed numerically to simulate the interaction of the shear production (SP), buoyancy flux (BP), and other energy components that contributes to the total optimal transient kinetic energy. These contributions affect the total energy by influencing the perturbation to extract kinetic energy (KE) from the mean flow. The decay of the amplification factor resulted from the positive amplification of both BP and SP, while the growth is attributed to the negative and positive amplification of BP and SP, respectively. The optimal SP is positively amplified, implying that there is the possibility of constant linear growth. These findings agree with the linear growth rate for increasing values of Grashof number.

Evaluation of Droplets Breakup Models of Emulsification Process in a Liquid–Liquid Taylor–Couette Flow

Evaluation of Droplets Breakup Models of Emulsification Process in a Liquid–Liquid Taylor–Couette Flow

Process intensification of a catalytic-wall Taylor-Couette reactor through unconventional modulation of its angular speed

This study proposes a non-conventional operation of a Taylor-Couette reactor by introducing a periodic flow perturbation as a modulated angular velocity of its inner cylinder. The applicability of this type of waveform has yet to be explored in Taylor-Couette reactors as a means to enhance mass transfer phenomena in reactive systems. A multiphysics numerical study was carried out considering a mass-transfer limited system with a catalytic reaction at the outer cylinder boundary of the reactor while applying different modulations of the inner cylinder angular speed. Results showed that a modular signal could yield conversions similar to a constant angular-speed operation. However, the use of modulating signals brings two essential benefits. First, it enhances mass transfer, which yields higher conversions by dynamically changing the flow patterns in the reactor. This improvement is demonstrated and discussed in terms of a dynamic mixing index, which accounts for the formation and abrupt disruption of Taylor vortices in the reactor. Second and more importantly, this type of operation leads to an overall reduction in the average electrical power required to drive the system (∼25 % reduction). The present study opens the possibility of using intelligent control strategies to optimize reactions and intensify conventional systems with non-conventional operation modes.

Universal relative scaling of longitudinal structure functions in shear-dominated turbulence

Shear significantly influences turbulence in the energy-containing range of shear-dominated flows, and the longitudinal structure functions do not have a universal form as they do in homogeneous isotropic turbulence. Despite this, the relative scaling of structure functions exhibits universal sub-Gaussian behaviour in shear-dominated flows, in particular for turbulent boundary layers, channels and Taylor–Couette flows. Our investigation of a turbulent vertical buoyancy layer at $Pr = 0.71$ using direct numerical simulation shows this universality even in moderate-Reynolds-number buoyancy-driven but shear-dominated boundary layers. It is demonstrated that the universality is related to the energy density of the eddies, which attains a hierarchical equilibrium in the energy-containing range of shear-dominated turbulence. We conjecture that the universal sub-Gaussian behaviour of the energy density of the energy-containing range, which was considered to be non-trivial in prior studies, is related to the universal anomalous scaling exponents of the inertial subrange turbulence. Based on this conjecture, we propose a hypothesis that relates large-scale eddies and the intermittent dissipation field in shear-dominated turbulence, highlighting a relationship between large and small scales. A phenomenological model is also developed to predict the scaling, which is verified using data from a turbulent boundary layer, half-channel and vertical buoyancy layer at friction Reynolds numbers spanning four orders of magnitude. Excellent agreement is observed.

Water Impalement Resistance and Drag Reduction of the Superhydrophobic Surface with Hydrophilic Strips.

Superhydrophobic surfaces (SHS) offer versatile applications by trapping an air layer within microstructures, while water jet impact can destabilize this air layer and deactivate the functions of the SHS. The current work presents for the first time that introducing parallel hydrophilic strips to SHS (SHS-s) can simultaneously improve both water impalement resistance and drag reduction (DR). Compared with SHS, SHS-s demonstrates a 125% increase in the enduring time against the impact of water jet with velocity of 11.9 m/s and a 97% improvement in DR at a Reynolds number of 1.4 × 104. The key mechanism lies in the enhanced stability of the air layer due to air confinement by the adjacent three-phase contact lines. These lines not only impede air drainage through the surface microstructures during water jet impact, entrapping the air layer to resist water impalement, but also prevent air floating up due to buoyancy in Taylor-Couette flow, ensuring an even spread of the air layer all over the rotor, boosting DR. Moreover, failure modes of SHS under water jet impact are revealed to be related to air layer decay and surface structure destruction. This mass-producible structured surface holds the potential for widespread use in DR for hulls, autonomous underwater vehicles, and submarines.

Flow transitions and effective properties in multiphase Taylor–Couette flow

The properties of multiphase flows are challenging to measure, and yet effective properties are fundamental to modelling and predicting flow behaviour. The current study is motivated by rheometric measurements of a gas-fluidized bed using a coaxial rheometer in which the fluidization rate and the rotational speed can be varied independently. The measured torque displays a range of rheological states: quasistatic, dense granular flow behaviour at low fluidization rates and low-to-moderate shear rates; turbulent toroidal-vortex flow at high shear rates and moderate-to-high fluidization rates; and viscous-like behaviour with rate-dependent torque at high shear rates and low fluidization or at low shear rates and high fluidization. To understand the solid-like to fluid-like transitions, additional experiments were performed in the same rheometer using single-phase liquid and liquid–solid suspensions. The fluidized bed experiments are modelled as a Bingham plastic for low fluidization rates, and as a shear-thinning Carreau liquid at high fluidization rates. The suspensions are modelled using the Krieger–Dougherty effective viscosity. The results demonstrate that, by using the effective properties, the inverse Bingham number marks the transition from solid-like to viscous-flow behaviour; a modified gap Reynolds number based on the thickness of the shear layer specifies the transition from solid-like to turbulent vortical flow; and a gap Reynolds number distinguishes viscous behaviour from turbulent vortical flow. The results further demonstrate that these different multiphase flows undergo analogous flow transitions at similar Bingham or Reynolds numbers and the corresponding dimensionless torques show comparable scaling in response to annular shear.

Turbulent mixed convection in a horizontal cylindrical cavity with the off-lattice Boltzmann method

Large-eddy simulations (LES) of turbulent mixed convection in a horizontal concentric cylindrical cavity when cylinders are subjected to rotations are investigated in the present study. A characteristic-based finite-difference off-lattice Boltzmann method solver is employed to investigate whether rotation enhances heat transfer in a horizontal concentric cylindrical cavity configuration where buoyancy acts perpendicular to the cylinder axis. The study focuses on assessing heat transfer characteristics for different rotation intensities in both counter-rotating-cylinders (CRC) and outer-rotating-cylinder (ORC) configurations. The LES simulations are performed for Rayleigh number, Ra=106, and Reynolds numbers, Re=0,500,1000, and 5000 for air (Pr=0.71), employing a classical Smagorinsky sub-grid-scale (SGS) model. The numerical solver is verified using two benchmark problems, namely (i) Buoyancy-driven flow in a differentially heated air-filled cubic enclosure and (ii) turbulent Taylor-Couette flow. The CRC and ORC configurations are chosen for the three-dimensional simulations in this study, prompted by the notable improvement in heat transfer observed in their respective two-dimensional simulations. A thorough comparison between two- and three-dimensional mixed convection results reveals substantial influence from secondary flow in the axial direction. While three-dimensional CRC simulations exhibit a modest increase in heat transfer at Re=500, the ORC configuration experiences a decrease in heat transfer rate as the Reynolds number increases, contrasting with two-dimensional simulation trends. The CRC configuration enhances the axial motion of the flow with Reynolds number, and its magnitude is larger than that of the ORC configuration and natural convection. A significant difference between two-dimensional and three-dimensional mixed convection simulation results is observed because of considerable secondary flow. Furthermore, the flow and thermal characteristics of turbulent mixed convection are investigated using streamlines, isotherms, Q-criterion, Nusselt number variation along the cylinders, anisotropy analysis, turbulent kinetic energy (TKE), dissipation, contours of radial and axial components of velocity.

Synthesis of Monodisperse Silica Nanospheres and Titania Beads from Rotating Cylinder System

Economic raw materials such as low-grade ethanol and Tetraethyl Ortosilicate (TEOS) were adopted as reactants for the synthesis of monodisperse silica nanospheres using a Taylor Vortex flow reactor with a rotating cylinder system. Factors affecting the mean diameter and size distribution of the silica grains were investigated by changing the rotating speed of the inner cylinder and the composition of reactants. In the case of titania beads, grain diameter size was adjusted by changing the amount of precursor (titanium (IV) isopropoxide, TIP), hexadecylamine (HDA) or Potassium Chloride (KCl) and reaction temperature as well as the rotation speed of the cylinder. As a demonstrative application, titania beads could be used for the removal of organic dye as a model contaminant dissolved in an aqueous medium using a photocatalytic decomposition system.

Taylor–Couette flow and heat transfer in an elliptical enclosure with a rotating inner cylinder

We numerically investigate Taylor–Couette flows within a system consisting of an elliptical outer cylinder and a rotating inner circular cylinder, with particular emphasis on the behavior of Taylor cells. The three-dimensional unsteady Navier–Stokes equations are solved under the assumption of axial periodicity. Also, a scalar transport equation is solved for the heat transfer. Our methodology employs a Fourier-spectral meshless discretization technique, which interpolates variables at scattered points using polyharmonic splines and appended polynomials. A pressure-projection algorithm achieves the time advancement of the flow equations. We present findings for an elliptical enclosure with an aspect ratio of two, examining a range of Reynolds numbers (Re) from subcritical to 300. Our analysis includes streamlines, axial velocity contours, pressure, vorticity, and temperature profiles. The results indicate that the flow remains steady up to Re≈300 before transitioning to an unsteady state at Re≈350.

Variational construction of tubular and toroidal streamsurfaces for flow visualization

Approximate streamsurfaces of a three-dimensional velocity field have recently been constructed as isosurfaces of the closest first integral of the velocity field. Such approximate streamsurfaces enable effective and efficient visualization of vortical regions in three-dimensional flows. Here we propose a variational construction of these approximate streamsurfaces to remove the limitation of Fourier series representation of the first integral in earlier work. Specifically, we use finite-element methods to solve a partial differential equation that describes the best approximate first integral for a given velocity field. We use several examples to demonstrate the power of our approach for three-dimensional flows in domains with arbitrary geometries and boundary conditions. These include generalized axisymmetric flows in the domains of a sphere (spherical vortex), a cylinder (cylindrical vortex) and a hollow cylinder (Taylor–Couette flow) as benchmark studies for various computational domains, non-integrable periodic flows (ABC and Euler flows) and Rayleigh–Bénard convection flows. We also illustrate the use of the variational construction in extracting momentum barriers in Rayleigh–Bénard convection.

Thermo-hydrodynamic lubrication and energy dissipation mechanism of a pump-turbine thrust bearing in load-rejection process

The dynamic behavior of the pump-turbine thrust bearing is important to the safety operation of the unit. This paper analyzed the lubrication and energy dissipation mechanism of pump-turbine thrust bearing during load-rejection based on the thermo-hydrodynamic model. The results show that the variation of the axial force and the maximum pad surface pressure is basically consistent with that of the inclination of the pad surface. The change of the friction loss is consistent with the change of rotational speed, while the change of pad surface temperature is affected by the combination of pad inclination and rotational speed. The chaotic flow in the oil tank is accompanied by different forms of vortices, such as Taylor vortices, vortex pairs, and Karman vortices, and results in a significant asymmetry in the pressure distribution. The flow in the bearing pad groove has an effect on the energy dissipation in the oil film. This paper provides a theoretical basis for the design and optimization of thrust bearings, and provides a reference for solving the problems of wear, oil mist, and other related problems of thrust bearings in engineering.

Hysteresis and ribbons in Taylor-Couette flow of a semidilute noncolloidal suspension

Hysteresis and ribbons in Taylor-Couette flow of a semidilute noncolloidal suspension

Research on the flywheel clearance flow and heat transfer characteristics in a horizontal canned motor reactor coolant pump

Flywheel is an important component for the canned motor reactor coolant pump (RCP). The flywheel of horizontal canned motor reactor coolant pump (HRCP) is different from RCP’s flywheel. Its clearance width or different thermal conductivity structure would affect the flow structure of flywheel clearance flow, and then affect temperature distribution of the clearance, which would further affects the life and stability of flywheel. Therefore, it is necessary to analyze the flow and heat transfer characteristics of flywheel clearance flow in detail. This paper applies a Detached Eddy Simulation (DES) method to study the flywheel clearance flow of HRCP. The results show that the DES method can capture the small-scale vortex structure in turbulence Taylor-Couette flow and then give a more refined flow field morphology. The authors further study the regularity of flywheel clearance flow with different thermal conductivity structure, and the results indicate that optimizing thermal conductivity structure can effectively reduce flywheel’s temperature. The authors also investigate the heat dissipation mechanism of different flywheel clearance structures.

Tricritical state and quasi-periodicity triggered by the non-linear elasticity in an Upper Convected Maxwell fluid confined between two co-oscillating cylinders about zero-mean

The effects of harmonically co-oscillating the inner and outer cylinders about zero mean rotation in a Taylor–Couette flow are examined numerically using Floquet theory, for the case where the fluid confined between the cylinders obeys the upper convected Maxwell model. Although stability diagrams and mode competition involved in the system were clearly elucidated recently by Hayani Choujaa et al. (2021) in weakly elastic fluids, attention is focused, in this paper, on the dynamic of the system at higher elasticity with emphasis on the nature of the primary bifurcation. In this framework, we are dealing with pure inertio-elastic parametric resonant instabilities where the elastic and inertial mechanisms are considered of the same order of magnitude. It turns out, on the one hand, that the fluid elasticity gives rise, at the onset of instability, to the appearance of a family of new harmonic modes having different axial wavelengths and breaking the spatio-temporal symmetry of the base flow: invariance in the axial direction generating the O(2) symmetry group and a half-period-reflection symmetry in the azimuthal direction generating a spatio-temporal Z2 symmetry group. On the other hand, new quasi-periodic flow emerging in the high frequency limit and other interesting bifurcation phenomena including bi and tricritical states are also among the features induced by the fluid elasticity. Lastly, and in comparison with the Newtonian configuration of this system, the fluid elasticity leads to a total suppression of the non-reversing flow besides emergence of instabilities with lower wavelengths. Such a comparison provides insights into the dynamics of elastic hoop stresses in altering the flow reversal in modulated Taylor–Couette flow.