Lab Manual 101L S24 (4)

pdf

School

University of Texas *

*We aren’t endorsed by this school

Course

101L

Subject

Physics

Date

Apr 3, 2024

Type

pdf

Pages

64

Uploaded by PresidentNeutronDog35

Report
1 PHYSICS 101L LABORATORY MANUAL K.W. Gentle J.A. Yeazell Department of Physics University of Texas at Austin January 2024 Copyright 2023, The University of Texas at Austin Note : The laboratory sessions do not begin immediately, you should consult the syllabus on Canvas for the schedule. You must prepare for each Lab. For Lab 0, you must read the Lab before joining your session of the class. You should always arrive prepared to begin work without further instruction. The course has no lectures; reading is required. If you do not arrive prepared, you will be unable to complete the Lab in the time available.
2 Getting Started The focus of this course is on the measurement and modeling of the motion of objects (Mechanics). You will design experiments and develop computational models that explore the key concepts in this area, e.g. conservation of momentum. The laboratory activities will start with a warmup exercise and then you will extend this to investigate a more complex phenomenon. This course has a few goals: 1) Foster a physics/science community. Find fellow physics/science majors to work with. Science and physics is often a product of teamwork. The development of colleagues and hopefully friends can make it easier, more interesting (and more fun) to learn and do physics. 2) Introduce you to attacking physics problems experimentally and numerically. It asks you to design experiments and introduces you to the computer modeling of physics problems. 3) Help you to communicate scientific results in a written format and introduce you to oral presentations of scientific topics. 4) Introduce you to finding scientific information. How to search databases for peer-reviewed articles on a topic. You are expected to read through the scheduled activity before the start of the lab session so that you are familiar with the topics and start thinking about the design of an experiment or of a computational model. With this preparation, you should be able to complete all the work for this course within the lab session (including any required scientific writing describing your experiment and results). The laboratory has Apple computers that will be used for data acquisition and can be used for computer modeling (you are welcome to use your own laptop for modeling). These Apple computers can be accessed in the usual way with your UTEID and associated password. Please save any data files you generate to UTBox or your favorite Cloud system (not to the local drives on the computers in the lab).
3 Table of Contents Lab 0. Measurement, Uncertainties, Accuracy, and Using Excel .............................................. 4 Lab 1. Scientific Literacy and Measurement Apparatus .......................................................... 16 Lab 2. Experiment on motion in 1D .......................................................................................... 23 Lab 3. Scientific Presentation .................................................................................................... 29 Lab 4. Modeling Motion ............................................................................................................. 33 Lab 5. Modeling Orbital Motion ................................................................................................. 39 Lab 6. Experiment on work, energy, and momentum ............................................................. 43 Lab 7. Computational Model of Collisions ................................................................................ 47 Lab 8. Experiment on Rotational Dynamics and Conservation of Angular Momentum ....... 52 Lab 9. Experiment on harmonic motion ................................................................................... 53 Lab 10. Computational Model of Waves ................................................................................... 58
Your preview ends here
Eager to read complete document? Join bartleby learn and gain access to the full version
  • Access to all documents
  • Unlimited textbook solutions
  • 24/7 expert homework help
4 Lab 0. Measurement, Uncertainties, Accuracy, and Using Excel This activity addresses a range of topics essential to experimental science. It covers: What is a measurement? What is the precision of a number or measurement? What is the uncertainty of a measurement? Error propagation: If you need to make more than one measurement to determine a more complicated quantity, how do you combine the uncertainties of each measurement to find the uncertainty of this complicated quantity? What is accuracy and how can you judge experimental results? How to use Excel to analyze data and to display experimental results graphically (including error bars). Accuracy and Uncertainty in Experimental Measurements Measurement is a central aspect of all science. Our ability to model the phenomena we observe is tied to the accuracy and confidence we have in those measurements. Here we describe how to determine the limits of experimental measurements. This topic is commonly known as “error analysis”. It is important to note that real errors or mistakes in a measurement are not included in this analysis. Real errors (also known as systematic errors) are things like using the measurement device improperly, forgetting to zero the device, the device is out of calibration, etc. If known these real errors should be eliminated, and the measurement repeated. The errors studied here typically stem from the basic accuracy of the measurement device and the random fluctuations in measurement due to noise. Precision and accuracy have two different meanings in error analysis. Precision is tied to number of digits to which a number is known or the significant digits in that number. Accuracy is how well a measurement compares with an accepted value. For example, you use a digital scale to measure a 1 kg standard mass. Your result is 1.0000252 kg. The difference from 1 kg is the accuracy (0.0000252). Repeated measurements reveal that the last digit fluctuates and this fluctuation leads to the precision of the measurement. The measurement has 3 significant digits and this is the precision. Precision is best defined by examples: for a measurement in decimal notation, the precision is the number of significant digits, e.g. 7.235 has 4 significant digits; 𝜋 or √2 has infinite precision, while the approximation 3.14 has only 3 significant digits; an integer 1, 2, … has infinite precision, but 2.0 only has 2 significant digits; in scientific notation the number 3.287 × 10 −5 has 4 significant digits. Note zeros to the right of the decimal point may act only as place holders and have no significance, e.g. 0.03 has only one significant digit while 0.0308 has 3 significant digits. The number of significant digits of a quantity that results from a multiplication or division is determined by the number with the least significant digits, e.g. ( 5.00 7.145 ) ( √3 2 )
5 The precisions of √3 and the integer 2 are infinite so the result of the calculation is determined the 3 significant digits in 5.00 . While doing the calculations you should keep more digits than are significant in the intermediate steps and only round to the correct number of significant digits at the end. The correct way to report the result of the above calculation is 0.606 . There are different rules for the significant figures resulting from addition and subtraction. The number with the least number of decimal places to the right of the decimal is the determining factor. If one has 5 places to the right of the decimal point while the other only has 3, then the sum or the difference is rounded to the 3 rd decimal place. For example, correctly rounded the following sum of 1.51475 + 0.322 = 1.836. Precision is closely tied to topic of uncertainty. The uncertainty in a measured value describes how well you know that value. Let’s say you me asure the length of a table with a ruler or meter stick and you find it be 2.154 m. Why did you only report 4 significant digits? On most rulers or meter sticks, the smallest division corresponds to a millimeter. The uncertainty in this measurement is typically defined as ±0.5 of this smallest division. All experimental measurements are reported as two numbers: the first number is the best estimate of the quantity being measured (usually the average of several measurements) and the second gives the bounds of the uncertainty. Exercise 1 : You and your partner should measure the length of a table twice with a meter stick and calculate the average. Then report your measurement below, ________________ ± ________________ The reported length should have a form like 2.154m ±0.0005m. Note that the uncertainty has the same units as the best estimate. This is known as absolute uncertainty. The uncertainty can also be given in terms of the ratio of uncertainty to the best estimate. This is known as the relative uncertainty. In the above example, the measurement of the table could also be written, 2.154m ± [0.0005/2.154] = 2.154m ± [0.0002] or 2.154m ± 0.02% The square brackets indicate it is a relative uncertainty.
6 Physical constants are reported with a specified precision or number of significant digits. There is an implied uncertainty tied to this precision. The implied minimum uncertainty is ±5 in the next decimal place. For example, a reported value of the universal Gravitational constant of ? = 6.67 × 10 −11 N ∙ m 2 /kg 2 has an uncertainty of ±0.005 × 10 −11 N ∙ m 2 /kg 2 . A measurement of a length as 2.154 ±0.0005 m does not mean that the actual length lies within ±0.0005m of the 2.154 m. Often these uncertainties are created on the assumption that ~2/3 of the time the actual value will lie in this range. If the table was measured with a more accurate and precise tool, you should not be too surprised if the result was outside the range. For example, interferometry can measure distances based on an accurately known wavelength of light and has a precision of smaller than the wavelength of the light. The measurement is a best estimate of the value tied to the measurement methods used. Error propagation The uncertainty of a quantity that is based on several measurements can be found if we know the uncertainties in each of the individual measurements of quantities needed to calculate this composite quantity, e.g. measurements of the length, width, and height of a block are needed to determine the volume of the block. Finding the uncertainty of a composite quantity is known as error propagation ”. Conventional error propagation analysis is a direct application of calculus using the Taylor Series expansion. That is, that the value of a function at small distance ? from a point ? 0 is, ?(? 0 + ?) = ?(? 0 ) + ? ?? ?? | ?=? 0 + ⋯ We are showing only the first two terms of the series since that is all we will need. The first term will be the best estimate of the measure and the second will be the uncertainty. Let’s assume we want to measure a composite quantity ? that is a function of a number of measurements of other quantities ? 𝑖 so that ? = ?(? 𝑖 ) . The uncertainty, ? , in ? is given by the second term of the Taylor series summed over all the measurements, ? 2 = ∑(∆ 𝑖 ) 2 ( ?? ?? 𝑖 ) 2 𝑖 where partial derivatives are used for each ? 𝑖 and the 𝑖 are the corresponding uncertainties for different measurments. The above equation for the uncertainty leads to the following convenient rules for propagating error:
Your preview ends here
Eager to read complete document? Join bartleby learn and gain access to the full version
  • Access to all documents
  • Unlimited textbook solutions
  • 24/7 expert homework help
7 1) Sums and Differences : If the desired quantity, ? , is the found by summing two (or more) measurements ? and ? (e.g. total length = sum of the two different lengths) then the absolute uncertainty ∆? in terms of the absolute uncertainties ∆? and ∆? is given by, (∆?) 2 = (∆?) 2 + (∆?) 2 + ⋯ (Note: the quantities MUST have the same units.) The uncertainty for quantities derived from differences also add using the exact same rule. Examples: (6.1cm ± 0.2cm) + (4.2cm ± 0.2cm) = 10.3cm ± √(0.2) 2 + (0.2) 2 cm = 10.3cm ± 0.3 cm (6.1cm ± 0.2cm) - (4.2cm ± 0.2cm) = 1.9cm ± √(0.2) 2 + (0.2) 2 cm = 1.9cm ± 0.3 cm The subtraction of quantities leads to an increased relative uncertainty in the result. In the second example, the relative uncertainty of the result is larger than either of quantities being subtracted. In the extreme case, where two nearly identical quantities are being subtracted, essentially all certainty can be lost (the relative uncertainty nears infinity)! Exercise 2 : Think back to exercise 1. Did you have to add together several measurements with the meter stick to get the length of the table? How would you report the measurement now? ________________ ± ________________ 2) Products and Quotients : If the desired quantity, ? , is the found by a product of two (or more) measurements ? and ? (e.g. area = product of length and width) then the relative uncertainty of ? is ∆?/?, where ? is the best estimate and ∆? is the absolute uncertainty and is given by, (∆?/?) 2 = (∆?/?) 2 + (∆?/?) 2 + ⋯ The exact same formula applies for quotients of measurements. Note ? and ? are best estimates and their absolute uncertainties are ∆? and ∆? , so the ratios are relative uncertainties. Example: (8.0N ± 1.0 N) / (2.0 m ± 0.2 m) = = 4.0 N/m ± [ √(1.0/8.0) 2 + (0.2/2.0) 2 ] = 4.0 N/m ± [0.16] 3) Quantities raised to Powers (not exponentials): If the desired quantity ? is given by raising the measurement ? to the power ? gives, ? = ? ? , then the relative uncertainty of the desired quantity is given by, (∆?/?) 2 = (|?|∆?/?) 2 Examples: (5.0cm ± 0.3cm ) 3 = (5.0cm ) 3 ± [ √3 2 (0.3/5.0) 2 ] =125 cm 3 ± [0.18] = 125 cm 3 ± 22.5cm 3 (9.0m 2 /s 2 ± 1.0m 2 /s 2 ) 1/2 = 3.0m/s ± [ √(0.5) 2 (1.0/9.0) 2 ] = 3.0m/s ± [0.055] = 3.0m/s ± 0.17m/s Exact numbers, like 𝝅 or numbers not from measurements, have no uncertainty.
8 Exercise 3 : You and your partner should each measure the length and diameter of a cylinder. From these measurements, calculate the volume and the uncertainty in the volume using error propagation. Make sure you use the rules on the previous page. Show your work and report your measurement of the volume, ________________ ± ________________ Checkpoint: Raise your hand and show the results of these three exercises to your TA. Make sure you can explain the steps of your calculation and justify any assumptions. If you have to wait for the TA, you may read on. How do you find the uncertainty for a given measurement? Often the measurement device has a given resolution or uncertainty. For example, an analytical balance may have an uncertainty of 1 𝜇g . There may be other times that the uncertainty of a measurement may not depend solely on the uncertainty of the instrument. Consider using a stopwatch to time an event. The time resolution of the stopwatch is typically at least 0.01 seconds. Unfortunately, a person’s reaction time is less certain. One way to estimate the uncertainty in the timing of the event is to repeatedly time the event and use either the standard deviation or the maximum deviation for the uncertainty. Before discussing accuracy, you will review some of the features of Excel which will simplify your calculations and statistical tests of your results. You will also use Excel to calculate standard deviations for sets of measurements. These can be used to plot error bars on a graph. Excel Summary (based on Office 365 education version) Open a new workbook in Excel and follow along to learn some of the basic ways we will be using Excel in this course. Entering an equation For example, say you wish to multiply the entries in B1 and C1 and put the product in D1 (see below). You can select B1 and type in a 2 and similarly put a 3 in C1. You can then select the D1 cell and type = B1*C1 and hitting Enter would give the result of 6. You could also select D1 and type = then click on the cell B1, type * , then click on the cell C1 and hit return to accomplish the same.
9 Automatically filling a column Click on the cell that contains the equation you wish to fill into the cells below it. Move the cursor to the bottom right-hand corner of the cell (black square), click and hold on that corner, and drag down the appropriate number of cells. For example, in A1 we typed a 1 . In A2 we typed =A1+1 and hit return. We then dragged the corner of A2 down through A6 with cell values being increased by 1 from the previous cell as shown in the adjacent screenshot. This is a convenient way to make an ? -coordinate for some function. See the following screenshots. In the last screenshot, B1 contains the equation =A1^2 that has also been filled down to B6 (i.e. ?(?) = ? 2 ) As seen in the above examples, Excel typically uses relative referencing when doing calculations. When you drag down a formula, Excel assumes you want to use a relative position for the data to
Your preview ends here
Eager to read complete document? Join bartleby learn and gain access to the full version
  • Access to all documents
  • Unlimited textbook solutions
  • 24/7 expert homework help
