Difference between revisions of "Angular Momentum Conservation"

From wwwelab
Jump to navigation Jump to search
m (Text replacement - "|left" to "|center")
 
(31 intermediate revisions by 3 users not shown)
Line 1: Line 1:
UNDER CONSTRUCTION
+
=Description of the Experiment=
 +
This control room allows the confirmation of angular momentum conservation by making collide a spinning disk with another one. Beyond that, the disk inertia momentum can be extrapolated based on the principles of conservation of energy.
 +
 
 +
<!-- Acho que este texto não está muito correcto. -->
 +
 
 +
 
 +
<div class="toccolours mw-collapsible mw-collapsed" style="width:420px">
 +
'''Links'''
 +
<div class="mw-collapsible-content">
 +
 
 +
*Video: rtsp://elabmc.ist.utl.pt/inertiadisks.sdp
 +
*Laboratory: Intermediate in [http://e-lab.ist.eu e-lab.ist.eu]
 +
*Control room: [unavaliable]
 +
*Level: ***
 +
 
 +
</div>
 +
</div>
 +
 
 +
{{#ev:youtube|X2OK3akY4GE|Slow motion video (12x) of the discs connecting.|center}}
  
=Description of the Experiment=
 
This control room allows the confirmation of angular momentum conservation by colliding a spinning disk with another. Moreover, the disk inertia momentum  can be extrapolated from energy conservation principles.
 
  
 
=Experimental Apparatus=
 
=Experimental Apparatus=
The experimental apparatus is based in a PC hard disk drive motor and its spinning disc with a mass of 115g., 13mm internal radius and 47mm external. A second disc with a mass of 69g and the same dimensions is held on top of it and can be dropped by a servo motor actuator.
+
The experimental apparatus is based on a PC hard disk drive motor and its spinning disc with a mass of 115g, 12.5mm internal radius and 47.5mm external. A second disc with a mass of 69g and the same dimensions is held on top of it and can be dropped by a servo motor actuator.
 +
 
 +
The motor of the apparatus can be used as a generator equipped with a switchable resistor, acting as an electromagnetic brake. The braking current and the voltage characteristic are measured allowing an accurate calculation of energy dissipation.
  
The apparatus' motor can be used as a generator equipped with a switchable resistor acting as an electromagnetic brake. The braking current&voltage characteristic is measured allowing a rigorous energy dissipation calculation.
 
  
 
=Protocol - Angular Momentum Conservation=
 
=Protocol - Angular Momentum Conservation=
  
[[File:Discos_velocidade_protocolo1.png|thumb|alt=|Figura 1: Angular velocity (rpm) as a function of time for a collision at 1000 rpm.]]
+
[[File:Discs_velocity_protocol1.png|thumb|alt=|Figure 1: Angular velocity (rpm) as a function of time for a collision at 1000 rpm.]]
  
The bottom disc is accelerated by the motor until it reaches a selected angular velocity. In this instant the motor is disconnected from supply and the disc allowed to rotate freely. When a certain pre-determined velocity is reached, the servo lets the suspended disc initially at rest fall on top of the rotating disc.
+
The bottom disc is accelerated by the motor until it reaches a selected angular velocity. At this point, the engine is disconnected from power and the disc is allowed to rotate freely. When a certain pre-determined velocity is reached, the servo lets the suspended disc initially at rest fall on top of the rotating disc.
  
 
Data taken from the experiment is given and plotted with the disc velocity in function of time.
 
Data taken from the experiment is given and plotted with the disc velocity in function of time.
  
'''Figure1''' is a plot of the results of an experient in which a servo lets the suspended disc fall when the disc below reaches 1000 rpm.
+
'''Figure1''' is a plot of the results of an experiment in which a servo lets the suspended disc fall when the disc below reaches 1000 rpm.
 +
 
 +
Doing a linear regression between the deceleration and the fall of the disc, it is possible to obtain the predicted velocity at any time. This gives us the rule of thumb of the friction deceleration related to angular velocity.
  
Doing a linear regression between the deceleration and fall of the disc, it is possible to obtain the predicted velocity at any time. This gives us the thumb rule for the friction deceleration related to angular velocity.
 
  
=Physics=
+
=Advanced Protocol - Moment of Inertia Evaluation=  
Using the following quantities:
+
 
 +
[[File:Discos_velocidade_protocolo2.png|thumb|alt=|Figure2: Rotational velocity as function of time after the electromagnetic breaking.|right|border|240px]]
 +
[[File:Discos_tensao_2fases.png|thumb|alt=|Figure3: Circuit schematic for voltage measurement.|right|border|240px]]
 +
[[File:Discos_tensao.png|thumb|alt=|Figure4: Voltage between two phases during breaking.|right|border|240px]]
 +
[[File:W_extrap.gif|thumb|alt=|Figure5: Extrapolation of w based on the initial slope of deceleration.|right|border|240px]]
 +
[[File:Discos_balanco_energetico.png|thumb|alt=|Figure6: Final energy balance showing electrical and mechanical component allowing you to get the total moment of inertia.|right|border|240px]]
 +
 
 +
The bottom disc is accelerated by the motor to a selected angular velocity. At that time the motor supply is disconnected and the disc is allowed to rotate freely. When a certain pre-determined velocity is reached, a relay puts each motor winding in parallel with a resistor  resistance similar to the windings of the motor  ('''Picture 3'''). These resistors will dissipate energy, acting as an electromagnetic brake. Voltage and velocity as functions of time are given in a table of results at the end of the session.
 +
 
 +
'''Figures 2''' and '''4''' are plots obtained through the table of results of an experiment in which the relay turns on when the rotating discs reach 1400 rpm.
 +
 
 +
Using the first data from '''Figure 2''', before the braking, to do a linear regression, you can get the angular deceleration caused by friction, assumed constant all the time, from the slope of the line. From this deceleration it will be possible to calculate  the instantaneous loss of angular momentum deferentially.
 +
 
 +
Between each speed acquisition an energy balance is done . The loss of total mechanical energy must be equal to the sum of losses by friction and electromagnetic breaking.
  
L - angular momentum
+
<math>  \Delta E_{mec} = \Delta E_{friction} + \Delta E_{elec}</math>
  
I - moment of inertia
+
The energy of a rotating body is <math> E_{rot}=\frac{I w^2}{2}</math> I being the moment of inertia. Then, the variation of mechanical energy between each acquisition will be:
  
&omega; - angular velocity
+
<math>\Delta E_{mec}=\frac{I(w_{n+1}^2-w_{n}^2)}{2}</math>
  
m - mass in rotation.
+
<math>w_{n+1}</math> and <math>w_{n}</math> being the angular velocity in two consecutive acquisitions.
  
For the angular momentum conservation:
+
Using the initial slope of deceleration <math>a</math> due to mechanical friction, one can extrapolate <math>w_{n+1}</math> in absence of electromagnetic breaking.
  
$L_i=L_f$
+
<math>w_{n+1}= w_{n} + a \Delta t</math>
  
$I_i \omega_i=I_f \omega_f$
+
Substituting this extrapolated <math>w_{n+1}</math> in the equation of total energy variation it is possible to calculate the energy dissipation due to mechanical friction:
  
$\frac{I_i}{I_f}=\frac{\omega_f}{\omega_i}$
+
<math>\Delta E_{friction}=\frac{I(w_{n}^2+2w_{n}a\Delta t + a^2\Delta t^2-w_{n}^2)}{2}</math>
  
$\frac{\frac{m_i\left (r_1^2+r_2^2 \right )}{2}}{\frac{m_f\left (r_1^2+r_2^2 \right )}{2}}=\frac{\omega_f}{\omega_i}$
+
<math>\Delta E_{friction}=\frac{I(2w_{n}a\Delta t + a^2\Delta t^2)}{2}</math>
  
$\frac{m_i}{m_f}=\frac{\omega_f}{\omega_i}$
+
A set of extrapolations of <math>w_{n+1}</math> can be seen in '''Figure 5'''.
  
The experimental results give:
+
The dissipated power is:
  
$\frac{\omega_f}{\omega_i}=\frac{623}{950}=0,656$
+
<math>P=VI=\frac{V^2}{R}</math>
  
while the predicted mass ratio is
+
The rms voltage across one winding is:
  
$\frac{m_i}{m_f}=\frac{115}{115+69}=0,625$
+
<math>V_{rms}=\frac{V_{measured}}{\sqrt{3}\sqrt{2}}</math>
  
Evaluating the accuracy:
+
In the setup used, the energy dissipates in 3 branches, so the power comes multiplied by 3.  Besides, each winding is in parallel with a resistor with the same resistance value <math>4,7\Omega</math>, which means the power equation will come multiplied by 2 and <math>R=4,7\Omega</math> will be used.
  
$\frac{\left|0,656-0,625\right|}{\left|0,625\right|}\times 100=4,9\%$
+
<math>P=3\times2\times\frac{V_{rms}^2}{R}=3\times2\times\left(\frac{V_{measured}}{\sqrt{3}\sqrt{2}}\right)^2\frac{1}{R}</math>
  
The speed ratio is different from the mass ratio by 4,9% which gives a good approximation for the angular momentum conservation.
+
<math>P=\frac{V^2}{R}</math>
  
Knowing the exact dimensions of the disks ($r_1=13mm, r_2=47mm$) and adding an error momentum on the equations one can infer an approximated value for the motor rotor momentum of inertia (or its mass knowing its average radius).
+
The energy dissipated will be:
  
$I_i \omega_i=I_f \omega_f$
+
<math>\Delta E_{ele}=P*\Delta t</math>
  
$\left (I_m + I_{Di}\right ) \omega_i=\left (I_m + I_{Df}\right ) \omega_f$
+
Where <math>\Delta t</math> is the time between acquisitions.
  
Solving in order to $I_m$
+
The ballance between each consecutive acquitition is summed in the end.
  
$I_m = \frac{I_{Df} \omega_f - I_{Di} \omega_i}{\omega_i-\omega_f}$
+
<math>Balance = \Delta E_{mec} - \Delta E_{friction} - \Delta E_{elec}</math>
  
=Advanced Protocol - Moment of Inertia Evaluation=
+
Finally the goal-seek function of Microsoft Excel is used to force the sum of balances to be 0 (zero) changing the value of I.
  
[[File:Discos_velocidade_protocolo2.png|thumb|alt=|Figura2: rotational speed]]
+
Using this method, we reached an experimental value of <math>1,274\times10^{-4} kg m^2</math>  for the moment of inertia.
[[File:Discos_tensao_2fases.png|thumb|alt=|Figura3: voltage measurement]]
 
[[File:Discos_tensao.png|thumb|alt=|Figura4: voltage between two phases]]
 
[[File:Discos_balanco_energetico.png|thumb|alt=|Figura5: energy balance]]
 
  
A disc with a total mass of 115g is accelerated by the hard-drive motor until it reaches a pre-defined angular velocity. At that time the motor is disconnected and the disc allowed to rotate freely. When a certain speed that the user defines previously is reached, a relay puts each motor winding in parallel with a resistor which resistance is the same as the motor's windings ('''Figure 3'''). These resistors will dissipate energy acting as an electromagnetic brake. Voltage and speed in funtion of time are given in a table of results in the end of the session.
+
'''Figure 6''' shows the energy of the disc in function of time, the energy lost by friction and electromagnetic breaking and the sum of all the energies that allows the verification of the conservation of energy throughout the experiment.
  
'''Figures 2''' and '''4''' are plots created in Microsoft Excel using the table of results of an experiment in which the relay turns on when the rotating discs reach 1400 rpm.
+
The disc is in fact a ring having an interior radius of 13mm and an exterior radius of 47mm, so its theoretical moment of inertia is:
  
Using the first data to do a linear regression, it is possible to get the speed of the motor over time if the relay didn't turn on.
+
<math>I=\frac{m\left(r_1^2+r_2^2\right)}{2}=\frac{0,115\left(0,013^2+0,047^2\right)}{2}=1,367\times 10^{-4}kg \; m^2</math>
  
Between each speed acquisition it is done an energy balance. The loss of total mechanical energy must be equal to the sum of losses by friction and electromagnetic breaking.
+
Evaluating the accuracy:
  
$\Delta E_{mec} = \Delta E_{fr} + \Delta E_{ele}$
+
<math>\frac{\left|1,274\times 10^{-4}-1,367\times 10^{-4}\right|}{\left|1,367\times 10^{-4}\right|}\times 100=6,8\%</math>
  
The energy of a rotating body is $E_{rot}=\frac{I w^2}{2}$ I being the moment of inertia. Then, the variation of mechanical energy between each acquisition will be:
+
This experiment gave results that differ by '''6,8%''' from the ones calculated theoretically.
  
$\Delta E_{mec}=\frac{I(w_{2exp}^2-w_{1exp}^2)}{2}$
+
=Physics=
 +
Using the following quantities:
  
$w_{2exp}$ and $w_{1exp}$ being the angular velocity in two consecutive acquisitions.
+
L - angular momentum
  
The loss of energy by friction will be:
+
I - moment of inertia
  
$\Delta E_{fr}=Iw_{exp}\left(w_{2wo/fr}-w_{1wo/fr}\right)$
+
&omega; - angular velocity
  
$w_{exp}$ is the angular velocity of the disc in that acquisition and $w_{2wo/fr}$ $w_{1wo/fr}$ are the extrapolated velocity of the disc in a no friction situation at the time of the acquition and the previous one respectively.
+
m - mass in rotation.
  
The dissipated power is:
+
For the angular momentum conservation:
  
$P=VI=\frac{V^2}{R}$
+
<math>L_i=L_f</math>
  
The rms voltage across one winding is:
+
<math>I_i \omega_i=I_f \omega_f</math>
  
$V_{rms}=\frac{V_{measured}}{\sqrt{3}\sqrt{2}}$
+
<math>\frac{I_i}{I_f}=\frac{\omega_f}{\omega_i}</math>
  
In the setup used the energy dissipates in 3 branches so the power comes multiplied by 3. Also, each winding is in parallel with a resistor with the same resistence value $4,7\Omega$, which means the power equation will come multiplied by 2 and $R=4,7\Omega$ will be used.
+
<math>\frac{\frac{m_i\left (r_1^2+r_2^2 \right )}{2}}{\frac{m_f\left (r_1^2+r_2^2 \right )}{2}}=\frac{\omega_f}{\omega_i}</math>
  
$P=3\times2\times\frac{V_{rms}^2}{R}=3\times2\times\left(\frac{V_{measured}}{\sqrt{3}\sqrt{2}}\right)^2\frac{1}{R}$
+
<math>\frac{m_i}{m_f}=\frac{\omega_f}{\omega_i}</math>
  
$P=\frac{V^2}{R}$
+
The experimental results give:
  
The energy dissipated will be:
+
<math>\frac{\omega_f}{\omega_i}=\frac{623}{950}=0,656</math>
  
$\Delta E_{ele}=P*\Delta t$
+
while the predicted mass ratio is
  
Where $\Delta t$ is the time between acquisitions.
+
<math>\frac{m_i}{m_f}=\frac{115}{115+69}=0,625</math>
  
In the end the ballance between each consecutive acquitition is summed.
+
Evaluating the accuracy:
  
$Balance = \Delta E_{mec} - \Delta E_{friction} - \Delta E_{ele}$
+
<math>\frac{\left|0,656-0,625\right|}{\left|0,625\right|}\times 100=4,9\%</math>
  
Finaly the goal-seek function of Microsoft Excel is used to force the sum of balances to be 0 (zero) changing the value of I.
+
The speed ratio is different from the mass ratio by 4,9%, which gives a good approximation of the angular momentum conservation.
  
Using this method, it was reached an experimental value of $1,274\times10^{-4}$ $kg$ $m^2$ for the moment of inertia.
+
Knowing the exact dimensions of the disks (<math>r_1=12,5mm, r_2=47,5mm</math>) and adding an error momentum to the equations, you can infer an approximated value of the motor rotor momentum of inertia (or its mass knowing its average radius).
  
'''Figure 5''' shows the energy of the disc in function of time, the energy lost by friction and electromagnetic breaking and the sum of all energies that allows the verification of the conservation of energy through all the experiment.
+
<math>I_i \omega_i=I_f \omega_f</math>
  
The disc is in fact a ring with interior radius 13mm and exterior 47mm so it's theoretical moment of inertia is:
+
<math>\left (I_m + I_{Di}\right ) \omega_i=\left (I_m + I_{Df}\right ) \omega_f</math>
  
$I=\frac{m\left(r_1^2+r_2^2\right)}{2}=\frac{0,115\left(0,013^2+0,047^2\right)}{2}=1,367\times 10^{-4}kg \; m^2$
+
Solving in for <math>I_m</math>
  
Evaluating the accuracy:
+
<math>I_m = \frac{I_{Df} \omega_f - I_{Di} \omega_i}{\omega_i-\omega_f}</math>
  
$\frac{\left|1,274\times 10^{-4}-1,367\times 10^{-4}\right|}{\left|1,367\times 10^{-4}\right|}\times 100=6,8\%$
 
  
This experiment gave results that differ by '''6,8%''' from the ones calculated theoreticaly.
+
=Links=
 +
*[[Conservação do Momento Angular | Portuguese version (Versão em Português)]]

Latest revision as of 14:12, 28 May 2015

Description of the Experiment

This control room allows the confirmation of angular momentum conservation by making collide a spinning disk with another one. Beyond that, the disk inertia momentum can be extrapolated based on the principles of conservation of energy.


Links

  • Video: rtsp://elabmc.ist.utl.pt/inertiadisks.sdp
  • Laboratory: Intermediate in e-lab.ist.eu
  • Control room: [unavaliable]
  • Level: ***


Experimental Apparatus

The experimental apparatus is based on a PC hard disk drive motor and its spinning disc with a mass of 115g, 12.5mm internal radius and 47.5mm external. A second disc with a mass of 69g and the same dimensions is held on top of it and can be dropped by a servo motor actuator.

The motor of the apparatus can be used as a generator equipped with a switchable resistor, acting as an electromagnetic brake. The braking current and the voltage characteristic are measured allowing an accurate calculation of energy dissipation.


Protocol - Angular Momentum Conservation

Figure 1: Angular velocity (rpm) as a function of time for a collision at 1000 rpm.

The bottom disc is accelerated by the motor until it reaches a selected angular velocity. At this point, the engine is disconnected from power and the disc is allowed to rotate freely. When a certain pre-determined velocity is reached, the servo lets the suspended disc initially at rest fall on top of the rotating disc.

Data taken from the experiment is given and plotted with the disc velocity in function of time.

Figure1 is a plot of the results of an experiment in which a servo lets the suspended disc fall when the disc below reaches 1000 rpm.

Doing a linear regression between the deceleration and the fall of the disc, it is possible to obtain the predicted velocity at any time. This gives us the rule of thumb of the friction deceleration related to angular velocity.


Advanced Protocol - Moment of Inertia Evaluation

Figure2: Rotational velocity as function of time after the electromagnetic breaking.
Figure3: Circuit schematic for voltage measurement.
Figure4: Voltage between two phases during breaking.
Figure5: Extrapolation of w based on the initial slope of deceleration.
Figure6: Final energy balance showing electrical and mechanical component allowing you to get the total moment of inertia.

The bottom disc is accelerated by the motor to a selected angular velocity. At that time the motor supply is disconnected and the disc is allowed to rotate freely. When a certain pre-determined velocity is reached, a relay puts each motor winding in parallel with a resistor resistance similar to the windings of the motor (Picture 3). These resistors will dissipate energy, acting as an electromagnetic brake. Voltage and velocity as functions of time are given in a table of results at the end of the session.

Figures 2 and 4 are plots obtained through the table of results of an experiment in which the relay turns on when the rotating discs reach 1400 rpm.

Using the first data from Figure 2, before the braking, to do a linear regression, you can get the angular deceleration caused by friction, assumed constant all the time, from the slope of the line. From this deceleration it will be possible to calculate the instantaneous loss of angular momentum deferentially.

Between each speed acquisition an energy balance is done . The loss of total mechanical energy must be equal to the sum of losses by friction and electromagnetic breaking.

[math] \Delta E_{mec} = \Delta E_{friction} + \Delta E_{elec}[/math]

The energy of a rotating body is [math] E_{rot}=\frac{I w^2}{2}[/math] I being the moment of inertia. Then, the variation of mechanical energy between each acquisition will be:

[math]\Delta E_{mec}=\frac{I(w_{n+1}^2-w_{n}^2)}{2}[/math]

[math]w_{n+1}[/math] and [math]w_{n}[/math] being the angular velocity in two consecutive acquisitions.

Using the initial slope of deceleration [math]a[/math] due to mechanical friction, one can extrapolate [math]w_{n+1}[/math] in absence of electromagnetic breaking.

[math]w_{n+1}= w_{n} + a \Delta t[/math]

Substituting this extrapolated [math]w_{n+1}[/math] in the equation of total energy variation it is possible to calculate the energy dissipation due to mechanical friction:

[math]\Delta E_{friction}=\frac{I(w_{n}^2+2w_{n}a\Delta t + a^2\Delta t^2-w_{n}^2)}{2}[/math]

[math]\Delta E_{friction}=\frac{I(2w_{n}a\Delta t + a^2\Delta t^2)}{2}[/math]

A set of extrapolations of [math]w_{n+1}[/math] can be seen in Figure 5.

The dissipated power is:

[math]P=VI=\frac{V^2}{R}[/math]

The rms voltage across one winding is:

[math]V_{rms}=\frac{V_{measured}}{\sqrt{3}\sqrt{2}}[/math]

In the setup used, the energy dissipates in 3 branches, so the power comes multiplied by 3. Besides, each winding is in parallel with a resistor with the same resistance value [math]4,7\Omega[/math], which means the power equation will come multiplied by 2 and [math]R=4,7\Omega[/math] will be used.

[math]P=3\times2\times\frac{V_{rms}^2}{R}=3\times2\times\left(\frac{V_{measured}}{\sqrt{3}\sqrt{2}}\right)^2\frac{1}{R}[/math]

[math]P=\frac{V^2}{R}[/math]

The energy dissipated will be:

[math]\Delta E_{ele}=P*\Delta t[/math]

Where [math]\Delta t[/math] is the time between acquisitions.

The ballance between each consecutive acquitition is summed in the end.

[math]Balance = \Delta E_{mec} - \Delta E_{friction} - \Delta E_{elec}[/math]

Finally the goal-seek function of Microsoft Excel is used to force the sum of balances to be 0 (zero) changing the value of I.

Using this method, we reached an experimental value of [math]1,274\times10^{-4} kg m^2[/math] for the moment of inertia.

Figure 6 shows the energy of the disc in function of time, the energy lost by friction and electromagnetic breaking and the sum of all the energies that allows the verification of the conservation of energy throughout the experiment.

The disc is in fact a ring having an interior radius of 13mm and an exterior radius of 47mm, so its theoretical moment of inertia is:

[math]I=\frac{m\left(r_1^2+r_2^2\right)}{2}=\frac{0,115\left(0,013^2+0,047^2\right)}{2}=1,367\times 10^{-4}kg \; m^2[/math]

Evaluating the accuracy:

[math]\frac{\left|1,274\times 10^{-4}-1,367\times 10^{-4}\right|}{\left|1,367\times 10^{-4}\right|}\times 100=6,8\%[/math]

This experiment gave results that differ by 6,8% from the ones calculated theoretically.

Physics

Using the following quantities:

L - angular momentum

I - moment of inertia

ω - angular velocity

m - mass in rotation.

For the angular momentum conservation:

[math]L_i=L_f[/math]

[math]I_i \omega_i=I_f \omega_f[/math]

[math]\frac{I_i}{I_f}=\frac{\omega_f}{\omega_i}[/math]

[math]\frac{\frac{m_i\left (r_1^2+r_2^2 \right )}{2}}{\frac{m_f\left (r_1^2+r_2^2 \right )}{2}}=\frac{\omega_f}{\omega_i}[/math]

[math]\frac{m_i}{m_f}=\frac{\omega_f}{\omega_i}[/math]

The experimental results give:

[math]\frac{\omega_f}{\omega_i}=\frac{623}{950}=0,656[/math]

while the predicted mass ratio is

[math]\frac{m_i}{m_f}=\frac{115}{115+69}=0,625[/math]

Evaluating the accuracy:

[math]\frac{\left|0,656-0,625\right|}{\left|0,625\right|}\times 100=4,9\%[/math]

The speed ratio is different from the mass ratio by 4,9%, which gives a good approximation of the angular momentum conservation.

Knowing the exact dimensions of the disks ([math]r_1=12,5mm, r_2=47,5mm[/math]) and adding an error momentum to the equations, you can infer an approximated value of the motor rotor momentum of inertia (or its mass knowing its average radius).

[math]I_i \omega_i=I_f \omega_f[/math]

[math]\left (I_m + I_{Di}\right ) \omega_i=\left (I_m + I_{Df}\right ) \omega_f[/math]

Solving in for [math]I_m[/math]

[math]I_m = \frac{I_{Df} \omega_f - I_{Di} \omega_i}{\omega_i-\omega_f}[/math]


Links