MATLAB code that solves and compiles the steady flow around a rotating sphere with variable viscosity of the form discussed in Miller et al (2020). The package contains code that solves the heated, rotating sphere boundary layer, the new model of the impinging region and the subsequent eruption of the radial jet.
All numerical methods are the same as those in the Steady Sphere Flow code. Full details can be seen in Smith et al (2023) [in draft].
To run the code, download all the files making sure the directory structure remains intact. Then simply enter BaseFlow(Re,Pr,λ) in the command window, where Re is the square root of the Reynolds number, Pr is the Prandtl number (~0.72 is the value for air) and λ is a sensitivity constant that represents a change in fluid (if set to 0 then the temperture independent case is solved). The data will be saved in a directory called Flows that will be created automatically. To see the results run the file figures.m.
