Engineering system and method of automatic controlling of the energy efficient mode of pipeline transporting process of fluid medium

Abstract

An engineering system, comprising a pumping system, a wave system, a system of measuring devices, and a dynamic controller, configured for generating pressure waves with controlled frequencies, and their propagating through turbulent flow in the suction and discharge pipelines of the pumping system. The engineering system performs a method of automatic controlling the energy efficient mode of pipeline transporting of fluid medium. The dynamic controller regulates and adjusts in real time optimal values of frequencies of generating pressure waves, providing minimum values of hydrodynamic resistance of turbulent flow in the discharge pipeline, for any average value of the flow velocity, compared to the hydrodynamic resistance of the turbulent flow having the same velocity, but not interacting with pressure waves. Automatic controlling the process of generating pressure waves with optimal frequencies by the dynamic controller, ensures the maximum energy efficiency of the process of pipeline transporting of fluid medium.

Claims

1. An engineering system for pipeline transporting process of fluid medium comprising systems: a pumping system, comprising: a pump, having an inlet and an outlet, a suction pipeline and a discharge pipeline, wherein the inlet is connected to the suction pipeline and the outlet is connected to the discharge pipeline, a motor, coupled to a pump, the pumping system is configured to pump flow medium through the suction pipeline, the pump, and the discharge pipeline; a wave system, comprising: a pressure waves generator, a motor, connected to the pressure waves generator, connecting pipes connecting the pressure waves generator parallel to the pump, to the suction pipeline and the discharge pipeline, the wave system is configured for generating pressure waves; a system for generating pressure waves, formed from the interconnected the pumping system and the wave system, configured in such a way, that during the operating pumping system and wave system, pressure waves are generating by pressure waves generator, and propagating through the medium flow in the connecting pipes, the suction pipeline, and the discharge pipeline; a system of measurement devices, comprising: measurement devices set on the elements of the system for generating pressure waves, connected to an input module of a dynamic controller, and configured for measuring current values of operating parameters of the turbulent flow, of the motor of pressure waves generator, and of the pressure waves; a dynamic controller, connected to the feedback of the system for generating pressure waves, comprising: a microprocessor, a memory, an input module, an output module, output-controlled devices, a programming device; the input module is connected to measurement devices; the microprocessor is interconnected to the input module, the memory, the programming device, and the output module; the output module is connected to the output-controlled devices; the dynamic controller is configured for providing in real-time the automatic controlling of energy efficient mode of the process of pipeline transporting of fluid medium through the pumping system.

2. The engineering system according to claim 1, wherein the ratio of the cross-sectional area of the connecting pipe to cross-sectional area of the discharge pipe is less than 0.04.

3. The engineering system according to claim 1, wherein each of connecting pipe of the pressure waves generator is equipped with a controlled valve.

4. The engineering system according to claim 3, wherein each control valve is connected to the input module of the dynamic controller.

5. A method of automatic controlling of energy efficient mode of the process of pipeline transporting of fluid medium performed by the engineering system, comprising the steps of: computing the optimal average values of frequencies of pressure waves interacting with turbulent flow in a pipe, providing the minimum average values of the hydrodynamic resistance coefficient for a discrete set of current average values of the flow velocities, covering an entire practical range; experimental checking and correcting of the computed minimum average values of the hydrodynamic resistance coefficient at average values of optimal frequencies of generated pressure waves, for the discrete set of said average values of the flow velocities; compiling the baseline file containing: a discrete set of said average values of the flow velocities; a discrete set of the average values of computed and experimentally corrected minimum values of the hydrodynamic resistance coefficient at average values of optimal frequencies of generated pressure waves for said discrete set of current average values of the flow velocities; uploading the baseline file in the memory of the dynamic controller; measuring current values of the flow velocity, current of the pressure waves generator motor, optimal frequencies of generated pressure waves, and pressures in the discharge pipeline and in the suction pipeline; transmitting the current values of said measured parameters to the input module; converting the transmitted current values of said measured parameters in the digital format by the input module; transmitting the digital format of said measured parameters to the microprocessor; providing receiving, recording, processing, calculating, and saving current values of the flow velocity, current of the pressure waves generator motor, optimal frequencies of generated pressure waves, and pressures in the discharge pipeline and in the suction pipeline, by the microprocessor; averaging the received current values of the flow velocity, current of the pressure waves generator motor, optimal frequency of generated pressure waves, pressures in the discharge pipeline and suction pipeline, and calculated parameters at each specified time interval; compiling the file of average parameters containing updated average values of the flow velocity, current of the pressure waves generator motor, optimal frequencies of pressure waves, and pressures in the discharge pipeline and in the suction pipeline; selecting the row in the baseline file, considering pressure waves action on the turbulent flow, in which the average value of the flow velocity is closest to the last updated average value of the flow velocity from the file of average parameters; transmitting average value of the current of pressure waves generator motor from the selected row of said baseline file to the output module; converting the received average value of the current of pressure waves generator motor into analog signal; transmitting the value of analog signal to the automatic controlled-device; changing the value of operating current of the motor of pressure waves generator motor to the transmitted optimal average value of the current, from the baseline file, corresponding to current value of the flow velocity in the discharge pipeline, optimal frequency of pressure waves, and minimum value of the hydrodynamic resistance coefficient; providing the energy efficient mode of the process of pipeline transporting of fluid medium.

6. The method according to claim 5 wherein the correcting amplitude of generated pressure waves is performed by variation of the valve opening angle.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

(1) To facilitate a complete understanding of the proposed invention, it is useful to simultaneously analyze the text of the Detailed Description and accompanying the proposed drawings.

(2) FIG. 1 is a schematic representation of the system of generating pressure waves 15, formed from the interconnected pumping system 1 and wave system 7.

(3) FIG. 2 is an architectural scheme of the engineering system for pipeline transportation of fluid medium in energy efficient mode.

(4) FIG. 3 is an example of schematic representation of the asymmetric compression pressure wave P.sup.+(t) and the rarefaction pressure wave P.sup.?(t), forming simultaneously in the pressure waves generator.

(5) FIG. 4 is a schematic representation of an example of the change of pressure vs frequency in the asymmetric rarefaction pressure wave P.sup.?(t).

(6) FIG. 5 is a schematic representation of the system of measurement devices, set on the elements of the engineering system, connected to a dynamic controller.

(7) FIG. 6 is a schematic representation of an example of the discrete function obtained as a result of measuring flow rate in the discharge pipeline with the updated time interval ?t. a), r=0; b), c), r=1.

(8) FIG. 7 is a schematic representation of an example of the discrete function obtained as a result of measuring frequencies of pressure waves corresponding to current values of the flow rate in the discharge pipeline, with the updated time interval ?t. a), r=0; b), c), r=1.

(9) FIG. 8 is a simplified block diagram of the dynamic controller.

(10) FIG. 9 is an example of graphical presentation of the head-flow curves of the pump and of the piping system.

(11) FIG. 10 is a schematic representation of flow domain, consisting of the volume of a portion of the discharge pipe 6, and conjugated with a volume of the portion of connection pipe 11.

(12) FIG. 11 is a schematic representation of the possible location of the average measured values of the current of the pump motor, a), current of the PWG motor, b), energy efficiency, c) relatively to average values: T.sub.Wj, T.sub.PWj, ?.sub.fWj, accordingly.

(13) FIG. 12, a-b is a flowchart illustrating a method of automatic controlling of the energy efficient mode of the process of pipeline transporting of fluid medium operation with the use of the engineering system of FIG. 2.

(14) FIG. 13 is a flowchart illustrating automatic correction of the amplitude of pressure waves providing the energy efficient mode during the process of pipeline transportation of fluid medium.

DETAILED DESCRIPTION OF THE INVENTION

(15) The current invention provides an engineering system and method of automatic controlling of the energy efficient mode of pipeline transporting process of fluid medium.

(16) The detailed description offered below does not consider all the details of each of the systems that comprise the engineering system, but sufficient, for understanding the operation of the engineering system. In the detailed description of the method, some nonessential actions are also omitted, which do not affect the essence and sequence of the steps of the practical realization of the method.

(17) In modern pumping systems the pipeline transportation of fluid media is carried out by a pump that provides a stationary pressure gradient along the discharge pipeline.

(18) Attempts to decrease the hydrodynamic resistance of the turbulent flow in a pipeline with the use of a source of additional perturbations of the flow velocity or pressure, such as artificial turbulators installed on the inner surfaces of the pipes that have a local effect. Therefore, they did not find industrial applications.

(19) The use of negative modulation of the flow-forming energy makes it possible to generate the controlled pressure waves of small amplitudes and low frequencies, and ensure their propagation through the flow of the transported medium in long pipelines.

(20) The interaction of pressure waves of suitable amplitudes and frequencies with the turbulent medium flow in the pipeline, ensures a decrease in hydrodynamic resistance. Therefore, the development of a method for automatic controlling of the energy efficiency of the process of pipeline transporting of flow medium by controlling amplitudes and frequencies of generated pressure waves is fundamental for industrial applications.

(21) To facilitate a complete understanding of the proposed invention, it is useful to study simultaneously the text of the Detailed Description, accompanied by the proposed drawings.

(22) The current invention provides an engineering system and a method of automatic controlling of the energy efficient mode of pipeline transporting process of fluid medium. The engineering system, performing the energy efficient mode, is designed for pipeline transportation of fluid medium and comprises the following interconnected systems: a pumping system; a wave system; a system for generating pressure waves; a system of measurement devices; a dynamic controller.

(23) FIG. 2 is an architectural scheme of an engineering system for pipeline transporting of fluid medium in an energy efficient mode. The main system in the pumping system 1, whose structure is shown in FIG. 1. The pumping system comprises a pump 2, a motor 3, a shaft 4, a suction pipeline 5 connected to an inlet of the pump, and a discharge pipe 6 connected to an outlet of the pump. The pump 2, connecting shaft 4, and motor 3, form a pump unit.

(24) A wave system 7, FIG. 1 comprising a pressure waves generator 8, a motor 9, a shaft 10 connecting the motor with the PWG, and connecting pipes 11 and 12 with control valves 13 and 14, connected to the discharge and suction pipelines of the system 1. The PWG 8, motor 9 and connecting shaft 10 form a wave unit.

(25) The wave system is designed to provide the set dynamic law of change of connecting surface over time, S.sub.C(t) during a period, T. The law S.sub.C(t) defines the law of change of pressure over time in the generated wave, P.sup.?(t).

(26) Calculated values of parameters of the generated waves depend on: physical properties of the fluid medium (density ?, viscosity ?, compressibility ?, temperature t.sub.f); operating parameters (flow velocity ?, suction pressure P.sub.sp, and discharge pressure P.sub.dp on the pump); geometrical parameters (diameters of the suction and discharge pipeline and of connecting pipes).

(27) Based on these parameters, the geometric shapes, and sizes of profiles of the openings in the walls of coaxially located cylindrical elements of the PWG are selected and calculated. During rotation of the inner cylinder relative to the outer cylinder, the intersection of the surfaces of the specified openings is carried out (details, U.S. Pat. No. 10,233,952). As a result of such intersecting, a connecting surface is formed, S.sub.C(t). It connects two volumes belonging to PWG: inner, with the suction pressure, P.sub.S and outer, with the discharge pressure, P.sub.D. During the frontal time, t.sub.F the value of connecting surface increases from zero to maximum value, of the given law, S.sub.F(t): t=0 . . . t.sub.F. During the back time, t.sub.B the value of connecting surface decreases from maximum value to zero, according to a given law, S.sub.B(t): t=t.sub.F . . . T. The law of change of the connecting surface S.sub.F(t) and S.sub.B(t) over time, pressures P.sub.S and P.sub.D, and a speed of the PWG motor, ? define the law of change of pressure over time in the generated pressure waves, P.sup.?.sub.F(t) and P.sup.?.sub.B(t), t=0 . . . t.sub.F . . . T.

