User Tools

Site Tools


gtspring2009:2dturbulence

Differences

This shows you the differences between two versions of the page.

Link to this comparison view

Both sides previous revision Previous revision
Next revision
Previous revision
gtspring2009:2dturbulence [2009/02/24 08:13]
paprocki
gtspring2009:2dturbulence [2010/02/02 07:55] (current)
Line 5: Line 5:
  
 ====== The experiment ====== ====== The experiment ======
 +The following are recurrence plots for runs with DC current at 3V, 6V, 9V, 12V and 15V respectively:​
 +
 +{{:​gtspring2009:​recplot_3v.jpg|}} ​ {{:​gtspring2009:​recplot_6v.jpg|}} ​ {{:​gtspring2009:​recplot_9v.jpg|}} ​
 +{{:​gtspring2009:​recplot_12v.jpg|}} ​ {{:​gtspring2009:​recplot_15v.jpg|}}
 +
 +The following videos correspond to the above plots: ​
 +
 +3V: [[http://​www.youtube.com/​watch?​v=HRbT2zgkMkg]]
 +6V: [[http://​www.youtube.com/​watch?​v=BtSvv9C9HJI&​feature=channel_page]] ​
 +9V: [[http://​www.youtube.com/​watch?​v=44kXAITPM1A&​feature=channel]] ​
 +12V: [[http://​www.youtube.com/​watch?​v=s2VSBRYrpIU&​feature=channel]] ​
 +15V: [[http://​www.youtube.com/​watch?​v=SanCC1bl6hY&​feature=channel]]
 +
 +
 +
 +{{:​gtspring2009:​paprocki1.jpg?​24}} Here's a description of the experiment, currently in progress:
 +
 +In general, we are investigating what is going on with forced two-dimensional turbulence. ​ We have a thin layer of electrically conducting fluid on top of a bed of strong magnets, with glass separating the fluid and the magnets. ​ Currently, the fluid is salt water. ​ Then, we drive current through the water using two large electrodes hooked up to a power source and a function generator. ​ The only signal we've tested thus far is a sine wave.  The current goal is to find a regime of flow somewhere between periodic and turbulent motion.
 +
 +Once we have enough preliminary data to know what we're actually looking for and what the parameters we'll be using are, we'll switch out the salt water for sulfuric acid, and make sure to keep the molarity consistent. ​ We'll also be using a voltage to current converter so that we'll be keeping the current consistent rather than the voltage, which is what is more important.
 +
 +We have a few different ways to look at data graphically right now.  The first one is the velocity field data.  Using particle image velocimetry software, we plug in a movie of the experiment into MATLAB and it spits out velocity field data, giving us an idea of which areas of the experiment have fluid moving in what direction. ​ From velocity field data, we can make recurrence plots, which are what is plastered all over this blog.  It's a scalar field that takes in the velocity field for two times, T + t and t, and finds the difference. ​ The formula is Sum over velocity field of |U(T+t)-U(t)|. ​ This will give a range of values, which is then assigned to different colors, and these colors are plotted on the recurrence plot.  We can also make recurrence plots for individual vectors in the vector field. ​ Recurrence plots appear to be "​missing"​ the upper triangle of the graph - this is intentional,​ as there isn't actually any data in that region to be looked at.  The horizontal axis is the t axis, and the vertical axis is the T axis.  The units are in seconds. ​ This allows us to spot periodicity better than a formula like |U(T)-U(t)|. ​ The last way we are currently looking at data are time series, which are plots of the magnitude of a given vector or entire vector field over time.
 +
 +//​[[jon.paprocki@gatech.edu|Jon Paprocki]] 2009-03-26 2:18 PM//
  
 {{:​gtspring2009:​paprocki1.jpg?​24}} Explanation on what we're actually doing COMING SOON *cue movie preview music* {{:​gtspring2009:​paprocki1.jpg?​24}} Explanation on what we're actually doing COMING SOON *cue movie preview music*
  
-Short, temporary explanation:​ Put a thin layer of salt water on top of some strong magnets and drive current through it.  Then, subject slave computer to hours of number crunching.+Short, temporary explanation:​ Put a thin layer of salt water on top of some strong magnets and drive current through it.  Then, subject slave computer to hours of number crunching. ​//​[[jon.paprocki@gatech.edu|Jon Paprocki]] 2009-02-24 11:13 AM// 
 + 
 +{{gtspring2009:​pc.jpg?​12 }} Can you start writing up the complete description of the experiment, in preparation for you term project report? Looking at it from Copenhagen I have no clue what is being measured (except that it looks cute)  --- //​[[predrag.cvitanovic@physics.gatech.edu|Predrag Cvitanovic]] 2009-02-11 13:49//
  
-//​[[jon.paprocki@gatech.edu|Jon Paprocki]] 2009-02-24 11:13 AM// 
  
 ===== Movies ===== ===== Movies =====
  
-{{:​gtspring2009:​paprocki1.jpg?​24}} Movies of the first and third recurrence plots for the first search for something between periodic and turbulence motion posted to [[http://​www.youtube.com/​user/​schatzturbulence|Schatz turbulence]]. //​[[jon.paprocki@gatech.edu|Jon Paprocki]] 2009-02-24 11:00 AM//+{{:​gtspring2009:​paprocki1.jpg?​24}} Here are movies of the .0625, .25, and .5 Hz runs mentioned below on the 03/26 post: 
 + 
 +[[http://​www.youtube.com/​watch?​v=DmY8NLv28Mo&​feature=channel_page|.0625 Hz video]] 
 +[[http://​www.youtube.com/​watch?​v=WJ0cWlpfT0c&​feature=channel_page|.0625 Hz velocity fields]] 
 +[[http://​www.youtube.com/​watch?​v=w3At3yzWO8w&​feature=channel_page|.25 Hz video]] 
 +[[http://​www.youtube.com/​watch?​v=dL9NhrgiWTQ&​feature=channel_page|.25 Hz velocity fields]] 
 +[[http://​www.youtube.com/​watch?​v=nanXO1dMXaQ&​feature=channel_page|.5 Hz video]] 
 +[[http://​www.youtube.com/​watch?​v=rr4smnXC65c&​feature=channel_page|.5 Hz velocity fields]] 
 + 
 +After these movies I'm going to switch to uploading them here - I meant to do it for this batch but forgot and uploaded them to youtube out of habit. 
 + 
 +//​[[jon.paprocki@gatech.edu|Jon Paprocki]] 2009-03-26 1:53 PM// 
 + 
 +{{:​gtspring2009:​paprocki1.jpg?​24}} Movies of the first and third recurrence plots for the first search for something between periodic and turbulence motion posted to [[http://​www.youtube.com/​user/​schatzturbulence|Schatz turbulence]]. It looks messy with lots of large conglomerates of glass spheres stuck together. ​ We did 22 runs in a row and didn't want to change out the salt water since that would probably change the conductivity as well.  It is also helpful since it shows that the PIV algorithm doesn'​t completely break when there is a lot of junk floating around.//​[[jon.paprocki@gatech.edu|Jon Paprocki]] 2009-02-24 11:00 AM//
  
 {{:​gtspring2009:​paprocki1.jpg?​24}} Movies of the hybrid runs posted below are now up at [[http://​www.youtube.com/​user/​schatzturbulence|Schatz turbulence]]. ​ //​[[jon.paprocki@gatech.edu|Jon Paprocki]] 2009-02-13 5:27 PM// {{:​gtspring2009:​paprocki1.jpg?​24}} Movies of the hybrid runs posted below are now up at [[http://​www.youtube.com/​user/​schatzturbulence|Schatz turbulence]]. ​ //​[[jon.paprocki@gatech.edu|Jon Paprocki]] 2009-02-13 5:27 PM//
Line 21: Line 59:
  --- //​[[dborrero@gatech.edu|Daniel B.]] 2009-02-03 11:​20// ​  --- //​[[dborrero@gatech.edu|Daniel B.]] 2009-02-03 11:​20// ​
  
-===== Plots =====+===== Plots =====  
 +{{gtspring2009:​paprocki1.jpg?​23}} Things have been quiet on this page lately, but I haven'​t been slacking. ​ We've been searching for a regime between periodic and turbulent flow, and here is the interesting stuff I've found in the last month. 
 + 
 +These are recurrence plots of runs at 15 volts, with frequencies of .0625 Hz, .125 Hz, .25 Hz, .5 Hz, .75 Hz, and 1 Hz respectively. 
 + 
 +{{:​gtspring2009:​research_projects:​schatz_lab:​d0625f_15v.jpg?​220}}{{:​gtspring2009:​research_projects:​schatz_lab:​d125f_15v.jpg?​220}}{{:​gtspring2009:​research_projects:​schatz_lab:​d25f_15v.jpg?​220}}{{:​gtspring2009:​research_projects:​schatz_lab:​d5f_15v.jpg?​220}}{{:​gtspring2009:​research_projects:​schatz_lab:​d75f_15v.jpg?​220}}{{:​gtspring2009:​research_projects:​schatz_lab:​1f_15v.jpg?​220}} 
 + 
 +You can see a very clear jump from vaguely semi-periodic to completely periodic motion between .25 and .5 Hz.  .0625, .125, and .25 Hz all appear to be just plain turbulent to the naked eye, but the recurrence data suggests otherwise. 
 + 
 +Also, I'm giving a talk on this project for the undergrad research symposium on April 1st in room 320 of the student center at 1:40 PM, if you're interested in something more detailed than this ramshackle blog. 
 + 
 +EDIT: movies of these posted above in movies section 
 + 
 +//​[[jon.paprocki@gatech.edu|Jon Paprocki]] 2009-03-26 11:50 AM// 
 + 
 +{{:​gtspring2009:​paprocki1.jpg?​24}} Recurrence plots for three higher voltages, with corresponding time series for the entire vector field. ​ These are also at a new, slightly higher frequency. ​ I don't know what the frequency was at before, but now the period of forcing is 3.92 seconds (for an entire cycle, so thats the time to go from +X to -X and back to +X). 
 + 
 +{{:​gtspring2009:​research_projects:​schatz_lab:​search23.jpg?​220}} {{:​gtspring2009:​research_projects:​schatz_lab:​search24.jpg?​220}} {{:​gtspring2009:​research_projects:​schatz_lab:​search25.jpg?​220}} \\ 
 + 
 +{{:​gtspring2009:​research_projects:​schatz_lab:​search23time.jpg?​220}}  
 +{{:​gtspring2009:​research_projects:​schatz_lab:​search24time.jpg?​220}}  
 +{{:​gtspring2009:​research_projects:​schatz_lab:​search25time.jpg?​220}} \\ 
 {{:​gtspring2009:​gibson.png?​24}} Gotcha. That makes sense. If you expand out %%(sin(t+T) - sin(t))^2%% you'll have some %%cos(t)%% and %%sin(t)%% mutliplicative factors that'​ll make r(t,T) oscillate in t.  {{:​gtspring2009:​gibson.png?​24}} Gotcha. That makes sense. If you expand out %%(sin(t+T) - sin(t))^2%% you'll have some %%cos(t)%% and %%sin(t)%% mutliplicative factors that'​ll make r(t,T) oscillate in t. 
  
Line 136: Line 196:
 --pivDifArray.m -> takes in a correlation.mat file and outputs a difArray.mat that essentially is just converting the structure array of correlation.mat into a more usable array containing only the velocity data.  This makes it easier to use vectorized code.  This isn't the original use of the function so it isn't very aptly named...\\ ​ --pivDifArray.m -> takes in a correlation.mat file and outputs a difArray.mat that essentially is just converting the structure array of correlation.mat into a more usable array containing only the velocity data.  This makes it easier to use vectorized code.  This isn't the original use of the function so it isn't very aptly named...\\ ​
 ---pivRecurrencePlot.m -> takes in a difArray.mat file and outputs a recurrencePlot.mat file containing recurrence data.\\ ​ ---pivRecurrencePlot.m -> takes in a difArray.mat file and outputs a recurrencePlot.mat file containing recurrence data.\\ ​
 +---PivRecurrencePlotNew.m -> uses the new formula for recurrence plots (doesn'​t square the difference vectors); outputs a recurrencePlotNew.mat file.\\
 +----piv1DRecSlicet.m -> grabs a "​slice"​ of the recurrence plot at a given time along the t axis and plots it.\\
 +----piv1DRecSliceTT.m -> grabs a slice of the recurrence plot at a given time along the T axis and plots it.\\
 +----minFinderY.m -> "​filters"​ out a recurrence plot by blotting out all the points that aren't minima along the Y axis, to some sensitivity.\\
 +----minFinderX.m -> same as minFinderY, but along the X axis.\\
 ---pivRecurrencePlotOne.m -> same as pivRecurrence,​ but only does it for one specified point and ouputs a recurrenceDataLoc###​.mat file\\ ​ ---pivRecurrencePlotOne.m -> same as pivRecurrence,​ but only does it for one specified point and ouputs a recurrenceDataLoc###​.mat file\\ ​
 ----pivRecurrenceImage.m -> takes in a recurrencePlot.mat or recurrenceDataLoc###​.mat file and outputs an image of it.\\  ----pivRecurrenceImage.m -> takes in a recurrencePlot.mat or recurrenceDataLoc###​.mat file and outputs an image of it.\\ 
Line 143: Line 208:
 8-o Checked out Dustin and Jon Paprocki'​s permission problems on the wiki. Enabled gt2009 permissions for Jon P.; 8-o Checked out Dustin and Jon Paprocki'​s permission problems on the wiki. Enabled gt2009 permissions for Jon P.;
 you should be able to start your own blog now. Jon P, your name is correct, but do you like Zappa?// John Gibson 2009-02-04 12:12// you should be able to start your own blog now. Jon P, your name is correct, but do you like Zappa?// John Gibson 2009-02-04 12:12//
- 
-{{:​gtspring2009:​paprocki1.jpg?​24}} To clarify, I can read/edit a few random pages on the wiki, but a lot of things I do not appear to have access to (such as starting my own mini-blog). //​[[jon.paprocki@gatech.edu|Jon Paprocki]] 2009-02-03 11:33// 
  
 :-P John G. could you please give Jon Paprocki (he's the undergrad working on this project) access to the blog? Thanks -- //​[[dborrero@gatech.edu|Daniel B.]] 2009-02-03 11:31// :-P John G. could you please give Jon Paprocki (he's the undergrad working on this project) access to the blog? Thanks -- //​[[dborrero@gatech.edu|Daniel B.]] 2009-02-03 11:31//
gtspring2009/2dturbulence.1235492027.txt.gz · Last modified: 2009/02/24 08:13 by paprocki