Materials issues in thermal modeling of thin film electrocaloric solid-state refrigerators
expand article infoGunnar Suchaneck, Linda Felsberg, Gerald Gerlach
‡ Technische Universität Dresden, Dresden, Germany
Open Access


Materials properties affecting EC device operation are discussed based on an analytically tractable model of a layered EC refrigerator. Special attention was paid to thermal and interface thermal resistances. Estimates of the average cooling power of a stacked MEMS-based EC refrigerator were made.


electrocaloric cooling, heat transfer, thermal resistance, heat transfer fluids, cooling power

1. Introduction

The electrocaloric (EC) effect is a reversible temperature change of a material which results from an adiabatic application of an electric field. EC refrigeration is an environment-friendly thermal energy conversion technology. Electric fields required for the EC refrigeration cycle can be supplied much easier and less expensively than the high magnetic fields required for the magnetocaloric (MC) refrigeration. Moreover, electrical energy for EC cooling can be provided by stationary or mobile solar cells, and by electric vehicle batteries. This opens up completely new possibilities for an environment-friendly industrialization of developing countries [1].

The extraordinary dielectric response of relaxor ferroelectrics makes them promising for application in electrocaloric refrigerators [2]. Nowadays, the most studied and best performing EC materials are polyvinylidene fluoride (·CH2-CF2·)n terpolymers (P(VDF–TrFE–CFE)) and irradiated copolymers (P(VDF–TrFE)) as well as solid solutions of (1-x)Pb(Mg1/3Nb2/3)O3-xPbTiO3 (PMN–PT) with x < 0.3 [3].

A giant EC temperature change, ΔTEC = 45 K, was obtained for Pb0.88La0.08Zr0.65Ti0.35O3(PLZT) thin films on a Pt/TiO2/SiO2/Si substrate at an electric field of 125 V/μm [4]. Also, lead-free perovskite relaxors BaZrxTi1-xO3(BZT) provide a ΔTEC value sufficient for practical cooling applications over a broad temperature range [5].

An EC refrigerator, i.e. a heat pump, is able to transport thermal energy against a temperature gradient from the load to the EC layer and from here to the heat sink. It consists of a polarizable EC material contacted by electrodes and thermal boundaries to the load and the heat sink. The thermal connections have to be opened and closed appropriately as the layer is heated or cooled. Heat is transferred from the load or to the heat sink either (i) via controlled heat switches as well as uncontrolled thermal rectifiers, or (ii) by pumping a gaseous or liquid heat transfer agent through the solid refrigerant [1].

The EC response of a multilayer ceramic capacitor (MLCC) under real environment and operational conditions, but without thermal connections, was considered in [6]. The two-dimensional transient heat transfer model included the EC effect as a source term and accounted for radiative and convective heat losses at the interfaces. Earlier, cooling powers for MLCCs was modelled by introducing thermal capacitances and bulk thermal resistances of the capacitor sublayers [7]. Later, this approach was extended to a finite element analysis of six material combinations, electrode-to-EC-layer-thickness ratio and volume of the MLCC terminals [8]. A simplified one-dimensional lumped model in [9] assumed an on-state thermal contact resistance of the thermal switch of 10-6 m2K/W. A two-dimensional model of a parallel plate EC refrigerator with a heat transfer fluid circulating between the plates was solved numerically by the finite element method (FEM) in [10]. The model accounted for the coupled effects of temperature and fluid velocity fields. During the polarization/depolarization processes, thermal insulation was applied on all the walls. The presence of the cold/hot heat exchangers during the fluid flow process was modelled by fixing the temperatures on the boundaries. Different EC materials (PVDF-based relaxors, PMN-PT and PLZT) have been tested with this model.

In [11], FEM was used to evaluate a micro-electromechanical system (MEMS)-technology refrigeration system based on the EC effect in a P(VDF–TrFE–CFE) terpolymer. The cooling device included two diaphragm actuators fabricated in the plane of a silicon wafer, which drive a heat transfer fluid back and forth across of 5-layer stacks of terpolymer that are placed between them. Distributions of temperature, pressure, and flow velocity of the heat transfer fluid (Galden HT-70) in the x-direction in the element as well as time-averaged cooling power density and Carnot coefficients of performance were calculated. The hot and cold ends had a constant temperature. The remaining surfaces were considered as thermally insulated bearing significant overestimations. An achievable cooling power density of 3 W/cm2 was reported at a temperature span of 15 K and a coefficient of performance of 31% Carnot.

In [12], a thermoelectric/EC heat pump was modeled using FEM. The proposed device consists of thermoelectric layers acting as thermal switches laminated with EC layers. Ideal thermal contact between the EC and the thermoelectric layers was assumed, i.e. the contact thermal resistance was not taken into account. Any heat losses from the EC material as well as radiation and convection losses were neglected. Thus, the obtained cooling power density in the order of 1 W/cm2 is strongly overestimated.