(28) A compression wave is generated in a similar way.

(29) A system for generating pressure waves comprising the interconnected pumping system 1 and the wave system 7, schematically presented in the FIG. 2. PWG 8 is connected parallel to the pump 2 through connecting pipes, 11 and 12, to the discharge 6 and suction 5 pipelines. Such configuration is designed to provide the generating pressure waves, and to propagate through the flow, along the connecting pipes 11 and 12, and onward, along suction and discharge pipelines, 5 and 6. During the process of pipeline transporting of fluid medium, in the suction and the discharge pipelines, average pressures <P.sub.S> and <P.sub.D>, and flow velocities, <?.sub.S> and <?.sub.D>, are established correspondingly. Thus, the necessary condition for generating pressure wave is the inequality: P.sub.D>P.sub.S.

(30) Physical meaning and design choice of the slope angles ? and ? between connecting pipes 12 and 11 and the suction and discharge pipes, 5 and 6 respectively, are presented by B. Pierre (Pressure Waves in Pipelines and Impulse Pumping. Physical Principles, Model Development and Numerical Simulation, Department of Petroleum Engineering and Applied Geophysics, Norwegian University of Science and Technology, Trondheim, December 2009).

(31) The process of generating waves is carried out as follows: suppose, that at the initial moment of time the value of the connecting surface is zero, S.sub.C(t)=0. When the motor 3 of a pump 2 is turned on, this system operates as a traditional pumping system 1 for transporting of fluid medium through the suction pipeline 5, the pump 2, and the discharge pipeline 6. When the PWG 8 is turned on, connecting surface becomes greater than zero, changes over time, linking volumes of medium with different pressures: inner volume with the suction pressure, P.sub.S and outer volume with the discharge pressure, P.sub.D. As a result of connecting volumes, an exchange of pressure impulse between them occurs. This process is accompanied by the simultaneous formation of compression and rarefaction pressure waves. The compression pressure wave P.sup.+(t) propagates through the connecting pipe 12 in the suction pipeline 5; the rarefaction pressure wave P.sup.?(t) propagates through the connecting pipe 11 in the discharge pipeline 6. A compression pressure wave, P.sup.+(t) is forming relatively to the suction pressure, P.sub.S, a rarefaction wave, P.sup.?(t) is forming relatively to the discharge pressure, P.sub.D, FIG. 3. Pressure waves P.sup.?(t) and P.sup.+(t) are generated only as a result of negative modulation of the flow-forming energy (details, U.S. Pat. No. 8,573,896). Physically, this means that a small part of the potential energy of the compressed flow, is converted into the wave energy, using the pump 2 and PWG 8, and is introduced back into the discharge flow. The turbulent flow generated this way is a perturbed weak wave flow.

(32) When the waves frequency changes, at constant pressures P.sub.S and P.sub.D between zones, the wave amplitude changes too. Indeed, the estimated value of the exchange of pressure impulse, passing perpendicularly to connecting surface S.sub.C(t) is:

(33) $\begin{array}{cc}{I}_p?{?}_0^T\left[P_D-P_S\right]*\frac{S_C\left(t\right)}{S_{\max }},& \left(6\right)\end{array}$ where S.sub.C(t) is the law of change of connecting surface over time during the period T, S.sub.C(max) is the maximum value of the connecting surface, corresponding to the maximum value of the pressure amplitude in the generated wave. FIG. 4 shows an example of a qualitative graphical interpretation of a rarefaction wave. Indeed, the amplitude of the generated pressure waves is inversely proportional to frequency, ?: P.sup.?(t.sub.F)?1/?. At the same time, the waveform also changes in each time interval, [0 . . . t.sub.F] and [t.sub.F . . . T].

(34) FIG. 4 clearly shows that, for example, a rarefaction wave with a frequency ?.sub.2, the total pump energy compression, [0 . . . P.sub.D] is represented as the sum of two components: a constant component, [0 . . . P.sup.?(t.sub.F2)] and a wave component, [P.sup.?(t.sub.F2) . . . P.sub.D]. The analogic energy structure is maintained for generated waves of any frequencies. Examples are waves with frequencies ?.sub.1 and ?.sub.3. The constant components of the pressure are respectively: [0 . . . P.sup.?(t.sub.F1)] and [0 . . . P.sup.?(t.sub.F3)], and of the wave components are [P(t.sub.F1) . . . P.sub.D] and [P.sup.?(t.sub.F3) . . . P.sub.D].

(35) The dynamic structure of turbulence and boundary layer of a medium flow of a given flow rate, <G.sub.j>??G(t), or flow velocity, <?.sub.i>???, interacting with pressure waves only of certain parameters, <P.sup.?>??P.sup.?(t), <?>???(t), <?>???(t), causes the physical effects described above that lead to decrease in the hydrodynamic resistance coefficient, ?.sub.W. The sign < > means the average value of parameter over a given averaging time interval ?t, containing a sufficient number of periods of pressure waves; ?G(t), ?P.sup.?(t), ?f(t), and ??(t) are dynamic deviations from the average values of the indicated parameters, during operating engineering system for pipeline transporting of fluid medium.

(36) Automatic control of the feedback between current value of the state of turbulent flow and parameters of generated waves is performed by the dynamic controller. The essence of this control is providing the continuous controlling frequency of generating pressure waves providing the decrease of hydrodynamic resistance of the current value of the flow rate.

(37) A system of measurement devices comprising measurement devices, installed on the elements of the engineering system, and connected to the dynamic controller, FIG. 5. This system is designed for measuring the current values of: the state parameters of fluid medium, a flow velocity, measured by a flow meter, an inlet pressure on a suction pipeline, pressures at the inlet and outlet of the discharge pipeline, a flow temperature; parameters of the pump unit components-current of the pump motor, speed of the pump impeller; parameters of the wave system-current of the motor of PWG, speed of the shaft of PWG; frequency of the generating pressure waves; the pressure waves parameters propagating in the discharge pipeline-amplitude, frequency, waveform.

(38) Each device is configured to measure the current value of parameter and to display on the screen the digital value of measured parameter in the standard or engineering units. For example: pressure can be displayed in Pa, Bar, atm; flow rate can be displayed in GPM, m.sup.3/s, m.sup.3/h.

(39) Speed of the pump impeller can be measured directly, by tachometer 25, FIG. 5, or calculated by using pump power characteristic and measured current of the pump motor by amperemeter 26.

(40) To ensure high accuracy of the frequency of the generated waves, in the wave system 7 a brushless motor 9 is used. It also smooths out the unevenness of the torque that occurs during operation of the PWG 8. The angular speed of the shaft 10 of the PWG 8 is determined by the torque-speed characteristic of the motor 9, measured by motor current, or by tachometer 27, FIG. 5. The system of measuring devices, the diagram of which is shown in FIG. 5, can be used in a hydrodynamic laboratory loop, for example, to adjust the parameters of generated pressure waves, compare and correct simulated results with experimental ones.

(41) Accuracy of measurement devices is an important requirement. Usually, digital measurement devices in the control system are used and most of them are already operating at any pumping station, except for measurement devices of the wave system. If among the measurement devices there are analog once, then the input module of dynamic controller converts analog signals into digital ones.

(42) During operation of the engineering system of the process of pipeline transporting of fluid medium, each measurement device shows instantaneous value of the current parameter, for example, flow rate G, pressure P.sub.D at the inlet of the discharge pipeline. All measurement devices are set to the same interval of time for updating values of measured parameters. In a particular case, the updated time interval of measured parameter can be changed. For example, for detailed analysis of the change of pressure over time in the wave, during a period, on a pressure sensor display, a shorter updated time interval is set.

(43) Power of the pump motor can be calculated using measured value of the current and standard voltage. This power can also be measured by watt meter. In this case, a watt meter is connected to the motor circuit. Usually, in a pumping stations both of these measurement devices are connected and installed on the control panel.

(44) The current values of the measured parameters change over time, under the influence of weak pressure and flow rate trends that exist in long pipelines, as well as more rough changes in flow rate caused by users. And in the case of the interaction of a turbulent flow with pressure waves, the dynamics and structure of the boundary layer are changed; the spectrum of turbulent fluctuations of the flow velocity and pressure, as well as macroscopic measured parameters, for example, the flow rate G and pressure P.sub.D, also are changed. These structural-dynamic changes in the turbulent flow, lead to change in turbulent viscosity and consequently, in the hydrodynamic resistance coefficient, ?.sub.W and the energy efficiency, E.sub.fW. The process of real-time controlling of the structural changes or dynamic state of the turbulent flow is determined by the frequency and amplitude of pressure waves interacting with the flow. Practically, this is carried out automatically by the PWG motor, based on the software data processing the measured current values of the flow parameters: flow velocity ?, measured by a flowmeter and calculated by ?=G/S; discharge P.sub.D and suction P.sub.S pressures, and hydrodynamic resistance coefficient, ?.sub.W.

(45) Usually, the updated time interval, ?t of the measured parameter(s), for example, of the flow rate, is quite short (of the order of tens ms) compared with the time interval ?.sub.W of responses of the output-controlled device and of the wave unit. Therefore, to control the wave unit using dynamic controller 17 with the output-controlled device 20, it is necessary to increase the duration of the time of the supplying current to the motor 9 of the PWG to the order of ?.sub.W. This requirement also relates to averaging time interval of the basic measured independent argument, flow rate. For satisfying the above specified condition, consider the next relation between time intervals: [m ?t]>?.sub.W, where m is an integer. The discrete function G.sub.Wj (i ?t), i=1, 2 . . . m, measured at the points ?t, 2 ?t . . . m ?t is converted into a continuous function on the same averaging time interval [m ?t], as a result of the simple averaging. FIG. 6:

(46) $<G_j>?{\stackrel{?}{G}}_j=\left(1/m\right)*{.Math.}_{i=1}^{i=m}{G}_j\left(i*\mathrm{?t}\right).$ The averaging time interval [m?t], or number m, is determined from practical considerations. It depends on the flow rate in the discharge pipeline and on the frequency of generated waves. The sign < > and x bar means the average value of the indicated parameter, over the given averaging time interval [m?t]. Values of the discrete function G.sub.Wj (i ?t) at the points i=1, 2 . . . m, form a scenario. Examples of scenarios: a)-c) in the FIG. 6. Differences in the values of the functions ?G.sub.j[(rm+i)?t]=<G.sub.j>?G.sub.j[(r m+i)?t]at points (r m+i)?t, i=1 . . . m, r=0, 1, 2 . . . N, are the deviations of G.sub.j [(rm+i)?t] from the average value. Thus, the flow rate in any moment of time t.sub.i, is defined as:
G.sub.1[(rm+i)?t]=<G.sub.j>??G.sub.j[(rm+i)?t].(7) t.sub.i?[(rm+i)?t], j=1, 2 . . . n, r=0, 1, 2 . . . N is the sequence of numbers of the averaging time intervals [m?t]. A starting of time count is performed from some reference point of counting. For example: i=1, 2 . . . m. r=0: t=1?t, 2?t . . . m ?t; ?G.sub.j(1?t), ?G.sub.j (2*?t) . . . ?G i (m ?t). r=1:t=(m+1)?t, (m+2)?t . . . (m+m)?t; ?G.sub.j [(m+1)?t], ?G.sub.j[(m+2)?t] . . . ?G.sub.j[(m+m)?t]. r=2: t=(2m+1)?t, (2m+2)?t . . . (2m+m)?t; ?G.sub.(j?1)[(2m+1)?t], ?G.sub.(j?1) [(2m+2)?t] . . . ?G.sub.(j?1)[(2m+m)?t].Math.m<N.