10 calculate the new result. For example, in the above screenshots, the cell A3 contains the formula =A2+1 and the cell B6 contains the formula =A6^2 . Absolute referencing Suppose you have a certain value (say a constant like ? ) that will be used in many calculations, and you want to type it in once, and then reference it after that. Putting dollar signs in front of the row and the column of the cell with this constant tells Excel to not use relative referencing. To use an absolute reference to the cell C1 in a formula, instead of typing C1 you type $C$1. In the screen shot below, a constant value of 9.8 was typed into C1 and used in the formula in column B which now represents the function ?(?) = 9.8? 2 . Standard deviation Standard deviation describes of how far a set of measurements varies from the mean or average of the measurement. In cells D1 thru D5 type in 5 random numbers that are close to 10, e.g. 9.875 (you might imagine these are 5 measurements of something). In D6, type =average(D1:D5) and hit enter to get the average of these measurements . In D6, type =stdev(D1:D5) and hit enter to get the standard deviation . Change some of the numbers in D1 thru D5 and see how the average and standard deviation changes. It is desirable that an experiment yield a standard deviation that is similar in size to the estimated of the uncertainty in an experiment found via the error propagated from the known uncertainties of the measurement devices. If the standard deviation is much larger than the estimated uncertainty it is likely there is unknown systematic error, i.e. a mistake. If the standard deviation is much smaller than the estimated uncertainty then you should be suspicious of the results. If these numbers represented measurements of a 10kg mass on a scale with a known uncertainty of 0.001 kg, what would you think about this set of measurements (Was it a good experiment? Should you redo it? Should you be suspicious of the results?)
11 Accuracy and Judging Scientific Results A peer-reviewed scientific article has been reviewed by other scientists who do research in the same area. They use their experience and knowledge of the topic to judge its suitability to be published. They also expect to see statistical arguments about the experimental results in the paper. One common test of the validity of the results is known as a chi-squared ( 𝜒 2 ) test. It is used to determine how the observed experimental results differ from the expected results. The assumption that the expected results are true is known as the “null hypothesis”. In equation form, it calculates the square of the difference between a particular observed result ( ? 𝑖 ) and the expected result ( ? 𝑖 ) divided by the expected result. It sums this over all the different possible results, 𝜒 2 = ∑ (? 𝑖 − ? 𝑖 ) 2 ? 𝑖 ? 𝑖=1 (Eq. 1) Large values of 𝜒 2 indicate a large deviation from the expected result or null hypothesis. The likelihood of finding a particular 𝜒 2 in a set of measurements can be found in tables or calculated in Excel. This likelihood is known as the ? -value. A ? -value close to 1 indicates that the measured results support the expected results, conversely a small ? -value suggests the expected result or null hypothesis may need to be rejected. In this exercise, you are going to look at the likelihood that a die is fair. You are checking a six- sided die (so ? = 6 ) and you assume it is fair (the null hypothesis) so that each face is just as likely to appear. The die is rolled ? times you would expect each face to appear ? 𝑖 = ?/6 times. Of course, if ? was a relatively small number of rolls the observed number of times a face ( ? 𝑖 ) will appear will deviate from this expected outcome. The ? -values (or likelihood that the observed or measured results confirm the expected results) are a function of both the value of 𝜒 2 and the number of degrees of freedom of the system. In this example, the number of degrees of freedom is 5 since if you know what happens with 5 of the faces of the die the last face is determined. Similarly, for a flipped coin, the number of degrees of freedom is 1. Below is a table of ? -values for different degrees of freedom and 𝜒 2 values. Table of 𝒑 -values as a function of degrees of freedom (rows) and 𝝌 𝟐 values (columns). d.f. / 𝜒 2 0.1 0.5 1 2 4 8 16 1 0.7518 0.4795 0.3173 0.1573 0.0455 0.0047 0.0001 2 0.9512 0.7788 0.6065 0.3679 0.1353 0.0183 0.0003 3 0.9918 0.9189 0.8013 0.5724 0.2615 0.0460 0.0011 4 0.9988 0.9735 0.9098 0.7358 0.4060 0.0916 0.0030 5 0.9998 0.9921 0.9626 0.8491 0.5494 0.1562 0.0068 For example, a 𝜒 2 = 1 for a die experiment would support that the die was fair at the level of 96.26%. However, the same 𝜒 2 = 1 for a flipped coin experiment, would support that the coin was fair only at the level of 31.73% (you should be suspicious). Exercise 4 : Your TA will give you a data set where the die you are testing for fairness was rolled 1500 times. The expected number of times each face should appear is 250 or ? 𝑖 = 250 . The
12 number of times each face did occur in these 1500 rolls will be given to you by your TA. Using Excel, calculate 𝜒 2 for these measurements (i.e., implement Eq. 1). Once you have found it, you can find the ? -value with the Excel function CHISQ.DIST.RT( 𝜒 2 ,df). The first argument of this function is the value of 𝜒 2 that you found and the second is the degrees of freedom for the die. 𝜒 2 = ________________ ? -value =________________ What is the likelihood that the die tested was fair? Should you be suspicious? Checkpoint: Raise your hand and show your spreadsheet and your results to your TA. Make sure you can explain the steps of your calculation and justify any assumptions. If you have to wait for the TA, you may read on. Excel Graphing After creating, entering, or importing the data you want to plot, drag the mouse to select the data you wish to graph. After selecting the data, click on the Insert tab and find the Chart section (see below). There is tiny picture of a ? − ? graph. If you click on it, there are several possible types of graphs. Below we chose the simp lest “X - Y (Scatter)” with just points. We also selected, under Quick Layout, Layout 1. We deleted the standard label for the ? axis and replaced it with “ ? ” and similarly on the ? axis we typed in ?(?) .
Your preview ends here
Eager to read complete document? Join bartleby learn and gain access to the full version
  • Access to all documents
  • Unlimited textbook solutions
  • 24/7 expert homework help
13 Fitting your data to a line or other function If you want to add a regression line, the “+ box” next to the graph contains a check box for a “Trendline” click that and a best fit line will be drawn on the graph. The arrow next to the Trendline box offer s you More Options. Click “More Options" and near the bottom of the resulting menu, click on the “display the equation” for the fitted line and “display R - squared value” (a measure of how good the fit is; closer to one is better). Often you may wish to plot experimental data vs. a theoretical prediction. The expected graph would be a line with the slope of one (assuming data and theory agree). The linear fit to this line and the correlation are a measure of how well the theoretical prediction fits the data. Note there are also other functions available to fit your data in Excel. y = 68.6x - 91.467 R² = 0.9583 -100 0 100 200 300 400 0 2 4 6 8 f(x) x
14 Uncertainty (Error) Bars To add error bars to the graph, the “+ box” next to the graph contains a check box for “Error Bars”. Check the box and some rather arbitrary error bars appear on each data point. Click on one of the vertical error bars on the graph. An Error Bars menu opens. If all the uncertainties are the same for each data point, choose fixed value and enter the uncertainty in the adjacent box. If the uncertainty varies with each point you can choose Custom. Click on the box “Specify Values” and a small dialog box opens (see below). The example below adds error bars to the ? data points (in column B) using the uncertainties given in column D. Click on the box below “Positive Error Value” and highlight D1 through D6. Then do the same for the box below “Negative Error Value”. The ? error bars are done in a similar fashion using the data in Column C. Exercise 5: Generate a set of data for an object moving in a way described by the formula ?(𝑡) = 𝑡 2 ? −3? where 𝑡 is the time in seconds and the position ? is meters. Create the time data over the range from 0 to 3 seconds at intervals of 0.1s. Assume the ? position can be measured with a resolution of 1cm (same value for all data points). Create a graph of ? vs 𝑡 with error bars. See details below for creating a properly formatted graph. Checkpoint: Raise your hand and show your properly formatted graph to your TA. It should have the following features: The ? and y axes should be labeled with a descriptive and units, e.g. for the ? axis “time (s)” and for the ? - axis “position (m)”. The font size should be reasonable (try 14 pt and then adjust).
15 The scales of the ? and ? axes should be such that the data fills the graph window well. The font size of the numbers should match that of the labels. Vertical error bars should be displayed for each point. You may assume the time is known with high resolution and so no horizontal error bars are needed. There is no lab report due for this experiment. If you have to wait for the TA to check your graph , begin reading Part 1 of next week’s activity.
Your preview ends here
Eager to read complete document? Join bartleby learn and gain access to the full version
  • Access to all documents
  • Unlimited textbook solutions
  • 24/7 expert homework help
16 Lab 1. Scientific Literacy and Measurement Apparatus In this lab, you will find and select a peer-reviewed scientific article that your group will read and later, during Lab 3, give an oral presentation on the article. This lab will also give an overview of the available measurement devices in this laboratory and how they function. Part 1: Finding Scientific Articles The UT library system has subscriptions to a wide number of scientific journals. They can be accessed at https://www.lib.utexas.edu/find-borrow-request . Click on Journals> . Type in the journal name (see screenshot below with example of American Journal of Physics). Click on Available Online and you will get the following result. You can search for specific articles or Click on the box with the arrow in it to browse current articles.
17 Your group will be using this general method today to find an article to do a short presentation in a later lab session. The goal is two-fold: to help you become familiar with accessing peer-reviewed scientific articles and to give you a chance to do an oral presentation on a scholarly topic. There are a few peer-reviewed journals that have articles that are meant for a more general audience, i.e. more easily understood by a non-specialist. These include the American Journal of Physics, Nature, and Science. Your task is to find an article in one of these that catches your group’s interest and during the next two weeks prepare a 10 minute talk discussing that article (what it was studying, what it found, how it found it, and its significance). The talk should use Powerpoint slides and each member of the group should talk equally. An example of what the slides should look like for such a talk can be found on Canvas. Checkpoint: Raise your hand and show the article you have chosen to the TA. Be prepared to explain why you have chosen it. You may proceed to the other part of this activity if you must wait for the TA. Part 2: Using the measurement apparatus Measuring distance In lab 0 you used a meter stick or ruler to measure length. This will be our standard way of measuring distance and position. Note the uncertainty associated with this measurement is half of the smallest division on the ruler or meter stick e.g. for a ruler with 1mm divisions the uncertainty is 0.5mm. Typical ruler uncertainty: 0.5mm
18 Measuring time Watches and smartphones typically come with a stopwatch app. Find it and familiarize yourself with how you stop, start, and reset it. The fundamental timing circuits on these devices are highly precise (uncertainties of less than 10 −6 s). However, they typically only report to the 1/100 th of a second and your reflexes on stopping and starting the stopwatch are even worse. Determining stopwatch uncertainty Each member of the group should start their stopwatch and try to stop it exactly on 2 seconds. Repeat this 5 times. Put the resulting exact times for all the trials by the group into a spreadsheet and calculate the standard deviation. Record it below. Stopwatch uncertainty :__________________ Checkpoint: When you have finished, raise your hand in the Zoom app to let the TA know you are ready to be checked out. Make sure you can explain how you did the measurement and justify your result. While you wait you can proceed to the next section. This lab has general apparatus (lab stands, rods, string, etc.) and specialized sensors. They are found in the cabinets in the laboratory room. Instructions for using the different sensors are below. The Pasco Motion Sensor is an ultrasonic (Sonar) probe. It transmits a pattern of ultrasonic pulses which strikes an object. The transmitter is also used as the receiver of the reflected ultrasonic pulses. The software compares the reflected pattern to the transmitted pattern to determine the distance from the sensor to the object. The specifications of the sensor are below: Range 0.15 to 8 m Resolution 1.0mm Maximum Sample Rate 50 Hz Pasco Motion Sensor
Your preview ends here
Eager to read complete document? Join bartleby learn and gain access to the full version
  • Access to all documents
  • Unlimited textbook solutions
  • 24/7 expert homework help
