Method for operating a variable displacement compressor and system with a variable displacement compressor

Abstract

A system and a method for operating a variable-displacement compressor are provided wherein a volume of a compression chamber of at least one cylinder is varied during a compression cycle, wherein during the cycle the at least one cylinder is alternately coupled to a low pressure side and a high pressure side of the compressor. A first thermodynamic state of the low pressure side and a second thermodynamic state of the high pressure side can be estimated and a current compression ratio

Description

BACKGROUND AND SUMMARY

The invention relates to a method for operating a variable displacement compressor and a system with a variable displacement compressor.

Variable displacement compressors are known in the art. Variable displacement compressors of the swash-plate type are widely used in automotive air conditioners. Such compressors include a displacement control mechanism, which operates to maintain the exit pressure of an evaporator or the suction pressure of a compressor associated with the exit pressure at a predetermined target value. The exit pressure of an evaporator or the suction pressure of a compressor associated with the exit pressure reflects the magnitude of the cooling load. The variable control mechanism controls the exit pressure of the evaporator or the suction pressure to adjust the inclination angle of the swash plate of the compressor which adjusts the displacement of the compressor. A compressor used in vehicles is usually driven by the vehicle engine. Being a main consumer of engine power, the load of the compressor to the engine has to be reduced under exceptional conditions, e.g. during acceleration or moving uphill or when starting the engine.

U.S. Pat. No. 4,578,960 B2 teaches a displacement control method for a variable displacement compressor which is incorporated in a refrigerant circuit of a vehicle air conditioning system. The displacement of the compressor is switched dependent on the temperature of the refrigerant. A control valve is provided for changes of the displacement of the compressor by changing the valve opening degree. The displacement is switched to a value between 30-50% of the maximum displacement if the temperature of the refrigerant drops below −1° C. If the temperature of the refrigerant is above 3° C., the compressor is operated at its maximum displacement.

The U.S. Pat. No. 6,484,520 B2 discloses a displacement control apparatus, wherein a variable displacement compressor of a swash plate type is employed. A drive shaft of the compressor is driven by a vehicle engine. A control valve changes the displacement of the compressor. An air conditioner controller produces a torque setting signal to the compressor controller which changes the valve opening and subsequently the displacement of the compressor. Thus, the compressor is controlled according to the torque.

It is desirable to provide a reliable and fast method for operating a variable displacement compressor. It is also desirable to provide a reliable system for operating a compressor. It is also desirable to provide a suitable computer program as well as a computer program product.

According to a first aspect of the invention a method for operating a variable-displacement compressor is provided, wherein a volume of a compression chamber of a cylinder is varied during a compressor cycle, also called compression cycle, wherein during the cycle said cylinder is alternately coupled to a low pressure side and a high pressure side of said compressor. The method comprises the steps of estimating a first thermodynamical state of said low pressure side and a second thermodynamical state of said high pressure side, estimating a current compression ratio of said compressor, calculating said displacement as a function of said compression ratio, estimating a compressor torque depending on the calculated displacement, and providing said compressor torque to a control unit for operating the compressor or for operating an engine coupled to the compressor. The control unit can operate the compressor directly or indirectly. For instance, the control unit can be an engine control unit where the engine drives the compressor. Knowing the displacement is useful when estimating the compressor torque, particularly when the compressor is employed in an air conditioning system of a vehicle. If the compressor torque is known, the control unit, such as an electronic control unit (ECU) or an engine control unit of a vehicle or the like, can easily compensate for the torque, thus making the control of the respective component, for instance the engine, easier. The drivability of the vehicle can also be improved. Particularly, the drivability is improved in the way that jerk problems are decreased. Other benefits gained by better compressor torque estimation which can be sent to an engine management system of a vehicle are decreased emissions of the engine as well as a decreased fuel consumption of the vehicle.

Additionally, the preferred method can easily be transferred to other systems. When employed into an air conditioning system, the method allows for better controlling of an evaporator coupled to the air conditioning system, particularly for better controlling of the evaporator temperature. Control of the evaporator can mean e.g. control of a fluid temperature, e.g. air temperature directly downstream the evaporator between the evaporator and the compressor and/or control of a coolant pressure between the evaporator and the compressor, which pressure can preferably correlate with the air temperature. Preferably, the evaporator can also be controlled by a control unit. For instance, the compressor control and torque estimation can be placed in the same control unit as the engine management software, or, for instance, a climate control unit can control the evaporator. A fluid which is compressed in the compressor can be a gas or a liquid.

Other than non-physical models, which depend strongly on measured data, the proposed method is based on a physical model of the compressor. As a physical model, the method is rather independent of the ambient environment of the compressor. A transfer of the method to another system can be done easily and fast, which is time and cost efficient. Knowing the actual displacement of the compressor is not only favorable for estimating the compressor torque but also for easier controlling an engine of an engine driven compressor, for instance when the compressor is coupled to a vehicular air conditioning system. Problems with a wind-up of regulator integrated parts of an evaporator, which can arise if the displacement is not known, can be avoided. Such problems can occur if the compressor is already at its maximum displacement and the requested cooling power is higher than can be provided by the compressor.

Preferably, the compressor is a swash-plate type compressor. However, other variable displacement compressor types are possible.

Only few, preferably three, signals are needed for the modeling on which the preferred method is based. The additional parameters for the model can be easily determined by thermodynamic property tables of the fluid being compressed in the compressor, such as a common refrigerant 1,1,1,2-tetrafluorethane, also known under the trade name R-134a, for example. This is reasonable, if the fluid cannot thermodynamically be treated as an ideal gas. If a fluid is used which can thermodynamically be treated as ideal gas, the parameters can be computed by thermodynamic formulas containing constant thermal capacities as known in the art and which can be found in thermodynamic textbooks.

The state on the low pressure side of the compressor can be determined with the help of an ambient temperature, particular an air temperature, downstream after an evaporator connected to the compressor inlet. The state on the high pressure side can be determined with the help of a measured pressure and the assumption that the compressor work is isentropic with a very high efficiency

