Lab 3 Graphing Processes Modeled by Exponential Functions

Download as docx, pdf, or txt
Download as docx, pdf, or txt
You are on page 1of 5

Lab 3: Graphing Processes Modeled by Exponentials

Save this document to your own computer. As you work through the lab, fill in tables and
copy and paste results from MATLAB as needed.

Many physical processes can be modeled using the exponential function. Some examples
include charging or discharging a capacitor, population growth or decay, heating and cooling,
radioactive decay, and chemical reaction rates.

A. Decomposition of Hydrogen Peroxide

First order chemical reactions can be modeled using exponential functions. Hydrogen Peroxide
(H2O2) decomposes as a 1st order reaction into water and oxygen gas:

2 H 2 O2 →2 H 2 O+O2

The concentration of hydrogen peroxide decreases exponentially according to the following


equation:
−kt
C ( t )=C 0 e

C(t) = Concentration at time t (M or mols/L)


C0 = Initial Concentration (M)
k = Reaction rate (s-1 )

The reaction rate, k, depends on the temperature according to Arrhenius’ Equation (also
exponential):

k = A e−E / (RT)
a

A = Frequency Factor (s-1 )


Ea = Activation Energy (J/mol)
R = Ideal Gas Constant = 8.314 (J/(mol*K))
T = Absolute Temperature (K)

1. Assume Ea = 7.5*104 J/mol and A = 1011 s-1. Complete the following table to calculate the
reaction rate for each of the given temperatures. When completed, check your values with
your T.A. The value of k should be calculated to five places behind the decimal point – use
the command format long in MATLAB. Remember e-3t is exp(-3*t) in MATLAB.

Temperature (oF) Temperature (oC) Absolute Temperature (K) Reaction Rate, k (s-1)
50 10 283.15 0.0015
77 25 298.15 0.0072
104 40 313.15 0.0308

2. Assume the initial concentration, Co of the hydrogen peroxide is 1 M. In MATLAB, on the


same graph, plot the concentration, C, of the hydrogen peroxide vs. time for all three
temperatures. The concentration should be on the y-axis and time should be on the x-axis.

Hints:
 Remember in MATLAB, e-3t is entered as exp(-3*t)
 You will need to think about a good time range for your plot. A good rule of thumb
is to remember that e-5 is reasonably close to zero.

3. Add the following to your plot then copy and paste your plot into this document.
 Label the x-axis as time (s)
 Label the y-axis as Concentration (M)
 Add a title: Decomposition of Hydrogen Peroxide
 Add a legend that shows the three temperatures in oF
4. Use your plot and the Data Cursor Tool to complete the following table.

Temperature (oF) 50oF 77oF 104oF


Time required for concentration to drop to ½ of the initial concentration 460s 96s 22.8s
Time required for concentration to drop to ¼ of the initial concentration 923.5s 192s 45.35s
Time required for concentration to drop to 10% of the initial 1530s 310s 75s
concentration

B. RC Circuits

The figure below shows a simple RC (Resistor Capacitor) circuit. When the switch is in position
P1, the capacitor will charge at an exponential rate to a voltage equal to the battery voltage (5 V
in this circuit). When the switch is in position P2, the capacitor will discharge at an exponential
rate to 0 V. While the capacitor is discharging, it will act as a temporary battery to light up the
LED (light emitting diode). The time it takes for the capacitor to charge or discharge depends on
the size of the resistors and the capacitor, specifically on the product, RC, referred to as a time
constant.

Capacitor Charging

The equation that describes how a capacitor charges over time in an RC circuit is:
−t /τ
V c (t)=E(1−e )

Vc(t) = Capacitor voltage at time t (V)


E = Voltage of the battery or source (V)
τ is the RC time constant (s)

The time constant, τ, is simply the product of the resistance and the capacitance:

τ =R ∙C
τ = RC time constant (s)
R = Resistance in ohms (Ω)
C = Capacitance in farads (F)
1. Calculate the time constant for each of the resistances shown in the table below.
Note: To calculate τ in seconds, capacitance must be in farads (F) and resistance must be in
ohms (Ω). 1 µF = 10-6 F and 1 kΩ = 103 Ω.
Hint: In MATLAB, you can enter the value 4.7*10-6 as 4.7*10^-6 or as 4.7e-6

Resistanc Capacitance τ (s)


e
1 kΩ 4.7µF 4.7 ms
2 kΩ 4.7µF 9.4 ms
10 kΩ 4.7µF 47 ms

2. Assume the battery voltage is 5V. On the same graph in MATLAB, plot the capacitor
voltage for all three resistance values. Determine a good time range.

3. Add the following to your plot then copy and paste your plot into this document.
 Label the x-axis as time (s)
 Label the y-axis as Capacitor Voltage (V)
 Add a title: Capacitor Charging in RC Circuit
 Add a legend that shows the three resistance values
4. Using the data cursor tool, determine how long it takes the capacitor voltage to reach 63.2%
and 98.2% of its final value (final value is battery voltage). Enter results in the table below:

Resistance Capacitanc Time to Reach Time to Reach


e 63.2% of 5V 98.2% of 5V
1 kΩ 4.7µF 4.7 milliseconds 0.019 seconds
2 kΩ 4.7µF 9 milliseconds 0.038 seconds
10 kΩ 4.7µF 46 milliseconds 0.19 seconds

Questions:

1. How is the time required to reach 63.2% of the final voltage related to the time constant,
τ?
A: the time required to reach 63.2% of the final voltage is equivalent to the time
constant.

2. Approximately how many time constants does it take the capacitor to reach 98.2% of its
final voltage?
A: it takes 4.0425 seconds of time constants for the capacitor to reach 98.2% of its final
voltage.

You might also like