arxiv: v1 [physics.ed-ph] 24 Apr 2018

Transcription

1 Computer simulations combined with experiments for calculus-based Physics laboratory course arxiv: v1 [physics.ed-ph] 24 Apr 2018 Sergey V. Samsonau a Princeton International School of Mathematics and Science, Princeton, NJ, USA Keywords: Physics, Education, Laboratory, Experiment, Simulation 1. Abstract This paper is presenting a set oflaboratoryclassesto be taught as a part ofa 1-year calculus-based Physics class. It is composed out of 7 modules designed to bring together experiments and computer simulations (numerical simulations). Each module uses both simulations and experiments to address a phenomena under study, and lasts for 3 weeks (21 weeks total for the whole set). Wolfram Mathematica is used for computer simulations. Topics are: Motion with a drag, Pendulum with an arbitrary amplitude, Magnus effect, Centripetal force acting on a pendulum, A ball rolling from paths of various shapes, Doppler shift and Fourier transformation, Equipotential lines. These laboratory classes were taught for 2 years at Princeton International School of Mathematics and Science. 2. Introduction In the recent decades computer simulations (numerical simulations) became an important way to do research in many areas. Complexity of many real system does not allow to obtain a model written through analytic equations. In such situations numerical simulations provide a good (if not the only one) way to do address them. Computer simulations allow people to perform weather forecasting[1], understand dynamics of star clusters[2], enhance particle physics research [3], allow to perform simulations of chemical interactions [4], as well as address many other systems. The question about educating people in the area of simulations was raised as early as in 1965 [5]. However up to now simulation did not become a part of a standard curriculum neither for the 1st or 2nd year of Physics curriculum of a typical US Physics program [6, 7, 8]. Up to know educators are working on finding a good way to include simulations in educational programs, as well as on evaluation of an efficiency of such an inclusion [9, 10, 11, 12, 13]. address: ssamsonau@gradcenter.cuny.edu (Sergey V. Samsonau)

2 A typical United States sequence of teaching Physics includes a 1-year algebrabased Physics class followed by a 1-year calculus-based class. And while it is definitely early to introduce this topic in algebra-based class, calculus-based class is well suited to do exactly that. In fact, many students may be making their decision about continuing their education in Physics or not at that level, and thus it is important to provide them with a full picture of what science is like, students shall learn that there are at least 3 ways one can do science: theory, experiment, and computer simulations. Many universities offer classes on computer simulations. However usually these are theoretical classes and don t use any experiments. We believe that introducing numerical simulations in a combination with experiments allows students to see and feel a direct connection between those, as well as understand better limitations and strengths of simulations. In such settings a design of simulations and experiments go hand in hand, and should be developed to be consistent with each other. Moreover, to the best of our knowledge there is no published set of laboratory classes available, designed to introduce students into the topic of use of computer simulations (with real experiments at the same time) during a standard calculus-based Physics class. This set of labs is designed in a paradigm of open inquiry labs, where there is no one definite answer, and students essentially perform a small research projects. Many decisions have to be made to reach a final result, such as setting up a model, designing experiment, choosing/making equipment for an experiment, determining a way to evaluate errors, and also importantly the way to compare results of simulation to results of real experiments. The results students obtain vary significantly depending on the decisions students make. The benefits of open-inquiry(and research-oriented) teaching approach are described elsewhere [14, 15, 16]. These labs were developed with a limited time for lab work in mind, giving students enough time to read about new concepts and implement code and experiments. This set was taught for two years in AP Physics Class in the following settings: 40 minutes once a week of class time (plus time allocated by students for work outside of classroom), 3 weeks per module. It is important to note, that most of the students who were doing this labs took Honors Physics class with an extensive open-inquiry laboratory classes, and were taught to handle open ended tasks, design experiment, analyze errors, and take responsibility for their own results. Completed work is submitted by students in a lab article format (instead of lab report ), which is build to be as close as possible to a typical research paper. The details about lab article format will be presented in discussed elsewhere. Wolfram Mathematica was chosen as a language for computer because it provides an easy to understand, follow (and modify) way to perform fairly good numerical simulations without a need of teaching students a full scale programming class. However, it is up to an instructor to give some references to those students who would like to learn more about deeper details and methods involved. The code provided bellow can be made smaller (in some cases), but it was delivered in this way for a better readability by students. While Mathemat- 2