Additionally, to estimate the torque, the engine speed has to be known. It is not necessary to know the mass flow or the actual swash-plate angle, which parameters are difficult to determine in praxis. Favorably, the method can be used quite general, which means that the estimation of the compressor torque does not depend on the ambient in which the compressor is utilized. For example, it can be utilized in various vehicles or different engines. Further, there is no influence on the model in different climes, so the model can be applied to compressors which are used in an air conditioning system in different countries all around the world. If a compressor of similar type is used, there is no need to build a new model. Developing a new model is very time consuming. Therefore, a general model such as used in the preferred method is cost-efficient and time-saving.

Only a few input signals are necessary. Thus, the model is easy to analyze. Furthermore, only a few sensors are necessary which allow for an economic air conditioning system, particularly in high-volume series which are typical in vehicle production.

Favorably, said first thermodynamical state on the low pressure side is extracted at least from a temperature, characteristic of a fluid temperature at an inlet of the compressor. Preferably, the temperature is an air temperature measured at an evaporator being in fluid connection with said inlet. From said temperature an input pressure value of said fluid can be determined. It might be possible to combine one of the temperature and pressure with another appropriate measure, for instance air flow or the like, to get an improved estimate.

A specific enthalpy, a specific volume, a specific internal energy and a specific entropy of said fluid at said inlet can be determined, depending on said temperature value and said pressure value. These parameters can be extracted from thermodynamic property tables for the fluid, e.g. the refrigerant being compressed in the compressor.

Favorably, said second thermodynamical state on the high pressure side is extracted from a pressure value of the fluid at an outlet on the high pressure side of said compressor.

Preferably, the pressure value is measured at the outlet of the compressor. An entropy at said outlet is set equal or at least approximately equal to an entropy at the inlet. An enthalpy can be determined depending on said pressure value and said entropy. An isentropic efficiency can also be calculated.

An outlet temperature value of said fluid at said outlet is determined depending on said enthalpy and said pressure value.

Finally, a specific enthalpy, a specific volume and a specific internal energy of said fluid at said outlet can be determined, depending on said temperature value and said pressure value. The parameters can be drawn from thermodynamic property tables describing the thermodynamic properties of the fluid compressed by the compressor.

Said displacement of the compressor can easily be determined as a product of a volumetric efficiency and said compression ratio and a constant characterizing a constant volume in the size of the minimum displacement.

Said compressor torque can be extracted from an energy balance on the compressor.

The torque can be extracted from said calculated displacement, said specific enthalpy at said inlet, said specific enthalpy at said outlet, said specific energy at said inlet and said specific energy at said outlet. The torque can be compensated for operating said compressor efficiently, particularly in a vehicle.

According to another aspect of the invention, a system comprising a variable displacement compressor is provided, wherein an evaporator is in fluid connection with an inlet of said compressor and a condenser in fluid connection with an outlet of said compressor. The compressor is coupled to an operation unit capable of operating said compressor according to a method according to an aspect of the present invention, providing a compressor torque depending on a calculated compressor displacement.

Preferably, a temperature sensor is coupled to the evaporator. By this, an input temperature of the fluid being compressed in the compressor, particularly a refrigerant such as tetrafluorethane, can be measured. Favorably, a pressure on the high pressure side of the compressor can be measured with a pressure sensor which is coupled to said outlet of said compressor. The compressor is driven by a motor, such as an electric motor, or an engine, particularly a vehicle engine, and the engine speed can easily be measured with an appropriate speed sensor. With these signals, the displacement of the compressor can be calculated and a torque of the compressor estimated.

In a preferred embodiment, the compressor is in fluid connection with an automotive based air conditioning system.

According to another aspect of the invention, a computer program is proposed, comprising a computer program code adapted to perform a method for use in the preferred method according to any of the aspects described above, when said program is run on a programmable microcomputer.

The computer program can favorably be adapted to be downloaded to a control unit or one of its components when run on a computer which is connected to the internet.

Preferably, a computer program product stored on a computer readable medium, comprising a program code for use in a method according to anyone of the aspects of the preferred method on a computer.

According to another aspect of the invention, a computer program product stored on a computer readable medium is proposed, comprising a program code for use in a method comprising at least the steps of estimating a first thermodynamical state of said low pressure side and a second thermodynamical state of said high pressure side, estimating a current compression ratio of said compressor, calculating a compressor displacement as a function of said compression ratio, estimating a compressor torque depending on said calculated displacement, and providing said compressor torque to a control unit.

BRIEF DESCRIPTION OF THE DRAWINGS

The present invention together with the above-mentioned and other advantages may best be understood from the following detailed description of the embodiments, but not restricted to the embodiments, wherein is shown schematically:

FIGS. 1 a-1d various states during an ideal compression cycle;

FIG. 2 an ideal pressure-volume diagram according to FIGS. 1a-1d;

FIG. 3 schematically a preferred air conditioning system with a sketch of a preferred compressor;

FIG. 4 a temperature-entropy diagram according to the ideal compression cycle in FIG. 2;

FIG. 5 enthalpy-entropy curves for an isentropic compressor and a real compressor;

FIG. 6 a top view of the compressor in FIG. 3; and

FIG. 7 a summary of a modeling on which a preferred method is based.

DETAILED DESCRIPTION

In the drawings, equal or similar elements are referred to by equal reference numerals.

The following propositions from the thermodynamics have been used in the modeling on which the preferred method is based:

Proposition 1: The thermodynamic state of a compressible system can be completely specified if two independent thermodynamic properties are known.

Proposition 2: Energy flow entering the system minus energy flow leaving the system is equal to the change in the total energy of the system, i.e.

$\begin{array}{cc}{\stackrel{.}{E}}_{\mathrm{in}}-{\stackrel{.}{E}}_{\mathrm{out}}={\left(\frac{dE}{dt}\right)}_{\mathrm{system}}& \left(\mathrm{e1}\right)\end{array}$

The terms in the energy balance contain energy forms according to,

