{"id":5519,"date":"2020-12-18T05:59:57","date_gmt":"2020-12-18T10:59:57","guid":{"rendered":"https:\/\/blogs.sw.siemens.com\/simulating-the-real-world\/?p=5519"},"modified":"2026-03-27T09:00:31","modified_gmt":"2026-03-27T13:00:31","slug":"chaotic-fluid-dynamics-part-3-wheels","status":"publish","type":"post","link":"https:\/\/blogs.sw.siemens.com\/simulating-the-real-world\/2020\/12\/18\/chaotic-fluid-dynamics-part-3-wheels\/","title":{"rendered":"Chaotic Fluid Dynamics Part 3 \u2013 Wheels"},"content":{"rendered":"\n<p>Rayleigh\u2013B\u00e9nard convection involves the transport of heat from a hot bottom surface to a cold top surface. If the temperature difference between top and bottom is big enough, fluid rolls up into counter-rotating convection cells that transport the heat. Under certain conditions these convection cells flip-flop chaotically between rotating in a clockwise and counter-clockwise directions. Ed Lorenz&#8217;s famous &#8216;strange attractor&#8217; is a graphical description of this behaviour:<\/p>\n\n\n<div class=\"wp-block-image\">\n<figure class=\"aligncenter size-large is-resized\"><img decoding=\"async\" src=\"https:\/\/upload.wikimedia.org\/wikipedia\/commons\/1\/13\/A_Trajectory_Through_Phase_Space_in_a_Lorenz_Attractor.gif\" alt=\"\" style=\"width:270px;height:270px\"\/><figcaption class=\"wp-element-caption\">Lorenz Attractor. Convection intensity vs Temperature variation vs Temperature difference for Rayleigh\u2013B\u00e9nard convection when Prandtl number = 10, Rayleigh number is 28 times greater than the critical value and aspect ratio of model is 2:1.<\/figcaption><\/figure><\/div>\n\n\n<h2 class=\"wp-block-heading has-large-font-size\">A Pedagogic Example<\/h2>\n\n\n\n<p>The <a href=\"https:\/\/en.wikipedia.org\/wiki\/Malkus_waterwheel\" target=\"_blank\" rel=\"noreferrer noopener\">Malkus waterwheel<\/a> is a fluid-mechanical model that intends to capture the same physical mechanisms as Rayleigh\u2013B\u00e9nard convection. A tap\/faucet supplies water at the top of the wheel, filling freely swinging cups positioned on the wheel periphery. Each cup has a hole in the bottom allow the liquid to slowly leak out:<\/p>\n\n\n\n<figure class=\"wp-block-video aligncenter\"><video controls src=\"https:\/\/videos.mentor-cdn.com\/mgc\/videos\/2500\/4b7bb472-96eb-4067-837c-908f3f6322fc-en-US-video.mp4\"><\/video><figcaption class=\"wp-element-caption\">Malkus waterwheel schematic adapted from https:\/\/en.wikipedia.org\/wiki\/Malkus_waterwheel<\/figcaption><\/figure>\n\n\n\n<p>As the cups are filled at the top, the wheel starts to move. If it moves too quickly, heavy cups rotating upwards tend the slow the wheel down. This results in cups being filled with more water at the top that then speeds the wheel up. The wheel flips between spinning in a clockwise and counter-clockwise direction. The red curve in the above animation is the center of gravity of the wheel, its path describes an attractor of the system.<\/p>\n\n\n\n<p>As a pedagogic aid, the Malkus waterwheel is often used as an example of a chaotic system. From a thermo-fluids perspective, the added water is akin to &#8216;adding cold&#8217; at the top, and the holes in the cups that let that water out can be viewed as &#8216;adding heat&#8217;.<\/p>\n\n\n\n<p>The Double Pendulum is also used as an example of what, on the surface appears a quite simple momentum based system, exhibits quite beautiful chaotic behaviour.<\/p>\n\n\n\n<figure class=\"wp-block-embed is-type-video is-provider-youtube wp-block-embed-youtube wp-embed-aspect-4-3 wp-has-aspect-ratio\"><div class=\"wp-block-embed__wrapper\">\n<iframe loading=\"lazy\" title=\"Double Pendulum Animation\" width=\"640\" height=\"480\" src=\"https:\/\/www.youtube.com\/embed\/BLsQexP_UYE?feature=oembed\" frameborder=\"0\" allow=\"accelerometer; autoplay; clipboard-write; encrypted-media; gyroscope; picture-in-picture; web-share\" referrerpolicy=\"strict-origin-when-cross-origin\" allowfullscreen><\/iframe>\n<\/div><\/figure>\n\n\n\n<h2 class=\"wp-block-heading has-large-font-size\">A Thermosyphon Equivalent<\/h2>\n\n\n\n<p>So can CFD to be used to model a Malkus waterwheel? It&#8217;s quite a complex fluid-solid interaction application involving filling, sloshing and rotational dynamics, something that <a href=\"https:\/\/www.plm.automation.siemens.com\/global\/en\/products\/simcenter\/STAR-CCM.html\" target=\"_blank\" rel=\"noreferrer noopener\">Simcenter STAR-CCM+<\/a> would be well suited to. However I&#8217;m going to take a slightly different approach, further abstract the application into an equivalent closed loop thermosyphon and use <a href=\"https:\/\/www.plm.automation.siemens.com\/global\/en\/products\/simcenter\/floefd-siemens-nx.html\" target=\"_blank\" rel=\"noreferrer noopener\">Simcenter FLOEFD for NX<\/a> to simulate its behaviour:<\/p>\n\n\n<div class=\"wp-block-image\">\n<figure class=\"aligncenter size-large\"><img loading=\"lazy\" decoding=\"async\" width=\"668\" height=\"576\" src=\"https:\/\/blogs.sw.siemens.com\/wp-content\/uploads\/sites\/56\/2020\/12\/SimcetnerFLOEFD_Waterwheel.png\" alt=\"\" class=\"wp-image-5525\" srcset=\"https:\/\/blogs.sw.siemens.com\/wp-content\/uploads\/sites\/56\/2020\/12\/SimcetnerFLOEFD_Waterwheel.png 668w, https:\/\/blogs.sw.siemens.com\/wp-content\/uploads\/sites\/56\/2020\/12\/SimcetnerFLOEFD_Waterwheel-600x517.png 600w\" sizes=\"auto, (max-width: 668px) 100vw, 668px\" \/><figcaption class=\"wp-element-caption\">Simcenter FLOEFD for NX model of a closed-loop thermosyphon<\/figcaption><\/figure><\/div>\n\n\n<p>The loop is filled with water. Heat is extracted at a constant rate from a small volume of the fluid at the top of the loop (equivalent to water being added at the top of the Malkus waterwheel). Heat is added at a constant rate through the rest of the wheel (equivalent to the holes in the cups that let water out).<\/p>\n\n\n\n<p>In which direction will the water flow? Let&#8217;s find out! (Spoiler alert &#8211; BOTH)<\/p>\n\n\n\n<figure class=\"wp-block-video aligncenter\"><video controls src=\"https:\/\/videos.mentor-cdn.com\/mgc\/videos\/5400\/d0ab5bab-79d8-406d-b8c9-c8d7b0402bb1-en-US-video.mp4\"><\/video><figcaption class=\"wp-element-caption\">Temperature (red hot, blue cold) and Velocity over time. The graph shows the mass flow rate at the bottom of the tube; positive = clockwise, negative = counter-clockwise<\/figcaption><\/figure>\n\n\n\n<p>Water at the top is cooled down due to the heat extraction, water everywhere else heats up. The cold water drops down the inside of the tube until such time as it forces the flow in one direction, reinforced by the heated water pushing upwards on the other side. Round it spins until it slows down due to cold water being pushed up too far, the flow stalls, temperature drops quickly at the top (e.g. cup is filled with more water at the top when the wheel slows down), then it starts spinning again in a new direction. Sometimes the same direction, sometimes the opposite direction.<\/p>\n\n\n\n<p>The flow is chaotic in that the number of spins in a given direction, and how often it changes rotation direction, are aperiodic:<\/p>\n\n\n\n<figure class=\"wp-block-video aligncenter\"><video controls src=\"https:\/\/videos.mentor-cdn.com\/mgc\/videos\/5400\/5777ce3c-d217-40e5-a1d2-6116e97e0b63-en-US-video.mp4\"><\/video><figcaption class=\"wp-element-caption\">Aperiodic rotational behaviour<\/figcaption><\/figure>\n\n\n\n<h2 class=\"wp-block-heading has-large-font-size\">Strangely Attractive<\/h2>\n\n\n\n<p>An attractor, plotted in so called &#8216;phase space&#8217;, is a way of collapsing a description of the behaviour of a system down into a bounded form. More useful than a time series plot that would go on for ever. It involves some integral properties of the system plotted against each other. Each point equating to a moment in time, joining the points up then shows a line that scribes a path.<\/p>\n\n\n\n<p>For this thermosyphon application, plotting the mass flow rate at the bottom of the pipe against the bulk average temperature at that plane creates a somewhat angry looking attractor:<\/p>\n\n\n<div class=\"wp-block-image\">\n<figure class=\"aligncenter size-large\"><img loading=\"lazy\" decoding=\"async\" width=\"1024\" height=\"683\" src=\"https:\/\/blogs.sw.siemens.com\/wp-content\/uploads\/sites\/56\/2020\/12\/ThermosyphonAttractor-1024x683.png\" alt=\"\" class=\"wp-image-5531\" srcset=\"https:\/\/blogs.sw.siemens.com\/wp-content\/uploads\/sites\/56\/2020\/12\/ThermosyphonAttractor-1024x683.png 1024w, https:\/\/blogs.sw.siemens.com\/wp-content\/uploads\/sites\/56\/2020\/12\/ThermosyphonAttractor-600x400.png 600w, https:\/\/blogs.sw.siemens.com\/wp-content\/uploads\/sites\/56\/2020\/12\/ThermosyphonAttractor-768x512.png 768w, https:\/\/blogs.sw.siemens.com\/wp-content\/uploads\/sites\/56\/2020\/12\/ThermosyphonAttractor-900x600.png 900w, https:\/\/blogs.sw.siemens.com\/wp-content\/uploads\/sites\/56\/2020\/12\/ThermosyphonAttractor.png 1197w\" sizes=\"auto, (max-width: 1024px) 100vw, 1024px\" \/><figcaption class=\"wp-element-caption\">Thermosyphon Attractor, mass flow rate vs. bulk average temperature<\/figcaption><\/figure><\/div>\n\n\n<p>Adding a 3rd property, the maximum temperature difference in the system at each time point, creates this 3D attractor:<\/p>\n\n\n\n<figure class=\"wp-block-video aligncenter\"><video controls src=\"https:\/\/videos.mentor-cdn.com\/mgc\/videos\/2500\/d03dc968-910b-4cbc-a8ea-619be5e78ba9-en-US-video.mp4\"><\/video><figcaption class=\"wp-element-caption\">3D Thermosyphon attractor, mass flow rate vs. temperature vs. maximum temperature difference<\/figcaption><\/figure>\n\n\n\n<p>If the system exhibited a fixed steady state behaviour its attractor would simply be a single point. Periodic time varying behaviour would be some kind of loop. An attractor of a chaotic system is called a &#8216;strange attractor&#8217; in that the line never crosses itself, it keeps on going for ever. If it did intersect itself then it would have achieved an identical state to as it did previously and so go on to repeat that same (periodic) behaviour. A strange attractor is a fractal in that its dimension is fractional. A loop would have a dimension of 1. The Lorenz attractor has a dimension of 2.06, somewhere between a plane and a 3D solid and nope, I still haven&#8217;t got my head around that either.  <\/p>\n\n\n\n<p>Despite the fact that a meshed based CFD approach has been used to derive the above attractor, and it is of a system that is twice abstracted from classic Rayleigh\u2013B\u00e9nard convection, it is heartening to see the similarity between it and the Lorenz attractor. Similar attractors can be observed in a number of differing fields, all going to show that there is a commonality to chaotic behaviour that underpins all physics (because there is only one physics, despite what you might have heard elsewhere).<\/p>\n\n\n\n<p><strong>Chaotic Fluid Dynamics Blog Series:<\/strong><\/p>\n\n\n\n<p><blockquote class=\"wp-embedded-content\" data-secret=\"A9XPkAnAZu\"><a href=\"https:\/\/blogs.sw.siemens.com\/simulating-the-real-world\/2020\/12\/08\/chaotic-fluid-dynamics-part-1-rabbits\/\">Chaotic Fluid Dynamics Part 1 &#8211; Rabbits<\/a><\/blockquote><iframe loading=\"lazy\" class=\"wp-embedded-content\" sandbox=\"allow-scripts\" security=\"restricted\" style=\"position: absolute; visibility: hidden;\" title=\"&#8220;Chaotic Fluid Dynamics Part 1 &#8211; Rabbits&#8221; &#8212; Simulating the Real World\" src=\"https:\/\/blogs.sw.siemens.com\/simulating-the-real-world\/2020\/12\/08\/chaotic-fluid-dynamics-part-1-rabbits\/embed\/#?secret=8evC5GQu6L#?secret=A9XPkAnAZu\" data-secret=\"A9XPkAnAZu\" width=\"600\" height=\"338\" frameborder=\"0\" marginwidth=\"0\" marginheight=\"0\" scrolling=\"no\"><\/iframe><\/p>\n\n\n\n<p><blockquote class=\"wp-embedded-content\" data-secret=\"TPOKjLps8A\"><a href=\"https:\/\/blogs.sw.siemens.com\/simulating-the-real-world\/2020\/12\/10\/chaotic-fluid-dynamics-part-2-butterflies\/\">Chaotic Fluid Dynamics Part 2 \u2013 Butterflies<\/a><\/blockquote><iframe loading=\"lazy\" class=\"wp-embedded-content\" sandbox=\"allow-scripts\" security=\"restricted\" style=\"position: absolute; visibility: hidden;\" title=\"&#8220;Chaotic Fluid Dynamics Part 2 \u2013 Butterflies&#8221; &#8212; Simulating the Real World\" src=\"https:\/\/blogs.sw.siemens.com\/simulating-the-real-world\/2020\/12\/10\/chaotic-fluid-dynamics-part-2-butterflies\/embed\/#?secret=pr82BMd5xq#?secret=TPOKjLps8A\" data-secret=\"TPOKjLps8A\" width=\"600\" height=\"338\" frameborder=\"0\" marginwidth=\"0\" marginheight=\"0\" scrolling=\"no\"><\/iframe><\/p>\n\n\n\n<p><blockquote class=\"wp-embedded-content\" data-secret=\"dkJZ2W9Fws\"><a href=\"https:\/\/blogs.sw.siemens.com\/simulating-the-real-world\/2020\/12\/18\/chaotic-fluid-dynamics-part-3-wheels\/\">Chaotic Fluid Dynamics Part 3 \u2013 Wheels<\/a><\/blockquote><iframe loading=\"lazy\" class=\"wp-embedded-content\" sandbox=\"allow-scripts\" security=\"restricted\" style=\"position: absolute; visibility: hidden;\" title=\"&#8220;Chaotic Fluid Dynamics Part 3 \u2013 Wheels&#8221; &#8212; Simulating the Real World\" src=\"https:\/\/blogs.sw.siemens.com\/simulating-the-real-world\/2020\/12\/18\/chaotic-fluid-dynamics-part-3-wheels\/embed\/#?secret=OsnkwEY172#?secret=dkJZ2W9Fws\" data-secret=\"dkJZ2W9Fws\" width=\"600\" height=\"338\" frameborder=\"0\" marginwidth=\"0\" marginheight=\"0\" scrolling=\"no\"><\/iframe><\/p>\n\n\n\n<p><blockquote class=\"wp-embedded-content\" data-secret=\"EQqU1kmXff\"><a href=\"https:\/\/blogs.sw.siemens.com\/simulating-the-real-world\/2021\/01\/04\/chaotic-fluid-dynamics-part-4-finding-feigenbaum\/\">Chaotic Fluid Dynamics Part 4 &#8211; Finding Feigenbaum<\/a><\/blockquote><iframe loading=\"lazy\" class=\"wp-embedded-content\" sandbox=\"allow-scripts\" security=\"restricted\" style=\"position: absolute; visibility: hidden;\" title=\"&#8220;Chaotic Fluid Dynamics Part 4 &#8211; Finding Feigenbaum&#8221; &#8212; Simulating the Real World\" src=\"https:\/\/blogs.sw.siemens.com\/simulating-the-real-world\/2021\/01\/04\/chaotic-fluid-dynamics-part-4-finding-feigenbaum\/embed\/#?secret=pwZilcrSai#?secret=EQqU1kmXff\" data-secret=\"EQqU1kmXff\" width=\"600\" height=\"338\" frameborder=\"0\" marginwidth=\"0\" marginheight=\"0\" scrolling=\"no\"><\/iframe><\/p>\n","protected":false},"excerpt":{"rendered":"<p>Rayleigh\u2013B\u00e9nard convection involves the transport of heat from a hot bottom surface to a cold top surface. If the temperature&#8230;<\/p>\n","protected":false},"author":2666,"featured_media":5520,"comment_status":"open","ping_status":"open","sticky":false,"template":"","format":"standard","meta":{"spanish_translation":"","french_translation":"","german_translation":"","italian_translation":"","polish_translation":"","japanese_translation":"","chinese_translation":"","footnotes":""},"categories":[1],"tags":[537,534,536],"industry":[],"product":[],"coauthors":[542],"class_list":["post-5519","post","type-post","status-publish","format-standard","has-post-thumbnail","hentry","category-news","tag-chaos","tag-simcenter-floefd","tag-thermosyphon"],"featured_image_url":"https:\/\/blogs.sw.siemens.com\/wp-content\/uploads\/sites\/56\/2020\/12\/Wheel1.png","_links":{"self":[{"href":"https:\/\/blogs.sw.siemens.com\/simulating-the-real-world\/wp-json\/wp\/v2\/posts\/5519","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/blogs.sw.siemens.com\/simulating-the-real-world\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/blogs.sw.siemens.com\/simulating-the-real-world\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/blogs.sw.siemens.com\/simulating-the-real-world\/wp-json\/wp\/v2\/users\/2666"}],"replies":[{"embeddable":true,"href":"https:\/\/blogs.sw.siemens.com\/simulating-the-real-world\/wp-json\/wp\/v2\/comments?post=5519"}],"version-history":[{"count":4,"href":"https:\/\/blogs.sw.siemens.com\/simulating-the-real-world\/wp-json\/wp\/v2\/posts\/5519\/revisions"}],"predecessor-version":[{"id":6187,"href":"https:\/\/blogs.sw.siemens.com\/simulating-the-real-world\/wp-json\/wp\/v2\/posts\/5519\/revisions\/6187"}],"wp:featuredmedia":[{"embeddable":true,"href":"https:\/\/blogs.sw.siemens.com\/simulating-the-real-world\/wp-json\/wp\/v2\/media\/5520"}],"wp:attachment":[{"href":"https:\/\/blogs.sw.siemens.com\/simulating-the-real-world\/wp-json\/wp\/v2\/media?parent=5519"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/blogs.sw.siemens.com\/simulating-the-real-world\/wp-json\/wp\/v2\/categories?post=5519"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/blogs.sw.siemens.com\/simulating-the-real-world\/wp-json\/wp\/v2\/tags?post=5519"},{"taxonomy":"industry","embeddable":true,"href":"https:\/\/blogs.sw.siemens.com\/simulating-the-real-world\/wp-json\/wp\/v2\/industry?post=5519"},{"taxonomy":"product","embeddable":true,"href":"https:\/\/blogs.sw.siemens.com\/simulating-the-real-world\/wp-json\/wp\/v2\/product?post=5519"},{"taxonomy":"author","embeddable":true,"href":"https:\/\/blogs.sw.siemens.com\/simulating-the-real-world\/wp-json\/wp\/v2\/coauthors?post=5519"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}