Viscoelastic Thermovibrational Flow Driven by Sinusoidal and Pulse (Square) Waves

: The present study aims to probe the role of an inﬂuential factor heretofore scarcely considered in earlier studies in the ﬁeld of thermovibrational convection, that is, the speciﬁc time-varying shape of the forcing used to produce ﬂuid motion under the effect of an imposed temperature gradient. Towards this end, two different temporal proﬁles of acceleration are considered: a classical (sinusoidal) and a pulse (square) wave. Their effects are analyzed in conjunction with the ability of a speciﬁc category of ﬂuids to accumulate and release elastic energy, i.e., that of Chilcott–Rallison ﬁnitely extensible nonlinear elastic (FENE-CR) liquids. Through solution of the related governing equations in time-dependent, three-dimensional, and nonlinear form for a representative set of parameters (generalized Prandtl number Pr g = 8, normalized frequency Ω = 25, solvent-to-total viscosity ratio ξ = 0.5, elasticity number ϑ = 0.1, and vibrational Rayleigh number Ra ω = 4000), it is shown that while the system responds to a sinusoidal acceleration in a way that is reminiscent of modulated Rayleigh–Bénard (RB) convection in a Newtonian ﬂuid (i.e., producing a superlattice), with a pulse wave acceleration, the ﬂow displays a peculiar breaking-roll mode of convection that is in between classical (un-modulated) RB in viscoelastic ﬂuids and purely thermovibrational ﬂows. Besides these differences, these cases share important properties, namely, a temporal subharmonic response and the tendency to produce spatially standing waves.