Recently, an EC device was numerically modeled which comprises multi-layers of EC and heat storage material films separated by thermal switches changing the contact thermal conductance by forming a larger gap [13]. Calculation were carried out for best-case parameters: electric field E = 180 V/µm, EC temperature change ΔTEC = 24 K, negligible interface thermal resistances, and Fourier number Fo = 1. This results in an average cooling power of 2.4 W/cm2 and an average heat transfer efficiency of η = 10% for an operational frequency of f = 100 Hz and a temperature difference between the hot and the cold side of the system of ΔThc = 20 K. Note, that Fo represents the ratio of the heat conduction rate to the rate of thermal energy storage in a solid [14]. Complete heat absorption or rejection steps were realized in EC device prototypes for Fo ~ 8…10 [11, 15].

In this work, we derive an analytically tractable model of a layered EC refrigerator, we identify and discuss materials properties affecting EC device operation and we illustrate the impact of these properties on device operation by means of our simplified model.

2. Physical fundamentals

2.1. Heat diffusion equation

The variation of temperature within the EC material is described by the heat transfer equation. For the sake of simplicity, we consider in the following an one-dimensional model. In the case of linear heat flow, the parabolic heat diffusion equation of a solid comprising a heat source is given by [16]


where Θ = TT0 is the temperature difference to the environment, α = κ/c the thermal diffusivity, κ the thermal conductivity, and c the volumetric specific heat assumed to be constant (c (E,T) ≈ c), V the volume, A the area and d the thickness of the EC layer, and h the heat transfer coefficient at the boundary describing heat loss to the environment. Introducing a bulk thermal resistance Rth,b = (κA)-1 and a thermal capacity Cth,b = cA, both per unit length, eq. (1) transforms into


where Q˙′ and h′ are the values of Q˙ and h per unit length, respectively. Since the Biot number Bi = hd/κ (characterizing the ratio of the thermal resistance of the EC materials bulk to that at the boundary) will be less than 0.1 up to values of d in the order of 100 μm, the temperature of the EC element during heat transfer remains nearly constant, i.e. the ∂2Θ/∂x2 = 0. This enables a lumped system approximation [14] simplifying eq. (2) to


where Rth,i = ARth,i = 1/h is the area-specific thermal resistance at the interface, Cth = cd the thermal capacitance of the EC layer per unit area, and Φ the heat flux through the interface. Eq. (3) represent the well-known bolometer equation [17].

The driving force of heat flux is the temperature difference Θ caused by charging or discharging the EC element which itself represents a dielectric capacitor. When charging a capacitor possessing a dielectric permittivity ε under isothermal conditions, the heat dQ delivered by a small change of the electric field dE is given by [18]


Since the EC temperature change is usually small, ΔΘEC << Θ, the corresponding thermal flux Φ(t) = dQ (t)/(Adt) might be written as


where <Τ> is the average temperature of the cooling cycle and where the field dependence of the temperature coefficient of dielectric permittivity 1/ε·∂ε/∂T is neglected. Note, that above the temperature Tm of maximum dielectric permittivity, 1/ε·∂ε/∂T is negative. The rise time of an output signal of a voltage pulse generator is mainly dependent on the switching speed of the power amplifier, while the fall time is determined by the resistor-capacitor (RC) circuitry at the generator output and the load. Here, we will not consider the physical origin of power amplifier switching. Assuming for the rising pulse


with τe = RC the time constant of the electronic circuit, an approximate solution of eq. (3) is given by


where τth = Rth,iCth and where the temperature coefficient of the dielectric permittivity was assumed to be constant. For very short rise times of the electric field, τe→0, eq. (7) reduces to


A pulse decay given by


yields for t´ = tt0


transforming in terms of t′ which for τe → 0 is equivalent to eq. (8) with opposite sign. A flexible EC cooling device with electrostatic actuation reported in [15] was driven with a high-voltage circuitry possessing an electrical time constant of about 140 ms. On the other hand, an analysis of eqs. (7) and (10) suggests time constants of less than 50 ms for an efficient device operation.

For sake of simplicity, we describe for a moment the thermodynamic cycles of the EC element by a periodic temperature change with a fundamental angular frequency ω = 2π/τc, where τc is the cycle time. Since the parabolic heat diffusion equation (1) has the same structure as the equation of an electrical RC-transmission line, multilayer problems are best treated by the matrix methods commonly used in electrical engineering [19]. In our case, temperature plays the role of a voltage and thermal flux the role of an electrical current. Each point, we shall describe by two quantities, the temperature Θ and the heat flux Φ. Denoting Θ and Φ as the temperature and thermal flux of a slab at the face z = 0 and Θ′ and Φ′ as their values at the face z = d, the following relations yield [16]

Θ′ = LΘ + MΦ,

Φ′ = NΘ + OΦ, (11)


L = cosh(kd), M = –(κk)-1sinh(kd),

N = –κksinh(kd), O = cosh(kd), (12)