(47) FIG. 6 is an example of the graphical interpretation of the definition (7): a) r=0; b), c) r=1 from the start of counting time. Obviously, scenarios b) and c) cannot be performed simultaneously. After the scenario a) is performed, either scenario b) or c), one of them can be performed. +?G(t) is defined in the points t i, in which, G(t)>(<G.sub.j>) in the averaging time interval [m?t], r=0, 1. Points above the lines <G.sub.j> and <G.sub.j+1> on the FIG. 6. ??G(t) is defined in the points t in which G(t.sub.i)?(<G.sub.j>) on the averaging time interval [m?t], r=0, 1. Points below the lines <G.sub.j> and <G.sub.j+1> on the FIG. 6.

(48) To assist with the understanding of the discussion of the present invention, the following definitions are provided: a controlled system, is the turbulent flow in a discharge pipeline, described by the state parameters: flow velocity, temperature, turbulence structure. a control function, is a pressure wave of the given amplitude and frequency, generated by PWG, and propagated through the turbulent flow in a discharge pipeline.

(49) The need to use the flow velocity in the process of automatic control of energy efficient mode is obviously. A pressure wave interacts with a turbulent flow in the discharge pipeline, which is uniquely described by the flow velocity, not by the flow rate. It is the flow velocity that determines structure and dynamic state of the boundary layer and intensity of turbulent pulsations, which interact with the pressure waves and lead to decrease in the normal components of pulsations, change the specter of turbulent pulsations, anizotropization of turbulence, i.e., to partial laminarization. Ultimately, all these phenomena lead to decrease of hydrodynamic resistance.

(50) However, the main parameter of any pipeline transportation system is the flow rate. This is a basic technological parameter and easily measured with high accuracy. Therefore, in the proposed invention, all measured parameters depend on the flow rate, explicitly or implicitly. However, the physical parameters that underlie the interaction of pressure waves with turbulent flow leading to a change in the state of the controlled system should be considered the flow velocity, both in the description of the patent and in the claims.

(51) The structure of the turbulent flow, and hence, the hydrodynamic resistance coefficient changes when pressure waves interact with the turbulent flow in the pipeline. The frequency of pressure waves controls the internal structure of the turbulent flow. At a certain frequency, the hydrodynamic resistance coefficient, corresponding to current value of the flow rate, reaches a minimum. This frequency is called optimal and corresponds to the energy efficient mode of the pipeline transporting process of flow medium.

(52) Presentation of the measured flow rate in the form (7), is the first step towards the implementing the process of automatic controlling of the energy efficient mode of pipeline transporting process of flow medium. Measuring the flow rate in real time, means the measuring of the state parameters of turbulent flow. State parameters define the structure of turbulent flow and boundary layer, which determine the value of the hydrodynamic resistance coefficient, ?.sub.W. The remaining measured and computed energy and dynamic parameters of the turbulent flow, presented in TABLE 2, are considered as functions of the flow rate. FIG. 7 shows an example of such a functional relationship: the frequency of generated pressure waves, as a function of the flow rate. This discrete function is obtained as a result of measuring frequencies of pressure waves, corresponding to current values of the flow rate in the discharge pipeline, with the updated time interval ?t, a), r=0; b), c), r=1.

(53) Power of the PWG motor can be calculated using measured value of the current and voltage, usually low constant voltage. This power can also be measured by wattmeter. In this case in the circuit the of the PWG motor a watt meter is connected.

(54) Below, we will consider only such average functions of the averaging time interval [m ?t], which are obtained as a result of measuring instantaneous values of parameters: flow rate, waves frequency, speed of the pump impeller etc. The value of ratio [m ?t]/?.sub.W acquires an important role in the process of automatic controlling of the energy efficient mode.

(55) The measured current values of parameters are transmitted to the input module of dynamic controller.

(56) A dynamic controller 17, FIGS. 2, 5 is a system, FIG. 8, comprising a microprocessor 33, a memory 34, a programming device 35, that can be a computer, an input module 18, an output module 19, and output-controlled devices 20-22. The dynamic controller is an interconnected system. It is connected to the system of measurement devices, FIG. 5 and to the system of generating pressure waves 15, FIG. 1. Also, it can be connected to a computer 23 or other interface device, 24. The dynamic controller 17 is a feedback controller connected to a controlled system to control the current state of turbulent flow in the discharge pipeline 6 by variation of the frequency of generated pressure waves by changing the speed of the PWG motor 9. All of current values of measured parameters are transmitted by measurement devices to the input module 18. The input module converts all measured parameters into digital format and transmitted to the microprocessor 33. The microprocessor processes the received data and, in accordance with the instructions, transmits it to the output module in digital format. The output module 19 converts all received digital signals into analog and transmits them to the output-controlled devices 20-22. Thus, the microprocessor receives data from the measuring devices, processes them in accordance with the program instructions, and performs controlling the state of the turbulent flow in the discharge pipeline using output-controlled devices 20-22.

(57) A simplified block diagram of the dynamic controller is presented in the FIG. 8. The display, power supply and its connections with the microprocessor, memory, input module, and output module are not shown. The dynamic controller 17 is designed to provide the real-time automatic controlling of the state parameters of the fluid medium transporting in the discharge pipeline in the energy efficient mode. Such mode is maintained by continuous controlling of the speed of PWG motor which determines the frequency of generated waves, interacting with fluid medium in a discharge pipeline. The microprocessor 33 receives from the input module 18 current values of measured parameters in digital format, processes them, and performs in real-time the energy optimization, by minimizing the value of the hydrodynamic resistance coefficient for the current value of the flow rate, by adjusting frequency of generated pressure waves. The mentioned optimization is carried out using a controlled output device 20 of dynamic controller, by controlling the speed, ? of the PWG motor 9.

(58) The dynamic controller 17 can be used for the pump control, for example, with the constant pump RPM using a controlled-device 22, without operating the wave system.

(59) The dynamic controller is configured also for manual operating mode. This means that the user can manually change or fix the current of the pump motor, current of the motor of the PWG, controlling the frequency of generated waves, and valve opening angle of the control valves.

(60) Before proceeding to the description of the method of energy efficient mode of pipeline transporting process of fluid medium, let us consider the database that is the basis for development of the program to implement the method of automatic controlling of this mode.

(61) The essence of this database is to simulate the turbulent motion of a fluid medium using the Navier-Stokes equation, under the action of the constant pressure and of the perturbed pressure by pressure waves, in the flow domain belonging to a straight horizontal pipeline section with a given constant cross section, FIG. 10. The purpose of the simulation is to compute the optimal frequencies of pressure waves that provide the minimum value of the hydrodynamic resistance coefficient of a turbulent flow at the current value of velocity. Modeling is carried out for a set of discrete values of flow velocities belonging to the entire practical range, {?.sub.j}, j=1, 2 . . . n. ?.sub.j=G.sub.j/S. As a result, we obtain two sets of solutions, for the constant pressure, P.sub.D=const: profile of the flow velocity, ?.sub.j(r), 0?r?R, R is the pipe radius; pressure P.sub.i on the inlet cross-section S.sub.i of the flow domain; pressure P.sub.O on the outlet cross-section S.sub.O of the flow domain; average value of the kinetic energy of turbulent motion; radial and longitudinal distributions of the turbulent kinetic energy; the wall shear stress; hydrodynamic resistance coefficient, ?.sub.j; and pressure losses on the flow domain; for a perturbed pressure by pressure waves P.sup.?(t), FIG. 3: profile of the average flow velocity, ?.sub.Wj(r), 0?r?R; average value of the pressure P.sub.Wi on the inlet cross-section S.sub.i of the flow domain; average value of the pressure P.sub.Wo on the outlet cross-section S.sub.O of the flow domain; average value of the wave's frequency providing minimum value of the hydrodynamic resistance coefficient, ?.sub.optj (?.sub.Wi); average value of the kinetic energy of turbulent motion; radial and longitudinal distributions of the turbulent kinetic energy; average value of the wall shear stress; average value of the hydrodynamic resistance coefficient, ?.sub.Wj (?.sub.optj(?.sub.Wj)); average value of the pressure losses on the flow domain.

(62) These two sets constitute the database for the baseline file which is used for the comparative analysis of the states of controlled system when the control function is changing.

(63) Next, the energy efficiency E.sub.f and ?.sub.fW are calculated and compared. Modeling is carried out using a sufficiently accurate computational tool, for example, Large Eddy Simulation (LES), which is widely used in Computational Fluid Dynamics (CFD). Calculations are carried out over the entire practical range of the flow velocity and consequently, of the Reynolds numbers.

(64) Obtaining the Baseline Data for the Turbulent Flow Under the Action of the Constant Discharge Pressure.

(65) The flow domain of horizontal cylindrical form, FIG. 10 is limited by inner cylindrical surface S.sub.F, inlet surface S.sub.i and outlet surface S.sub.O, both perpendicular to longitudinal axis of symmetry. The turbulent motion of the fluid medium is simulated under the action of the constant pressure over the entire practical range of the flow velocity, {?.sub.j}, j=1, 2 . . . n; the corresponding values of the hydrodynamic resistance coefficient, ?.sub.J are computed. The discrete set of the values {?.sub.1, ?.sub.2 . . . ?.sub.n} uniformly covers the entire practical range of the flow velocity, that is ?.sub.j+1??.sub.j=??=const, j=1, 2 . . . n. In the simulation process the flow velocity is used. Flow rate is directly related to the cross-sectional area of the flow domain (pipe) and can be measured easily and accurately. The discrete set of the values of the computed flow velocities {?.sub.1, ?.sub.2 . . . ?.sub.n} is equivalent to the discrete set of the values of the measured flow rates {G.sub.1, G.sub.2 . . . G.sub.n}. And G.sub.j=?.sub.j*S, G.sub.j+1?G.sub.j=? G=const1. If LES is performed on another computer, then the results are uploaded in the programming device 35. The simulation results are verified experimentally.

(66) Experimental verification of the hydrodynamic resistance coefficient, ? obtained as a result of the LES, is carried out as follows: a horizontal experimental section of the pipeline with a length L and a diameter D is selected, at the distance of more than 35 pipe diameters downstream from the conjugation surface S.sub.D, FIG. 10. Pressure sensors are installed at the ends of the experimental section of pipeline. The pump is turned on, the flow rate is fixed, the pressures at the ends of experimental section of the pipeline are measured, and the pressure losses ?P.sub.? are computed. The flow rate G.sub.j is adjusted discretely with a constant step, ? G.sub.j; j=1, 2 . . . n. Similarly, the pressure losses ?P.sub.?j are computed. The experimental results of the pressure losses are compared with theoretical results, according to the Dacy-Weisbach equation. Theoretical and experimental results are corrected.