3 ica has a great and extensive help with many interactive examples, some teachers and students may find useful to use books introducing Wolfram Mathematica [17, 18]. Also, taking into account how many students spend a lot of time playing computer games, a reference on use of computer simulations for game development can be used as an additional motivation [19, 20]. 3. Description of Laboratories In many of the modules students are responsible for designing their own experiments and are allowed to use all available equipment in the physics laboratory. As these assignments have a research component, students are in fact encouraged to develop their own devices in addition to equipment available in physics laboratory. In case there is a strong engineering program, they can benefit from resources available there as well. Of course, students can ask for an assistance from an instructor, but the instructor should not do the work for students and give them enough freedom and challenge to explore and learn. Thus in summary students do: Develop their own experiments to address a specified phenomena Develop (or adjust) computer simulations Compare the results of experiments with the ones produced by computer simulations Redesign experiments or simulations if needed Files with a code presented bellow can be downloaded from this link 1 in the Wolfram Mathematica notebook format Motion with a drag Task with potential subtasks Study motion with a drag Model free fall motion with a drag Determine which drag type to use (turbulent vs laminar flow) Make corresponding plots Design and perform experiment, to verify model. Experiment can contain modifications to address separate issues under study 1 3

4 Distribution of work by days Day one. An instructor makes an intro to Mathematica: basic operations, plotting, solving equations, etc. Ask students to describe free fall with kinematic equation. Plot it using an analytical equation. After that an instructor gives a basic example of use of NDSolve with free fall in no air resistance situation. Compare this solution to kinematic equation. Then compare to experiment (drop a heavy ball). Ask students to learn (at home) how to use NDSolve to solve differential equations. They shall find and read general theory related to Euler and Runge- Kutta methods. Homework: At home, students will make Mathematica program to integrate an equation of motion with a drag. They should (read and) decide which drag to use: laminar or turbulent. Students also design an experiment in which they can measure effect of a drag. For example, one can use water in high and wide graduated cylinder (we use 2000mL Polymethylpentene cylinder from FLINN scientific). Then we can either drop a small object to sink there (dense objects will do this too fast), or produce bubbles to observe how they rise (bubbles can be produced by putting a foam under water) Day two. Verify that students have a working code for situation with drag. If they do, they can do an experiment. Homework: analyze experimental data. Think whether the designed and performed experiment allows to observe effects of drag. If not, think about another design. Also, enter into Mathematica program parameters of a real experiment they are using and compare results. Day three. After they compare results between numerical simulation and real experiment at home they may need to redesign experiment and redo it Code: starter code Prepare path to save files in current directory SetDirectory [ NotebookDirectory [ ] ] ; Use a kinematic equation g = 9.8; ykin[ t, y0, v0 ]:= y0 + v0 t + g tˆ2/2; Plot [ykin[t, 0, 0], {t, 0, 2}, AxesLabel > { t, s, y, m }]; Use a differential equation s = NDSolve [{y [ t]==g, y[0]==0, y [0]==0}, y, {t, 0, 5}]; Plot [ Evaluate [y[t ]/. s], {t, 0, 3}, PlotRange > All, AxesLabel > { t, s, y, m }] Show together p1 = Plot [ykin[t, 0, 0], {t, 0, 3}, PlotStyle > {Black, Dotted}, AxesLabel > { t, s, y, m }, ImageSize > 300, PlotLabel > Kinematic ]; p2 = Plot [ Evaluate [y[t ]/. s ],{t,0,3}, PlotStyle > {Gray, Dashed}, AxesLabel > { t, s, y, m }, ImageSize > 300, 4

5 PlotLabel > NDSolve ]; p3 = Show[p1, p2, PlotLabel > Two plots shown together ]; p4 = Deploy@Grid[{{p1, p2, p3}}, Dividers > Gray, Spacings > {2, 2}] Export[ p4. eps, p4]; y,m Kinematic y,m NDSolve y,m Two plots shown together t,s t,s t,s Code: add drag Students have to read about laminar flow, turbulent flow, and Reynolds number k = 1; ( Students have to think how to choose this constant ) rho = 1; ( Students have to read and choose an appropriate value for air or water ) ; Cd = 1; ( Students have to read and choose this drag constant ) A = Pi (0.01) 2; ( cross section area ) m = 0.01; LaminarDrag[v ]:= Sign [v] k Abs[v] A TurbulentDrag[v ]:= Sign [v] 1/2 rho vˆ2 Cd A Fg = m g; Tmax = 0.5; ( Solve ) sfree = NDSolve [{y [ t]==fg /m, y[0]==0, y [0]==0}, y, {t,0,tmax}]; sl = NDSolve [{y [ t]==(fg + LaminarDrag[y [ t ]]) /m, y[0]==0, y [0]==0}, y, {t,0,tmax}]; st = NDSolve [{y [ t]==(fg + TurbulentDrag[y [ t ]]) /m, y[0]==0, y [0]==0}, y, {t,0,tmax}]; ( import actual data ) data = Import[ sample. csv ]; ( use //Table to see it as a table ) pdata = ListPlot [data, PlotStyle > Black ]; ( Plot ) py = Plot [{ Evaluate [y[t ]/.sfree ], Evaluate [y[t ]/.sl], Evaluate [y[t ]/.st]}, {t, 0, Tmax}, 5

6 PlotRange > All, PlotLegends > { y No Drag, y Laminar, y Turbulent }, AxesLabel > { t, s, y, m }, PlotLabel > Coordinate ] pv = Plot [{ Evaluate [y [ t ]/.sfree ], Evaluate [y [ t ]/.sl], Evaluate [y [ t ]/.st]}, {t, 0, Tmax}, PlotRange > All, PlotLegends > { v No Drag, v Laminar, v Turbulent }, AxesLabel > { t, s, y, m }, PlotLabel > Velocity ] pp = Deploy@Grid[{{pY, pv}}, Dividers > Gray, Spacings > {2, 2}] Export[ pp. eps,pp]; y,m Coordinate v,m Velocity t,s t,s No Drag No Drag -0.6 Laminar Drag Turbulent Drag -2 Laminar Drag Turbulent Drag Code: error evaluation and analysis datatime = data[[ All, 1]]; ( time values in data ) datacoordinate = data[[ All, 2]]; ( coordinate values in data ) Prediction = Evaluate [y[ datatime ]/.sl] [ [ 1 ] ] ; ( corresponding prediction from simulation ) ( Standart Deviation of a difference ) StDev = Standar ddevi ati on [ Prediction datacoordinate]; ( Maximum Percentage Error ) MaxPercentageError = Max[100 ( Prediction datacoordinate) /(0.5 ( Prediction + datacoordinate)) ]; ( add error bar to plot ) Needs[ ErrorBarPlots ] ( load package ) InstError = 0.04; ( set up an instrumental error ) error = ConstantArray[InstError, Length [data ] ]; ( add instrumental error ) ( students can extend this approach to include accidental error ) witherror=transpose[{data[[ All, 1]], data[[ All, 2]], error }]; pdatanew = Show[ ListPlot [data], ErrorListPlot[withError, PlotStyle > Black ] ] ; ( Plot ) pyerror = Show[ Plot[{ Evaluate [y[t ]/.sfree ], Evaluate [y[t ]/.sl], Evaluate [y[t ]/.st]}, {t, 0, Tmax}, PlotRange > All, PlotLegends > { No Drag, Laminar Drag, 6

7 Turbulent Drag }, AxesLabel > { t, s, y, m }, ImageSize > 400, PlotLabel > Coordinate ], pdatanew]; ( prepare separate plot for errors ) erp1 = Plot [{InstError}, {x, 0, Tmax}, PlotStyle > Gray ]; erp2 = Plot[{ InstError}, {x, 0, Tmax}, PlotStyle > Gray ]; perror = ListPlot [ Transpose[{dataTime, datacoordinate Prediction}], PlotLabel > Error, AxesLabel > { t, s, Prediction Actual, m }, ImageSize > 400, PlotStyle > Black ]; perrort = Show[pError, erp1, erp2]; pperror = Deploy@Grid[{{pYError, perrort}}, Dividers > Gray, Spacings > {2, 2}] Export[ pperror. eps,pperror ]; y,m Coordinate Prediction - Actual, m Error t,s No Drag 0.02 Laminar Drag TurbulentDrag t,s Pendulum with an arbitrary amplitude Task with potential subtasks Study motion of a real pendulum. The simple formula for simple pendulum is valid in an approximation of small amplitude, but what will happen if the angle is larger? Model motion with any amplitude Study influence of mass Study influence of drag Make corresponding plots Make animation of a motion Design and perform experiment, to verify model. Experiment can contain modifications to address separate issues under study 7