19 Quick Start 1. Connect the Motion Sensor to your PASPORT 850 Universal Interface 2. Turn on the interface and make sure it is connected to the computer. 3. Start the Capstone Software and Click one of the available templates. The selected layout appears in the Display Area. 4. Click <Select Measurement> for each display. 5. Select the measurement in the drop-down list e.g. Motion Sensor. 6. Place an object in front of the sensor at least 15 cm away. 7. Click “Record”. 8. Move the object in a straight line directly away from or toward the sensor (recall the results may not be reliable if the object get closer than 15 cm). Click “Stop”. 9. A graph of the Position vs. time is usually the default. By clicking the label of the vertical axis you may select Velocity or Acceleration from the Measurements Menu. Exercise 1: Find the Help tab at the top of the Pasco Capstone window. Use it to figure out how to export your data and then open it in Excel. Make a plot of the position vs. time in Excel (make sure it is properly formatted as described in the previous lab). Checkpoint: When you feel you have learned how to use this sensor let the TA know you are ready to be checked out. You can begin reviewing how to use the next sensor while you wait. Pasco Rotary Motion Sensor
20 The Pasco Rotary Motion Sensor measures the angular motion of a rotating axle. A wheel with fine angular divisions is attached to the axle and light is reflected off the wheel. The variation of the light as the markings pass is converted into angular position data. This wheel and light sensor package is known as an optical encoder. Note the linear motion of a string passing over the pulley can also be found to a high precision. The specifications of the sensor are below: Three-step Pulley: 10, 29 and 48 mm diameter Resolution: ±0.09° or 0.00157 radian or 4,000 divisions/revolution o Linear resolution for 10 mm pulley is 7.85 × 10 −6 m. Maximum Rotation Rate: 30 revolutions per second Optical Encoder: Bidirectional, indicates direction of motion. Quick Start 1. Connect the Rotary Motion Sensor to your PASPORT 850 Universal Interface 2. Turn on the interface and make sure it is connected to the computer. Start the Capstone Software. 3. Start the Capstone Software and Click one of the available templates. The selected layout appears in the Display Area. 4. Click <Select Measurement> for each display. 5. Select the measurement in the drop-down list e.g. Rotary Motion Sensor. 6. Click or press ‘Record’ to begin recording data. Give the pulley assembly on the sensor a spin. Once it has stopped spinning, Click “Stop”. 7. A graph of the Angular Position vs. time is usually the default. By clicking the label of the vertical axis you may select Angular Velocity or Angular Acceleration from the Measurements Menu. Checkpoint: When you feel you have learned how to use this sensor let the TA know you are ready to be checked out. You can begin reviewing how to use the next sensor while you wait. The Pasco Photogate is an infrared light source mounted in one arm of the U-bracket aimed at a photodetector in the other arm. If something passes between the arms the resulting voltage signal falls from high voltage to low voltage when the light blocked. It can be used to measure the time it takes for an object to pass through the gate (e.g. time from front edge to back edge). Often it is Pasco Photogates
21 used in pairs of photogates to measure the time between something passing through first one gate and then the other. It can also be used for oscillatory motion. Time resolution (the fastest the voltage can fall from high to low) as an edge rapidly blocks the light is 0.1 ms. This is our fastest sensor in the lab. It is, however, a bit more complicated to setup in the software. Quick Start 1. Attach the Photogate Head to a support rod or lab stand. 2. Turn on the interface and make sure it is connected to the computer. Start the Capstone Software. 3. Click on Hardware Setup near the top of the left column. 4. You will be using one of the Digital Channels 1 2 3 4 on the left. 5. Plug the cable from your photogate into Channel 1 into the and select that channel on the computer. 6. From the pop-up menu, select the Photogate. (If no menu appears, you can make the selection from the Add sensor/interface button) 7. When a Photogate is selected, a new icon will appear on the left in the tools palette (looks like photogate) 8. Click on this new icon. A time wizard menu will expand (there are many options). There are 6 steps to the wizard. Just Click next to proceed to step 2. Make sure the Photogate to be used is checked. Click next. In step 3, Choose Pendulum Timer. In step 4, Period and Speed should be checked. In step 5, enter the width of the object you are going to pass through the gate (e.g. pencil or ruler).and select Setup Timer. In step 6, you can choose the name of the timer. Click Finish. 9. Once the timer is set up, prepare the display area. Select Table from the right column of icons and drag it into the white display area. Note that Table has a row of tool icons at its top. They are used to control the Table. Your Table needs three columns. Click on the icon at the left end of the row to add a column. 10. Select a measurement for each column by clicking on the entry box. For the left column, Select Time(s). For the remaining column(s), select Period(s), Speed (m/s), 11. Select Continuous Mode. 12. Select Record to begin data acquisition. Try to move your object periodically back and forth through the gate. Click Stop. Checkpoint: When you feel you have learned how to use this sensor let the TA know you are ready to be checked out. Make sure you can explain how it is finding the period and speed. You can begin reviewing how to use the next sensor while you wait.
Your preview ends here
Eager to read complete document? Join bartleby learn and gain access to the full version
  • Access to all documents
  • Unlimited textbook solutions
  • 24/7 expert homework help
22 Using the Sensors on your SmartPhone There are a number of free Apps that give you access to the sensors on your smartphone. Some good examples are Physics Toolbox Sensor Suite and Phyphox. Depending on your type of phone (Android or Iphone) you may find one better than the other. They both provide 1) access to data from the accelerometers in your phone; 2) an audio generator; 3) an audio oscilloscope; and 4) A spectrum analyzer. Primarily, we will be using the accelerometers. The accelerometers on a wide variety of phones have similar specifications, i.e. the resolution of the acceleration measurement is 0.0025? , where ? is the gravitational acceleration on the surface of the earth. In this review of this apparatus, we will only look at the accelerometer. Download the appropriate App and open it. Find the g-Force meter or Acceleration or whatever it is called on the app you chose. You may wish to look at a couple different apps and choose whichever you and your group finds easiest to use. Start the accelerometer. Move the phone around a bit to see what happens. Explore the different settings options. Clear the data. Let’s try accelerating the phone in only one dimension. Place it on the surface of your lab table and slide it back and forth in a direction parallel to a side of the phone (it may be helpful to use a meter stick as a guide rail). After practicing this motion, start the accelerometer (for Toolbox tap the + button to start recording data) then start the motion. Stop the recording the data and then export it as an Excel file via email to yourself (or if you prefer to all work on the same file also to you partners). Open the data in Excel and plot the magnitude of the acceleration vs. time. You will want to plot just the data from when you start moving the phone till when it is stopped moving. Make sure you all have roughly the same plot. Checkpoint: When you have created this plot (make sure it is properly formatted as described in the previous lab), let the TA know you are ready to be checked out. Make sure you can explain what is occurring in the plot. What is being plotted? How do significant features in the plot correspond to certain motions of the phone? Discuss with your partners: How you can use your acceleration vs. time data to obtain velocity vs. time data? Position vs. time data? In Excel, how can you integrate a function? There is no lab report for this lab.
Your preview ends here
Eager to read complete document? Join bartleby learn and gain access to the full version
  • Access to all documents
  • Unlimited textbook solutions
  • 24/7 expert homework help
23 Lab 2. Experiment on motion in 1D In this activity, you will measure the motion of a particle in one dimension. Warmup: Modeling Motion with Calculus and Motion Matching The functional representations of the position, velocity, and acceleration of a particle ( ?(𝑡), 𝑣(𝑡), and ?(𝑡) ) can all be connected by simple calculus. The time derivative of the position is the velocity, ? ?? (?(𝑡)) = 𝑣(𝑡) , or conversely the integral with respect to time of the velocity is the position, ?(𝑡) = ∫ 𝑣(𝑡′)?𝑡′ ? 0 . The second derivative with respect to time of the position is the acceleration, ? 2 ?? 2 (?(𝑡)) = ?(𝑡) . The integral with respect to time of the acceleration is the velocity, 𝑣(𝑡) = ∫ ?(𝑡′)?𝑡′ ? 0 . Exercise 1: The derivative of a function at a point is the slope of that function at that point. Below we have a graph of the position vs. time of a particle. On a separate piece of paper, copy this position vs. time graph and then below it graph the velocity of the particle vs. time, and below that graph the acceleration of the particle vs. time (in a column making sure to align the time axes). 0 0.05 0.1 0.15 0 0.5 1 1.5 2 2.5 position (m) time (s)
Your preview ends here
Eager to read complete document? Join bartleby learn and gain access to the full version
  • Access to all documents
  • Unlimited textbook solutions
  • 24/7 expert homework help
24 Exercise 2: In the graph below, you are given the acceleration of a particle as a function of time. On a separate piece of paper, copy this acceleration vs. time graph and then below it graph the velocity of the particle vs. time, and below that graph the position of the particle vs. time (in a column making sure to align the time axes). Recall that the integral of a function is the area under the curve; e.g. the velocity at a time 𝑡 is the area under this acceleration curve from zero to 𝑡 . Warmup experiment : You will try to create some particular ?(𝑡), 𝑣(𝑡), and ?(𝑡) . The idea is to gain some kinesthetic sense of what these relations mean. Secure a Pasco Motion Sensor (Sonar type) to a lab rod and plug the sensor into the interface. Load the data acquisition program. Adjust the angle of the transmitter/receiver so that you can move a book (or lab manual) over some distance (recall the sensor works over a range of 0.15m to 8m from the sensor). Move the book and make sure the sensor is recording the motion of the book. Try to do the following (for the checkpoint you will be asked to explain the results to the TA for one of these): 1. Move the book at constant velocity towards the sensor and then away from the sensor. Examine the associated ?(𝑡) and ?(𝑡) . You should be able to explain the results to your TA. 0 0.05 0.1 0.15 0 0.5 1 1.5 2 2.5 acceleration (m/s 2 ) time (s)
Your preview ends here
Eager to read complete document? Join bartleby learn and gain access to the full version
  • Access to all documents
  • Unlimited textbook solutions
  • 24/7 expert homework help
25 2. Move the book at constant acceleration towards the sensor and then away from the sensor. Examine the associated ?(𝑡) and 𝑣(𝑡) . You should be able to explain the results to your TA. 3. Hold the book stationary for a duration of 2 seconds at each of the following three positions from the face of the sensor (0.2m, 0.5m, 1.0m). The transitions between the positions should be rapid so that the position curve ?(𝑡) resembles steps. Examine the associated 𝑣(𝑡) and ?(𝑡) . You should be able to explain the results to your TA. Checkpoint: Show your graph for the 3 rd case to your TA and be able to explain the experimental results to your TA. If you are waiting for the TA, proceed to the next section and begin to design your experiment. Motion of object with constant acceleration Design an experiment where an object moves with a constant acceleration. The constant value must be between 2m/s 2 and 5m/s 2 . You can choose to measure that motion with any of the available sensors (Sonar Motion Sensor, Rotary Motion Sensor, Photogates, Accelerometer on phone). The experiment must safe for both you and your sensor. For example, slamming your phone into a hard surface could break it and make the experiment unrepeatable or trying to run with an acceleration of 5m/s 2 is not safe in the lab (hitting walls or other students). Have your TA okay your design for safety before performing the experiment. The experiment and analysis should: Measure the motion over a significant time span (~ 1 second) where the acceleration appears to be constant and has a value in the desired range. Find, during that time span, the average acceleration and the standard deviation of the acceleration from that average value. Calculate, from the uncertainties in your measurements, the uncertainty in the acceleration. Look back at Lab 1 for the uncertainties of the various measuring devices and refer to Lab 0 (page 7) for how to calculate the uncertainty in acceleration from the uncertainties of the measuring devices that you used. Compare this to the standard deviation found above. Is the standard deviation larger or smaller than the calculated uncertainty? Use your measurements, to plot ?(𝑡) using Excel and calculate and include vertical Error Bars on the graph. Label the axes with the quantity on that axis and the appropriate units, e.g. for the horizontal axis, time (s) . Finally, compare your results with those of other lab groups. Are your standard deviations smaller/larger? Are your uncertainties smaller/larger (calculated via error propagation)? Are your error bars smaller/larger?
Your preview ends here
Eager to read complete document? Join bartleby learn and gain access to the full version
  • Access to all documents
  • Unlimited textbook solutions
  • 24/7 expert homework help
26 The abbreviated lab report due for this experiment should have: an Abstract (200 words or less) summarizing your experiment. Its first sentence should describe what the experiment is designed to do. The next two sentences should highlight key points in the design of your experiment. The last should give the most significant experimental result(s). A Results section giving the average acceleration and the standard deviation of the acceleration from that average value. It should contain the uncertainty in the acceleration and how that uncertainty was determined. and the plot of acceleration vs. time with error bars. The plot should include a linear fit to the data (displaying the equation for the linear fit). The plot should be properly formatted as described in Lab 0, Exercise 5. A Conclusions section stating the value of the constant acceleration and the uncertainty of that value. Compare the standard deviation of the acceleration to the uncertainty in the acceleration that you calculated. Is the object undergoing constant acceleration to the limits of your measuring device? What does the slope of the linear fit tell you about the data? Comment on any variation from constant acceleration you see in your data and suggest possible causes. Print out your lab report and exchange it with another group and have them check that your report contains all the required elements. If you need to wait, make preparations for next week’s presentation (see reminder below). Make any necessary revisions to the report and turn it in to your TA. The grading rubric for this lab report is found on the following page. REMINDER: Next week you will be presenting the findings in the paper your group selected in Lab 1. A sample Powerpoint presentation can be found on Canvas. A copy of the feedback form that will be used to critique the presentation can also be found on Canvas. In preparation, practice your presentation in front of friends or family or you may video and watch yourself to polish your delivery. Make sure that all parts of your group’s presentation fits in the 10 minute time window. For example, for a group of two, each person should talk about 5 minutes. Make sure the parts of the talk fit together smoothly. Your TA will have created a UT Box folder and shared it with you. Please upload your completed presentation to that folder. Please include one of your lastnames in the filename so you can easily identify it.
Your preview ends here
Eager to read complete document? Join bartleby learn and gain access to the full version
  • Access to all documents
  • Unlimited textbook solutions
  • 24/7 expert homework help
