Lumped-element model: Difference between revisions

Content deleted Content added
Link suggestions feature: 3 links added.
 
(29 intermediate revisions by 21 users not shown)
Line 1:
{{short description|Simplification of a physical system into a network of discrete components}}
{{Technical|date=August 2019}}
{{Refimprove|date=August 2023}}
[[File:Ohm's Law with Voltage source TeX.svg|right|thumb|Representation of a lumped model made upconsisting of a voltage source and a resistor.]]
 
The '''lumped-element model''' (also called '''lumped-parameter model''', or '''lumped-component model''') simplifiesis thea description[[idealization (philosophy of thescience)|simplified]] behaviourrepresentation of spatially distributeda [[physical systemssystem]] intoor circuit that assumes all components are concentrated at a [[Topologysingle (electricalpoint circuits)|topology]]and consistingtheir ofbehavior discretecan entitiesbe thatdescribed approximateby theidealized behaviourmathematical ofmodels. theThe distributedlumped-element model simplifies the system underor certaincircuit assumptionsbehavior description into a [[Topology (electrical circuits)|topology]]. It is useful in [[electrical network|electrical systems]] (including [[electronics]]), mechanical [[multibody system]]s, [[heat transfer]], [[acoustics]], etc. This is in contrast to [[distributed parameter system]]s or models in which the behaviour is distributed spatially and cannot be considered as localized into discrete entities.
 
Mathematically speaking, theThe simplification reduces the [[State space (controls)|state space]] of the system to a [[counting number|finite]] [[dimension]], and the [[partial differential equation]]s (PDEs) of the continuous (infinite-dimensional) time and space model of the physical system into [[ordinary differential equation]]s (ODEs) with a finite number of parameters.
 
== Electrical systems ==
<!-- started as a copy from https[[Special:PermanentLink//en.wikipedia.org/w/index.php?title=Lumped_matter_discipline&oldid=544612541]] -->
 
=== Lumped-matter discipline ===
The '''lumped-matter discipline''' is a set of imposed assumptions in [[electrical engineering]] that provides the foundation for '''lumped-circuit abstraction''' used in [[Network analysis (electrical circuits)|network analysis]].<ref>Anant Agarwal and Jeffrey Lang, course materials for 6.002 Circuits and Electronics, Spring 2007. MIT OpenCourseWare ([http://ocw.mit.edu/courses/electrical-engineering-and-computer-science/6-002-circuits-and-electronics-spring-2007/video-lectures/6002_l1.pdf PDF]), [[Massachusetts Institute of Technology]].</ref> The self-imposed constraints are:
# The change of the magnetic flux in time outside a conductor is zero. <math display="block">\frac{\partial \Phi_B} {\partial t} = 0</math>
 
1.# The change of the magnetic fluxcharge in time outsideinside aconducting conductorelements is zero. <math display="block">\frac{\partial q} {\partial t} = 0</math>
3.# Signal timescales of interest are much larger than [[propagation delay]] of [[electromagnetic waves]] across the lumped element.
::<math>\frac{\partial \phi_B} {\partial t} = 0</math>
 
2. The change of the charge in time inside conducting elements is zero.
::<math>\frac{\partial q} {\partial t} = 0</math>
 
3. Signal timescales of interest are much larger than propagation delay of [[electromagnetic waves]] across the lumped element.
 
The first two assumptions result in [[Kirchhoff's circuit laws]] when applied to [[Maxwell's equations]] and are only applicable when the circuit is in [[steady state (electronics)|steady state]]. The third assumption is the basis of the lumped-element model used in [[Network analysis (electrical circuits)|network analysis]]. Less severe assumptions result in the [[distributed-element model]], while still not requiring the direct application of the full Maxwell equations.
 
=== Lumped-element model ===
The lumped-element model of electronic [[Electrical network|circuitcircuits]]s makes the simplifying assumption that the attributes of the circuit, [[Electrical resistance|resistance]], [[capacitance]], [[inductance]], and [[Gain (electronics)|gain]], are concentrated into idealized [[electrical component]]s; [[resistor]]s, [[capacitor]]s, and [[inductor]]s, etc. joined by a network of perfectly [[Electrical conduction|conducting]] wires.
 
The lumped-element model is valid whenever <math>L_c \ll \lambda</math>, where <math>L_c</math> denotes the circuit's characteristic length, and <math>\lambda</math> denotes the circuit's operating [[wavelength]]. Otherwise, when the circuit length is on the order of a wavelength, we must consider more general models, such as the [[distributed-element model]] (including [[transmission line]]s), whose dynamic behaviour is described by [[Maxwell's equations]]. Another way of viewing the validity of the lumped-element model is to note that this model ignores the finite time it takes signals to propagate around a circuit. Whenever this propagation time is not significant to the application the lumped-element model can be used. This is the case when the propagation time is much less than the [[period (physics)|period]] of the signal involved. However, with increasing propagation time there will be an increasing error between the assumed and actual phase of the signal which in turn results in an error in the assumed amplitude of the signal. The exact point at which the lumped-element model can no longer be used depends to a certain extent on how accurately the signal needs to be known in a given application.
Line 29 ⟶ 26:
Real-world components exhibit non-ideal characteristics which are, in reality, distributed elements but are often represented to a [[first-order approximation]] by lumped elements. To account for leakage in [[capacitor]]s for example, we can model the non-ideal capacitor as having a large lumped [[resistor]] connected in parallel even though the leakage is, in reality distributed throughout the dielectric. Similarly a [[wire-wound resistor]] has significant [[inductance]] as well as [[Electrical resistance|resistance]] distributed along its length but we can model this as a lumped [[inductor]] in series with the ideal resistor.
 
== Thermal systems ==
A '''lumped-capacitance model''', also called '''lumped system analysis''',<ref>{{cite book | last = Incropera |author2=DeWitt |author3=Bergman |author4=Lavine | title = Fundamentals of Heat and Mass Transfer | url = https://archive.org/details/fundamentalsheat00incr_869 | url-access = limited | edition = 6th | year = 2007 | isbn = 978-0-471-45728-2 | publisher = John Wiley & Sons | pages = 260–261[https://archive.org/details/fundamentalsheat00incr_869/page/n267 260]–261}}</ref> reduces a [[thermal system]] to a number of discrete “lumps” and assumes that the [[temperature]] difference inside each lump is negligible. This approximation is useful to simplify otherwise complex [[differential equation|differential]] heat equations. It was developed as a mathematical analog of [[electrical capacitance]], although it also includes thermal analogs of [[electrical resistance]] as well.
 
The lumped-capacitance model is a common approximation in transient conduction, which may be used whenever [[heat conduction]] within an object is much faster than heat transfer across the boundary of the object. The method of approximation then suitably reduces one aspect of the transient conduction system (spatial temperature variation within the object) to a more mathematically tractable form (that is, it is assumed that the temperature within the object is completely uniform in space, although this spatially uniform temperature value changes over time). The rising uniform temperature within the object or part of a system, can then be treated like a capacitative reservoir which absorbs heat until it reaches a steady thermal state in time (after which temperature does not change within it).
Line 36 ⟶ 33:
An early-discovered example of a lumped-capacitance system which exhibits mathematically simple behavior due to such physical simplifications, are systems which conform to ''Newton's law of cooling''. This law simply states that the temperature of a hot (or cold) object progresses toward the temperature of its environment in a simple exponential fashion. Objects follow this law strictly only if the rate of heat conduction within them is much larger than the heat flow into or out of them. In such cases it makes sense to talk of a single "object temperature" at any given time (since there is no spatial temperature variation within the object) and also the uniform temperatures within the object allow its total thermal energy excess or deficit to vary proportionally to its surface temperature, thus setting up the Newton's law of cooling requirement that the rate of temperature decrease is proportional to difference between the object and the environment. This in turn leads to simple exponential heating or cooling behavior (details below).
 
=== Method ===
To determine the number of lumps, the [[Biot number]] (Bi), a dimensionless parameter of the system, is used. Bi is defined as the ratio of the conductive heat resistance within the object to the [[Convection (heat transfer)|convective heat transfer]] resistance across the object's boundary with a uniform bath of different temperature. When the [[thermal resistance]] to heat transferred into the object is larger than the resistance to heat being [[diffused]] completely within the object, the Biot number is less than 1. In this case, particularly for Biot numbers which are even smaller, the approximation of ''spatially uniform temperature within the object'' can begin to be used, since it can be presumed that heat transferred into the object has time to uniformly distribute itself, due to the lower resistance to doing so, as compared with the resistance to heat entering the object.
 
If the Biot number is less than 0.1 for a solid object, then the entire material will be nearly the same temperature, with the dominant temperature difference will bebeing at the surface. It may be regarded as being "thermally thin". The Biot number must generally be less than 0.1 for usefully accurate approximation and heat transfer analysis. The mathematical solution to the lumped-system approximation gives [[Newton's law of cooling]].
 
A Biot number greater than 0.1 (a "thermally thick" substance) indicates that one cannot make this assumption, and more complicated [[heat transfer]] equations for "transient heat conduction" will be required to describe the time-varying and non-spatially-uniform temperature field within the material body.
Line 59 ⟶ 56:
The equations describing the three heat transfer modes and their thermal resistances in steady state conditions, as discussed previously, are summarized in the table below:
 
{| class="wikitable" style="margin:1em auto; text-align:center;"
<center>
{| class="wikitable" style="text-align:center;"
|+Equations for different heat transfer modes and their thermal resistances.
|-
Line 77 ⟶ 73:
|Radiation
|<math>\dot{Q}=\frac{T_{\rm surf}-T_{\rm surr}}{\left ( \frac{1}{h_rA_{\rm surf}} \right )}</math>
|<math>\frac{1}{h_rA}</math>, where<br /><math>h_r= \epsilon \sigma (T_{\rm surf}^{2}+T_{\rm surr}^{2})(T_{\rm surf}+T_{\rm surr})</math>
|-
|}
</center>
 
In cases where there is heat transfer through different media (for example, through a [[composite material]]), the equivalent resistance is the sum of the resistances of the components that make up the composite. Likely, in cases where there are different heat transfer modes, the total resistance is the sum of the resistances of the different modes. Using the thermal circuit concept, the amount of heat transferred through any medium is the quotient of the temperature change and the total thermal resistance of the medium.
 
As an example, consider a composite wall of cross-sectional area <math>A</math>. The composite is made of an <math>L_1</math> long cement plaster with a thermal coefficient <math>k_1</math> and <math>L_2</math> long paper faced fiber glass, with thermal coefficient <math>k_2</math>. The left surface of the wall is at <math>T_i</math> and exposed to air with a convective coefficient of <math>h_i</math>. The right surface of the wall is at <math>T_o</math> and exposed to air with convective coefficient <math>h_o</math>.
 
<!-- Deleted image removed: [[Image:Thermal Circuits2.jpg|thumb|center|500px| {{pufc|Thermal Circuits2.jpg}}]] -->
 
Using the thermal resistance concept, heat flow through the composite is as follows:
<math display="block">\dot{Q}=\frac{T_i-T_o}{R_i+R_1+R_2+R_o}=\frac{T_i-T_1}{R_i}=\frac{T_i-T_2}{R_i+R_1}=\frac{T_i-T_3}{R_i+R_1+R_2}=\frac{T_1-T_2}{R_1}=\frac{T_3-T_o}{R_0}</math>
where <math>R_i=\frac{1}{h_iA}</math>, <math>R_o=\frac{1}{h_oA}</math>, <math>R_1=\frac{L_1}{k_1A}</math>, and <math>R_2=\frac{L_2}{k_2A}</math>
 
==== Newton's law of cooling ====
<math>\dot{Q}=\frac{T_i-T_o}{R_i+R_1+R_2+R_o}=\frac{T_i-T_1}{R_i}=\frac{T_i-T_2}{R_i+R_1}=\frac{T_i-T_3}{R_i+R_1+R_2}=\frac{T_1-T_2}{R_1}=\frac{T_3-T_o}{R_0}</math>
{{Main|Newton's Lawlaw of Coolingcooling}}
 
'''Newton's law of cooling''' is an [[empirical relationship]] attributed to English physicist [[Isaac Newton|Sir Isaac Newton]] (1642 - 17271642–1727). This law stated in non-mathematical form is the following:
where
 
<math>R_i=\frac{1}{h_iA}</math>, <math>R_o=\frac{1}{h_oA}</math>, <math>R_1=\frac{L_1}{k_1A}</math>, and <math>R_2=\frac{L_2}{k_2A}</math>
 
====Newton's law of cooling====
{{Main|Newton's Law of Cooling}}
 
'''Newton's law of cooling''' is an empirical relationship attributed to English physicist [[Isaac Newton|Sir Isaac Newton]] (1642 - 1727). This law stated in non-mathematical form is the following:
 
{{Quotation|The rate of heat loss of a body is proportional to the temperature difference between the body and its surroundings.}}
 
Or, using symbols:
:<math display="block">\text {Rate of cooling} \sim \!\, \Delta T</math>
 
An object at a different temperature from its surroundings will ultimately come to a common temperature with its surroundings. A relatively hot object cools as it warms its surroundings; a cool object is warmed by its surroundings. When considering how quickly (or slowly) something cools, we speak of its ''rate'' of cooling - how many degrees' change in temperature per unit of time.
:<math>\text {Rate of cooling} \sim \!\, \Delta T</math>
 
An object at a different temperature from its surroundings will ultimately come to a common temperature with its surroundings. A relatively hot object cools as it warms its surroundings; a cool object is warmed by its surroundings. When considering how quickly (or slowly) something cools, we speak of its ''rate'' of cooling - how many degrees' change in temperature per unit of time.
 
The rate of cooling of an object depends on how much hotter the object is than its surroundings. The temperature change per minute of a hot apple pie will be more if the pie is put in a cold freezer than if it is placed on the kitchen table. When the pie cools in the freezer, the temperature difference between it and its surroundings is greater. On a cold day, a warm home will leak heat to the outside at a greater rate when there is a large difference between the inside and outside temperatures. Keeping the inside of a home at high temperature on a cold day is thus more costly than keeping it at a lower temperature. If the temperature difference is kept small, the rate of cooling will be correspondingly low.
 
As Newton's law of cooling states, the rate of cooling of an object - whether by [[thermal conduction|conduction]], [[ConvectiveConvection (heat transfer)|convection]], or [[thermal radiation|radiation]] - is approximately proportional to the temperature difference Δ''T''. Frozen food will warm up faster in a warm room than in a cold room. Note that the rate of cooling experienced on a cold day can be increased by the added convection effect of the [[wind]]. This is referred to as [[wind chill]]. For example, a wind chill of -20&nbsp;°C means that heat is being lost at the same rate as if the temperature were -20&nbsp;°C without wind.
 
==== Applicable situations ====
This law describes many situations in which an object has a large thermal capacity and large conductivity, and is suddenly immersed in a uniform bath which conducts heat relatively poorly. It is an example of a thermal circuit with one resistive and one capacitative element. For the law to be correct, the temperatures at all points inside the body must be approximately the same at each time point, including the temperature at its surface. Thus, the temperature difference between the body and surroundings does not depend on which part of the body is chosen, since all parts of the body have effectively the same temperature. In these situations, the material of the body does not act to "insulate" other parts of the body from heat flow, and all of the significant insulation (or "thermal resistance") controlling the rate of heat flow in the situation resides in the area of contact between the body and its surroundings. Across this boundary, the temperature-value jumps in a discontinuous fashion.
 
Line 120 ⟶ 107:
In such a situation, the object acts as the "capacitative" circuit element, and the resistance of the thermal contact at the boundary acts as the (single) thermal resistor. In electrical circuits, such a combination would charge or discharge toward the input voltage, according to a simple exponential law in time. In the thermal circuit, this configuration results in the same behavior in temperature: an exponential approach of the object temperature to the bath temperature.
 
==== Mathematical statement ====
Newton's law is mathematically stated by the simple first-order differential equation:
:<math display="block"> \frac{d Q}{d t} = - h \cdot A(T(t)- T_{\text{env}}) = - h \cdot A \Delta T(t)\quad </math>
 
where
:*''Q'' is thermal energy in [[joule]]s
:*''h'' is the [[heat transfer coefficient]] between the surface and the fluid
:*''A'' is the surface area of the heat being transferred
:*''T'' is the temperature of the object's surface and interior (since these are the same in this approximation)
:*''T''<sub>env</sub>'' is the temperature of the environment
:*Δ''T''(''t)'') = ''T''(''t)'') - ''T''<sub>env</sub>'' is the time-dependent thermal [[gradient]] between environment and object
 
Putting heat transfers into this form is sometimes not a very good approximation, depending on ratios of heat conductances in the system. If the differences are not large, an accurate formulation of heat transfers in the system may require analysis of heat flow based on the (transient) heat transfer equation in nonhomogeneous or poorly conductive media.
Line 138 ⟶ 124:
If the entire body is treated as lumped-capacitance heat reservoir, with total heat content which is proportional to simple total [[heat capacity]] <math>C</math>, and <math>T</math>, the temperature of the body, or <math>Q = C T</math>. It is expected that the system will experience [[exponential decay]] with time in the temperature of a body.
 
From the definition of heat capacity <math>C</math> comes the relation <math>C = dQ/dT</math>. Differentiating this equation with regard to time gives the identity (valid so long as temperatures in the object are uniform at any given time): <math>dQ/dt = C (dT/dt)</math>. This expression may be used to replace <math>dQ/dt</math> in the first equation which begins this section, above. Then, if <math>T(t)</math> is the temperature of such a body at time <math>t</math>, and <math>T_\text{env}</math> is the temperature of the environment around the body:
:<math display="block"> \frac{d T(t)}{d t} = - r (T(t) - T_{\mathrmtext{env}}) = - r \Delta T(t)\quad </math>
 
where <math>r = hA/C</math> is a positive constant characteristic of the system, which must be in units of <math>s^{-1}</math>, and is therefore sometimes expressed in terms of a characteristic [[time constant]] <math>t_0</math> given by: <math>rt_0 = 1/t_0r = -\Delta T(t)/(dT(t)/dt)</math>. Thus, in thermal systems, <math>t_0 = C/hA</math>. (The total [[heat capacity]] <math>C</math> of a system may be further represented by its mass-[[specific heat capacity]] <math>c_p</math> multiplied by its mass <math>m</math>, so that the time constant <math>t_0</math> is also given by <math>mc_p/hA</math>).
:<math> \frac{d T(t)}{d t} = - r (T(t) - T_{\mathrm{env}}) = - r \Delta T(t)\quad </math>
 
where
 
<math>r = hA/C</math> is a positive constant characteristic of the system, which must be in units of <math>s^{-1}</math>, and is therefore sometimes expressed in terms of a characteristic [[time constant]] <math>t_0</math> given by: <math>r = 1/t_0 = \Delta T/(dT(t)/dt)</math>. Thus, in thermal systems, <math>t_0 = C/hA</math>. (The total [[heat capacity]] <math>C</math> of a system may be further represented by its mass-[[specific heat capacity]] <math>c_p</math> multiplied by its mass <math>m</math>, so that the time constant <math>t_0</math> is also given by <math>mc_p/hA</math>).
 
The solution of this differential equation, by standard methods of integration and substitution of boundary conditions, gives:
: <math display="block"> T(t) = T_{\mathrm{env}} + (T(0) - T_{\mathrm{env}}) \ e^{-r t}. \quad </math>
 
: <math> T(t) = T_{\mathrm{env}} + (T(0) - T_{\mathrm{env}}) \ e^{-r t}. \quad </math>
 
If:
Line 155 ⟶ 136:
 
then the Newtonian solution is written as:
: <math display="block"> \Delta T(t) = \Delta T(0) \ e^{-r t} = \Delta T(0) \ e^{-t/t_0}. \quad </math>
 
: <math> \Delta T(t) = \Delta T(0) \ e^{-r t} = \Delta T(0) \ e^{-t/t_0}. \quad </math>
 
This same solution is almost immediately apparent if the initial differential equation is written in terms of <math>\Delta T(t)</math>, as the single function to be solved for.
:<math display="block"> \frac{d T(t)}{d t} = \frac{d \Delta T(t)}{d t} = - \frac{1}{t_0} \Delta T(t)\quad </math>
'
:<math> \frac{d T(t)}{d t} = \frac{d \Delta T(t)}{d t} = - \frac{1}{t_0} \Delta T(t)\quad </math>
 
=== Applications ===
This mode of analysis has been applied to [[forensic science]]s to analyze the time of death of humans. Also, it can be applied to [[HVAC]] (heating, ventilating and air-conditioning, which can be referred to as "building climate control"), to ensure more nearly instantaneous effects of a change in comfort level setting.<ref>Heat Transfer - A Practical Approach by ''Yunus A Cengel''</ref><!-- need a better reference-->
 
== Mechanical systems ==
The simplifying assumptions in this ___domain are:
* all objects are [[rigid body|rigid bodies]];
* all interactions between rigid bodies take place via [[kinematic pair]]s (''joints''), [[spring (device)|springsprings]]s and [[dashpot|dampers]].
 
== Acoustics ==
In this context, the lumped-component model extends the distributed concepts of [[Acousticacoustic theory]] subject to approximation. In the acoustical lumped-component model, certain physical components with acoustical properties may be approximated as behaving similarly to standard electronic components or simple combinations of components.
 
*A rigid-walled cavity containing air (or similar compressible fluid) may be approximated as a [[capacitor]] whose value is proportional to the volume of the cavity. The validity of this approximation relies on the shortest wavelength of interest being significantly (much) larger than the longest dimension of the cavity.
Line 181 ⟶ 160:
A simplifying assumption in this ___domain is that all heat transfer mechanisms are linear, implying that radiation and convection are linearised for each problem.
 
Several publications can be found that describe how to generate lumped-element models of buildings. In most cases, the building is considered a single thermal zone and in this case, turning multi-layered walls into lumped elements can be one of the most complicated tasks in the creation of the model. The dominant-layer method is one simple and reasonably accurate method.<ref>Ramallo-González, A.P., Eames, M.E. & Coley, D.A., 2013. Lumped Parameter Models for Building Thermal Modelling: An Analytic approach to simplifying complex multi-layered constructions. Energy and Buildings, 60, pp.174-184.</ref> In this method, one of the layers is selected as the dominant layer in the whole construction, this layer is chosen considering the most relevant frequencies of the problem. In his thesis,<ref>Ramallo-González, A.P. 2013. Modelling Simulation and Optimisation of Low-energy Buildings. PhD. University of Exeter.</ref>
 
Lumped-element models of buildings have also been used to evaluate the efficiency of domestic energy systems, by running many simulations under different future weather scenarios.<ref>Cooper, S.J.G., Hammond, G.P., McManus, M.C., Ramallo-Gonzlez, A. & Rogers, J.G., 2014. Effect of operating conditions on performance of domestic heating systems with heat pumps and fuel cell micro-cogeneration. Energy and Buildings, 70, pp.52-60.</ref>
 
== Fluid systems ==
Lumped-element models of buildings have also been used to evaluate the efficiency of domestic energy systems, by running many simulations under different future weather scenarios.<ref>Cooper, S.J.G., Hammond, G.P., McManus, M.C., Ramallo-Gonzlez, A. & Rogers, J.G., 2014. Effect of operating conditions on performance of domestic heating systems with heat pumps and fuel cell micro-cogeneration. Energy and Buildings, 70, pp.52-60.
Fluid systems can be described by means of [[Lumped parameter cardiovascular model|lumped-element cardiovascular models]] by using voltage to represent pressure and current to represent flow; identical equations from the electrical circuit representation are valid after substituting these two variables. Such applications can, for example, study the response of the human cardiovascular system to [[ventricular assist device]] implantation.<ref>Farahmand M, Kavarana MN, Trusty PM, Kung EO. "Target Flow-Pressure Operating Range for Designing a Failing Fontan Cavopulmonary Support Device" IEEE Transactions on Biomedical Engineering. DOI: 10.1109/TBME.2020.2974098 (2020)</ref>
</ref>
 
== See also ==
* [[Isomorphism#System isomorphism|System isomorphism]]
* [[Model order reduction]]
 
== References ==
{{Reflist}}
 
Line 196 ⟶ 178:
* [http://simulation.uni-freiburg.de/downloads/ims IMTEK Mathematica Supplement (IMS)], the Open Source IMTEK Mathematica Supplement (IMS) for lumped modelling
 
[[Category:Conceptual modellingmodels]]
[[Category:Mechanics]]
[[Category:Acoustics]]
[[Category:Electronic circuits]]
[[Category:Electronic design]]