Jump to content

Time constant: Difference between revisions

From Wikipedia, the free encyclopedia
Content deleted Content added
No edit summary
Tag: Reverted
No edit summary
 
(25 intermediate revisions by 15 users not shown)
Line 1: Line 1:
{{Short description|Constant in physics and engineering}}
In [[physics]] and [[engineering]], the '''better luck next time''', usually denoted by the [[Greek language|Greek]] letter {{mvar|τ}} (tau), is the [[parameter]] characterizing the response to a step input of a first-order, [[LTI system theory|linear time-invariant]] (LTI) system.<ref name= Lipták>{{cite book |title=Instrument Engineers' Handbook: Process control and optimization |author=Béla G. Lipták |url=https://books.google.com/books?id=pPMursVsxlMC&pg=PA100 |page=100 |isbn=978-0-8493-1081-2 |year=2003 |publisher=CRC Press |edition=4}}</ref><ref group="note">Concretely, a first-order LTI system is a system that can be modeled by a single [[Ordinary differential equation|first order differential equation]] in time. Examples include the simplest single-stage electrical [[RC circuit#Series circuit|RC circuit]]s and [[RL circuit#Series circuit|RL circuit]]s.</ref> The better luck next time is the main [[characteristic unit]] of a first-order LTI system.
{{Technical|date=February 2023}}

In [[physics]] and [[engineering]], the '''time constant''', usually denoted by the [[Greek language|Greek]] letter {{mvar|τ}} (tau), is the [[parameter]] characterizing the response to a step input of a first-order, [[LTI system theory|linear time-invariant]] (LTI) system.<ref name= Lipták>{{cite book |title=Instrument Engineers' Handbook: Process control and optimization |author=Béla G. Lipták |url=https://books.google.com/books?id=pPMursVsxlMC&pg=PA100 |page=100 |isbn=978-0-8493-1081-2 |year=2003 |publisher=CRC Press |edition=4}}</ref><ref group="note">Concretely, a first-order LTI system is a system that can be modeled by a single [[Ordinary differential equation|first order differential equation]] in time. Examples include the simplest single-stage electrical [[RC circuit#Series circuit|RC circuit]]s and [[RL circuit#Series circuit|RL circuit]]s.</ref> The time constant is the main [[characteristic unit]] of a first-order LTI system. It gives speed of the response.


In the time domain, the usual choice to explore the time response is through the [[step response]] to a [[Heaviside step function|step input]], or the [[impulse response]] to a [[Dirac delta function]] input.<ref name=Wie>{{cite book |title=Space vehicle dynamics and control |url=https://archive.org/details/spacevehicledyna00wieb_0 |url-access=registration |page=[https://archive.org/details/spacevehicledyna00wieb_0/page/100 100] |author=Bong Wie |isbn=978-1-56347-261-9 |year=1998 |publisher=American Institute of Aeronautics and Astronautics}}</ref> In the frequency domain (for example, looking at the [[Fourier transform]] of the step response, or using an input that is a simple sinusoidal function of time) the time constant also determines the [[Bandwidth (signal processing)|bandwidth]] of a first-order time-invariant system, that is, the frequency at which the output signal power drops to half the value it has at low frequencies.
In the time domain, the usual choice to explore the time response is through the [[step response]] to a [[Heaviside step function|step input]], or the [[impulse response]] to a [[Dirac delta function]] input.<ref name=Wie>{{cite book |title=Space vehicle dynamics and control |url=https://archive.org/details/spacevehicledyna00wieb_0 |url-access=registration |page=[https://archive.org/details/spacevehicledyna00wieb_0/page/100 100] |author=Bong Wie |isbn=978-1-56347-261-9 |year=1998 |publisher=American Institute of Aeronautics and Astronautics}}</ref> In the frequency domain (for example, looking at the [[Fourier transform]] of the step response, or using an input that is a simple sinusoidal function of time) the time constant also determines the [[Bandwidth (signal processing)|bandwidth]] of a first-order time-invariant system, that is, the frequency at which the output signal power drops to half the value it has at low frequencies.
Line 7: Line 10:
Time constants are a feature of the [[lumped system analysis]] (lumped capacity analysis method) for thermal systems, used when objects cool or warm uniformly under the influence of convective cooling or warming.<ref name=NATO>{{cite book |title=Physically-based Modelling and Simulation of Climate and Climatic Change |editor=Michael E. Schlesinger |agency=NATO |chapter-url=https://books.google.com/books?id=oDsv-wPWdkUC&pg=PA63 |page=627 |author=GR North |chapter=Lessons from energy balance models |edition=NATO Advanced Study Institute on Physical-Based Modelling |publisher=Springer |year=1988 |isbn=978-90-277-2789-3}}</ref>
Time constants are a feature of the [[lumped system analysis]] (lumped capacity analysis method) for thermal systems, used when objects cool or warm uniformly under the influence of convective cooling or warming.<ref name=NATO>{{cite book |title=Physically-based Modelling and Simulation of Climate and Climatic Change |editor=Michael E. Schlesinger |agency=NATO |chapter-url=https://books.google.com/books?id=oDsv-wPWdkUC&pg=PA63 |page=627 |author=GR North |chapter=Lessons from energy balance models |edition=NATO Advanced Study Institute on Physical-Based Modelling |publisher=Springer |year=1988 |isbn=978-90-277-2789-3}}</ref>


Physically, the time constant represents the elapsed time required for the system response to decay to zero if the system had continued to decay at the initial rate, because of the progressive change in the rate of decay the response will have actually decreased in value to {{math|1&thinsp;/&thinsp;''e'' ≈ 36.8%}} in this time (say from a step decrease). In an increasing system, the time constant is the time for the system's [[step response]] to reach {{math|1 − 1&thinsp;/&thinsp;''e'' ≈ 63.2%}} of its final (asymptotic) value (say from a step increase). In radioactive decay the time constant is related to the [[decay constant]] ('''λ'''), and it represents both the mean lifetime of a decaying system (such as an atom) before it decays, or the time it takes for all but 36.8% of the atoms to decay. For this reason, the time constant is longer than the [[half-life]], which is the time for only 50% of the atoms to decay.
Physically, the time constant represents the elapsed time required for the system response to decay to zero if the system had continued to decay at the initial rate, because of the progressive change in the rate of decay the response will have actually decreased in value to {{math|1&thinsp;/&thinsp;''e'' ≈ 36.8%}} in this time (say from a step decrease). In an increasing system, the time constant is the time for the system's [[step response]] to reach {{math|1 − 1&thinsp;/&thinsp;''e'' ≈ 63.2%}} of its final (asymptotic) value (say from a step increase). In radioactive decay the time constant is related to the [[decay constant]] ('''λ'''), and it represents both the mean lifetime of a decaying system (such as an atom) before it decays, or the time it takes for all but 36.8% of the atoms to decay. For this reason, the time constant is longer than the [[half-life]], which is the time for only 50% of the atoms to decay.


==Differential equation==
==Differential equation==
{{Main|LTI system theory}}
{{Main|Linear time-invariant system}}
First order LTI systems are characterized by the differential equation
First order LTI systems are characterized by the differential equation
:<math> \tau \frac{dV}{dt} + V = f(t) </math>
<math display="block"> \tau \frac{dV}{dt} + V = f(t) </math>


where {{mvar|τ}} represents the [[exponential decay]] constant and {{mvar|V}} is a function of time {{mvar|t}}
where {{mvar|τ}} represents the [[exponential decay]] constant and {{mvar|V}} is a function of time {{mvar|t}}
:<math> V = V(t). </math>
<math display="block"> V = V(t). </math>
The right-hand side is the ''forcing function'' {{math|''f''(''t'')}} describing an external driving function of time, which can be regarded as the system ''input'', to which {{math|''V''(''t'')}} is the ''response'', or system output. Classical examples for {{math|''f''(''t'')}} are:
The right-hand side is the ''forcing function'' {{math|''f''(''t'')}} describing an external driving function of time, which can be regarded as the system ''input'', to which {{math|''V''(''t'')}} is the ''response'', or system output. Classical examples for {{math|''f''(''t'')}} are:


The [[Heaviside step function]], often denoted by {{math|''u''(''t'')}}:
The [[Heaviside step function]], often denoted by {{math|''u''(''t'')}}:
:<math>u(t)=\begin{cases} 0, & t < 0 \\ 1, & t \ge 0 \end{cases} </math>
<math display="block">u(t) = \begin{cases} 0, & t < 0 \\ 1, & t \ge 0 \end{cases} </math>
the [[Dirac delta function|impulse function]], often denoted by {{math|''δ''(''t'')}}, and also the sinusoidal input function:
the [[Dirac delta function|impulse function]], often denoted by {{math|''δ''(''t'')}}, and also the sinusoidal input function:
:<math> f(t) = A \sin(2 \pi f t) </math>
<math display="block"> f(t) = A \sin(2 \pi f t) </math>
or
or
:<math> f(t) = A e^{j \omega t }, </math>
<math display="block"> f(t) = A e^{j \omega t }, </math>
where {{mvar|A}} is the amplitude of the forcing function, {{mvar|f}} is the frequency in Hertz, and {{math|1=''ω'' = 2''π f''}} is the frequency in radians per second.
where {{mvar|A}} is the amplitude of the forcing function, {{mvar|f}} is the frequency in Hertz, and {{math|1=''ω'' = 2''π f''}} is the frequency in radians per second.


===Example solution===
===Example solution===
An example solution to the differential equation with initial value {{math|''V''<sub>0</sub>}} and no forcing function is
An example solution to the differential equation with initial value {{math|''V''<sub>0</sub>}} and no forcing function is
:<math> V(t) = V_0 e^{-t / \tau} </math>
<math display="block"> V(t) = V_0 e^{-t / \tau} </math>


where
where
:<math> V_0 = V(t=0) </math>
<math display="block"> V_0 = V(t=0) </math>


is the initial value of {{mvar|V}}. Thus, the response is an exponential decay with time constant {{mvar|τ}}.
is the initial value of {{mvar|V}}. Thus, the response is an exponential decay with time constant {{mvar|τ}}.
Line 38: Line 41:


Suppose
Suppose
<math display="block"> V(t) = V_0 e^{-t / \tau}. </math>
<math display="block"> V(t) = V_0 e^{-t / \tau}. </math>


This behavior is referred to as a "decaying" exponential function. The time {{mvar|τ}} (tau) is referred to as the "time constant" and can be used (as in this case) to indicate how rapidly an exponential function decays.
This behavior is referred to as a "decaying" exponential function. The time {{mvar|τ}} (tau) is referred to as the "time constant" and can be used (as in this case) to indicate how rapidly an exponential function decays.
Line 49: Line 52:


# Let <math>t=0</math>; then <math>V = V_0 e^0</math>, and so <math>V=V_0</math>
# Let <math>t=0</math>; then <math>V = V_0 e^0</math>, and so <math>V=V_0</math>
# Let <math>t= \tau</math>; then <math>V=V_0 e^{-1} \approx 0.37 V_0 </math>
# Let <math>t= \tau</math>; then <math>V=V_0 e^{-1} \approx 0.37 V_0 </math>
# Let <math>V = f(t) = V_0 e^{-t / \tau}</math>, and so <math display="inline"> \lim_{t \to \infty}f(t) = 0 </math>
# Let <math>V = f(t) = V_0 e^{-t / \tau}</math>, and so <math display="inline"> \lim_{t \to \infty}f(t) = 0 </math>
# Let <math>t=5 \tau</math>; then <math>V = V_0 e^{-5} \approx 0.0067V_0 </math>
# Let <math>t=5 \tau</math>; then <math>V = V_0 e^{-5} \approx 0.0067V_0 </math>
Line 59: Line 62:
Suppose the forcing function is chosen as sinusoidal so:
Suppose the forcing function is chosen as sinusoidal so:


:<math> \tau \frac{dV}{dt} + V = f(t) = Ae^{j \omega t }. </math>
<math display="block"> \tau \frac{dV}{dt} + V = f(t) = Ae^{j \omega t }. </math>


(Response to a real cosine or sine wave input can be obtained by taking the real or imaginary part of the final result by virtue of [[Euler's formula]].) The general solution to this equation for times {{math|''t'' ≥ 0 s}}, assuming {{math|1=''V''(''t'' = 0) = ''V''<sub>0</sub>}} is:
(Response to a real cosine or sine wave input can be obtained by taking the real or imaginary part of the final result by virtue of [[Euler's formula]].) The general solution to this equation for times {{math|''t'' ≥ 0 s}}, assuming {{math|1=''V''(''t'' = 0) = ''V''<sub>0</sub>}} is:


:<math>\begin{align}
<math display="block">\begin{align}
V(t) &= V_0e^{-t/\tau} + \frac{Ae^{-t/\tau}}{\tau}\int_0^t \, dt'\ e^{t'/\tau} e^{j\omega t'} \\
V(t) &= V_0e^{-t/\tau} + \frac{Ae^{-t/\tau}}{\tau}\int_0^t \, dt'\ e^{t'/\tau} e^{j\omega t'} \\[1ex]
&= V_0 e^{-t/\tau} + \frac{1 / \tau}{j\omega +1/\tau} A\left( e^{j \omega t} - e^{-t/\tau}\right).
&= V_0 e^{-t/\tau} + \frac{\frac 1 \tau}{j\omega + \frac 1 \tau} A\left( e^{j \omega t} - e^{-t/\tau}\right).
\end{align}</math>
\end{align}</math>


For long times the decaying exponentials become negligible and the [[steady-state]] solution or long-time solution is:
For long times the decaying exponentials become negligible and the [[steady-state]] solution or long-time solution is:


:<math> V_{\infty}(t) = \frac{1/\tau}{j\omega +1/\tau}Ae^{j \omega t}.</math>
<math display="block"> V_{\infty}(t) = \frac{1/\tau}{j\omega +1/\tau}Ae^{j \omega t}.</math>


The magnitude of this response is:
The magnitude of this response is:
:<math> |V_{\infty}(t)| = A\frac{1}{\tau\left(\omega^2 +(1/\tau)^2\right)^{1/2}} = A \frac{1}{\sqrt{1+(\omega \tau)^2 }}.</math>
<math display="block"> |V_{\infty}(t)| = A\frac{1}{\tau\left(\omega^2 +(1/\tau)^2\right)^{1/2}} = A \frac{1}{\sqrt{1+(\omega \tau)^2 }}.</math>
By convention, the bandwidth of this system is the frequency where {{math|{{!}}''V''<sub>∞</sub>{{!}}<sup>2</sup>}} drops to half-value, or where {{math|1=''ωτ'' = 1}}. This is the usual [[Bandwidth (signal processing)|bandwidth]] convention, defined as the frequency range where power drops by less than half (at most &minus;3&nbsp;dB). Using the frequency in hertz, rather than radians/s ({{math|1=''ω'' = 2''πf''}}):
By convention, the bandwidth of this system is the frequency where {{math|{{!}}''V''<sub>∞</sub>{{!}}<sup>2</sup>}} drops to half-value, or where {{math|1=''ωτ'' = 1}}. This is the usual [[Bandwidth (signal processing)|bandwidth]] convention, defined as the frequency range where power drops by less than half (at most &minus;3&nbsp;dB). Using the frequency in hertz, rather than radians/s ({{math|1=''ω'' = 2''πf''}}):


:<math> f_\mathrm{3dB} = \frac {1}{2 \pi \tau}. </math>
<math display="block"> f_\mathrm{3dB} = \frac {1}{2 \pi \tau}. </math>


The notation {{math|''f''<sub>3dB</sub>}} stems from the expression of power in [[decibels]] and the observation that half-power corresponds to a drop in the value of {{math|{{!}}''V''<sub>∞</sub>{{!}}}} by a factor of 1/√2 or by 3 decibels.
The notation {{math|''f''<sub>3dB</sub>}} stems from the expression of power in [[decibels]] and the observation that half-power corresponds to a drop in the value of {{math|{{!}}''V''<sub>∞</sub>{{!}}}} by a factor of 1/2 or by 3 decibels.


Thus, the time constant determines the bandwidth of this system.
Thus, the time constant determines the bandwidth of this system.
Line 87: Line 90:
Suppose the forcing function is chosen as a step input so:
Suppose the forcing function is chosen as a step input so:


:<math> \frac{dV}{dt} + \frac{1}{\tau} V = f(t) = A u(t), </math>
<math display="block"> \frac{dV}{dt} + \frac{1}{\tau} V = f(t) = A u(t), </math>


with {{math|''u''(''t'')}} the Heaviside step function. The general solution to this equation for times {{math|''t'' ≥ 0 s}}, assuming {{math|1=''V''(''t'' = 0) = ''V''<sub>0</sub>}} is:
with {{math|''u''(''t'')}} the [[unit step function]]. The general solution to this equation for times {{math|''t'' ≥ 0 s}}, assuming {{math|1=''V''(''t'' = 0) = ''V''<sub>0</sub>}} is:


:<math> V(t) = V_0 e^{-t/\tau} + A \tau \left( 1 - e^{-t/\tau}\right).</math>
<math display="block"> V(t) = V_0 e^{-t/\tau} + A \tau \left( 1 - e^{-t/\tau}\right).</math>


(It may be observed that this response is the {{math|''ω'' → 0}} limit of the above response to a sinusoidal input.)
(It may be observed that this response is the {{math|''ω'' → 0}} limit of the above response to a sinusoidal input.)


The long-time solution is time independent and independent of initial conditions:
The long-time solution is time independent and independent of initial conditions:
:<math> V_{\infty} = A \tau. </math>
<math display="block"> V_{\infty} = A \tau. </math>


The time constant remains the same for the same system regardless of the starting conditions. Simply stated, a system approaches its final, steady-state situation at a constant rate, regardless of how close it is to that value at any arbitrary starting point.
The time constant remains the same for the same system regardless of the starting conditions. Simply stated, a system approaches its final, steady-state situation at a constant rate, regardless of how close it is to that value at any arbitrary starting point.


For example, consider an electric motor whose startup is well modelled by a first-order LTI system. Suppose that when started from rest, the motor takes {{sfrac|1|8}} of a second to reach 63% of its nominal speed of 100 RPM, or 63 RPM&mdash;a shortfall of 37 RPM. Then it will be found that after the next {{sfrac|1|8}} of a second, the motor has sped up an additional 23 RPM, which equals 63% of that 37 RPM difference. This brings it to 86 RPM&mdash;still 14 RPM low. After a third {{sfrac|1|8}} of a second, the motor will have gained an additional 9 RPM (63% of that 14 RPM difference), putting it at 95 RPM.
For example, consider an electric motor whose startup is well modelled by a first-order LTI system. Suppose that when started from rest, the motor takes {{sfrac|1|8}} of a second to reach 63% of its nominal speed of 100 RPM, or 63 RPM&mdash;a shortfall of 37 RPM. Then it will be found that after the next {{sfrac|1|8}} of a second, the motor has sped up an additional 23 RPM, which equals 63% of that 37 RPM difference. This brings it to 86 RPM&mdash;still 14 RPM low. After a third {{sfrac|1|8}} of a second, the motor will have gained an additional 9 RPM (63% of that 14 RPM difference), putting it at 95 RPM.


In fact, given ''any'' initial speed {{nowrap|''s'' ≤ 100 RPM,}} {{sfrac|1|8}} of a second later this particular motor will have gained an additional {{nowrap|0.63 × (100 − ''s'') RPM.}}
In fact, given ''any'' initial speed {{nowrap|''s'' ≤ 100 RPM,}} {{sfrac|1|8}} of a second later this particular motor will have gained an additional {{nowrap|0.63 × (100 − ''s'') RPM.}}
Line 110: Line 113:
[[File:Series RC resistor voltage.svg|thumb|right|230px|Inductor voltage step-response.]]
[[File:Series RC resistor voltage.svg|thumb|right|230px|Inductor voltage step-response.]]


In an [[RL circuit]] composed of a single resistor and inductor, the time constant ''<math>\tau</math>'' (in [[second]]s) is
In an [[RL circuit]] composed of a single resistor and inductor, the time constant <math>\tau</math> (in [[second]]s) is
:<math> \tau = \frac{ L }{ R } </math>
<math display="block"> \tau = \frac{ L }{ R } </math>


where ''R'' is the [[Electrical resistance|resistance]] (in [[ohm]]s) and ''L'' is the [[inductance]] (in [[Henry (unit)|Henrys]]).
where ''R'' is the [[Electrical resistance|resistance]] (in [[ohm]]s) and ''L'' is the [[inductance]] (in [[Henry (unit)|henrys]]).


Similarly, in an [[RC circuit]] composed of a single resistor and capacitor, the time constant <math>\tau</math> (in seconds) is:
Similarly, in an [[RC circuit]] composed of a single resistor and capacitor, the time constant <math>\tau</math> (in seconds) is:
:<math> \tau = R C </math>
<math display="block"> \tau = R C </math>


where ''R'' is the resistance (in [[ohm]]s) and ''C'' is the [[capacitance]] (in [[farad]]s).
where ''R'' is the resistance (in [[ohm]]s) and ''C'' is the [[capacitance]] (in [[farad]]s).
Line 122: Line 125:
Electrical circuits are often more complex than these examples, and may exhibit multiple time constants (See [[Step response]] and [[Pole splitting]] for some examples.) In the case where [[Negative feedback amplifier|feedback]] is present, a system may exhibit unstable, increasing oscillations. In addition, physical electrical circuits are seldom truly linear systems except for very low amplitude excitations; however, the approximation of linearity is widely used.
Electrical circuits are often more complex than these examples, and may exhibit multiple time constants (See [[Step response]] and [[Pole splitting]] for some examples.) In the case where [[Negative feedback amplifier|feedback]] is present, a system may exhibit unstable, increasing oscillations. In addition, physical electrical circuits are seldom truly linear systems except for very low amplitude excitations; however, the approximation of linearity is widely used.


In digital electronic circuits another measure, the [[fanout of 4|FO4]] is often used. This can be converted to time constant units via the equation <math>5\tau = \text{FO4}</math>.<ref>{{cite book|doi=10.1109/ACSSC.2003.1292037|chapter=Logical effort of carry propagate adders|title=The Thirty-Seventh Asilomar Conference on Signals, Systems & Computers, 2003|pages=873–878|year=2003|last1=Harris|first1=D.|last2=Sutherland|first2=I.|isbn=0-7803-8104-1|s2cid=7880203}}</ref>
In digital electronic circuits another measure, the [[fanout of 4|FO4]] is often used. This can be converted to time constant units via the equation <math>5\tau = \text{FO4}</math>.<ref>{{cite book| doi=10.1109/ACSSC.2003.1292037| chapter = Logical effort of carry propagate adders| title=The Thirty-Seventh Asilomar Conference on Signals, Systems & Computers, 2003 | pages=873–878| year=2003| last1=Harris|first1=D.| last2=Sutherland|first2=I.| isbn=0-7803-8104-1| s2cid=7880203}}</ref>


===Thermal time constant===
===Thermal time constant===
Line 129: Line 132:
Time constants are a feature of the [[lumped system analysis]] (lumped capacity analysis method) for thermal systems, used when objects cool or warm uniformly under the influence of [[Convection (heat transfer)|convective cooling or warming]]. In this case, the heat transfer from the body to the ambient at a given time is proportional to the temperature difference between the body and the ambient:<ref name=Lewis>{{cite book |title=Fundamentals of the finite element method for heat and fluid flow |author1=Roland Wynne Lewis |author2=Perumal Nithiarasu |author3=K. N. Seetharamu |url=https://books.google.com/books?id=dzCeaWOSjRkC&pg=PA151 |page=151 |isbn=978-0-470-84789-3 |year=2004 |publisher=Wiley}}</ref>
Time constants are a feature of the [[lumped system analysis]] (lumped capacity analysis method) for thermal systems, used when objects cool or warm uniformly under the influence of [[Convection (heat transfer)|convective cooling or warming]]. In this case, the heat transfer from the body to the ambient at a given time is proportional to the temperature difference between the body and the ambient:<ref name=Lewis>{{cite book |title=Fundamentals of the finite element method for heat and fluid flow |author1=Roland Wynne Lewis |author2=Perumal Nithiarasu |author3=K. N. Seetharamu |url=https://books.google.com/books?id=dzCeaWOSjRkC&pg=PA151 |page=151 |isbn=978-0-470-84789-3 |year=2004 |publisher=Wiley}}</ref>


:<math> F = hA_s \left( T(t) - T_a\right ), </math>
<math display="block"> F = hA_s \left( T(t) - T_a\right ), </math>


where ''h'' is the [[heat transfer coefficient]], and ''A''<sub>s</sub> is the surface area, ''T'' is the temperature function, i.e., ''T''(''t'') is the body temperature at time ''t'', and ''T''<sub>a</sub> is the constant ambient temperature. The positive sign indicates the convention that ''F'' is positive when heat is ''leaving'' the body because its temperature is higher than the ambient temperature (''F'' is an outward flux). If heat is lost to the ambient, this heat transfer leads to a drop in temperature of the body given by:<ref name=Lewis/>
where ''h'' is the [[heat transfer coefficient]], and ''A''<sub>s</sub> is the surface area, ''T'' is the temperature function, i.e., ''T''(''t'') is the body temperature at time ''t'', and ''T''<sub>a</sub> is the constant ambient temperature. The positive sign indicates the convention that ''F'' is positive when heat is ''leaving'' the body because its temperature is higher than the ambient temperature (''F'' is an outward flux). If heat is lost to the ambient, this heat transfer leads to a drop in temperature of the body given by:<ref name=Lewis/>


:<math> \rho c_p V \frac {dT}{dt} = -F, </math>
<math display="block"> \rho c_p V \frac {dT}{dt} = -F, </math>


where ''ρ'' = density, ''c''<sub>p</sub> = [[specific heat]] and ''V'' is the body volume. The negative sign indicates the temperature drops when the heat transfer is ''outward'' from the body (that is, when ''F'' > 0). Equating these two expressions for the heat transfer,
where ''ρ'' = density, ''c''<sub>p</sub> = [[specific heat]] and ''V'' is the body volume. The negative sign indicates the temperature drops when the heat transfer is ''outward'' from the body (that is, when ''F'' > 0). Equating these two expressions for the heat transfer,


:<math> \rho c_p V \frac {dT}{dt} = -hA_s \left( T(t) - T_a \right). </math>
<math display="block"> \rho c_p V \frac {dT}{dt} = -hA_s \left( T(t) - T_a \right). </math>


Evidently, this is a first-order LTI system that can be cast in the form:
Evidently, this is a first-order LTI system that can be cast in the form:


:<math>\frac {dT}{dt} +\frac {1}{\tau} T = \frac{1}{\tau} T_a, </math>
<math display="block">\frac {dT}{dt} +\frac {1}{\tau} T = \frac{1}{\tau} T_a, </math>


with
with


:<math>\tau = \frac{\rho c_p V}{hA_s}.</math>
<math display="block">\tau = \frac{\rho c_p V}{hA_s}.</math>


In other words, larger masses ''ρV'' with higher heat capacities ''c''<sub>p</sub> lead to slower changes in temperature (longer time constant ''τ''), while larger surface areas ''A''<sub>s</sub> with higher heat transfer ''h'' lead to more rapid temperature change (shorter time constant ''τ'').
In other words, larger masses ''ρV'' with higher heat capacities ''c''<sub>p</sub> lead to slower changes in temperature (longer time constant ''τ''), while larger surface areas ''A''<sub>s</sub> with higher heat transfer ''h'' lead to more rapid temperature change (shorter time constant ''τ'').
Line 151: Line 154:
Comparison with the introductory [[#Differential_equation|differential equation]] suggests the possible generalization to time-varying ambient temperatures ''T''<sub>a</sub>. However, retaining the simple constant ambient example, by substituting the variable Δ''T'' ≡ (''T − T''<sub>a</sub>), one finds:
Comparison with the introductory [[#Differential_equation|differential equation]] suggests the possible generalization to time-varying ambient temperatures ''T''<sub>a</sub>. However, retaining the simple constant ambient example, by substituting the variable Δ''T'' ≡ (''T − T''<sub>a</sub>), one finds:


:<math>\frac {d\Delta T}{dt} +\frac {1}{\tau} \Delta T = 0. </math>
<math display="block">\frac {d\Delta T}{dt} +\frac {1}{\tau} \Delta T = 0. </math>


Systems for which cooling satisfies the above exponential equation are said to satisfy [[Newton's law of cooling]]. The solution to this equation suggests that, in such systems, the difference between the temperature of the system and its surroundings Δ''T'' as a function of time ''t'', is given by:
Systems for which cooling satisfies the above exponential equation are said to satisfy [[Newton's law of cooling]]. The solution to this equation suggests that, in such systems, the difference between the temperature of the system and its surroundings Δ''T'' as a function of time ''t'', is given by:
Line 159: Line 162:
where Δ''T''<sub>0</sub> is the initial temperature difference, at time ''t'' = 0. In words, the body assumes the same temperature as the ambient at an exponentially slow rate determined by the time constant.
where Δ''T''<sub>0</sub> is the initial temperature difference, at time ''t'' = 0. In words, the body assumes the same temperature as the ambient at an exponentially slow rate determined by the time constant.


===Time constants in neuroscience===
===Time constants in Biophysics===
In an excitable cell such as a [[Myocyte|muscle]] or [[neuron]], the time constant of the [[membrane potential]] <math>\tau</math> is
In an excitable cell such as a [[Myocyte|muscle]] or [[neuron]], the time constant of the [[membrane potential]] <math>\tau</math> is
:<math>\tau = r_m c_m</math>
<math display="block">\tau = r_m c_m</math>


where ''r''<sub>m</sub> is the resistance across the membrane and ''c''<sub>m</sub> is the [[capacitance]] of the membrane.
where ''r''<sub>m</sub> is the resistance across the membrane and ''c''<sub>m</sub> is the [[capacitance]] of the membrane.
Line 168: Line 171:


The time constant is used to describe the rise and fall of membrane voltage, where the rise is described by
The time constant is used to describe the rise and fall of membrane voltage, where the rise is described by
:<math> V(t) = V_\textrm{max} \left(1 - e^{-t /\tau}\right) </math>
<math display="block"> V(t) = V_\textrm{max} \left(1 - e^{-t /\tau}\right) </math>


and the fall is described by
and the fall is described by
:<math> V(t) = V_\textrm{max} e^{-t /\tau} </math>
<math display="block"> V(t) = V_\textrm{max} e^{-t /\tau} </math>


where [[voltage]] is in millivolts, time is in seconds, and <math>\tau</math> is in seconds.
where [[voltage]] is in millivolts, time is in seconds, and <math>\tau</math> is in seconds.


''V''<sub>max</sub> is defined as the maximum voltage change from the [[resting potential]], where
''V''<sub>max</sub> is defined as the maximum voltage change from the [[resting potential]], where
:<math> V_\textrm{max} = r_m I </math>
<math display="block"> V_\textrm{max} = r_m I </math>


where ''r''<sub>m</sub> is the resistance across the membrane and ''I'' is the membrane current.
where ''r''<sub>m</sub> is the resistance across the membrane and ''I'' is the membrane current.
Line 188: Line 191:
===Exponential decay===
===Exponential decay===
{{Further|Exponential decay}}
{{Further|Exponential decay}}
In [[exponential decay]], such as of a [[radioactive decay|radioactive]] isotope, the time constant can be interpreted as the [[mean lifetime]]. The [[half-life]] ''T''<sub>''HL''</sub> is related to the exponential time constant <math>\tau</math> by
In [[exponential decay]], such as of a [[radioactive decay|radioactive]] isotope, the time constant can be interpreted as the [[mean lifetime]]. The [[half-life]] ''T''<sub>HL</sub> or ''T''<sub>1/2</sub> is related to the [[exponential decay constant]] <math>\tau</math> by
:<math> T_{HL} = \tau \ln 2. </math>
<math display="block"> T_{1/2} = T_\text{HL} = \tau \ln 2. </math>
The reciprocal of the time constant is called the [[decay constant]], and is denoted {{nowrap|<math>\lambda = 1/\tau</math>.}}
The reciprocal of the time constant is called the [[decay constant]], and is denoted {{nowrap|<math>\lambda = 1/\tau</math>.}}


Line 198: Line 201:


==See also==
==See also==
{{div col}}
*[[RC time constant]]
*[[RC time constant]]
*[[Cutoff frequency]]
*[[Cutoff frequency]]
*[[Exponential decay]]
*[[Exponential decay]]
*[[Lead-lag compensator]]
*[[Lead–lag compensator]]
*[[Length constant]]
*[[Length constant]]
*[[Rise time]]
*[[Rise time]]
Line 208: Line 212:
*[[Impulse response]]
*[[Impulse response]]
*[[Step response]]
*[[Step response]]
*[[Transition time]]
*[[Settling time]]
*[[Settling time]]
{{div col end}}


==Notes==
==Notes==
Line 227: Line 231:
[[Category:Physical constants]]
[[Category:Physical constants]]
[[Category:Neuroscience]]
[[Category:Neuroscience]]
[[Category:Durations]]

Latest revision as of 03:17, 22 June 2024

In physics and engineering, the time constant, usually denoted by the Greek letter τ (tau), is the parameter characterizing the response to a step input of a first-order, linear time-invariant (LTI) system.[1][note 1] The time constant is the main characteristic unit of a first-order LTI system. It gives speed of the response.

In the time domain, the usual choice to explore the time response is through the step response to a step input, or the impulse response to a Dirac delta function input.[2] In the frequency domain (for example, looking at the Fourier transform of the step response, or using an input that is a simple sinusoidal function of time) the time constant also determines the bandwidth of a first-order time-invariant system, that is, the frequency at which the output signal power drops to half the value it has at low frequencies.

The time constant is also used to characterize the frequency response of various signal processing systems – magnetic tapes, radio transmitters and receivers, record cutting and replay equipment, and digital filters – which can be modelled or approximated by first-order LTI systems. Other examples include time constant used in control systems for integral and derivative action controllers, which are often pneumatic, rather than electrical.

Time constants are a feature of the lumped system analysis (lumped capacity analysis method) for thermal systems, used when objects cool or warm uniformly under the influence of convective cooling or warming.[3]

Physically, the time constant represents the elapsed time required for the system response to decay to zero if the system had continued to decay at the initial rate, because of the progressive change in the rate of decay the response will have actually decreased in value to 1 / e ≈ 36.8% in this time (say from a step decrease). In an increasing system, the time constant is the time for the system's step response to reach 1 − 1 / e ≈ 63.2% of its final (asymptotic) value (say from a step increase). In radioactive decay the time constant is related to the decay constant (λ), and it represents both the mean lifetime of a decaying system (such as an atom) before it decays, or the time it takes for all but 36.8% of the atoms to decay. For this reason, the time constant is longer than the half-life, which is the time for only 50% of the atoms to decay.

Differential equation

[edit]

First order LTI systems are characterized by the differential equation

where τ represents the exponential decay constant and V is a function of time t The right-hand side is the forcing function f(t) describing an external driving function of time, which can be regarded as the system input, to which V(t) is the response, or system output. Classical examples for f(t) are:

The Heaviside step function, often denoted by u(t): the impulse function, often denoted by δ(t), and also the sinusoidal input function: or where A is the amplitude of the forcing function, f is the frequency in Hertz, and ω = 2π f is the frequency in radians per second.

Example solution

[edit]

An example solution to the differential equation with initial value V0 and no forcing function is

where

is the initial value of V. Thus, the response is an exponential decay with time constant τ.

Discussion

[edit]

Suppose

This behavior is referred to as a "decaying" exponential function. The time τ (tau) is referred to as the "time constant" and can be used (as in this case) to indicate how rapidly an exponential function decays.

Here:

  • t is time (generally t > 0 in control engineering)
  • V0 is the initial value (see "specific cases" below).

Specific cases

[edit]
  1. Let ; then , and so
  2. Let ; then
  3. Let , and so
  4. Let ; then

After a period of one time constant the function reaches e−1 = approximately 37% of its initial value. In case 4, after five time constants the function reaches a value less than 1% of its original. In most cases this 1% threshold is considered sufficient to assume that the function has decayed to zero – as a rule of thumb, in control engineering a stable system is one that exhibits such an overall damped behavior.

Relation of time constant to bandwidth

[edit]
An example response of system to sine wave forcing function. Time axis in units of the time constant τ. The response damps out to become a simple sine wave.
Frequency response of system vs. frequency in units of the bandwidth f3dB. The response is normalized to a zero frequency value of unity, and drops to 1/√2 at the bandwidth.

Suppose the forcing function is chosen as sinusoidal so:

(Response to a real cosine or sine wave input can be obtained by taking the real or imaginary part of the final result by virtue of Euler's formula.) The general solution to this equation for times t ≥ 0 s, assuming V(t = 0) = V0 is:

For long times the decaying exponentials become negligible and the steady-state solution or long-time solution is:

The magnitude of this response is: By convention, the bandwidth of this system is the frequency where |V|2 drops to half-value, or where ωτ = 1. This is the usual bandwidth convention, defined as the frequency range where power drops by less than half (at most −3 dB). Using the frequency in hertz, rather than radians/s (ω = 2πf):

The notation f3dB stems from the expression of power in decibels and the observation that half-power corresponds to a drop in the value of |V| by a factor of 1/2 or by 3 decibels.

Thus, the time constant determines the bandwidth of this system.

Step response with arbitrary initial conditions

[edit]
Step response of system for two different initial values V0, one above the final value and one at zero. Long-time response is a constant, V. Time axis in units of the time constant .

Suppose the forcing function is chosen as a step input so:

with u(t) the unit step function. The general solution to this equation for times t ≥ 0 s, assuming V(t = 0) = V0 is:

(It may be observed that this response is the ω → 0 limit of the above response to a sinusoidal input.)

The long-time solution is time independent and independent of initial conditions:

The time constant remains the same for the same system regardless of the starting conditions. Simply stated, a system approaches its final, steady-state situation at a constant rate, regardless of how close it is to that value at any arbitrary starting point.

For example, consider an electric motor whose startup is well modelled by a first-order LTI system. Suppose that when started from rest, the motor takes 1/8 of a second to reach 63% of its nominal speed of 100 RPM, or 63 RPM—a shortfall of 37 RPM. Then it will be found that after the next 1/8 of a second, the motor has sped up an additional 23 RPM, which equals 63% of that 37 RPM difference. This brings it to 86 RPM—still 14 RPM low. After a third 1/8 of a second, the motor will have gained an additional 9 RPM (63% of that 14 RPM difference), putting it at 95 RPM.

In fact, given any initial speed s ≤ 100 RPM, 1/8 of a second later this particular motor will have gained an additional 0.63 × (100 − s) RPM.

Examples

[edit]

Time constants in electrical circuits

[edit]
Capacitor voltage step-response.
Inductor voltage step-response.

In an RL circuit composed of a single resistor and inductor, the time constant (in seconds) is

where R is the resistance (in ohms) and L is the inductance (in henrys).

Similarly, in an RC circuit composed of a single resistor and capacitor, the time constant (in seconds) is:

where R is the resistance (in ohms) and C is the capacitance (in farads).

Electrical circuits are often more complex than these examples, and may exhibit multiple time constants (See Step response and Pole splitting for some examples.) In the case where feedback is present, a system may exhibit unstable, increasing oscillations. In addition, physical electrical circuits are seldom truly linear systems except for very low amplitude excitations; however, the approximation of linearity is widely used.

In digital electronic circuits another measure, the FO4 is often used. This can be converted to time constant units via the equation .[4]

Thermal time constant

[edit]

Time constants are a feature of the lumped system analysis (lumped capacity analysis method) for thermal systems, used when objects cool or warm uniformly under the influence of convective cooling or warming. In this case, the heat transfer from the body to the ambient at a given time is proportional to the temperature difference between the body and the ambient:[5]

where h is the heat transfer coefficient, and As is the surface area, T is the temperature function, i.e., T(t) is the body temperature at time t, and Ta is the constant ambient temperature. The positive sign indicates the convention that F is positive when heat is leaving the body because its temperature is higher than the ambient temperature (F is an outward flux). If heat is lost to the ambient, this heat transfer leads to a drop in temperature of the body given by:[5]

where ρ = density, cp = specific heat and V is the body volume. The negative sign indicates the temperature drops when the heat transfer is outward from the body (that is, when F > 0). Equating these two expressions for the heat transfer,

Evidently, this is a first-order LTI system that can be cast in the form:

with

In other words, larger masses ρV with higher heat capacities cp lead to slower changes in temperature (longer time constant τ), while larger surface areas As with higher heat transfer h lead to more rapid temperature change (shorter time constant τ).

Comparison with the introductory differential equation suggests the possible generalization to time-varying ambient temperatures Ta. However, retaining the simple constant ambient example, by substituting the variable ΔT ≡ (T − Ta), one finds:

Systems for which cooling satisfies the above exponential equation are said to satisfy Newton's law of cooling. The solution to this equation suggests that, in such systems, the difference between the temperature of the system and its surroundings ΔT as a function of time t, is given by:

where ΔT0 is the initial temperature difference, at time t = 0. In words, the body assumes the same temperature as the ambient at an exponentially slow rate determined by the time constant.

Time constants in Biophysics

[edit]

In an excitable cell such as a muscle or neuron, the time constant of the membrane potential is

where rm is the resistance across the membrane and cm is the capacitance of the membrane.

The resistance across the membrane is a function of the number of open ion channels and the capacitance is a function of the properties of the lipid bilayer.

The time constant is used to describe the rise and fall of membrane voltage, where the rise is described by

and the fall is described by

where voltage is in millivolts, time is in seconds, and is in seconds.

Vmax is defined as the maximum voltage change from the resting potential, where

where rm is the resistance across the membrane and I is the membrane current.

Setting for t = for the rise sets V(t) equal to 0.63Vmax. This means that the time constant is the time elapsed after 63% of Vmax has been reached

Setting for t = for the fall sets V(t) equal to 0.37Vmax, meaning that the time constant is the time elapsed after it has fallen to 37% of Vmax.

The larger a time constant is, the slower the rise or fall of the potential of a neuron. A long time constant can result in temporal summation, or the algebraic summation of repeated potentials. A short time constant rather produces a coincidence detector through spatial summation.

Exponential decay

[edit]

In exponential decay, such as of a radioactive isotope, the time constant can be interpreted as the mean lifetime. The half-life THL or T1/2 is related to the exponential decay constant by The reciprocal of the time constant is called the decay constant, and is denoted .

Meteorological sensors

[edit]

A time constant is the amount of time it takes for a meteorological sensor to respond to a rapid change in a measure, and until it is measuring values within the accuracy tolerance usually expected of the sensor.

This most often applies to measurements of temperature, dew-point temperature, humidity and air pressure. Radiosondes are especially affected due to their rapid increase in altitude.

See also

[edit]

Notes

[edit]
  1. ^ Concretely, a first-order LTI system is a system that can be modeled by a single first order differential equation in time. Examples include the simplest single-stage electrical RC circuits and RL circuits.

References

[edit]
  1. ^ Béla G. Lipták (2003). Instrument Engineers' Handbook: Process control and optimization (4 ed.). CRC Press. p. 100. ISBN 978-0-8493-1081-2.
  2. ^ Bong Wie (1998). Space vehicle dynamics and control. American Institute of Aeronautics and Astronautics. p. 100. ISBN 978-1-56347-261-9.
  3. ^ GR North (1988). "Lessons from energy balance models". In Michael E. Schlesinger (ed.). Physically-based Modelling and Simulation of Climate and Climatic Change (NATO Advanced Study Institute on Physical-Based Modelling ed.). Springer. p. 627. ISBN 978-90-277-2789-3. {{cite book}}: Unknown parameter |agency= ignored (help)
  4. ^ Harris, D.; Sutherland, I. (2003). "Logical effort of carry propagate adders". The Thirty-Seventh Asilomar Conference on Signals, Systems & Computers, 2003. pp. 873–878. doi:10.1109/ACSSC.2003.1292037. ISBN 0-7803-8104-1. S2CID 7880203.
  5. ^ a b Roland Wynne Lewis; Perumal Nithiarasu; K. N. Seetharamu (2004). Fundamentals of the finite element method for heat and fluid flow. Wiley. p. 151. ISBN 978-0-470-84789-3.
[edit]