Use of heated corrugations for propulsion

Abstract We show that a heating pattern applied to a surface with a corrugation pattern generates a propulsive effect. The same heating pattern applied to a smooth surface creates propulsion through nonlinear thermal streaming associated with pitchfork bifurcation. The combination of groove and heating patterns generates thermal drift, representing a forced response whose magnitude changes with the relative position of both patterns. Thermal drift is always present, while thermal streaming requires sufficiently intense heating. When both effects are active, a change in the relative position of these patterns produces rapid changes in the magnitude of propulsion, resulting in the formation of limit points. The strength of propulsion increases with a decrease in the Prandtl number and with an addition of uniform heating.


Introduction
The effectiveness of fluid transportation depends on the availability of propulsive devices.These devices take mechanical power from a source and convert it into a propulsive effect.As a wide range of phenomena is utilized for energy conversion, the resulting devices can be classified for the convenience of this discussion as either concentrated or distributed systems.Examples of concentrated propulsion include fans, pumps, jets, propellers, etc. whose analysis and designs attract plenty of attention.Concentrated systems are found in biology, e.g.fish propulsion using flapping foils and bird propulsion using a different form of flapping foils, and have become a subject of intense research in recent years (Wu 2011;Floryan, Burren & Smits 2018;Mannam & Krishnankutty 2019;Smits 2019).Distributed propulsion relies on the interaction between a solid surface and the adjacent fluid and is frequently used by small-scale biological organisms, e.g.snail locomotion (Chan et al. 2005;Lee et al. 2008), cilia and flagella (Taylor 1951;Blake & Sleigh 1974;Katz 1974;Brennen & Winet 1977;Lauga 2016), etc.Its use can also be found in large systems, e.g.bowel movements, although this is less common.Analysis and development of distributed systems for technical applications have been focused on the peristaltic effect, starting with the work of Shapiro and Fung and their co-workers (Fung & Yih 1968;Shapiro et al. 1969;Jaffrin & Shapiro 1971).Ali et al. (2020) provide a recent review.These analyses considered long slow waves, such as those encountered in biological systems.Current results demonstrate that fast, short peristaltic waves are more effective in providing a route to novel applications (Floryan et al. 2021;Floryan & Haq 2022;Haq & Floryan 2022).
Concentrated and distributed propulsions can be used simultaneously to achieve a performance not accessible to any of them individually.Active flow control can be viewed as a form of distributed propulsion, where improved performance is achieved by spending energy on reducing friction (Hoepffner & Fukagata 2009).There is a vast amount of literature on active control, but without referring to it as the development of distributed propulsion systems (Cattafesta & Sheplak 2011).The objectives of active control were narrowly defined as a search for net energy savings, i.e. energy savings larger than the energy cost required to achieve them, rather than a search for propulsion augmentation techniques.Such savings are not achievable in smooth conduits modified by wall transpiration (Bewley 2009); indications are that the situation can differ in non-smooth conduits (Floryan 2023).Results in the literature demonstrate that proper surface topography modifications reduce resistance (Walsh 1983;Fukagata et al. 2009;Mohammadi & Floryan 2013a,b, 2015;Moradi & Floryan 2013).It is also possible to produce chaotic stirring at an energy cost less than required to maintain the same flow rate in a smooth conduit (Gepner & Floryan 2020).
Forcing (actuation) patterns imposed on the contact surface between fluid and solid create distributed propulsion.These actuations reduce the wall shear acting on the fluid and increase its velocity relative to the solid boundary.Numerous concepts have been investigated, e.g.plasma-based actuators (Inasawa et al. 2013) and piezo-driven actuators (Fukunishi & Ebina 2001).See also the review provided by Cattafesta & Sheplak (2011).More recent contributions were focused on the thermal drift created by the interaction of groove and heating patterns (Abtahi & Floryan 2017;Inasawa et al. 2021) through the pattern interaction effect (Floryan & Inasawa 2021) on pumping driven by thermal drift (Floryan, Haq & Panday 2022), on propulsion driven by nonlinear thermal streaming (Floryan, Panday & Aman 2022) and on nonlinear streaming created by distributed wall transpiration (Jiao & Floryan 2021a,b).
Body forces distributed within the fluid volume can create distributed propulsion.Heating patterns produce a propulsive force and reduce viscous friction simultaneously (Hossain et al. 2012;Floryan & Floryan 2015;Hossain & Floryan 2016).Utilizing available waste heat reduces the energy cost required by such heating.While groove patterns generally increase friction (unless one uses unique shapes), a judicious combination of heating and groove patterns can reduce pressure losses (Hossain & Floryan 2016;Floryan, Shadman & Hossain 2018).Sound represents another form of volume actuation, but little is known about its effectiveness (Kato et al. 1997).
This paper describes the use of heated grooves to generate distributed propulsion.The model problem used in the analysis is described in § 2. Section 3 discusses the thermal drift effect and provides its essential characterization.Section 4 reviews the properties of nonlinear thermal streaming.Propulsion, when both effects are simultaneously active, is discussed in § 5. Section 6 provides a summary of the main conclusions.

