WO2013064840A1 - A method of controlling a dynamic physical system that exhibits a chaotic behaviour - Google Patents

A method of controlling a dynamic physical system that exhibits a chaotic behaviour Download PDF

Info

Publication number
WO2013064840A1
WO2013064840A1 PCT/GB2012/052739 GB2012052739W WO2013064840A1 WO 2013064840 A1 WO2013064840 A1 WO 2013064840A1 GB 2012052739 W GB2012052739 W GB 2012052739W WO 2013064840 A1 WO2013064840 A1 WO 2013064840A1
Authority
WO
WIPO (PCT)
Prior art keywords
control
rate
controlled
term
control term
Prior art date
Application number
PCT/GB2012/052739
Other languages
French (fr)
Inventor
Tjeerd Victor Siebe Maria OLDE SCHEPER
Andrew Robert CARNELL
Original Assignee
Oxford Brookes University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Priority to GB1119036.0 priority Critical
Priority to GBGB1119036.0A priority patent/GB201119036D0/en
Application filed by Oxford Brookes University filed Critical Oxford Brookes University
Publication of WO2013064840A1 publication Critical patent/WO2013064840A1/en

Links

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B13/00Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
    • G05B13/02Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
    • G05B13/04Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F02COMBUSTION ENGINES; HOT-GAS OR COMBUSTION-PRODUCT ENGINE PLANTS
    • F02DCONTROLLING COMBUSTION ENGINES
    • F02D41/00Electrical control of supply of combustible mixture or its constituents
    • F02D41/0002Controlling intake air
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F03MACHINES OR ENGINES FOR LIQUIDS; WIND, SPRING, OR WEIGHT MOTORS; PRODUCING MECHANICAL POWER OR A REACTIVE PROPULSIVE THRUST, NOT OTHERWISE PROVIDED FOR
    • F03DWIND MOTORS
    • F03D7/00Controlling wind motors
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F24HEATING; RANGES; VENTILATING
    • F24FAIR-CONDITIONING; AIR-HUMIDIFICATION; VENTILATION; USE OF AIR CURRENTS FOR SCREENING
    • F24F11/00Control or safety arrangements
    • F24F11/30Control or safety arrangements for purposes related to the operation of the system, e.g. for safety or monitoring
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B13/00Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
    • G05B13/02Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
    • G05B13/04Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
    • G05B13/047Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators the criterion being a time optimal performance criterion
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E10/00Energy generation through renewable energy sources
    • Y02E10/70Wind energy
    • Y02E10/72Wind turbines with rotation axis in wind direction

Abstract

A method of controlling a dynamic physical system comprising a plurality of variable quantities. A model of the system comprising a plurality of variables representing the variable quantities, and a plurality of respective rate equations that describe the rate of change of the variables, is obtained. A control term in at least one rate equation from the plurality of rate equations is identified. A rate control function is derived from, for at least one of the variables in the rate equation, the proportion of the variable to the growth rate of the rate equation, and the rate control function is applied to the control term to provide a stabilised control term. The dynamic physical system is then controlled by modifying at least one of the quantities represented by the variables in the control term, so that the control term derived from the modified quantities is substantially the same as the stabilised control term.

Description

