Time constant

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.

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 \approx 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 \approx 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

First order LTI systems are characterized by the differential equation

\tau {\frac {dV}{dt}}+V=f(t)

where $τ$ represents the exponential decay constant and $V$ is a function of time $t$

V=V(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)$ :

u(t)={\begin{cases}0,&t<0\\1,&t\geq 0\end{cases}}

the impulse function, often denoted by $δ (t)$ , and also the sinusoidal input function:

f(t)=A\sin(2\pi ft)

or

f(t)=Ae^{j\omega t},

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

An example solution to the differential equation with initial value $V 0$ and no forcing function is

V(t)=V_{0}e^{-t/\tau }

where

V_{0}=V(t=0)

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

Discussion

Suppose

V(t)=V_{0}e^{-t/\tau }.

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)
$V 0$ is the initial value (see "specific cases" below).

Specific cases

Let $t=0$ ; then $V=V_{0}e^{0}$ , and so $V=V_{0}$
Let $t=\tau$ ; then $V=V_{0}e^{-1}\approx 0.37V_{0}$
Let $V=f(t)=V_{0}e^{-t/\tau }$ , and so ${\textstyle \lim _{t\to \infty }f(t)=0}$
Let $t=5\tau$ ; then $V=V_{0}e^{-5}\approx 0.0067V_{0}$

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

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

f 3dB

. 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:

\tau {\frac {dV}{dt}}+V=f(t)=Ae^{j\omega t}.

(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 \geq 0 s$ , assuming $V (t = 0) = V 0$ is:

{\begin{aligned}V(t)&=V_{0}e^{-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 {\frac {1}{\tau }}{j\omega +{\frac {1}{\tau }}}}A\left(e^{j\omega t}-e^{-t/\tau }\right).\end{aligned}}

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

V_{\infty }(t)={\frac {1/\tau }{j\omega +1/\tau }}Ae^{j\omega t}.

The magnitude of this response is:

|V_{\infty }(t)|=A{\frac {1}{\tau \left(\omega ^{2}+(1/\tau )^{2}\right)^{1/2}}}=A{\frac {1}{\sqrt {1+(\omega \tau )^{2}}}}.

By convention, the bandwidth of this system is the frequency where $| V \infty | 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$ ):

f_{\mathrm {3dB} }={\frac {1}{2\pi \tau }}.

The notation $f 3dB$ stems from the expression of power in decibels and the observation that half-power corresponds to a drop in the value of $| V \infty |$ 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

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

\tau

.

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

{\frac {dV}{dt}}+{\frac {1}{\tau }}V=f(t)=Au(t),

with $u (t)$ the Heaviside step function. The general solution to this equation for times $t \geq 0 s$ , assuming $V (t = 0) = V 0$ is:

V(t)=V_{0}e^{-t/\tau }+A\tau \left(1-e^{-t/\tau }\right).

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

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

V_{\infty }=A\tau .

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

Time constants in electrical circuits

Capacitor voltage step-response.

Inductor voltage step-response.

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

\tau ={\frac {L}{R}}

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 $\tau$ (in seconds) is:

\tau =RC

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 $5\tau ={\text{FO4}}$ .^[4]

Thermal time constant

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]

F=hA_{s}\left(T(t)-T_{a}\right),

where h is the heat transfer coefficient, and A_s is the surface area, T is the temperature function, i.e., T(t) is the body temperature at time t, and T_a 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]

\rho c_{p}V{\frac {dT}{dt}}=-F,

where ρ = density, c_p = 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,

\rho c_{p}V{\frac {dT}{dt}}=-hA_{s}\left(T(t)-T_{a}\right).

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

{\frac {dT}{dt}}+{\frac {1}{\tau }}T={\frac {1}{\tau }}T_{a},

with

\tau ={\frac {\rho c_{p}V}{hA_{s}}}.

In other words, larger masses ρV with higher heat capacities c_p lead to slower changes in temperature (longer time constant τ), while larger surface areas A_s 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 T_a. However, retaining the simple constant ambient example, by substituting the variable ΔT ≡ (T − T_a), one finds:

{\frac {d\Delta T}{dt}}+{\frac {1}{\tau }}\Delta T=0.

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:

\Delta T(t)=\Delta T_{0}e^{-t/\tau },

where ΔT₀ 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

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

\tau =r_{m}c_{m}

where r_m is the resistance across the membrane and c_m 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

V(t)=V_{\textrm {max}}\left(1-e^{-t/\tau }\right)

and the fall is described by

V(t)=V_{\textrm {max}}e^{-t/\tau }

where voltage is in millivolts, time is in seconds, and $\tau$ is in seconds.

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

V_{\textrm {max}}=r_{m}I

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

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

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

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

In exponential decay, such as of a radioactive isotope, the time constant can be interpreted as the mean lifetime. The half-life T_HL or T_1/2 is related to the exponential decay constant $\tau$ by

T_{1/2}=T_{\text{HL}}=\tau \ln 2.

The reciprocal of the time constant is called the decay constant, and is denoted $\lambda =1/\tau$ .

Meteorological sensors

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.

Notes

↑ 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

↑ 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.
↑ Bong Wie (1998). Space vehicle dynamics and control. American Institute of Aeronautics and Astronautics. p. 100. ISBN 978-1-56347-261-9.
↑ 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. NATO. p. 627. ISBN 978-90-277-2789-3.
↑ 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.
1 2 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.

External links

This article is issued from Wikipedia. The text is licensed under Creative Commons - Attribution - Sharealike. Additional terms may apply for the media files.

[2] 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.

[Lipták-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.

[Wie-3] Bong Wie (1998). Space vehicle dynamics and control. American Institute of Aeronautics and Astronautics. p. 100. ISBN 978-1-56347-261-9.

[NATO-4] 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. NATO. p. 627. ISBN 978-90-277-2789-3.

[5] 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.

[Lewis-6] 1 2 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.

[1]

[note 1]

[2]

[3]

[4]

[5]