27
Your preview ends here
Eager to read complete document? Join bartleby learn and gain access to the full version
  • Access to all documents
  • Unlimited textbook solutions
  • 24/7 expert homework help
28 Rubric for Lab 2 Report If they have all the required elements in their report, their grade should be 100. The abbreviated lab report should have: an Abstract (200 words or less) summarizing your experiment. (if over 200 words minus 5) o First sentence should describe what the experiment is designed to do. (10 pts) o The next sentences should highlight key points in the design of experiment. (10 pts) o The last sentences should give the most significant experimental result(s). For example, “The acceleration was found to be constant to within 5% of 3 m/s 2 . This variation was larger than the uncertainty in the measurement”. (10 pts) A Results section giving o the average acceleration (5 pts) o the standard deviation of the acceleration from that average value. (5 pts) o It should contain the uncertainty in the acceleration and how that uncertainty was determined. (15pts) o A correctly formatted plot (10pts) of acceleration vs. time with error bars (5pts). The plot should include a linear fit to the data (displaying the equation for the linear fit). (5pts) A Conclusions section with o the value of the constant acceleration (5pts) o the uncertainty of that value (5pts) o Comparison of the standard deviation of the acceleration to the uncertainty in the acceleration that you calculated. Is it constant acceleration to the limits of the measurement? (5pts) o What the slope of the linear fit says about the data (slope can indicate general increase or decrease of acceleration) (5pts) o Comment on any variation from constant acceleration you see in your data and suggest possible causes. (5pts) If no variation is present give these 5 pts.
Your preview ends here
Eager to read complete document? Join bartleby learn and gain access to the full version
  • Access to all documents
  • Unlimited textbook solutions
  • 24/7 expert homework help
29 Lab 3. Scientific Presentation During this lab session: Your group will be presenting the findings in the paper you selected in Lab 1. o You will be expected to be able to answer a couple of questions from the audience. When you are not giving your presentation, you will be providing feedback on the other groups presentations. A copy of the feedback form is at the end of this activity. You are expected to ask a question on some of the talks. Your TA will use the same feedback form as a rubric to grade your presentation. Each check box is worth 10 points. Modelling Physics In the next activity, you will build a computational model of a particle’s motion. We will use Python as our programming language and, in particular, Vpython which offers 3D visualization. To prepare for modeling the motion, you will look at representing physical quantities with vectors (arrows for velocity or force) in this language. Our use of Vpython is focused on precisely how to model the physics of the problem. If you would like to learn Python in more depth, there are many good tutorials available on the web e.g., https://docs.python.org/3/tutorial/ Using Vpython and Vectors Web VPython is a browser-based version of Vpython which has all the functions needed for this modeling (all the instructions will assume that you are using this version). Go to vpython.org on your browser and Click on “Using VPython without installing any software”. Under the “ Web VPython: Log in, write, run, share, export programs stored in the cloud section click on “Read more” and follow the instructions for getting started with Web Vpython (making a free account, making a test program, etc.). Once you finished these preliminaries, proceed to the next section.
Your preview ends here
Eager to read complete document? Join bartleby learn and gain access to the full version
  • Access to all documents
  • Unlimited textbook solutions
  • 24/7 expert homework help
30 Vectors and Vector Operations Position Vectors Position vectors describe where something is with relation to the origin of a coordinate system. They can be represented pictorially by arrows which point from the origin to the position of an object (see figure). Type the following two lines into a new program in Web VPython . They will load the visualization library and define the size, location and color of the graphics window where you can draw your vectors and anything else you would like to model. from visual import * scene = display(width=600, height=600, center=vector(0,0,0), background=color.white) The display window or scene is initially oriented with the positive ? axis pointing to right, the positive ? axis pointing to the top of the page, and the positive ? axis pointing towards you out of the page. The center of the scene is at the origin. VPython works directly with vectors and defines a vector in terms of an ordered triple (x, y, z). The next lines of code define the origin of the coordinate system to be (x=0, y=0, z=0) and also define two vectors A and B . These ordered triples define the locations of the tips of the arrows for the vectors A and B . Type these lines of code into your program below the your first two lines. # The origin and tips of vectors A and B origin=vector(0,0,0) A=vector(1,1,0) B=vector(1.5,-1.5,0) The arrow function is used to create arrows to represent these position vectors. The tails of the arrows are both at the origin of the coordinate system. This is indicated by the pos variable in each arrow. The axis variables contain the locations of the tips of the position vectors. You can choose to make the color of the arrows different, e.g. blue and green. The width of arrow is defined by the variable shaftwidth . Note anything after a # is a comment and does no operation. It does provide helpful information regarding what the code is meant to do. # Draw A and B on the coordinate system Aarrow=arrow(pos=origin, axis=A, color=color.blue, shaftwidth = 0.1) Barrow=arrow(pos=origin, axis=B, color=color.green, shaftwidth = 0.1) Add these lines of code, run the program and you should see two arrows. Check that the positions of the vectors make sense. x y origin
Your preview ends here
Eager to read complete document? Join bartleby learn and gain access to the full version
  • Access to all documents
  • Unlimited textbook solutions
  • 24/7 expert homework help
31 Using the 3D Capabilities of VPython When you first run the program, the scene shows you the x-y plane. Your viewpoint is looking toward the origin along the z-axis (the negative z-axis goes into the computer screen). You can change that with your mouse. Holding down both buttons and moving the mouse up zooms your viewpoint inwards. Similarly, moving the mouse down zooms your viewpoint outwards. Holding down the right mouse button alone and moving the mouse in any direction rotates your viewpoint around the origin. Relative Position Vectors How do you get from B to A (i.e., from the tip of vector B to the tip of vector A )? A vector that has its tail at B and its tip at A is known as a relative position vector. This relative position vector is the difference, r = A - B. Implement code that will draw this new vector, r = A - B . You will need two lines after the Barrow line. The first sets r equal to A-B ; r = A - B . The second is an arrow function for r e.g. rarrow=arrow( … ) filling in the appropriate characteristics in these parentheses. Note the tail of the arrow representing this new vector r must sit on the tip of the arrow representing B. Run your code and make sure the three arrows are oriented correctly. Checkpoint: Ask your TA to check that your results look correct. How would your results change if you added the two vectors, A + B? You can try this if you have to wait for the TA. Make sure everyone in the group has a copy of this program. You will be using it in the warmup next week.
Your preview ends here
Eager to read complete document? Join bartleby learn and gain access to the full version
  • Access to all documents
  • Unlimited textbook solutions
  • 24/7 expert homework help
32 Feedback form for Presentation Names of Presenters ____________________________________________________________ Required Content Title slide (Title of Article and Journal, Authors of Article, Names of Presenters) Introductory slide(s) (Explains significance of the main topic of the Article) Procedure slide(s) (How study was accomplished; sketches of apparatus, etc.) Results slide(s) (Presents the main results of paper, including graphs, tables, etc.) Conclusion slide (Explains significance of results and indicates what can be studied in the future). Presentation or Delivery Clear presentation (understandable; scientific jargon defined) Appropriate speed (spoke at pace that could be followed) Fit time well (used 10 minutes of length well) Spoke to audience (Made eye contact) Answered questions clearly Constructive Feedback on Presentation (One thing that would make the talk better)
Your preview ends here
Eager to read complete document? Join bartleby learn and gain access to the full version
  • Access to all documents
  • Unlimited textbook solutions
  • 24/7 expert homework help
33 Lab 4. Modeling Motion In this activity, you will build a computational model of a particle’s motion. We will use Python as our programming language and, in particular, Vpython which offers 3D visualization. In the warmup, to prepare for modeling the motion, you will look at representing physical quantities with vectors (arrows for velocity or force) in this language. It is focused on precisely how to model the physics of the problem. If you would like to learn Python in more depth, there are many good tutorials available on the web e.g., https://docs.python.org/3/tutorial/ Warmup: Making Particles move with Vpython Web VPython, a browser-based version of Vpython, can be found at webvpython.org Login to the account you created last week and open the program you created. You will be adding code to make an object move. Adding an Object to Move and Making it Move The most common object you will work with is a sphere. To use a sphere to represent a yellow ball sitting at the tip of the B vector the code is, ball=sphere(pos=B, radius=0.2, color=color.yellow) Add this line into your program after the lines that you had written last week. The variable pos is where the sphere is located and is a vector. The radius of the ball is 0.2 m. To define the velocity of the ball you need a magnitude (or speed) and a direction. If you want the ball to go from the tip of B to the tip of A , you want it to move in the direction of the r vector. We can define that direction as a unit vector where we take r and divide it by its magnitude (gets rid of the units and leaving only direction and a magnitude of 1). Conveniently, there is a function mag() which gives the magnitude of a vector, e.g. mag(r) gives the magnitude of r. In code, unitr=r/mag(r) # this is the unit vector of r If the speed is equal to 0.1 m/s then the ball’s velocity can be written in code as ball.v=0.1*unitr Add the above lines of code below the line where you defined ball. .
Your preview ends here
Eager to read complete document? Join bartleby learn and gain access to the full version
  • Access to all documents
  • Unlimited textbook solutions
  • 24/7 expert homework help
34 You will use a while loop to make the ball move. The code you need is described here and in the flow chart to the right. The code itself is given below. Please take time to understand it since we will use this type of coding often. This code first sets values for the time step, dt , and for the time we want the ball to stop moving, tend . It also sets the time variable, t , to zero. The line while t<tend: starts a loop where all the indented lines of code that follow it will be repeated until t > tend . Each time through the loop the current position of the ball ( ball.pos) is found by adding the ball’s velocity times dt to the previous value of ball’s position. Then the current time t is calculated by adding dt to the previous value of t . Be careful: code like t=t+dt is not an algebraic equation! It means add dt to the old value of the time to get the current time. Since there are no more indented lines, the while statement is executed next and the loop continues until t>tend . Note that the while statement ends in a colon and that all the subsequent lines have exactly the same indentation … coding requires you to pay attention to such details! #Making the ball move from B to A dt=0.005 # a small time step i.e. 0.005 s tend=mag(r)/mag(ball.v) # tend is the time it takes for a ball to travel the distance from B to A t=0.0 # start time t at zero seconds while t<tend: #loop to do the following indented code over and over until t > tend rate(1000) #executes code only at 1000Hz (slows the motion enough to be seen on screen) ball.pos=ball.pos+ball.v*dt #updates the ball position after dt has passed t=t+dt #increase time t by dt Vpython allows you to directly access the position vector of the ball using the variable ball.pos Once you have added all these lines of code, run the program and you should see the ball move from the tip of B to the tip of A. The program ends when the ball reaches A since there are no lines of code after the while loop. Note we calculated tend so the ball would stop at the tip of A. You will use your warmup code below so please do not delete it! Set the values of the: time step, dt , time to end the loop, tend , and set time t =zero While statement if t<tend enter loop If t>tend leave loop Rate statement slows the motion so it is visible Multiply the velocity of ball by the time step and add it to the previous value of the ball position to get the current position Add dt to t to get the current time. Next line of code after loop. Not indented
Your preview ends here
Eager to read complete document? Join bartleby learn and gain access to the full version
  • Access to all documents
  • Unlimited textbook solutions
  • 24/7 expert homework help
35 Modeling motion in the presence of a constant acceleration Exercise 1: Constant acceleration In the warmup, you moved a ball in steps by multiplying a constant velocity by a small time step dt . If there is an acceleration present the velocity will change with time, so you will have to update both the position of the ball and velocity of the ball at each time step. For example, in the warmup code you would first need to define the acceleration vector above the while loop. Let’s assume it is accelerating from the tip of B to the tip of A at 0.1 m/s 2 , so the code is ball.a= 0.1*unitr #vector for acceleration of ball Inside the loop, directly after the ball.pos line, add ball.v=ball.v+ball.a*dt #update the ball velocity after dt has passed Once you have added these lines to the appropriate places in your warmup code, run the program and you should see the ball accelerate along r and keep going (note the scene automatically zooms out to keep the ball in view and so the ball and the arrows grow smaller). Make a copy of your code and call it Projectile on Moon. Exercise 2: Projectile motion on a moon (in vacuum; ? = 𝒈/𝟔 ) You will be modeling the launch of a ball on a moon and finding how far it has gone when it strikes the surface (the range). There are a number of things you may remove from your copied code, e.g. the vectors A , B , and the vectors r and unitr. Make sure to keep the first two lines of the code and the definition of ball as a sphere. You can keep the while loop but since you do not know when it will strike the surface set tend=50000 (a presumably long enough time). You will need to write some additional code to model the launch a ball from the surface of the moon (see below). You will vary the angle of launch to find the result which produces the maximum range. Changes or additions to the code: The ball’s initial position should be the origin. Change this in the definition of ball. The initial velocity of ball.v has a speed 𝑣 0 = 2m/s and an angle, 𝜃 , that will be varied. You can define the velocity vector as the triple (𝑣 0 cos 𝜃, 𝑣 0 sin 𝜃, 0) . In Vpython the sin() and cos() functions assume the angle of launch 𝜃 is in radians. The constant pi is predefined in Vpython, e.g. the code cos(pi/3) would equal 1/2.
Your preview ends here
Eager to read complete document? Join bartleby learn and gain access to the full version
  • Access to all documents
  • Unlimited textbook solutions
  • 24/7 expert homework help
