{"id":5518,"date":"2016-12-13T14:55:48","date_gmt":"2016-12-13T19:55:48","guid":{"rendered":"http:\/\/pages.vassar.edu\/magnes\/?p=5518"},"modified":"2016-12-13T14:55:48","modified_gmt":"2016-12-13T19:55:48","slug":"viscous-drag-and-motion","status":"publish","type":"post","link":"https:\/\/pages.vassar.edu\/magnes\/2016\/12\/13\/viscous-drag-and-motion\/","title":{"rendered":"Viscous Drag and Motion"},"content":{"rendered":"<p><b>Overview<\/b><\/p>\n<p><span style=\"font-weight: 400\">This project aimed to simulate the viscous drag forces that occur when a ball falls in fluids with different densities and viscosities, resulting in different Reynolds numbers. In this case, there are two graphs produced. The first shows the ball falling in the medium in phase-space. The second shows the vertical distance traveled by the ball in time. The \u00a0ball falling in the fluid with the larger Reynold number should take a \u00a0longer time to fall because the drag forces on the ball are greater. Fluids include both liquids and gases. <\/span><\/p>\n<p><span style=\"font-weight: 400\">Reynolds numbers determine the type of flow the fluid will make. When Re&lt;&lt;1 the fluid follows Stoke&#8217;s Law. When Re&gt;&gt;1 the fluid does not,it follows the Navier-Stokes equations. The Reynolds number is a proportion of the inertial forces to the viscous forces, Re = inertial\/viscous.<\/span><\/p>\n<p><b>Motivations<\/b><\/p>\n<p><span style=\"font-weight: 400\">While I studied abroad, I took a fluid dynamics class that covered the Navier-Stokes equations and complicated geometries of objects moving in fluids. The most basic geometry is a sphere in 3D and a circle\/disk in 2D. I wanted to be able to visualize and create a code that could interpret various inputs for fluids and produce a 2D plot of the motion. It proved more difficult to do than initially thought because I do not have a solid background in fluid dynamics, just that one course. As a result, I referenced derivations for the Navier-Stokes, but most were more complicated than the 2D system I wanted to model. I referenced my notes from the course too and <\/span><a href=\"https:\/\/youtu.be\/51-6QCJTAjU?t=17m39s\"><span style=\"font-weight: 400\">YouTube videos from MIT<\/span><\/a><span style=\"font-weight: 400\">. <\/span><\/p>\n<p><b>Figures (click on the images to make them larger)<\/b><\/p>\n<p><a href=\"http:\/\/pages.vassar.edu\/magnes\/files\/2016\/12\/water_vs_honey.jpg\"><img loading=\"lazy\" decoding=\"async\" class=\"aligncenter wp-image-5519\" src=\"http:\/\/pages.vassar.edu\/magnes\/files\/2016\/12\/water_vs_honey-300x149.jpg\" alt=\"water_vs_honey\" width=\"615\" height=\"305\" srcset=\"https:\/\/pages.vassar.edu\/magnes\/files\/2016\/12\/water_vs_honey-300x149.jpg 300w, https:\/\/pages.vassar.edu\/magnes\/files\/2016\/12\/water_vs_honey-768x381.jpg 768w, https:\/\/pages.vassar.edu\/magnes\/files\/2016\/12\/water_vs_honey-1024x507.jpg 1024w, https:\/\/pages.vassar.edu\/magnes\/files\/2016\/12\/water_vs_honey-624x309.jpg 624w\" sizes=\"auto, (max-width: 615px) 100vw, 615px\" \/><\/a><\/p>\n<p><span style=\"font-weight: 400\">Figure 1. The parameters of this plot are those of water and honey at room temperature (20\u00b0 C). The viscosities are 1.0 (Pa s) and 1.36 (Pa s) and the densities are 0.7978 (kg\/L) and 7.20 (kg\/L) respectively. The ball travels slower in the honey, Re<\/span><span style=\"font-weight: 400\">2<\/span><span style=\"font-weight: 400\">. <\/span><\/p>\n<p><a href=\"http:\/\/pages.vassar.edu\/magnes\/files\/2016\/12\/randomvalues.jpg\"><img loading=\"lazy\" decoding=\"async\" class=\"aligncenter wp-image-5520\" src=\"http:\/\/pages.vassar.edu\/magnes\/files\/2016\/12\/randomvalues-300x177.jpg\" alt=\"randomvalues\" width=\"615\" height=\"363\" srcset=\"https:\/\/pages.vassar.edu\/magnes\/files\/2016\/12\/randomvalues-300x177.jpg 300w, https:\/\/pages.vassar.edu\/magnes\/files\/2016\/12\/randomvalues-768x453.jpg 768w, https:\/\/pages.vassar.edu\/magnes\/files\/2016\/12\/randomvalues-1024x604.jpg 1024w, https:\/\/pages.vassar.edu\/magnes\/files\/2016\/12\/randomvalues-624x368.jpg 624w\" sizes=\"auto, (max-width: 615px) 100vw, 615px\" \/><\/a><\/p>\n<p><span style=\"font-weight: 400\">Figure 2. The parameters for this plot are viscosities of 32 (Pa s) \u00a0and 400 (Pa s), densities of 20 (kg\/L) and 75 (kg\/L), and kinematic viscosities of 1.6 and 5.33, respectively. The Reynolds numbers are displayed on the plots. The ball moves slower in the larger Reynolds number fluid.<\/span><\/p>\n<p><b>Fluid Mechanics Physics<\/b><\/p>\n<p><a href=\"http:\/\/pages.vassar.edu\/magnes\/files\/2016\/12\/Screen-Shot-2016-12-13-at-14.44.04.png\"><img loading=\"lazy\" decoding=\"async\" class=\"aligncenter wp-image-5522 size-large\" src=\"http:\/\/pages.vassar.edu\/magnes\/files\/2016\/12\/Screen-Shot-2016-12-13-at-14.44.04-1024x164.png\" alt=\"screen-shot-2016-12-13-at-14-44-04\" width=\"625\" height=\"100\" srcset=\"https:\/\/pages.vassar.edu\/magnes\/files\/2016\/12\/Screen-Shot-2016-12-13-at-14.44.04-1024x164.png 1024w, https:\/\/pages.vassar.edu\/magnes\/files\/2016\/12\/Screen-Shot-2016-12-13-at-14.44.04-300x48.png 300w, https:\/\/pages.vassar.edu\/magnes\/files\/2016\/12\/Screen-Shot-2016-12-13-at-14.44.04-768x123.png 768w, https:\/\/pages.vassar.edu\/magnes\/files\/2016\/12\/Screen-Shot-2016-12-13-at-14.44.04-624x100.png 624w, https:\/\/pages.vassar.edu\/magnes\/files\/2016\/12\/Screen-Shot-2016-12-13-at-14.44.04.png 1806w\" sizes=\"auto, (max-width: 625px) 100vw, 625px\" \/><\/a><\/p>\n<p><span style=\"font-weight: 400\">The above equations are the Navier-Stokes equations in 3D for cartesian coordinates. For the 2D version that I used, the z-direction components are removed. These equations are only used for when Re&gt;&gt;1, for high Reynolds numbers. Because the Navier-Stokes equations can only be solved for very specific cases I simplified the process by using the drag equation, <\/span><\/p>\n<p><a href=\"http:\/\/pages.vassar.edu\/magnes\/files\/2016\/12\/Screen-Shot-2016-12-13-at-14.45.17.png\"><img loading=\"lazy\" decoding=\"async\" class=\"alignnone size-full wp-image-5523\" src=\"http:\/\/pages.vassar.edu\/magnes\/files\/2016\/12\/Screen-Shot-2016-12-13-at-14.45.17.png\" alt=\"screen-shot-2016-12-13-at-14-45-17\" width=\"284\" height=\"52\" \/><\/a><\/p>\n<p><span style=\"font-weight: 400\">and then solved Newton\u2019s second law for the vertical and horizontal directions, and formulated the kinematics equations. Buoyancy forces were ignored. This is a different case than than when Re&lt;&lt;1. <\/span><\/p>\n<p><span style=\"font-weight: 400\">When Re&lt;&lt;1, the system follows viscous drag using Stoke\u2019s Law<\/span><\/p>\n<p><a href=\"http:\/\/pages.vassar.edu\/magnes\/files\/2016\/12\/Screen-Shot-2016-12-13-at-14.46.07.png\"><img loading=\"lazy\" decoding=\"async\" class=\"alignnone size-full wp-image-5524\" src=\"http:\/\/pages.vassar.edu\/magnes\/files\/2016\/12\/Screen-Shot-2016-12-13-at-14.46.07.png\" alt=\"screen-shot-2016-12-13-at-14-46-07\" width=\"192\" height=\"36\" \/><\/a><\/p>\n<p><span style=\"font-weight: 400\">From here you can solve Newton\u2019s second law and formulate the kinematic equations again. Stoke\u2019s Law only works for these cases, otherwise Navier-Stokes and the drag equation are used. Buoyancy forces were ignored here too.\u00a0<\/span><\/p>\n<p><b>Results<\/b><\/p>\n<p><span style=\"font-weight: 400\">Computational fluid dynamics (CFD) is a growing and complicated field of both mathematics and physics. The program I developed is part of the general sphere of CFD but not close to the type of programs are used and developed. CFD is used to simulate airplane designs, model the way fluids flow around objects. My program has some bugs because of the simplifications I used for the equations. By using an adapted version of the Navier-Stokes equations, not all parameter values work correctly with the model. For the many combinations of inputs, however, the larger Reynolds number fluid produces a slower movement in the fluid. <\/span><\/p>\n<p><b>Conclusion<\/b><\/p>\n<p><span style=\"font-weight: 400\">In the future, I would like to fix the bugs caused by the simplifications and possibly create a program that uses the Navier-Stokes equation in its full 3D form. If possible, it would be useful to have a GUI for the program because that would allow the user to more easily see the differences in the motion when the viscosities and densities of the fluids are changed. Perhaps adding more than two fluids to the model as well as making the visualization as a movie and not stagnant plots.<\/span><\/p>\n<p><b>References<\/b><\/p>\n<p>http:\/\/soliton.ae.gatech.edu\/labs\/windtunl\/classes\/hispd\/hispd06\/ns_eqns.html<\/p>\n<p>http:\/\/www.personal.psu.edu\/wzl113\/Lesson%20Plan.htm#GF<\/p>\n<p>http:\/\/scienceworld.wolfram.com\/physics\/StokesVelocity.html<\/p>\n<p>https:\/\/www.grc.nasa.gov\/www\/k-12\/airplane\/nseqs.html<\/p>\n<p>https:\/\/www.grc.nasa.gov\/www\/k-12\/airplane\/reynolds.html<\/p>\n","protected":false},"excerpt":{"rendered":"<p>Overview This project aimed to simulate the viscous drag forces that occur when a ball falls in fluids with different densities and viscosities, resulting in different Reynolds numbers. In this case, there are two graphs produced. The first shows the ball falling in the medium in phase-space. The second shows the vertical distance traveled by [&hellip;]<\/p>\n","protected":false},"author":3475,"featured_media":0,"comment_status":"open","ping_status":"closed","sticky":false,"template":"","format":"standard","meta":{"footnotes":""},"categories":[128934],"tags":[],"class_list":["post-5518","post","type-post","status-publish","format-standard","hentry","category-molly"],"_links":{"self":[{"href":"https:\/\/pages.vassar.edu\/magnes\/wp-json\/wp\/v2\/posts\/5518","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/pages.vassar.edu\/magnes\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/pages.vassar.edu\/magnes\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/pages.vassar.edu\/magnes\/wp-json\/wp\/v2\/users\/3475"}],"replies":[{"embeddable":true,"href":"https:\/\/pages.vassar.edu\/magnes\/wp-json\/wp\/v2\/comments?post=5518"}],"version-history":[{"count":4,"href":"https:\/\/pages.vassar.edu\/magnes\/wp-json\/wp\/v2\/posts\/5518\/revisions"}],"predecessor-version":[{"id":5527,"href":"https:\/\/pages.vassar.edu\/magnes\/wp-json\/wp\/v2\/posts\/5518\/revisions\/5527"}],"wp:attachment":[{"href":"https:\/\/pages.vassar.edu\/magnes\/wp-json\/wp\/v2\/media?parent=5518"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/pages.vassar.edu\/magnes\/wp-json\/wp\/v2\/categories?post=5518"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/pages.vassar.edu\/magnes\/wp-json\/wp\/v2\/tags?post=5518"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}