where dD = (2D/ω)1/2 = (Dτc/π) is the penetration depth of the thermal oscillation. Note, that for dielectric thin and thick films with d < 100 µm we get L = O ≈ 1, M = –d/κ, and N = –κk2d for ω = 1 Hz since D is in the order of 10-6 m2/s. The value of M represents the area-specific thermal resistance of the film bulk.

Eq. (11) may be considered as matrix equation


In eq. (14), Θ and Φ must be multiplied by their time dependence exp(iωt) which is omitted here entirely. It should be again included when real or imaginary parts have to be taken at the end of the calculation. Assuming perfect thermal contact between the faces of the slabs a comprising n sublayers, the rth being of thickness dr, thermal conductivity κr and thermal diffusivity Dr with Θr and Φr and Θ′r and Φ′r at its left-hand and right-hand faces, respectively, recursive application of eq. (14) gives for a one-dimensional heat flow across several sublayers


Given any two of Θr, Φr, Θ′r and Φ′r the other two can be found using this method. The multiplication of matrices in eq. (15) can be carried out successively. However, deriving explicit formulas for slabs of n layers is very cumbersome [20–22].

In the presence of thermal contact resistances between the slabs or at the surfaces, they may be expressed also in matrix notation and included into eq. (15) [16]. For example, in case of a thermal interface resistance Rth,i1 between the first and second slab we find


and, correspondingly,


The final result of this calculation is are linear relations between the temperature and fluxes Θ1, Θ′n, Φ1, Φ′n at the two surfaces of the composite slab. The surface conditions will provide two more equations so that all four quantities Θ1, Θ′n, Φ1, Φ′n will be determined. Recursive application of this matrix scheme yields for a one-dimensional heat flow across j slabs


where M is the heat transfer matrix. Thus, the matrix method allows numerically a very easy analysis of multilayer structure by multiplying numerical matrices.

The relative thickness of the EC layer in prototype EC cooling devices amounts to 0.30 to 0.35 dD [11, 15]. In accordance with the above made approximation of a thermally thin system, we consider in the following a stack of j films satisfying dj << dD. In this case, we obtain:


2.2. Electrocaloric material

According to eq. (8), the EC temperature change given by ΔTEC = Θ(t = 0) is defined by the volumetric specific heat c and the temperature coefficient dε/dT of dielectric permittivity. Since the values of c at room temperature approach the high temperature limit, it is expected that they do not vary strongly among various pyroelectric materials [23]. Therefore, dε/dT is a key parameter. Relaxor ferroelectrics possess above Tm large values of 1/ε·∂ε/∂T as high as a few 102 K-1 (Table 1).

Temperature coefficient of dielectric permittivity above Tm for BZT20 and selected relaxor ferroelectrics.

Refrigerant ε(Tm) dε/dT, K-1 Ref.
BaZr0.2Ti0.8O3 10600 -100 [5]
-87 [24]
BaZr0.25Ti0.75O3 6950 -72.5 [25]
18940 -377 [26]
BaZr0.3Ti0.7O3 33400 -560 [27]
BaZr0.35Ti0.65TiO3 11550 -81.5 [28]
BaSn0.24Ti0.76TiO3 16000 -30 [29]
Ba0.20Pb0.80ZrO3 11080 -340 [30]
0.9PbMg1/3Nb2/3O3-0.1PbTiO3 11400 -130 [31]

Recently, BaZr0.2Ti0.8O3 (BZT20) was demonstrated to be a promising material for EC application [5, 24]. BZT20 is a composition near the invariant critical point, where four different phases, including cubic, tetragonal, orthorhombic, and rhombohedral phases coexist. It is located slightly below the crossover to relaxor behavior [28, 32] where an intermediate state between ferroelectric and relaxor behavior is expected. On the other hand, BZT20 thick films show the shift of the dielectric permittivity peak towards higher temperatures with increasing frequency characteristic for relaxor ferroelectrics [24].

EC refrigerants are evaluated by a material criterion which characterizes the efficiency of the physical cooling process and, therefore, is independent on the performance of different thermodynamic cycles [1]:


Here, εε0·E2·tanδ is the non-recoverable electrical loss with tanδ the loss tangent, and cΔTEC the heat transferred from the load to the heat sink within one refrigeration cycle with c the specific heat of the cooling element and ΔTEC the EC temperature change. At this point, the tanδ comes into play as an additional material parameter since – following eq. (8) – the dielectric permittivity cancels out. It has to been noted that tanδ of eq. (20) needs to be investigated in the frequency range of EC device operation (0.1 to 10 Hz). In the presence of a ferroelectric hysteresis, it mainly consists of losses due to minor hysteresis loops described by an apparent loss tangent [33]:


where P–(E) and P+(E) are the field dependence of dielectric polarization at the upper and lower branch of the minor hysteresis loop, respectively. Obviously, tanδ depends significantly on materials synthesis. Values of tanδ < 0.07 provide refrigerant efficiencies of Φ > 0.95 exceeding the ones of any other solid state cooling technology.

2.3. Thermal interface resistance