36 Define the acceleration of the ball, ball.a=vector(0.0,-9.8/6.0,0.0). Make sure these variables that are based on ball are defined after ball is defined. You should define an arrow, varrow , representing the velocity of the ball. Its tail should remai n on the ball and the arrow’s direction and magnitude should equal the changing velocity. The statement defining varrow should be before the while loop. However, you will have to update this arrow’s characteristics varrow.pos and varrow.axis inside the time while loop. You should represent the surface of the moon with the following code that uses the box function. Place this code on a line after the scene statement near the top of the program, surface= box(pos=vec(0,0,0), length=25, height=.1, width=0) This puts the center of the box at the origin and makes the length 25 meters. You may need to change the size and position of the surface so you can better observe the trajectory. You can find the range by finding the ball’s position when it hits the surf ace. The ? position of the ball can be found by appending a .x to the vector ball.pos , i.e. ball.pos.x (the components of any vector in Vpython can be called by appending .x , .y , or .z to it). You will need an appropriate test for ending the while loop when the ball hits the surface. One way is to use an if statement to see if the ? -position of the ball has passed through the moon’s surface. If th e condition is true, you set the time to greater than the end of the while loop so the while loop will end. Any indented code after an if statement is executed if the test is true if ball.pos.y < 0.0: #if the ball is below the surface do next line, otherwise skip it. t=tend + 1 # increase t to after the end time so loop ends when it hits ground when you add these lines at the end of the while loop the first line will be at the same indentation as the previous lines, but the second line will be indented even further. There are other ways to code this condition. You could use a similar logical test ( ball.pos.y > 0.0 ) in the while statement. You should add a trail or trajectory by adding two characteristics to the sphere representing the ball . That is after the characteristics already typed inside the parentheses for this sphere, type , make_trail=True, trail_type="points" Print the range (or ? -position of the ball when it hits the ground), use the following code AFTER the while loop (it must NOT be indented). range= str(ball.pos.x) ## converts variable to string for printing in label label( pos=ball.pos+vec(0,-1.0,0), text=range ) ## prints range A string variable is in a code format (ASCII) that can be printed. The ordinary format of a numerical variable is not printable. Run the program. Vary the angle until you find the angle that produces the maximum range.
Your preview ends here
Eager to read complete document? Join bartleby learn and gain access to the full version
  • Access to all documents
  • Unlimited textbook solutions
  • 24/7 expert homework help
37 Checkpoint: Run your program with the angle that produces the maximum range and show it to your TA. Explain your code and the results to your TA. If you are waiting for the TA, make a new copy of your program and begin using it in the next section. Exercise 3: Projectile motion on the earth (air resistance and 𝒈 ) Now you will model the effect of air resistance on the projectile motion on the surface of the earth. Pushing the air out of the way results in a velocity dependent force in the opposite direction of the motion. The acceleration due to the air is often modeled with ? ?𝑖? = −?𝑣 ⃗(𝑡) where ? is a constant and for this problem ? = 0.15 s −1 . Modif y the acceleration in your code to include the change to earth’s acceleration of gravity and to include this ? ?𝑖? . This net acceleration now changes with velocity so it will need to be updated as the velocity changes with time. That is, the acceleration must be updated within the time loop. Add an arrow representing the acceleration of the ball which now changes with time. Increase the launch speed so the trajectory is of a reasonable size (try 𝑣 0 = 20m/s). You may also find it convenient to change the size and center of the surface (box function) and the center of the scene. Run the program. For the given value of ? , vary the angle until you find the angle that produces the maximum range. How does the angle which produces the maximum range change when you double ? ? For a fixed angle of 𝜋/4 , vary the value of ? and qualitatively observe how the trajectory changes. Record this data and your observations below. Angle for maximum range=_________________ ; Range=____________________________ Angle for max. range ( 𝟐? )=_________________ ; Range=____________________________ For a fixed angle of 𝝅 𝟒 , describe what happens to the trajectory as you vary the value of ?. Checkpoint: Run your program with the angle that produces the maximum range for ? = 0.14 s −1 and show it to your TA. Be able to explain your results and your observations to your TA. If you are waiting for the TA, make a new copy of your program and begin using it in the next section.
Your preview ends here
Eager to read complete document? Join bartleby learn and gain access to the full version
  • Access to all documents
  • Unlimited textbook solutions
  • 24/7 expert homework help
38 Exercise 4: Wing as projectile on the earth (air resistance, lift, and 𝒈 ) In your remaining time, assume the projectile has the shape of a wing so that it experiences lift as it moves through the air. The lift force is also velocity dependent and its magnitude varies with the speed squared. You can model the magnitude of the lift acceleration as |? ?𝑖?? |= ?|𝑣 ⃗(𝑡)| 2 , where ? is a constant (what are the units of ? ?). Its direction is perpendicular to the motion (hopefully, the ? component is positive). One way to create a vector perpendicular to another is to use the cross-product. For example, consider the case where the wing is moving in the positive ? - direction and it’s oriented/shaped so that the lift will make the wing ri se in the positive ? - direction (like an airplane taking off). We can write this 𝑣 as |𝑣 ⃗|? ̂ and if we take the cross product of this with a vector in the negative ? direction, then (|𝑣 ⃗|? ̂) × (– ?̂) = |𝑣 ⃗|? ̂ , where ? ̂, ? ̂, and are unit vectors in the respective directions. If the cross product is unfamiliar, refresh your memory by looking at your textbook or googling it. Fortunately, Vpython has a built-in cross product function cross(ball.v, vector(0,0,-1)) will give you a vector perpendicular to ball.v with the same magnitude as ball.v, where ball.v is the velocity of the ball. So the lift acceleration vector can be written as c*mag(ball.v)* cross(ball.v, vector(0,0,-1)) . Add this code for the lift acceleration vector to the line of code where you define your net acceleration. Find the angle that gives the maximum range for a value of ? = 0.015 m -1 and for ? = 0.15 s −1 . Angle for maximum range=_________________ ; Range=____________________________ Change the size of the constant c to see its effect on the motion. Do some of the trajectories remind you of throwing a paper airplane? Describe or sketch them. Checkpoint: Run your program with the angle that produces the maximum range for ? = 0.015 m -1 , ? = 0.15 s −1 and show it to your TA. Be able to explain your results (why the trajectory looks like it does) and your observations regarding what happens as you increase c to your TA. If you are waiting for the TA, make sure that each lab partner has a copy of all the programs you have created. You will need them for next week’s lab and future labs! There is no Lab Report for this lab.
Your preview ends here
Eager to read complete document? Join bartleby learn and gain access to the full version
  • Access to all documents
  • Unlimited textbook solutions
  • 24/7 expert homework help
39 Lab 5. Modeling Orbital Motion In this lab, you will build a computational model of a spaceship traveling in the vicinity of a planet and later also an artificial moon. The forces on the spaceship will depend upon its position with respect to the position of the planet and later also with respect to the position of the moon. Newton’s law of gravity gives the attractive force bet ween two objects with masses ? and ? as ? = −? ?? |?⃗| 2 where ?⃗ is the relative position vector between the two objects and is the associated unit vector of ?⃗ and ? is the Gravitational Constant. The magnitude of the force varies with the masses ( ?, ? ) and the distance between them |?⃗| . The direction of the force is determined by the unit vector, ?̂ = ?⃗ /| ?⃗ |. Recall last week’s code, unitr = r/mag(r) # this is the unit vector of r Warmup: A spaceship orbiting a planet Make a new program giving it a new name like orbit . Start it with the same two lines as last lab, from visual import * scene = display(width=600, height=600, center=vector(0,0,0), background=color.white) You are welcome to copy and paste these lines from last week’s progra m and it would be helpful to refer to that program as you build a new program to fit the new problem described below. Represent the planet as a sphere at the origin and in your code give it the mass of 3.0 × 10 24 kg and give it a radius of ? 𝑃 = 3.0 × 10 6 m. Represent the spaceship as a sphere and give it a mass of 10,000 kg. Choose a radius that makes it small but visible (e.g. ? 𝑃 /10). Place the ship at a distance of 4 ? 𝑃 from the origin on the positive ? axis. Some suggestions for how to do this are given below. In your constants section, define the gravitational constant as G and a constant for the radius of the planets and spaceship. Something like, G=6.67e-11 # Gravitational Constant Prad=3.e6 # Radius of planet Define spheres to represent the planet and the ship (note characteristics of an object like ship.m must come after the definition of the object). planet=sphere(pos=origin, radius=Prad, color=color.green) ship=sphere(pos=vec(4.0*Prad,0,0), radius=Prad/10., color=color.yellow, make_trail=True, trail_type="points", interval=10000)
Your preview ends here
Eager to read complete document? Join bartleby learn and gain access to the full version
  • Access to all documents
  • Unlimited textbook solutions
  • 24/7 expert homework help
40 planet.m=3.e24 # Mass of planet ship.m=10000 # Mass of ship ship.v= 2500.*vec(0,1,0) # initial velocity of ship Note that the trail for the ship is now points at every 10,000 th position of the ship (this speeds up the running of the code.). We have also defined an initial velocity for the ship of 2500 m/s in the positive ? direction. You will need to place the code for the ship’s acceleration in a while loop for time , since this acceleration will depend upon the position of the spaceship, the acceleration must be updated as the ship moves. Note we will assume the planet is so massive that its motion is not influenced by its attraction to the ship (or later to the moon), i.e. the planet remains at the origin. In the while loop, you define a relative position vector r that points from the planet to the ship ( r must be in the loop since it changes with time) and the ship acceleration, r= ship.pos-planet.pos ship.a = - G*planet.m/(mag(r))**2 * r/mag(r) Orbital times are long so increase the time duration of your while loop to at least 50,000 seconds. You will also want to increase the rate of execution so you don’t have to observe this in real time. For example, the beginning of your while loop could look like, while t<50000: #loop to do the following indented code for 50,000 s rate(1000000) #executes at 1,000,000 Hz Use a time step, dt=0.2s, as your initial choice (you will change this later). Once you have finished this code, run the program. The spaceship should make an elliptical orbit passing through the planet. The code doesn’t know the spaceship crashed. Try changing the initial velocity and see how the orbit changes. Also try increasing the time step and see what happens. As a self-check, you should analytically calculate the velocity needed for a circular orbit by setting the magnitude of the centripetal force equal to the magnitude of the gravitational force, ?𝑣 2 |?⃗| = ? ?? |?⃗| 2 then check that this velocity creates a circular orbit in your program. Exercise 1: Place a moon in circular orbit about a planet In the warmup, you placed a spaceship in various orbits about a planet. Now add a moon to your program. It should be 1/60 th the mass of the planet. It should be in a circular orbit with a radius of 25 planet radii from the center of the planet. You can use the analytical solution for
Your preview ends here
Eager to read complete document? Join bartleby learn and gain access to the full version
  • Access to all documents
  • Unlimited textbook solutions
  • 24/7 expert homework help
41 the velocity for a circular orbit to determine the initial velocity of the moon (see above). Solve for the initial velocity of the moon and place it in the program. Make all the necessary additions and run the program and make sure that the moon moves as expected. Exercise 2: Turning on the interaction between the spaceship and the moon The moon also exerts a force on the spaceship. Add the code necessary to model the acceleration due to the moon’s force. Note the moon is not at a fixed position and you will have to calculate another relative position vector between the moon and the ship (call it rmoon ). Give the ship an initial velocity of 3000 m/s in the positive ? direction. Run the program and make sure that the spaceship still moves as it did in the warmup. You will have to use your mouse controls to zoom in on the spaceship. Increase the mass of the moon until there is an observable change to the orbit of the ship over the 50,000s duration of the evolution. ? ???? = ______________________ (for observable change) Checkpoint: Show your code and results to your TA. If you need to wait for the TA, make a copy of your code (in this new code make sure the mass of the moon is reset to 1/60 th of the planet). Start thinking about how you would add code that would give thrust to your ship. Exercise 3: Giving the spaceship thrust and launching it to the moon Set the initial velocity of the spaceship so that it is in circular orbit around the planet. Assume the ship’s rocket engine has a thrust that can produce a constant acceleration of 2.5 m/s 2 for a set period of time. That is, you turn it on for so many seconds and while it is on the ship has an additional acceleration in the direction of the velocity of the ship, i.e. 𝑣̂ , with a magnitude of 2.5 m/s 2 . Write code to implement this acceleration. You choose when to turn on and off the engine. The if statement in Vpython is helpful for this. Any indented lines of code after the if statement will be executed if the logical test is passed (otherwise these indented lines will be skipped). For example, you can define a variable on which is ordinarily zero, but is 1.0 between certain times, on=0.0 # engine switch is usually in off position if 1000< t <2000: #set the on variable to 1.0 between 1000s and 2000s on=1.0 # engine switch is now in on position So you can add an engine acceleration vector to your ship.a code. If this engine acceleration vector is multiplied by the variable on , then for most of the while loop the engine will be off, but it will be turned on between the times of 1000s and 2000s.
Your preview ends here
Eager to read complete document? Join bartleby learn and gain access to the full version
  • Access to all documents
  • Unlimited textbook solutions
  • 24/7 expert homework help