(67) In the process of experimentation, the error of all measuring devices is determined. The final corrected results are presented in TABLE 1.

(68) TABLE-US-00001 TABLE 1 Baseline data. P.sub.D = const. j G.sub.j ?.sub.j Re.sub.j Ps.sub.j P.sub.Dj I.sub.j ?.sub.j N.sub.Pj ?.sub.j V.sub.j E.sub.fj 1 G.sub.1 ?.sub.1 Re.sub.1 Ps.sub.1 P.sub.D1 I.sub.1 ?.sub.1 N.sub.1 ?.sub.1 V.sub.1 E.sub.f1 2 G.sub.2 ?.sub.2 Re.sub.2 Ps.sub.2 P.sub.D2 I.sub.2 ?.sub.2 N.sub.2 ?.sub.2 V.sub.2 E.sub.f2 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . q G.sub.q ?.sub.q Re.sub.q Ps.sub.q P.sub.Dq I.sub.q ?.sub.q N.sub.q ?.sub.q V.sub.q E.sub.fq . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . n G.sub.n ?.sub.n Re.sub.n Ps.sub.n P.sub.Dn I.sub.n ?.sub.n N.sub.n ?.sub.n V.sub.n E.sub.fn

(69) These results constitute the baseline data (reference data) for the turbulent motion of the fluid medium in a pipeline, under the action of stationary pressure, P.sub.D=const. In essence, the baseline data at constant pressure is the table function in which, the measured independent argument is the flow rate, G.sub.j, j=1, 2 . . . n. This is the basic measurable technological parameter. The remaining parameters presented in the columns of the TABLE 1 are functions, dependent on the flow rate. Moreover, some of them are measured, and some are computed using values of measured parameters. For example, the flow velocity, ?.sub.j=G.sub.j/S, S is the cross-sectional area of the flow domain (pipe); S=S.sub.i=S.sub.o, FIG. 10; Reynolds number, Re.sub.j=A*G.sub.j, A=4/(?.sup.Dv); D=2R is the pipe diameter; V.sub.j is the volume of the pumped fluid, during the given time interval, for example, ??.sub.G; Ef.sub.j is the energy efficiency. Subscript q is integer, q=n/2.

(70) This is the baseline file, considering the action of constant pressure on the turbulent flow in the discharge pipeline. It is developed by programming device 35, is uploaded and is stored in the memory 34 of the dynamic controller 17.

(71) It is obvious, that the main measurable flow parameter, flow rate G.sub.j (TABLE 1), is easily generalized to other similar flows, using the Reynolds similarity criterion, Re.sub.j=4 G/?D?=?.sub.jD/?. For example, if the kinematic viscosity of the main flow and of a similar flow are equal, then from the condition of similarity of flows it follows: ?.sub.j/D.sub.j=?.sub.?/D.sub.?. This condition can be generalized to flows with other pipe diameters and the entire practical range of the flow velocity. Subscript ? relates to parameters of the similar flow. Pressures P.sub.S and P.sub.D, as well as parameters of the pump, Q j and the pump motor, I.sub.j, N.sub.Pj, clearly depend on the flow rate, diameter and length of the pipeline, and also on its configuration, bends and elevation. Parameters Psi, P.sub.Dj, I.sub.j, ?.sub.j, N.sub.Pj, G.sub.j, are intrinsic parameters of the pumping system. They serve to determine the energy and power parameters of the pumping system, for example, calculating the energy efficiency E.sub.f, of the pipeline transporting process of flow medium, using (1). They also serve to monitor the operating dynamics and stability of the pumping system. Therefore, intrinsic parameters of the pumping system are added as data points to the baseline file, TABLE 1.

(72) The increase of energy efficiency of perturbed turbulent flow by pressure waves is carried out not only explicitly, due to the decrease of hydrodynamic resistance, but also implicitly, due to the change in flow rate. FIG. 9 shows an example of the graphical presentation of the head-flow curve of the pump, curve 1, and of the system, curve 2 relating to unperturbed flow with the volumetric flow rate G. When the PWG is turned on, the interaction of pressure waves with turbulent flow leads to decrease of hydrodynamic resistance of the system. Due to internal changes in turbulent flow, under the influence of pressure disturbances, the system head-flow curve increases curvature and takes the position 3; 4 is the static head, H.sub.S. Pressure losses ?P.sub.?W are decreased: ?P.sub.?W<?P.sub.?. This decrease is accompanied by increase in the flow rate, from G.sub.? to G.sub.?W. Thus, the full contribution to the increase of energy efficiency is also determined by the increase in flow rate, ?G.sub.?W=G.sub.?W?G.sub.?.

(73) Obtaining the Baseline Data for the Turbulent Flow Under the Action of Perturbed Discharge Pressure by Pressure Waves.

(74) The purpose of this baseline data is to compute the optimal values of frequencies, <?.sub.opt j>.sub.opt and amplitudes, P.sup.?(t.sub.F) opt of pressure waves, interacting with the turbulent medium flow, and to provide the minimum values of the hydrodynamic resistance coefficient, ?.sub.W compared with the hydrodynamic resistance coefficient, A of the same discharge pipeline, of the same current values of the flow velocity, but under the action of the constant pressure. To achieve this purpose, the turbulent motion of the medium flow in the discharge pipeline, interacting with pressure wave, P.sup.?(t), FIG. 3, is simulated, by using LES. During the simulation, the wave frequency and amplitudes vary. Calculations are carried out over the entire practical range of the flow velocity and consequently Reynolds numbers.

(75) The configuration of flow domain is selected in such a way, that the physical picture of the simulated interaction of pressure waves with the turbulent flow in the pipeline is as close as possible to the real case scenario.

(76) Flow domain, FIG. 10, consists of the volumes of portions of the discharge pipe 6 and of the connection pipe 11, connected through a conjugated surface, S.sub.J. Surfaces bounding the flow domain are: an inlet S.sub.I and an outlet S.sub.O, the inner surface of the discharge pipe S.sub.F, a conjugated surface S.sub.I, an inner surface of connecting pipe, S.sub.CP, and a surface of cross-sectional area of connecting pipe, S.sub.W. Pipe volumes communicate only through the conjugated surface S.sub.J. Therefore, the boundary conditions on the surface S.sub.W essentially set the boundary conditions on the surface S.sub.J, through which the wave energy is transmitted to the turbulent flow in the discharge pipeline. And besides: S.sub.I=S.sub.O, S.sub.W/S.sub.D<0.04.

(77) To identify the effects of decreasing, and further minimizing the hydrodynamic resistance of the flow domain, operating boundary conditions are used. On the surfaces bounding the flow domain the following boundary conditions are inserted: on the S.sub.I, a mass flow rate; on the S.sub.O, pressure; on the surfaces S.sub.F and S.sub.CP, the no-slip conditions of the flow velocity; on the surface SW, the pressure wave action, P.sub.S.sup.?(t). Pressure P.sub.S.sup.?(t) is changed by the wave law relatively to the discharge pressure P.sub.D, FIG. 3. During the LES simulating process, the boundary conditions on the surface S w are set, by User Defined Function. This function completely defines the rarefaction wave P.sub.S.sup.?(t) on the surface S.sub.W: the form, amplitude, and frequency. Boundary conditions on the surface S.sub.W essentially set the boundary conditions on the surface S.sub.J, through which the wave energy is transmitted to the turbulent flow in the discharge pipeline.

(78) The essence of the LES simulation is to compute the optimal frequency, ?.sub.opt j=<?.sub.optj>???.sub.j corresponding to current value of the flow velocity ?.sub.Wj [(rm+i)?t]=<?.sub.Wj>???.sub.Wj[(rm+i)?t)] or flow rate, G.sub.Wj[(rm+i)?t]=<G.sub.Wj>T ?G.sub.Wj [(rm+i)?t)], which provides the minimum value of the hydrodynamic resistance coefficient, ?.sub.W=?.sub.Wmin, i.e. ?.sub.Wj/?.sub.j=min, and ?.sub.Wj=?.sub.wj.sup.min, j=1, 2 . . . n. In this case the following inequality is satisfied:
[?.sub.j?(?.sub.wj.sup.min)]<?.sub.?.(8) where ?.sub.? is the minimum value of the difference in the hydrodynamic resistance coefficients, valid for any value of the flow rate, G.sub.Wj:G.sub.Wj=G.sub.j.

(79) Note, that inequality (8) is satisfied in a certain range of the flow rate changes. In practice, this means that there is a certain range of waves frequencies changes corresponding to a certain range of flow rate changes in which inequality (8) is satisfied. Under the action of pressure waves, the stabilized state of turbulent flow is provided. This is a fundamental principle in the process of automatic controlling of dynamic systems, that the turbulent flow is interacting with pressure waves.

(80) It is this correspondence between the frequency ?.sub.opt j and the flow velocity, ?.sup.wj, that plays a central role in the process of automatic controlling of the energy efficient mode of pipeline transporting process of fluid medium. This process is carried out in real-time by the dynamic controller 17 with the use of its feedback between the state parameters of the turbulent flow in the discharge pipeline, ?.sub.wj and the frequency of generated pressure waves, ?.sub.optj. Therefore, the table function ?.sub.optj(?.sub.Wj) (TABLE 2), is obtained as a result of the LES simulation. Then it is corrected experimentally. The final data ?.sub.Wj, ?.sub.optj(?.sub.Wj), ?.sub.min Wj form the functional dependences, ?.sub.optj(?.sub.Wj) and ?.sub.min Wj(?.sub.optj(?.sub.Wj)) which is the baseline file, considering the action of pressure waves on the turbulent flow.

(81) If LES is performed on another computer, then the results are uploaded into the programming device 35.

(82) The simulation results are verified and corrected on the laboratory hydrodynamic loop. Boundary conditions used on the surfaces S.sub.I and S.sub.W of the flow domain during the process of LES simulation are practically realized as follows: 1. the fluid flow and pressure applied to the S.sub.i surface is provided by a pump; 2. the frequency of the generated pressure waves, set on the surface S.sub.W, FIG. 10, is provided by the rotation speed of the pressure waves generator, which is controlled by the current in the motor 9, FIG. 1; 3. the amplitude and form of generated pressure waves, set on the surface S.sub.W, are determined by the shapes and sizes of the profiles of the surfaces of the openings in the walls of coaxially located cylindrical elements, providing connection between zones with different pressures P.sub.sp and P.sub.dp of the flow medium, during operation of the PWG (details U.S. patent Ser. No. 10/233,952).

(83) If the amplitude of the generated pressure waves obtained as a result of the LES is not enough to provide the minimum values of hydrodynamic resistance ?.sub.wj.sup.min and of energy efficiency, ?.sub.fW on the laboratory hydrodynamic loop, then it is corrected by non-essential variation of the valve opening angle ? of the controlled valve, for example 13, relatively to the default value ?.sub.D.

(84) Experimentally, the hydrodynamic resistance coefficients ?.sub.j and ?.sub.wj.sup.min are not measured directly, but implicitly, by measuring pressures at the ends of the experimental section of the pipeline, pump power or current in the pump motor. Energy efficiency E.sub.f and ?.sub.fW are calculated using measured parameters: flow rate, power or current in the pump motor, over a given averaging time interval.