At thin film interfaces, the highest value of R"th,i measured at room temperature to date is 1.2·10-5 m2K/W for Bi/Hydrogen-terminated diamond. The thermal interface resistance at room temperature of other metal/dielectric combinations falls within a relatively narrow range, 3.3·10-8 m2K/W < R"th,i < 12·10-8 m2K/W [34].

The minimum possible R"th,iph-ph produced by a purely harmonic process involving two phonons (one phonon on each side of the interface) is given by the phonon radiation limit which in the limit of high enough temperatures, where the Dulong-Petit law holds, is given by


Here, k is the Boltzmann constant, fmax = kΘ/h the vibrational cutoff frequency of the metal, h the Planck constant, Θ the Debye temperature, and vD,


the Debye velocity of a two dimensional dielectric with vl and vt the average longitudinal and transverse sound velocities, respectively. A compilation of all data used for calculations is given in Table 2. In result, eq. (22) yields a phonon-phonon contribution to the interface resistance of R"th,iph-ph = 4.1·10-10 m2K/W.

Material parameters used for the calculation of the interface thermal resistance.

Material Parameter Value Ref.
Ni Θ 495 K [35]
vl 6040 m/s [36]
ρ 8.902 g/cm3 [36]
M 58.69 g/mol [36]
l ph 1.2 nm [37]
c e 0.3174 MJ/m3K [38]
Cu τe-ph 1…4 ps [39]
W ~0.5 ps [40]
BaTiO3 vl 6860 m/s [41]
vt 3870 m/s
κ 2.61 W/mK [42]
BaSr0.3Ti0.7O3 εb 420 [43]
BaSr0.48Ti0.52O3 380 [44]
BaSr0.55Ti0.45O3 566 [45]
BaSr0.65Ti0.35O3 488 [46]
Pt/BaSr0.3Ti0.7O3/Pt Ci 0.032 F/m2 [43]
Pt/BaSr0.48Ti0.52O3/Pt difif 0.076 [44]
Pt/BaSr0.55Ti0.45O3/Pt 0.048 [45]
Pt/BaSr0.65Ti0.35O3/Pt 0.174 [46]

While phonons dominate the heat conduction in dielectrics, electrons dominate it in metals. Hence, an energy transfer occurs between electrons and phonons during heat transport across the metal-dielectric interface. In result, two thermal resistances - volumetric electron-phonon (R"the-ph) in the metal and interfacial phonon (R"th,iph-ph) – are in series. The former one is given by [47]


where Re is the electron cooling rate or electron-to-phonon energy transfer per unit volume, Re = ce/τ with ce the electron heat capacity per unit volume and τe-ph is the relaxation time characterizing electron-phonon energy loss or electron cooling. The value κph is the phonon thermal conductivity,


with ρ the density, NA the Avogadro constant, M the molar mass, and lph the phonon mean-free-path of the metal. Assuming τe-ph = 2 ps (cf. values of W and Cu in Table 2), we arrive at R"the-ph = 8.6·10-10 m2K/W. Thus, the interfacial thermal resistance in the order of 10-9 m2K/W calculated by simplified models is much lower than the corresponding experimental data.

From a microscopic point of view, all surfaces in contact have deviations from their idealized geometry. Because of these imperfections, two bodies in contact touch actually only at a few discrete points. The heat flux close to the interface is then constricted in these micro-contact regions [48].

The thermal contact resistance is defined as the ratio between the temperature drop at the interface and the heat flux crossing the interface. Similarly, electrical capacitance C is the ratio of the total charge on a conductor to its electric potential (potential difference to a reference electrode). Both capacitance and thermal conductance are affected by edge effects. Therefore, there is a close relation between Rth,i and the interface capacitance Ci given by [49]


where εb is the bulk dielectric permittivity. The area-related value Rth,i = ARth,i is then given by


with di the thickness and εi the dielectric permittivity of the interface layer. With regard to the values of Ci and εb in Pt/BaSrxTi1-xO3/Pt structures (cf. Table 2), Rth,i amounts to 1…5·10-8 m2K/W crossing the lower limit of the thermal interface resistance for thin film metal/ dielectric contacts reported in [34].

In macroscopic cooling systems, the thermal interface resistance is generally negligible since it is in series with other thermal resistances exhibiting much larger dimensions and, thus, much larger resistance. Here, heat is transferred from the load or to the heat sink (i) via thermal switches or (ii) by a gaseous or liquid heat transfer agents [1]. A compilation of typical interface resistances of solid-solid and liquid-solid thermal contact resistances is given in Table 3.

Thermal contact resistances of solid-solid and liquid-solid configurations.