42 You goal is to turn on the engine’s acceleration is such a way that the ship reaches the moon. A fixed value of the time step, dt , can make the model take a very long time to run. When the spaceship is moving very slowly a longer dt can be used. For example, after entering the while loop define dt in terms of the velocity of the ship, dt=100.0/mag(ship.v) Make sure that your choice of constant in dt still gives the same evolution if it is decreased by a factor of 2, e.g. decrease the 100.0 above to 50.0 and the evolution should be unchanged. Run the program. Vary the duration of the engine’s acceleration and when it occurs to put the ship on a trajectory to hit the moon. Crashing is okay … for this model. Checkpoint: Show your code and results to your TA. If you have to wait, start reading the next section and plan the problem you would like to address. Exercise 4: Different orbital trajectory problems Now that you have your program, turn it to solving a problem of your choice. Some examples are: 1) putting the ship into orbit around the moon; 2) Use the gravitational interaction with the moon to slingshot the ship to a different trajectory at a higher speed; 3) Loop around the moon and return to an orbit around the planet; 4) Starting from the ship in orbit around the moon use your rocket engine to put the ship into an orbit around the planet; 5) Trying a planet with two moons. You are also free to explore your own problem. Checkpoint: Explain your problem and the results to your TA. Show what you had to do in your code to implement this. If you are waiting for the TA, make sure that each lab partner has a copy of all the programs you have created. You will need them for future labs! There is no Lab Report for this lab.
Your preview ends here
Eager to read complete document? Join bartleby learn and gain access to the full version
  • Access to all documents
  • Unlimited textbook solutions
  • 24/7 expert homework help
43 Lab 6. Experiment on work, energy, and momentum In this activity, you will find the energy and momentum of various systems and determine how your results align with the laws for conservation of energy and for conservation of momentum. Warmup (Kinematics and Work-Kinetic Energy Theorem): You will observe the fall of a golf ball and its first and second bounces off the floor. You will need a golf ball, a stopwatch (probably your phone), and a ruler or meter stick. The experiment, as written below, requires at least 2 people: One person holds the ball high enough that the golf will bounce to a little less than 1 meter from the floor. This will ease the measurement of maximum height. This person drops the golf ball and also tries to mark how high above the floor, y max, the ball bounces on the first bounce. Another person starts the stopwatch when the ball first strikes the floor and stops it when the ball strikes the floor again. This is twice the time required to go from the floor to the maximum height after the bounce (in the approximation that we neglect air resistance). Repeat this 5 times and record the floor to maximum height times and the maximum heights that occurred in an Excel spreadsheet (any trial in which a real error or mistake occurs sh ould be discarded and the trial repeated, e.g. stopwatch didn’t stop when pressed). Calculate the average time and average height and also the standard deviation for both of these measurements. Summarize your results in a table (example below). Averages t floortomax = y max = Standard deviation ( t floortomax )= ( y max )= Use kinematic equations for constant acceleration (i.e. 𝑣 ? = 𝑣 𝑖 − ?𝑡 , … ) with your experimental observations to determine the average initial speed of the ball just after it bounces off the floor. Use the average values t floortomax and y max . 𝑖?𝑖𝑡𝑖?? 𝑣 = ___________ m/s (from kinematic equations) Use error propagation to calculate the uncertainty in this initial velocity found from a time measurement. Use the standard deviation, ( t floortomax ), for the uncertainty of your time measurement. Refer to Lab 0 (pg. 7) for formula for error propagation. ?(𝑖?𝑖𝑡𝑖?? 𝑣) = ___________ m/s (from Newton’s second law)
Your preview ends here
Eager to read complete document? Join bartleby learn and gain access to the full version
  • Access to all documents
  • Unlimited textbook solutions
  • 24/7 expert homework help
44 Use conservation of energy (kinetic energy at floor equals potential energy at maximum height) to calculate the initial speed of the ball just after it bounces off the floor. Use the average value y max . initial v = ___________ m/s (from conservation of energy ) Use error propagation to calculate the uncertainty in this initial velocity found from a height measurement. Use the standard deviation, ( y max ), for the uncertainty of your height measurement. ?(𝑖?𝑖𝑡𝑖?? 𝑣) = ___________ m/s (from conservation of energy ) Checkpoint: Compare these two ways of measuring the initial velocity. Which do you believe is more accurate? Be prepared to explain why. Was momentum conserved for the golf ball in its collision with the floor? Was energy conserved in the collision? Again, be prepared to explain why. Is there some system where momentum and/or energy is conserved? If you are waiting for the TA, proceed to the next section and begin to design your experiment. Conservation of Momentum and Conservation of Energy Design an experiment where an object interacts with another object (e.g., golf ball colliding with another golf ball) and you measure the velocity of the objects before and after the interaction. No net external force should act upon these objects . You will be trying to determine the answers to the following key questions in your experiment: Is momentum conserved in the system of the two objects to within the uncertainty of your measurements? Is energy conserved in the system of the two objects to within the uncertainty of your measurements? You can use any of the available measuring devices you explored in Lab 2. Design the experiment and discuss it with your TA. There is also a variety of other apparatus available: air tables (pucks or light objects can move with little friction); tracks to guide rolling balls; pendulum accessories for the rotary sensor. The experiment must be safe for both you and your sensors. Have your TA okay your design for safety . The experiment and analysis should: Measure the velocities of the objects immediately before and after the interaction. From your measurements, derive the momenta and energies before and after the interaction. Calculate, from the uncertainties in your measurements, the uncertainties in your momenta and your energies. Determine if momentum and/or energy is conserved in this interaction.
Your preview ends here
Eager to read complete document? Join bartleby learn and gain access to the full version
  • Access to all documents
  • Unlimited textbook solutions
  • 24/7 expert homework help
45 The abbreviated lab report due for this experiment should have: an Abstract (200 words or less) summarizing your experiment. Its first sentence should describe what the experiment is designed to do. The next two sentences should highlight key points in the design of your experiment. The last sentence or two should give the most significant experimental result(s). For example, it should indicate the degree to which you were able to determine if momentum and/or energy was conserved in the experiment. That is, how does any change in momentum compare with the experimental uncertainty in the momentum? Similarly, make the comparison for the change in energy and the uncertainty of the energy. A procedure that lays out the experimental design in detail. A diagram or sketch of the apparatus is essential. A conclusion that discusses your experimental results. Explain how you determined if energy and/or momentum were conserved. Discuss any details of your experimental design that may have affected your ability to make these conservation determinations. Print out your lab report and exchange it with another group and have them check that your report contains all the required elements. If you need to wait for another group to finish look at nex t week’s lab and begin your preparations (coding). Make any necessary revisions and turn it in to your TA. The rubric your TA will use to grade this report is below. Rubric for Report for Lab 6. If they have all the required elements in their report, their grade should be 100. The abbreviated lab report due for this experiment should have: an Abstract (200 words or less) summarizing your experiment. (if over 200 words minus 5 pts.) o First sentence should describe what the experiment is designed to do. (10 pts) o The next sentences should highlight key points in the design of experiment. (10 pts) o The last sentences should give the most significant experimental result(s). For example, “momentum was found to be conserved within the measurement precisio n of the apparatus” (10 pts) A Procedure section giving o A detailed description of the experiment (15 pts) o A sketch of the apparatus. (20 pts) A Conclusions section with o Stating whether momentum was found to be conserved and why. (10 pts)
Your preview ends here
Eager to read complete document? Join bartleby learn and gain access to the full version
  • Access to all documents
  • Unlimited textbook solutions
  • 24/7 expert homework help
46 o Stating whether energy was found to be conserved and why. (10 pts) o Comment on any improvements that could be made in your experimental design to increase the precision of your results. (15 pts.)
Your preview ends here
Eager to read complete document? Join bartleby learn and gain access to the full version
  • Access to all documents
  • Unlimited textbook solutions
  • 24/7 expert homework help
47 Lab 7. Computational Model of Collisions In this lab, you will build a computational model of the collision of two charged objects. The Coulomb force is conservative (conserves energy) and has the form, ? = ? ?? |?⃗| 2 Where ? is a constant that is approximately 8 .987552 × 1 0 9 N · m 2 · C −2 , ?⃗ is the relative position vector between the two objects and is the associated unit vector of ?⃗. There are strong similarities between the Coulomb force and the gravitational force. There is one major difference, when the charges ? and ? have the same sign the force is repulsive (gravity is always an attractive force). In this activity, the two charges will have the same sign so the charges will never come in contact, but instead bounce away from each other due to the Coulomb force. Warmup: 1-D Collision As you write this new program your spaceship program from Lab 5 will be a helpful guide. However, there are a few significant differences: You will need to replace the gravitational law with the Coulomb law. Give the objects charges with the same sign (e.g. both positive) so the force is repulsive (not attractive). The masses of the two objects will be similar. In the previous program the mass of the planet was so much larger that we assumed it didn’t move. You will need to have code for the mo vement of both objects. Fortunately, Newton’s 3 rd law makes the magnitude of the force the same on both objects and makes the directions the opposite (one Force vector is the negative of the other). In your new program, let’s define your constants, vari ables, and objects in the following way for ease of debugging. Your two objects are ball1 and ball2 and each is represented by a sphere. Define these variables as follows, k=8.987552e9 # Coulomb Constant ball1rad=0.1 # Radius of ball one ball2rad=0.2 # Radius of ball two ball1= sphere(pos=vec(0,0,0), radius=ball1rad, color=color.green, make_trail=True, trail_type="points", interval=100) # position is the origin ball2= sphere(pos=vec(-5,0,0), radius=ball2rad, color=color.red, make_trail=True, trail_type="points", interval=100) # position is -5 meters on the x-axis ball1.m=1.0 # Mass of ball 1 ball2.m=3.0 # Mass of ball 2 ball1.q=5e-6 # Charge of ball 1 in Coulombs ball2.q=5e-6 # Charge of ball 2 in Coulombs
Your preview ends here
Eager to read complete document? Join bartleby learn and gain access to the full version
  • Access to all documents
  • Unlimited textbook solutions
  • 24/7 expert homework help
48 Let ball1 be initially stationary, but give ball2 an initial velocity of 1 m/s in the positive ? direction, e.g. ball2.v=vec(1,0,0) Modify your while loop for time with a reasonable duration, time step, and rate of execution. It may be easiest to check this without the accelerations, i.e. try it first with ball1.a and ball2.a equal to vec(0,0,0) and see if the motion makes sense. Implement the accelerations of each ball and update their motions in the while loop. Run the program and decide if you have correctly modeled the behavior of two colliding balls. As always check that you are using a sufficiently small time step. The behavior should not change if the time step is reduced by half. Try changing the initial velocities and/or the masses of the objects and see what happens. In particular, check what happens when the two objects have the same mass. Checkpoint: Show your code and results to your TA, in particular show the equal mass result. Make sure you can explain the motion of the particles in terms of conservation of momentum and energy. If the force between the particles was not the Coulomb force, but instead some other force that varies with one over the distance cubed, ? = ? ?? |?⃗| 3 What would you need to change in your program to model the collision with this new force? Discuss this with your group and explain it to your TA. If you need to wait for the TA, make a copy of your code and begin modifying it to address the problem of the next section, collisions in two dimensions. Collisions in 2 and 3 dimensions Since your code is written in terms of vectors, it should handle collisions in higher dimensions. Change the initial position of ball2 to slightly above the ? -axis but still give it an initial velocity parallel to the ? -axis, e.g. ball2= sphere(pos=vec(-5,0.1,0), radius=ball2rad, color=color.red) ball2.v=vec(1,0,0) Run the program and record the initial and final momenta for the system of ball1 and ball2 . Similarly, record the initial and final kinetic energies. There is a table on the following page to record these.
Your preview ends here
Eager to read complete document? Join bartleby learn and gain access to the full version
  • Access to all documents
  • Unlimited textbook solutions
  • 24/7 expert homework help