(85) Inequality (8) imposes the following condition on energy efficiency:
[E.sub.fj?(?.sub.fwj.sup.min)]<?.sub.E,(9) where ?.sub.E is the lower acceptable limit of energy efficiency.

(86) Inequality (9) is based on experimental measured parameters and is used for laboratory testing of simulated results. It is also used by the program for providing the proposed method of automatic controlling of the energy efficient mode in an industrial environment, by the use of current measured parameters of the system for generating pressure waves. AE is set practically. Its estimated value is based on the configuration of the pipeline, physical properties of the flow medium, and of the flow rate.

(87) The sketch of laboratory hydrodynamic loop is similar to one represented on the FIG. 5, with non-essential changes. Instead of the dynamic controller 17 the simplified manual control panel is used. In addition, the suction pipeline is equipped with a valve for manual flow control; the experimental section is a part of the discharge pipeline, equipped with pressure sensors at the beginning and at the end; connecting pipes 11, 12 of the PWG, equipped with pressure meters, for detailed analysis of the law of change of pressure over time in the generated waves. During experimental testing of simulation results, the flow rate, current of the pump motor, and current of the motor of PWG, defining the frequency of generated waves, are varied manually. By default the waves amplitude, is set by the connecting surface S.sub.C(t), and can be non-essential corrected manually, by discrete varying of the valve opening angle ? of the control valve, for example, 13. For a more detailed analysis of the law of pressure change in the wave, a shorter updated time interval can be set on the pressure sensors. The optimal value of frequency ?.sub.opt Wj at a fixed flow rate G.sub.Wj ensures the minimum value of the hydrodynamic resistance coefficient, ?.sub.minWj, which is calculated according to the measured values of pressure at the ends of the experimental section of pipeline. Experiments are carried out for entire practical range of values of the flow velocities, {?.sub.Wj}, j=1, 2 . . . n. ?.sub.Wj=G.sub.Wj/S, same as for the baseline data for the constant pressure: G.sub.Wj=G.sub.j, j=1, 2 . . . n. G.sub.Wj+1?G.sub.Wj=?G.sub.W=?G=const1, FIG. 6. Entire practical range means that the flow velocity in a pipeline of a given diameter practically varies from the minimum to the maximum possible values. For a more accurate calculation of the optimal waves frequency, ?.sub.opt W and further experimental correction, it is required, that the intervals of the flow rate ?G.sub.W be relatively small, FIG. 6. The same requirement relates to the intervals of the flow velocity. In addition, the maximum and minimum deviations of two consecutive values of the flow rate must match, i.e.
G.sub.Wj+1??G.sub.Wj+1=G.sub.Wj+?G.sub.Wj,
?G.sub.Wj+1=?G.sub.Wj=?G.sub.W=const1, j=1, 2 . . . n.(10)

(88) Equality (10) is a fundamental condition in terms of the program for implementation of the method of automatic controlling of energy efficient mode of pipeline transporting process of the flow medium. The condition of cross-linking (10) at the ends of any two successive intervals ?G.sub.Wj and ?G.sub.W+1 is required. An example of this condition, for scenarios a) and c), G.sub.Wj+?G.sub.W=G.sub.Wj+1??G.sub.W is presented geometrically on the FIG. 6, the dotted line. Cross accent, x is the measured instantaneous values of the flow rate, ?G(i?t). Example: c), double line, G.sub.Wj??G[(m+1)?t], i=1.

(89) If the flow rate in the pipeline changes significantly, for example, for technological needs, then the average values of the flow rate change over time until another stationary state is stabilized, FIG. 6, c); k>1. To this continuous range of flow rates corresponds a continuous range of optimal waves frequencies, that provides the minimum value of the hydrodynamic resistance coefficient. However, while the deviations of the flow rate from any average value are the same, ??G.sub.W, (10), the corresponding deviations of the frequencies from the average value, ??.sub.opt Wj and ??.sub.optWj+1, are different. This is due to nonlinear dependence between the optimal frequencies of pressure waves and flow rate. Any two successive values optimal pressure wave frequencies of this nonlinear dependence (function) are expressed by the condition:
?.sub.optWj+?.sub.fWj??.sub.optWj+1???.sub.Wj+1
??.sub.Wj???.sub.Wj+1, j=1, 2 . . . n.(11)

(90) This means, that for different optimal frequencies, the intervals of deviation from the average value differ. This is schematically represented in FIG. 7. Cross accent, x is the measured instantaneous value of the frequency, ??(i?t). Example: double line, ?.sub.optWj??? [(m+1)?t], i=1.

(91) Thus, the results of the LES are corrected experimentally.

(92) In the process of experimental correction of simulated parameters in laboratory conditions and further in industrial environment, the number of experimental points m is much greater than indicated in FIGS. 6, 7.

(93) Experimental laboratory correcting of computed parameters is carried out as follows: if the diameter of the pipeline is equal to diameter of the simulated flow domain, then, the function ?.sub.optWj(G.sub.Wj) is corrected. And ?.sub.optWj (G.sub.Wj) is equivalent to ?.sub.optWj (?.sub.Wj); wherein G.sub.Wj=G.sub.j. The frequency is varied in the vicinity of the computed optimal value, at which the minimum value of ?.sub.min Wj is reached. This is the value of optimal frequency, ?.sub.optWj. To achieve the optimal value of frequency ?.sub.opt Wj, if necessary, the amplitude of the pressure waves is non-essential corrected by varying the valve opening angle ? of the controlled valve, for example 13. Other words, the function ?.sub.op wj (?.sub.wj) is corrected in the entire practical range of the flow velocity and the nonlinear functional dependence ?.sub.opt Wj (?.sub.Wj) is built.

(94) The deviations of the measured instantaneous values of parameters from the average are due to the dynamic changes of turbulent flow, actual process of generating waves, as well as the accuracy of measuring devices. The deviations of the measured instantaneous values of parameters from the average, for example, of the flow rate, ?G(i ?t), over an averaging time interval [m ?t], are located in a fairly wide range. A part of deviations belong to the interval, ?G(i?t)?[G.sub.Wj??G.sub.W, G.sub.Wj+?G.sub.W] adjacent to average value G.sub.Wj, FIG. 6. Other parts are located outside the interval ?G.sub.W, ?G(i ?t).Math.[G.sub.Wj??G.sub.W, G.sub.Wj+?G.sub.W]. They are located within the interval ?G.sub.W, which is adjoined to the average values of G.sub.Wj?1 or G.sub.Wj+1 and are satisfactorily approximated as:
1. G.sub.Wj, if ?G(i?t)?[G.sub.Wj??G.sub.W,G.sub.Wj+?G.sub.W];
2. G.sub.Wj?1, if ?G(i?t)<G.sub.Wj??G.sub.W;
3. G.sub.Wj+1, if ?G(i?t)>G.sub.Wj+?G.sub.W.(12)

(95) By substituting G.sub.W with ?.sub.W/S in expression (12), we obtain similar expressions for the average values of the flow velocity.

(96) Thus, the experimentally measured average values of G.sub.WX cover the entire practical range of flow rate, and of the flow velocity. Wherein, any measured value of current flow rate G.sub.WX is satisfactorily approximated by one of the conditions (12). Similar approximations are also valid for the flow velocity ?.sub.wX, control function ?.sub.opt and any other measured parameter.

(97) Such a complementary approach, theoretical and experimental, is aimed to defining the corrected control function ?.sub.opt and the controlled dynamic system, which is a turbulent flow in the discharge pipeline in the state with the minimum value of the hydrodynamic resistance coefficient, ?.sub.min W.

(98) All measured and calculated parameters are recorded in the baseline file. Computed and experimental results are corrected. The final corrected results are presented in the TABLE 2; for any j, j=1, 2 . . . n: G.sub.Wj=G.sub.j, or ?.sub.Wj=?.sub.j.

(99) These results constitute the baseline data for the turbulent motion of the fluid medium in a discharge pipeline, interacting with pressure waves. Part of the baseline data is obtained by recording parameters from measurement devices, the next part has been obtained as a result of computer simulation, and the last part is computed. Each computed parameter is linear related to some measured parameter. For example, the flow velocity, ?.sub.W:?.sub.W=K G.sub.W, K=1/S, S is the pipe cross-sectional area, where the flowmeter is installed.

(100) Average measured and computed parameters presented in the TABLE 2 are recorded in the baseline file. Further, all n rows in the baseline file are sorted in ascending order of the average values of the flow rate G.sub.Wj: G.sub.W1<G.sub.W2< . . . <G.sub.Wn. The analogical chain of inequalities is valid for the values of average values of the flow velocity ?.sub.Wj:?.sub.w1<?.sub.w2< . . . <?.sub.wn. Besides, the interval between two consecutive average values of the flow rate, interacting with pressure waves, is the same, and is equal to ?G.sub.W: G.sub.W2=G.sub.W1+?G.sub.W; G.sub.W3=G.sub.W2+?G.sub.W=G.sub.W1+2?G.sub.W, . . . , G.sub.Wp=G.sub.W1+p*?G.sub.W. p=1, 2 . . . n. If G.sub.W1=0, then, G.sub.Wp=p*?G.sub.W. p=0, 1, 2 . . . n.

(101) The proposed method of automatic control of energy efficient mode uses only average values of measured parameters over an averaging time interval [m?t]. Sources and magnitude of deviations of measured values from the average, will not be analyzed and considered.

(102) Considering the action of pressure waves on the turbulent flow, a baseline file is created by using programming device 35, which is uploaded and stored in the memory 34 of the dynamic controller 17. The program, developed for implementation of the method of automatic controlling of the energy efficient mode is installed in the microprocessor 33. The program processes the received data from measurement devices, and can use any parameter of the baseline files for calibrating, correcting, and comparing the measured parameters with the parameters of the baseline files.

(103) The data contained in the baseline files is used during the process of automatic controlling of energy efficient mode.

(104) In TABLE 2, more columns can be added. For example: the column comparing the energy efficiency E.sub.fj and ?.sub.fWj during a day, a week, a month; column of the cost of pipeline transportation of a volume of fluid medium, in energy efficient mode and in the normal mode during a day, a week, a month; column (E.sub.fj??.sub.fj) provides an estimated quantity of saved energy.

(105) Experimental studies carried out on a laboratory hydrodynamic loop show that the power component N.sub.e in balance equation (4) is about ?0.1% of the pump power and changes slightly with changes in flow rate in the discharge pipeline.

(106) The data from TABLE 2 is obtained as a result of simulation of the turbulent motion of flow medium in the flow domain of the given sizes and configuration; the optimal frequencies are computed, for which the hydrodynamic resistance coefficient reaches the minimum value for any value of the flow rate belonging to entire practical range.

(107) The main measured parameter is the flow rate, G.sub.Wj, that is equivalent to calculated flow velocity ?.sub.Wj=G.sub.Wj/S and can be generalized to other similar flows using the hydrodynamic similarity criteria. The control function, ?.sub.opt j, is inextricably linked to the level of turbulence of the flow, and is determined by the flow velocity. Current in the PWG motor, T.sub.Wj, its speed, ?.sub.j, power, N.sub.PWj, and the valve opening angle ?.sub.s, are the intrinsic parameters of the wave system. The intrinsic parameters of the pumping system are: current, T.sub.PWj, speed, ?.sub.Wj, and power, N.sub.Wj of the pump motor, pressures, P.sub.SWj and P.sub.DWj, amplitude of pressure waves, P.sub.j, and energy efficiency ?.sub.fWj. All of them depend on the flow rate, diameter, length, and configuration of the pipeline. The intrinsic parameters of the wave system, which determine the control function, reflect the dynamic state of the controlled system, determined by the flow velocity. The intrinsic parameters of the pumping and wave systems are used by the program for continuously providing the energy efficient mode of the process of pipeline transporting of fluid medium and to monitor the stability of the operation of the wave and pumping systems. Therefore, intrinsic parameters of the pumping and wave systems are added as data points to the baseline file considering the turbulent motion of the flow medium under the action of pressure waves.