Contact Rth,c, m2K/W Parameter Ref.
Si/Si ~2·10-4 0.169 MPa [50]
Cu/Si ~1·10-4
Cu-In/Si 2.72·10-5
Cu-PCM1)/Si 5.22·10-5
Cu-CNT/Si ~2.5·10-5
Cu-PCM1)-CNT/Si ~3.4·10-5
poly-Si/SiNx 1.7·10-5 [51]
Cu/Cu (sheets) 1.59·10-5 3.6 MPa [52]
1.05·10-5 14.4 MPa
Al/Al (sheets) 4.76·10-5 3.6 MPa
4.00·10-5 14.4 MPa
1.09·10-6 K = 110 (0.1N) [53]
Si/Hg 2.6·10-6 224 (1N)
CNT columns/Si 5.7·10-5 - (0.1 N)
CNT columns/Si 2.04·10-5 8.4 (1N)
CNT columns-Au/Si 4.16·10-5 4.4
CNT/Si 4.6 (0.1N)
CNT/Si 8.87·10-5 27.8 (1N)
CNT/CNT 4.23·10-4 6.1
CNT-Polyimide ~2·10-3 2) [15]
Hybrid solid-liquid (Cu-H2O)-PTFE/Si 1.3·10-5 d H2O = 75 µm [54]
Si-H2O-PTFE/Glass 1.0·10-5 d H2O = 6.1 µm [55]
2.0·10-5 11.6 µm
3.2·10-5 19.3 µm
Cu/PDMS 3.25·10-5 0.46 MPa [56]
ZnO, BN/PDMS ~9·10-6
Carbon fiber/epoxy resin 0.3-1.8·10-5 [57]
Pt/SiO2 3.2·10-9 d SiO2 = 26, 440 nm [58]
Si/SiO2 < 1.7·10-9 d SiO2 = 26 nm
Al/SiO2 0.7-1·10-8 d Al = 50-120 nm [59]
Si/SiO2 6-8·10-8 d SiO2 = 110-518 nm
Al/epoxy resin 4.6·10-4 [60]
Glass/epoxy resin 8.6·10-4
Grease 0.6-2·10-5 [61]
Galistan 7.7·10-6 0.172 MPa [62]
Galistan 2.8·10-5 0.284 MPa [63]

2.4. Heat transfer fluids

In an active EC regenerator configuration, liquid or gaseous heat transfer agents are pumped through the EC refrigerant. They absorb heat from the bed and transfer it to the cold sink. Thereby, an electric field is applied solely to the refrigerant. The flow channels are usually rectangular ducts requiring a two-dimensional flow analysis. On the other hand, in a one-dimensional model, the determination of a heat transfer coefficient h (the inverse of h represents the thermal interface resistance) describing the heat transfer between the fluid and the solid is sufficient. The heat transfer across a boundary is defined by the Nusselt number Nu – the ratio of convective to conductive heat transfer normal to the considered boundary:


with dH the hydraulic diameter. For steady state, laminar flow of an incompressible liquid with constant physical properties in a rectangular duct of constant cross-sectional area characterized by its aspect ratio AR (the ratio of the-smaller-to-the larger duct size), the value of Nu is approximately given by [64]:

Nu = 8.235(1 – 2.0421AR + 3.0853AR2

2.4765AR3 + 1.0578AR4 – 0.1861AR5). (29)

At Bi < 0.1, the temperature of the EC element during heat transfer remains nearly constant (cf. section 2.1). Here, a duct high of z = 0.6 mm and AR = 0.15 yield thermal resistances of


The thermal resistance of the fluid, Rth,f, is determined by the inverse of the product of its mass flow rate dm/dt and its specific heat at constant pressure cp [16]. The total thermal resistance of the fluid, Rth,f, for a given fluid velocity v then yields [65]


For calculation, an upper limit of the fluid velocity of v = 0.1 m/s was chosen providing a back and forth fluid motion of 1 cm at a cycle time of τc = 0.2 s. This results in vc >> h leading to Rth,fRth,s-f. Values of Rth,s-f are compiled in Table 4. Also, the thermal penetration depth at ω = 1 Hz, dD,1Hz, is listed. It allows a quick estimation of dD at a given cycle time by

Physical properties of heat transfer fluids, thermal resistance at the solid-fluid interface Rth,s-f for laminar flow, and values of dD,1Hz.

Heat transfer fluid κ, W/mK c, MJ/m3K Rth,s-f , 10-3m2K/W dD ,1Hz, mm
Silicone oil, 20 cSt 0.142a 1.52a 1.36 0.30
Water 0.606b 4.19b 0.32 0.38
HT 70 0.07c 1.62 2.76 0.21


Summing up, the thermal resistance at the interface to the heat transfer fluid limits the heat transfer of active EC regenerators.

3. Electrocaloric device and cooling power

A MEMS-based EC element with a stack of three cooling cells is illustrated in Figure 1. Multilayer ceramic thin films with interdigitated electrodes replicating multilayer ceramic capacitors (MLCCs) [7] are positioned on a thin SiNx membrane inside a Si frame produced by silicon wet etching. For thermal insulation, the inclined walls of the silicon die are covered by thin SiO2 film as adhesion layer and a polydimethylsiloxane (PDMS) film as thermal barrier. The silicon dies are stacked using PDMS interlayers consisting of fluid channels. The bulk thermal resistances of the constituent layers are compiled in Table 5. A contact thermal resistance of the BZT/Pt interface amounting Rth,c = 5·10-8 m2K/W was calculated by means of eq. (27). For the Pt/SiNx interface, a value of Rth,c = ~3·10-9 m2K/W obtained in literature for Pt/SiO2 [58] was used.

