ISSN ONLINE(2278-8875) PRINT (2320-3765)
N.Visali1, S.Niranjan2
|
Related article at Pubmed, Scholar Google |
Visit for more related articles at International Journal of Advanced Research in Electrical, Electronics and Instrumentation Engineering
Electrification of ship propulsion system with hybrid power systems created interest in modeling of all marine vessels. Reduction in the emission of green house gases, fuel consumption and efficient energy management with improved system dynamic performance is possible with hybrid power systems. With penetration of power electronic converters into power systems, ship board direct current distribution systems offer further more advantages such as space, weight saving and flexible arrangement of equipment. On board marine vessel consists of synchronous generator, energy storage system, fuel cell, photo voltaic array and induction motor drive. Power electronic converters are modeled through detailed model. Mechanical components are modeled with mathematical analysis approach. Power sharing among the proposed hybrid power system is observed in MATLAB/Simulink software.
Keywords |
||||||||||||||||||||||||||||||||||||||||||||||||||||||
Simulation, DC distribution systems, detailed modeling, PVA module, and vector (or) field oriented control. | ||||||||||||||||||||||||||||||||||||||||||||||||||||||
I. INTRODUCTION |
||||||||||||||||||||||||||||||||||||||||||||||||||||||
Electrical installations are present in every marine vessel ranging from powering of communication and monitoring systems, running of different motors to high power electrical installation for electric propulsion [1]-[2]. With the possibility to control electrical motors with variable speed drives in a large power range with compact, reliable and cost competitive solutions, the use of electrical propulsion has emerged in new application areas. Electrification of propulsion system brought us advantages like improved life cycle cost by reduced fuel consumption and maintenance, less propulsion noise and vibrations, Less space consuming etc[3]-[5]. | ||||||||||||||||||||||||||||||||||||||||||||||||||||||
Now a day‘s all electrical ships are enabled by hybrid power systems [6]. Hybrid power systems use low calorie renewable energy resources to generate electrical energy [7]-[9]. In this paper hybrid power system incorporates a set of primary and auxiliary energy sources in order to provide the propulsion power and at the same time, to energize the ship board electric loads. A gas or diesel energy source comes under primary sources which are capable to drive both propulsion and also ship loads. Fuel cell, photo voltaic, and energy storage systems comes under auxiliary sources which participate only in improving ship dynamic performance and also ship load system. | ||||||||||||||||||||||||||||||||||||||||||||||||||||||
Models are derived for different electrical and mechanical elements including the synchronous generator–rectifier system, inverter, dc/dc converters, diesel engine, propeller, and ship hydrodynamics. Power electronic converters are modeled in detailed modeling to give an effective solution and to improve speed of simulation. A simulation platform is developed in MATLAB/Simulink for system-level studies of hybrid electric ships. As a case study, marine vessel is simulated in different modes of operation. The simulation results of a power sharing control among two diesel generators, a fuel cell unit, photo voltaic array and an energy storage system show practical utility of such a simulation tool in system studies associated with design, evaluation, power management. | ||||||||||||||||||||||||||||||||||||||||||||||||||||||
II. SYSTEM CHARACTERISTICS AND OVERVIEW |
||||||||||||||||||||||||||||||||||||||||||||||||||||||
The single line diagram of an overall system is shown in the above Fig. 1. Modeling of system-level analysis of ship dc distribution power systems requires specific features imposed by the system characteristics. First of all, there are mechanical and electrical components in such a system with extremely different dynamics. This difference ranges from a very small time constant of a power electronic switch, which is on the order of nanoseconds to a large time constant of ship hydrodynamics, which is tens of seconds. Detailed switching models of power electronic converters are considered hear. Second, for system level studies, in which large variations occur in system states, small signal or linearized averaging is invalid, and large signal models must be developed instead. Third, the interactions among connected elements should be included in order for the system model to be realistic. In other words, each element should receive the reaction of the connected elements, as soon as it affects them. | ||||||||||||||||||||||||||||||||||||||||||||||||||||||
III. MODELING OF POWER ELECTRONIC CONVERTERS |
||||||||||||||||||||||||||||||||||||||||||||||||||||||
It is the power electronics, area that brought a great change in technology of electric power industry. Compared to rotating commutation static conversion is more efficient and reliable. The DC-DC converters are very prominent for energy management in shipboard hybrid power systems [10]. These converters are used to shift voltage from one level to another among the dc distribution system. Here we use two types of DC-DC converters, a boost converter to step up the voltage generated from fuel cell module and a boost type bidirectional converter for battery cell module. | ||||||||||||||||||||||||||||||||||||||||||||||||||||||
A. DC-DC Boost Converter | ||||||||||||||||||||||||||||||||||||||||||||||||||||||
The step up dc to dc converter is basically called as boost converter [10]. Topological circuit diagram of full bridge boost converter is shown in above fig.2. Here we can observe that when controllable switch (Sc) is in âÃâ¬Ãâon‘ state, current flows through inductor and it charges. When switch is in âÃâ¬Ãâoff‘ state the energy stored inductor and supply voltage together gets voltage boost. The amount of boosting voltage depends on duty cycle of controllable switch. | ||||||||||||||||||||||||||||||||||||||||||||||||||||||
B. DC-DC Bi-Directional Converter | ||||||||||||||||||||||||||||||||||||||||||||||||||||||
Most of the existing bidirectional dc-dc converters fall into the generic circuit structure illustrated in the above fig. 3 which is characterized by a current fed or voltage fed on one side. Based on the placement of the auxiliary energy storage, the bidirectional dc-dc converter can be categorized into buck and boost type [12]. The buck type is to have energy storage placed on the high voltage side, and the boost type is to have it placed on the low voltage side [11]. The high frequency transformer based system is an attractive one to obtain isolation between the source and load sides, since transformer can isolate the voltage sources and provide the impedance matching between them. | ||||||||||||||||||||||||||||||||||||||||||||||||||||||
C. Synchronous Machine-Rectifier System | ||||||||||||||||||||||||||||||||||||||||||||||||||||||
Nearly all of the electrical power used throughout the world is generated by the synchronous machines driven either by steam or hydraulic turbine or any combustion engines. The synchronous machine is the principle means of converting mechanical energy to electrical energy [13]. Here we prefer d-q model to eliminate the dependency of inductance on rotor position. Since rotor field, d and q-axis are already aligned in rotor coordinates it is only stator voltages, currents and flux linkages that have to be transformed to the rotor coordinates [14]. This transformation is done through park‘s transformation, valid for voltages, currents, and flux linkages. | ||||||||||||||||||||||||||||||||||||||||||||||||||||||
As a case study synchronous machine rectifier system is modeled in detailed model. The response to a step change in the higher loads (0.13Ω to 0.2Ω) is also shown in the above fig. 4. Similarly for a step change in the lower loads (2.8Ω to 1Ω) is observed here in the above fig. 5. Modeling of the synchronous machine in the detailed model gives us better precision than any other average value models with constant and also variable parameters. | ||||||||||||||||||||||||||||||||||||||||||||||||||||||
D. Modeling of induction motor drive | ||||||||||||||||||||||||||||||||||||||||||||||||||||||
Induction motor drives have been the work horses for variable speed applications ranging from fractional horse power to multi-megawatt. Modeling of induction motor drive is done in vector control method, though scalar control method is simple to implement, due to its inherent coupling effect gives sluggish response and the system is prone to instability. A vector controlled induction motor drive operates like a separately exited dc motor drive [14]. | ||||||||||||||||||||||||||||||||||||||||||||||||||||||
IV. MODELING OF PHOTO VOLTAIC ARRAY |
||||||||||||||||||||||||||||||||||||||||||||||||||||||
Photo voltaic solar cells (PVSC) are emerging as a renewable green energy power source. Use of solar energy gives reduces fuel consumption and emission of green house gases. Combined series/parallel combinations of photo voltaic solar cells, which are usually represented by a simplified equivalent circuit model shown in below fig. 6 or by equation 3. | ||||||||||||||||||||||||||||||||||||||||||||||||||||||
Where βT = 0.004 and γT = 0.06 for the cell used Ta=200c is the ambient temperature during cell testing. This is used to obtain the modified model of cell for another ambient temperature Tx. If solar irradiation (Sc) increases from Sx1 to Sx2 the cell operating temperatures and currents also increases to new values. Thus this change due to solar irradiation can be implemented using correction factors. | ||||||||||||||||||||||||||||||||||||||||||||||||||||||
Where Sc is the original solar irradiation level and Sx is the new level. CSI, CSI are correction factors for cell output voltages and currents and ΔTC is the change in the temperature with αS slope change in the operating temperature. By using all correction factors new output voltage VCX and photo current Iphx are obtained as follows: | ||||||||||||||||||||||||||||||||||||||||||||||||||||||
Here above fig. 7 is the sub model of the photo voltaic array block in above fig. 8. | ||||||||||||||||||||||||||||||||||||||||||||||||||||||
V. MODELING OF MECHANICAL SYSTEMS |
||||||||||||||||||||||||||||||||||||||||||||||||||||||
1. Propeller Modeling | ||||||||||||||||||||||||||||||||||||||||||||||||||||||
Propeller can be modeled by implementing relations between torque, thrust and speed. For fixed propeller these relationships can be expressed as given below eq. 10 and eq.11 [16]. | ||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
||||||||||||||||||||||||||||||||||||||||||||||||||||||
2. Ship Hydro Dynamics | ||||||||||||||||||||||||||||||||||||||||||||||||||||||
The vessel hydrodynamic model could be complicated when dynamic positioning is purposed [17]. For power system analysis in integrated power and propulsion systems, a one dimensional forward-motion model suffices. | ||||||||||||||||||||||||||||||||||||||||||||||||||||||
Therefore, the ship can be treated as an inertial mass with a resistive drag, on which the propeller thrust Tp applies. The resistive drag is proportional to the square of ship speed Vs. | ||||||||||||||||||||||||||||||||||||||||||||||||||||||
3. Diesel Engine | ||||||||||||||||||||||||||||||||||||||||||||||||||||||
Diesel engine dynamics can be modeled in different levels of complexity, depending on the application. In this paper, it is approximated by a time delay τ and a time constant τC [16] as represented in (14). | ||||||||||||||||||||||||||||||||||||||||||||||||||||||
Where Tm is the generated torque, Ky is the torque constant, and Y is the fuel index (governor setting). The time delay is half the period between consecutive cylinder firings, which can be calculated by eq. 15 [37]. | ||||||||||||||||||||||||||||||||||||||||||||||||||||||
Where nm is the engine rotational speed in [rps] and N is the number of cylinders. The time constant is calculated by eq.16.In order to model the interaction between the diesel engine and the synchronous generator, a speed control loop is established between speed governor and synchronous generator. | ||||||||||||||||||||||||||||||||||||||||||||||||||||||
VI. SYSTEM SIMULATION AND RESULTS |
||||||||||||||||||||||||||||||||||||||||||||||||||||||
The below fig. 9 represents overview of the system implemented in MATLAB/Simulink. Interfaces among elements are bidirectional to represent interactions among connected elements. This facilitates the interface among different types of electrical, mechanical, and electromechanical elements. Parameters for different parts of the system are given in the Appendix. To represent the simulation results, a sailing profile including high speed (140 rad/sec), moderate speed (70 rad/se), and low speed (30 rad/sec) is studied in this section with and without solar module. | ||||||||||||||||||||||||||||||||||||||||||||||||||||||
Fig. 10 (a), (b) shows per unit mechanical power and terminal voltage of generator set-1. Fig. 11 (a), (b) shows per unit mechanical power and terminal voltage of and generator set-2. | ||||||||||||||||||||||||||||||||||||||||||||||||||||||
Fig. 12(a) shows the variation in the stator current of motor used for propulsion. Fig. 12(b) represents ship actual speeds and reference speeds of propulsion motor. Fig. 12(c) represents propeller reference torque and electromagnetic torque of induction motor over the simulation time. | ||||||||||||||||||||||||||||||||||||||||||||||||||||||
In the first cycle, the ship is accelerating from zero to full speed, which examines the power system response to high stress and also due to high inertia during start current drawn during this cycle of operation is high. The propeller accelerates much faster than the ship, as dynamics of the propulsion system is much faster than ship motion so during first cycle there is a deviation of the actual speed from reference speed, but in later cycles of operation both speeds are same with small deviations. In reference electromagnetic torque and actual electromagnetic torque of the propulsion induction motor there is minor deviation between these two torques which is due to load variation during first cycle is high i.e. from zero to maximum torque value, so variation between torques is more in first cycle of operation. In remaining two stages of operation better match is obtained between reference torque and actual torque. | ||||||||||||||||||||||||||||||||||||||||||||||||||||||
Fig. 13 to fig.16 shows dc currents of different units that are connected to the dc grid. During high speeds 140 rad/sec, the propulsion system demands highest amount of power. Therefore, both of the diesel generators, the fuel cell, and the energy storage supply power to the grid during the first cycle of simulation. The propulsion power highly reduces with the reduction of propeller rotational speed and ship speed. | ||||||||||||||||||||||||||||||||||||||||||||||||||||||
Reduction of the ship speed by half leads to reduction in the propeller speed. At this speed range 70 rad/sec, one of the generators (G2) shuts down and the power is mainly supplied by the remaining generator (G1) and the fuel cell. At slow speed operation (30 rad/sec), the generator G1 also shuts down, and the power is mainly supplied by the fuel cell, and battery. | ||||||||||||||||||||||||||||||||||||||||||||||||||||||
Fig. 17 shows the dc bus voltage which is maintained constant with small droop. Fig. 18 shows total load current drawn by overall system. The battery state of charge (SOC) is maintained within certain limits. This is done by starting an idle diesel generator when SOC reaches the lower boundary, and by shutting down an active diesel generator when SOC exceeding the upper limit. Fig. 18 shows SOC, voltage, and current of the battery. Fig.20 shows load current of ship auxiliary and hotel load rated at 175hp, which is modeled by a constant dc motor load. | ||||||||||||||||||||||||||||||||||||||||||||||||||||||
B. With solar module under operationThis addition of solar module provides ultra-low emission mode of operation, which is important when ship is approaching/leaving harbors at low speeds. Improvement in voltage profile is observed with addition of additional solar module. Some of important waveforms are shown below. | ||||||||||||||||||||||||||||||||||||||||||||||||||||||
Fig. 22 shows the dc bus voltage which is maintained constant with small droop we can observe the improvement in voltage profile with additional hybrid power system. Fig. 23 to fig.27 shows dc currents of different units that are connected to the dc grid. Since there is no change in load pattern the electromagnetic torque, thrust, auxiliary load current and total load current drawn will remains same as shown in the above case without solar module under operation. | ||||||||||||||||||||||||||||||||||||||||||||||||||||||
The control methods for power sharing and voltage regulation of all-electric ships with dc distribution systems can be tested and optimized by this simulation program. The derived models and the simulation platform can be reconfigured to assist in design, implementation, and control of experimental setups of dc power systems for all-electric ships. | ||||||||||||||||||||||||||||||||||||||||||||||||||||||
VII. CONCLUSION |
||||||||||||||||||||||||||||||||||||||||||||||||||||||
A simulation platform is developed using the derived models of different components presented in the paper. Detailed modeling method is used to model all power electronic converters without neglecting system dynamic behavior. We have observed the exact system behavior under transient conditions. The simulation results for a sailing profile of an all-electric ship and how power sharing is done between different mechanical and electrical components with different dynamics. A significant reduction in emission of green house gases and use of non-renewable energy (diesel) is done by adding a green energy source so that better energy management is done and improvement in voltage profile is obtained. With this green power sources ultra low emission of effluent gases like near shore is obtained. | ||||||||||||||||||||||||||||||||||||||||||||||||||||||
Figures at a glance |
||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
||||||||||||||||||||||||||||||||||||||||||||||||||||||
References |
||||||||||||||||||||||||||||||||||||||||||||||||||||||
|