A METHOD OF CONTROLLING A DYNAMIC PHYSICAL SYSTEM THAT
EXHIBITS A CHAOTIC BEHAVIOUR
Field of the Invention
The present invention concerns the control of dynamic physical systems. More particularly, but not exclusively, the invention concerns the control of dynamic physical systems that exhibit chaotic behaviour.
Background of the Invention
Chaotic dynamic physical systems are dynamic physical systems that exhibit behaviour that is deterministic in a mathematical sense (the behaviour is precisely determined by the state at any particular time), but nevertheless
unpredictable over time due to the sensitivity of the system on the state. Due to this unpredictability, when attempting to control a dynamic physical system it is generally considered undesirable to allow it to enter a chaotic state, as this limits the ability to control the system.
However, avoiding chaotic states may mean restricting the system to inefficient states. For example, in a chemical reaction, while they are unpredictable the chaotic states may make use of resources more efficiently than the stable, non-chaotic states. Further, in systems beyond a certain size chaotic behaviour may be inevitable, and so chaotic behaviour can be avoided only by restricting to very small-scale systems. For this reason, the control of chaotic systems has been of some interest.
A dynamic physical system can be represented by a plurality of variables, each representing a quantity of the system that varies over time. (For example, if the system is a chemical process, the variables could represent the quantities of the different chemicals or their
concentrations.) The system can then be modelled by a plurality of rate equations for each of the variables, which describe how the quantities vary over time based on the current state.
As is well known, the instantaneous state of a system can be considered to be a point moving around a "state space", where the dimensions of the state space represent the different variables of the system, and the position of the point at a given time is determined by the values of the variables at that time. The state of an example Rossler system (a well-known simple dynamic physical system that exhibits chaotic behaviour) as a point moving through state space is shown in Figure 1.
It is known that for dynamic physical systems in ranges where chaotic behaviour occurs, the global behaviour of the system is dominated by a subset of the rate equations, which tend to be the non-linear parts of the system that allow the quantities of the system to grow at an exponential rate. The local rate of expansion is proportional to the local behaviour of each of the variables. A known method of controlling a chaotic system is the OGY (Ott, Grebogi and Yorke) method, as first described in E. Ott, C. Grebogi, and J. A. Yorke, "Controlling chaos", Phys. Rev. Lett. 64 1196 (1990). The OGY model is based upon the properties of chaotic systems discussed above. In the OGY method, a model of the dynamic physical system is obtained, and analysed to identify an unstable periodic orbit around which the system cycles is identified. The system is then controlled by applying small, pre-determined changes to the variables of the system, in proportion to the local behaviour of the variables, in order to keep the system in or near the chosen orbit.
Another known method of controlling a chaotic system is the Pyragas continuous control method. Similarly to the OGY method, an unstable periodic orbit around which the system cycles is identified. However, in this method adjustments are made to the system in accordance with pre-determined time delays, which need to be carefully matched with the dynamics of the system to allow successful control.
It is a disadvantage of the both OGY method and the Pyragas method that they require the system in operation to be analysed in order to identify a particular unstable periodic orbit, so that the properties of that orbit can be used to implement the method. Further, neither method allows control of systems in all cases.
The present invention seeks to mitigate the above- mentioned problems. Alternatively or additionally, the present invention seeks to provide an improved method of controlling a dynamic physical system, in particular a dynamic physical system that exhibits chaotic behaviour.
Summary of the Invention
In accordance with a first aspect of the invention there is provided a method of controlling a dynamic physical system, wherein the system comprises a plurality of variable quantities, the method comprising the steps of:
obtaining a model of the system, the model comprising a plurality of variables representing the variable quantities, and a plurality of respective rate equations that describe the rate of change of the variables;
identifying a control term in at least one rate equation from the plurality of rate equations;
deriving a rate control function from, for at least one of the variables in the rate equation, the proportion of the variable to the growth rate of the rate equation;
applying the rate control function to the control term to provide a stabilised control term; and
controlling the dynamic physical system by modifying at least one of the quantities represented by the variables in the control term, so that the control term derived from the modified quantities is substantially the same as the stabilised control term.
The control term is a term that in at least some states of the system dominates the control of the variable whose behaviour is defined by the rate equation. As discussed above, in ranges where chaotic behaviour occurs, the local rate of expansion of a variable is proportional to the local behaviour of the variables in the control term. The rate control function is derived from the proportions of the variables, and so acts to stabilise the rate equation. The stabilisation of the model given by the rate control function can be applied to the system itself by modifying the quantities represented by the variables accordingly, thus allowing the system to be stabilised. The control term will in most cases be a term that would under certain conditions expand exponentially if the control provided by the invention were not applied.
Unlike the known methods described above, the method of this aspect of the invention does not require the model to be analysed when in operation, so that the properties of a particular unstable orbit can be identified. Rather, only the rate equations defining the model need to be considered, so that the rate control equation can be derived.
Preferably, the control term includes the variable whose behaviour is described by the rate equation. Such a term usually allows exponential expansion of the system, as the feedback of the value of the variable on its own growth rate is what causes the exponential growth.
Preferably, the proportion qx of a variable x to the growth rate is given by the equation:
where μχ is a constant. Preferably again, the rate control function is of the form: where qxi to qxn are the proportions of the variables Xi to xn to the growth rate, and f and ξ are scalars. The scalars may f and ξ may advantageously be varied so as to stabilise the system into a pre-determined chosen orbit. The rate control function may include a biasing term, in the form of a constant term added to the original rate control function above .
Preferably, the at least one rate equation describes exponential growth of its respective variable. This is because a rate equation allowing exponential growth will tend to dominate the behaviour of the system in at least some ranges. Preferably, the control term contributes to the growth of the respective variable of the at least one rate equation.
The method is preferably implemented using a computer or computing device. Alternatively the method may be implemented using an electronic circuit. The dynamical physical system may be controlled by a control system, which preferably comprises a computer or computing device.
Alternatively the control system may comprise an electronic circuit. The control system preferably measures one or more quantities of the system using one or more input devices, sensors or the like. The values of quantities of the system may be calculated by the control system from the
measurements . The control system may calculate the
modification required to the controlled quantity of the system, and modify the quantity using one or more control devices . In accordance with a second aspect of the invention there is provided a dynamic physical system controlled according to any of the methods described above. The dynamical physical system may comprise a control system, which preferably comprises a computer or computing device. Alternatively the control system may comprise an electronic circuit. The control system may be arranged to measure one or more quantities of the system using one or more input devices, sensors or the like. The control system may be arranged to calculated values of quantities of the system from the measurements. The control system may be arranged to calculate the modification required to the controlled quantity of the system, and to modify the quantity using one or more control devices.
In the first aspect of the invention, the dynamic physical system may be a bioreactor. The bioreactor may be a bioethanol fermentor. Advantageously, the control term is D(CSf-C3), where Cs is the substrate concentration, CSf is a forcing term, and D is a constant. In such a case,
preferably the rate control function is:
>¾ where
Alternatively, the control term may be -DCP, where Cp is the product concentration and D is a constant. This allows the system to be controlled by controlling the concentration of the product, or in other words by removing the product from the system as it is produced. In such a case, preferably the rate control function is:
where
Ce
CI*
Cx is the biomass concentration and Ce is the active
component fraction.
In accordance with a third aspect of the invention there is provided a bioreactor controlled according to any of the methods described above.
In accordance with a fourth aspect of the invention there is provided a bioethanol fermentor controlled
according to any of the methods described above.
In the first aspect of the invention, the dynamic physical system may be a wind turbine. Preferably the wind turbine is controlled by varying the pitch of the blades of the turbine. Advantageously, the control term is Tr, the aerodynamic torque of the wind turbine. In that case, he rate control function may be:
where
The control term may be the wind speed. Alternatively, the control term may be the power output of the turbine generator. The wind turbine may be controlled by a control function that is a combination of a control functions derived from control terms for the wind speed and the power output of the turbine generator respectively.
The wind turbine may be controlled by controlling the torque of the power generator of the wind turbine.
Alternatively and/or additionally, the wind turbine may be controlled by controlling the angular velocity of the power generator of the wind turbine. The control term may be the angular velocity of the rotational input of the power generator. Alternatively and/or additionally, the control term may be the torque of the rotational input of the power generator .
In accordance with a fifth aspect of the invention there is provided a wind turbine controlled according to any of the methods described above.
In the first aspect of the invention, the dynamic physical system may be an HVAC system. The control term may be the volumetric flow-rate ^¾ of the exhaust air. In that case, the rate control function may be:
where : i v, +
Alternatively, the control term may be the volumetric flow-rate of the conditioned supply air. In that case, the rate control function may be: where :
Ή
In the first aspect of the invention, the dynamic physical system may be a semiconductor laser system. The semiconductor laser system may be of a non-linear laser ring cavity type. Advantageously, the control term may be the amount of dissipation within the cavity. Alternatively, the semiconductor laser system may be a laser operating with delayed optical feedback. Advantageously, the control term may be the feedback strength of the system. The
semiconductor laser system may be a physical optoelectronic device with a feedback loop. Advantageously, the control term may be the voltage at the output of the low-pass filter.
In accordance with a sixth aspect of the invention there is provided a laser system controlled according to any of the methods described above.
In accordance with a seventh aspect of the invention there is provided a non-linear laser ring cavity laser system controlled according to any of the methods described above .
In accordance with an eight aspect of the invention there is provided a laser operating with delayed optical feedback according to any of the methods described above.
In accordance with an ninth aspect of the invention there is provided a physical optoelectronic device laser system with a feedback loop according to any of the methods described above.
In the first aspect of the invention, the dynamic physical system may be an internal combustion engine.
The internal combustion engine may be controlled by controlling the pressure of the air intake of a cylinder of the internal combustion engine. In that case, the control term may be the drive-side torque acting against crankshaft rotation. Alternatively, the control term may be the crankshaft angular velocity of the internal combustion engine. Advantageously, both the torque of the cylinder against the rotation of the crankshaft and the crankshaft angular velocity are control terms.
Alternatively or additionally, the internal combustion engine may be controlled by controlling the pressure of the fuel injection line of a cylinder of the internal combustion engine. In that case, the control term may be the mass of oxygen in the cylinder.
Alternatively, the control term may be the pressure of the exhaust chamber of the internal combustion engine.
Alternatively, the control term may be the mass of fuel in a cylinder of the engine. Alternatively, the control term may be the ignition time of a cylinder of the engine.
In accordance with a tenth aspect of the invention there is provided an internal combustion engine controlled according to any of the methods described above. The internal combustion engine may be controlled by controlling the operation at least two cylinders of the internal combustion engine. Advantageously, the at least two cylinders are controlled using different control terms. Description of the Drawings
Embodiments of the present invention will now
described by way of example only, with reference to
accompanying schematic drawings of which:
Figure 1 is a graph showing the behaviour of a Rossler system;
Figure 2 is a flowchart describing the method of a first embodiment of the invention;
Figure 3a is a graph showing the change over time in the biomass concentration in a conventional bioethanol fermentor;
Figure 3b is a graph showing the change over time in the product concentration;
Figure 3c is a graph showing the change over time in the substrate concentration;
Figure 3d is a graph showing the change over time in the active component concentration;
Figure 4 is a graph showing the substrate concentration against the product concentration in the bioethanol fermentor;
Figure 5a is a graph showing the change over time in the biomass concentration in a bioethanol fermentor in which the control method of an embodiment of the present invention is applied at time 1000; Figure 5b is a graph showing the change over time in the product concentration; Figure 5c is a graph showing the change over time in the substrate concentration;
Figure 6 is a graph showing the change over time in the control function in the bioethanol fermentor; Figure 7 is a graph showing the substrate concentration against the product concentration in the bioethanol fermentor;
Figure 8 is a graph showing the results of the method on the product concentration as the system undergoes several parameter changes
Figure 9 is a graph showing the change in windspeed over time in a weather simulation;
Figure 10 is a graph showing the change over time in the angular velocity of a typical wind turbine and a wind turbine according to an embodiment of the invention;
Figure 11 is a graph showing the change over time in the rotor torque of the wind turbines;
Figure 12 is a graph showing the change over time in the shaft twisting of the wind turbines;
Figure 13 is a graph showing the change over time in the power output of the wind turbines;
Figure 14 is a graph showing the change in the control function in the wind turbine of the embodiment of the invention;
Figure 15 is a schematic diagram of a laboratory HVAC
system;
Figure 16 shows the behaviour of a typical non-linear
laser ring cavity laser system and a non-linear laser ring cavity laser system according to an embodiment of the invention;
Figure 17 shows the behaviour of the laser systems of
Figure 16 when the control is periodically turned on and off;
Figure 18 is a phase diagram showing the behaviour of a typical laser operating with delayed optical feedback and a laser operating with delayed optical feedback according to an embodiment of the invention;
Figure 19 shows the behaviour of the laser systems of
Figure 18 over time;
Figure 20 shows the behaviour of the laser system of the invention when the pumping current is periodically varied;
Figure 21 shows the output power over time of a typical physical optoelectronic device laser system with a feedback loop and a physical optoelectronic device laser system with a feedback loop according to an embodiment of the invention;
Figure 22 shows the output voltage of the laser systems; Figure 23 is a cylinder of an internal combustion engine; Figure 24 shows the change over time of angular velocity of the crankshaft of the internal combustion of the invention;
Figure 25 shows the change over time of the control
function based on the angular velocity of the crankshaft of the internal combustion of the invention; Figure 26 shows the change over time of the control function based on the torque of the cylinder against the rotation of the crankshaft of the internal combustion of the invention;
Figure 27 shows the change over time of the control
function based on the mass of oxygen in the cylinder of the internal combustion of the invention .
Detailed Description
A first embodiment of a method of the present invention is now described with reference to the flow chart of Figure 2. In the present embodiment, the method is being used to control a dynamic physical system comprising a set of variable quantities.
First, a model for the system is obtained (step 10) . The system is modelled by a set of variables that represent the variable quantities of the system, and the model comprises a set of rate equations for the variables that describe the rate of change of the quantities based on its current state. The model may be obtained, for example, by studying the system to derive its properties, in other words deriving the model from the observed properties of the system. Alternatively, the system may already have been studied and so a model is already available.
Next, a control term for at least one of the rate equations is identified (step 11) . As described above in relation to known methods of controlling chaotic behaviour, in chaotic states the global behaviour of a system tends to be dominated by non-linear parts of the system that allow quantities grow at an exponential rate. This implies that the growth of a variable in such a state will often be determined by a term including the variable itself, as this feedback of the value of the variable on its own growth rate causes the exponential growth. Thus, a control term for a rate equation for a particular variable will often be a term of the equation that includes the variable itself. However, it is not necessary (or even desirable) that the control term contains the particular variable; the control term can be any term within the rate equation that is liable to expand at an exponential rate.
A rate control function is then derived from the rate equation (step 12) . In particular, the rate control function is derived using the proportions of the variables that contribute to the growth rate of the rate equation. This is because, again as discussed above in relation to known methods of control of chaos, the local rate of expansion is proportional to the local behaviour of each of the variables.
The rate control function is then applied to the control term (step 13) . This provides a stabilised control term for the rate equation, in which the rate control function limits the control term in proportion to the divergence rate of the variable. Thus, if the model of the system is modified by using the stabilised control term in place of the original control term, this should give a model in which the system is stabilised. In some embodiments the rate control function may depend on one or more scalars. A scalar may apply to the rate control function as a whole, in which case it simply defines the strength of the control applied to the system. The scalar must not be too small, in which case the control is not sufficient to prevent exponential growth, or too large, in which case it would materially affect the overall behaviour of the system (for example by preventing chaotic behaviour altogether) . A suitable value for such a scalar may be determined experimentally, for example; any value that acts to keep the system stable is sufficient.
Alternatively or additionally, a scalar may apply only to the proportions of the variables in the rate control function, and not to the rate control function as a whole. Such a scalar will have a more precise effect on the behaviour of the rate control function beyond simply determining whether it successfully stabilises the system or not. For example, by varying the scalar it may be possible to select different orbits into which the system can be stabilised. Similarly to above, a suitable value for such a scalar may be determined experimentally. In practice, a scalar applied to the rate control function as a whole may need to be adjusted to take account of different values for a scalar applied to the proportions only.
Thus, a stabilised model of the system has been provided. In order to stabilise the system itself, the quantities whose variables are present in the control term are adjusted so that the control term derived from the modified quantities is the same (or at least substantially the same) as the stabilised control term (step 14) . In other words, the rate control function identifies what the result of calculating the control term using the values of the quantities should be in order to stabilise the system, and thus gives a result which the quantities should be varied so as to achieve.
To give a simple example, if a growth term consists of a single variable x, and at a given point in time the rate control function has a value of 0.5, this indicates that the quantity represented by the variable x should be reduced by one half in order to keep the system stable. If on the other hand the growth term comprises more than one variable, the quantities they represent can be varied in any suitable way such that the rate control function is satisfied. This means that the quantity that is most easily to adjust can be varied, for example.
Example 1
An embodiment in which the method of the invention is applied to an example system is now described. The system is the Rossler system, a well-known dynamic physical system that exhibits chaotic behaviour, which is defined by the following rate equations:
Considering third rate equation for the variable z, the growth of the variable z is given by the terms zx and β/α, and so the control term is taken to be the term zx (as the constant term β/α cannot be varied) . The proportion of each variable in the control term, to the growth rate is given by the quotients :
Figure imgf000021_0001
where ux and uz are constants.
These quotients are then used to derive a rate control function σ for the variable z, as follows:
Figure imgf000021_0002
where f and ξ are scalars as discussed above. (The scalar f is used to set the overall level of control applied to the system, while the scalar ξ is used to stabilise the system to different orbits.)
The rate control functions are then applied to the control term in rate equation, to give the following modified control term: ff( 'X.f z) ZX
which can be substituted in the rate equation as follows:
Figure imgf000021_0003
This gives a stabilised system, and as described above the rate control function together with the control term defines a target for varying the quantities represented by z and x in order to stabilise the system.
Example 2
An embodiment in which the method of the invention is applied to a bioreactor is now described; that is, a
biochemical process to manufacture a desired product by means of a biochemical reaction. Bioreactors and generic models thereof are described in Michael A. Henson.
Exploiting cellular biology to manufacture high-value products. IEEE Control Systems Magazine, pages 54-62, August 2006.
In particular, the method of the invention as applied to a bioethanol fermentor is now described. A bioethanol fermentor is a system that produces ethanol by fermenting biomass such as waste agricultural material (sugar, fruit, grains, potatoes etc.) . The biomass is fermented using a microorganism such as Zymomonas mobilis . A bioethanol fermentor is described in I. M. Jobses, G. T. Egberts, K. C. Luyben, and J. A. Roels, "Fermentation kinetics of Zymomonas mobilis at high ethanol concentrations: Oscillations in continuous cultures", Biotechnology and Bioengineering, 28(6) :868 - 877, 1986. Also described therein is a model of the bioethanol fermentor. Similarly, a bioethanol fermentor and its model are described in M. E. E. Abashar and S. S. E. H. Elnashaie, "Dynamic and chaotic behavior of periodically forced fermentors for bioethanol production. Chemical
Engineering Science", 65(16) :4894 - 4905, 15 August 2010. However, in this case the model of the bioethanol fermentor includes a sinusoidal forcing term, which describes a periodic adjustment to one of the quantities of the system, as described below.
A model of the system with forcing term is given by the following rate equations:
d .t> CJS c«< (k-s - ¾2 CF - k3 C - 0 C«
il t KS
Figure imgf000023_0001
for the active component concentration Ce, the product concentration Cp, the substrate concentration Cs, and the biomass concentration Cx of the system. μ is the specific growth rate of the system, given by the equation:
where max and Ks are constants. CSf is the forcing term, the periodic feed concentration given by the equation:
Figure imgf000023_0002
where Cso, A and ω are constants and t is time. A defines the amount of forcing applied, and ω the frequency of the forcing. (A model of the unforced system can be obtained by simply setting A to zero.) The other values in the
equations are constants.
An example of the chaotic behaviour of the system is shown by Figures 3a to 3d, which are graphs showing the changes in the biomass concentration, product concentration, substrate concentration and active component concentration respectively over time, and by Figure 4, which is a graph of the substrate concentration against the product
concentration .
Considering the rate equation for the substrate concentration Cs, the control term for this equation is the term D(CSf-C3) . The proportion of the variable Cs to the growth rate of the rate equation is given by the quotient:
C& + iii
From this equation the rate control function o(Cs) is derived :
This is applied to the control term to give a modified rate equation with stabilised control term as follows: i m ... L.
-v,.. ,/
This then shows how the quantities represented by the variables making up the control term can be adjusted in order to stabilise the system. In other words, in this particular case the rate control function defines how to adjust the substrate concentration Cs of the system in order to keep the system stabilised.
Considering instead the rate equation for the product concentration Cp, the control term for this equation is taken to be the term -DCP. This is because, in this case, it can be seen that production of the product Cp is inhibited by the presence of the product itself, and it is desired to control the system by removing the product as it is produced. Considering then the rate equation that contribute, the variables that contribute to the growth of the rate equation are the variables Cx and Ce, and their proportions to the growth rate of the rate equations are given by the equations :
i. x
€\·> 4- iia
From these equations the rate control function o(Cp) is derived :
This is applied to the control term to give a modified rate equation with stabilised control term as follows:
Figure imgf000025_0001
While in this embodiment control of both Cs and Cp has been described in order to control the system, in practice control of only one is required. However, control of the system could be achieved by controlling both variables (in practice controlling their underlying quantities, of course) , if desired.
The effect of the applying the control method is shown in Figures 5a to 5c, which are graphs showing the changes in the biomass concentration, product concentration and substrate concentration respectively over time, in which the control method is applied at time 1000. Similarly, Figure 6 is a graph showing the control function over time, as applied from time 1000. Figure 7 is a graph of the
substrate concentration against the product concentration showing the effect of applying the control method. Finally, Figure 8 is a graph showing the results of the method on the product as the system undergoes several parameter changes; at every lOOOh the forcing parameter A is increased. As can be seen, this makes the system go through several
oscillations, but nevertheless the method stabilises the system independent of the forcing term into stable
oscillations. (The dotted line indicates the yield of ethanol as mean extraction concentration at each forcing parameter change, and notably as can be seen from the right- hand side of the graph the yield can be increased with higher controlled chaotic parameters.)
Example 3
An embodiment in which the method of the invention is applied to a wind turbine power generator is now described. A wind turbine comprises a number of blades arranged around a rotor shaft, which convert kinetic energy in the wind into rotational movement. The rotor shaft is connected
mechanically to a generator, which converts the rotational movement into electricity. There exist many types of wind turbine design, each with their own particular
characteristics. In the following embodiment, the wind turbine is a variable speed, Horizontal Axis Wind Turbine (HAWT) . However, the invention is equally applicable to other designs of wind turbine.
Typical wind turbine control methods are primarily dependant on the wind speed, as measured by an anemometer located at or near the top of the turbine structure. The turbine control is typically split into three regions, primarily governed by the prevailing wind speed but also by the rotational frequency of the generator shaft and the mechanical and electrical limits of the equipment. These control regions are defined as follows:
Region 1: In region 1, the generator is decoupled from the rotor shaft (i.e. rotation of the rotor shaft is not transmitted to the generator) . If the wind speed is deemed too low to start the turbine the blades are pitched so as to produce minimal aerodynamic torque, so as to reduce stress on the blades while they are stationary. When the wind speed rises above a pre-determined value (based on the type of turbine) , the blades are pitched to the angle which provides maximal aerodynamic torque, with the result that the wind acts to rotate the turbine. Once the generator shaft has reached a pre-determined angular velocity (again based on the type of turbine), the generator control torque is enabled, coupling the rotor shaft to the generator. At this stage region 2 power production is entered.
Region 2: The wind speed in region 2 is less than that required to produce the rated power of the generator and the generator control torque is adjusted in an attempt to track optimal power capture for the current rotor angular velocity and wind speed. Region 3: The wind speed in region 3 is at least enough for the wind turbine to produce rated power output. The generator control torque is held at the rated torque. The control system is now responsible for maintaining the turbine at rated power by altering the rotor blade pitch to control the aerodynamic torque acting on the blades, and thereby controlling the amount of power extracted from the wind. The standard, non-proprietary, method of control used to alter the blade pitch and maintain rated power in this region is proportional-integral-derivative (PID) control, as described in Bossanyi, E. 1987, Adaptive Pitch Control for a 250 kW Wind Turbine, Proceedings of the British Wind Energy Conference, pp. 85-92; and Boukhezzar, B. and Siguerdid ane , H. Nonlinear Control of a Variable-Speed Wind Turbine Using a Two-Mass Model, IEEE Transactions on energy conversion, vol26, No.l. 2011.
The progression from region to region described above assumes a relatively smooth and steady change in wind speed. However, a particularly difficult scenario from the
perspective of a control mechanism is that of gusting wind conditions. During large and variable wind gusts the wind turbine control mechanism must be able to maintain the angular velocity of the rotor assembly below a rated value and minimise mechanical stresses. Under such conditions, and without an appropriate control mechanism as provided by the present invention, power production must be halted.
It is known to use a Proportional Integral Derivative (PID) controller to as a method of control in many
industries, including for the control of wind turbines.
However, PID control is linear in nature, whereas a wind turbine is a non-linear dynamical system, as wind speed is inherently chaotic. The performance shortcomings of PID control can be mitigated using feedback techniques such as neural networks and fuzzy control, but these add complexity to the tuning process and can require high precision and sampling rates. Further, PID control must be adjusted for different operating conditions in order to ensure that thresholds are not exceeded, and the inability of PID control to work throughout the full spectrum of operating conditions can itself be a reason for requiring the shutdown of a wind turbine in highly gusty wind conditions.
A model of the wind turbine is given by the following equations, taken from Eisenhut, C., Krug, F., Schram, C. and Klockl, B. 2007, Wind-Turbine Model for System Simulations Near Cut-in Wind Speed, IEEE Transactions on energy
conversion, vol. 22, No. 2; and Soltani, M., Wisniewski, R., Brath, .P and Boyd, S, Load Reduction of Wind Turbines Using Receding Horizon Control, Proceedings IEEE Multi-Conference on Systems and Control, Denver, September 2011:
lrti>v TT— Κβθ— Ββθ— Βτ ν
Figure imgf000029_0001
θ ^ \ -
Figure imgf000029_0002
where : 71
P
π
9 2Ν3λ
and where Cp is approximated by a non-linear function of λ and β, from Muhando, E.B., Senjyu, T . , Urasaki, N., Yona, A., Funabashi, T, Robust Predictive Control of Variable- Speed Wind Turbine Generator by Self-Tuning Regulator 2007 IEEE :
Figure imgf000030_0001
with :
and Ci = 0.5176, c = 116, C3 = 0.4, C4 = 5, C5 = 21 and 0.0068.
The following terms are used in the above equations and below :
T and u subscripts denote either rotor or generator
the blade pitch angle. A value of :::: 0 is used here.
the blade pitch angle: 90
Bo damping coefficient,, chosen to be 3<?5 Nrrt rad"
B,- and Bif Represents frictsve forces. 28 and 0.2 respectively
Power coefficie t
Optimal rotor power coefficient
and I« Inertia. 3es and 30 respectively kg
Stiffness constant of the drive-train assembly, 1 sr' Hm rad"
Shaft torsion rad
Tip-speed ratio and is given by t»T /V
Optimal tip-speed ratio
Gearing of the drivetrain. M=80
and Angular velocity of the rotor and generator respectively rats s "
Air density, taken here to be a constant 1.22521 kg n"!
Power output of the turbine generator W
Radius of the blade assembly, set h ere to .213 m
Generator torque Urn
Aerodynamic torque acting on turbine Hm Generator time constant
Biade pitch time constant
Wind speed as measured by the turbine anemometer m s'"
In order to control the wind turbine during gusty conditions, the aerodynamic torque Tr is adjusted in order to keep the amount of power extracted from the wind within the rated limit of the turbine. Considering the rate equation lr >r r the control term for this equation is simply Tr . From this, the following rate control function is derived :
where : and applying this to the existing model gives:
/v-:.'.v ···· <7: i : .. ϊ γ h'J* ···· ! ,;>■■■■ />,
In practice, the application of this control function to the aerodynamic torque would be used to determine how the blade pitch β should be adjusted.
The effect of applying the control method is shown with reference to Figures 9 to 14. Figure 9 shows a simulated wind speed data over a two day period. For the simulation, the constants of the control function were chosen as follows: f = 0.74, μΤΓ = 1.6e6 and ξ = -4.6. As can be seen, the wind speed on the second day was highly variable. For this example, the power which the turbine must not exceed is chosen to be ~400 KW, with a rated generator rotational frequency of 2350 rpm and therefore a rated generator torque of 1619 Nm. The wind speed at which this rated power occurs is around 10.31ms-1.
Figures 10 to 13 show the angular velocity of the generator, rotor torque, shear stress and power output of the turbine, without the control function being used shown by the variable dotted line, and with the control function being used shown by the variable solid line. (In practice if the turbine were to operate as shown by the dotted line it would suffer catastrophic failure, and so failsafe measures would act to prevent any movement of the blades as all, giving the flat dashed lines shown in the Figures.) Figure 14 shows the behaviour of the control function. As can be seen, over the first day during which the control function is not used (so maintains a value of 1), the systems behave identically. On day two, where the gusty conditions mean the control function is used, it can be seen that control function acts to keep the operation of the turbine within safe limits, while as shown in Figure 13 in particular allowing a significant amount of power to be produced by the turbine. In the simulation the turbine was able to produce ~4.88MWh over the second day, while without the control function the turbine would essentially be unusable on the second day so no power would be produced.
In a related embodiment, the control functions are derived by considering the wind turbine system to be a combination of two distinct systems working together. The first system is the turbine rotor assembly, which needs to be limited in cases of high and/or gusty wind. The second system is the generator system, which needs to be controlled to operate within its rated power limits.
Control functions as follows are derived:
Figure imgf000033_0001
where : v ···.·.·. V3
V'3 + μ.νϊ
These are applied to the model of the wind turbine as follows :
Figure imgf000034_0001
so the torque Tg and angular velocity cog of the generator are controlled by the control function derived from the angular velocity cog of the generator; and the blade pitch is controlled as follows:
, 1 .. „ where :
so the blade pitch β is controlled by the combination of the control functions based on the wind speed V and the power output of the turbine generator P.
Other controls that may be used are as follows:
<3 u.v f« ) :::: fwT e5^ ^ ;:
where :
Figure imgf000034_0002
Simulations of the model incorporating the above control functions, particularly in contrast to models incorporating PID control, provided a more controlled power output and so were able to operate at a higher mean power output while still ensuring that power spikes would not exceed the rated output of the generator. Further, the drive train of the turbine had a lower mean level of stress, and lower maximum stress peak amounts.
Example 4
An embodiment in which the method of the invention is applied to a Heating Ventilation and Air Conditioning (HVAC) system is now described. HVAC systems may be found
operating in a large number of different scenarios, such as offices, automobiles, factories, laboratories, refrigeration systems and life-support/protective suits.
A schematic drawing of a laboratory HVAC system is shown in Figure 15. The drawing shows a laboratory 100 and adjacent space 101; heat is conducted between the two as shown by the arrows 109. The laboratory 100 contains a fume hood 102. An actuator 109 draws air from the fume hood to a fume exhaust 105, causing air to be drawn from the
laboratory 100 into the fume hood 102. A further actuator 108 draws air from the laboratory 100 to a general exhaust 104. Finally, supply air 103 is blown into the laboratory 100 by an actuator 107, in the process being heated by a heating coil 106. A model for a laboratory HVAC system, derived from Osman, A., Mitchell, J.W. & Klein, S.A, Development of a Simulator for Laboratory HVAC System
(http : / /www . inive . org/members_area/medias /pdf/ Inive /clima200 0/1997/P274.pdf) is as follows:
Figure imgf000036_0001
where the terms used are as follows:
Figure imgf000036_0002
In some laboratory spaces the pressure of the air in the laboratory, P, will be maintained at a lower value than that of the adjacent spaces, Pa(j. This will help to prevent the flow of any harmful substances out laboratory.
Alternatively, in a space such as a clean-room, this relationship is reversed as in this case it is desirable to ensure that no contaminant can enter the laboratory from the adjacent spaces.
The control system for such an HVAC system would typically be responsible for maintaining specific variables such as temperature, pressure differential (Pa(j - P) , volumetric air-flow rates, general exhaust among others, within well-defined ranges corresponding to safety and/or comfort .
Considering the first equation, the control term is *'¾ (the volumetric flow-rate of the exhaust air) , giving a control function:
Figure imgf000037_0001
where :
Considering instead the second equation, the control term is (the volumetric flow-rate of the conditioned supply air), giving a control function:
(7 where :
(¾ + Ui:!)
The control functions can then be applied to the control terms in the first and second equation respectively, in order to determine how the exhaust air and/or conditioned supply air should be controlled.
Example 5
An embodiment in which the method of the invention is applied to control chaotic behaviour of lasers is now described. Semiconductor lasers are used in many
applications, for example in optical communication networks, DVD players and many other applications. The dynamic behaviour of lasers can be extremely complex, in particular under the influence of external perturbations such as feedback or delay coupling.
A first type of semiconductor laser system is a nonlinear laser ring cavity, which can be modelled using an Ikeda map as described in K. Ikeda, Multiple-valued
stationary state and its instability of the transmitted light by a ring cavity system, Optics Communications,
30 ( 2 ) : 257 -2 61 , 197 9 . By plotting the real and imaginary parts of the complex equation:
Figure imgf000038_0001
with x = Re(z) and y = Im(z), the Ikeda map can be described as follows, where typically a = 1 , b = 0.9, κ = 0 . 4 and η =
6 . 0 :
Figure imgf000039_0001
a -t- 1?{χ¾ cos φ --- \} n s φ
b('x¾sin«t> -ΐ-ιΐί,εοε φ)
Considering the equations xn+i and yn+i , the control term is φ, which represents the amount of dissipation within the cavity, giving a control function of :
x li i
i χ il l ^.... f t e ~ι.η * ίί .... ί -,Λ ' ν + χ→' ^ '
where
Figure imgf000039_0002
giving :
, ! ; .. . . ,.2
H« } ¾ 5
Taking μ = 5 and ξ = 1, Figure 16 shows the
uncontrolled Ikeda map, with on top two different stabilised orbits, a five period orbit indicated by open circles, and a fifteen period orbit indicated by stars. The orbit to which the system is stabilised will depend upon the proximity of the system to the orbit when control is enabled. By turning the control off and on again, different orbits may be found. To demonstrate this, Figure 17 shows a simulation in which control is enabled after 1000 timesteps and is turned off again after a further 1500 timesteps. This pattern is then repeated four times with the period of 2500 timesteps. It can be seen from Figure 17 that the system stabilises subsequently into a period 15 orbit and then a period 5 orbit, the period 15 orbit and the period 5 orbit again. A second type of semiconductor laser system is a laser operating with delayed optical feedback. A model of this type of system, the Lang-Kobayashi model, is described in T Sano, Antimode dynamics and chaotic itinerancy in the coherence collapse of semiconductor lasers with optical feedback, Physical Review A, 50 (3) : 2719-2726, 1994, and is as follows : φ(1
t) - 2y y P(t ·- t}?{t) cos (Δφ -f ? )
Figure imgf000040_0001
-ΔΝ - P{t K - Δ ) + ΔΙ where P is the photon number, φ the slowly varying part of the optical phase, and the deviation of the carrier number from the threshold value Nth of the solitary laser ΔΝ = N - Nth. ω is the optical angular frequency, which is assumed to be zero. κ = 1000 is the damping constant for the photon number, and AJ the pumping current deviation from the threshold. The feedback strength γ is given by:
Figure imgf000040_0002
where R is the power reflectivity of the external mirror, r that of the laser facets, Tc the round trip time of the cavity and η the coupling ratio. Other parameter values are a = 6, β = 10"5, n = 1, Nth = 10
Considering the equations for the changes in P or φ, the control term is taken to be γ; this is also a practical choice as it can be controlled in a physical implementation, The control function is σ then given by:
Φ
Φ + μψ
o ?, φ :— f e" ^ q*
where f = 6 and ξ = -1, giving the new term:
ri
~V
Figure 18 shows the phase space representation of the carrier number of the system versus the phase difference, with the uncontrolled system shown underneath and the controlled system shown in black on top. Figure 19
similarly shows the change in carrier number over time, where control is enabled at time = 200au. In both it can be seen that the control function substantially stabilises the behaviour of the system. Further, figure 20 shows the behaviour of the behaviour of the system over time where every 500 time steps the pumping current AJ is increased (from -10 up to +1) . It can be seen that the system goes through short transients but stabilises quickly every time. The horizontal line (e.g. from 500 to 1000) means that the system is in controlled steady state. The thick black bars mean that the system is oscillating at a controlled stable but high frequency. Because the control method of the invention is effective in controlling the chaotic system independent of parameter changes which may cause
bifurcations, it can be used to maintain the stability of the laser with different amount of the pumping current AJ. A third type of semiconductor laser system is a
physical optoelectronic device with a feedback loop. A model of this type of system, the Blakely model, is
described in JN Blakely, L Illing, and Daniel J Gauthier,
High-speed chaos in an optical feedback system with flexible timescales, Quantum Electronics, IEEE, 40 ( 3 ) : 299-305 , 2004; and in Y. G. Zheng and Z. H. Wang. Stability and Hopf bifurcations of an optoelectronic time-delay feedback system, Nonlinear Dynamics, 57(1-2) : 125-134, September 2008. Its main features are delay feedback with additional low- pass and high-pass filters to give it a band-pass
characteristic. The Blakely model is described by the equations :
Figure imgf000042_0001
where v(t) is the voltage at the output of the low-pass filter, p(t) is the laser output which relates to the voltage at the output of the high-pass filter, Vdet(t) is the voltage output of the (non-linear) photodiode, po = 26 is the emission power, T0 = 19.1 is the timed delay in the feedback loop, ii = 0.66 is the low-pass filter time
constant, ih = 22 is the high-pass filter time constant, γ = 0.0053 is the system amplification by feedback strength, k = 4.8 is the voltage to power conversion strength, a = 1.89 determines the sensitivity of the interferometer and β = 0.8 is the fringe visibility. Considering the equation for the change in time of p, the control term of this equation is taken to be v(t), voltage at the output of the low-pass filter. This is also practical for physical implementation, as feedback strength is not easily modifiable. The control function and modified rate equation are then:
V
< > ··—:
v
-—
sip. ) ^ fεΕ *v
pit) ' -{pit}— ) ., , , ,
____ ,,,, : _ . j v| j o†p, j where f = 20 and ξ = -1.
Figures 21 and 22 show the change of the output power and voltage over time, for the uncontrolled system in grey and the controlled system in black. The control is
implemented at time = 200au. The system is immediately controlled into a stable 4-orbit that greatly reduces the needed output power to maintain the stability of the system, thereby making it more efficient and economical, the reduction in needed output being approximately 56%.
Example 6
An embodiment in which the method of the invention is applied to an internal combustion engine is now described. An internal combustion engine comprises a crankshaft which is driven by a plurality of cylinders . An exemplary cylinder is shown in Figure 23. The cylinder 200 comprises a cylinder cavity 201 within which a piston 202 can move. The piston 202 drives a piston rod 203 which in turn rotates the crankshaft 204. The cylinder 200 has a fuel intake 205 though which a fuel injector injects fuel (for example octane) into the cylinder cavity 201, and an air intake 206 through which air is taken into the cylinder cavity 201. As will be well understood by the skilled person, the cylinder 200 has "four-stroke" operation comprising an intake stroke, a compression stroke, a power stroke and an exhaust stroke.
In conventional engines various parameters such as amount of fuel to be injected are controlled by means of lookup tables, which use properties of the engine such as current speed, torque, temperature, acceleration and the like in order to provide pre-determined ideal values for the parameters. The ideal values are typically determined by testing an engine throughout its range of performance at the product development stage. However, while this allows a certain level of performance to be guaranteed, it does not allow the optimal performance of which the engine may be capable. In the internal combustion engine of the present embodiment of the invention, the various parameters of the engine are controlled using control functions based on the properties of the engine, as described below.
The operation of the engine is modelled as follows. The height h of the cylinder head with respect to the crankshaft angle A is described by:
*¾ _
li reosf A.) - (l* -- r^si^iA }2
where r is the radius of the crankshaft 204 and 1 is the length of the arm 203. The mass-flow equations for the fuel and other components are as follows:
r
Figure imgf000045_0001
where typical modelling assumptions have been made regarding heat loss (zero heat loss during the combustion process) and ideal gas behaviour. Dissociation of other combustion species is not specifically modelled, rather it has been assumed that excess fuel is not combusted fully or partially and is represented as octane in the combustion products.
The crank dynamics are then described by:
Iclb — Tc Κβθ - Ββθ - α>, αν
0
N
while the torque on the crank Tc due to the force acting on the piston head as a result of the pressure inside the cylinder cavity is described by: rsiriiA)
φ— arcs in (
I
Figure imgf000046_0001
F1 = Fcos(d)]
Figure imgf000046_0002
Tc = rFc
where A is the crank-angle and is in the range [0, 4n) radians, F is the force exerted on the piston-head by the gas mixture, Fi is the force transmitted down the piston rod, Fc is the tangential force exerted on the crankshaft at the interface between the crankshaft and the piston rod and Tc is the torque exerted on the crankshaft by the piston.
Combustion is modelled using the octane combustion equations described in Flagan, Richard C. & Seinfeld, John H. Fundamentals of air pollution engineering. Prentice-Hall, Inc., Englewood Cliffs, New Jersey (1988) . In a
stoichiometric situation (there is exactly the amount of air required to burn all the fuel) or off-stoichiometric situation (more air than required to burn all the fuel) , for one mole of fuel the reaction equation is:
s 4- x(C¾ + 3.76N2)→ 8CO24- H204- (3.76χ)Ν24- ίχ 12.5j£>2
In a rich air/fuel mixture (not enough air required to burn all the fuel) again for one mole of fuel the reaction equati is:
C*HiR +
Figure imgf000046_0003
where χ = 12.5 for complete stoichiometric combustion. It is assumed that there is no partial burning of octane during rich conditions, and that octane that is not combusted is simply present in the reactants as pure octane. It is also assumed that air composition is 80% Nitrogen and 20% Oxygen.
The heat produced during this reaction is derived from a heat balance equation of this reaction equation, assuming constant-volume specific heat capacity values of the
constituent species. The temperature change due to
combustion is approximated by calculating the change in temperature using an averaging of specific heat values for all substances involved (taken at a reasonable temperature of approx. 500K) . This model is representative enough of the dynamics involved in the combustion process and the assumptions do not adversely affect the results the fine details of the combustion process do not significantly affect the general behaviour of the system. For example, the modelling of dissociation of species such as carbon monoxide, hydrogen and nitrogen oxide would not
significantly alter the macroscopic dynamics of the system.
Using this model, the following control functions are derived based on the drive-side torque acting against crankshaft rotation, the crankshaft angular velocity and the mass of oxygen in the cylinder (which could in practice be estimated using mass flow sensors) :
Figure imgf000047_0001
where : r μ , .
The control functions are then used to control the pressure of the air intake and the pressure of the fuel injection line as follows:
Figure imgf000048_0001
( L - MI;
!,J ,J
In other words, the control function based on the mass of oxygen in the cylinder is used to control the pressure of the fuel injection line, while the control functions based on the drive-side torque acting against crankshaft rotation and the crankshaft angular velocity are used to control the pressure of the air intake.
Other control functions that could be used are as follows :
where
m <·
Figure imgf000048_0002
The terms used in the above are as follows
A Crank-angle. [0, An) radians (four-stroke mechanism) rad
Relative atomic trass of carbon: 12,0107
Relative atomic trass of hydrogen: 1.00784
Relative atomic trass of nitrogen: 14.0067
Relative atomic trass of oxygen: 15.9994
Area of the air-infst orifice: *t0.007&i rrr Area of the exhaust orifice: :?i0.0075 m2
«r Area of the fuel-Injector orifice: O.0012 m2 c a d subscripts denote either crankshaft or drive-side
"Average" specific heat (assuming specific heat values for substances at ··· 500K) J mol- ' K" 1 Damping coefficient Nm rad" 1 sec "
Orifice coefficient of ffow for octane, set to 1 for simplification
Q„ Orifice coefficient of ffow for oxygen, set to 1 for simplification
C!;j Orifice coefficient of ffow for nitrogen, set to 1 for simplification
Orifice coefficient of ffow for carbon dioxide, set to 1 for simplification
Orifice coefficient of ffow for water, set to 1 for simplification
<pc and φ, Coefficients of friction. 0.05 and 0 1 respectively
an Enthalpy of octane combustion (assume typical enthalpy of formation values) J root" '
Appro*, temperature change due to combustion (assume c?*9): dTemp dH/cf* K
F Force exerted on cylinder head due to the pressure of gasses In the cylinder cavity N
Force transmitted down crank rod. Accounts for cylinder/rod bearing interface N Fc Force exerted on crankshaft. Translation of Ft to tangent of crankshaft surface N
I..- and ld Inertia of the crank and drive-shaft respectively kg m2
Stiffness constant of the drive-train assembly Nm rad"
Length of the cylinder connecting rod: 0.25 m
Mass of fuei {octane} kg
Mass of oxygen kg
Mass of nitrogen kg
Mass of water kg
Mass of carbon dioxide kg
Number of moles of octane in reaction products mol umber of moles of octane in reaction reactants mol Number of moles of oxygen in reaction products mol Number of moles of oxygen in reaction reactants mol Number of moles of nitrogen in reaction products mol
Number of moles of nitrogen in reaction products
Number of moles of nitrogen in reaction reactants
Number of moles of water in reaction products
Number of moles of water in reaction reactants
Number of moles of carbon dioxide in reaction product
Number of moles of carbon dioxide in reaction reactants
Figure imgf000049_0001
Gearing of the drivetrain: Scalar value
Engine power output
Pressure of the air intake
Pressure in cylinder cavity ( During combustion, APCisi calculated using dTemp)
Pressur of the exhaust chamber
Pressure of the fuel injection line
VP Partial pressure of octane
pp Partial pressure of oxygen
VP Partial pressure of nitrogen
pp Partial pressure of carbon dioxide
PP Partial pressure of gaseous water cut and: α*ά Anguiar velocity of the crank and drive-side respectively rad s
Pf Fuel (octane) density, taken here to be a constant 745 kg m
PuAtmospheric oxygen density, taken here to be a constant 1.43 kg ητ pil; Atmospheric nitrogen density, taken here to be a constant 1 ,251 kg ητ
Pi «2 Carbon diox de density kg ητ
Water density kg ητ
Radius of cylinder head: 0.05 m
T radius of the crank arm: 0.05 m
¾fA Boolean function, 1 if cylinder is in the air-intake phase, zero otherwise
¾' Boolean function, 1 if cylinder is in the combustion phase, zero otherwise
% Boolean function, 1 if cylinder is in the exhaust phase, zero otherwise
Wv Boolean function, 1 if cylinder is m the injection phase, zero otherwise
ip Engine power time constant s
'7,· Crank torque Nm fa Drive-side (Drive-shaft) torque acting against crankshaft rotation Nm
Θ Drive-train torsion rad
V Volume of cylinder cavity (varies with cylinder position) m3
The crankshaft angular velocity of the engine
controlled in the embodiment is shown in Figure 24, where the external torque demand on the system is increased substantially at times 5xl05 and 10x10s. Instead of
stalling, it can be seen that the control functions adjust the fuel and air injection pressures to maintain a stable system. The variations of the control functions over time are shown in Figures 25 to 27.
The skilled person will appreciate that in other embodiments each cylinder of the engine could be controlled using a different set of control functions to control different parameters of the cylinder.
Whilst the present invention has been described and illustrated with reference to particular embodiments, it will be appreciated by those of ordinary skill in the art that the invention lends itself to many different variations not specifically illustrated herein. In particular, it will be appreciated that the invention will have application to chaotic system found in many different areas of technology. By way of example only, certain possible variations will now be described.

Claims

Claims
1. A method of controlling a dynamic physical system, wherein the system comprises a plurality of variable quantities, the method comprising the steps of:
obtaining a model of the system, the model comprising a plurality of variables representing the variable quantities, and a plurality of respective rate equations that describe the rate of change of the variables;
identifying a control term in at least one rate equation from the plurality of rate equations;
deriving a rate control function from, for at least one of the variables in the rate equation, the proportion of the variable to the growth rate of the rate equation;
applying the rate control function to the control term to provide a stabilised control term; and
controlling the dynamic physical system by modifying at least one of the quantities represented by the variables in the control term, so that the control term derived from the modified quantities is substantially the same as the stabilised control term.
2. A method as claimed in claim 1, wherein the control term includes the variable whose behaviour is described by the rate equation.
3. A method as claimed in claim 1 or claim 2, wherein the proportion qx of a variable x to the growth rate is given by the equation: X
where μχ is a constant.
4. A method as claimed in any preceding claim, wherein the rate control function is of the form:
where qxi to qxn are the proportions of the variables xl to xn to the growth rate, and f and ξ are scalars.
5. A method as claimed in claim 4, where f and ξ are varied so as to stabilise the system into a pre-determined orbit.
6. A method as claimed in any preceding claim, wherein the at least one rate equation describes exponential growth of its respective variable.
7. A method as claimed in any preceding claim, wherein the control term contributes to the growth of the respective variable of the at least one rate equation.
8. A dynamic physical system controlled according to any of the methods of claims 1 to 7.
9. A method as claimed in any of claims 1 to 7, wherein the dynamic physical system is bioreactor.
10. A method as claimed in claim 9, wherein the bioreactor is a bioethanol fermentor.
11. A method as claimed in claim 9, wherein the control term is D(CS f-C3), where Cs is the substrate concentration, CS f is a forcing term, and D is a constant.
12. A method as claimed in claim 11, wherein the rate control function is:
where :
13. A method as claimed in claim 9, wherein the control is -DCP, where Cp is the product concentration and D is constant .
14. A method as claimed in claim 11, wherein the rate control function is:
where :
Figure imgf000054_0001
Cx is the biomass concentration and Ce is the active
component fraction.
15. A bioreactor controlled according to the method of claim 9.
16. A bioethanol fermentor controlled according to any of the methods of claims 10 to 14.
17. A method as claimed in any of claims 1 to 7, wherein the dynamic physical system is a wind turbine.
18. A method as claimed in claim 17, wherein the wind turbine is controlled by varying the pitch of the blades of the turbine.
19. A method as claimed in claim 17 or 18, wherein the control term is Tr, the aerodynamic torque of the wind turbine .
20. A method as claimed in claim 19, wherein the rate control function is:
where :
21. A method as claimed in claim 17 or 18, wherein the control term is the wind speed.
22. A method as claimed in claim 17 or 18, wherein the control term is the power output of the turbine generator.
23. A method as claimed in any of claims 17 to 22, wherein the wind turbine is controlled by controlling the torque of the power generator of the wind turbine.
24. A method as claimed in claim 17 or 23, wherein the wind turbine is controlled by controlling the angular velocity of the power generator of the wind turbine.
25. A method as claimed in claim 17, 23 or 24, wherein the control term is the angular velocity of the rotational input of the power generator.
26. A method as claimed in claim 17, 23 or 24, wherein the control term is the torque of the rotational input of the power generator.
27. A wind turbine controlled according to any of the methods of claim 17 to 26.
28. A method as claimed in any of claims 1 to 7, wherein the dynamic physical system is an HVAC system.
29. A method as claimed in claim 28, wherein the control term is the volumetric flow-rate ^f; of the exhaust air.
30. A method as claimed in claim 29, wherein the rate control function is: Jr where
31. A method as claimed in claim 28, wherein the control term is the volumetric flow-rate of the conditioned supply air.
32. A method as claimed in claim 31, wherein the rate control function is:
where : ifc, = —
{ V}: - - ½¥ )
33. A method as claimed in any of claims 1 to 7, wherein the dynamic physical system is a semiconductor laser system.
34. A method as claimed in claim 33, wherein the
semiconductor laser system is of a non-linear laser ring cavity type.
35. A method as claimed in claim 34, wherein the control term is the amount of dissipation within the cavity.
36. A method as claimed in claim 27, wherein the semiconductor laser system is a laser operating with delayed optical feedback.
37. A method as claimed in claim 36, wherein the control term is the feedback strength of the system.
38. A method as claimed in claim 33, wherein the
semiconductor laser system is a physical optoelectronic device with a feedback loop.
39. A method as claimed in claim 38, wherein the control term is the voltage at the output of the low-pass filter.
40. A semiconductor laser system controlled according to the method of claim 33.
41. A non-linear laser ring cavity laser system controlled according to the methods of any of claims 34 to 35.
42. A laser operating with delayed optical feedback
controlled according to the methods of any of claims 36 to 37.
43. A physical optoelectronic device laser system with a feedback loop controlled according to the methods of any of claims 38 to 39.
44. A method as claimed in any of claims 1 to 7, wherein the dynamic physical system is an internal combustion engine.
45. A method as claimed in claim 44, wherein the internal combustion engine is controlled by controlling the pressure of the air intake of a cylinder of the internal combustion engine.
46. A method as claimed in claim 45, wherein the control term is the drive-side torque acting against crankshaft rotation of the internal combustion engine.
47. A method as claimed in claim 44 or 45, wherein the control term is the crankshaft angular velocity of the internal combustion engine.
48. A method as claimed in any of claims 44 to 47, wherein the internal combustion engine is controlled by controlling the pressure of the fuel injection line of a cylinder of the internal combustion engine.
49. A method as claimed in claim 48, wherein the control term is the mass of oxygen in the cylinder.
50. An internal combustion engine controlled according to the method of any of claims 44 to 49.
51. An internal combustion engine as claimed in claim 50, wherein the internal combustion engine is controlled by controlling the operation at least two cylinders of the internal combustion engine.
52. An internal combustion engine as claimed in claim 51, wherein the at least two cylinders are controlled using different control terms.
PCT/GB2012/052739 2011-11-03 2012-11-05 A method of controlling a dynamic physical system that exhibits a chaotic behaviour WO2013064840A1 (en)