Figure 1.

(a) Electrocaloric element: 1, 2 – outer electrodes, 3 – EC MLC, 4 – SiNx membrane, 5 – Si wafer, 6 – cavity sidewall coating and (b) stack of EC three cells.

Bulk thermal resistances of the constituent layers of the EC device in Figure 1.

Layer κ, W/mK d, µm Rth, 10-3m2K/W Ref.
Thermal insulation 10 [69]
Ni electrodes 90.7 200×2 4.41·10-2 [7]
BZT20 layer 3.4 200×6.5 0.5 [5]
Pt electrode 71.6 0.1 1.4·10-6 [36]
SiNx membran 10 0.2 2.0·10-5 [70]
PDMS coating 0.2 1000 5 [71]

The thermodynamic cycle of an EC refrigeration process in such a cooler design includes four steps [11, 72]: (i) Adiabatic polarization: the EC material is heated by applying the electric field; (ii) Heat rejection: while the electric field is on, the fluid absorbs the heat from the EC material as it flows from the cold side to the hot side rejecting the heat carried by the fluid to the hot side (heat sink); (iii) Adiabatic depolarization: the EC materials is cooled by turning off the electric field; (iv) Heat absorption: the fluid moves back from the hot side to the cold side, releases heat to the EC material and absorbs heat from the cold side (load). Thus, by periodic cycling, heat is extracted from load and released to the heat sink. Here, the fluid reaches its lowest temperature in the cold side.

For a cycle time τc = 2mCthRth proportional to the thermal time constant of the EC device, the cycle-averaged cooling power of a single cell thermally isolated at the top is given by [9, 72]:


whith ΔΘ* ≈ ΔTEC and


Here, TI denotes the thermal insulation layer. Neglecting the small bulk and interface thermal resistances, assuming Rth,bTI → ∞ – equivalent to zero temperature gradient at the boundary between the cells - and taking into account, that the heat flux is divided between the heat transfer fluid and its PDMS environment, we arrive at


Taking water as heat transfer fluid, assuming ΔTEC = 5 K and m = 2, the cycle-averaged cooling power of a single cell amount to 180 mW/cm2. Following Table 3, the corresponding <q˙″> values for silicone oil and HT-70 amount to 51.5 and 25.7 mW/cm2, respectively. This is comparable to an experimental value of 30 mW/cm2 reported recently for a cooling device using a flexible electrocaloric (EC) polymer film and an electrostatic actuation mechanism [15]. The dependence of the averaged cooling power of a single cell on cycle time is depicted in Figure 2.

Figure 2.

Cycle-averaged cooling power of a single EC cell in dependence on cycle time.

For a cell stack, only the bottom and top cells should be thermally insulated to the environment by a layer of low thermal conductivity. Here, MEMS-structures consisting of vacuum cavities [69] or highly porous materials find application. At the top and bottom surfaces, an additional heat flux divider appears which lower the cooling power of the cell by a factor of


Figure 3 illustrated the dependence of F on the number of stacked cells assuming that a vertical temperature gradient is present only in the top and bottom cells of the stack.

Figure 3.

Efficiency factor versus number of stacked cells.

The estimations made in this work illustrate that a solid state EC refrigerator with a cooling power of up to 2 W/cm2 can be realized using a stack of 10 single MEMS-based cells.


A simplified analytical model of solid state EC refrigerators allows analysis of the materials impact on device operation, enables the indentification of the key material parameters and makes it possible to predict the device performance. For a stack of MEMS-based EC cells in an active regenerator mode, cycle-averaged cooling powers of up to 2 W/cm2 were estimated.