(108) Simulating the turbulent motion of the flow medium under the action of pressure waves in the flow domain with specific configuration, sizes, and the experimental testing and specifying the obtained results, represent the methodological steps of practical realization of the proposed method. The dimensions of the flow domain are arbitrary. They do not impose any restrictions on the proposed invention.

(109) The obtained results are scalable and can be generalized to pumping and wave systems of other sizes and powers using hydrodynamic, energy, and geometric similarity criteria. Usually, sizes of the flow domain are chosen in a practical range of pipe diameters 6, FIG. 1.

(110) TABLE-US-00002 TABLE 2 Baseline data. Pressure waves action. j G.sub.Wj ?.sub.Wj Re.sub.Wj ?.sub.optj ?.sub.Wj ?.sub.j N.sub.PWj P.sub.SWj P.sub.DWj 1 G.sub.W1 ?.sub.W1 Re.sub.W1 ?.sub.opt1 ?.sub.W1 ?.sub.1 N.sub.W1 P.sub.SW1 P.sub.DW1 2 G.sub.W2 ?.sub.W2 Re.sub.W2 ?.sub.opt2 ?.sub.W2 ?.sub.2 N.sub.W2 P.sub.SW2 P.sub.DW2 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . q G.sub.Wq ?.sub.Wq Re.sub.Wq ?.sub.optq ?.sub.Wq ?.sub.q N.sub.Wq P.sub.SWq P.sub.DWq . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . n G.sub.Wn ?.sub.Wn Re.sub.Wn ?.sub.optn ?.sub.Wn ?.sub.n N.sub.Wn P.sub.SWn P.sub.DWn ? ?.sub.PWj ?.sub.Wj N.sub.Wj ?.sub.Wj V.sub.Wj E.sub.fj ?.sub.fWj E.sub.fj-?.sub.fWj ? ?.sub.PW1 ?.sub.W1 N.sub.W1 ?.sub.W1 V.sub.W1 E.sub.f1 ?.sub.fW1 E.sub.fj-?.sub.fWj ? ?.sub.PW2 ?.sub.W1 N.sub.W2 ?.sub.W2 V.sub.W2 E.sub.f2 ?.sub.fW2 E.sub.f2-?.sub.fW2 . . . . . . . . . . . . . . ? ?.sub.PWq ?.sub.Wq N.sub.Wq ?.sub.Wq V.sub.Wq E.sub.fq ?.sub.fWq E.sub.fq-?.sub.fWq . . . . . . . . . . . . . . . . ? ?.sub.PWn ?.sub.Wn N.sub.Wn ?.sub.Wn V.sub.Wn E.sub.fn ?.sub.fWn E.sub.fn-?.sub.fWn

(111) The main results of the simulation are the computed optimal values of frequencies and amplitudes of pressure waves, providing minimum value of the hydrodynamic resistance coefficient for any value of the flow velocity, belonging to the entire practical range. They are generalized to all fluid media with the same, or very close physical properties: density, dynamic viscosity, compressibility. Therefore, when using dynamic controllers in industrial pumping stations to provide an energy efficient mode of the process of pipeline transporting of fluid media with the above-mentioned physical parameters, these results are uploaded directly into the memory of the dynamic controllers. And intrinsic parameters of the pumping and wave systems are generated during their operation, then are added in the baseline file during the compiling. They can be taken as estimated reference data for calculation of the energy efficiency of real pipeline transporting process.

(112) Therefore, main results can be used for pipelines of any configuration. Wherein, if the initial wave amplitude is not optimal, it can be non-essential corrected by the opening angle ? of the control valve 13, to provide the energy efficient mode. The program used for the automatic controlling in real-time of the energy efficient mode of the process of pipeline transporting of fluid medium, provides a continuously optimal value of the control function, by processing the feedback data of the dynamic state of the controlled system, which is determined by the flow velocity. The control function is maintained in the optimal mode, ?.sub.opt providing minimum value of the hydrodynamic resistance coefficient, ?.sub.minWj (?.sub.optj), and maximum energy efficiency, ?.sub.fWj for any values of the current flow velocity, ?.sub.Wj in the discharge pipeline.

(113) In the process of automatic control of the energy efficient mode, only a part of the measured parameters is used, presented in Table 2 and FIG. 5. For example, the speed of the pressure wave generator motor, which determines the frequency of the generated pressure waves, is determined uniquely, both by using a tachometer 27 and by using an amperemeter or wattmeter. In the proposed method, the current in the PWG motor, ?w is used in the process of automatic control of the energy efficient mode. Similarly, the current in the pump motor, T.sub.PW is used in this process. However, the use of the currents T.sub.W and T.sub.PW for this purpose does not in any way limit the use of the proposed method.

(114) In the engineering system for pipeline transportation of fluid medium, FIG. 1, a method of automatic controlling of energy efficient mode of pipeline transporting process of fluid medium is carried out as follows: The inlet section of the suction pipeline 5, FIG. 1 is immersed in a reservoir filled with fluid medium, in FIG. 1 it is not shown. The pump 3 and the system of measurement devices, FIG. 5 are turned on. Pumping system 1 is operated and the flow medium is transported through the suction pipeline 5, pump 2, and the discharge pipeline 6. A stationary pressure P.sub.S and a flow velocity ?.sub.S are established on the inlet cross section of the suction pipeline S.sub.i. The pressure P.sub.D and the flow velocity ?.sub.D are established on the outlet cross section of the discharge pipeline S.sub.O. The pump pressures P.sub.sp and P.sub.dp are established accordingly, on the inlet and on the outlet of the pump, providing the flow rate G, FIG. 5. For multistage centrifugal pumps, S.sub.S=S.sub.D, therefore ?.sub.S=?.sub.D.

(115) Pump operation mode is carried out by a dynamic controller 17. For example, a constant speed of the pump impeller or a constant current in the pump motor.

(116) The PWG motor 9 is turned on and the pressure waves generator 8 starts to work, generating pressure waves. Speed of the motor is set to default which, usually, may not be optimal. The generated compression pressure wave propagates from the PWG through the connecting pipe 12 in the suction pipeline 5, FIG. 1. The rarefaction pressure wave propagates from the PWG through the connecting pipe 11 in the discharge pipeline 6. The interaction of the rarefaction pressure wave with the flow in the discharge pipeline 6, changes the structure of turbulence and boundary layer 36, FIG. 10. The hydrodynamic resistance coefficient, ?.sub.W decreases. As follows from the head-flow curve of the system, presented in FIG. 9, the flow rate is increased from G.sub.? to G.sub.?w.

(117) Measurement devices 25-32 are measuring instantaneous values of the: speed of the pump motor, by tachometer 25; current of the pump motor, by amperemeter 26; speed of the PWG motor, by tachometer 27; pressure on the inlet of the suction pipeline, by pressure meter 28; flow rate of the fluid medium through the pumping system, by flowmeter 29; discharge pressure on the beginning of the discharge pipeline, by pressure meter 30; temperature of the fluid medium, by thermometer 31; current of the motor of the PWG, by amperemeter 32.

(118) Speed of the pump motor and of the PWG motor can be calculated by the standard voltage supplied to each motor and the measured current.

(119) Values of measured parameters are determined by the internal dynamic state of the turbulent flow. Weak changes in the values of the turbulent flow parameters, the trend, occur in the vicinity of the average value of the measured quantity, for example, of the flow rate, G.sub.W. More significant changes of the flow rate (for example, by the operator), lead to the change of the average value. Thus, the turbulent flow as a controlled system is constantly in a dynamic state. Instantaneous values of all measured parameters are transmitted to the input module 18, FIG. 8. Then, the input module 18 converts the received data into digital format and transmits it to the microprocessor 33. The microprocessor 33 records all received current values of measured parameters into the database file, for example, Excel spreadsheet, in accordance with the sequence presented in TABLE 2. The program is built in such a way, that every m rows of the recorded current instantaneous values of the measured parameters are averaged. The result of averaging is recorded in the row (m+1). Each row with average values of parameters is also recorded in a separate file, for example, Excel spreadsheet. Let us name this file the file of average parameters. The arrangement of parameters in the columns of this file corresponds to TABLE 2. Rows with average values of parameters are written sequentially: M=1, 2 . . . M.sub.T. This is the basic file, based on which, the program performs all steps of the method of automatic controlling of the energy efficient mode of the process of pipeline transporting of fluid medium. To provide the specified method, the program links with all types of files: baseline file, database file, file of average parameters, and uses any needed parameters during the computing process. The program can also select and display on the screen of the dynamic controller any parameter from each of the above-mentioned types of files.

(120) Note. Below we propose one of the schemes of the computational process to ensure an energy efficient mode, using the averaging time interval [m ?t]. However, for this purpose, one can use the other averaging time intervals that are multiples of [m ?t], for example, several averaging time intervals [m?t]. In addition, during the computational process you can use data from any cell of the last recorded row of average values of measured parameters in the file of average parameters, regardless of the number of averaging time intervals used.

(121) Let us consider, how the automatic control of an energy-efficient mode of the process of transporting a flow medium during the transition of the flow from one stationary state to another stationary state, with an insignificant change in the flow rate is performed.

(122) Further, in the text all average measured parameters are accompanied by the subscript wx.

(123) Let us assume that the process of pipeline transporting of fluid medium is carried out in the optimal mode. At the same time, the current values of the average parameters correspond to the last updated row (M?1) in the file of average parameters: G.sub.WXM-1, ?.sub.WXM?1, Re.sub.WXM?1, ?.sub.optWXM-1 . . . ?.sub.f WXM?1.

(124) In the file of average parameters, after an averaging time interval [m?t], the following set of average parameter values is recorded, in accordance with the row M: G.sub.WXM, ?.sub.WXM, Re.sub.WXM, ?.sub.opt WXM . . . ?.sub.fWXM.

(125) The search for the control function, ?.sub.optWXM, corresponding to current state of the controlled system, i.e. the turbulent flow in the discharge pipeline with the flow state parameter G.sub.WXM is a dynamic process. This dynamic process, performed by the program, is started as soon as row M is recorded in the file of average parameters. The implementation of this search is performed successively by the steps.

(126) Let the measured average value of the flow rate be G.sub.WXM. Dividing G.sub.WXM by the interval ?G.sub.W we obtain:

(127) $\begin{array}{cc}\frac{{\stackrel{?}{G}}_{wxM}}{?{\stackrel{?}{G}}_w}=p_M,p_M=\left[p_M\right]+r_M& \left(13\right)\end{array}$ p.sub.M is the quotient; [p.sub.M] is the integer part of the quotient, [p.sub.M]=0, 1, 2 . . . n, r.sub.M is the remainder, a fractional part of the quotient, 0?r.sub.M<1.

(128) Comparison of the remainder r.sub.M and choice of the control function.