8 Distribution of work by days Students are free to schedule their work on their own. An instructor will provide starter code on the first and second days. Day one. Starter code for pendulum with any amplitude. Using this code students can start developing model, think how to obtain period, etc. Optional: students can find period automatically (by finding maximum of phi(t)) and do plot of period vs amplitude. Also instructor gives a code to generalize drag force using Piecewise function of Mathematica. Some students may choose to add drag force to their equations. Day two. An instructor helps to write code to create animations (shows some examples). Also show briefly results when drag is taken into account. Don t give the code to produce this result. Day three. Students choose how to use this time Code: starter code Small amplitude SetDirectory [ NotebookDirectory [ ] ] ; L = 0.5; g = 9.81; T = 2 Pi Sqrt [L/Abs[g]] Use differential equation (any amplitude) m = 0.1; Ft[phi ]:= m Abs[g] Sin [phi] MomInertia = m L2 s = NDSolve [{phi [ t]==(ft[phi[t ]] L)/MomInertia, phi[0]==pi/2, phi [0]==0}, phi, {t, 0, 10}]; p1 = Plot [{ Evaluate [phi[t ]/. s ], Evaluate [phi [ t ]/. s ], Evaluate [phi [ t ]/. s ]}, {t, 0, 10}, PlotRange > All, PlotLegends > { phi, phi, phi }] Export[ pendulum p1. eps, p1] 8