This work was supported by the German Research Foundation (DFG) within the Priority Program «Ferroic Cooling» (SPP1599, project B6). The authors thank R. Liebschner for providing Figure 1b.


  • 4 Lu S.G., Rožič B., Zhang Q.M., Kutnjak Z., Li Xinyu, Furman E., Gorny L. J., Lin M., Malič B., Kosec M., Blinc R., Pirc R. Organic and inorganic relaxor ferroelectrics with giant electrocaloric effect. Appl. Phys. Lett. 2010; 97(16): 162904.
  • 5 Qian Xiao-Shi, Ye Hui-Jian, Zhang Ying-Tang, Gu Haiming, Li Xinyu, Randall C.A., Zhang Q.M. Giant electrocaloric response over a broad temperature range in modified BaTiO3 ceramics. Adv. Funct. Mater. 2014; 24(9): 1300.
  • 6 Smith N.A.S., Rokosz M.K., Correia T.M. Experimentally validated finite element model of electrocaloric multilayer ceramic structures. J. Appl. Phys. 2014; 116(4): 044511.
  • 7 Kar-Narayan S., Mathur N.D. Predicted cooling powers for multilayer capacitors based on various electrocaloric and electrode materials. Appl. Phys. Lett. 2009; 95(24): 242903.
  • 8 Crossley S., McGinnigle J.R., Kar-Narayan S., Mathur N.D. Finite-element optimisation of electrocaloric multilayer capacitors. Appl. Phys. Lett. 2014; 104(8): 082909.
  • 9 Ju Y.S., Solid-state refrigeration based on the electrocaloric effect for electronics cooling. J. Electron. Packag. 2010; 132(4): 041004.
  • 10 Aprea C., Greco A., Maiorino A., Masselli C. A comparison between different materials in an active electrocaloric regenerative cycle with a 2D numerical model. Int. J Refrig. 2016; 69: 369.
  • 11 Guo D., Gao J., Yu Y.-J., Santhanam S., Slippey A., Fedder G.K., McGaughey A.J.H., Yao S.C. Design and modeling of a fluid-based micro-scale electrocaloric refrigeration system. Int. J. Heat Mass Transf. 2014; 72: 559.
  • 12 Feng D., Yao S.-C., Zhang T., Zhang Q. Modeling of a smart heat pump made of laminated thermoelectric and electrocaloric materials. J. Electron. Packag. 2016; 138(4): 041004.
  • 13 Hirasawa S., Kawanami T. Shirai K. Electrocaloric refrigeration using multi-layers of electrocaloric material films and thermal switches. Heat Transfer Eng. 2018; 39(12): 1091–1099.
  • 14 Bergmann T.L., Lavine A.S., Incropera F.P., Dewitt D.P. Fundamentals of Heat and Mass Transfer, 7th edition. Hoboken (NJ): John Wiley & Sons, 2011, 1050 p.
  • 15 Ma R., Zhang Z., Tong K., Huber D., Kornbluh R., Ju Y. S., Pei Q. Highly efficient electrocaloric cooling with electrostatic actuation. Science. 2017; 357(6356): 1130.
  • 16 Carslaw H.S., Jaeger J.C. Conduction of Heat in Solids, 2nd Edition. Oxford: Oxford Science Publications, 1959, 510 p.
  • 17 Putley E.H. The pyroelectric detector. in Semiconducton and Semimetals. Vol. 5. R.K. Willardson and A.C. Beer (Eds.). New York: Academic Press, 1970, pp. 259–285.
  • 18 Kubo R. Thermodynamics. An Advanced Course with Problems and Solutions. Amsterdam; New York: North-Holland Pub. Co, 1968, Ch. 31, chapter 3 § 31, p. 186.
  • 19 Adby P.R. Applied Circuit Theory: Matrix and Computer Methods (Ellis Horwood Series in Electrical and Electronic Engineering). London: Ellis Horwood Ltd., 1980.
  • 22 Samoilov V.B., Yoon Y.S. Frequency response of multilayer pyroelectric sensors. IEEE Trans. Ultrason. Ferroelectr. Freq. Control. 1998; 45(5); 1246.
  • 24 Ye H.-J., Qian X.-S., Jeong D.-Y., Zhang S., Zhou Y., Shao W.-Z., Zhen L., Zhang Q.M. Giant electrocaloric effect in BaZr0.2Ti0.8O3 thick film. Appl. Phys. Lett. 2014; 105(15): 152908.
  • 25 Kell R.C., Hellicar N.J. Structural transitions in barium titanate-zirconate transducer materials. Acta Acustica united with Acustica. 1956; 6(2): 235–245.
  • 27 Yu Z., Ang C., Guo R., Bhalla A.S. Ferroelectric-relaxor behavior of Ba(Ti0.7Zr0.3)O3 ceramics. J. Appl. Phys. 2002; 92(5): 2655.
  • 28 KleemannW., Miga S., Dec J., Zhai J. Crossover from ferroelectric to relaxor and cluster glass in BaTi1-xZrxO3 (x = 0.25–0.35) studied by non-linear permittivity. Appl. Phys. Lett. 2013; 102(23): 232907.
  • 29 Lei C., Bokov A.A., Ye Z.-G. Ferroelectric to relaxor crossover and dielectric phase diagram in the BaTiO3–BaSnO3 system. J. Appl. Phys. 2007; 101(8): 084105.
  • 31 Hirose S., Usui T., Crossley S., Nair B., Ando A., Moya X., Mathur N.D. Progress on electrocaloric multilayer ceramic capacitor development. APL Mater. 2016; 4(6), 064105.
  • 36 Lide D.R. (Ed.). CRC Handbook of Chemistry and Physics. Boca Raton (FL, USA): CRC Press. 2005. 2660 p.
  • 38 Stoner E.C. VI. The specific heat of nickel. The London, Edinburgh, and Dublin Philosophical Magazine and Journal of Science. Series 7. 1936; 22(145): 81–106.
  • 39 Fujimoto J.G., Liu J.M., Ippen E.P. Femtosecond laser interaction with metallic tungsten and nonequilibrium electron and lattice temperatures. Phys. Rev. Lett. 1984; 53(19): 1837.
  • 44 Park W.Y., Hwang C.S. Film-thickness-dependent Curie-Weiss behavior of (Ba,Sr)TiO3 thin-film capacitors having Pt electrodes. Appl. Phys. Lett. 2004; 85(22): 5313.
  • 45 Lee B.T., Hwang C.S. Influences of interfacial intrinsic low-dielectric layers on the dielectric properties of sputtered (Ba,Sr)TiO3 thin films. Appl. Phys. Lett. 2000; 77(1): 124.
  • 46 Wang J., Zhang T., Zhang B., Jiang J., Pan R., Ma Z. Interfacial characteristic of (Ba,Sr)TiO3 thin films deposited on different bottom electrodes. J. Mater. Sci: Mater. Electron. 2009; 20(12): 1208.
  • 47 Majumar A., Reddy P. Role of electron-phonon coupling in thermal conductance of metal-nonmetal interfaces. Appl. Phys. Lett. 2004; 84(23): 4768.
  • 48 Mantelli M.B.H., Yovanovich M.M. Thermal Contact Resistance, ch. 16. In: Spacecraft Thermal Control Handbook. Vol. I: Fundamental Technologies. El Segundo (CA): The Aerospace Press, 2002, 836 p.
  • 49 Iossel’ Yu.Ya., Kochanov E.S., Strunsky M.G. Rasshet elektrisheskoi emkosti [Calculation of Electrical Capacity]. Leningrad: Energoisdat,. 1981, 288 p. (In Russ.)
  • 51 Song W.-B., Sutton M.S., Talghader J.J. Thermal contact conductance of actuated interfaces. Appl. Phys. Lett. 2002; 81(7): 1216.
  • 54 Jia Y., Ju Y.S. Solid-liquid hybrid thermal interfaces for low-contact pressure thermal switching. J. Heat Transfer. 2014; 136(7): 074503.
  • 55 Cha G., Ju Y.S. Reversible thermal interfaces based on microscale dielectric liquid layers. Appl. Phys. Lett. 2009; 94(21): 211904.
  • 56 Xu Y., Luo X., Chung D.D.L. Sodium silicate based thermal interface material for high thermal contact conductance. J. Electron. Packag. 1999; 122(2): 128.
  • 57 Kumar K., Ayyagari N., Fisher T.S. Effects of graphene nanopetal outgrowths on internal thermal interface resistance in composites. ACS Appl. Mater. Interfaces. 2016; 8(10): 6678.
  • 58 Kimling J., Philippi-Kobs A., Jacobsohn J., Oepen H.P., Cahill D.G. Thermal conductance of interfaces with amorphous SiO2 measured by time-resolved magneto-optic Kerr-effect thermometry. Phys. Rev. B. 2017; 95(18): 184305.
  • 59 Zhu J., Tang J.D., Wang W., Liu J., Holub K.W., Yang R. Ultrafast thermoreflectance techniques for measuring thermal conductivity and interface thermal conductance of thin films. J. Appl. Phys. 2010; 108(9): 094315.
  • 60 Pietrak K., Wiśniewski T.S., Kubiś M. Application of flash method in the measurements of interfacial thermal resistance in layered and particulate composite materials. Thermochimica Acta. 2017; 654: 54.
  • 61 Narumanchi S., Mihalic M., Kelly K., Eesley G. Thermal interface materials for power electronics applications. 11th Intersociety Conf. Thermal and Thermomechanical Phenomena in Electronic Systems (ITHERM 2008). Orlando, FL (USA), 2008, 10 pp.
  • 62 Khounsary A.M., Chojnowski D., Assoufid L., Worek W.M. Thermal contact resistance across a copper-silicon interface. Proc. SPIE. Vol. 3151, High heat flux and synchrotron radiation beamlines. San Diego (CA, USA), 1997, pp. 45–51.
  • 63 Asano M., Ogata J., Yosinaga Y. Quantitative evaluation of contact thermal conductance in a vacuum as a result of simulating the effect of cooling. SPIE Proc. Vol. 1739, High heat flux engineering. San Diego (CA, USA), 1993, pp. 652–656.
  • 64 Shah R.K., London A.L. Laminar Flow Forced Convection in Ducts. New York: Academic Press, 1978, p. 205.
  • 65 Takács G., Szabó P.G., Bognár Gy. Modelling of the flow-rate dependent partial thermal resistance of integrated microscale cooling structures. Microsyst. Technol. 2017; 23(9): 4001.
  • 69 Koehler R., Bruchhaus R., Pitzer D., Primig R., Schreiter M., Wersing W., Winkler B., Gerlach G., Hofmann G., Heß N. Pyroelectric thin film presence detector arrays with micromachined pixels. Integr. Ferroelectrics. 2002; 44(1): 77.
  • 70 Zhang X., Grigoropoulos C.P. Thermal conductivity and diffusivity of freestanding silicon nitride thin films. Rev. Sci. Instrum. 1995; 66: 1115.