(129) Let us assume that [p.sub.M]=j. Then, depending on the value of r.sub.M, for any j, the average value of the flow rate G.sub.WXM can be located at different distances relatively to G.sub.Wj. For example, black circles in FIG. 6. The conditions proposed below determine all mentioned possible locations of the average values of the flow rate G.sub.WXM and are represented schematically in FIG. 6.
1. r.sub.M?0.5 and |G.sub.Wj?G.sub.WXM|<?G.sub.W, then G.sub.WXM is approximated by the average value G.sub.Wj;
2. 0.5?r.sub.M<1 and G.sub.WXM<G.sub.Wj??G.sub.W, then G.sub.WXM is approximated by the average value G.sub.Wj?1;
3. 0.5?r.sub.M<1 and G.sub.WXM>G.sub.Wj+?G.sub.W, then G.sub.WXM is approximated by the average value G.sub.Wj+1.(14)

(130) Conditions similar to (14) are satisfied by the average values of the measured current in the pump motor, current in the motor of pressure waves generator, and of the calculated energy efficiency. They are obtained by replacing the average flow rate in (14) with the average current in the pump motor, T.sub.WXM, the average current in the motor of the pressure waves generator, T.sub.PWXM, which determines the frequency of the generated pressure waves, ?.sub.opt WXM, and the average value of energy efficiency, ?.sub.fWXM. A schematic representation of the average values of these updated parameters, taken from the file of average parameters, is shown on the corresponding axes on the in FIG. 11. a) possible location of the average values of the current of the pump motor, T.sub.WXM, relatively to average value T.sub.Wj, marked as x; b) possible location of the average values of the current of the motor of PWG, T.sub.PWXM, relatively to average value T.sub.PWj, marked as *; c) possible location of the average values of the energy efficiency, ?.sub.fWxM relatively to average value E.sub.fWj marked as +.

(131) The average values of the current of the pump motor, T.sub.PW and of the energy efficiency, ?.sub.fW directly depend on the average value of the flow rate, G.sub.W. Therefore, insignificant changes in flow rate affect proportionally the current in the pump motor and the energy efficiency, as also on the intervals ?ID.sub.W and ?E.sub.fW. However, the average value of the current in the motor of pressure waves generator, T.sub.PW and value of the interval ?I.sub.PW depends on the average value of the flow rate, while the optimal frequencies of pressure waves weakly depend on average value of the flow velocity, and therefore, on the average value of the flow rate.

(132) Practical execution of all steps in a logical sequence, which implementing a method of automatically controlling of energy efficient mode of the process of pipeline transporting of fluid medium, with any changes in flow rate, is described in detail below, flowchart of which is presented in FIG. 12. Separate numbering and naming of blocks in FIG. 12 is missing, since each block describes the performed action in sufficient detail. The combined block 37 can be considered as a separate block. The baseline file can be directly uploaded into the memory of the dynamic controller 17. In this case, the dynamic controller can be used in automatic control systems of energy efficient mode of the process of pipeline transportation of fluid media with the same viscosity. The remaining blocks presented in the flowchart on the FIG. 12, together with the baseline file, are consequently executed by the program, to provide the energy efficient mode. The combined block 37 can be built in the separate computing application and can be used to create the baseline file comprising results of simulation of the interaction of turbulent flow with pressure waves in a pipe, providing minimum hydrodynamic resistance for optimal values of frequencies and amplitudes of pressure waves.

(133) The word parameters used in the blocks in FIG. 12 covers measured parameters, G.sub.WJ, I.sub.Wj, I.sub.PWj, P.sub.SWj, P.sub.DWj and calculated parameters, each of which, uses at least one measured parameter. For example, velocity of the flow, ?.sub.Wj=G.sub.Wj/S. In the process of computing certain steps, in separate blocks, the only specific, specified parameters are used. The valve opening angle, ? is not measured or calculated parameter. The amplitude of generated pressure waves is corrected by changing ?. All Measured and calculated parameters are recorded in the Table 2. Each row of averaged parameter values, over an averaging time interval, is also recorded in the file of average parameters.

(134) The software implementation of the process of automatic control of energy efficient mode is shortly described below.

(135) Let us consider, how the automatic control of the energy efficient mode of the process of transporting a flow medium during the transition from one stationary state to another stationary state, with an insignificant change in the flow rate is performed.

(136) Let us use equality (13).