$\begin{array}{cc}{\stackrel{.}{E}}_{\mathrm{in}}={\stackrel{.}{Q}}_{\mathrm{in}}+{\stackrel{.}{W}}_{\mathrm{in}}+\stackrel{.}{m}\left(h_{\mathrm{in}}+{\mathrm{gz}}_{\mathrm{in}}+\frac{{\upsilon }_{\mathrm{in}}^2}{2}\right)& \left(\mathrm{e2}\right)\\ {\stackrel{.}{E}}_{\mathrm{out}}={\stackrel{.}{Q}}_{\mathrm{out}}+{\stackrel{.}{W}}_{\mathrm{out}}+\stackrel{.}{m}\left(h_{\mathrm{out}}+{\mathrm{gz}}_{\mathrm{out}}+\frac{{\upsilon }_{\mathrm{out}}^2}{2}\right)& \left(\mathrm{e3}\right)\\ {\left(\frac{dE}{dt}\right)}_{\mathrm{system}}={\left(\frac{d\left(\mathrm{mu}\right)}{dt}\right)}_{\mathrm{system}}+{\left(\frac{d\left(\mathrm{mgz}\right)}{dt}\right)}_{\mathrm{system}}+{\left(\frac{d\left(m\frac{{\upsilon }^2}{2}\right)}{dt}\right)}_{\mathrm{system}}+{\left(\frac{d\left(I\frac{{\omega }^2}{2}\right)}{dt}\right)}_{\mathrm{system}}& \left(\mathrm{e4}\right)\end{array}$

For clarity, the parameters used in the above-quoted and in the following equations (e1) -(e31) are listed below:

D	displacement	[m3]
E	energy	[J]
Ė	power	[W]
g	gravitational acceleration	[m/s2]
h	specific enthalpy	[J/kg]
I	inertia	[kgm2]
k	gear ratio
K	volume factor	[m3]
m	mass	[kg]
{dot over (m)}	mass rate	[kg/s]
M	compressor torque	[Nm]
n	number of cylinders in the compressor
P	pressure	[Pa]
Q	heat	[J]
{dot over (Q)}	heat power	[W]
s	specific entropy	[kJ/kgK]
t	time	[s]
T	temperature	[° C.]
Tevap	air temperature after evaporator	[° C.]
Tsat	saturation temperature at low pressure	[° C.]
Tsh	superheat of tetrafluorethane before compression	[° C.]
u	specific internal energy	[J/kg]
v	specific volume	[m3/kg]
V	volume	[m3]
ωcomp	compressor speed	[rad/s]
ωengine	engine speed	[rad/s]
W	work	[J]
{dot over (W)}	work power	[W]
z	height related to a specific reference level	[m]
ηs	isentropic efficiency
ηvol	volumetric efficiency
ρ	density	[kg/m3]
υ	speed	[m/s]

with the indices

in  related to compressor inlet (12 in FIG. 3)
out related to compressor outlet (16 in FIG. 3)

The method is described for a compressor (10 in FIG. 3) where a fluid is compressed for an automotive air conditioning system. As fluid a refrigerant tetrafluorethane, also known as R-134a, is preferred, which can easily be liquefied at high pressure.

An ideal compression cycle undergoes four processes that are distinguished by four important states illustrated in FIGS. 1a-1d. FIG. 2 illustrates a complete compression cycle in an ideal pressure-volume (PV)-diagram known from thermodynamics. Point 1a refers to FIG. 1a, point 1b to FIG. 1b, point 1c to FIG. 1c and point 1d to FIG. 1d.

In a cylinder 30, a piston 32 can move up and down along a cylinder axis. Depending on the axial position of the piston 32, a volume of a compression chamber 34 above the piston 32 is varied. In the compression chamber 34 the fluid is contained. The compression chamber 34 is opened or closed by a valve 36. In FIG. 1a, the fluid has been taken in at the low pressure side of the compressor and the compression chamber 34 has its maximum volume (point 1a in FIG. 2) at minimum pressure (P_in in FIG. 2). In this position the pressure valve 36 is closing (indicated by an arrow in FIG. 1a) and the piston 32 starts to move upwards. As a result, the volume of the compression chamber 34 decreases with increasing pressure in the compression chamber 34. When the pressure in the cylinder 30 gets equal (FIG. 1b; point 1b in FIG. 2) to a predetermined maximum pressure (P_out in FIG. 2), the valve 36 opens (indicated by an arrow in FIG. 1b) and the fluid squeezes out, and the piston 32 continues to go up thereby further decreasing the volume of the compression chamber 34 at constant maximum pressure (P_out in FIG. 2) until the valve 36 closes (indicated by an arrow in FIG. 1c). The cylinder 32 has now reached its minimum volume 34 (FIG. 1c; point 1c in FIG. 2). When the piston 32 goes downwards, the volume of the compression chamber 34 increases and the fluid expands until the pressure gets equal to the predetermined minimum pressure (P_in in FIG. 2) at the low pressure side of the compressor (FIG. 1d; point 1d in FIG. 2). Now the valve 36 is opening again (indicated by an arrow in FIG. 1d) and the piston 32 continues to go down thereby further increasing the volume of the compression chamber 34 with constant minimum pressure (P_in in FIG. 2) until the maximum volume of the compression chamber 34 is reached (FIG. 1a, point 1a in FIG. 2) and the compression cycle is closed. A double-headed arrow symbolizes the variable displacement of the compressor D in FIG. 1a, which, for instance, can be varied by varying an inclination angle and/or a linear movement of a swash plate (not shown).

FIG. 3 depicts schematically a preferred air conditioning system 58 comprising a compressor 10 (simplified illustration) with an inlet 12 (suction side) on a low pressure side 24 and an outlet 16 (discharge side) on a high pressure side 26. The inlet 12 forms a suction pressure zone and the outlet 16 forms a discharge pressure zone. Preferably the air conditioning system 58 is a vehicle air conditioning system.

The details of the compressor 10 are well known to a person skilled in the art and therefore not described explicitly. Upstream of the inlet 12 an evaporator 50 is arranged which is in operational fluid connection with the inlet 12. A temperature sensor 54 is mounted at the evaporator 54. Downstream of the outlet 16 a condenser 52 is provided in fluid connection with the outlet 16. A pressure sensor 56 is coupled to the outlet 16. Bypass conduits 14 and 18 are connected to the inlet 12 and the outlet 16, respectively. Unit 60 is a control unit controlling evaporator 50 and condenser 52. In the evaporator 50, the fluid to be compressed in the compressor 10 is evaporated, and in the condenser 52, the compressed fluid is cooled after compression. Preferably, unit 60 is an engine control unit operating the vehicle engine (not shown).

