Title: Linear stability of the flow of a second order fluid past a wedge

Citation
Cracco M, Davies C, PhillipsTN (2020). Linear stability of the flow of a second order fluid past a wedge. Cardiff University. http://doi.org/10.17035/d.2020.0112209117


Access Rights: Data is provided under a Creative Commons Attribution (CC BY 4.0) licence
Access Method: Click to email a request for this data to opendata@cardiff.ac.uk

Cardiff University Dataset Creators

Dataset Details
Publisher: Cardiff University
Date (year) of data becoming publicly available: 2020
Coverage start date: 02/04/2019
Coverage end date: 01/07/2020
Data format: .m
Software Required: Matlab
Estimated total storage size of dataset: Less than 100 megabytes
Number of Files In Dataset: 12
DOI: 10.17035/d.2020.0112209117

Description

The linear stability analysis of Rivlin-Ericksen fluids of second order is investigated for boundary layer flows, where a semi-infinite wedge is placed symmetrically with respect to the flow direction. Second order fluids belong to a larger family of fluids called Order fluids, which is one of the first classes proposed to model departures from Newtonian behaviour. Second order fluids can model non-zero normal stress differences, which is an essential feature of viscoelastic fluids. The linear stability properties are studied for both signs of the elasticity number K, which characterises the nonNewtonian response of the fluid. Stabilisation is observed for the temporal and spatial evolution of two-dimensional disturbances when K > 0, in terms of increase of critical Reynolds numbers and reduction of growth rates, whereas the flow is less stable when K < 0. By extending the analysis to three-dimensional disturbances, we show that a positive elasticity number K destabilises streamwise independent waves, while the opposite happens for K < 0. We show that, as for Newtonian fluids, the nonmodal amplification of streamwise independent disturbances is the most dangerous mechanism for transient energy growth which is enhanced when K > 0 and diminished when K < 0.

The following files generate Figures 2-15 and Tables I-II in the paper:


‘meanFlow’ generates the mean flow profiles for a given value of γ and a chosen range of K, then saves the
results to .csv files. The .csv files are then loaded and used by the following files.


'velProfiles' generates figures 2(a)-(f): Velocity profile and relative variation with respect to the Newtonian
profile for increasing and decreasing values of the parameter K; (a), (b) γ = 0 (flat plate); (c), (d) γ = 0.5
(flow past a wedge); (e), (f) γ = 1 (stagnation flow).


‘eigenvalues’ generates figures 3(a)-(b): Comparison between Newtonian and non-Newtonian eigenspectrum
for the temporal problem for flow over a flat plate (γ = 0) with α* = 0.179, Re = 580 and (a) K = 0.03, (b) K
= −0.03. The least damped eigenvalues are those in the grey circle.


‘growthRatesSpatial’ generates figure 4(b): Spatial growth rates for flow over a flat plate (γ = 0) and Re =
580. Newtonian and non-Newtonian cases with K = 0.01, −0.05.


‘growthRatesTemporal’ generates figure 4(a) and 5(a)-(d). Figure 4(a): Temporal growth rates for flow over a
flat plate (γ = 0) and Re = 580, Newtonian and non-Newtonian cases with K = ±0.01. Figure 5: Temporal
growth rates for a flow past a wedge, Newtonian and non-Newtonian cases; (a) γ = 0.5, Re = 10000, K =
3e−04; (b) γ = 1 (stagnation point), Re = 27000, K =1e−04 ; (c) γ = 1.2, Re = 27000, K = 1e−04 ; (d) γ =
−0.14 (inflection point), Re = 300, K = 0.05.


‘neutralCurves’ generates figures 6(a)-(e): Temporal neutral curves in the Newtonian and non-Newtonian
cases. (a) γ = 0 (flat plate), K̃ = ±1e+03 ; (b) γ = 0.5, K̃ = ±1e+04 , x0 = 1; (c) γ = 1 (stagnation point), K̃ =
±2.5e+04 , x0 = 1; (d) γ = 1.2, K̃ = ±5e+04 , x0 = 1; (e) γ = −0.14 (inflection point), K̃ = ±100; x0 = 1.
‘neutralCurves’ can also be used, with appropriate modifications, to calculate the critical Reynolds numbers
of Table I.