(137) 1. If r.sub.M=0, then p.sub.M=[p.sub.M], the program automatically assigns to the subscript j the value j=[p.sub.M], activates the row, j=[p.sub.M] from the baseline file for the turbulent motion of the fluid medium in a discharge pipeline, interacting with pressure waves. When G.sub.WXM belongs to the interval G.sub.WXM ?[G.sub.W[pM]??G.sub.W, G.sub.W[pM]+?G.sub.W], it is approximated by the value G.sub.W[pM], in accordance with the 1.sup.st condition of (14). Further, the average value of the current T.sub.W[pM] from the corresponding cell is transmitted in the digital format to the output module 19, where it is converted into analog signal and is transmitted to the automatic output controlled-device 20, which changes the optimal value of operating current of the motor 9 at the current value received from the output module 19, corresponding to optimal value of the frequency ?.sub.opt W[pM], at which the minimum value of the hydrodynamic resistance coefficient is achieved, {dot over (?)}.sub.WpXM=?.sub.WpXM min (?.sub.opt W [pM]). If ?.sub.WpXM satisfies the inequality
|?.sub.WpXM??.sub.W[pM]|<?.sub.?W,(15) then the wave system 7 operates at these optimal parameters values during the averaging time interval [m?t]. ?.sub.?W is the maximum permissible value of the difference of hydrodynamic resistance coefficients. Inequality (8) imposes the following condition on energy efficiency:
|?.sub.fWPXM?(?.sub.fW[PM].sup.in|<?.sub.EX,(16) ?.sub.EX is the maximum permissible value of the difference of the energy efficiencies. Physical conditions (15) provide practical conditions (16). The control of fulfillment conditions of the energy efficient mode, is performed by last block of the flowchart in FIG. 12, b. This is carried out by the variation of the valve opening angle ?, proposed in the flowchart of FIG. 13.

(138) If ?.sub.WpXM differs significantly from ?.sub.W[pM], i.e. (?.sub.WpXM??.sub.W[pM])>?.sub.?W, and therefore, inequality (16) also changes sign, then the amplitude of generated pressure waves is corrected by changing the valve opening angle N of the control valve 13, until the condition (16) is satisfied. For example, the valve opening angle is discretely increasing relatively to a default value, ?.sub.D with the constant step ??. If at the same time the ?.sub.WpXM is increasing, then the direction of the valve opening angle changes in reverse; the valve opening angle ? is decreasing. Moreover, after each step of decreasing the opening angle, ?.sub.s=?.sub.0?s*(??), the energy efficiency, ?.sub.fWPXM is compared using the inequality (16); s=1, 2 . . . S, S is the number of the discrete steps of the change of the opening angle ?; ?.sub.0 is the default value of the valve opening angle. As soon as the condition (16) is satisfied, the valve opening angle is no longer changed and the wave system 7 is working during the averaging time interval [m ?t] in the optimal mode, with optimal parameters of the row j: j=[p.sub.M] of the baseline file, providing the energy efficient mode.

(139) 1.1. If 0<r.sub.M?0.5, and G.sub.WXM belongs to the interval G.sub.WXM?[G.sub.W[pM]??G.sub.W, G.sub.W[pM]+?G.sub.W], then j=[p.sub.M] and G.sub.WXM is approximated by G.sub.W[pM], in accordance with the 1.sup.st condition of (14). The row j=[p.sub.M] in the baseline file remains activated. The current T.sub.PW[pM] in the PWG motor 9 remains the same and, consequently, the optimal frequency ?.sub.optW[pM] of the generated waves also does not change. Thus, when the average values of the flow rate are located within the interval ??G.sub.W relatively to the average value G.sub.W[pM] (weak trend of the stationary flow), the control function ?.sub.opt w [.sub.p M] is in a dynamic state, slightly changing in the vicinity of the optimal value, ?.sub.optW[pM], the 1.sup.st condition of (14), but being applied to optimal frequency.

(140) The state of the controlled system with the flow rate G.sub.W[pM] whose control function is determined by ?.sub.optW[pM] is quasi-stationary, or formally stationary, during the average time interval [m?t]; while the energy efficiency reaches the minimum value, ?.sub.fW[PM].sup.min.

(141) 2. If 0.5<r.sub.M<1, and the average current value of the flow rate satisfies the inequality, G.sub.W[pM]+r.sub.M<G.sub.W[pM]??G w, then G.sub.W[pM]+rM is approximated by the 2.sup.nd condition of (14), G.sub.W[pM]?1. The program automatically assigns to the subscript j the value: j=[p.sub.M]?1 and activates the row j=[p.sub.M]?1 from the baseline file. Further, the average value of the current T.sub.PW[p1]?1 from the corresponding cell is transmitted to the output module 19, then, is converted into analog signal, and is transmitted to the output-controlled device 20. It sets new optimal value of the current T.sub.PW[p1]?1 in the PWG motor 9. Speed of the PWG motor is adjusted to the new value of the current, T.sub.PW[pM]?1, and therefore the optimal frequency of the generated pressure waves is adjusted to its new value ?.sub.opt W[pM]?1, at which the minimum value of the energy efficiency is reached, ?.sub.fWxM=?.sub.fW.sup.min. If E.sub.fW [pM]?1 satisfies the inequality (16), then the wave system 7 operates at these optimal parameter values during the averaging time interval [m ?t]. If inequality (16) is not satisfied, the amplitude of the generated pressure waves is corrected by varying the valve opening angle ? of the control valve 13, until the condition (16) is satisfied. Performing of this procedure is described above.

(142) 3. If 0.5<r.sub.M<1, and the average current value of the flow rate satisfies the inequality, G.sub.W[pM]+rM>G.sub.W[pM]+?G.sub.W, then G.sub.W[pM]+rM is approximated by the 3.sup.rd condition of (14), G.sub.W[pM]+1. The program automatically assigns to the subscript j the value: j=[pM]+1 and activates the row j=[p.sub.M]+1 from the baseline file. Further, the average value of the current T.sub.PW[pM+1 from the corresponding cell is transmitted to the output module 19, then it is transmitted to the output-controlled device 20. It sets new optimal value of the current T.sub.PW[pM]+1 in the PWG motor 9. The PWG motor speed is adjusted to the new current values, T.sub.PW[pM]+1, and therefore the optimal frequency of the generated pressure waves is changed to its new value ?.sub.opt W[pM]+1, at which the minimum value of the energy efficiency is reached, ?.sub.fWxM=?.sub.fWXM.sup.min. If ?.sub.fW[pM]+1 satisfies the inequality (?.sub.fWxM??.sub.fW [pM]+1)<?.sub.EX, then the wave system 7 operates at these optimal parameters values during the average time interval [m ?t]. If inequality (16) is not satisfied, the wave amplitude is corrected by varying the valve opening angle ? of the control valve 13 by controlled-device 21, until condition (16) is satisfied. Performing of this procedure is described above.

(143) Such insignificant adjustment of the values of control function between the average value ?.sub.opt W[pM] and neighboring values ?.sub.opt W[pM]?1 and ?.sub.opt W[pM]+1 is usually realized in a quasi-stationary flow, when the flow rate smoothly oscillates (trend) between the neighboring values of G.sub.W[pM]?1 and G.sub.W[pM]+1 relatively to average value G.sub.W[pM].

(144) Similar relationships are obtained for average values of flow velocities, when replacing the flow rate G to ?*S in equality (13). Therefore, the proposed analysis in terms of the average values of the flow rate is also valid in terms of the average values of the flow velocity.

(145) Let us assume, that the current values of the flow rate in the discharge pipeline equal to G.sub.WXA and is carried out in an energy efficient mode. Suppose, that as a result of technological needs, the flow rate is changed. The new value of the steady flow becomes G.sub.WXB. The change in flow rate from G.sub.WXA to G.sub.WXB is a transitional process. It is carried out in the time interval ?t.sub.x, called a transition period, which depends on the value of the difference G.sub.WXB?G.sub.WXA and its velocity. Moreover, a period of the transient process is much more than the averaging time interval: ?t.sub.x>>[m?t]. The program automatically performs all actions, in order to provide the energy efficient mode, using only parameters received from the measurement devices. The average values of the measured flow rate, during the transition period ?t.sub.x, change discretely, and after each averaging time interval [m?t] a new set of transitional values of the flow rate is generated: {G.sub.WXh}, h=A, A+1, A+2 . . . B?1, B; or, {G.sub.WXA, G.sub.WXA+1, G.sub.WXA+2 . . . G.sub.WXB?1, G.sub.WXB}.

(146) Assume, that the first measured value of the transient flow rate is G.sub.WXA+1. Dividing by the interval ?G.sub.W we obtain:

(147) $\begin{array}{cc}\frac{{\stackrel{?}{G}}_{WXA+1}}{?{\stackrel{?}{G}}_W}=p_{A+1},p_{A+1}=\left[p_{A+1}\right]+r_{A+1},h=A+1.& \left(17\right)\end{array}$ The structure of the division result (17) is similar to (13) but in different notations. p.sub.A+1 is the quotient, p.sub.A+1=[p.sub.A+1], 0<r.sub.A+1<1. Therefore, the steps of searching and selecting the value of the control function, ?.sub.opt WXh, and transmitting the value of the optimal current to the output module 19 of the dynamic controller 17, are similar with the described above comparisons 1-4 of the fractional part of the quotient, 0?r.sub.h<1.

(148) If r.sub.A+1=0, then p.sub.A+1=[p.sub.A+1], and the subscript j is assigned a value [p.sub.A+1]: j=[p.sub.A+1]. G.sub.WXA+1?[G.sub.W[pA+1]??G.sub.W, G.sub.W[pA+1]+?G.sub.W] in accordance with the 1.sup.st condition of (14).

(149) If 0?r.sub.A+1?0.5 and G.sub.WXA+1?[G.sub.W[pA+1]??G.sub.W, G.sub.W[pA+1]+?G.sub.W] then the subscript j is assigned a value [p.sub.A+1]: j=[p.sub.A+1]. For both cases, r.sub.A+1=0 and 0?r.sub.A+1?0.5, G.sub.WXA+1 is approximated by G.sub.W[pA+1].

(150) The program automatically activates the row j=[p.sub.A+1] from the baseline file, and performs sequentially the following steps: 1. transmitting value of the current T.sub.W[pA+1] from the corresponding cell to the output module 19; 2. transmitting value of the current T.sub.W[pA+1] from the output module 19 to the output-controlled device 20; 3. re-setting the optimal value of the current of the PWG motor 9; 4. providing the generating pressure waves of optimal frequency, ?.sub.optW[pA+1]; 5. reaching the minimum value of the hydrodynamic resistance coefficient, {circumflex over (?)}.sub.W[pA+1]min; 6. checking the condition (16); 7. correcting the amplitude of pressure waves by varying the valve opening angle ?, of the control valve, for example, 13, by controlled-device 21, if necessarily. 8. providing the energy efficient mode of the process pipeline transporting of fluid medium.

(151) Detailed of performing steps 1-8 is described above.

(152) If 0.5?r.sub.A+1<1, and the average current value of the flow rate satisfies the inequality, G.sub.WX[pA+1]+r.sub.A1+1<G.sub.W[pA+1]??G.sub.W, then G.sub.WX[pA+1]+r.sub.A+1 is approximated by G.sub.W[pA+1]?1, the 2.sup.nd condition of (14). The program automatically assigns to the subscript j the value: j=[p.sub.A+1]?1, activates the row j=[p.sub.A+1]?1 from the baseline file, and performs sequentially the above steps 1-8.

(153) If 0.5<r.sub.A+1<1, and the average current value of the flow rate satisfies the inequality, G.sub.WX[pA+1]+rA+1>G.sub.W[pA+1]+?G.sub.W, then G.sub.WX[pA+1]+rA+1 is approximated by the 3.sup.rd condition of (14), G.sub.W[pA+1]+1. The program automatically assigns to the subscript j the value: j=[p.sub.A+1]+1, activates the row j=[p.sub.A+1]+1 from the baseline file, and performs sequentially the above steps 1-8.

(154) Similar computational steps are performed automatically by the program for the entire set of transitional values of flow rate: {G.sub.WXA, G.sub.WXA+1, G.sub.WXA+2 . . . G.sub.WXB-1, G.sub.WXB}. When the transitional values of flow rate G.sub.WXAh get closer to the next stationary value G.sub.WXB, h=A, A+1 . . . B?1, B, the difference between two successive values G.sub.WX[Ah]?G.sub.WX[Ah]?1 becomes quite small. The flow state G.sub.WXh at h=B is stationary. At the same time, the computational steps for choosing the optimal frequencies, ?.sub.opt W [Ah] are carried out sequentially, as described above.

(155) Providing the energy efficient mode is performed using a program, that analyzes, compares, computes, and continuously adjusts the value of the control function, ?.sub.opt, to the current state of the controlled system.

(156) To compare the energy efficiency of the action of pressure waves on the process of pipeline transporting of fluid medium, formula (1) is used. Comparative results calculated using the intrinsic parameters of the pumping system (the action of the constant pressure on a turbulent flow) and of the intrinsic parameters of the pumping and wave systems (the action of pressure waves on a turbulent flow), presented in the TABLE 2.

(157) Thus, in the proposed engineering system, the method of automatic controlling of the energy efficient mode of the process of pipeline transporting of fluid medium is provided in real time in both, stationary and non-stationary turbulent flows. In other words, the proposed energy efficient mode provides in real-time the optimal value of the control function for both the stationary state of the controlled system and for any transient states.

(158) The proposed engineering system for pipeline transportation of fluid medium and method of automatic controlling of the energy efficient mode can be successfully used for pipeline transportation of technological media, brines, different heterogeneous media, whose solid particles are sedimented on the inner surface of a pipe, and form a solid phase layer. The formation of such a layer leads to the decrease in the cross-sectional area of the flow, an increase in power of the pump motor, and a decrease in flow rate. The energy costs of operating such pipelines become economically unprofitable.

(159) The use of the system of generating pressure waves 15 in relation to the transportation of such media has the following advantage: reducing the rate of sedimentation of solid particles on the inner surface of the pipe, due to longitudinal sign alternating acceleration of the elementary fluid particles in the field of the pressure wave; increasing the life of pipeline.

(160) Prevention of decrease in sedimentation of solid particles is ensured not only by the frequency of pressure waves, but also by their form and amplitude. Even a small increase in the service life of pipes in industrial conditions is very important from economic standpoint. For example, in salt production factories, certain technological processes are accompanied by pumping brines at a temperature of 203? F. After about six months of operation, the cross-sectional area of the flow in this pipeline is significantly reduced. The specific energy consumption is increased and pipeline operation becomes unprofitable. These pipes are replaced with new ones.

(161) A baseline file similar to the one described in the method is used in control systems for certain industries. For control systems in such industries, the baseline file is similar to the one presented in the described method. A separate column in the baseline file must be added, to visualize the dynamics of specific energy consumption on a daily or even hourly basis. As soon as this parameter reaches the upper limit, determined by the energy efficiency, the dynamic controller sends an alert signal to the operator. The pumping system can be stopped automatically or by the operator.

(162) The modern filtration systems, for example, wastewater filtration systems, are operating only under the action of the constant pressure. The main problem in the modern wastewater filtration technology is the relatively quick filter clogging and filter cake forming. Consequently, the filter hydrodynamic resistance is increased. The pump power, and specific energy consumption to maintain the same flow rate, are also increased. Frequent filters replacements are relatively expensive. Therefore, at least a slight extension of the filter life, with acceptable energy efficiency, is economically profitable.

(163) The proposed engineering system and method of automatic controlling of the energy efficient mode can be successfully used in the wastewater filtration technology. In such case, in engineering system FIG. 2, a filter is installed behind the conjugated surface S.sub.J downstream, FIG. 9 (filter is not shown).

(164) Longitudinal sign-alternative wave action on the filter cake, maintains it in unstable, loosened state; the state with more porosity. This state increases the performance of the filter and filtration efficiency. In addition, periodical action of pressure to the impurities within the filter itself, prevents their coagulation and settling on the pore walls, or solid filter particles, i.e., slowing down the filter clogging. Both these pressure wave actions lead to an increase in the life of the filter.

(165) For some filtration systems the PWG can be connected parallel to the filter.

(166) For filtration processes, the wave system 7 must generate pressure waves of such form and amplitude as to minimize the hydrodynamic resistance of the current state of the controlled system, new filter, and the filter cake. The optimal form, amplitude, and frequency of pressure waves for a given filter material and types of wastewater impurities is determined as a result of modeling and further experimental correction. This data forms the basis of the baseline file. Essentially, it is a table function ?.sub.opt fj(?, P.sup.?(t), P.sub.fj), where P.sub.fj is the average value of the inlet pressure on the filter; this is the main technological parameter of the filtering process; ?.sub.opt fj is the optimal value of the control function of the filter system: filter and filter cake; P.sup.?(t) is the law of change of pressure versus time in the generated wave. As applied to the filtration process, the TABLE 2 is modified. Some parameters are removed, for example, Reynolds number, average flow velocity. In the second column, the flow rate, G.sub.Wj is replaced to the main technological parameter of the filtering process, P.sub.fj. Further, all n rows in the baseline file are sorted in ascending order of the average values of the inlet pressure on the filter P.sub.fj:P.sub.f1<P.sub.f2< . . . <P.sub.fn. Besides, the interval between two consecutive average values of the inlet pressure is the same and is equal to ?P.sub.fW. The energy efficiency of the filtration process is determined by formula (1), where the denominator is the volume of filtered liquid over the time interval ??.sub.G.

(167) It will be understood that the proposed engineering system with such dynamic controller, with an installed program and baseline file of technological parameters of the microprocessor and memory is provided by way of example but not limitation, and embodiment of the present invention can be used in engineering systems for pipeline transporting of fluid media with more complicated properties, and in other technological processes with changeable parameters, for the purpose of automatic controlling of energy efficient mode. Moreover, the control function does not necessarily have to be of the optimal value of frequency of generated pressure waves. For some technological applications, the form of pressure wave or amplitude may be optimal. Ultimately, the amplitude-frequency characteristics of the wave are important. It is the dynamics of pressure changes in the wave that determines the energy efficiency of the wave's impact on the controlled system. Performing such action on the controlled systems by pressure waves of various forms and amplitudes is ensured by the frequency of the generated waves, which are automatically controlled by the speed of pressure waves generator.