The compressor 10 contains several cylinders 30a, 30b, 30c, 30d and the same number of pistons 32a, 32b, 32c, 32d that all can rotate around a shaft 22 that transfers the torque. The cylinders 30a-30d are supplied with suitable valves (not shown) located at the inlet 12 and at the outlet 16. The valves are opened and closed in accordance with the compression cycle as depicted in FIG. 1 and FIG. 2.

The shaft 22 is coupled to the vehicle engine (not shown) and driven by the engine rotation. The refrigerant, which enters the compressor 10 as fluid through the inlet 12 enters the respective cylinder 30a . . . 30d at the inlet 12. In FIG. 3, as an example, cylinder 30a receives the fluid on the low pressure side 24. Then the cylinder 30a is rotated by the shaft 22 to the outlet 16. When the shaft 22 rotates, the piston 32a glides with its foot 38a on a swash plate 20 which is inclined compared with the axis of shaft 22. The swash plate 20 is used as a cam plate forcing the piston 32a upwards during the rotational movement of the cylinder 30a from the low pressure side 24 to the high pressure side 26 of the compressor 10 and downwards during the rotational movement of the cylinder 30a from the high pressure side 26 to the low pressure side 24 of the compressor 10.

This means that the inclined swash plate 20 converts rotation into a reciprocation of each piston 32a, 32b, 32c, 32d, which in turn alters correspondingly the volume of the compression chamber 34a, 34b, 34c, 34d for the fluid. For instance, the position of the cylinder 30a relative to the inclined swash plate 20 defines the stroke of its piston 32a, i.e. the displacement D of the compressor chamber 34a. Taking a situation as depicted in FIG. 3 where the piston 32a is at its lower inversion position with the maximum volume of the compression chamber 34a as a starting point, during rotation of the cylinder 30a the volume of the compression chamber 34a continuously decreases in the cylinder 30a and the fluid is continuously compressed until a predetermined maximum pressure when the valve at the outlet 16 opens and the fluid leaves the compression chamber 34a through the outlet 16 on the high pressure side 26 of the compressor 10. Continuing the rotational movement of the cylinder 30a the volume of the compression chamber 34a further decreases at constant maximum pressure until the compression chamber 34a has reached its minimum volume (and the piston 32a its upper inversion position). Further continuing the rotational movement of the cylinder 30a the volume of the compression chamber 34a increases again in the cylinder 30a, and the fluid is expanded thereby decreasing the pressure of the fluid until a predetermined minimum pressure when the valve at the inlet 12 opens and additional fluid is taken in into the compression chamber 34a through the inlet 12 on the low pressure side 26 of the compressor 10. Further continuing the rotational movement of the cylinder 30a the volume of the compression chamber 34a further increases at constant minimum pressure until the compression chamber 34a has reached its maximum volume (and the piston 32a has reached its lower inversion position as depicted in FIG. 3). This movement is recurrently repeated for each of the cylinders 30a-30d, i.e. each of the cylinders 30a-30d is subject to this rotation of the cylinder 30a-30d around the shaft 22 and caused by the rotation of the shaft 22 and the corresponding reciprocal linear movement of its piston 32a-32d caused by the inclined swash plate 20.

The inclination of the swash plate 20 can be changed by pivoting it around its joint 20a, e.g. dependent on operation conditions of the air conditioning system 58 the displacement of the compressor 10 can be altered.

Only three signals need to be used as input signals for the modelling of the behaviour of the compressor:

• • pressure P_out downstream after compressor 10,
  • air temperature T_evap downstream after evaporator 50,
  • rotational engine speed ω_engine, wherein the compressor 10 is driven by the vehicle engine via shaft 22. It is a favourable benefit of the invention that it is not necessary to use a mass flow measure or a displacement measure, other than what has been used in art of air condition engineering.

The compressor torque M is computed in three steps. Firstly, the thermodynamic states are determined on both low and high pressure sides 24, 26 of the compressor 10. Then the displacement D is calculated and finally the torque M is computed by determining an energy balance on the compressor 10.

As the typical refrigerants, for instance the preferred refrigerant tetrafluorethane, usually cannot be treated as an ideal gas, and their thermodynamic properties usually cannot be computed by formulas containing constant thermal capacities, thermodynamic property tables for said refrigerants (as for instance tetrafluorethane) are preferably being used. Such tables can be found for example in various textbooks for the respective fluid used in the compressor 10. If two independent properties are given, a third arbitrary property can be determined from these tables according to proposition 1 as described before.

The sensor 54 measuring the ambient, particularly air temperature T_evap after the evaporator 50 is used to determine the state upstream of the compressor 10. Since the fluid cools the air through heat exchange in the evaporator 50, the air and the fluid have approximately the same temperature downstream the evaporator 50 i.e.T_in=T_evap  (e5)

On its way to the compressor 10 the temperature T_in of the fluid changes only marginally. This small change can be neglected.

By approximating the pressure at the low pressure side 24 upstream of the compressor 10 as the saturation pressure, this property is determined. This can be done even though there is some superheating going on, as shown in FIG. 4.

FIG. 4 depicts functional dependencies of temperature T vs. entropy s (T-s-diagram) known from thermodynamics for a compression process as indicated in FIGS. 1a-1d. Curve A represents a saturation curve, while curves B and C represent curves at constant pressure, i.e. minimum pressure P_in (curve C) at the inlet 12 of the compressor 10 and maximum pressure P_out (curve B) at the outlet 16 of the compressor 10, as described in FIG. 3.

In the saturated case, as indicated by curve A, the temperature T increases monotonically with increasing entropy s until a maximum temperature T is reached at a maximum of the curve A and drops monotonically behind the maximum with increasing entropy s. On the saturation curve A all the thermodynamic properties can be determined if one of them is known.