9 p Code: animation, period Animation x[ t ]:= L Sin [phi[t ]] y[ t ]:= L Cos[phi[t ]] MaxTime =2; AnF := Animate[ Graphics [ {Red, PointSize > Large, Point[{x[t], y[t]}/. sol ], Green, Thick, Line[{{0,0}, {x[t], y[t]}}/. sol ]}, PlotRange > {{ 1.1 L, 1.1 L}, { 1.1 L, 0}}, Axes > True, Ticks > False ], {t, 0, MaxTime, MaxTime/100}, AnimationRunning > False ] AnF/. sol > s ( Animation ) Compare period at different amplitudes s2 = NDSolve [{phi [ t]==(ft[phi[t ]] L)/MomInertia, phi[0]==pi/2, phi [0]==0}, phi, {t, 0, 10}]; s6 = NDSolve [{phi [ t]==(ft[phi[t ]] L)/MomInertia, phi[0]==pi/6, phi [0]==0}, phi, {t, 0, 10}]; s10 = NDSolve[{phi [ t]==(ft[phi[t ]] L)/MomInertia, phi[0]==pi/10, phi [0]==0}, phi, {t, 0, 10}]; p2 = Plot [{ Evaluate [phi[t ]/.s2], Evaluate [phi[t ]/. s6], Evaluate [phi[t ]/.s10]}, {t, 0, 10}, PlotRange > All, PlotLegends > { Pi/2, Pi/6, Pi/10 }] Export[ pendulum p2. eps,p2] 9

10 Pi/2 Pi/6 Pi/ Animate AnF/. sol > s10 ( Animation ) 3.3. Magnus effect Task with potential subtasks Study motion affected by Magnus effect. In this case the motion should be described using differential equations for two coordinates (x and y). Model motion in 2 dimensions which is not affected by the Magnus effect. Students have to remember that Drag force direction will change, according to the change in velocity direction. Model motion affected by Magnus effect. Students have to remember that direction of both Drag and Magnus forces will change, according to the change in velocity direction. Study influence of object s radius (if cylindrical object is used) Study influence of rotation speed Make corresponding plots Design and perform experiment, to verify model. Experiment can contain modifications to address separate issues under study Distribution of work by days Students are free to schedule their work on their own. Starter code can be shown as an example on a screen, but not given to students (allowing students to write their own implementation). Assistance can be given to those who needs help. Note: Students will design their own experiments. For initial demonstration of an effect an instructor can use a paper cylinder (from a sheet of paper) rolling 10

11 from an incline. Different radii cylinders can be used. Alternatively, two paper cups connected on the bottom with a tape, can be launched with a rubber band, as can be found on some of YouTube videos Code: Magnus effect 2D Motion with drag SetDirectory [ NotebookDirectory [ ] ] ; rho = 1.2; ( kg/mˆ3, air density ) R = ; ( radius ) Cd = 1.17; ( drag coefficient ) TurbulentDrag[v ]:= (1/2) rho vˆ2 Cd A ; ( if turbulent drag is used ) A = Pi Rˆ2; ( crossection area ) m = 0.01 ; ( 10g ) g = 9.8; V[vx, vy ]:= Sqrt[vx2 + vy2]; CosAlpha[vx, vy ]:= vx/v[vx, vy]; SinAlpha[vx, vy ]:= vy/v[vx, vy]; Tmax = 1.5; sfree = NDSolve [{ m x [ t]==0, m y [ t]==m g, x[0] == 0, x [0]==5, y[0]==0, y [0]==0}, {x, y}, {t, 0, Tmax}]; st = NDSolve [{ m x [ t]==turbulentdrag[v[x [ t], y [ t ]]] CosAlpha[x [ t ], y [ t ]], m y [ t]==m g + TurbulentDrag[V[x [ t], y [ t ]]] SinAlpha[x [ t ], y [ t ]], x[0]==0, x [0]==5, y[0]==0, y [0]==0}, {x, y}, {t, 0, Tmax}]; pf = ParametricPlot [ Evaluate [{x[t], y[t]}/.sfree ], {t, 0, Tmax}, AxesLabel > { x coordinate, m, y coordinate, m }, PlotLabel > Trajectory of motion in 2D with drag, AspectRatio > Automatic, PlotStyle > Blue ]; pd = ParametricPlot [ Evaluate [{x[t], y[t]}/.st], {t, 0, Tmax}, AxesLabel > { x coordinate, m, y coordinate, m }, PlotLabel > Trajectory of motion in 2D, AspectRatio > Automatic, PlotStyle > {Dashed, Red}]; Show[pF, pd]; Add Magnus effect (using cylinder as an object) L = 0.2; ( cylinder length ) r = 0.02; ( cylinder radius ) 11