Introduction
Thermally driven flows represent a vast category of phenomena with extensive background applications in various technological realms and industrial areas. In nature (where they are omnipresent), such processes are typically driven by density differences, which, in turn, are due to gradients of temperature. An increase in temperature typically makes a fluid lighter; as a result, in a steady gravitational field, relatively warm parcels of fluid tend to rise, whereas cold ones sink. Thermally driven fluid motion, however, is not an exclusive prerogative of steady acceleration fields. Time-varying (with variable sign and/or direction) body forces can indeed support similar density-driven mechanisms and excite other fascinating forms of thermal convection [1,2].
Such apparently innocuous observation is indeed at the basis of a well-known dichotomy in the literature, that is, the distinction between the classical fluid motion of thermogravitational nature and the so-called "thermovibrational" convection. The latter is called in this way as the simplest approach to have a fluid undergoing an acceleration that changes continuously sign in time is to "shake" the system containing the fluid itself. Constant amplitude "vibrations" with a given frequency result in a constant amplitude acceleration that varies in time with the same frequency of the parent vibrations and has localized minima and maxima of equal strength (absolute value). In other words, constant amplitude vibrations can produce an acceleration signal that displays no mean value (i.e., the related time-averaged value over a period of the vibrations is zero [3]).
Fluid motion induced by constant-amplitude vibrations may be regarded as the natural counterpart of buoyancy flow produced by steady gravity (see, e.g., in [4]). In fact, over the years this realization has led to the development of two branches of investigation running in parallel in the literature. These have shown that both variants are extremely sensitive to the relative inclination of the (steady or oscillatory) acceleration and the dominant temperature gradient. If they are parallel (e.g., concurrent gravity and temperature gradient), fluid motion requires that a threshold is exceeded in terms of a proper characteristic number (the so-called Rayleigh number). If they are not parallel (a special case being the condition in which they are perfectly perpendicular), then fluid motion is enabled regardless of the "position" of the considered case in the space of parameters (see, e.g., in [5] and references therein).
In the present study, we will concentrate on the first paradigm, i.e., the angle between acceleration and temperature gradient being exactly zero.
In these circumstances, when the acceleration is constant in time, the spatio-temporal behavior of the flow, which emerges from an initially quiescent state as the control parameter is increased, strongly depends on the nature of the considered fluid. For Newtonian fluids (where the only stresses present in the fluid are of a purely frictional nature), convection typically emerges in the form of steady parallel rolls. Time dependence is typically excited through a secondary bifurcation of the Hopf type (this concept being connected to the existence of well-defined instability regions in the space of parameters, e.g., the Busse balloon for the infinite layer case [6][7][8][9][10][11]). For viscoelastic fluids (special liquids where the existence of stresses that have an elastic nature is also permitted), the primary flow is intrinsically oscillatory (i.e., the first bifurcation is a Hopf bifurcation) if the elasticity of the fluid is sufficiently high; in turn, this behavior can be immediately put in relation with the well-known notion of "overstability", i.e., the process by which the critical threshold is significantly lowered with respect to the equivalent Newtonian case [12][13][14][15][16][17][18][19][20].
Unlike thermogravitational fluid motion, however, the thermovibrational convection being considered in the present work is oscillatory by definition. Using a phenomenological perspective, the emerging flow may be seen as the combination of a "linear response" (fluid velocity and temperature oscillating with the same frequency of the forcing) and a nonlinear contribution. The properties of the latter are essentially determined by the equations governing the dynamics, which are intrinsically nonlinear per se (the nonlinearity residing essentially in the convective terms of the momentum and energy equations for standard Newtonian fluids, additional nonlinearity being brought in by the additional equations for the elastic stresses for the viscoelastic fluids). In general, the nonlinear contribution consists of various harmonics (oscillatory modes with frequencies that are integer multiples of the forcing one) and a "stationary" component, i.e., a value obtained by time-averaging the fluid velocity (or temperature) over the period of the forcing, which is not zero. Various studies conducted under the assumption of Newtonian fluid have shown that the relative importance of the linear (fluctuating) response and the time-averaged contribution depends essentially on the amplitude of the vibrations and their frequency; the aforementioned mean value typically becomes significant as the frequency grows if the amplitude of vibrations is small (this resulting in the so-called Gershuni's regime of fluid motion [3,[21][22][23][24][25][26][27]).
Vice versa, the fluctuating component is dominant if relatively small frequencies are considered [28,29]. Investigations carried out for viscoelastic fluids in this condition have led to the conclusion that the concept of overstability is still applicable [30]. Moreover, it has been shown that in terms of spatio-temporal behavior, three-dimensional (3D) standing waves are the preferred outcome for these fluids (i.e., mechanisms of convection where localized plumes periodically invert their sense of motion, giving rise to extremely ordered lattice-like patterns [30]).
The equivalent 3D studies conducted for thermogravitational flows in viscoelastic liquids have revealed quite different dynamics where elongated rolls are still possible. It has been found that the parallel convective structures typical of steady Rayleigh-Bénard (RB) convection are replaced by a peculiar mechanism of oscillations where rolls continuously break and reform in time with a new orientation in space (which is essentially perpendicular to the original one [31,32]).
In such a context, the present study aims to probe the role of an influential factor that before now has been scarcely considered, that is, the specific (temporal) shape of the forcing used to produce thermovibrational convection. All the studies cited in this introduction were indeed obtained assuming for the forcing a canonical waveform, that is, a classical sinusoidal acceleration signal. The main objective of the present numerical investigation is the assessment of the modifications induced in the patterning behavior of a viscoelastic fluid when classical sinusoidal vibrations are replaced with a "square wave". This can be seen as a non-sinusoidal (yet periodic) forcing where the amplitude alternates (yet at a constant frequency) between fixed minimum and maximum values, with the same duration at minimum and maximum and instantaneous transitions between such extrema (also known as pulse waves). Such a change in the waveform is expected to have a remarkable impact on the dynamics of viscoelastic fluids given the ability of these liquids to accumulate and release energy according to intrinsic temporal scales [30].

Mathematical Model
Following the same rationale underpinning [30], a shallow infinite layer of aspect ratio AR = w/ = 15 is considered, where w represents the extension in the x and z directions (square basis), and accounts for the depth of the layer. To mimic an infinite domain, cyclic (periodic) boundary conditions are applied along the sides of the layer, while no slip walls play the role of top and bottom boundaries. Such solid walls are differentially heated in such a way that a temperature gradient is established within the liquid. The related difference of temperature reads ∆T = T h − T c , where the subscripts h and c denote the "hot" and "cold" sides, respectively. A scheme of the domain is depicted in Figure 1. Proper casting of the balance equations in compact form requires a model is introduced for the buoyancy force. In this regard (although in this case the acceleration acting in the system is not steady), the Oberbeck-Boussinesq approximation can still be considered valid. This must be associated with an acceleration that is an harmonic function of time, namely, H(t * ) (the symbol * indicates that the quantity is dimensional). This function depends on the "type" of the vibrations that are exciting the system. The classical way to model harmonic vibrations for thermovibrational flows (see, e.g., in [1,[3][4][5]29]) is to consider a sinusoidal displacement s(t * ) = b sin(ωt * )n, where b is the amplitude, ω is the angular frequency (ω = 2π f ),n is the direction of the vibrations. In such a case, taking the second derivative of s(t * ) yields

Cyclic Boundary Conditions
where γ = −bω 2n represents the amplitude of the acceleration itself.
It is straightforward to observe that, if in place of a sinusoidal acceleration, another type of forcing is needed, the sin function in Equation (1) must be replaced with the desired waveform. For the case of a pulse (square) wave, this reads [33] H(t * ) = γ tanh(10 sin(ωt * )) (2) For the sake of simplicity, hereafter the generic harmonic acceleration will be referred to as H(t * ) = γh(t * ).
As explained to a certain extent in the introduction, the considered fluid is assumed to be viscoelastic. In line with our preceding work [5,31], this property is modeled in the framework of the so-called finitely extensible nonlinear elastic model in Chilcott-Rallison form, generally known as FENE-CR [34,35]. This specific paradigm stems from the idea that a viscoelastic fluid can be seen as a mixture of a Newtonian solvent (having dynamic viscosity η s ) and a polymer solute (having dynamic viscosity η p ) with molecules roughly behaving as couples of beads inter-linked by a spring (dumbbell model). The resulting viscoelastic fluid has a total viscosity η 0 = η s + η p . It is common practice to define two mutually dependent dimensionless numbers related to the viscosity, i.e., the solvent-to-total viscosity ratio ξ = η s /η 0 and the viscosity ratio This specific model allows to couple the extra viscoelastic stress tensorτ * with the momentum equation. This new quantity is evaluated with a new transport equation.
To summarize, in the frame of this approach, the fluid behavior is governed by the continuity, momentum, energy, and viscoelastic stress tensor balance equations, which in dimensional form read as follows: where t * is the time, u * is the velocity, T * is the temperature, p * is the pressure, ρ is the density of the fluid, β T is the thermal expansion coefficient, α is the thermal diffusivity, and λ is the so-called relaxation time. The function f (tr(τ * )) is a quantity related to the possible deformation of the polymeric molecule, which in the FENE-CR is limited and modelled through the the so-called finite extensibility parameter (L 2 ); it reads The balance equations can then be nondimensionalized using , α/ , 2 /α, α/ 2 , ρα 2 / 2 , ∆T, and ρν s α/ 2 (where ν s = η s /ρ) as scaling factors for the length, velocity, time, frequency, pressure, temperature, and extra-stress tensor, respectively. This operation generates the non-dimensional groups that fully characterize the problem under consideration, namely, the Prandtl number for the Newtonian solvent Pr = ν s /α, a generalized version of this parameter for the viscoelastic fluid Pr g = Pr/ξ, the non-dimensional frequency Ω = 2 ω/α, the elasticity number ϑ = λα/ 2 , ξ and ζ (both previously defined), and the vibrational Rayleigh number: where ν 0 = η 0 /ρ is the total kinematic viscosity. The balance equations in their dimensionless form therefore can be cast in condensed form as

Numerical Method
The computational platform openFOAM ® is employed here to solve numerically the set of balance equations previously defined. In particular, the related segregated timemarching procedure directly relies on the PISO (Pressure Implicit with Splitting of Operator) algorithm to couple velocity and pressure fields, and the Rie and Chow scheme [36] to avoid checkerboarding of the latter one.
The balance equations have been discretized over a uniform Cartesian grid having 200 × 35 × 200 cells (see in [30] for the complete grid independence analysis). All the diffusive and convective terms have been evaluated though a second-order accurate central difference scheme with the exception of the convective term of Equation (12), where the Minmod scheme has been used to improve the algorithm stability. Furthermore, to mitigate the effect of possible numerical singularities, the numerical solver implements Equation (10) as [37]: As the reader will realize by comparing Equations (10) and (13), an extra diffusive term Prζ∇ 2 u has been added to both sides of the momentum equation. Although the two equations are equivalent from a mathematical standpoint, in the numerical framework, the two extra diffusive terms can be integrated using different approaches (namely, one is integrated implicitly while the other one explicitly). This technique generates quantitatively negligible numerical diffusion, which, however, improves the stability and the robustness of the solver, as witnessed by earlier studies (see in [5,[29][30][31]38] for further details about the code validation and its application to the case of both Newtonian and viscoelastic fluids).

Results
We focus on a layer of viscoelastic (FENE-CR) fluid with Pr g = 8, ξ = 0.5, and ϑ = 0.1 subjected to either sinusoidal or square vibrations having frequency Ω = 25. For both cases, the resulting vibrational Rayleigh number is Ra ω = 4000. Furthermore, the model parameter L 2 is set to 200 as in our previous studies (see in [5,31]). Let us recall that the FENE-CR can adequately represent highly elastic solutions known as "Boger fluids" (see, e.g., in [39]), able to retain an essentially constant viscosity over a wide range of shear rates. Relevant examples are represented by a class of water-based polymer dilute solutions at ambient or moderate temperatures, e.g., water between 25°C and 50°C with limited amount of a polymer such a PAM, PEG, PEO, PVP, Xanthan Gum, etc., for which the Prandtl number would be similar to that considered in the present work (Pr g ∼ = 8) and ξ < 1. The rheological parameters (Pr g and ξ) considered here are almost identical to those examined by Li and Khayat [39], who assumed a Boger fluid with Pr g ∼ = 7 (and ξ varying in the range between 0 and 0.79). Assuming λ ∼ = 10 −3 s (a typical realistic value for small polymer concentrations) and thermal diffusivity α ∼ = 1.5 × 10 −7 m 2 /s over this range of temperatures, these values would correspond to ϑ = 0.1 setting 0.04 mm as the distance between the plates.

Forcing and Related System Response
Given the complexity of the considered subject, before starting to deal with the structure and evolution of the fluid dynamic field, for the convenience of the reader, in this subsection the system response is initially assessed in terms of "localized" or "global" parameters (by which drawing initial conclusions on the behavior of the system is relatively straightforward).
In particular, in order to highlight differences and similarities between the two "shapes" of the acceleration signal, we concentrate on the time evolution of two representative quantities, namely, the axial fluid velocity (the velocity component perpendicular to the solid walls) recorded by a virtual probe located in the center of the layer, i.e., (7.5, 0.5, 7.5), and an integral quantity, that is the Nusselt number calculated on the hot (bottom) plate as: where n plate in the unit vector normal to the boundary, and A is the area of the heated plate. Along these lines, Figure 2 provides the evolution in time of the two aforementioned characteristic thermal and fluid dynamic quantities together with the forcing. It can be seen that, when the external force (in red) h(t) > 0, i.e., the acceleration has a positive direction (with respect to the unit vectorn, see Figure 1), the effect of the vibration is stabilizing, i.e., the vibrations act as a damper for the convective instability. By contrast, when the external force h(t) < 0, the vibrations act in a way that is similar to the steady gravitational acceleration in the classical problem of Rayleigh-Bénard convection. Therefore, in this time interval vibrations are destabilizing, i.e., they promote convective motion. These simple arguments also clarify the reason why the Nusselt number and the forcing are out of phase. The Nusselt number is largely increased in the stage where vibrations play a destabilizing role and tends to the unit value (corresponding to quiescent conditions) as their effect is reverted.
Given this premise, it is easy to realize that, due the peculiar shape of the forcing, the switch between stabilizing and destabilizing phases can happen in different ways: the smooth transition, typical of the sinusoidal wave is taken over by an abrupt change in the acceleration direction for the case with the pulse wave. Moreover, as a result of the different shape of the forcing signal, the pulse wave might be expected to exert a stronger influence on the fluid.
In order to prove the last statement, it is worth comparing Nu(t) for the two cases. As quantitatively substantiated by the data, the maximum and the time average value of the Nusselt number are lower for the case with sinusoidal acceleration (Nu max = 1.88, Nu av = 1.18) and bigger for the pulse one (Nu max = 2.17, Nu av = 1.30). A justification for this trend can be elaborated in its simplest form on the basis of the argument that while the sinusoidal acceleration attains its maximum intensity for a single instant in time, the square acceleration maintains a constant maximum intensity for half of the period of oscillations.
The axial velocity signal can also be used to get useful insights into these dynamics. Interestingly, in the stabilizing semi-period, both signals (hereafter v S (t) and v P (t) for the sinusoidal and pulse waves respectively) reach a zero value at certain time. However, a careful analysis of their shape reveals that while v P (t) is stabilized to a null value in a certain interval of the stabilizing period, v S (t) changes continuously in this interval, i.e., dv S /dt = 0 in every time sub-interval (while it can locally reach 0). In the light of these initial findings, the reader may expect that when the the fluid is excited by pulse waves, it will display a more involved behavior in the destabilizing phase if compared with the case of sinusoidal forcing for the same value of Ra ω (we will come back to this interesting concept later).
As a concluding task for this initial analysis, we also calculated the frequency of the Nu(t) and v(t) signals. In agreement with the earlier results by Boaro and Lappa [30] and Lyubimova and Kovalevskaya [33], obtained with a nonlinear numerical approach and linear stability analysis, respectively, the flow field exhibits a half-subharmonic response to both types of wave, i.e., the frequency of v(t) is half of the frequency of the acceleration, that is Ω/2 (interestingly as shown in Figure 3, the related frequency spectrum displays three dominant frequencies, namely Ω/2, Ω + Ω/2 = 3Ω/2 and 2Ω + Ω/2 = 5Ω/2, where the second and the third components are produced by the obvious nonlinear interference between the flow sub-harmonic response and the forcing). Yet in accordance with the work in [30], Nu(t) shows a simple harmonic response (its frequency is exactly Ω). Here, we select t 0 = 17.94 for the case of sinusoidal acceleration and t 0 = 11.14 for the pulse acceleration (see Figure 2). We also take advantage of this subsection to introduce another parameter useful for the understanding of these phenomena, that is, the ratio between the relaxation time of the polymer λ and the period of the external forcing T ω = 2π/ω [5,30]. It is easy to verify that this parameter depends on two non-dimensional numbers defined in the preceding section and reads Σ = ϑΩ/2π. As in our simulation both ϑ and Ω are fixed, Σ = 0.4.

Patterning Behaviour and 3D Evolution
In this section (entirely dedicated to the analysis of the 3D behaviour of the considered flows), first we examine the 3D isosurfaces of the axial velocity for the case of sinusoidal shaking for t ∈ I 1 , see Figure 4. At the time t 0 (Figure 4a), the flow displays an ordered convective structure consisting of groups of parallel rolls each oriented along a given preferential direction. The evolution of this pattern is quite simple. As time progresses, the intensity of the axial velocity becomes smaller (Figure 4b) until the rolls change their sense of rotation (Figure 4c). This pulsation of the velocity field is a common solution for viscoelastic Rayleigh-Bénard systems and it is generally referred to in the literature as standing wave (see in [1] and reference therein).
Although after the inversion of the sense of rotation the intensity of the axial velocity tends to increase, it never reaches again the same value that it had at t = t 0 (Figure 4d). Indeed, the acceleration is changing direction, and therefore the stabilizing effect of the vibrations starts to influence the amplitude of convection. As evident in Figures 4e-h, interestingly, the rolls keep pulsating. However, this pulsation is no longer a result of a variation in the direction of the buoyancy force, rather it is driven by the residual elastic energy stored in the polymer molecules. This concept will be discussed afterwards.
Here, we limit ourselves to pointing out that, from a phenomenological point of view, in the initial part of the interval I 2 , Figure 5a, even if the acceleration is in its stabilizing phase, the intensity of the rolls increases again: the rolls, initially disperse in the domain, start to grow in intensity until they merge together (Figure 5a-d). It is also worth noticing that after the last inversion of sense, visible in Figure 4a-d, the rolls have not changed again their sense of rotation. However, as now the fluid motion is sustained only by the elastic energy stored in the fluid in the preceding destabilizing phase and not by the buoyancy force, the intensity of the velocity field becomes smaller as the elastic energy is dissipated through viscosity effects (compare Figures 4d and 5d).
As the flow evolves, the acceleration changes it sign and starts to play again a destabilizing role, thereby promoting the rise of the velocity field (see Figure 5f,g). Now, the buoyancy force that drives the fluid motion, along with the elastic energy, leads to a situation similar to the initial one (t = t 0 ), both for intensity and shape of the rolls, but with a different sense of rotation as evident in Figure 5h.
The second part of the oscillation period (I 3 and I 4 ) evolves in similar way to the one just described, and for the sake of brevity, we do not describe it in detail.
Rather, we concentrate on the companion case, that is, the situation with the pulse acceleration (see Figure 6 for the behavior over the interval I 1 ). Remarkably, it can be seen that the rolls are no longer aligned along a particular direction. Instead, they form randomly organized structures. The elongated shape that characterizes the case with sinusoidal vibration is now replaced by a mix between short longitudinal and rounded rolls as depicted in Figure 6a. Recalling the concepts introduced in the earlier subsection, this scenario might be seen as a realization of the increased complexity of the flow expected when the acceleration follows a pulse wave time evolution.
Although the spatio-temporal behavior differs in the two cases under analysis, however the leading mechanism of evolution remains basically the same, namely, a standing wave. Indeed, the velocity field initially undergoes a shrinkage in intensity (Figure 6b,c) until eventually the fluid changes its sense of rotation (Figure 6d,e). It is worth observing that also in this case the overall intensity of the velocity field after the first inversion is smaller than that in the initial frame (compare Figure 6a-e).
Despite these similarities, however, a closer look at the first part of the pattern evolution (Figure 6a-c) reveals that a secondary mechanism is at play in this case, namely, a breaking roll mode of convection [31,32]. This specific behavior has been already observed as a typical mode of convection in classical Rayleigh-Bénard viscoelastic flows. Its dynamics basically consists of convective rolls that break and reform but in a direction that is perpendicular to the original one. To visualize this process, as an example, the reader may consider the first roll in the corner of the layer in Figure 6a. It is easy to observe that, as the pattern evolves, this rolls tends to merge with the closest roll of the same sign (Figure 6b,c).
Nonetheless, this peculiar mechanism is mediated by the predominant standing wave and by the fact that in the meantime the vibrations became stabilizing for the flow field (and therefore, the fluid motion is no longer driven by buoyancy but only by the elastic energy stored by the polymer molecules).
When the vibrations become stabilizing, the standing wave is still the dominant solution. In this case the fluid has stored enough energy to continue pulsating and to invert the sense of rotation two other times (Figure 6e-g).
As the interval I 2 is entered, the vibrations are more or less at half of their stabilizing phase. In this stage, the acceleration has a constant value (as opposed to a sinusoidal acceleration, which changes continuously in time). This constant stabilizing acceleration has the effect to dampen and almost annihilate the convective structures. This becomes more evident when the patterning behaviors for sinusoidal and pulse forcing are compared in I 2 (compare Figures 5 and 7). As the system enters the second part of I 2 , the acceleration suddenly changes its direction. As a result, convective instabilities are newly excited and the velocity field starts increasing its intensity. The spatial arrangement of the rolls is the same as at the beginning of I 1 but with an inversion of the sense of rotation (as depicted in Figure 7e-h).
The dynamics in I 3 and in I 4 are relatively similar to the one just described and for this reason we omit to report them. However, it is worth highlighting that the breaking rolls mechanism is still effective over I 3 .
As a concluding remark for this subsection, one last point should be clarified. As we have illustrated on the basis of a purely heuristic approach, a key factor driving the fascinating dynamics of the system under analysis is its ability to retain elastic energy and release it when the external time dependent forcing does not play a destabilizing role. A classical way to estimate the local elastic energy in a viscoelastic fluid is to calculate the trace of the extra-stress tensor tr(τ) (this quantity being proportional to the local elastic energy) [40].
In Figure 8, we report the evolution of tr(τ), hereafter we will simply call this quantity "elastic energy", along a fixed line parallel to the x direction (belonging to the mid-height plane, i.e., 0 < x < 15, y = 0.5, and z = 7.5) as a function of time (over the interval 2T Ω ). On this basis, it is easy to see that moving from a sinusoidal to a pulse wave acceleration increases the amount of elastic energy stored by the polymer molecules. Further, note that the elastic energy field exhibits a harmonic response to the periodic forcing, regardless of its shape (waveform).

Discussion and Conclusions
The results presented in this study should be regarded as an extension of previous work of the present authors on these subjects. On one hand, we have considered conditions close to those examined by Boaro and Lappa [30], i.e., thermovibrational convection in an infinite layer of viscoelastic fluid having Pr g = O(10), and with a similar frequency of the external forcing (set as sinusoidal wave in Boaro and Lappa [30]). On the other hand, a FENE-CR fluid with L 2 = 200 having Pr g = 8, ξ = 0.5, and ϑ = 0.1 has been assumed (as in Lappa and Boaro [31], where attention was paid to classical Rayleigh-Bénard convection, i.e., the case with acceleration constant in time and space).
In the present work, the parameter Σ has been fixed to 0.4 (a value smaller than 1), which implies that in a period of vibrations, the polymer molecules are potentially able to relax and go back to their initial configuration. This is the feature allowing the system to release elastic energy when it enters the stabilizing phase of the vibrations (where such energy is progressively dissipated by viscous effects, similarly to what happens in an harmonic oscillator, see Figure 8).
As illustrated in the results section, in the case of sinusoidal forcing, the periodic excitation and dampening of convective modes (mediated by the accumulation and release of elastic energy) can generate complex-ordered structures.
Following the classification made by Rogers and coworker to categorize complex ordered structures arising from modulated Rayleigh-Bénard convection in Newtonian fluids [41][42][43][44], here, the observed patterns might be classified as roll superlattices. The roll superlattice is a subharmonic pattern emerging in states that are slightly supercritical. Once the Rayleigh number is increased and exceeds a certain threshold (which, in turn, depends on the amplitude of vibrations), these planforms can undergo a transition to a square superlattice (SQS) state [41].
Given these premises, it is interesting to compare the results of the present study (ϑ = 0.1), with those reported in [30] (considering in particular the cases with ϑ = 0.06 and 0.24). In that work, "more" supercritical solutions were examined and the obtained superlattice-like structures were akin to the SQS observed in the modulated Rayleigh-Bénard problem. This indicates that the general conclusion that an increase in the critical parameter can turn a roll superlattice into an SQS also holds for thermovibrational convection in a viscoelastic fluid.
Moving forward to the case with pulse wave acceleration, it can be concluded that the ordered patterns seen for the sinusoidal forcing are taken over by a distribution of rolls that is more similar to the one observed in classical RB convection in viscoelastic fluids. In particular, the initial distribution of the rolls is reminiscent of the flow pattern that arises in a shallow viscoelastic liquid bridge (see [31] where the following case was considered: A = 0.17, Ra = 2000, ϑ = 0.1, where A = height/diameter and Ra represent the classical Rayleigh number). In [31], the flow was found to display a peculiar breaking-roll evolution. In the present case, this mechanism (sustained by buoyancy forces) vanishes as soon as the system enters its stabilizing phase (the differences between the two systems can obviously be found in the different cause-and-effect relationships enabled by steady and time-varying accelerations).
For the sake of completeness, we wish to report that the ability of modulation to cause a departure from the tendency of classical RB convection to produce steady parallel rolls has also been observed in Newtonian fluids where in place of mechanical forcing (vibrations), thermal forcing was implemented (in terms of time-dependent top and/or bottom temperature, see Hohenberg and Swift [45] and Meyer et al. [46]). We take those studies, where such thermal modulations were proved to cause a competition between the standard parallel rolls and hexagonal cells through a mechanism akin to that effective in non-OB (Oberbeck-Boussinesq) systems, as a cue for a possible extension of the present research.
Another exciting prospect for the future is to investigate the vibrationally modulated problem for additional (non-unit) values of the parameter Σ, for which the relaxation time of the polymer molecules is smaller or larger than the period of the forcing. This factor is also expected to contribute significantly to the number of modes that can be excited in the fluid and alter the mechanisms depicted above, thereby influencing the emerging patterns. Funding: This work has been supported by the UK Space Agency (STFC grants ST/S006354/1, ST/V005588/1 and ST/W002256/1) in the framework of the "particle vibration" (T-PAOLA) project.

Data Availability Statement:
Publicly available datasets were analyzed in this study. These data can be found in the pure repository of the University of Strathclyde.

Conflicts of Interest:
The authors declare no conflicts of interest.