49 The values of the velocities at the start of the program should be used to calculate the initial initial momenta and energies. When do you record the final velocities or when is the collision over? We started with ball1 and ball2 separated by a large distance (5 meters) so there would be little repulsion between the particles. The final velocities and momenta should be taken when ball1 and ball2 are again separated by 5 meters. Adjust the length of the time loop to make this so. You will need to print out the four final momentum components (make sure you change the position where each prints so you can read and identify them). From these momenta, you can also calculate the final energies, e.g. ? = |? ⃗| 2 /2? . There is example code below the table demonstrating how to print out one of the components. ? 1?,𝑖 = _____________ ? 1?,𝑖 = _____________ ? 1,𝑖 = _____________ ? 2?,𝑖 = _____________ ? 2?,𝑖 = _____________ ? 2,𝑖 = _____________ ? 1?,? = _____________ ? 1?,? = _____________ ? 1,? = _____________ ? 2?,? = _____________ ? 2?,? = _____________ ? 2,? = _____________ p1xf=ball1.m*ball1.v.x p1xfstring= str(p1xf) ## converts final ball1 momentum to string for printing label( pos= vec(0,0.5,0), text= p1xfstring ) ## prints above origin You will need to print out the four final momentum components (make sure you change the position where each prints so you can read and identify them). How could you change the initial conditions to make the collision a 3D problem? Try it out and see if you are correct. Checkpoint: Show your code and results to your TA. Be sure you can explain how your results show that energy and momentum are conserved. If you need to wait for your TA, read ahead about the Rutherford experiment and consider how you could model this with your code. Modeling the Rutherford Experiment The Rutherford experiment is one of the crucial experiments of modern atomic physics which revealed the basic structure of the atom. It found that the atom consists of a small but heavy, positively charged nucleus surrounded by light negatively charged electron(s). This cloud of electrons is large and determines the volume of an atom or in a solid the spacing between nuclei. The radius of the nucleus is on the order of 10 −15 m while the diameter of the atom is roughly 10 −10 m, still the nucleus contains nearly all the mass of the atom. So a solid can be thought of as containing very localized heavy masses and largely empty space.
Your preview ends here
Eager to read complete document? Join bartleby learn and gain access to the full version
  • Access to all documents
  • Unlimited textbook solutions
  • 24/7 expert homework help
50 In the experiment, a beam of alpha particles was fired at a gold foil and the directions the alpha particles were scattered were measured. For this model, you can ignore the electrons of the gold foil by assuming they are spread uniformly in the foil so will not affect the path of the alpha particles. Alpha particles are positively charged and are ~2% of mass of gold nucleus. The radius of an alpha particle is smaller than that of a gold nucleus by about a factor of seven. In this program, we will use the same radius for both for ease of visibility (you will find that the actual radius is not significant in this problem). The characteristics of a gold nucleus and an alpha particle can be coded as, goldrad=7.e-15 # Radius of gold nucleus alpharad=7.e-15 # Radius of alpha particle gold= sphere(pos=vec(0,0,0), radius=goldrad, color=color.green, make_trail=True, trail_type="points",interval=100) alpha= sphere(pos=vec(-200*goldrad,0,0), radius=alpharad, color=color.red, make_trail=True, trail_type="points",interval=100) gold.m=197.*1.67e-27 # Mass of gold nucleus alpha.m=4.*1.67e-27 # Mass of alpha gold.q=79.*1.6e-19 # Charge of gold alpha.q=2.0*1.6e-19 # Charge of alpha You may assume the nucleus remains fixed at the origin and an alpha particle travels towards it on a path parallel to the ? -axis at high speeds (but non-relativistic speeds). alpha.v=vec(1e7,0,0) # velocity of alpha particle A beam of alpha particles is made up of alpha particles with the same initial velocities, but different initial positions, i.e. different displacements from the ? -axis. In your code, You will need to modify your time steps and duration of while loop to handle the large changes in orders of magnitude of sizes, speeds, charges, etc. As a first estimate for dt , try dt = 1e-23 s. Why is it so small? How could you guess it? How long does it take the alpha particle to travel the radius of the gold atom? Model a beam of alpha particles by taking 101 alpha particles each with a different ? -offset from the ? -axis. Start with alpha.pos=vec(-200*goldrad, -200*goldrad,0) and the increase by 4 gold radii for the next particle and so on. Adding another while loop around the while time loop to change these initial conditions would be a good way to accomplish this. Below is a portion of code that nests a two particle collision loop (like you wrote above) inside another while loop that changes the initial conditions of the alpha particle and then repeats the collision code 101 times. j=0
Your preview ends here
Eager to read complete document? Join bartleby learn and gain access to the full version
  • Access to all documents
  • Unlimited textbook solutions
  • 24/7 expert homework help
51 while j < 101: #outer loop which sets initial positions, velocities of the alpha particles j=j+1 # counter for the number of the alpha particle, e.g. 1 st , 2 nd , 3 rd , alpha.pos=vec(-200*goldrad,-(200-4*j)*goldrad,0) # different for each alpha r=alpha.pos-gold.pos # initial relative position vector alpha.v=vec(1e7,0,0) # resets it to same initial value for each alpha alpha.a= k*gold.q*alpha.q/alpha.m/mag(r)/mag(r)*r/mag(r) #initial alpha acceleration t=0 # start time at zero seconds for each alpha particle while t<29e-20: #inner time loop for collision of one alpha particle with gold atom rate(1000000) alpha.pos=alpha.pos+alpha.v*dt #update the alpha position after dt has passed alpha.v=alpha.v+alpha.a*dt #update the alpha velocity after dt has passed r=alpha.pos-gold.pos # new relative position vector alpha.a= k*gold.q*alpha.q/alpha.m/mag(r)/mag(r)*r/mag(r) #new alpha acceleration t=t+dt #increase time by dt Make sure each alpha particle leaves a trail so you can see the distribution of trajectories. Checkpoint: Show your code and results to your TA. Is it clear why this set of trajectories would lead you to believe the nucleus was very small? What would happen if the charge of the nucleus was spread uniformly over a very large area? What would you have to change in your code if you broke up the nucleus into a thousand smaller particles (each 1/1000 of the charge and 1/1000 of the mass) spread over a large area? What assumptions would change? Discuss this new model of the atom in your group and then present it to your TA.
Your preview ends here
Eager to read complete document? Join bartleby learn and gain access to the full version
  • Access to all documents
  • Unlimited textbook solutions
  • 24/7 expert homework help
52 Lab 8. Experiment on Rotational Dynamics and Conservation of Angular Momentum In this activity, you will measure the rotational motion of a system. Rotational dynamics is quite like the dynamics of linear motion. For linear motion the linear momentum is conserved. For rotational motion angular momentum is conserved. As a warmup to the analyses needed for this experiment, answer the following questions, 1) There is an object bolted to the axle of a pulley. The axle is vertical. A string is wrapped around the pulley (at a radius, ? ) and leaves it horizontally going over a second pulley whose axle is horizontal. A mass ? is attached to the other end of the string and allowed to fall (see picture on following page). The tension in the string exerts a torque, 𝜏 , on the object-pulley system. Sketch the forces acting on the falling mass. Write an equation that describes the acceleration, ? , of the mass ? and solve for the tension, 𝑇 , in the string. Assume both ? and ? (gravitational acceleration) are positive. 2) Assume you measure the angular acceleration, 𝛼 , of the object-pulley system due to the falling mass in the question described above. Using Newton's second law for rotary motion ( 𝜏 = 𝐼𝛼) and the result from above, calculate the moment of inertia, 𝐼 , for the object-pulley system. Write an equation for 𝐼 in terms of ? , ? , ? , and 𝛼 . Recall that 𝜏 = 𝑇? for this setup.
Your preview ends here
Eager to read complete document? Join bartleby learn and gain access to the full version
  • Access to all documents
  • Unlimited textbook solutions
  • 24/7 expert homework help
53 Warmup (Rotational Dynamics of Spinning Disk): You will angularly accelerate the pulley on the Rotary Motion sensor with a string is wrapped around the pulley. The other end of the string is attached to a falling known mass, ? . See picture below. There should not yet be an aluminum disk attached to the rotary motion sensor (only the 3 step pulley assembly). The string is wrapped around the largest pulley in the pulley assembly. Recall from Lab 1 o Three-step Pulley: 10, 29 and 48 mm diameters o Resolution: ±0.09° or 0.00157 radian or 4,000 divisions/revolution Record the rotational motion ( 𝜃(𝑡) , ?(𝑡) , and 𝛼(𝑡) ) using the Pasco sensor and interface. Use the angular acceleration data to calculate the moment of inertia of the rotary motion sensor system. Use the problems you solved above to guide your analysis. That is, moment of inertia should be calculated from measurements of the mass of the falling mass, radius of the pulley (where string is wrapped; note this is given in specifications of rotary sensor), the gravitational acceleration ( ?) , and the average angular acceleration 𝛼 ̅ . Use the standard deviation from this average value 𝛼 ̅ to estimate the uncertainty in 𝛼 . Record the moment of inertia below for rotary motion sensor system. Calculate and include the uncertainty of the moment of inertia. 𝐼 ?????? = _______________________________ Now attach the aluminum disk with the thumb screw to make a disk rotary motion sensor system. Repeat your measurements and perform similar analyses for this system. Record
Your preview ends here
Eager to read complete document? Join bartleby learn and gain access to the full version
  • Access to all documents
  • Unlimited textbook solutions
  • 24/7 expert homework help
54 the moment of inertia below for disk-rotary motion sensor system. Calculate and include the uncertainty of the moment of inertia. 𝐼 ?𝑖??,?????? = _______________________________ Checkpoint: Was angular momentum conserved for the disk-rotary motion sensor system in this experiment? How about if you include the falling mass in the system? If not, is there a system where angular momentum is conserved? If you are waiting for the TA, proceed to the next section and begin to design your experiment. Conservation of Angular Momentum During a Brief Interaction Design an experiment where a rotating object interacts with another object (rotating or not rotating) and you measure the angular velocity of the objects before and after the interaction. You would like the angular momentum to be conserved in the system of the two objects to within the uncertainty of your measurements. You will also be measuring whether energy is conserved in the system of the two objects to within the uncertainty of your measurements. Let’s assume that one of the objects is the disk-rotary sensor setup from the warm-up. The following are some questions to help you design your experiment (discuss them with other groups and TA if you like): 1) If you want angular momentum to be conserved should there be an external torque on the system? How could you spin the disk-rotary sensor platform to satisfy this? 2) Is there a way to connect two disk-rotary sensors by a string to enable the desired interaction? 3) Could you drop something onto the disk-rotary sensor and have the desired interaction? 4) If a cart is attached by a string to a stationary disk-rotary sensor and it rolls past the rotary sensor causing the string to suddenly go taut, could the cart-disk-rotary sensor system display the desired angular momentum conserving interaction? How would you measure the linear velocity of the cart? The experiment must be safe for both you and your sensors. Have your TA okay your design before performing the experiment . The experiment and analysis should: Measure the angular velocities and/or linear velocities of the objects immediately before and after the interaction. From your measurements, derive the angular momenta and energies before and after the interaction. Calculate, from the uncertainties in your measurements, the uncertainties in your angular momenta and your energies.
Your preview ends here
Eager to read complete document? Join bartleby learn and gain access to the full version
  • Access to all documents
  • Unlimited textbook solutions
  • 24/7 expert homework help
55 Determine if angular momentum and/or energy is conserved in this interaction. The lab report due for this experiment should have: an abstract (200 words or less) summarizing your experiment. The abstract should highlight key points in the design of your experiment and in the experimental results. In particular, it should indicate the degree to which you were able to determine if angular momentum and/or energy was conserved in the experiment. That is, how does any change in angular momentum compare with the experimental uncertainty in the angular momentum? Similarly, make the comparison for the change in energy and the uncertainty of the energy. An introduction where the significance of the conservation of angular momentum is described. It should contain at least one peer-reviewed journal article where this conservation of angular momentum plays a key role. A sentence or two summarizing this article should be included. A procedure that lays out the experimental design in detail. A diagram or sketch of the apparatus is essential. A conclusion that discusses what you have found from the experimental results. Explain how you determined if energy and/or angular momentum were conserved. Discuss any details of your experimental design that may have affected your ability to make these conservation determinations. Print out your lab report and exchange it with another group and have them check that your report contains all the required elements. If you need to wait for another group to finish look at next week’s lab and begin your preparations for this experiment on harmonic motion. Make any necessary revisions and turn it in to your TA. The rubric your TA will use to grade this report is below. Rubric (If they have all the required elements in their report, their grade should be 100.) The abbreviated lab report due for this experiment should have: an Abstract (200 words or less) summarizing your experiment. (if over 200 words minus 5) o First sentence should describe what the experiment is designed to do. (10 pts) o The next sentences should highlight key points in the design of experiment. (10 pts) o The last sentences should give the most significant experimental result(s). For example, “angular momentum was found to be conserved within the measurement precision of the apparatus” (10 pts) An Introduction o The significance of conservation of angular momentum should be described. (5 pts)
Your preview ends here
Eager to read complete document? Join bartleby learn and gain access to the full version
  • Access to all documents
  • Unlimited textbook solutions
  • 24/7 expert homework help
56 o It should contain at least one peer-reviewed journal article where this conservation of angular momentum plays a key role. A brief summary of this article should be included. (15pts) A Procedure section giving o A detailed description of the experiment (10 pts) o A sketch of the apparatus. (10 pts) A Conclusions section with o Stating whether angular momentum was found to be conserved and why. (10 pts) o Stating whether energy was found to be conserved and why. (10 pts) o Discuss any details of your experimental design that may have affected your ability to make these conservation determinations. Comment on any improvements that could be made in your experimental design to increase the precision of your results. (10 pts.) Lab 9. Experiment on harmonic motion In this lab, you will explore systems that display harmonic motion (sinusoidal oscillations). Warmup (harmonic and anharmonic motion of a pendulum): We’ve mounted a Rotary Motion sensor on a lab rod so that it’s axle is horizontal. One end of the a rod is attached to the axle at one end and a mass to the other end (as shown in diagram below). For small, initial angular displacements (~ 10 degrees), record the rotational motion ( 𝜃(𝑡) , ?(𝑡) , and 𝛼(𝑡) ) using the Pasco sensor and interface. . o Where do the maxima of ? and 𝛼 occur? o What is the period and frequency of the oscillation (use several cycles of data)? 𝑇 = ____________________ ? = _______________________
Your preview ends here
Eager to read complete document? Join bartleby learn and gain access to the full version
  • Access to all documents
  • Unlimited textbook solutions
  • 24/7 expert homework help