The area representing superheating is located to the right of the saturation curve A, and in this area the compression of the fluid in the compressor 10 is done. Superheating (sometimes referred to as boiling retardation, or boiling delay) is the phenomenon in which a liquid is heated to a temperature higher than its standard boiling point, without actually boiling. This can be caused by rapidly heating a homogeneous substance while leaving it undisturbed (so as to avoid the introduction of bubbles at nucleation sites).

For the fluid at minimum pressure P_in the temperature T increases with increasing entropy s until a first intersection point of curves C and A at low entropy s and is constant beneath the saturation curve A. In this region, the pressure P_in is constant. For higher entropy values, the curve C intersects curve A at a second intersection point at a higher entropy s and starts to rise again.

Curve B represents the fluid at maximum pressure P_out which shows qualitatively the same behavior at a higher temperature level. Due to the higher temperature T for the maximum pressure P_out, the intersection points between curves A and B are located in between the intersection points of curves A and C. The region where pressure P_out and temperature T are constant is accordingly smaller than for the minimum pressure P_in.

The amount of superheat T_sh, which is indicated in FIG. 4, is around 10 K, and represents the increase in temperature T when compressing the fluid during the upward movement of the piston 32 (see for instance FIGS. 1a-1d). As mentioned above, the compression of the fluid is done in the superheated area located to the right of the saturation curve A. The compression process experienced by the fluid between the minimum pressure P_in at the inlet 12 to the maximum pressure P_out at the outlet 16 of the compressor 10 is indicated by a line between curves B and C, which line is slightly curved towards higher entropy s, indicating not only a rise in temperature T but also a slight increase in entropy s during compression of the fluid.

The thermodynamic state is completely determined by a single property if the fluid is a saturated gas, and in this case it is sufficient to know the saturation temperature T_sat given byT_sat=T_in−T_sh  (e6)

The pressure P_in at the inlet 12 is then determined from tables listing saturated fluid properties (as for instance for a refrigerant as tetrafluorethane (R-134a) properties),

$\begin{array}{cc}{T}_{\mathrm{sat}}\underset{\underset{\underset{R-134a}{\mathrm{saturated}\mathrm{fluid}.e.g.}}{\mathrm{tables}}}{\underbrace{\Rightarrow }}{P}_{\mathrm{in}}& \left(\mathrm{e7}\right)\end{array}$

FIG. 4 is reflecting this relation, wherein on the right edge of the saturation curve A there is only one pressure value for P_in corresponding to a specific temperature T at the intersection point of curves A and C.

With two known properties, the remaining ones are found in property tables covering the fluid in the compressor 10, for example superheated tetrafluorethane.

$\begin{array}{cc}\begin{array}{c}{T}_{\mathrm{in}}\\ p_{\mathrm{in}}\end{array}\}\underset{\underset{\underset{\underset{\underset{e.g.R-134a}{\mathrm{fluid},}}{\mathrm{superheated}}}{\mathrm{tables}}}{︸}}{\Rightarrow }\begin{array}{c}{h}_{\mathrm{in}}\\ v_{\mathrm{in}}\\ u_{\mathrm{in}}\\ s_{\mathrm{in}}\end{array}\}& \left(\mathrm{e8}\right)\end{array}$

Typically, compressors work almost isentropic and this can be used to determine the state at the outlet 16 of the compressor 10. Two possible compressor processes I-2s and I-2 are shown in FIG. 5 in an h-s diagram indication the functional dependence of the specific enthalpy h from the specific entropy s known in thermodynamics. The upper curve reflects the enthalpy-dependence on the entropy s of the maximum pressure P_out at the outlet 16 and the lower curve the enthalpy-dependence on the entropy s of the minimum pressure P_in at the inlet 12 of the compressor 10. Both curves increase monotonically in specific enthalpy h with increasing specific entropy s, revealing a higher specific enthalpy h for the higher pressure P_out compared to the lower pressure P_in.

Two paths, I-2s and I-2, are drawn in showing the compression process from minimum pressure P_in to maximum pressure P_out. The left path I-2s going from the lower curve P_in at a lower enthalpy h to the upper curve P_out with a higher enthalpy at constant entropy s is for an isentropic compressor 10, the right path I-2 is for a real compressor 10. If the compression is done without any increase in entropy s, which is hard to achieve in reality, the process path would look like the I-2s line. As a result of losses by friction, the entropy s will increase somewhat making the path look more like I-2 in FIG. 5, where an increase of the enthalpy h by moving from P_in to P_out results in an increase in entropy s in the fluid. The latter path I-2 corresponds to the real compressor process.

The compressor work is equal to the enthalpy difference over the compressor 10 when potential and kinetic energy are neglected. If the isentropic efficiency n_s is defined as isentropic compressor work divided by real compressor work this results in

$\begin{array}{cc}{\eta }_s=\frac{h_{2s}-h_1}{h_2-h_1}& \left(\mathrm{e9}\right)\end{array}$ with indices according to the states in FIG. 5. Indices 1 and 2 correspond to the indices in and out used elsewhere.

Equation (e8) is used to determine the entropy s₁ at the compressor inlet 12. The entropy s_2s is defined according to FIG. 5 and thuss_2s=s₁  (e10)

The pressure P_out at the outlet 16 is measured with the pressure sensor 56, and with two properties at 2s, the enthalpy h at that point can be determined,

$\begin{array}{cc}\begin{array}{c}{s}_{2s}\\ p_{\mathrm{out}}\end{array}\}\underset{\underset{\underset{\underset{\underset{R-134a}{\mathrm{fluid}.e.g.}}{\mathrm{superheated}}}{\mathrm{tables}}}{︸}}{\Rightarrow }{h}_{2s}& \left(\mathrm{e11}\right)\end{array}$

Another way to write equation (e9) is

$\begin{array}{cc}{h}_2=\frac{h_{2s}-h_1}{{\eta }_s}+h_1& \left(\mathrm{e12}\right)\end{array}$

It is appropriate to put n_s=0.9 for this type of compressor 10. Finally the temperature T_out is determined through the property tables,