Priority Applications (2)

Application Number Priority Date Filing Date Title
GB1119036.0 2011-11-03
GBGB1119036.0A GB201119036D0 (en) 2011-11-03 2011-11-03 A method of controlling a dynamic physical system

Applications Claiming Priority (5)

Application Number Priority Date Filing Date Title
US14/356,110 US9740180B2 (en) 2011-11-03 2012-11-05 Method of controlling a dynamic physical system that exhibits a chaotic behaviour
KR1020147015156A KR20140100946A (en) 2011-11-03 2012-11-05 A method of controlling a dynamic physical system that exhibits a chaotic behaviour
CN201280065939.5A CN104303115A (en) 2011-11-03 2012-11-05 A method of controlling a dynamic physical system that exhibits a chaotic behaviour
EP12788245.4A EP2774006A1 (en) 2011-11-03 2012-11-05 A method of controlling a dynamic physical system that exhibits a chaotic behaviour
JP2014539408A JP6161618B2 (en) 2011-11-03 2012-11-05 Control method of dynamic physical system showing chaotic behavior

Publications (1)

Publication Number Publication Date
WO2013064840A1 true WO2013064840A1 (en) 2013-05-10

Family

ID=45375783

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/GB2012/052739 WO2013064840A1 (en) 2011-11-03 2012-11-05 A method of controlling a dynamic physical system that exhibits a chaotic behaviour