57 The period is independent of the initial angle as long as the initial angles are small. If you can’t recall why, check your textbook. How much does the period change for large initial angles? Repeat the above experiment and analyses for large, initial angular displacements (90 degrees or more). o How do the graphs of the rotational motion ( 𝜃(𝑡) , ?(𝑡) , and 𝛼(𝑡) ) compare? Is the motion harmonic? o What is the period and frequency of the oscillation (use several cycles of data)? 𝑇 = ____________________ ? = _______________________ o Checkpoint: Did the period change? By what percent? Could you tell by just looking at the two sets that the 𝜽(?) data was different? Why was it easiest to see a difference in the 𝜶(?) data? See if you can argue why 𝜶(?) looks nearly constant at the peaks. Be prepared to show your graphs and explain your results and conclusions to the TA. If you are waiting for the TA, proceed to the next section and begin to design your experiment. Explore Damped Oscillators Design an experiment to add damping to the pendulum system used in the warm-up. The amount of damping should be able to be varied so that you can explore a range of damping from underdamped to critical damping. Refer to your textbook for definitions of underdamped and critical damping, if you do not have your text available you may consult the website below for a review. http://hyperphysics.phy-astr.gsu.edu/hbase/oscda.html The experiment must be safe for both you and your sensors. Have your TA okay your design before performing the experiment. The experiment and analysis should: Measure the position of the oscillator as a function of time for different levels of damping. Adjust the damping so that you can display the following 2 cases: o Weakly damped: Envelope of oscillations decays exponentially showing a few cycles before coming to a stop. o Strongly damped: No oscillations, just decay
Your preview ends here
Eager to read complete document? Join bartleby learn and gain access to the full version
  • Access to all documents
  • Unlimited textbook solutions
  • 24/7 expert homework help
58 Make two graphs of the motion ( 𝜽 𝒗?. ?) of the damped pendulum displaying the two of the above cases, i.e. one graph for each case. Try to fit your weakly damped experimental data to a function of the form 𝐴? −?? cos ?𝑡 . Create another column in Excel with this function and graph it along with your data. Adjust the parameters 𝐴, ? , ??? ? to make them match. Use your data to help you make good guesses for the parameters, e.g. use the time between peaks to give you a guess for the period and from that ? . Using absolute addressing for the parameters is helpful, e.g. if a parameter is in the cell G2, reference it in the formula as $G$2 then the graph will update immediately when you change this parameter. What is the frequency of the oscillation? What is the decay constant? ? = _____________________________ ; ? = _______________________________ Try to fit your strongly damped experimental data to a decaying exponential. What is the decay constant? Excel has a trendline function that you can use here. ? = _______________________________ The lab report due for this experiment should have: an abstract (200 words or less) summarizing your experiment. The abstract should highlight key points in the design of your experiment and in the experimental results. In particular, it should describe the means of damping. It should give the experimental results for the weakly damped case and the strongly damped case. An introduction where the significance of harmonic motion is described. It should contain at least one peer-reviewed journal article where harmonic motion or oscillators play a key role in modeling a physical system. A brief summary of this article should be included to give a purpose and background to your experiment. A procedure that lays out the experimental design in detail. A diagram or sketch of the apparatus is essential. A results section that includes the two graphs (weakly damped and strongly damped cases). The fit curve for each case should also be plotted on the associated graph. Describe what each graph is displaying and include the fit parameters for the curves in your description. A conclusion that discusses what you have found from the experimental results. Discuss any modifications of your experimental design that you might make to improve your ability to adjust the damping.
Your preview ends here
Eager to read complete document? Join bartleby learn and gain access to the full version
  • Access to all documents
  • Unlimited textbook solutions
  • 24/7 expert homework help
59 Print out your lab report and exchange it with another group and have them check that your report contains all the required elements. If you need to wait for another group to finish look at next week’s lab and begin your preparations for modeling waves. Make any necessary revisions and turn it in to your TA. The rubric your TA will use to grade this report is on the following page.
Your preview ends here
Eager to read complete document? Join bartleby learn and gain access to the full version
  • Access to all documents
  • Unlimited textbook solutions
  • 24/7 expert homework help
60 Rubric (If they have all the required elements in their report, their grade should be 100.) The abbreviated lab report due for this experiment should have: an Abstract (200 words or less) summarizing your experiment. (if over 200 words minus 5) o First sentence should describe what the experiment is designed to do. (5 pts) o The next sentences should highlight key points in the design of experiment. (5 pts) o The last sentences should give the most significant experimental result(s). For example, “in the case of weak damping, we found an oscillation frequency of 3 Hz and a decay constant of 0.1 s -1 ” (5 pts) An Introduction o The significance of using damped harmonic motion to model physical systems should be described. (5 pts) o It should contain at least one peer-reviewed journal article where a phenomena is modeled as a damped harmonic oscillator. A brief summary of this article should be included to give a purpose and background to your experiment. (10pts) A Procedure section giving o A detailed description of the experiment (10 pts) o A sketch of the apparatus. (10 pts) A Results section with properly formatted graphs of o Weakly damped oscillation, with fit function (10 pts) o Strongly damped oscillation, with fit function (10 pts) o Description of each graph including the fit parameters (10 pts) A Conclusions section with o Stating what type of damping was observed. (10 pts) o Discuss any modifications of your experimental design that you might make to improve your ability to adjust the damping. (10 pts.)
Your preview ends here
Eager to read complete document? Join bartleby learn and gain access to the full version
  • Access to all documents
  • Unlimited textbook solutions
  • 24/7 expert homework help
61 Lab 10. Computational Model of Waves In this activity, you will build a computational model of waves propagating in a medium. The medium will be represented by many spheres connected by springs with the same spring constant ? . For example, imagine a line of these spheres (one dimensional). The force on the 𝑖𝑡ℎ sphere depends upon its position and the positions of its neighbors, ? 𝑖 = −?(? 𝑖 − ? 𝑖−1 ) − ?(? 𝑖 − ? 𝑖+1 ) Convince yourself that this force is zero when the 𝑖𝑡ℎ sphere is at the same distance from each neighbor. Displacing one of these spheres will lead to a force trying to restore it to its original position. However, this force will also affect its neighbors leading to a disturbance or wave. Warmup: 1-D Longitudinal Waves In this warmup, you will be dealing with many similar objects or atoms. It is convenient to introduce the concepts of class and list . A list is an array holding objects and associated variables. A class can be thought of as a function that creates a list of objects. Below is a example program that defines a set of atoms arranged in a line and each having a random velocity. Note that you can access the variables associated with an object in terms of the position of the object in the list, e.g. the velocity of the first atom c.atoms[0].v from visual import * N = 11 # N of atoms in set spacing = 1 # distance between atoms atom_radius = 0.25*spacing # radius of atom scene.center = 0.5*(N-1)*vector(1,0,0) # center the scene on the middle atom class set: # define a class of objects to represent set of atoms; def __init__(self, N, atom_radius, spacing ): # self is the variable that is equal to the class set, N, radius, spacing self.atoms = [ ] #defines list storing the objects for x in range(0,N,1): #loop defining objects ... spheres as atoms atom = sphere(radius=atom_radius, color = color.red) # all atoms same radius and color atom.pos = vector(x,0,0)*spacing #initial position of an atom on x-axis atom.v = vec.random() # give each atom a random velocity self.atoms.append( atom ) #add each atom to list c = set(N, atom_radius, spacing) # the variable c is assigned to the set; in the class it is self dt=.005 #time step t=0 #time while t<1000: #time loop i=0 #list position for first atom rate(100) while i<N: #updating position for each atom c.atoms[i].pos=c.atoms[i].pos+c.atoms[i].v*dt i=i+1 t=t+dt
Your preview ends here
Eager to read complete document? Join bartleby learn and gain access to the full version
  • Access to all documents
  • Unlimited textbook solutions
  • 24/7 expert homework help
62 Read through the example program and make sure you understand the purpose of each line. You do not need to retype the example program there is a copy ( atomsinrow.py ) you may download on Canvas . Run it. You should see eleven atoms initially on the ? -axis that then start to move in different directions. Make a copy of this program and you are now going to add the spring forces between the atoms. Assume the ends of the chain of atoms are fixed so they do not evolve in time. That is, you do not need to calculate a force for the first or last atoms or update their velocities or positions. Set all the atom.v=vec(0,0,0) , not random. Your new program will need a few things, You will need to define the mass of each atom and spring constant. To start use a value of one for the mass m and the spring constant k . In the class set, you will need to define a variable atom.force=vec(0,0,0) The forces on the atoms change with time so inside the time loop, you will need a loop to calculate all the forces on the non-fixed atoms. Before the position update loop add another similar loop. It should look something like, i=1 while i<N-1: c.atoms[i].force = “ you insert the code for the force on the ith atom ” i=i+1 After the force loop, you will need to change the position update loop to use the forces calculated in the previous loop to first update the velocity of each nonfixed atom and then the position of each non-fixed atom. If you run the program, you should see a set of stationary atoms (all the net forces are zero). To make something happen, you can change the initial position of one of the atoms. For example, you can move the fixed first atom closer to the second atom. Insert something like the following immediately before your time loop. c.atoms[0].pos=c.atoms[0].pos+vec(0.5*spacing,0,0) This will create a disturbance that propagates down the chain. Run the program and decide if this is what you see. Remove the above statement and inside the class definition let’s arrange the initial atoms so they are initially more closely spaced on the right end of the chain (or higher density of atoms on the right end). That is, atom.pos=vector((N-1)*sin(.5*pi*x/(N-1)),0,0)*spacing
Your preview ends here
Eager to read complete document? Join bartleby learn and gain access to the full version
  • Access to all documents
  • Unlimited textbook solutions
  • 24/7 expert homework help
63 Set spacing to 1. Run it and see if the motion of the atoms makes sense to you. Checkpoint: Show your code and results to your TA. Where on the chain is the atom that is displaced the most, that is travels the largest distance back and forth (maxima)? Where on the chain are atoms that are displaced the least (minima or node)? Where does the density of atoms vary the most/least (that is the spacing between atoms varies the most/least)? What is the relation between the variations in the displacement and the variations in density? If you need to wait for the TA, make a copy of your code and begin modifying it to address the problem of the next section, transverse waves. Transverse Waves Since your code is written in terms of vectors, it should handle transverse waves with essentially no modification. Instead of creating a disturbance by displacing the atoms along the ? -axis, you will displace the atoms in the ? -direction to create the disturbance. Make a sketch below of what the atoms and springs would look like if you displaced the central atom half a spacing in the positive y direction. Which way is the net force on that central atom? In your code, you will need to reset the ? positions of the atoms back to evenly spaced in the class definition. Also, Increase the number of atoms to 101 Try displacing only the central atom by 5 spacings in the +y direction. (Which is the central atom?). What happens? Does it remind you of anything? Change the initial positions of atoms so that they will cause the atoms to oscillate as if they were a string oscillating in its fundamental mode. What function should you use to describe the initial y positions? How about for the second harmonic?
Your preview ends here
Eager to read complete document? Join bartleby learn and gain access to the full version
  • Access to all documents
  • Unlimited textbook solutions
  • 24/7 expert homework help
64 Checkpoint: Show your code and results to your TA. Be sure you can show both the fundamental and the second harmonic. If you need to wait for your TA, make a copy of your code and begin to modify it for the next activity. Modeling driven waves on a string If the ? position of the first atom oscillates sinusoidally it will drive waves on the string. The ? position of the first atom should have the form, ?(𝑡) = 𝐴 sin(2𝜋?𝑡) Where 𝐴 is an amplitude and ? is the frequency of this driving oscillation. Modify your program to make this happen. Implement the function, ?(𝑡) = 𝐴 sin(2𝜋?𝑡), in code. Where should you place this line in your program? o Recall that the position variable for the first atom in code is denoted c.atoms[0].pos o Remember that this variable is a vector! The ? position should remain unchanged, you want oscillations only in the ? direction. Find the frequencies that drive the fundamental and the second harmonic. ? ??????????? = ____________________ ; ? ?????? = ____________________ Checkpoint: Show your code and results to your TA. Be sure you can show the drive causing the string to vibrate in the fundamental and also show the case for driving the second harmonic. Why does the amplitude of the standing waves grow rapidly? There is one final activity for this course. It is a short quiz called Modeling Physics. You may discuss it with your group before you enter your answers, but you should each individually enter your answers on Canvas. Your grade on this quiz will NOT be based on these entries. Once you have completed your entries, please go over your answers with your TA. Your grade will be based on participation and the grade will be entered manually. There is no lab report for this activity.
Your preview ends here
Eager to read complete document? Join bartleby learn and gain access to the full version
  • Access to all documents
  • Unlimited textbook solutions
  • 24/7 expert homework help