$\begin{array}{cc}\begin{array}{c}{h}_2\\ p_{\mathrm{out}}\end{array}\}\underset{\underset{\underset{\underset{R-134a}{\mathrm{fluid}.e.g.}}{\mathrm{superheated}}}{\mathrm{tables}}}{\underbrace{\Rightarrow }}{T}_{\mathrm{out}}& \left(\mathrm{e13}\right)\end{array}$

Remaining property values can be collected from the tables similar to equation (e8),

$\begin{array}{cc}\begin{array}{c}{T}_{\mathrm{out}}\\ p_{\mathrm{out}}\end{array}\}\underset{\underset{\underset{\underset{\underset{R-134a}{\mathrm{fluid}.e.g.}}{\mathrm{superheated}}}{\mathrm{tables}}}{︸}}{\Rightarrow }\begin{array}{c}{h}_{\mathrm{out}}\\ v_{\mathrm{out}}\\ u_{\mathrm{out}}\end{array}\}& \left(\mathrm{e14}\right)\end{array}$

The compression ratio of the fluid in the compressor 10 is used to estimate the variable displacement D. This ratio is defined as the quotient of the fluid density at the outlet 16 and the fluid density at the inlet 12 of the compressor 10.

FIG. 6 is used to relate the compressor 10 in FIG. 3 to the ideal compression cycle described in FIG. 2 and FIG. 1a-1d. It shows a view of the compressor 10 as depicted in FIG. 3 from above. A specific cylinder 30a-30d going through a compression cycle can be studied by comparing the Figures. The valves 36a-36d are closing at Ia and Ic (FIG. 2) at maximum and minimum volume of the compression chambers 34a-34d, respectively. The valves 36a-36d are opening at Ib and Id (FIG. 2) whose location vary somewhat as a result of different high and low pressures in the high and low pressure sides 26, 24 of the compressor 10. The volume of the compression chamber 34a-34d decreases linearly in the direction of an axis (not shown) going from Ia to Ic (FIG. 2). The rotating cylinders 30a-30d are exposed to the high pressure below line 28 and exposed to the low pressure above the line 28 in the FIG. 6. Line 28 indicates a physical barrier between the low and the high pressure sides 24, 26 of the compressor 10.

It is assumed that just after the cylinder valve 36 (FIGS. 1a-1d) has been shut in Ia (FIG. 2), the valve 36 is exposed to the high pressure and a half revolution later, just after Ic (FIG. 2), the valve 36 is exposed to the low pressure again. The locations of the states Ib and Id vary and depend on existing pressures. State Ib occurs later, and with smaller volume the larger P_out is, while state Id occurs earlier, and with smaller volume the larger P_in is.

Letting V_i denote the cylinder volume of the compression chamber 34a-34d at state i (with i corresponding to Ia, Ib, Ic, Id, see FIGS. 1a-1d and FIG. 2) and n be the number of cylinders 30a-30d in the compressor 10, the displacement D can be written asD=(V_a−V_d)n  (e15)

The mass m_in entering a cylinder 30a-30d at the compressor inlet 12 between Id and Ia (FIG. 2), and the mass m_out leaving the same cylinder 30a-30d later at the compressor outlet 16 between Ib and Ic (FIG. 2), arem_in=(V_a−V_d)ρ_in  (e16)m_out=(V_b−V_c)ρ_out  (e17)

where ρ denotes density of the fluid. Since there is no accumulation of mass these two masses m_in, m_out must be equal,

$\begin{array}{cc}{V}_a-V_d=\frac{{\rho }_{\mathrm{out}}}{{\rho }_{\mathrm{in}}}\left(V_b-V_c\right)& \left(\mathrm{e18}\right)\end{array}$

Swapping densities for specific volumes v and using (e15) and (e18) give

$\begin{array}{cc}D=\frac{v_{\mathrm{in}}}{v_{\mathrm{out}}}\left(V_b-V_c\right)n& \left(\mathrm{e19}\right)\end{array}$

The specific volumes are given by equations (e8) and (e14), n is known and V_c is constant since a change in the inclination of the swash plate 20 does not change the location of the piston 32a-32d closest to the joint 20a of the plate 20, see FIG. 3.

Holding the compressor rotational speed ω at a constant level, V_b is a function of pressure as well as of the inclination of the swash plate 20. The function is not exactly known but can be estimated. Assume the inclination of the swash plate 20 has just increased by pivoting it downward around its joint 20a. If the valve 36 at Ib (FIG. 2) is opening at the same location as before V_b would increase. On the other hand, a larger inclination implies larger compression ratio and higher pressure which causes the valve 36 to open later than before, i.e. Ib moves towards Ic (FIG. 2). The conclusion of this is that V_b increases as a result of increasing inclination, but decreases as a result of increasing high pressure. For that reason, V_b is approximated as a constant if the compressor rotational speed ω is constant, and this means that V_b is a function of compressor rotational speed ω only.

Merging the constants and V_b in equation (e19) to η_volK gives

$\begin{array}{cc}D={\nu }_{\mathrm{vol}}K\frac{v_{\mathrm{in}}}{v_{\mathrm{out}}}& \left(\mathrm{e20}\right)\end{array}$ where K is a constant volume in the size of the minimum displacement D,

$\begin{array}{cc}K=\frac{\left(V_b-V_c\right)n}{{\eta }_{\mathrm{vol}}}& \left(\mathrm{e21}\right)\end{array}$ and η_vol is the volumetric efficiency.

If the inclination of the swash plate 20 is constant, equation (e21) contains constants only apart from η_vol and V_b. These variables are directly proportional with respect to each other according to the equation. The valve 36 (FIGS. 1a-1b) is opening in Ib (FIG. 2), but the opening process is taking some time and a high compressor rotational speed ω would mean that the cylinder 30a-30d will rotate a longer distance during this time. This way V_b is decreasing which means that η_vol is decreasing as well. The volumetric efficiency η_vol can be assumed to decrease linearly from 1 to 0.5 with increasing compressor rotational speed ω.

The main equation of the preferred model is based on the principle of energy conservation presented in equation (e1). In its complete form it looks like