Country Status (7)

Country Link
US (1) US9740180B2 (en)
EP (1) EP2774006A1 (en)
JP (1) JP6161618B2 (en)
KR (1) KR20140100946A (en)
CN (1) CN104303115A (en)
GB (1) GB201119036D0 (en)
WO (1) WO2013064840A1 (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112363394A (en) * 2020-11-02 2021-02-12 中国科学院力学研究所 Ventilation supercavitation form control method based on machine learning

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6331750B2 (en) * 2014-06-23 2018-05-30 三菱自動車工業株式会社 Engine control device
US10344615B2 (en) * 2017-06-22 2019-07-09 General Electric Company Method and system for schedule predictive lead compensation
EP3492772B1 (en) * 2017-11-29 2022-02-09 Dr. Ing. h.c. F. Porsche AG Method for vibration attenuation in an automotive drivetrain, corresponding data processing apparatus, computer program, and data carrier signal
JP6765394B2 (en) 2018-05-22 2020-10-07 京セラ株式会社 Measuring device and control method of measuring device
CN108873690B (en) * 2018-08-29 2021-06-01 东北大学 Trajectory tracking method of second-order strict feedback chaotic system
CN108845494B (en) * 2018-08-29 2021-06-01 东北大学 Second-order strict feedback chaotic projection synchronization method
CN108931917B (en) * 2018-09-04 2021-06-01 东北大学 Three-order strict feedback chaotic projection synchronization method
CN109062054B (en) * 2018-09-04 2021-06-01 东北大学 Three-order strict feedback chaotic track tracking method
CN113109391B (en) * 2021-05-21 2022-03-04 中国计量大学 Adiabatic reaction dynamic specific heat capacity measuring method based on compensation power modulation

Family Cites Families (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB8329218D0 (en) * 1983-11-02 1983-12-07 Ffowcs Williams J E Reheat combustion system for gas turbine engine
CA2132194A1 (en) * 1993-01-19 1994-08-04 Robert R. Klevecz Method for control of chaotic systems
JP4490000B2 (en) * 2001-06-19 2010-06-23 本田技研工業株式会社 Air-fuel ratio control device for internal combustion engine
DE60236394D1 (en) * 2001-09-20 2010-07-01 Honda Motor Co Ltd Control device for multipurpose engine
JP3880861B2 (en) * 2002-01-22 2007-02-14 本田技研工業株式会社 Air-fuel ratio control device for internal combustion engine
US20060165138A1 (en) * 2005-01-21 2006-07-27 Alexander Kachanov Frequency doubling of semiconductor lasers to generate 300-600 nm light
JP2008529062A (en) * 2005-01-21 2008-07-31 プレジデント アンド フェローズ オブ ハーバード カレッジ オフィス オブ テクノロジー アンド トレードマーク ライセンシング System and method for providing a tunable optical parametric oscillator laser system providing dual frequency output for nonlinear vibrational spectroscopy and microscopy
JP4337868B2 (en) * 2006-12-12 2009-09-30 トヨタ自動車株式会社 Control device for gas turbine engine
US20080277486A1 (en) * 2007-05-09 2008-11-13 Johnson Controls Technology Company HVAC control system and method
ES2382631B1 (en) * 2009-09-03 2013-05-03 Gamesa Innovation & Technology, S.L. METHODS AND SYSTEMS OF AIRCRAFT CONTROL
BR112013022994A2 (en) * 2011-03-14 2016-12-06 Dow Brasil Sudeste Ind Ltda method for controlling foaming in an aqueous alcoholic fermentation broth
US8793004B2 (en) * 2011-06-15 2014-07-29 Caterpillar Inc. Virtual sensor system and method for generating output parameters

Non-Patent Citations (19)

* Cited by examiner, † Cited by third party
Title
ATAEI M ET AL: "Control of chaos in permanent magnet synchronous motor by using optimal Lyapunov exponents placement", PHYSICS LETTERS A, NORTH-HOLLAND PUBLISHING CO., AMSTERDAM, NL, vol. 374, no. 41, 13 September 2010 (2010-09-13), pages 4226 - 4230, XP027260826, ISSN: 0375-9601, [retrieved on 20100820] *
BOSSANYI, E.: "Adaptive Pitch Control for a 250 kW Wind Turbine", PROCEEDINGS OF THE BRITISH WIND ENERGY CONFERENCE, 1987, pages 85 - 92
BOUKHEZZAR, B.; SIGUERDIDJANE, H.: "Nonlinear Control of a Variable-Speed Wind Turbine Using a Two-Mass Model", IEEE TRANSACTIONS ON ENERGY CONVERSION, vol. 26, no. 1, 2011
DONGMEI WU ET AL: "Non-feedback Control on Chaotic Logistic System", INFORMATION PROCESSING, 2009. APCIP 2009. ASIA-PACIFIC CONFERENCE ON, IEEE, PISCATAWAY, NJ, USA, 18 July 2009 (2009-07-18), pages 215 - 217, XP031505134, ISBN: 978-0-7695-3699-6 *
E. OTT; C. GREBOGI; J. A. YORKE: "Controlling chaos", PHYS. REV. LETT., vol. 64, 1990, pages 1196, XP002926846, DOI: doi:10.1103/PhysRevLett.64.1196
EISENHUT, C.; KRUG, F.; SCHRAM, C.; KLOCKL, B.: "Wind-Turbine Model for System Simulations Near Cut-In Wind Speed", IEEE TRANSACTIONS ON ENERGY CONVERSION, vol. 22, no. 2, 2007, XP011184737, DOI: doi:10.1109/TEC.2006.875473
FLAGAN; RICHARD C.; SEINFELD; JOHN H.: "Fundamentals of air pollution engineering", 1988, PRENTICE-HALL, INC.
FRADKOV A L ET AL: "Control of chaos: Methods and applications in engineering", ANNUAL REVIEWS IN CONTROL, PERGAMON, AMSTERDAM, NL, vol. 29, no. 1, 1 January 2005 (2005-01-01), pages 33 - 56, XP027657942, ISSN: 1367-5788, [retrieved on 20050101] *
I. M. JOBSES; G. T. EGBERTS; K. C. LUYBEN; J. A. ROELS: "Fermentation kinetics of Zymomonas mobilis at high ethanol concentrations: Oscillations in continuous cultures", BIOTECHNOLOGY AND BIOENGINEERING, vol. 28, no. 6, 1986, pages 868 - 877
JN BLAKELY; L ILLING; DANIEL J GAUTHIER: "High-speed chaos in an optical feedback system with flexible timescales, Quantum Electronics", IEEE, vol. 40, no. 3, 2004, pages 299 - 305
JUNJI OHTSUBO: "Feedback Induced Instability and Chaos in Semiconductor Lasers and Their Applications", OPTICAL REVIEW, SPRINGER, BERLIN, DE, vol. 6, no. 1, 1 January 1999 (1999-01-01), pages 1 - 15, XP019353663, ISSN: 1349-9432, DOI: 10.1007/S10043-999-0001-Z *
K. IKEDA: "Multiple-valued stationary state and its instability of the transmitted light by a ring cavity system", OPTICS COMMUNICATIONS, vol. 30, no. 2, 1979, pages 257 - 261, XP024490117, DOI: doi:10.1016/0030-4018(79)90090-7
M. E. E. ABASHAR; S. S. E. H. ELNASHAIE: "Dynamic and chaotic behavior of periodically forced fermentors for bioethanol production", CHEMICAL ENGINEERING SCIENCE, vol. 65, no. 16, 15 August 2010 (2010-08-15), pages 4894 - 4905, XP027130791
MICHAEL A. HENSON: "Exploiting cellular biology to manufacture high-value products", IEEE CONTROL SYSTEMS MAGAZINE, August 2006 (2006-08-01), pages 54 - 62
MUHANDO, E.B.; SENJYU, T.; URASAKI, N.; YONA, A.; FUNABASHI, T: "Robust Predictive Control of Variable-Speed Wind Turbine Generator by Self-Tuning Regulator", IEEE, 2007
PHYSICAL REVIEW A, vol. 50, no. 3, 1994, pages 2719 - 2726
SOLTANI, M.; WISNIEWSKI, R.; BRATH, .P; BOYD, S: "Load Reduction of Wind Turbines Using Receding Horizon Control", PROCEEDINGS IEEE MULTI-CONFERENCE ON SYSTEMS AND CONTROL, DENVER, September 2011 (2011-09-01)
Y. G. ZHENG; Z. H. WANG: "Stability and Hopf bifurcations of an optoelectronic time-delay feedback system", NONLINEAR DYNAMICS, vol. 57, no. 1-2, September 2008 (2008-09-01), pages 125 - 134, XP019685820
YAZDANPANAH AM ET AL: "Adaptive control of chaos in nonlinear chaotic discrete -time systems", PHYSICS AND CONTROL, 2005. PROCEEDINGS. 2005 INTERNATIONAL CONFERENCE SAINT PETERSBURG, RUSSIA AUG. 24-26, 2005, PISCATAWAY, NJ, USA,IEEE, 24 August 2005 (2005-08-24), pages 913 - 915, XP010839637, ISBN: 978-0-7803-9235-9 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112363394A (en) * 2020-11-02 2021-02-12 中国科学院力学研究所 Ventilation supercavitation form control method based on machine learning
CN112363394B (en) * 2020-11-02 2022-04-15 中国科学院力学研究所 Ventilation supercavitation form control method based on machine learning

Also Published As

Publication number Publication date
EP2774006A1 (en) 2014-09-10
GB201119036D0 (en) 2011-12-14
CN104303115A (en) 2015-01-21
US9740180B2 (en) 2017-08-22
JP6161618B2 (en) 2017-07-12
US20140309755A1 (en) 2014-10-16
KR20140100946A (en) 2014-08-18
JP2014532931A (en) 2014-12-08

Similar Documents

Publication Publication Date Title
US9740180B2 (en) Method of controlling a dynamic physical system that exhibits a chaotic behaviour
JP6248124B2 (en) Renewable energy compatible gas turbine and control method thereof
Lescher et al. Robust gain scheduling controller for pitch regulated variable speed wind turbine
WO2016038672A1 (en) Power generation system and power generation method
Corradini et al. An observer-based blade-pitch controller of wind turbines in high wind speeds
CN109296500B (en) Maximum wind energy capture method based on robust control theory
Prats et al. Improving transition between power optimization and power limitation of variable speed, variable pitch wind turbines using fuzzy control techniques
Burkart et al. Nonlinear control of wind turbines: An approach based on switched linear systems and feedback linearization
Dang et al. Nonlinear model predictive control (NMPC) of fixed pitch variable speed wind turbine
Ranjan et al. Parabolic trough solar–thermal–wind–diesel isolated hybrid power system: active power/frequency control analysis
Bratcu et al. Optimal control of wind energy conversion systems: From energy optimization to multi-purpose criteria-A short survey
Dang et al. Model predictive control for maximum power capture of variable speed wind turbines
Rubio et al. Control of two electrical plants
Bagherieh et al. Utilization of blade pitch control in low wind speed for floating offshore wind turbines
Gosk Model predictive control of a wind turbine
Dang et al. A multi-objective optimal nonlinear control of variable speed wind turbine
CN106712615B (en) The control method of Gas Turbine Generating Units in energy mix electric system
Gionfra et al. Hierarchical control of a wind farm for wake interaction minimization
Sirigu et al. Wind generator innovative blade design: variable twist and start-up control
Pachauri et al. Pitch angle controlling of wind turbine system using proportional-integral/fuzzy logic controller
Sahin et al. Performance comparison of two turbine blade pitch controller design methods based on equilibrium and frozen wake assumptions
Freeman et al. Direct model-reference adaptive control of variable speed horizontal-axis wind turbines
Nejadkhaki et al. Integrative control and design framework for an actively variable twist wind turbine blade to increase efficiency
Ofualagba et al. The modeling and dynamic characteristics of a variable speed wind turbine
Cavanini et al. A Control Strategy for Variable-Speed Variable-Pitch Wind Turbines within the regions of partial-and full-load operation without wind speed feedback

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 12788245

Country of ref document: EP

Kind code of ref document: A1

ENP Entry into the national phase

Ref document number: 2014539408

Country of ref document: JP

Kind code of ref document: A

WWE Wipo information: entry into national phase

Ref document number: 14356110

Country of ref document: US

NENP Non-entry into the national phase

Ref country code: DE

WWE Wipo information: entry into national phase

Ref document number: 2012788245

Country of ref document: EP

ENP Entry into the national phase

Ref document number: 20147015156

Country of ref document: KR

Kind code of ref document: A