12 omega=10; ( assume it is a constant during motion, for a first approximation ) ( Students have carefully consider sign of omega ) G = 2 Pi r 2 omega; Magnus[v ]:= L rho v G sm = NDSolve [{ m x [ t]==turbulentdrag[v[x [ t], y [ t ]]] CosAlpha[x [ t ], y [ t ]] + Magnus[V[x [ t ], y [ t ] ]] CosAlpha[x [ t ], y [ t ]], m y [ t]==m g + TurbulentDrag[V[x [ t], y [ t ]]] SinAlpha[x [ t ], y [ t ]] Magnus[V[x [ t ], y [ t ] ]] SinAlpha[x [ t ], y [ t ]], x[0] == 0, x [0]==5, y[0]==0, y [0]==0}, {x, y}, {t, 0, Tmax}]; pm = ParametricPlot [ Evaluate [{x[t], y[t]}/.sm], {t, 0, Tmax}, AxesLabel > { x coordinate, m, y coordinate, m }, PlotLabel > Trajectory of motion in 2D, and Magnus effect, AspectRatio > Automatic, PlotStyle > Dotted]; pt = Legended[Show[pF, pd, pm], LineLegend[{Blue, {Dashed, Red}, {Dotted, Green}}, { Free, Drag, Magnus }]] Export[ Mangus plot. eps, pt]; 12

13 y coordinate, m T omotionin2dwithdrag x coordinate, m Free Drag M Centripetal force and Pendulum Task with potential subtasks Study a force acting on a string suspension point produced due to a pendulum motion. In this lab students will compare results from an experiment, a computer simulation, and from an analytical expression. Make a numerical model. Study the influence of various factors. Obtain an analytical expression for force acting on a force sensor (first, assuming no air drag, and then with air drag). Make corresponding plots. Design and perform experiment, to verify model. Experiment can contain modifications to address separate issues under study. Compare results produced by three approaches: analytical, numerical, experimental. 13

14 Distribution of work by days Students are free to schedule their work on their own. No starter code is given this time. Note: For experimental part students can use any force sensor they wish. In our lab we use a force sensor from Vernier company, with a data interface to Vernier Logger Pro software. Students can import the collected data into Mathematica using Import command and plot it on the same graph with points produced by simulation. Routing of a thread holding the bob of a pendulum fromtheforcesensorshouldbedoneinsuchaway,thatthereisonlylongitudinal pull on a force sensor Code No new code needed for this lab A ball rolling from paths of various shapes Task with potential subtasks Given tracks of irregular shapes, students will create a numerical model of a motion of a ball rolling along that tracks. This introduces an idea of use of numerical simulation to describe motion of an object given conditions we can t describe analytically. For example, motion of a space ship moving near the planets, large asteroids, and other objects. Make image of tracks, and obtain coordinates describing their shapes. Make a numerical model of a ball rolling down the tracks. Make corresponding plots. Compare results: numerical simulation, experimental Distribution of work by days Students are free to schedule their work on their own. An instructor will provide starter code on the first day. Note: We use PSworks Conservation of Energy Track from Flinn Scientific, by anyother setup with severaldifferent paths can be used. Ball rollsdown each path with different dependence of coordinate and velocity form time. Students can use a video camera in a smartphone to obtain experimental dependence of those values and then compare them with prediction of a model. In a simplified approach, only a time of the whole motion can be used to compare prediction of model to observation. In order to make a numerical model, students will have to extract x-y coordinates combination describing the shape of each track. This can be done by making photo of a track, importing it into Mathematica, and then using Coordinates tool available under the image to extract coordinates. Subsequent interpolation is required to obtain a function, derivative of which will be calculated. Moving average smoothing can be used before interpolation. Alternatively, video processing tools of Mathematica can be used to extract the information about trajectory and a dependence of coordinates from time automatically (the code for that part is provided below). 14