$\begin{array}{cc}{\stackrel{.}{Q}}_{\mathrm{in}}+{\stackrel{.}{W}}_{\mathrm{in}}+{\stackrel{.}{m}}_{\mathrm{in}}\left(h_{\mathrm{in}}+\frac{{\upsilon }_{\mathrm{in}}^2}{2}+{\mathrm{gz}}_{\mathrm{in}}\right)-{\stackrel{.}{Q}}_{\mathrm{out}}-{\stackrel{.}{W}}_{\mathrm{out}}-{\stackrel{.}{m}}_{\mathrm{out}}\left(h_{\mathrm{out}}+\frac{{\upsilon }_{\mathrm{out}}^2}{2}+{\mathrm{gz}}_{\mathrm{out}}\right)=\frac{d\left(m_{\mathrm{system}}\left(u_{\mathrm{system}}+\frac{{\upsilon }_{\mathrm{system}}^2}{2}+{\mathrm{gz}}_{\mathrm{system}}\right)+\frac{I_{\mathrm{system}}{\omega }_{\mathrm{system}}^2}{2}\right)}{dt}& \left(\mathrm{e22}\right)\end{array}$ where the system is taken as the compressor 10. For a preferred embodiment of the system according to the invention a number of simplifications are done,

• • Adiabatic process. {dot over (Q)}_in={dot over (Q)}_out=0
  • Potential and kinetic flow energy are neglected

$h>>\frac{{\upsilon }^2}{2}+\mathrm{gz}.$

• • The system does not do any work on the surrounding. {dot over (W)}_out=0
  • No accumulation of mass in the system. {dot over (m)}_in={dot over (m)}_out={dot over (m)}
  • The position of the compressor is fixed, thus no change in kinetic or potential energy.

$\begin{array}{cc}\Rightarrow \frac{d\left(m_{\mathrm{system}}\left(\frac{{\upsilon }_{\mathrm{system}}^2}{2}+{\mathrm{gz}}_{\mathrm{system}}\right)\right)}{dt}=0& \left(\mathrm{e23}\right)\end{array}$

The first simplification is relatively gross since there is some heat generated in the compressor 10. The energy balance according to equation (e22) now turns into

$\begin{array}{cc}{\stackrel{.}{W}}_{\mathrm{in}}+\stackrel{.}{m}\left(h_{\mathrm{in}}-h_{\mathrm{out}}\right)={\left(\frac{d\left(\mathrm{mu}+\frac{I{\omega }_{\mathrm{comp}}^2}{2}\right)}{dt}\right)}_{\mathrm{system}}& \left(\mathrm{e24}\right)\end{array}$

Since the provided power is equal to the compressor rotational speed times the compressor torque M, equation (e24) becomes

$\begin{array}{cc}{w}_{\mathrm{comp}}M+\stackrel{.}{m}\left(h_{\mathrm{in}}-h_{\mathrm{out}}\right)={\left(\frac{d\left(\mathrm{mu}+\frac{I{\omega }_{\mathrm{comp}}^2}{2}\right)}{dt}\right)}_{\mathrm{system}}& \left(\mathrm{e25}\right)\end{array}$

The mass flow {dot over (m)} is a function of compressor rotational speed ω_comp, displacement D and density ρ according to

$\begin{array}{cc}\stackrel{.}{m}=\frac{{\omega }_{\mathrm{comp}}}{2\pi }D{\rho }_{\mathrm{in}}=\frac{{\omega }_{\mathrm{comp}}}{2\pi }D\frac{1}{v_{\mathrm{in}}}& \left(\mathrm{e26}\right)\end{array}$ where ω_comp is the compressor rotational speed in radians per second and D is the displacement in cubic meters per revolution.

While the compressor 10 is connected to the engine a constant gear ratio is achieved according toω_comp=kω_engine  (e27)where ω_engine is the engine rotational speed and k is the gear ratio. Both are known.

By solving equation (e25) for M and using the mass flow from equation (e26) and the compressor rotational speed ω_comp from equation (e27) the torque M can be expressed as

$\begin{array}{cc}M=\frac{D}{2\pi v_{\mathrm{in}}}\left(h_{\mathrm{out}}-h_{\mathrm{in}}\right)+\frac{1}{k{\omega }_{\mathrm{engine}}}{\left(\frac{d\left(\mathrm{mu}+\frac{{\mathrm{Ik}}^2{\omega }_{\mathrm{engine}}^2}{2}\right)}{dt}\right)}_{\mathrm{system}}& \left(\mathrm{e28}\right)\end{array}$

Replacing the displacement D with equation (e20) gives

$\begin{array}{cc}M=\frac{K}{2\pi v_{\mathrm{out}}}\left(h_{\mathrm{out}}-h_{\mathrm{in}}\right)+\frac{1}{k{\omega }_{\mathrm{engine}}}{\left(\frac{d\left(\mathrm{mu}+\frac{{\mathrm{Ik}}^2{\omega }_{\mathrm{engine}}^2}{2}\right)}{dt}\right)}_{\mathrm{system}}& \left(\mathrm{e29}\right)\end{array}$ where the state properties are given by equations (e8) and (e14).

m_system is the mass of fluid in the compressor 10 and equals the displacement D multiplied by the density at the compressor inlet 12 from equation (e20) yields

$\begin{array}{cc}{m}_{\mathrm{system}}=D\frac{1}{v_{\mathrm{in}}}={\eta }_{\mathrm{vol}}K\frac{1}{v_{\mathrm{out}}}& \left(\mathrm{e30}\right)\end{array}$

The specific internal energy u can not be determined since the fluid in the compressor 10 is not in thermal equilibrium. An approximation is made using the average value of the internal energy at the inlet 12 and outlet 16 of the compressor 10.

$\begin{array}{cc}{u}_{\mathrm{system}}=\frac{u_{\mathrm{in}}+u_{\mathrm{out}}}{2}& \left(\mathrm{e31}\right)\end{array}$

