COMPUTATIONAL FLUID DYNAMICS MODEL

This Computation Fluid Dynamic (CFD) model of a sugar beet clamp has been developed using the OpenFOAM open source software. The table below gives links to posts that describe the development of the model.

NB: owing to this work being under development for publication, the publication on this website of the details of the model development have stopped. The final model count is closer to 80, with many of them run multiple times. The expected manuscript submission date is end 2022. At the latest, the paper will be in the thesis this work is a part of, with a defense date set for 26 May 2023. That said…

The github repository for these models are public. You need an OpenFOAM installation to run the model, ParaView to inspect the results, and a fair bit of specialised knowledge to modify the model (properly) …

The physical model is fairly idealised, but also has a fairly good resolution and is built on some sound experimental data.

Some of the key features of the model include:

• separate momentum and energy equation solvers. The momentum equation is solved first, with the resulting fields passed to the energy equation solver. This fits the available data (15 minute interval averages), and allows the slower/ small dt momentum equation to be solved over a shorter period, and the fast / large dt energy equation to be solved over the entire 15 minute period, all while the whole system runs as quickly as possible.
• separate energy equations for fluid and solid. We have good data for the fluid phase (air), but it is the solid phase we really care about. Also, it is possible that the temperature loggers deep in the experimental clamp we had are actually in the solid phase: they were driven in after the construction of the clamp. Over much of the domain, the temperature of the two phases seems very similar.
• inlet and outlet swap. By having a symmetrical domain, changes in the wind direction can be accounted for by just choosing the appropriate mesh for the given 15 minute interval step.
• cover region. This is pretty broad (15 cm) to accommodate a coarse grid. The actual cover is 0.11 cm thick.
• Ts at the clamp edge. To ensure no odd conduction in the solid phase (Ts) at the edge of the clamp region from outside of the clamp, the Ts field in the open region is allowed to follow the Ts in the clamp region.