15 Code: starter code Let say you get a list of point from the Coordinates tool ( generating a list of point to be used as an example ) a = Table[{ i, iˆ2 3 i Random[]}, {i, 1, 3, 0.2}]; ListPlot [a, AspectRatio > Automatic, PlotStyle > {PointSize [Medium], Red}] Do an interpolation b = Interpolation [a]; ( b now is a function of x ) Show[ Plot [b[x], {x, 1, 3}], ListPlot [a, PlotStyle > {PointSize [Medium], Red}], AspectRatio > Automatic ] Find an angle of the curve Der = b ; ( derivative, or tan(alpha) ) Plot[{b[x], Der[x]}, {x, 1, 3}, AspectRatio > Automatic ] angle [x ]:= ArcTan [Der[x]]/ Degree ; Plot [angle [x], {x, 1, 3}] D&'()*+,-. 6 Data I!"#$% Angle / Now students write and solve system of equations. Then they can use ND- Solve 15

16 Code: automatic Image processing Use a static, preferably monotone background to produce a video (students have to use slow motion option of their smartphones to have sufficient time resolution). Import the video. S e t D i r e c t o r y [ N o te b o o k D i r e c t o r y [ ] ] ; p = Import [ v i d e o.mov, I m a g e L i s t ] ; ( p = Import [ v i d e o. mp4, { AVI, I m a g e L i s t } ] ; ) Length [ p ] Display selected frames Show [ p [ [ 1 ] ], Im a g e S i z e > ] ; Show [ p [ [ 1 0 ] ], Im a g e S i z e > ] ; Show [ p [ [ 2 7 ] ], Im a g e S i z e > ] ; Deploy@Grid[ { { p1, p2, p3} }, Dividers > Gray, Spacings > { 2, 2 } ] S e t D i r e c t o r y [ N o te b o o k D i r e c t o r y [ ] ] ; Export [ r o l l i n g d o w n f r a m e s. e p s, pp ] ; p1 p2 p3 pp = = = = Frame 1000 Frame 1100 Frame 4679 Extract positions BallPositions = ImageFeatureTrack [ p [ [ 1 ( t a k e s a l o n g time to compute ) ;; 2 7 ] ], MaxFeatures > 1 ] ; ppos = L i s t P l o t [ BallPositions, Im a g e S i z e > 4 0 0, P l o t L a b e l > P o s i t i o n s ] Remove duplicates in x and do interpolation, apply smoothing, interpolate, plot ( Remove r e c o r d s which have th e same x ) ( add frame number f i r s t, to keep i n f o r m a t i o n about time ) PositionsWithTime = MapThread [ F l a t t e n /@ Thread [{##}] &, { BallPositions, Range [ 1, Length [ BallPositions ] ] } ] [ [ A l l, 1 ] ] ; Lenght b e f o r e removal o f d u p l i c a t e s PositionsWithTime ] ] i s : <> T o S tr i n g [ Length [ ( remove d u p l i c a t e s ) PositionsWithTimeNoDublicates = Union [ PositionsWithTime, SameTest > ( # 1 [ [ 1 ] ] == # 2 [ [ 1 ] ] &) ] ; Lenght a f t e r removal o f d u p l i c a t e s i s : <> T o S tr i n g [ Length [ PositionsWithTimeNoDublicates ] ] ( s e p a r a t e back i n t o time and c o o r d i n a t e s ) time = PositionsWithTimeNoDublicates [ [ A l l, 3 ] ] ; coord = PositionsWithTimeNoDublicates [ [ A l l, { 1, 2 } ] ] ; 16

17 pcoord = ListPlot [coord, ImageSize > 400, PlotLabel > Positions, no duplicates ] Do Interpolation ( if smooting needed smothed = MovingAverage[ coord, 2] ; ListPlot [smothed] ff= Interpolation [smothed ]; ) ff6 = Interpolation [coord]; Minx = Min[coord[[ All, 1 ]]]; Maxx = Max[coord[[ All, 1 ]]]; pfit = Plot [ ff6 [x], {x, Minx, Maxx}, ImageSize > 400, PlotLabel > Fit ] Get derivative ( Derivative ) pder = Plot [ ff6 [x], {x, Minx, Maxx}, ImageSize > 400, PlotLabel > Derivative ] Deploy@Grid[{{pPos, pfit, pder}}, Dividers > Gray, Spacings > {2, 2}] P:;<=>?@A RSU Fit bcd efghijklmn NOQ _`a JKL BCE FGH \]^ VWX YZ[ Please note, an alternative code for image processing is available for download at this link Doppler shift and Fourier transformation Task with potential subtasks Study Doppler effect and use it to produce a sound of an object moving with a given velocity. Along the way students will get familiar with Fourier transformation and idea of alternative representations of signals: using time domain vs using frequency domain. Make a program producing a sound form a moving object given a sound from a standing object. Record sound of a signal of a standing car

18 Produce sound of moving car using program Record sound of a signal of a moving car Compare sound generated numerically with real sound. It is an interesting taskforstudentstothinkwhatisthegoodwaytocomparepredictedsignal and recorded signal, as they have to produce an objective way to compare those, including error estimation, etc. Make corresponding plots Distribution of work by days Day one. Concept of Doppler shift and Fourier transformation will be explained, starter code will be given. Some examples with tuning fork or other instruments can be given. Student will work at home on understanding the underlying theory and adjusting program if needed. Day two. An instructor will produce a sound using a signal with his/her car, and students will record it with their smartphones. Then they will have to produce the shifted sound. Day three. An instructor will drive the same car at miles/hour near students, who will be recording the sound Code: Fourier Transformation, starter code ClearAll [ Global ] record1 = Sound [SoundNote[ G, 1, Violin ]] Periodogram[ record1] 1s 18

19 Shifting sound by a fixed frequency res=audiofrequencyshift [ record1, Quantity[1000, Hertz ]] Periodogram[{record1, res}, PlotRange > {{0, 4000}, All}, ScalingFunctions > Absolute ] Shifting sound using a given coefficient new = AudioPitchShift[record1, 1.5] 19

20 Periodogram[{new, record1}, PlotRange >{{0,4000}, All}, ScalingFunctions > Absolute ] Do audio recording and simply drag the file into Mahtematica. Record2 = (drag and drop audio file ) At this stage students have enough knowledge on how to work with sound to complete the task. Some students may also benefit from functions helping to extract data in a form of a list. Code for such functions is presented bellow. ( get Amplitude data in time domain ) AudioData[new] AudioData[new] [ [ 1 ] ] ; AudioData[new] [ [ 2 ] ] ; ListPlot [AudioData[new ] [ [ 1 ] ] ] ; ListPlot [AudioData[new ] [ [ 2 ] ] ] ; ( squared magnitude of the discrete Fourier transform (power spectrum) of list. ) PerData = PeriodogramArray[ new] ListLinePlot [PerData, PlotRange > {{0, 4000}, All }] ( apply filter if needed ) newfiltered = BandpassFilter[ new, {400, 1000}] PerDataFiltered = PeriodogramArray[ newfiltered]; ListLinePlot [PerDataFiltered, PlotRange > {{0, 4000}, All }] 20

21 3.7. Equipotential lines lab Task with potential subtasks Use numerical simulations to produce an image of equipotential lines, given the shape of electrodes. Compare to experiment. Using Mathematica calculate position of equipotential lines, given interesting (to explore) combination of electrodes in a plane. Note: a relatively simple way to construct irregularly shaped electrodes with Mathematica, is to combine ellipses (see an example in the code) Print an image produced by Mathematica using printer On a conductive paper make the same electrodes arrangement using conductive ink, or conductive tape. On a printed paper make small holes for voltmeter probes along the equipotential lines. Align this paper with the conductive paper. Apply voltage to electrodes on conductive paper. Using a voltmeter measure the voltage in the premade holes. Evaluate how precise numerical predictions are Distribution of work by days Day one: starter code and first calculations. Day two: Preparing conductive paper with electrodes. Day three: Final measurements Code: starter code boundaries={ y 0, ( boundaries [[1]] ) y 1,( boundaries [[2]] ) 0 x, ( boundaries [[3]] ) 2.5+x, ( boundaries [[4]] ) 1/200 ( 0.2+x)2 (y 0.5)2, ( boundaries [[5]] ) 1/200 ( 1+x)2 (y 0.5)2,( boundaries [[6]] ) 1/5 ((x 2)Cos[ Pi/6]+(y 0.38)Sin [ Pi/6])2/0.252 ((x 2)Sin [ Pi/6]+(y 0.38)Cos[Pi/6]) 2/0.052,( boundaries [[7]] ) 1/5 ((x 2)Cos[ Pi/6]+(y 0.62)Sin[ Pi/6]) 2/0.252 ((x 2)Sin[ Pi/6]+(y 0.62)Cos[ Pi/6]) 2/0.052,( boundaries [[8]] ) 1/1.25 ((x 2.2)Cos[ Pi/2]+(y 0.5)Sin [ Pi/2]) 2/0.252 ((x 2.2)Sin [ Pi/2]+(y 0.5)Cos[ Pi/2]) 2/0.032 ( boundaries [[9]] ) }; DefinedRegion = ImplicitRegion [And@@(#<=0&/@boundaries), {x, y}]; 21

22 Show[ RegionPlot[ DefinedRegion], ContourPlot[ Evaluate [ Thread [ boundaries ==0]], {x, 0.1, 5}, {y, 0.1, 2}, ContourStyle > {Purple, Purple, Purple, Purple, Blue, Red, Green, Green, Green }], AspectRatio > Automatic ] Specify voltage on boundaries and solve ( If condition for a boundary not given, potential can have any value there ) sol=ndsolvevalue[{laplacian[u[x, y], {x, y}]==0, {DirichletCondition[u[x, y]==5, boundaries[[5]]==0], DirichletCondition[u[x, y]==0, boundaries[[6]]==0], DirichletCondition[u[x, y]== 2, boundaries[[7]]==0], DirichletCondition[u[x, y]== 2, boundaries[[8]]==0], DirichletCondition[u[x, y]== 2, boundaries[[9]]==0]}}, u, {x, y}, Element[{x, y}, DefinedRegion ] ]; ContourPlot[ sol [x, y], Element[{x, y}, DefinedRegion], Mesh > None, ColorFunction > Electric Potential Map, Contours > 18, AspectRatio > Automatic ]

23 4. Conclusion Abovewepresentedthe complete setof7laboratorymodulesreadytobe implemented in advanced high schools or college settings. Of course, this sequence is not set in stone an adjustment can/should be made by an instructor to adjust for a particular group of students. This is our hope, that this paper will serve as a basis for development of many other laboratory courses combining experiment with computer simulations. It would be interesting to perform a study involving many students to obtain statistical data on pedagogical outcomes. 5. Acknowledgements We express gratitude to students of Princeton International School of Mathematics and Science for their effort, open minded approach, and interesting ideas about implementation of some of the topics. We appreciate the support of administration and faculty of Princeton International School of Mathematics and Science. References [1] Martin J. Miller and Piotr K. Smolarkiewicz. Predicting weather, climate and extreme events. Journal of Computational Physics, 227(7): , mar [2] Sverre J. Aarseth. Computer simulations of star cluster dynamics. Vistas in Astronomy, 15(C):13 37, jan [3] F Karsch and E Laermann. Numerical simulations in particle physics. Reports on Progress in Physics, 56(11): , nov [4] Jacob R. Gissinger, Benjamin D. Jensen, and Kristopher E. Wise. Modeling chemical reactions in classical molecular dynamics simulations. Polymer, 128: , oct [5] Robert L. Egbert. The Role of Computer Simulation in Education, [6] David Halliday. Fundamentals of physics. Wiley, [7] Douglas C. Giancoli. Physics: Principles with applications. Pearson, [8] CollegeBoard. The College Board AP Physics C Descriptions, [9] Nico Rutten, Wouter R. Van Joolingen, and Jan T. Van Der Veen. The learning effects of computer simulations in science education. Computers and Education, 58(1): , [10] Bat-Sheva Eylon and Et al. Computer Simulations as Tools for Teaching and Learning: Using a Simulation Environment in Optics. Journal of Science Education and Technology, 5(2):93 110, jun

24 [11] Elizabeth Gates Bradley and Brittany Kendall. A Review of Computer Simulations in Teacher Education. Journal of Educational Technology Systems, 43(1):3 12, sep [12] Cândida Sarabando, José P. Cravino, and Armando A. Soares. Contribution of a Computer Simulation to Students Learning of the Physics Concepts of Weight and Mass. Procedia Technology, 13: , jan [13] Ersin Bozkurt and Aslan Ilik. The effect of computer simulations over students beliefs on physics and physics success. Procedia - Social and Behavioral Sciences, 2(2): , jan [14] Desmond H. Murray, Sherine Obare, and James Hageman. Early Research: A strategy for inclusion and student success. In ACS Symposium Series, volume 1231, pages 1 32, jan [15] Erik R. Mohlhenrich, Sergey V. Samsonau, and Roxanne P. Spencer. Integrating science through authentic research in secondary schools. In 2018 IEEE Integrated STEM Education Conference (ISEC), pages IEEE, mar [16] Clark A Chinn and Betina A Malhotra. Epistemologically Authentic Inquiry in Schools: A Theoretical Framework for Evaluating Inquiry Tasks. Inc. Sci Ed, 86: , [17] Stephen Wolfram. An Elementary Introduction to the Wolfram Language. Wolfram Media, Incorporated, [18] Cliff Hastings, Kelvin Mischo, and Michael Morrison. Hands-on start to Wolfram Mathematica: and programming with the Wolfram language. Wolfram Media, [19] David M. Bourg, Kenneth. Humphreys, and Bryan. Bywalec. Physics for game developers. O Reilly Media, [20] Grant. Palmer. Physics for game programmers. Apress,