As a summary, FIG. 7 depicts the modeling on which the preferred method is based. FIGS. 2, 4 and 5 depict basic thermodynamical relations for a compression cycle depicted in FIG. 1. By measuring three signals T_evap, P_out, and ω_engine. Knowing T_evap and P_out of the fluid, the parameters T_in, P_in and T_out of the fluid can be derived from equations (e5), (e7) and (e13). From the latter three parameters by applying equations (e8) and (e14) the specific volumes v_in and v_out, as well as the specific entropies h_in,h_out and energies u_in,u_out are derived. Applying the equation (e20) to the specific volumes v_in and v_out yields the actual displacement D of the compressor 10, and applying equation (e29) to the displacement D, the specific entropies h_in, h_out and specific energies u_in,u_out, the torque M of the compressor 10 is calculated. The compressor rotational speed 107_comp can easily be derived from the engine rotational speed ω_engine with equation (e27).

By calculating the actual displacement D, the actual temperature of the evaporator 50 (FIG. 3) can be controlled better as well as the control current of the control unit 60 for controlling the evaporator 50 (FIG. 3).

The invention is not restricted to the preferred embodiments described above. For instance, it is possible to use a variable displacement compressor without a swash plate or with a rotating swash plate instead of rotating cylinders and the like. The compressor displacement can be varied not only by variation of the inclination of the swash plate, influencing the active stroke of the cylinders 30a-30d, but also by linear movement of the swash plate away from or towards the cylinders 30a-30d, which influences the maximum shift of the pistons 32a-32d.

Claims

1. A method for operating a variable-displacement compressor, wherein a volume of a compression chamber of at least one cylinder is varied during a compression cycle, wherein during the cycle the at least one cylinder is alternately coupled to a low pressure side and a high pressure side of the compressor, comprising estimating a first thermodynamic state of the low pressure side using a temperature (Tevap)characteristic of a temperature (Tin) of a fluid at an inlet of the compressor and estimating a second thermodynamic state of the high pressure side using a pressure value (Pout) of the fluid at an outlet on high pressure side of the compressor,estimating a current compression ratio (specific volume in/specific volume out

2. The method according to claim 1, wherein the characteristic temperature (Tevap) is an ambient temperature measured at an evaporator in fluid connection with the inlet of the compressor.

3. The method according to claim 2, wherein from the characteristic temperature (Tevap) an input pressure value (Pin) of the fluid is determined.

4. The method according to claim 3, wherein a specific enthalpy (hin), a specific volume (vin), a specific internal energy (uin) and a specific entropy (sin) of the fluid at the inlet are determined, depending on the fluid temperature value (Tin)and the pressure value (Pin) of the fluid.

5. The method according to claim 2, wherein the pressure value of the fluid (Pout) is measured at the outlet.

6. The method according to claim 5, wherein an entropy (s2s) at the outlet is set equal or at least approximately equal to an entropy (s1) at the inlet.

7. The method according to claim 5, wherein an enthalpy (h2s) is determined depending on the pressure value (Pout) of the fluid and the entropy (S2s).

8. The method according to claim 7, wherein an outlet, temperature value (Tout) of the fluid at the outlet is determined depending on the enthalpy (h2s) and the pressure value (Pout) of the fluid.

9. The method according to claim 5, wherein a specific enthalpy (hout), a specific volume (vout)and a specific internal energy (uout) of the fluid at the outlet are determined depending on the temperature value (Tout) and the pressure value (Pout) of the fluid.

10. The method according to claim 1, wherein the displacement (D) is determined as a product of a volumetric efficiency (nvol) and the compression ratio

11. The method according to claim 1, wherein the compressor torque (M) is extracted from an energy balance over the compressor.

12. The method according to claim 1, wherein the torque (M) is extracted from the calculated displacement (D), a specific enthalpy (hin) at the inlet, a specific enthalpy (hout) at the outlet, a specific energy (uin) at the inlet and a specific energy (uout) at the outlet.

13. The method according to claim 1, wherein the torque (M) is compensated for operating an engine coupled to the compressor.

14. A system comprising a variable displacement compressor,an evaporator in operational connection with an inlet of the compressor,a condenser in operational connection with an outlet of the compressor,wherein the compressor is coupled to a control unit arranged to operate the compressor according to a method that causes a compressor torque to be provided depending on a calculated compressor displacement, wherein, in the compressor a volume of a compression chamber of at least one cylinder is varied during a compression cycle and during the cycle the at least one cylinder is alternately coupled to a low pressure side and a high pressure side of the compressor, the method comprisingestimating a first thermodynamic state of the low pressure side using a temperature (Tevap) characteristic of a temperature (Tin) of a fluid at an inlet of the compressor and estimating a second thermodynamic state of the high pressure side using a pressure value (Pout) of the fluid at an outlet on the high pressure side of the compressor,estimating a current compression ratio (specific volume in/specific volume out

15. The system according to claim 14, wherein a temperature sensor is coupled to the evaporator.

16. The system according to claim 14, wherein a pressure sensor is coupled to the outlet of the compressor.

17. The system according to claim 14, wherein the compressor is in fluid connection with an automotive based air conditioning system.

18. A non-transitory computer program product comprising a computer program code adapted to perform a method of for use in a method according to claim 1 when the program is run on a programmable microcomputer.

19. The non-transitory computer program product according to claim 18 adapted to be downloaded to a control unit or one of its components when run on a computer which is connected to the internet, wherein the control unit is coupled to a system comprising a variable displacement compressor, an evaporator in operational connection with an inlet of the compressor, a condenser in operational connection with an outlet of the compressor, wherein the compressor is coupled to the control unit, the control unit being capable of operating the compressor according to the method according to claim 1 , providing a compressor torque (M) depending on a calculated compressor displacement (D).

20. A non-transitory computer program product stored on a computer readable medium, comprising a program code for use in a method according to claim 1 on a computer.

21. A non-transitory computer program product stored on a computer readable medium, comprising a program code for use in a method comprising at least the steps of estimating a first thermodynamic state of the low pressure side using a temperature (Tevap) characteristic of a temperature (Tin) of a fluid at an inlet of the compressor and estimating a second thermodynamic state of the high pressure side using a pressure value (Pout) of the fluid at an outlet on the high pressure side of the compressor,estimating a current compression ratio (specific volume in/specific volume out