Melbourne to Perth

Even longer than last time! The return trip took 20 days to go around TAS, up the east coast, across the top end and down the west coast to complete one full lap.

Keeping the car (and myself) running in the October heat of the top end, after already 10 days on the road, was a lot more grueling than I expected. But I made it!

The schedule I developed was robust enough to allow for a few detours & side-trips, and the cost estimate was accurate enough that overall spending was only 8% over-budget, with fuel being 3% under!

There are some incredibly beautiful parts of this country, which I’m already planning to go back to. And there are some parts where it’s a miracle anything lives there at all.

Pictured: Mt Wellington, Hobart; Hunter Valley, NSW; NT/WA Border

Scotty dogs

Developed with a coworker as a Christmas present for his wife – one for the garden, one for the desk

Frisbee flight

Flow vis of frisbee boundary layer separation (flow is from the top down) – from Potts, J. 2005

Frisbees are super interesting! They generate lift in a way similar to a wing, but because they are so short in a transverse direction, 3d effects dominate. The easiest way to keep stable flight with this sort of shape (aside from gyroscopic stabilisation i.e. throwing it) is to ensure a turbulent boundary layer – which is why frisbees have ridges! The ridges serve to trip the boundary layer and keep it under control, maintaining a more steady lift distribution. Frisbees without ridges are much more sensitive to angle of attack, and harder to throw properly.

I wanted to model this boundary layer separation – and I still have a lot of learning to do! Below is a Q contour animation of my latest model, using LES with the smallest elements on the order of the Taylor microscale. The frisbee was set up at a 5 degree angle of attack with a 15 m/s headwind, which is identical to the image above. The geometry does not have ridges, but it does have a rotating wall boundary condition (counter-clockwise when viewed from above). The weird upstream artifacts are due to a poor mesh from snappyHexMesh – the sharp transition between cell sizes interferes with the LES algorithm, which uses the cell size as a spatial filter – I’ll see if I can get a structured mesh set up.

The aerodynamic coefficients given by this model are C_l ~ 0.20 and C_d ~ 0.10 – the drag coefficient is close to that measured by Potts but the lift coefficient should be closer to 0.4. This likely has a lot to do with the poor mesh.

Vortex shedding

A quick study looking at the Strouhal number for a medium-length cylinder – above is an animation of Q isosurfaces, coloured by vertical velocity. I ran some 2d preliminary models, followed by a 3d model with a structured mesh that I then refined. The Reynolds number was 51,355, which according to data from Achenbach (1968) is in a transitional range for the coefficient of drag:

My results gave C_d = 0.73, which is higher than for a long cylinder the above chart indicates – likely due to the three-dimensionality of the flow around the ends.

A plot of C_l vs time clearly shows the oscillations of the Karman vortex sheet, at a frequency of about 37 Hz. This gives a Strouhal number of about 0.12, which is not far off experimental values of 0.18-0.50.

Also visible in the plot is that after about 0.4s, the oscillations become unstable and appear to ‘beat’. The RMS value does remain pretty constant at around 0.4, but I’m curious whether this is due to end effects or maybe just the integration schemes used!

Return to BFS flow

I’ve been working to improve my CFD skills, and have set up a workstation running OpenFOAM – what better case to practice on than one I’ve already done before?

I wasn’t entirely happy with my results from the last attempt, thinking most of the error was due to an improper mesh, so I ran the same case using a better mesh, and an assortment of turbulence models in 2d and 3d.

Streamlines in the 2d BFS flow

The reattachment length was found by evaluating the wall shear stress along the bottom wall of the downstream section – where this is zero is the reattachment point. I found the standard k-\epsilon model did not reach mesh convergence at any reasonable point, but both k-\omega and realisable k-\epsilon did, giving x_R/H ~ 4.8 and 5.6 respectively. Both of these are still below the real steady-state value of 7.0, but are much more promising than last time!

I used k-\omega for the 3d case as well – it was developed specifically for internal flows and is supposedly the best RANS model for BFS flow, but it looks like I need to work on my calibration as I got an x_R/H of 4.5! This may also be because I imposed a symmetry condition on the centre of the duct – the flow has been reported as two-dimensional along this plane but there are transient 3d effects to take into account.

Three-dimensionality of the near-wall flow

Perth to Melbourne

That was a long drive! 9 days along the south coast.

It was smooth sailing except for a few hiccups – thankfully I brought some tools and the MX5 workshop manual or it would have been a lot longer than 9 days…


Sheet Metal Guillotine

This was designed for a machine design unit.

Requirements were to cut 3 mm mild steel sheet, up to 2400 mm wide, and be powered by one hydraulic cylinder while still maintaining a straight & level cut.

Several iterations led to a design that could be CNC machined from one single 20 x 1800 x 6000 plate, with a few extra small parts.