‘stab3DNewt’ generates figures 7(a)-(b): Contour plot for the temporal growth rate, ωi, in the Newtonian
case (K = 0) for the flat plate (γ = 0). The red asterisk (∗) represents max ) represents max α,β ωi. The black line represents the neutral curve. (a) Re0 = 500; (b) Re0 = 1000.


‘stab3D’ generates figures 8-10. Figures 8(a),(b): Contour plots for ωi in the non-Newtonian cases for
the flat plate (γ = 0) and Re0 = 500. The red asterisk (∗) represents max ) represents maxα,β ωi. The black line represents the neutral curve. (a) K = −0.001; (b) K = 0.001. Figures 8(c),(d): Comparison of Newtonian (-) and non-Newtonian (- -) temporal growth rates for (c) α = 0.02; (d) β = 0.2.


Figures 9(a),(b): Contour plots for ωi in the non-Newtonian cases for the flat plate (γ = 0) and Re0 = 1000.
The red asterisk (∗) represents max ) represents maxα,β ωi . The black line represents the neutral curve. (a) K = −0.0001; (b) K = 0.0001. Figures 9(c),(d): Comparison of Newtonian (-) and non-Newtonian (- -) temporal growth rates for (c) α = 0.02; (d) β = 0.2.


Figures 10(a),(b): Contour plots for ωi in the non-Newtonian cases for the flow past a corner (γ = −0.14) and
Re0 = 150. The red asterisk (∗) represents max ) represents maxα,β ωi. The black line represents the neutral curve. (a) K = −0.003; (b) K = 0.003. Figures 10(c),(d): Comparison of Newtonian (-) and non-Newtonian (- -) temporal growth rates for (c) α = 0.02; (d) β = 0.2.


‘transientGrowthAlphaBeta’ generates figures 11-12. Figure 11: Contour plot of G max for γ = 0 (flat plate)
and Re0 = 1000. The black line indicates where an exponentially unstable mode exists. (a) K = 0; (b) K = 10
−4 ; (c) K = −1e−04. Figure 12: Contour plot of G max for γ = 0.5 and Re0 = 500. (a) K = 0; (b) K = 1e−04;
(c) K = −1e−04.


‘transientGrowthAlphaBeta’ can also be used, with appropriate modifications and change of scalings, to
generate Table II: Largest global optima for Reθ = 166 and Reθ = 385. The asterisk (∗) represents max ) indicates where an exponentially unstable mode exists and GΓ is calculated excluding the TS wave. The missing values indicate where an exponential unstable mode exists also as β→0.


‘transientGrowthAlphaWi’ generates figure 13(a)-(b): Ratio of non-Newtonian to Newtonian maximum
possible amplification for the flat plate γ = 0 and Re0 = 500, β = 0.6. (a) Gmax/Gmax,Newt ; (b) tmax/tmax,Newt .


‘optimalDistCompare’ generates figure 14(a)-(d): Comparison between Newtonian and non-Newtonian
optimal disturbances for γ = 0.5, Re0 = 500, α = 0.6, β = 0. (a) wall-normal, |v0|, and spanwise, |w0|, initial
velocities; (b) wall-normal, |vmax|, and spanwise, |wmax|, velocities at t = tmax ; (c) streamwise, |u0|, initial
velocities; (d) streamwise, |umax |, streamwise velocity at t = tmax .


‘optimalDist’ generates figure 15(a)-(d): Optimal disturbance for the stagnation point flow with γ = 1, Re0 =
500, α = 0.6, β = 0 and a non-Newtonian parameter K = 0.0001. (a),(b) disturbance at t = 0; (c),(d)
disturbance at t = tmax.


The figures can be obtained by specifying parameters (gamma – angle parameter, Re – Reynolds number
etc.) given in the captions of the figures.




Keywords

Linear stability

Related Projects

Last updated on 2020-31-07 at 10:16