Problem formulation
Consider two horizontal plates with a gap between them filled with a fluid, as shown in figure 1.The upper plate is smooth, while the lower one has sinusoidal corrugations, resulting in the gap geometry of the form where the subscripts L and U refer to the lower and upper plates, respectively, A * is the (peak to bottom) amplitude of the groove, α * is its wavenumber and stars identify dimensional quantities.The gap extends to ±∞ in the x*-direction, its mean opening is 2h * and it is periodic with wavelength λ * = 2π/α * .The gravitational acceleration g * is acting in the negative y*-direction.The fluid is assumed to have a thermal conductivity k * , specific heat c * , thermal diffusivity κ * = k * /ρ * c * , kinematic viscosity ν * , dynamic viscosity μ * , thermal expansion coefficient Γ * and variations of its density ρ * follow the Boussinesq approximation.Conditions required for the use of this approximation are discussed by Tritton (1977).Results of experiments for thermal conditions similar to those used in this analysis (Inasawa et al. 2019(Inasawa et al. , 2021;;Floryan & Inasawa 2021) demonstrate that the Boussinesq approximation captures well the fluid response.These experiments included qualitative (flow topology) and quantitative (velocity measurements) comparisons with the theoretical model.The upper plate is free to move within its plane, and the lower plate is stationary.The upper plate does not move in the case of an isothermal system, as no mechanism is available to generate a propulsive force.The introduction of patterned heating produces convection in the slot, and it is interesting to determine if this convection can force movement of the upper plate.We assume that the heating applied to the lower plate results in a simple sinusoidal temperature distribution while the upper plate is kept isothermal, i.e. (2.2a,b) where the subscripts 'mean' and 'per' refer to the mean and periodic parts, respectively, T * per is the peak-to-peak amplitude of the periodic component and Ω * denotes the phase shifts between the heating and groove patterns.We use T * U for reference with all material 980 A41-3 properties and define the relative temperature (2.3a,b) where Finally, we use half of the mean gap height h * as the length scale and κ * ν * /(g * Γ * h * 3 ) as the temperature scale to find the dimensionless expressions for the geometry and temperatures of the form (2.5a,b) where  (2.6c,d) where (u, v) are the velocity components in the (x, y) directions, respectively, scaled with U v * = ν * /h * as the velocity scale, p stands for the pressure scaled with ρ * U * 2 ν as the pressure scale and Pr = ν * /κ * denotes the Prandtl number.These equations are subject to the following boundary conditions: (2.7d-g) One must impose constraints in the form of to prevent any external pressure gradient from affecting the fluid movement.If convection generates a mean shear stress at the upper plate, then this plate accelerates until it reaches velocity U top , which eliminates this stress -all forces acting on the upper plate must balance out, so, since the only permitted force is due to shear, the mean shear must be zero.The periodic components of shear do not matter as they cancel out over one period.
The velocity of the upper plate U top is not known a priori -it is determined from the known solution during the postprocessing.Identification of the plate driving forces requires evaluation of the surface forces acting on the plates.The stress vector at the lower plate σ L is expressed as where , is the normal unit vector pointing outwards.The x-component of this vector has the form where (σ xv,L , σ xp,L ) denote the viscous and pressure contributions, respectively.The x-component of the total force F x,L acting on the fluid at the lower plate per unit length has the form where F xp,L and F xv,L stand for the pressure and viscous contributions, which must be equal and opposite.Similar expressions for the upper plate have the form and the total force must be null.System (2.6)-(2.8)must be solved numerically with sufficient accuracy to provide a data quality suitable for analysing the effects of geometric irregularities and heating and their interactions.Modelling the irregularities requires using techniques able to accurately handle many geometries with an acceptable cost.The spectrally accurate immersed boundary method accomplishes this, as it provides both the necessary geometric flexibility and accuracy and the ability to handle many geometries with a minimum labour cost.
We express the field equations in terms of the stream function defined in the usual manner, i.e. u = ∂ψ/∂y, v = −∂ψ/∂x, and normalized with condition ψ(y L ) = 0. Further, we capture the x-dependences of this function, as well as temperature θ , by expressing them as Fourier expansions based on the wavenumber α.Then we substitute these expansions into the field equations and separate the Fourier modes formulating nonlinear, ordinary differential equations for the modal functions.We express the modal functions using Chebyshev expansions and employ the Galerkin projection method to construct algebraic equations for the expansion coefficients.The irregular shape of the lower plate is accommodated using the spectrally accurate immersed boundary conditions concept (Szumbarski & Floryan 1999;Husain & Floryan 2008, 2010;Husain et al. 2009).The resulting constraints expressing the flow and boundary conditions along the lower corrugated plate and boundary conditions along the smooth upper plate are accommodated using the tau method.The nonlinear algebraic equations resulting from the discretization process are solved using a fixed-point iterative technique.The discretization process provided spectral accuracy with the absolute error controlled by the number of Chebyshev polynomials and Fourier modes used in the computations.We determined that using 15 Fourier modes and 50 Chebyshev polynomials provided a minimum of five digit accuracy for all quantities of interest.The reader can find algorithm details in Panday & Floryan (2021).
The heating activates different physical mechanisms.We begin the discussion by presenting these mechanisms in detail and then discuss the propulsion generated by these effects.

Thermal drift effect
We begin the discussion with a description of the thermal drift effect.If the upper plate does not experience any external resistance, viscous shear associated with convective motion may accelerate this plate until it reaches a velocity U top -the mean shear becomes zero under this condition.
Typical velocity and temperature fields are illustrated in figure 2. When the hot spots are located at the groove crests (Ω = 0), the flow topology is made of counter-rotating pairs of rolls with the fluid moving upwards above the hot spots (figure 2a).When the hot spots are located at the groove bottoms (Ω = π), the flow field has a similar structure (figure 2c).All other positions of the hot spots result in the formation of a net horizontal flow directed either to the right (0 < Ω < π) or to the left (π < Ω < 2π), as illustrated in figure 2(b,d), respectively.The flow direction is dictated by the relative position of the heating and groove patterns and represents a special form of thermal drift (Abtahi & Floryan 2017, 2018;Inasawa et al. 2021) driven by the pattern interaction effect (Floryan & Inasawa 2021).
The strength of the drift pulling the upper plate is quantified in terms of the velocity U top .Typical variations of U top as a function of the phase shift Ω illustrated in figure 3(a) demonstrate that positioning the hot spots either at Ω = π/2 or at Ω = 3π/2 produces the fastest plate movement.This speed increases proportionally to the groove amplitude, as illustrated in figure 3(b).
Figure 4(a) illustrates the distribution of the x-component of the pressure force σ xp,L acting on the fluid at the lower plate that drives the fluid in the horizontal direction.The form of the pressure field is dictated by the convection pattern, whose structure, in turn, is dictated by the heating pattern.As this pattern moves relative to the grooves, the projection of the pressure field onto the surface topography changes, producing σ xp,L of different magnitudes and directions -this projection is responsible for horizontal fluid movement (Abtahi & Floryan 2017).The convective motion generates an x-component of the shear stress σ xv,L at the grooved plate with large periodic x-variations (see figure 4b).The mean value of the shear stress F xv,L /λ must be equal to the mean value of the pressure force F xp,L /λ as the upper plate does not support any forces.
Uniform heating applied to the lower plate increases the plate velocity as the intensity of convection increases, while uniform cooling produces the opposite effect, as illustrated in figure 5(a).Reduction of the Prandtl number increases the strength of convection as mixing associated with convective motion is less effective in reducing horizontal temperature gradients (Hossain & Floryan 2013, 2022), leading to an increase of U top , while its reduction has the opposite effect, as illustrated in figure 5(b).
The stability of thermal drifts was assessed using the technique outlined in Appendix A. No unstable eigenvalues have been found, which suggests that thermal drifts are stable in the range of parameters of interest in this analysis.

Nonlinear thermal streaming effect
If the upper plate does not experience external resistance and the lower plate is smooth, sufficiently intense periodic heating may generate plate movement through nonlinear thermal streaming (Floryan 2023).The fluid motion has the form of counter-rotating The plate movement results from a system bifurcation, as illustrated in figure 7. The heating intensity, as measured by Ra per,cr , must reach a minimum level that strongly depends on the heating wavenumber -the minimum corresponds to Ra per,cr ≈ 900, it occurs at α ≈ 2.1, rapidly increasing when the heating wavelength is too short or too long (see figure 7).The value of Ra per,cr rapidly increases when the Prandtl number increases as convection reduces the horizontal temperature gradients required to initiate this effect (see figure 7b).An addition of uniform heating minimizes the intensity of periodic heating required to produce bifurcation.In contrast, constant cooling has the opposite effect (see figure 7c).

Combined effect
When the thermal drift and streaming effects are active, the combined effect occurs for a sufficiently high heating intensity.One can estimate the minimum heating required to activate the streaming for the grooved plate by placing hot spots either at the groove peaks or at the groove troughs, as this eliminates thermal drift.Results in figure 8 demonstrate relatively minor changes in critical conditions to those found in the case of a smooth plate.
The same results also illustrate different responses when the hot spots are placed at the groove troughs and the groove peaks.The combined system response depends on the relative position of the heating and groove patterns, as illustrated in figure 9(a).The value of U top varies as a function of phase difference Ω in a regular, nearly sinusoidal manner for Ra per = 500, which means that the flow response is driven by thermal drift and nonlinear streaming plays no role.Increase heating intensity to Ra per = 800, i.e. an increase of a factor of 1.6, results in a rise in the maximum possible U top by a factor of 4.7 while retaining its sinusoidal variations with Ω.When heating intensity is further increased to Ra per = 900, variations of U top remain regular, except when the phase difference is near Ω = π, which corresponds to hot spots near the groove troughs.Nonlinear streaming becomes activated, leading to the appearance of limit points in variations of U top as a function of Ω -for such conditions, a slight change of Ω can change the direction of plate movement.System response shows hysteresis -when Ω increases towards Ω = π and crosses the value of 3.23, U top jumps from U top = 1.53 to U top = −2.46.When Ω decreases back, U top does not jump back until the phase shift reaches a value of Ω = 3.06, when it jumps from U top = −1.67 to U top = 2.43.The hysteresis loop is relatively narrow, i.e. it occurs between Ω = 3.06 and Ω = 3.23, with U top being positive at the upper side of the loop and negative at the lower side.Increasing the heating intensity to Ra per = 1000 results in the formation of additional limit points near Ω = 0, which correspond to hot spots near the groove peaks with another hysteresis loop forming there and a widening of the hysteresis loop near Ω = π.Stability analysis (see Appendix A for details) demonstrates that the formation of limit points corresponds to the loss of stability of a particular branch of the bifurcation diagram (figure 9b).No stationary unstable states were found past the limit points, which suggests that only non-stationary states may exist under such conditions.
Evolution of U top as a function of Ra per near Ω = 0, which corresponds to hot spots near the groove peaks, is illustrated in figure 10(a).Using Ω = 0 results in a symmetric pitchfork bifurcation, i.e. the upper plate moves either to the right or to the left in the same manner.Placement of the hot spot very close to the groove peak and on its left side (Ω = 0.08) produces plate movement to the right, with U top increasing monotonically with Ra per , i.e. no bifurcation takes place.A separate solution branch corresponding to the plate moving to the left appears when Ra per > 850.This branch does not connect to the branch describing plate movement to the right -it ends at a limit point at Ra per = 956 with U top = −1.49.Moving the hot spot further away to the left (Ω = 0.52) produces a similar response but with U top varying with Ra per in a more regular manner and with a new solution branch forming at a larger Ra per , i.e. the limit point of the new branch is at Ra per = 1056 with U top = −2.71.Placing a hot spot on the right side of the groove peak (Ω = −0.08)forces plate movement to the right with a new solution branch corresponding to negative U top appearing at a limit point Ra per = 956 and U top = 1.49.Moving the hot spot further away (Ω = −0.52)moves the limit point to Ra per = 956 and U top = −1.49.
Variations of U top as a function of Ra per near Ω = π, which correspond to hot spots located near the groove troughs, are qualitatively similar to those observed for Ω near 0 (see figure 10b).Using Ω = π leads to a symmetric pitchfork bifurcation.Moving the hot spot to the left of the trough (Ω = 0.08) produces negative U top , which is a monotonic function of Ra per .A second branch, corresponding to positive U top , appears at Ra per = 896 with U top = 1.39.Moving the hot spots further away to the left (Ω = 0.52) produces more regular variations of U top as a function of Ra per , with the second branch originating at Ra per = 1000 with U top = 2.7.Moving the hot spot to the right produces a near mirror image of flow response with limit points being at Ra per = 896 and U top = −1.39 for Ω = −0.08 and Ra per = 1000 and U top = −2.7 for Ω = −0.52.
Comparisons of the data displayed in figures 8-10 show regular variations of U top as a function of Ra per for Ra per <∼ 800, where thermal drift dominates, and regular variations for Ra per >∼ 1100 where, depending on Ω, either nonlinear streaming or thermal drift dominates, producing a nearly identical system response.The system response exhibits somewhat unpredictable but significant changes of U top , including changes of direction for Ra per varying between these limits.
The effectiveness of thermal drift increases proportionally to the groove amplitude, as illustrated in figure 3(b).When thermal streaming is active, variations of U top with A become complex and change with Ω and Ra per as illustrated in figure 11. Figure 11(a) shows the system response for Ra per = 900, which is just above the critical heating of Uniform heating strengthens thermal drift (see figure 5a) and reduces the critical Rayleigh number required for the onset of thermal streaming (see figure 7c).Effects of adding uniform heating when both effects are active are illustrated in figure 12.In general, uniform cooling reduces U top and eliminates bifurcations that would occur with uniform cooling.Heating has the opposite effect, as it lowers the critical Rayleigh number and increases the achievable U top .
The type of fluid strongly affects the magnitude of U top .The propulsive effect relies on horizontal temperature gradients, which diminish with an increase of Pr as this leads to more intense convection resulting in a reduction of these gradients.The thermal streaming and drift are affected similarly, as illustrated in figures 5(b) and 7(b).The cumulative effect is shown in figure 13.Reduction of Pr reduces the critical Rayleigh number and increases the magnitude of U top .A slight increase of Prandtl number from Pr = 0.71 to Pr = 1 eliminates bifurcations and leads to a significant reduction of U top .
Stability analysis was carried out using the formulation described in Appendix A. No unstable eigenvalues have been found, which suggests that the flows generated by the combined thermal streaming and thermal drift reported here are stable.

Summary
Analysis of the propulsive effect generated by heated grooves has been carried out.The model problem consists of parallel horizontal plates with the upper plate free to move and the lower plate equipped with grooves and subject to periodic heating.Perfectly tuned sinusoidal groove and heating patterns were considered, i.e. the same wavenumber described them.Heating forced movement of the upper plate whose velocity was used to measure propulsion effectiveness.Two effects contribute to propulsion.The first effect is produced by an interaction of groove and heating patterns creating a thermal drift representing a forced system response.The strength of this effect depends on the relative position of both patterns.The upper plate movement reaches the maximum when hot spots are located approximately half-way between the groove peaks and troughs.It changes direction depending if the hot spot is on the right or left side of the groove peak.Thermal drift disappears when the hot spots are located at the groove peaks or the groove troughs.
The second effect relies on nonlinear thermal streaming and requires sufficiently intense heating for its activation.This effect represents system bifurcation.Both the thermal drift and the nonlinear streaming are active simultaneously and generate plate movement regardless of the relative position of both patterns.Similar to pure thermal drift, the fastest plate movement occurs for hot spots half-way between the groove peaks and troughs.A somewhat slower movement results when hot spots are located near the groove peaks and troughs.Under such conditions, the direction and velocity of plate movement may change drastically in response to a slight change in the relative position of both patterns.The formation of limit points and hysteresis loops was observed in response to changes in the relative position of both patterns.
The magnitude of propulsion increases with a decrease in the Prandtl number and expands with the addition of a uniform heating component.Despite extensive testing, investigating the stability of flows generated by the heating failed to identify any unstable eigenvalue.
to the disturbance equations of the form ∂ω where σ = σ r + iσ i is the complex frequency with σ i describing the rate of growth of disturbances and σ r describing their frequency and c.c. stands for the complex conjugate.
Here, G D (x, y), Ω D (x, y) and κ D (x, y) are the x-periodic amplitude functions accounting for the spatial flow modulations.The spatial distribution of disturbances is characterized by the spanwise (μ) and streamwise (δ) wavenumbers, with their amplitude functions being periodic in x.The former completely characterizes the overall system periodicity in the spanwise direction, while its periodicity in the longitudinal direction depends on the ratio of α and δ.The system can be aperiodic in this direction for an irrational ratio of α and δ (incommensurable system) and could be periodic but with wavelengths varying by several orders of magnitude (commensurable systems).Wavenumber α is a control parameter as its value characterizes heating and topography patterns and can be selected arbitrarily.The stability analysis requires consideration of all possible δ values for a given α, which gives rise to commensurability effects.It also requires consideration of all possible μ values, but this does not lead to commensurability effects.Here, we follow the formulation proposed by Floryan (1997Floryan ( , 2002)), which properly handles commensurability; see also Panday & Floryan (2023).
Substitutions of (A5) into (A3) lead to the partial differential equations for the amplitude functions.Since these functions are x-periodic, they can be expressed in terms of Fourier expansions of the form

Figure 1 .
Figure 1.Schematic diagram of the flow system.

Figure 3 .Figure 4 .
Figure 3. Variation of the upper plate velocity U top (a) as a function of the phase shift Ω and (b) as a function of the groove amplitude A for α = 2, Pr = 0.71, Ra per = 500, Ra uni = 0. Thin dotted lines in (a) show Ω giving rise to maximum U top and in (b) identify the character of the asymptotic variations for A → 0.

Figure 5 .Figure 6 .
Figure 5. Variations of the upper plate velocity U top as a function of the phase shift Ω for (a) Pr = 0.71, Ra uni = −100, −50, 0, 50, 100, and for (b) Ra uni = 0, Pr = 0.5, 0.71, 7.In both cases, α = 2, A = 0.05 and Ra p = 500.Thick dashed lines in (a) show variations of the maximum of |U top | as a function of Ω. Term 'multiplied by 5' means that the quantity of interest was multiplied by a factor of 5 before plotting.

Figure 9 .
Figure 9. (a) Variations of the upper plate velocity U top as a function of the phase shift Ω for α = 2, A = 0.025, Pr = 0.71, Ra uni = 0 and Ra per = 500, 800, 900, 1000.Circles identify limit points.The dotted lines provide the plate velocity achieved with smooth plates.(b) Variations of the disturbance amplification rate σ i (see Appendix A for definition) as a function of the phase shift Ω for flow conditions corresponding to red lines in (a).These disturbances are characterized by frequency σ r = 0 (see Appendix A for definition).

Figure 11 .
Figure11.Variations of the upper plate velocity U top as a function of the phase shift Ω for (a) Ra per = 900, A = 0, 0.001, 0.005, 0.02, and (b) for Ra per = 1000, A = 0, 0.005, 0.0093, 0.0094, 0.025.All results are for = 2, Ra uni = 0 and Pr = 0.71.Circles identify limit points.Black dashed lines provide the plate velocity achieved with smooth plates.Positions of limit points were determined using stability theory.Ra per = 880 for a smooth plate.No bifurcation occurs at A = 0.0004, but one occurs at A = 0.0006, producing wide hysteresis loops.Further increase of the amplitude to A = 0.005 dramatically narrows down the hysteresis loop.Further increase to A = 0.02 appears to change the magnitude of U top but does not affect the size of the hysteresis loop.Figure11(b) illustrates what happens when heating corresponds to Ra per = 1000.No bifurcation occurs for A = 0.005 and A = 0.0093, but one is formed when the amplitude reaches A = 0.0094.The results show that the amplitude must increase by a factor of nearly 10 to produce bifurcation when heating increases from Ra per = 900 to Ra per = 1000.The hysteresis loop is very narrow at A = 0.0094 but significantly widens when the amplitude increases to A = 0.025, which demonstrates a trend opposite to that observed for Ra per = 900.Uniform heating strengthens thermal drift (see figure5a) and reduces the critical Rayleigh number required for the onset of thermal streaming (see figure7c).Effects of adding uniform heating when both effects are active are illustrated in figure12.In general, uniform cooling reduces U top and eliminates bifurcations that would occur with uniform cooling.Heating has the opposite effect, as it lowers the critical Rayleigh number and increases the achievable U top .The type of fluid strongly affects the magnitude of U top .The propulsive effect relies on horizontal temperature gradients, which diminish with an increase of Pr as this leads to more intense convection resulting in a reduction of these gradients.The thermal streaming and drift are affected similarly, as illustrated in figures 5(b) and 7(b).The cumulative effect is shown in figure13.Reduction of Pr reduces the critical Rayleigh number and increases the magnitude of U top .A slight increase of Prandtl number from Pr = 0.71 to Pr = 1 eliminates bifurcations and leads to a significant reduction of U top .Stability analysis was carried out using the formulation described in Appendix A. No unstable eigenvalues have been found, which suggests that the flows generated by the combined thermal streaming and thermal drift reported here are stable.

Figure 12 .Figure 13 .
Figure 12.Variations of the upper plate velocity U top as a function of the phase shift Ω for α = 2, A = 0.025, Ra per = 1000, Pr = 0.71 and Ra uni = −100, −50, 0, 50, 100.Circles identify limit points.Dashed lines provide the plate velocity achieved with smooth plates.Positions of limit points were determined using stability theory.
-d)where j is the unit vector along the vertical y-direction.The relevant boundary conditions have the formy = 1: u D,per = ∂u D ∂y mean = v D = w D = θ D = 0, (A4a) y = y L : u D = v D = w D = θ D = 0. (A4b)System (A3)-(A4) defines a linear stability problem for the spatially modulated flow with modulations characterized by the wavenumber α.The disturbances are represented as[V D , ω D , θ D ](x, y, z, t) = [G D , Ω D , κ D ](x,y) e i(δx+μz−σ t) + c.c., (A5) y), g m η ( y), g m φ ( y)] e imαx + c.c., (A6b) κ D (x, y) = +N D m=−N D [g m θ ( y)] e imαx + c.c., stand for the modal functions and N D is the number of Fourier modes retained in the Fourier expansions used to represent disturbance quantities.The basic state flow quantities are also expressed as Fourier expansions, i.e. [u B , v B , θ B ](x, y) = denotes the number of Fourier modes retained in the representation of the basic state.Substitution of (A6)-(A7) into the partial differential equations for the amplitude functions, elimination of g m u ( y) and g m w ( y) and separation of Fourier modes leads to the ordinary differential equations for the modal functions for g is the uniform Rayleigh number measuring the intensity of the uniform heating component and Ra per = g * Γ * h * 3 T * per /(κ * ν * ) is the periodic Rayleigh number measuring the intensity of the periodic heating component.Determination of whether the heating can create a movement of the upper plate requires a solution of the field equations of the form