[go: up one dir, main page]

Next Issue
Volume 13, February-1
Previous Issue
Volume 13, January-1
 
 
energies-logo

Journal Browser

Journal Browser

Energies, Volume 13, Issue 2 (January-2 2020) – 220 articles

Cover Story (view full-size image): To cope with stringent emission rules, modern vessels shift towards hybrid/electric vessels by employing a battery energy storage system as a main part of their onboard power systems. Therefore, this research paper focuses on large-scale integration of renewables and proper sizing and location of battery energy storage systems in order to meet the required power for the vessels. View this paper.
  • Issues are regarded as officially published after their release is announced to the table of contents alert mailing list.
  • You may sign up for e-mail alerts to receive table of contents of newly released issues.
  • PDF is the official format for papers published in both, html and pdf forms. To view the papers in pdf format, click on the "PDF Full-text" link, and use the free Adobe Reader to open them.
Order results
Result details
Section
Select all
Export citation of selected articles as:
10 pages, 1577 KiB  
Article
Growth and Production of Lipids in Raphidocelis subcapitata Immobilized in Sodium Alginate Beads
by Amel Benasla and Robert Hausler
Energies 2020, 13(2), 506; https://doi.org/10.3390/en13020506 - 20 Jan 2020
Cited by 5 | Viewed by 4181
Abstract
The growth and production of lipids in the green microalga Raphidocelis subcapitata immobilized in alginate gel are studied. The beads are made from alginate (2% w/v) and CaCl2 (1% w/v). The dry weight, the concentration of [...] Read more.
The growth and production of lipids in the green microalga Raphidocelis subcapitata immobilized in alginate gel are studied. The beads are made from alginate (2% w/v) and CaCl2 (1% w/v). The dry weight, the concentration of cells, and the lipid content are determined after dissolution of the beads in a sodium phosphate buffer. The results show that variations in biomass do not reflect variations in the number of cells in R. subcapitata. Cells divide more rapidly (Gc = 3.45 ± 0.3 days) than biomass is produced (Gm = 4.1 ± 0.4 days) during the exponential growth phase. Therefore, the average mass of the immobilized cells decreases until it reaches its minimum at the end of the exponential phase. Thus, during the stationary phase, cell division ceases while biomass production continues, and the average mass of the immobilized cells increases. In the present study, it is shown that this increase is due to the accumulation of lipids following the depletion of nitrates and phosphates in the culture medium. A lipid content of 24.7 ± 2.5% (dcw) and a lipid productivity of LP = 29.8 ± 3.0 mg/L/day are recorded at the end of culture. These results suggest that immobilized R. subcapitata has promising potential for biodiesel production. Full article
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Growth of <span class="html-italic">R. subcapitata</span> immobilized in alginate beads over a period of 10 days (culture performed in batch and in duplicate): (<b>a</b>) variations in the number of cells; (<b>b</b>) variations in biomass; (<b>c</b>) variations of the average mass of cells. <span class="html-italic">t</span><sub>0</sub> <span class="html-italic">= d</span><sub>1</sub>.</p>
Full article ">Figure 2
<p>Production of total lipids in <span class="html-italic">R. subcapitata</span> immobilized in alginate beads, over a period of 10 days (culture performed in batch and in duplicate). (<b>a</b>) Nitrate consumption; (<b>b</b>) phosphate consumption; (<b>c</b>) variations in lipid content per dry weight of cells; (<b>d</b>) variations in lipid concentration per cell. <span class="html-italic">t</span><sub>0</sub> <span class="html-italic">= d</span><sub>1</sub>.</p>
Full article ">
16 pages, 4010 KiB  
Article
Control Strategies and Economic Analysis of an LTO Battery Energy Storage System for AGC Ancillary Service
by Bingxiang Sun, Xitian He, Weige Zhang, Yangxi Li, Minming Gong, Yang Yang, Xiaojia Su, Zhenlin Zhu and Wenzhong Gao
Energies 2020, 13(2), 505; https://doi.org/10.3390/en13020505 - 20 Jan 2020
Cited by 7 | Viewed by 3424
Abstract
With the rapid growth of renewable energy and the DC fast charge pile of the electric vehicle, their inherent volatility and randomness increase a power system’s unbalance of instantaneous power. The need for power grid frequency regulation is increasing. The energy storage system [...] Read more.
With the rapid growth of renewable energy and the DC fast charge pile of the electric vehicle, their inherent volatility and randomness increase a power system’s unbalance of instantaneous power. The need for power grid frequency regulation is increasing. The energy storage system (ESS) can be used to assist the thermal power unit so that a better frequency regulation result is obtained without changing the original operating mode of the unit. In this paper, a set of different charging/discharging control strategies of the lithium titanate battery (LTO) is proposed, which are chosen according to the interval of the State of energy (SOE) to improve the utilization rate of the ESS. Finally, the cost-benefit model of the ESS participating in automatic generation control ancillary service is established. Case analysis proves that after a 1.75 MWh ESS is configured for a 600 MW thermal power unit, Kp and D is increased from 1.42 to 6.38 and 2857 to 6895 MW. The net daily income is increased from 20,284 yuan to 199,900 yuan with a repayment period of 93 days. The results show that the control strategies and the energy configuration method can improve the performance and economic return of the system. Full article
(This article belongs to the Special Issue Thermal and Energy Management of Battery-Operated Systems)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>The setting point control process of automatic generation control (AGC) generator.</p>
Full article ">Figure 2
<p>The energy storage system (ESS) participates in AGC ancillary service.</p>
Full article ">Figure 3
<p>The curve of (<b>a</b>) State of Charge (SOC)–open circuit voltage (OCV) and (<b>b</b>) SOC–state of energy (SOE).</p>
Full article ">Figure 4
<p>Strategy of charging and discharging.</p>
Full article ">Figure 5
<p>The strategy of ESS when the thermal power unit is in reverse.</p>
Full article ">Figure 6
<p>The Negative Strategy of ESS.</p>
Full article ">Figure 7
<p>The simulation results of the 3 MWh energy storage system: (<b>a</b>) Power, (<b>b</b>) Storage, and (<b>c</b>) SOE.</p>
Full article ">Figure 8
<p>The generator follows the AGC instruction curve for 24 h.</p>
Full article ">Figure 9
<p>Analysis of output deviation and regulating energy (<b>a</b>) The output deviation between AGC command and units (MW) (<b>b</b>) Energy required for each regulation (kWh).</p>
Full article ">Figure 10
<p>The relationship between <span class="html-italic">K<sub>p</sub></span>, <span class="html-italic">D</span>, and energy of ESS.</p>
Full article ">Figure 11
<p>The relationships between system net profit, repayment period, and energy configuration.</p>
Full article ">Figure 12
<p>Data perspective of the output difference: (<b>a</b>) difference of power without ESS and output power of ESS and (<b>b</b>) difference of power without ESS and difference of power with ESS.</p>
Full article ">
19 pages, 1253 KiB  
Article
Taguchi Method and Numerical Simulation for Variable Viscosity and Non-Linear Boussinesq Effects on Natural Convection over a Vertical Truncated Cone in Porous Media
by Ken Ming Tu, Kuo Ann Yih, Fu I Chou and Jyh Horng Chou
Energies 2020, 13(2), 504; https://doi.org/10.3390/en13020504 - 20 Jan 2020
Cited by 3 | Viewed by 3121
Abstract
This study uses an optimization approach representation and numerical solution for the variable viscosity and non-linear Boussinesq effects on the free convection over a vertical truncated cone in porous media. The surface of the vertical truncated cone is maintained at uniform wall temperature [...] Read more.
This study uses an optimization approach representation and numerical solution for the variable viscosity and non-linear Boussinesq effects on the free convection over a vertical truncated cone in porous media. The surface of the vertical truncated cone is maintained at uniform wall temperature and uniform wall concentration (UWT/UWC). The viscosity of the fluid varies inversely to a linear function of the temperature. The partial differential equation is transformed into a non-similar equation and solved by Keller box method (KBM). Compared with previously published articles, the results are considered to be very consistent. Numerical results for the local Nusselt number and local Sherwood number with the six parameters (1) dimensionless streamwise coordinate ξ, (2) buoyancy ratio N, (3) Lewis number Le, (4) viscosity-variation parameter θ r , (5) non-linear temperature parameter δ 1 , and (6) non-linear concentration parameter δ 2 are expressed in figures and tables. The Taguchi method was used to predict the best point of the maxima of the local Nusselt (Sherwood) number of 3.8636 (5.1156), resulting in ξ (4), N (10), Le (0.5), θ r (−2), δ 1 (2), δ 2 (2) and ξ (4), N (10), Le (2), θ r (−2), δ 1 (2), δ 2 (2), respectively. Full article
(This article belongs to the Section J: Thermal Management)
Show Figures

Figure 1

Figure 1
<p>Flow model and physical coordinate system.</p>
Full article ">Figure 2
<p>The dimensionless temperature profile for two values of Le and <math display="inline"><semantics> <mrow> <msub> <mi mathvariant="sans-serif">θ</mi> <mi mathvariant="normal">r</mi> </msub> </mrow> </semantics></math>.</p>
Full article ">Figure 3
<p>The dimensionless concentration profile for two values of Le and <math display="inline"><semantics> <mrow> <msub> <mi mathvariant="sans-serif">θ</mi> <mi mathvariant="normal">r</mi> </msub> </mrow> </semantics></math>.</p>
Full article ">Figure 4
<p>The dimensionless temperature profile for two values of <math display="inline"><semantics> <mrow> <msub> <mi mathvariant="sans-serif">δ</mi> <mn>1</mn> </msub> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <msub> <mi mathvariant="sans-serif">δ</mi> <mn>2</mn> </msub> </mrow> </semantics></math>.</p>
Full article ">Figure 5
<p>The dimensionless concentration profile for two values of <math display="inline"><semantics> <mrow> <msub> <mi mathvariant="sans-serif">δ</mi> <mn>1</mn> </msub> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <msub> <mi mathvariant="sans-serif">δ</mi> <mn>2</mn> </msub> </mrow> </semantics></math>.</p>
Full article ">
29 pages, 2132 KiB  
Article
An All-At-Once Newton Strategy for Marine Methane Hydrate Reservoir Models
by Shubhangi Gupta, Barbara Wohlmuth and Matthias Haeckel
Energies 2020, 13(2), 503; https://doi.org/10.3390/en13020503 - 20 Jan 2020
Cited by 9 | Viewed by 3327
Abstract
The migration of methane through the gas hydrate stability zone (GHSZ) in the marine subsurface is characterized by highly dynamic reactive transport processes coupled to thermodynamic phase transitions between solid gas hydrates, free methane gas, and dissolved methane in the aqueous phase. The [...] Read more.
The migration of methane through the gas hydrate stability zone (GHSZ) in the marine subsurface is characterized by highly dynamic reactive transport processes coupled to thermodynamic phase transitions between solid gas hydrates, free methane gas, and dissolved methane in the aqueous phase. The marine subsurface is essentially a water-saturated porous medium where the thermodynamic instability of the hydrate phase can cause free gas pockets to appear and disappear locally, causing the model to degenerate. This poses serious convergence issues for the general-purpose nonlinear solvers (e.g., standard Newton), and often leads to extremely small time-step sizes. The convergence problem is particularly severe when the rate of hydrate phase change is much lower than the rate of gas dissolution. In order to overcome this numerical challenge, we have developed an all-at-once Newton scheme tailored to our gas hydrate model, which can handle rate-based hydrate phase change coupled with equilibrium gas dissolution in a mathematically consistent and robust manner. Full article
(This article belongs to the Special Issue Advances in Natural Gas Hydrates)
Show Figures

Figure 1

Figure 1
<p>Representation of the phases and components in a representative elementary volume (REV). For any pore-filling phase <math display="inline"> <semantics> <mrow> <mi>β</mi> <mo>=</mo> <mi>g</mi> <mo>,</mo> <mi>w</mi> <mo>,</mo> <mi>h</mi> </mrow> </semantics> </math>, phase saturation is defined as <math display="inline"> <semantics> <mrow> <msub> <mi>S</mi> <mi>β</mi> </msub> <mo>:</mo> <mo>=</mo> <mfrac> <msub> <mi>V</mi> <mi>β</mi> </msub> <msub> <mi>V</mi> <mi>p</mi> </msub> </mfrac> </mrow> </semantics> </math>. Total and apparent porosites are defined as <math display="inline"> <semantics> <mrow> <mi>ϕ</mi> <mo>:</mo> <mo>=</mo> <mfrac> <msub> <mi>V</mi> <mi>p</mi> </msub> <msub> <mi>V</mi> <mi>t</mi> </msub> </mfrac> </mrow> </semantics> </math>, and <math display="inline"> <semantics> <mrow> <msub> <mi>ϕ</mi> <mrow> <mi>e</mi> <mi>f</mi> <mi>f</mi> </mrow> </msub> <mo>:</mo> <mo>=</mo> <mfrac> <mrow> <msub> <mi>V</mi> <mi>p</mi> </msub> <mo>−</mo> <msub> <mi>V</mi> <mi>h</mi> </msub> </mrow> <msub> <mi>V</mi> <mi>t</mi> </msub> </mfrac> <mo>=</mo> <mfenced separators="" open="(" close=")"> <mn>1</mn> <mo>−</mo> <msub> <mi>S</mi> <mi>h</mi> </msub> </mfenced> <mi>ϕ</mi> </mrow> </semantics> </math>, respectively.</p>
Full article ">Figure 2
<p>Regional seismic profile across the western part of the Danube paleo delta in SW to NE direction, depicting the geological setting for Example 1 (<a href="#sec4dot1-energies-13-00503" class="html-sec">Section 4.1</a>). 2D RMCS line 09. Interpretation of the seismic data according to [<a href="#B50-energies-13-00503" class="html-bibr">50</a>]. BCL: buried channel-levee.</p>
Full article ">Figure 3
<p>Problem setting for Example 1 (<a href="#sec4dot1-energies-13-00503" class="html-sec">Section 4.1</a>). (<b>a</b>) The initial state of the system and identifies the corresponding gas hydrate stability zone (GHSZ). (<b>b</b>) The state of the system at <math display="inline"> <semantics> <mrow> <mi>t</mi> <mo>=</mo> <msub> <mi>t</mi> <mi>n</mi> </msub> <mo>&gt;</mo> <mn>0</mn> </mrow> </semantics> </math>, illustrating how the GHSZ shifts as a result of sedimentation over time.</p>
Full article ">Figure 4
<p>Numerical results for Example 1 (<a href="#sec4dot1-energies-13-00503" class="html-sec">Section 4.1</a>). The figure compares the NCP and the primary variable switching (PVS) schemes in terms of the cumulative CPU time required to solve the problem and the evolution of the time-step size during the simulation.</p>
Full article ">Figure 5
<p>Numerical results for Example 1 (<a href="#sec4dot1-energies-13-00503" class="html-sec">Section 4.1</a>). The figure shows snapshots of <math display="inline"> <semantics> <msub> <mi>S</mi> <mi>h</mi> </msub> </semantics> </math>, <math display="inline"> <semantics> <msub> <mi>S</mi> <mi>g</mi> </msub> </semantics> </math>, <math display="inline"> <semantics> <msubsup> <mi>χ</mi> <mi>w</mi> <mrow> <mi>C</mi> <msub> <mi>H</mi> <mn>4</mn> </msub> </mrow> </msubsup> </semantics> </math>, and <math display="inline"> <semantics> <msubsup> <mi>χ</mi> <mi>w</mi> <mi>c</mi> </msubsup> </semantics> </math> at <math display="inline"> <semantics> <mrow> <mi>t</mi> <mo>=</mo> <mn>22,500</mn> </mrow> </semantics> </math> years, that is, the time when the gas phase first appears; <math display="inline"> <semantics> <mrow> <mi>t</mi> <mo>=</mo> <mn>135,000</mn> </mrow> </semantics> </math> years, that is, the time up to which the PVS scheme solved in 240 CPU-hours; and <math display="inline"> <semantics> <mrow> <mi>t</mi> <mo>=</mo> <msub> <mi>t</mi> <mrow> <mi>e</mi> <mi>n</mi> <mi>d</mi> </mrow> </msub> <mo>=</mo> <mn>300,000</mn> </mrow> </semantics> </math> years. For <math display="inline"> <semantics> <mrow> <mi>t</mi> <mo>=</mo> <mn>22,500</mn> </mrow> </semantics> </math> years and <math display="inline"> <semantics> <mrow> <mi>t</mi> <mo>=</mo> <mn>135,000</mn> </mrow> </semantics> </math> years, the solutions of both NCP and PVS schemes is plotted for comparison.</p>
Full article ">Figure 5 Cont.
<p>Numerical results for Example 1 (<a href="#sec4dot1-energies-13-00503" class="html-sec">Section 4.1</a>). The figure shows snapshots of <math display="inline"> <semantics> <msub> <mi>S</mi> <mi>h</mi> </msub> </semantics> </math>, <math display="inline"> <semantics> <msub> <mi>S</mi> <mi>g</mi> </msub> </semantics> </math>, <math display="inline"> <semantics> <msubsup> <mi>χ</mi> <mi>w</mi> <mrow> <mi>C</mi> <msub> <mi>H</mi> <mn>4</mn> </msub> </mrow> </msubsup> </semantics> </math>, and <math display="inline"> <semantics> <msubsup> <mi>χ</mi> <mi>w</mi> <mi>c</mi> </msubsup> </semantics> </math> at <math display="inline"> <semantics> <mrow> <mi>t</mi> <mo>=</mo> <mn>22,500</mn> </mrow> </semantics> </math> years, that is, the time when the gas phase first appears; <math display="inline"> <semantics> <mrow> <mi>t</mi> <mo>=</mo> <mn>135,000</mn> </mrow> </semantics> </math> years, that is, the time up to which the PVS scheme solved in 240 CPU-hours; and <math display="inline"> <semantics> <mrow> <mi>t</mi> <mo>=</mo> <msub> <mi>t</mi> <mrow> <mi>e</mi> <mi>n</mi> <mi>d</mi> </mrow> </msub> <mo>=</mo> <mn>300,000</mn> </mrow> </semantics> </math> years. For <math display="inline"> <semantics> <mrow> <mi>t</mi> <mo>=</mo> <mn>22,500</mn> </mrow> </semantics> </math> years and <math display="inline"> <semantics> <mrow> <mi>t</mi> <mo>=</mo> <mn>135,000</mn> </mrow> </semantics> </math> years, the solutions of both NCP and PVS schemes is plotted for comparison.</p>
Full article ">Figure 6
<p>Numerical results for Example 1 (<a href="#sec4dot1-energies-13-00503" class="html-sec">Section 4.1</a>). Figure shows the process of hydrate dissociation → gas migration → hydrate reformation as a result of rising GHSZ between 60,000 years <math display="inline"> <semantics> <mrow> <mo>≤</mo> <mi>t</mi> <mo>≤</mo> </mrow> </semantics> </math> 120,000 years. The new gas hydrate layer grows using the methane gas supplied by the dissociating gas hydrate layer below.</p>
Full article ">Figure 6 Cont.
<p>Numerical results for Example 1 (<a href="#sec4dot1-energies-13-00503" class="html-sec">Section 4.1</a>). Figure shows the process of hydrate dissociation → gas migration → hydrate reformation as a result of rising GHSZ between 60,000 years <math display="inline"> <semantics> <mrow> <mo>≤</mo> <mi>t</mi> <mo>≤</mo> </mrow> </semantics> </math> 120,000 years. The new gas hydrate layer grows using the methane gas supplied by the dissociating gas hydrate layer below.</p>
Full article ">Figure 6 Cont.
<p>Numerical results for Example 1 (<a href="#sec4dot1-energies-13-00503" class="html-sec">Section 4.1</a>). Figure shows the process of hydrate dissociation → gas migration → hydrate reformation as a result of rising GHSZ between 60,000 years <math display="inline"> <semantics> <mrow> <mo>≤</mo> <mi>t</mi> <mo>≤</mo> </mrow> </semantics> </math> 120,000 years. The new gas hydrate layer grows using the methane gas supplied by the dissociating gas hydrate layer below.</p>
Full article ">Figure 7
<p>Problem setting for Example 2 (<a href="#sec4dot2-energies-13-00503" class="html-sec">Section 4.2</a>). Figure (<b>a</b>) highlights the essential features of the problem setting like the locations of the sea floor, the production well, the initial base of the GHSZ, and the initial hydrate distribution within the hydrate layer, and marks our domain of interest within the computational domain. Figure (<b>b</b>) identifies the relevant regions of the computational domain. <math display="inline"> <semantics> <mo>Ω</mo> </semantics> </math> denotes the computational domain, <math display="inline"> <semantics> <mrow> <msub> <mo>Ω</mo> <mi>H</mi> </msub> <mo>⊂</mo> <mo>Ω</mo> </mrow> </semantics> </math> denotes the hydrate layer, <math display="inline"> <semantics> <mrow> <msub> <mo>Ω</mo> <mi>I</mi> </msub> <mo>⊂</mo> <mo>Ω</mo> </mrow> </semantics> </math> denotes the domain of interest, and <math display="inline"> <semantics> <mrow> <msub> <mo>Ω</mo> <mi>H</mi> </msub> <mo>∩</mo> <msub> <mo>Ω</mo> <mi>I</mi> </msub> <mo>≠</mo> <mo>∅</mo> </mrow> </semantics> </math>. <math display="inline"> <semantics> <mrow> <mo>∂</mo> <msub> <mo>Ω</mo> <mrow> <mi>w</mi> <mi>e</mi> <mi>l</mi> <mi>l</mi> </mrow> </msub> </mrow> </semantics> </math> denotes the production well boundary, while <math display="inline"> <semantics> <mrow> <mo>∂</mo> <msub> <mo>Ω</mo> <mi>L</mi> </msub> </mrow> </semantics> </math> denotes the left boundary excluding the production well. <math display="inline"> <semantics> <mrow> <mo>∂</mo> <msub> <mo>Ω</mo> <mrow> <mi>s</mi> <mi>f</mi> </mrow> </msub> </mrow> </semantics> </math> denotes the top boundary corresponding to the sea floor. <math display="inline"> <semantics> <mrow> <mo>∂</mo> <msub> <mo>Ω</mo> <mi>R</mi> </msub> </mrow> </semantics> </math> and <math display="inline"> <semantics> <mrow> <mo>∂</mo> <msub> <mo>Ω</mo> <mi>B</mi> </msub> </mrow> </semantics> </math> denote the right and the bottom boundaries, respectively.</p>
Full article ">Figure 8
<p>Numerical results for Example 2 (<a href="#sec4dot2-energies-13-00503" class="html-sec">Section 4.2</a>). The figure shows snapshots of the quantities of interest (QoIs) (from left to right: <math display="inline"> <semantics> <msub> <mi>S</mi> <mi>h</mi> </msub> </semantics> </math>, <math display="inline"> <semantics> <msub> <mi>S</mi> <mi>g</mi> </msub> </semantics> </math>, <math display="inline"> <semantics> <msubsup> <mi>χ</mi> <mi>w</mi> <mrow> <mi>C</mi> <msub> <mi>H</mi> <mn>4</mn> </msub> </mrow> </msubsup> </semantics> </math>, and the GHSZ) within the domain of interest <math display="inline"> <semantics> <msub> <mo>Ω</mo> <mi>I</mi> </msub> </semantics> </math> at different times. Note that for GHSZ, a value of 1 indicates an unstable zone, and <math display="inline"> <semantics> <mrow> <mo>−</mo> <mn>1</mn> </mrow> </semantics> </math> indicates a stable zone.</p>
Full article ">Figure 9
<p>Numerical results for the reference test case of Example 2 (<a href="#sec4dot2-energies-13-00503" class="html-sec">Section 4.2</a>). The figure shows snapshots of the QoIs (from left to right: <math display="inline"> <semantics> <msub> <mi>S</mi> <mi>h</mi> </msub> </semantics> </math>, <math display="inline"> <semantics> <msub> <mi>S</mi> <mi>g</mi> </msub> </semantics> </math>, <math display="inline"> <semantics> <msubsup> <mi>χ</mi> <mi>w</mi> <mrow> <mi>C</mi> <msub> <mi>H</mi> <mn>4</mn> </msub> </mrow> </msubsup> </semantics> </math>, and the GHSZ) within the domain of interest <math display="inline"> <semantics> <msub> <mo>Ω</mo> <mi>I</mi> </msub> </semantics> </math> at <math display="inline"> <semantics> <mrow> <mi>t</mi> <mo>=</mo> <msub> <mi>t</mi> <mrow> <mi>e</mi> <mi>n</mi> <mi>d</mi> </mrow> </msub> <mo>=</mo> <mn>360</mn> </mrow> </semantics> </math> days. Note that for the GHSZ, 1 indicates an unstable zone, and <math display="inline"> <semantics> <mrow> <mo>−</mo> <mn>1</mn> </mrow> </semantics> </math> indicates a stable zone.</p>
Full article ">Figure 10
<p>Numerical result for Example 2 (<a href="#sec4dot2-energies-13-00503" class="html-sec">Section 4.2</a>). A comparison of the time-step size evolution for a random hydrate phase distribution and a homogeneous hydrate phase distribution.</p>
Full article ">
12 pages, 2716 KiB  
Article
Analysis of the Theoretical Performance of the Wind-Driven Pulverizing Aerator in the Conditions of Góreckie Lake—Maximum Wind Speed Method
by Ewa Osuch, Andrzej Osuch, Piotr Rybacki and Andrzej Przybylak
Energies 2020, 13(2), 502; https://doi.org/10.3390/en13020502 - 20 Jan 2020
Cited by 6 | Viewed by 2390
Abstract
The eutrophication of surface waters is a natural process; however, anthropogenic activities significantly accelerate degradation processes. Most lakes in Poland and in the world belong to the poor and unsatisfactory water quality class. It is therefore necessary to limit negative anthropogenic impacts and [...] Read more.
The eutrophication of surface waters is a natural process; however, anthropogenic activities significantly accelerate degradation processes. Most lakes in Poland and in the world belong to the poor and unsatisfactory water quality class. It is therefore necessary to limit negative anthropogenic impacts and introduce restoration methods, in particular those that are safe for the aquatic ecosystem. One of these is a pulverizing aeration Podsiadłowski method that uses only wind energy. The method allows for the moderate oxygenation of hypolimnion water, which maintains the oxygen conditions in the overlying water zone in the range of 0–1 mg O2·dm-1. The purpose of the work was to develop a new method of determining the efficiency of the aerator pulverization unit in the windy conditions of the lake. The method consists in determining the volumetric flow rates of water in the aerator pulverization unit, based on maximum hourly wind speeds. The pulverization efficiency in the conditions of Góreckie Lake was determined based on 6600 maximum hourly wind speeds in 2018. Based on the determined model, the theoretical performance of the machine was calculated, which in the conditions of Góreckie Lake in 2018 amounted to less than 79,000 m3 per year (nine months of the effective aerator operation). Full article
Show Figures

Figure 1

Figure 1
<p>Wind-driven pulverizing aerator located on Góreckie Lake.</p>
Full article ">Figure 2
<p>Scheme of the pulverization unit in the aerator. Note: 1, coagulant applicator; 2, pressuring chamber; 3, pulverizing circle; 4, suction chamber; 5, pressuring hoses; 6, suction hoses.</p>
Full article ">Figure 3
<p>Rotational speed model of the aerator pulverization wheel in relation to wind speed.</p>
Full article ">Figure 4
<p>Efficiency of the pulverization unit of the aerator in March.</p>
Full article ">Figure 5
<p>Efficiency of the pulverization unit of the aerator in April.</p>
Full article ">Figure 6
<p>Efficiency of the pulverization unit of the aerator in May.</p>
Full article ">Figure 7
<p>Efficiency of the pulverization unit of the aerator in June.</p>
Full article ">Figure 8
<p>Efficiency of the pulverization unit of the aerator in July.</p>
Full article ">Figure 9
<p>Efficiency of the pulverization unit of the aerator in August.</p>
Full article ">Figure 10
<p>Efficiency of the pulverization unit of the aerator in September.</p>
Full article ">Figure 11
<p>Efficiency of the pulverization unit of the aerator in October.</p>
Full article ">Figure 12
<p>Efficiency of the pulverization unit of the aerator in November.</p>
Full article ">Figure 13
<p>Efficiency of the pulverization unit of the aerator in 2018.</p>
Full article ">
20 pages, 5864 KiB  
Article
Dynamic Modeling of a Parallel-Connected Solid Oxide Fuel Cell Stack System
by Chien-Chang Wu and Tsung-Lin Chen
Energies 2020, 13(2), 501; https://doi.org/10.3390/en13020501 - 20 Jan 2020
Cited by 6 | Viewed by 2943
Abstract
This study proposes novel simulation methods to model the power delivery function of a parallel-connected solid-oxide-fuel-cell stack system. The proposed methods are then used to investigate the possible thermal runaway induced by the performance mismatch between the employed stacks. A challenge in this [...] Read more.
This study proposes novel simulation methods to model the power delivery function of a parallel-connected solid-oxide-fuel-cell stack system. The proposed methods are then used to investigate the possible thermal runaway induced by the performance mismatch between the employed stacks. A challenge in this modeling study is to achieve the same output voltage but different output current for each employed stack. Conventional fuel-cell models cannot be used, because they employ fuel flow rates and stack currents as the input variables. These two variables are unknown in the parallel-connected stack systems. The proposed method solves the aforementioned problems by integrating the fuel supply dynamics with the conventional stack models and then arranging them in a multiple-feedback-loop configuration for conducting simulations. The simulation results indicate that the proposed methods can model the transient response of the parallel-connected stack system. Moreover, for the dynamics of the power distribution, there exists an unstable positive feedback loop between employed stacks when the stack temperatures are low, and a stable negative feedback loop when the stack temperatures are high. A thermal runaway could be initiated when the dynamics of the stack temperature is slower than that of the current distribution. Full article
(This article belongs to the Section A5: Hydrogen Energy)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Schematic of the two fuel-cell (FC) stacks connected in parallel. The stack system is followed by a DC–DC converter to deliver power to the load. The green arrows represent the fuel supplying and exiting the stacks; the blue arrows represent electricity output; the red arrows represent the heat exchange between stacks.</p>
Full article ">Figure 2
<p>Proposed simulation method for the power-delivery operation of the parallel-connected stack system.</p>
Full article ">Figure 3
<p>Comparisons between the I–V curve data of a SOFC system obtained from three sources: experimental data excerpted from [<a href="#B31-energies-13-00501" class="html-bibr">31</a>], simulation data obtained from the static electrochemical model presented in this study, and simulation data obtained from the dynamic model stated in this study.</p>
Full article ">Figure 4
<p>Transient response of the SOFC system obtained when the two employed stacks are the same. The total requested power was 1500 W. Case 1: stacks have a small ohmic loss. Case 2: stacks have a large ohmic loss. (<b>a</b>) Stack output current, (<b>b</b>) stack voltage, (<b>c</b>) stack output power, and (<b>d</b>) stack temperature.</p>
Full article ">Figure 5
<p>Transient response of the stack voltages when the two employed FC stacks are the same. Case 1: stacks have a small ohmic loss. Case 2: stacks have a large ohmic loss. (<b>a</b>) Nernst voltage, (<b>b</b>) activation polarization, (<b>c</b>) ohmic polarization, and (<b>d</b>) concentration polarization.</p>
Full article ">Figure 6
<p>Transient response of the channel pressures and fuel utilizations when the two employed FC stacks are the same. Case 1: stacks have a small ohmic loss. Case 2: stacks have a large ohmic loss. (<b>a</b>) Anode channel pressure, (<b>b</b>) cathode channel pressure, (<b>c</b>) hydrogen utilization, and (<b>d</b>) oxygen utilization.</p>
Full article ">Figure 7
<p>Transient response of the FC system obtained when the two employed FC stacks have different ohmic properties. Case 1 includes the heat exchange between stacks, whereas case 2 excludes it. (<b>a</b>) Stack output current, (<b>b</b>) stack voltage, (<b>c</b>) stack output power, and (<b>d</b>) stack temperature. The system encounters a thermal runaway when there is no heat exchange between stacks.</p>
Full article ">Figure 8
<p>Transient response of the stack voltages when the two employed FC stacks have different ohmic properties. Case 1 includes the heat exchange between stacks, whereas case 2 excludes it. (<b>a</b>) Nernst voltage, (<b>b</b>) activation polarization, (<b>c</b>) ohmic polarization, and (<b>d</b>) concentration polarization. Stack voltage differs largely when there is no heat exchange between stacks.</p>
Full article ">Figure 9
<p>Transient response of the channel pressures and fuel utilizations when the two employed FC stacks have different ohmic properties. Case 1 includes the heat exchange between stacks, whereas case 2 excludes it. (<b>a</b>) Anode channel pressure, (<b>b</b>) cathode channel pressure, (<b>c</b>) hydrogen utilization, and (<b>d</b>) oxygen utilization. The channel pressures are almost the same, but the hydrogen utilization rate differs largely.</p>
Full article ">Figure 10
<p>Transient response of the parallel-connected stacks when the two employed stacks have different ohmic properties and a lower heat capacity. Case 1 includes the heat exchange between stacks, whereas case 2 excludes it. (<b>a</b>) Stack output current, (<b>b</b>) stack voltage, (<b>c</b>) stack output power, and (<b>d</b>) stack temperature. The steady-state performance of the two cases and stacks are almost the same.</p>
Full article ">Figure 11
<p>FC stack voltages vary with the stack temperature and output current.</p>
Full article ">
23 pages, 924 KiB  
Article
International Trade Disputes over Renewable Energy—the Case of the Solar Photovoltaic Sector
by Agnieszka Hajdukiewicz and Bożena Pera
Energies 2020, 13(2), 500; https://doi.org/10.3390/en13020500 - 20 Jan 2020
Cited by 14 | Viewed by 6552
Abstract
The development of the renewable energy industry is a priority of economic policies in many countries, since it is viewed as one of the key growth sectors in the economy, playing also a very important role in mitigating climate change. At the international [...] Read more.
The development of the renewable energy industry is a priority of economic policies in many countries, since it is viewed as one of the key growth sectors in the economy, playing also a very important role in mitigating climate change. At the international level, renewable energy is an issue of international cooperation but also an area of high trade tensions between countries. The main goal of this paper is to examine the nature and sources of recent trade disputes in the solar photovoltaic sector, which is the most dynamically growing sector in the green energy industry. In particular, the paper explores the links between the contemporary trade disputes and modern protectionism and between protectionist policies and practices and the export competitiveness in the growing sector of the economy. To achieve the aim of the study we explore in detail the WTO trade disputes over photovoltaic (PV) products, which occurred in the years 2007–2018. The products covered by the analysis were solar modules and cells classified under the HS code 854140. In our research we also used measures of descriptive statistics, hierarchical cluster analysis and revealed comparative advantage indexes. Our key results demonstrate the existence of links between protectionist policy causing trade conflicts and the export competitiveness. The research has also allowed us to identify problems of future studies concerning the association between trade protectionism and global value chains in the solar energy sector. Full article
(This article belongs to the Special Issue Political Economy of Energy Policies)
Show Figures

Figure 1

Figure 1
<p>Research conceptualization (research strategy).</p>
Full article ">Figure 2
<p>Share of analyzed countries in total exports and concentration of world’s photovoltaic panels exports in 2007–2018. Source: own elaboration based on the ITC Trade Map [<a href="#B84-energies-13-00500" class="html-bibr">84</a>].</p>
Full article ">Figure 3
<p>Revealed Comparative Advantage (<span class="html-italic">RCA</span>) index versus country’s exports share in the world exports of photovoltaic panels. Source: own elaboration.</p>
Full article ">
27 pages, 2573 KiB  
Article
Is Investing in Companies Manufacturing Solar Components a Lucrative Business? A Decision Tree Based Analysis
by Sebastian Klaudiusz Tomczak, Anna Skowrońska-Szmer and Jan Jakub Szczygielski
Energies 2020, 13(2), 499; https://doi.org/10.3390/en13020499 - 20 Jan 2020
Cited by 5 | Viewed by 3041
Abstract
In an era of increasing energy production from renewable sources, the demand for components for renewable energy systems has dramatically increased. Consequently, managers and investors are interested in knowing whether a company associated with the semiconductor and related device manufacturing sector, especially the [...] Read more.
In an era of increasing energy production from renewable sources, the demand for components for renewable energy systems has dramatically increased. Consequently, managers and investors are interested in knowing whether a company associated with the semiconductor and related device manufacturing sector, especially the photovoltaic (PV) systems manufacturers, is a money-making business. We apply a new approach that extends prior research by applying decision trees (DTs) to identify ratios (i.e., indicators), which discriminate between companies within the sector that do (designated as “green”) and do not (“red”) produce elements of PV systems. Our results indicate that on the basis of selected ratios, green companies can be distinguished from the red companies without an in-depth analysis of the product portfolio. We also find that green companies, especially operating in China are characterized by lower financial performance, thus providing a negative (and unexpected) answer to the question posed in the title. Full article
(This article belongs to the Special Issue Economics of Sustainable and Renewable Energy Systems)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Global electricity generation from renewable sources, 1990–2017. Data source: [<a href="#B1-energies-13-00499" class="html-bibr">1</a>].</p>
Full article ">Figure 2
<p>Electricity generation from renewable sources; share in the overall energy production for selected countries, 2010–2017 (in %). Data source—[<a href="#B1-energies-13-00499" class="html-bibr">1</a>,<a href="#B4-energies-13-00499" class="html-bibr">4</a>,<a href="#B5-energies-13-00499" class="html-bibr">5</a>,<a href="#B6-energies-13-00499" class="html-bibr">6</a>]. Target for Turkey, People’s Republic of China, India is for 2030, target for Chinese Taipei is for 2025 and target for Thailand is for 2036; * means the percentage of non-fossil energy in total electric power generation.</p>
Full article ">Figure 3
<p>Return on assets (ROA) values for red and green companies.</p>
Full article ">Figure 4
<p>Return on sales (ROS) values for red and green companies.</p>
Full article ">Figure 5
<p>Return on equity (ROE) values for red and green companies.</p>
Full article ">
19 pages, 6542 KiB  
Article
Measurements Based Analysis of the Proton Exchange Membrane Fuel Cell Operation in Transient State and Power of Own Needs
by Andrzej Wilk and Daniel Węcel
Energies 2020, 13(2), 498; https://doi.org/10.3390/en13020498 - 20 Jan 2020
Cited by 1 | Viewed by 3036
Abstract
Currently, fuel cells are increasingly used in industrial installations, means of transport, and household applications as a source of electricity and heat. The paper presents the results of experimental tests of a Proton Exchange Membrane Fuel Cell (PEMFC) at variable load, which characterizes [...] Read more.
Currently, fuel cells are increasingly used in industrial installations, means of transport, and household applications as a source of electricity and heat. The paper presents the results of experimental tests of a Proton Exchange Membrane Fuel Cell (PEMFC) at variable load, which characterizes the cell’s operation in real installations. A detailed analysis of the power needed for operation fuel cell auxiliary devices (own needs power) was carried out. An analysis of net and gross efficiency was carried out in various operating conditions of the device. The measurements made show changes in the performance of the fuel cell during step changing or smooth changing of an electric load. Load was carried out as a change in the current or a change in the resistance of the receiver. The analysis covered the times of reaching steady states and the efficiency of the fuel cell system taking into account auxiliary devices. In the final part of the article, an analysis was made of the influence of the fuel cell duration of use on obtained parameters. The analysis of the measurement results will allow determination of the possibility of using fuel cells in installations with a rapidly changing load profile and indicate possible solutions to improve the performance of the installation. Full article
Show Figures

Figure 1

Figure 1
<p>Diagram of the fuel cell stack measurement system.</p>
Full article ">Figure 2
<p>Fuel cell characteristics and operating area.</p>
Full article ">Figure 3
<p>Power of fuel cell auxiliary devices as a function of electric current.</p>
Full article ">Figure 4
<p>Power of fuel cell auxiliary devices as a function of excess air ratio.</p>
Full article ">Figure 5
<p>Fuel cell auxiliary power index as a function of current.</p>
Full article ">Figure 6
<p>Fuel cell characteristics obtained by means of measurements.</p>
Full article ">Figure 7
<p>Changes in fuel cell power at different inlet coolant temperature T<sub>1</sub> and different excess air ratio λ.</p>
Full article ">Figure 8
<p>Characteristics of the gross efficiency as a function of excess air ratio λ for T<sub>1</sub> = 54 °C.</p>
Full article ">Figure 9
<p>Characteristics of the net efficiency as a function of air excess ratio λ for T<sub>1</sub> = 54 °C.</p>
Full article ">Figure 10
<p>Fuel cell parameters changes at rapid increase of load.</p>
Full article ">Figure 11
<p>Voltage fluctuation rate ε<sub>V</sub> and hydrogen flow with time under step resistance change.</p>
Full article ">Figure 12
<p>Changes in fuel cell parameters during step increasing of the load.</p>
Full article ">Figure 13
<p>Changes in fuel cell parameters during step decreasing of the load.</p>
Full article ">Figure 14
<p>Voltage fluctuation rate ε<sub>V</sub> with time under load increase.</p>
Full article ">Figure 15
<p>Voltage fluctuation rate ε<sub>V</sub> with time under load decrease.</p>
Full article ">Figure 16
<p>Fuel cell characteristics V = f(I) obtained during many years of operation.</p>
Full article ">Figure 17
<p>Fuel cell characteristics P = f(I) obtained during many years of operation.</p>
Full article ">
15 pages, 1327 KiB  
Article
Parameter Estimation of Three Diode Photovoltaic Model Using Grasshopper Optimization Algorithm
by Omnia S. Elazab, Hany M. Hasanien, Ibrahim Alsaidan, Almoataz Y. Abdelaziz and S. M. Muyeen
Energies 2020, 13(2), 497; https://doi.org/10.3390/en13020497 - 20 Jan 2020
Cited by 97 | Viewed by 4220
Abstract
While addressing the issue of improving the performance of Photovoltaic (PV) systems, the simulation results are highly influenced by the PV model accuracy. Building the PV module mathematical model is based on its I-V characteristic, which is a highly nonlinear relationship. All the [...] Read more.
While addressing the issue of improving the performance of Photovoltaic (PV) systems, the simulation results are highly influenced by the PV model accuracy. Building the PV module mathematical model is based on its I-V characteristic, which is a highly nonlinear relationship. All the PV cells’ data sheets do not provide full information about their parameters. This leads to a nonlinear mathematical model with several unknown parameters. This paper proposes a new application of the Grasshopper Optimization Algorithm (GOA) for parameter extraction of the three-diode PV model of a PV module. Two commercial PV modules, Kyocera KC200GT and Solarex MSX-60 PV cells are utilized in examining the GOA-based PV model. The simulation results are executed under various temperatures and irradiations. The proposed PV model is evaluated by comparing its results with the experimental results of these commercial PV modules. The efficiency of the GOA-based PV model is tested by making a fair comparison among its numerical results and other optimization method-based PV models. With the GOA, a precise three-diode PV model shall be established. Full article
(This article belongs to the Section A2: Solar Energy and Photovoltaic Systems)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Three diode Photovoltaic (PV) model.</p>
Full article ">Figure 2
<p>Pseudo code of the Grasshopper Optimization Algorithm (GOA).</p>
Full article ">Figure 3
<p>Objective function convergence. (<b>a</b>) KC200GT; (<b>b</b>) MSX60.</p>
Full article ">Figure 4
<p>GOA-based simulation results and practical data of KC200GT module at various temperature conditions, G = 1000 W/m<sup>2</sup>. (<b>a</b>) I-V curves; (<b>b</b>) P-V curves.</p>
Full article ">Figure 5
<p>GOA-based simulation results and practical data of KC200GT module at various irradiation conditions, temperature = 25 °C (<b>a</b>) I-V curves; (<b>b</b>) P-V curves.</p>
Full article ">Figure 6
<p>GOA-based simulation results and practical data of MSX-60 module at various temperature conditions, G = 1000 W/m<sup>2</sup>. (<b>a</b>) I-V curves. (<b>b</b>) P-V curves.</p>
Full article ">Figure 7
<p>Absolute current error of the GOA-based PV model, MATLAB model, whale optimization algorithm (WOA) model, genetic algorithm (GA) model, and simulated annealing technique (SA) model. (<b>a</b>) KC200GT. (<b>b</b>) MSX-60.</p>
Full article ">
17 pages, 6431 KiB  
Review
Review of Time and Space Harmonics in Multi-Phase Induction Machine
by Vladimir Kindl, Radek Cermak, Zelmira Ferkova and Bohumil Skala
Energies 2020, 13(2), 496; https://doi.org/10.3390/en13020496 - 19 Jan 2020
Cited by 22 | Viewed by 7153
Abstract
Modern multiphase electric machines take advantage of additional degrees of freedom for various purposes, including harmonic current injection to increase torque per ampere. This new approach introduces a non-sinusoidal air gap flux density distribution causing additional technical problems and so the conventional assumptions [...] Read more.
Modern multiphase electric machines take advantage of additional degrees of freedom for various purposes, including harmonic current injection to increase torque per ampere. This new approach introduces a non-sinusoidal air gap flux density distribution causing additional technical problems and so the conventional assumptions need to be revised. The paper presents a methodology for synthesis of air gap magnetic field generated by a symmetrically distributed multiphase windings including the rotor field reaction due to the machine’s load. The proposed method is suitable either for single-layer or double layer windings and can be adopted either for full-pitched or chorded winding including slots effects. The article analyses the air gap flux density harmonic content and formulates conclusions important to multiphase induction motors. It also discusses effects of time harmonic currents and illustrates the principle of changing number of pole-pairs typical for harmonic currents being injected to increase torque. Full article
(This article belongs to the Special Issue Advances in Rotating Electric Machines)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>(<b>a</b>) Air gap field generated by 3-phase winding; field distribution (<b>left</b>), (<b>b</b>) frequency spectrum (<b>right</b>).</p>
Full article ">Figure 2
<p>(<b>a</b>) Air gap field generated by 5-phase winding; field distribution (<b>left</b>), (<b>b</b>) frequency spectrum (<b>right</b>).</p>
Full article ">Figure 3
<p>(<b>a</b>) Air gap field generated by 7-phase winding; field distribution (<b>left</b>), (<b>b</b>) frequency spectrum (<b>right</b>).</p>
Full article ">Figure 4
<p>(<b>a</b>) No-load flux density distribution inside the machine, (<b>b</b>) air gap flux density.</p>
Full article ">Figure 5
<p>(<b>a</b>) Flux density distribution inside the machine under load, (<b>b</b>) air gap flux density.</p>
Full article ">Figure 6
<p>(<b>a</b>) Comparison between no-load and loaded operational state; field distribution (<b>b</b>) from FEA, frequency spectrum (right).</p>
Full article ">Figure 7
<p>(<b>a</b>) Time dependency of bar current (<b>left</b>) and (<b>b</b>) its harmonic spectrum (<b>right</b>).</p>
Full article ">Figure 8
<p>(<b>a</b>) Final air gap flux density for motor operating under load condition; flux density distribution (<b>left</b>)—blue line inverted, (<b>b</b>) frequency spectrum (<b>right</b>).</p>
Full article ">Figure 9
<p>(<b>a</b>) Stator phases re-assembling, (<b>b</b>) mmf for the fundamental and the 3rd harmonics.</p>
Full article ">Figure 10
<p>Time harmonics acting in multiphase winding: (<b>a</b>) harmonics suitable for injection, (<b>b</b>) harmonics inappropriate for injecting.</p>
Full article ">Figure 11
<p>Time harmonics effect demonstration in three-phase machines.</p>
Full article ">Figure 12
<p>Time harmonics effect demonstration in five-phase machines.</p>
Full article ">Figure 13
<p>Time harmonics effect demonstration in seven-phase machines.</p>
Full article ">
20 pages, 3241 KiB  
Article
Research on Additional Control Technology Based on Energy Storage System for Improving Power Transfer Capacity of Multi-Terminal AC/DC System with Low Cost
by Zheng Wu, Laifu Li, Yubo Yuan, Xiaodong Yuan, Chenyu Zhang, Li Kong, Wei Pei and Wei Deng
Energies 2020, 13(2), 495; https://doi.org/10.3390/en13020495 - 19 Jan 2020
Viewed by 1973
Abstract
The multi-terminal AC/DC system will become one of the important forms of the future power grid. The negative impedance characteristic caused by the constant power load in the DC network will reduce the power transfer capacity between the terminals, especially when a grid [...] Read more.
The multi-terminal AC/DC system will become one of the important forms of the future power grid. The negative impedance characteristic caused by the constant power load in the DC network will reduce the power transfer capacity between the terminals, especially when a grid fault occurs in AC system at any terminal. Energy storage has played an important role in improving the stability of AC and DC systems. This paper proposes an additional control method based on an energy storage system to improve system power transfer capacity with low cost. The state space model of two-terminal AC/DC system is established, and the feedback laws for additional control are further designed by Lyapunov theory. Furthermore, the additional control strategies based on the energy storage system is built, without changing the existing control system of each control object. Finally, the corresponding system simulation model is established by Matlab/Simulink for analysis and verification. The research results show that the proposed additional control method is effective. The power transfer limitation can be overcome by only adding small damping energy with the stable DC voltages under large disturbances, and the power transfer capacity between the terminals can be significantly improved with low control cost. Full article
(This article belongs to the Section A1: Smart Grids and Microgrids)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Typical structure of the multi-terminal AC/DC system.</p>
Full article ">Figure 2
<p>Classic VSC control strategies.</p>
Full article ">Figure 3
<p>System structure and the equivalent circuit. (<b>a</b>) System structure under failure at one VSC; (<b>b</b>) Equivalent circuit.</p>
Full article ">Figure 4
<p>The power transfer boundary. (<b>a</b>) P<sub>bus</sub> from −180 kW to 180 kW; (<b>b</b>) P<sub>bus</sub> from −20 kW to 50 kW.</p>
Full article ">Figure 5
<p>The simulation result when <span class="html-italic">P<sub>s</sub></span> changes from 90kW to 95kW: (<b>a</b>) The DC bus voltage; (<b>b</b>) The DC-side voltage of VSC2.</p>
Full article ">Figure 6
<p>The simulation result when <span class="html-italic">P<sub>s</sub></span> changes from 90kW to 98kW: (<b>a</b>) The DC bus voltage; (<b>b</b>) The DC-side voltage of VSC2.</p>
Full article ">Figure 7
<p>Additional control strategy.</p>
Full article ">Figure 8
<p>The power transfer boundary after applying the additional control based on the ESS.</p>
Full article ">Figure 9
<p>Simulation results of <span class="html-italic">U<sub>dc</sub></span>: (<b>a</b>) <span class="html-italic">P<sub>s</sub></span> changes from 90 kW to 100 kW; (<b>b</b>) <span class="html-italic">P<sub>s</sub></span> changes from 90 kW to 120 kW; (c) <span class="html-italic">P<sub>s</sub></span> changes from 90 kW to 150 kW; (d) <span class="html-italic">P<sub>s</sub></span> changes from 90 kW to 180 kW.</p>
Full article ">Figure 10
<p>Simulation results of <span class="html-italic">U<sub>s</sub></span>: (<b>a</b>) <span class="html-italic">P<sub>s</sub></span> changes from 90 kW to 100 kW; (<b>b</b>) <span class="html-italic">P<sub>s</sub></span> changes from 90 kW to 120 kW; (<b>c</b>) <span class="html-italic">P<sub>s</sub></span> changes from 90 kW to 150 kW; (<b>d</b>) <span class="html-italic">P<sub>s</sub></span> changes from 90 kW to 180 kW.</p>
Full article ">Figure 11
<p>The curves of the additional power of the ESS: (<b>a</b>) <span class="html-italic">P<sub>s</sub></span> changes from 90 kW to 100 kW; (<b>b</b>) <span class="html-italic">P<sub>s</sub></span> changes from 90 kW to 120 kW; (<b>c</b>) <span class="html-italic">P<sub>s</sub></span> changes from 90 kW to 150 kW; (<b>d</b>) <span class="html-italic">P<sub>s</sub></span> changes from 90 kW to 180 kW.</p>
Full article ">Figure 12
<p>The performance comparison with the implementation of the addition control: (<b>a</b>) DC bus voltage; (<b>b</b>) DC-side voltage of VSC2; (<b>c</b>) The additional power of the ESS; (<b>d</b>) The steady-state recovering time.</p>
Full article ">Figure 13
<p>The sensitivity analysis of the weighting parameters: (<b>a</b>) The power transfer boundary under different values of k at <b>Q</b>[1,1]; (<b>b</b>) The power transfer boundary under different values of k at <b>Q</b>[2,2]; (<b>c</b>) The power transfer boundary under different values of k at <b>Q</b>[3,3]; (<b>d</b>) The power transfer boundary under different values of k at <b>Q</b>[4,4].</p>
Full article ">Figure 13 Cont.
<p>The sensitivity analysis of the weighting parameters: (<b>a</b>) The power transfer boundary under different values of k at <b>Q</b>[1,1]; (<b>b</b>) The power transfer boundary under different values of k at <b>Q</b>[2,2]; (<b>c</b>) The power transfer boundary under different values of k at <b>Q</b>[3,3]; (<b>d</b>) The power transfer boundary under different values of k at <b>Q</b>[4,4].</p>
Full article ">
27 pages, 1267 KiB  
Review
Internet of Things (IoT) and the Energy Sector
by Naser Hossein Motlagh, Mahsa Mohammadrezaei, Julian Hunt and Behnam Zakeri
Energies 2020, 13(2), 494; https://doi.org/10.3390/en13020494 - 19 Jan 2020
Cited by 482 | Viewed by 39725
Abstract
Integration of renewable energy and optimization of energy use are key enablers of sustainable energy transitions and mitigating climate change. Modern technologies such the Internet of Things (IoT) offer a wide number of applications in the energy sector, i.e, in energy supply, transmission [...] Read more.
Integration of renewable energy and optimization of energy use are key enablers of sustainable energy transitions and mitigating climate change. Modern technologies such the Internet of Things (IoT) offer a wide number of applications in the energy sector, i.e, in energy supply, transmission and distribution, and demand. IoT can be employed for improving energy efficiency, increasing the share of renewable energy, and reducing environmental impacts of the energy use. This paper reviews the existing literature on the application of IoT in in energy systems, in general, and in the context of smart grids particularly. Furthermore, we discuss enabling technologies of IoT, including cloud computing and different platforms for data analysis. Furthermore, we review challenges of deploying IoT in the energy sector, including privacy and security, with some solutions to these challenges such as blockchain technology. This survey provides energy policy-makers, energy economists, and managers with an overview of the role of IoT in optimization of energy systems. Full article
(This article belongs to the Special Issue Sustainable Energy Systems: Efficiency and Optimization)
Show Figures

Figure 1

Figure 1
<p>Energy supply chain.</p>
Full article ">Figure 2
<p>Diagram describing the components of an IoT platform.</p>
Full article ">Figure 3
<p>A centralized data connectivity in a smart city concept.</p>
Full article ">Figure 4
<p>Share of residential energy consumption.</p>
Full article ">Figure 5
<p>Applications of IoT in an integrated smart energy system.</p>
Full article ">
18 pages, 5059 KiB  
Article
Bio-Crude Production through Aqueous Phase Recycling of Hydrothermal Liquefaction of Sewage Sludge
by Ayaz A. Shah, Saqib S. Toor, Tahir H. Seehar, Rasmus S. Nielsen, Asbjørn H. Nielsen, Thomas H. Pedersen and Lasse A. Rosendahl
Energies 2020, 13(2), 493; https://doi.org/10.3390/en13020493 - 19 Jan 2020
Cited by 56 | Viewed by 4862
Abstract
Hydrothermal liquefaction (HTL) is a promising technology for the production of bio-crude. However, some unresolved issues still exist within HTL, which need to be resolved before its promotion on a commercial scale. The management of the aqueous phase is one of the leading [...] Read more.
Hydrothermal liquefaction (HTL) is a promising technology for the production of bio-crude. However, some unresolved issues still exist within HTL, which need to be resolved before its promotion on a commercial scale. The management of the aqueous phase is one of the leading challenges related to HTL. In this study, the sewage sludge has been liquefied at 350 °C with and without catalyst (K2CO3). Subsequently, aqueous phase recycling was applied to investigate the effect of recycling on bio-crude properties. Obtained results showed that the energy recovery in the form of bio-crude increased by 50% via aqueous phase recirculation, whereas nitrogen content in the bio-crude was approximately doubled after eight rounds of recycling. GCMS characterization of the aqueous phase indicated acetic acid as a major water-soluble compound, which employed as a catalyst (0.56 M), and resulted in a negligible increase in bio-crude yield. ICP-AES highlighted that the majority of the inorganics were transferred to the solid phase, while the higher accumulation of potassium and sodium was found in the aqueous phase via successive rounds of recycling. Full article
(This article belongs to the Special Issue Thermochemical Biorefining II)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Optimization of the aqueous phase via (<b>a</b>) different processes, (<b>b</b>) present study.</p>
Full article ">Figure 2
<p>Schematic diagram of the HTL process with the recycling of the aqueous phase.</p>
Full article ">Figure 3
<p>Product yield at sub and supercritical conditions.</p>
Full article ">Figure 4
<p>Effect of aqueous phase recycling on product yield.</p>
Full article ">Figure 5
<p>Effect of concentration of acetic acid on bio-crude yield.</p>
Full article ">Figure 6
<p>Van Krevelan diagram shows the pathways of liquefaction products.</p>
Full article ">Figure 7
<p>(<b>a</b>) Volatility curves of bio-crudes, and (<b>b</b>) DTG curves of bio-crudes.</p>
Full article ">Figure 8
<p>(<b>a</b>) Effect of aqueous phase recycling on the composition of bio-crudes, and (<b>b</b>) chromatograms of bio-crudes.</p>
Full article ">Figure 9
<p>Carbon recovery in HTL products via aqueous phase recycling.</p>
Full article ">Figure 10
<p>TOC, TN, and pH of aqueous phase.</p>
Full article ">Figure 11
<p>Distribution of inorganic elements at (<b>a</b>) C0, (<b>b</b>) C5, and (<b>c</b>) Cat-acetic acid.</p>
Full article ">
18 pages, 2924 KiB  
Article
Simplified Sensorless Current Predictive Control of Synchronous Reluctance Motor Using Online Parameter Estimation
by Ahmed Farhan, Mohamed Abdelrahem, Amr Saleh, Adel Shaltout and Ralph Kennel
Energies 2020, 13(2), 492; https://doi.org/10.3390/en13020492 - 19 Jan 2020
Cited by 30 | Viewed by 3643
Abstract
In this paper, a simplified efficient method for sensorless finite set current predictive control (FSCPC) for synchronous reluctance motor (SynRM) based on extended Kalman filter (EKF) is proposed. The proposed FSCPC is based on reducing the computation burden of the conventional FSCPC by [...] Read more.
In this paper, a simplified efficient method for sensorless finite set current predictive control (FSCPC) for synchronous reluctance motor (SynRM) based on extended Kalman filter (EKF) is proposed. The proposed FSCPC is based on reducing the computation burden of the conventional FSCPC by using the commanded reference currents to directly calculate the reference voltage vector (RVV). Therefore, the cost function is calculated for only three times and the necessity to test all possible voltage vectors will be avoided. For sensorless control, EKF is composed to estimate the position and speed of the rotor. Whereas the performance of the proposed FSCPC essentially necessitates the full knowledge of SynRM parameters and provides an insufficient response under the parameter mismatch between the controller and the motor, online parameter estimation based on EKF is combined in the proposed control strategy to estimate all parameters of the machine. Furthermore, for simplicity, the parameters of PI speed controller and initial values of EKF covariance matrices are tuned offline using Particle Swarm Optimization (PSO). To demonstrate the feasibility of the proposed control, it is implemented in MATLAB/Simulink and tested under different operating conditions. Simulation results show high robustness and reliability of the proposed drive. Full article
(This article belongs to the Special Issue Electric Machines and Drive Systems for Emerging Applications)
Show Figures

Figure 1

Figure 1
<p>(<b>a</b>) Coordinates of SynRM. (<b>b</b>) Flow chart of EKF.</p>
Full article ">Figure 2
<p>(<b>a</b>) Two-level voltage source inverter. (<b>b</b>) Voltage vectors.</p>
Full article ">Figure 3
<p>Conventional finite set current predictive control (FSCPC) of SynRM.</p>
Full article ">Figure 4
<p>Proposed FSCPC of SynRM.</p>
Full article ">Figure 5
<p>Simulation results of the estimated and actual variables of SynRM baesd on EKF for step changes in speed (from top): speed, speed error, position, position error, stator resistance, q-axis inductance, d-axis inductance, load torque, stator d-axis current, and stator q-axis current.</p>
Full article ">Figure 6
<p>Simulation results of the estimated and actual variables of SynRM baesd on EKF for low speeds and reversal speed (from top): speed, speed error, position, position error, stator resistance, q-axis inductance, d-axis inductance, load torque, stator d-axis current, and stator q-axis current.</p>
Full article ">Figure 7
<p>Simulation results of the proposed EKF at 25% and 50% change in the SynRM stator resistance <math display="inline"> <semantics> <msub> <mi>R</mi> <mi>s</mi> </msub> </semantics> </math> (from top to bottom): speed (<math display="inline"> <semantics> <mrow> <msub> <mi>ω</mi> <mrow> <mi>r</mi> <mo>,</mo> <mi>r</mi> <mi>e</mi> <mi>f</mi> </mrow> </msub> <mo>,</mo> <msub> <mi>ω</mi> <mi>r</mi> </msub> <mo>,</mo> <msub> <mover> <mi>ω</mi> <mo stretchy="false">^</mo> </mover> <mi>r</mi> </msub> </mrow> </semantics> </math>), resistance (<math display="inline"> <semantics> <mrow> <msub> <mi>R</mi> <mi>s</mi> </msub> <mo>,</mo> <msub> <mover accent="false"> <mi>R</mi> <mo stretchy="false">^</mo> </mover> <mi>s</mi> </msub> </mrow> </semantics> </math>), stator inductances (<math display="inline"> <semantics> <mrow> <msub> <mi>L</mi> <mrow> <mi>d</mi> <mi>q</mi> </mrow> </msub> <mo>,</mo> <msub> <mover> <mi>L</mi> <mo stretchy="false">^</mo> </mover> <mrow> <mi>d</mi> <mi>q</mi> </mrow> </msub> </mrow> </semantics> </math>), load torque (<math display="inline"> <semantics> <mrow> <msub> <mi>T</mi> <mi>L</mi> </msub> <mo>,</mo> <msub> <mover> <mi>T</mi> <mo stretchy="false">^</mo> </mover> <mi>L</mi> </msub> </mrow> </semantics> </math>), and estimated stator currents.</p>
Full article ">Figure 8
<p>Simulation results of the proposed EKF at 25% and 50% change in the SynRM stator inductances <math display="inline"> <semantics> <msub> <mi>L</mi> <mrow> <mi>d</mi> <mi>q</mi> </mrow> </msub> </semantics> </math> (from top): speed (<math display="inline"> <semantics> <mrow> <msub> <mi>ω</mi> <mrow> <mi>r</mi> <mo>,</mo> <mi>r</mi> <mi>e</mi> <mi>f</mi> </mrow> </msub> <mo>,</mo> <msub> <mi>ω</mi> <mi>r</mi> </msub> <mo>,</mo> <msub> <mover> <mi>ω</mi> <mo stretchy="false">^</mo> </mover> <mi>r</mi> </msub> </mrow> </semantics> </math>), stator inductances (<math display="inline"> <semantics> <mrow> <msub> <mi>L</mi> <mrow> <mi>d</mi> <mi>q</mi> </mrow> </msub> <mo>,</mo> <msub> <mover> <mi>L</mi> <mo stretchy="false">^</mo> </mover> <mrow> <mi>d</mi> <mi>q</mi> </mrow> </msub> </mrow> </semantics> </math>), resistance (<math display="inline"> <semantics> <mrow> <msub> <mi>R</mi> <mi>s</mi> </msub> <mo>,</mo> <msub> <mover> <mi>R</mi> <mo stretchy="false">^</mo> </mover> <mi>s</mi> </msub> </mrow> </semantics> </math>), load torque (<math display="inline"> <semantics> <mrow> <msub> <mi>T</mi> <mi>L</mi> </msub> <mo>,</mo> <msub> <mover> <mi>T</mi> <mo stretchy="false">^</mo> </mover> <mi>L</mi> </msub> </mrow> </semantics> </math>), and estimated stator currents.</p>
Full article ">Figure 9
<p>Simulation results of the proposed EKF at change in the load torque to 1.5 N.m (from top): speed (<math display="inline"> <semantics> <mrow> <msub> <mi>ω</mi> <mrow> <mi>r</mi> <mo>,</mo> <mi>r</mi> <mi>e</mi> <mi>f</mi> </mrow> </msub> <mo>,</mo> <msub> <mi>ω</mi> <mi>r</mi> </msub> <mo>,</mo> <msub> <mover> <mi>ω</mi> <mo stretchy="false">^</mo> </mover> <mi>r</mi> </msub> </mrow> </semantics> </math>), load torque (<math display="inline"> <semantics> <mrow> <msub> <mi>T</mi> <mi>L</mi> </msub> <mo>,</mo> <msub> <mover> <mi>T</mi> <mo stretchy="false">^</mo> </mover> <mi>L</mi> </msub> </mrow> </semantics> </math>), estimated quadrature current, resistance (<math display="inline"> <semantics> <mrow> <msub> <mi>R</mi> <mi>s</mi> </msub> <mo>,</mo> <msub> <mover> <mi>R</mi> <mo stretchy="false">^</mo> </mover> <mi>s</mi> </msub> </mrow> </semantics> </math>), and stator inductances (<math display="inline"> <semantics> <mrow> <msub> <mi>L</mi> <mrow> <mi>d</mi> <mi>q</mi> </mrow> </msub> <mo>,</mo> <msub> <mover> <mi>L</mi> <mo stretchy="false">^</mo> </mover> <mrow> <mi>d</mi> <mi>q</mi> </mrow> </msub> </mrow> </semantics> </math>).</p>
Full article ">Figure 10
<p>Simulation results for step change in the SynRM parameters of the proposed FSCPC and conventional FSCPC: (<b>a</b>) step change in the stator resistance <math display="inline"> <semantics> <msub> <mi>R</mi> <mi>s</mi> </msub> </semantics> </math> and (<b>b</b>) step change in the stator inductances <math display="inline"> <semantics> <msub> <mi>L</mi> <mrow> <mi>d</mi> <mi>q</mi> </mrow> </msub> </semantics> </math>.</p>
Full article ">Figure 11
<p>Simulation results for the comparison between proposed FSCPC and conventional FSCPC at rated speed 1500 rpm and rated torque 1 N.m.</p>
Full article ">
13 pages, 1092 KiB  
Article
Enhanced Biogas Production of Cassava Wastewater Using Zeolite and Biochar Additives and Manure Co-Digestion
by Chibueze G. Achi, Amro Hassanein and Stephanie Lansing
Energies 2020, 13(2), 491; https://doi.org/10.3390/en13020491 - 19 Jan 2020
Cited by 26 | Viewed by 7126
Abstract
Currently, there are challenges with proper disposal of cassava processing wastewater, and a need for sustainable energy in the cassava industry. This study investigated the impact of co-digestion of cassava wastewater (CW) with livestock manure (poultry litter (PL) and dairy manure (DM)), and [...] Read more.
Currently, there are challenges with proper disposal of cassava processing wastewater, and a need for sustainable energy in the cassava industry. This study investigated the impact of co-digestion of cassava wastewater (CW) with livestock manure (poultry litter (PL) and dairy manure (DM)), and porous adsorbents (biochar (B-Char) and zeolite (ZEO)) on energy production and treatment efficiency. Batch anaerobic digestion experiments were conducted, with 16 treatments of CW combined with manure and/or porous adsorbents using triplicate reactors for 48 days. The results showed that CW combined with ZEO (3 g/g total solids (TS)) produced the highest cumulative CH4 (653 mL CH4/g VS), while CW:PL (1:1) produced the most CH4 on a mass basis (17.9 mL CH4/g substrate). The largest reduction in lag phase was observed in the mixture containing CW (1:1), PL (1:1), and B-Char (3 g/g TS), yielding 400 mL CH4/g volatile solids (VS) after 15 days of digestion, which was 84.8% of the total cumulative CH4 from the 48-day trial. Co-digesting CW with ZEO, B-Char, or PL provided the necessary buffer needed for digestion of CW, which improved the process stability and resulted in a significant reduction in chemical oxygen demand (COD). Co-digestion could provide a sustainable strategy for treating and valorizing CW. Scale-up calculations showed that a CW input of 1000–2000 L/d co-digested with PL (1:1) could produce 9403 m3 CH4/yr using a 50 m3 digester, equivalent to 373,327 MJ/yr or 24.9 tons of firewood/year. This system would have a profit of $5642/yr and a $47,805 net present value. Full article
(This article belongs to the Special Issue Biogas for Rural Areas)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Cumulative CH<sub>4</sub> production based on volatile solids (VS) added to each reactor for cassava wastewater (CW) digested alone and co-digested with poultry litter (PL) shown in (<b>A</b>), co-digestion with dairy manure (DM) shown in (<b>B</b>), and co-digestion with zeolite (ZEO) and biochar (B-Char) at low and high concentrations (LC and HC) shown in (<b>C</b>) and (<b>D</b>), respectively.</p>
Full article ">Figure 2
<p>Cumulative methane (CH<sub>4</sub>) production based on volatile solids (VS) added to each reactor for cassava wastewater (CW) digested alone and co-digested with poultry litter (PL), dairy manure (DM), zeolite (ZEO), and biochar (B-Char) at low and high concentrations (LC and HC) at five time points (Days 9, 15, 20, 37, and 48) in the 48-day digestion period.</p>
Full article ">Figure 3
<p>Scale-up model for a cassava processing factory, with two 25 m<sup>3</sup> digesters plumbed in series to treat cassava wastewater and poultry manure with the utilization of the digestate for fertilizer.</p>
Full article ">
24 pages, 1485 KiB  
Article
Performance Assessment of Japanese Electric Power Industry: DEA Measurement with Future Impreciseness
by Toshiyuki Sueyoshi and Mika Goto
Energies 2020, 13(2), 490; https://doi.org/10.3390/en13020490 - 19 Jan 2020
Cited by 7 | Viewed by 3569
Abstract
This study examines the performance of Japanese electric power companies from 2003 to 2020. We use an observed data set from 2003 to 2015 and a forecasted data set from 2016 to 2020. The Japanese deregulation of the industry needs to be completed [...] Read more.
This study examines the performance of Japanese electric power companies from 2003 to 2020. We use an observed data set from 2003 to 2015 and a forecasted data set from 2016 to 2020. The Japanese deregulation of the industry needs to be completed by April 2020. As a method, this study uses data envelopment analysis (DEA) environmental assessment, which measures performance from a holistic perspective. This research adds a new analytical capability to the DEA-based assessment by including an analytical ability to handle an “imprecise” data set. We apply the proposed approach to investigate the performance of these companies before and after the disaster of Fukushima Daiichi nuclear power plant (11 March 2011). All electric power companies have suffered from business damage due to the nuclear disaster. The Japanese government has developed a policy scheme on how to recover from the huge handling costs resulting from the disaster. Nuclear energy has been long considered the most useful approach to handle climate change. However, many industrial nations have changed policy direction since the nuclear disaster. The Japanese government allocates the costs to not only Tokyo Electric Power Company, which produced the nuclear disaster, but also the other incumbent electric power companies that own nuclear power plants. Under the current Japanese scheme, financial conditions have been gradually recovering from the damage due to the managerial efforts and by indirectly allocating the expenditure to consumers and tax payers. Full article
(This article belongs to the Special Issue Data Envelopment Analysis (DEA) Applied to Energy and Environment)
Show Figures

Figure 1

Figure 1
<p>Range for supporting line on <span class="html-italic">x</span> and <span class="html-italic">g.</span></p>
Full article ">Figure 2
<p>Methodological comparison and forecasting.</p>
Full article ">Figure 3
<p>Ten electric power companies by service area. (Source: Electricity Review Japan, The Federation of Electric Power Companies of Japan 2012.)</p>
Full article ">Figure 4
<p>Annual shift of unified efficiencies of electric power companies.</p>
Full article ">Figure 5
<p>Annual shift of unified efficiencies of electric power companies: Pooled data from 2016 to 2020.</p>
Full article ">
13 pages, 4554 KiB  
Article
Electric Vehicle Battery Performance Investigation Based on Real World Current Harmonics
by Sid-Ali Amamra, Yashraj Tripathy, Anup Barai, Andrew D. Moore and James Marco
Energies 2020, 13(2), 489; https://doi.org/10.3390/en13020489 - 19 Jan 2020
Cited by 15 | Viewed by 4564
Abstract
Electric vehicle (EV) powertrains consist of power electronic components as well as electric machines to manage the energy flow between different powertrain subsystems and to deliver the necessary torque and power requirements at the wheels. These power subsystems can generate undesired electrical harmonics [...] Read more.
Electric vehicle (EV) powertrains consist of power electronic components as well as electric machines to manage the energy flow between different powertrain subsystems and to deliver the necessary torque and power requirements at the wheels. These power subsystems can generate undesired electrical harmonics on the direct current (DC) bus of the powertrain. This may lead to the on-board battery being subjected to DC current superposed with undesirable high- and low- frequency current oscillations, known as ripples. From real-world measurements, significant current harmonics perturbations within the range of 50 Hz to 4 kHz have been observed on the high voltage DC bus of the EV. In the limited literature, investigations into the impact of these harmonics on the degradation of battery systems have been conducted. In these studies, the battery systems were supplied by superposed current signals i.e., DC superposed by a single frequency alternating current (AC). None of these studies considered applying the entire spectrum of the ripple current measured in the real-world scenario, which is focused on in this research. The preliminary results indicate that there is no difference concerning capacity fade or impedance rise between the cells subjected to just DC current and those subjected additionally to a superposed AC ripple current. Full article
(This article belongs to the Section E: Electric Vehicles)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>BEV powertrain architecture.</p>
Full article ">Figure 2
<p>Real-world battery current measurement within 10 kHz resolution and current spectrum.</p>
Full article ">Figure 3
<p>Example test signal as described in <a href="#sec3dot2-energies-13-00489" class="html-sec">Section 3.2</a>, showing a DC signal with a superimposed AC cycle (for illustrative purposes).</p>
Full article ">Figure 4
<p>Experimental test rig description.</p>
Full article ">Figure 5
<p>Experimental procedure.</p>
Full article ">Figure 6
<p>Normalized capacity fade as a function of number of cycles for normal condition group (i.e., cells 1–3) versus AC ripple condition group (i.e., cells 3–5) and retained energy versus total energy throughput.</p>
Full article ">Figure 7
<p>The difference in temperature during discharge operation for all cells at fresh condition (<b>a</b>) and after 250 cycles (<b>b</b>).</p>
Full article ">Figure 8
<p>DC impedance increase as a function of number of cycles and total energy throughput for normal condition set group (i.e., cells 1–3) versus AC ripple condition set group (i.e., cells 4–6).</p>
Full article ">Figure 9
<p>Electrochemical impedance spectroscopy (EIS) plots, for (<b>a</b>) all fresh cells at 80% SoC, (<b>b</b>) Cell 1 and cell 5 at 80% SoC—fresh cells, (<b>c</b>) change of impedance spectrum with degradation after 250 cycles of cell1 and cell 5 at 80% SoC, (<b>d</b>) all fresh cells at 20% SoC, (<b>e</b>) Cell 1 and cell 5 at 20% SoC—fresh cells, (<b>f</b>) change of impedance spectrum with degradation after 250 cycles of cell 1 and cell 5 at 20% SoC.</p>
Full article ">Figure 10
<p>Bode plot shows impedance vs. frequency and phase vs. frequency for cells 1 (<b>black</b>) and 5 (<b>red</b>) at 80% SOC, for (<b>a</b>) fresh cells, and (<b>b</b>) after 250 cycles were performed.</p>
Full article ">
24 pages, 3534 KiB  
Article
Assessing the Techno-Economics and Environmental Attributes of Utility-Scale PV with Battery Energy Storage Systems (PVS) Compared to Conventional Gas Peakers for Providing Firm Capacity in California
by Sashwat Roy, Parikhit Sinha and Syed Ismat Shah
Energies 2020, 13(2), 488; https://doi.org/10.3390/en13020488 - 19 Jan 2020
Cited by 10 | Viewed by 5850
Abstract
The United States needs to add at least 20 GW of peaking capacity to its grid over the next 10 years, led by large-scale projects in California, Texas and Arizona. Of that, about 60% must be installed between 2023 and 2027, meaning that [...] Read more.
The United States needs to add at least 20 GW of peaking capacity to its grid over the next 10 years, led by large-scale projects in California, Texas and Arizona. Of that, about 60% must be installed between 2023 and 2027, meaning that the energy storage industry has more time to build an economic advantage by lowering costs and improving performance to compete with conventional gas peakers. In this paper, we assess the technical feasibility of utility-scale PV plus battery energy storage (PVS) to provide high capacity factors during summer peak demand periods using a target period capacity factor (TPCF) framework as an alternative to natural gas peakers. Also, a new metric called “Lifetime Cost of Operation” (LCOO) is introduced to provide a metric, focusing on the raw installation and operational costs of PVS technology compared to natural-gas fired peaker plants (simple cycle or conventional combustion turbine) during the target period window. The target period window is the time period during which it is valuable for power plants to provide firm capacity usually during early or late evening peak demand periods in the summer months (from April to September); a framework for which is increasingly being asked for by utilities in recent request for proposals (RFPs). A 50 MWAC PV system with 60 MW/240 MWh battery storage modelled in California can provide >98% capacity factor over a 7–10 p.m. target period with lower LCOO than a conventional combustion turbine natural gas power plant. Full article
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Electricity consumption in California for 2017 (residential and non-residential sectors) by county. Source: California Energy Commission (CEC, 2018) staff. <a href="http://ecdms.energy.ca.gov/" target="_blank">http://ecdms.energy.ca.gov/</a>.</p>
Full article ">Figure 2
<p>California’s peak to average electricity demand ratio has been generally rising in the past decade, particularly in Southern California, although it is highly dependent on yearly climate, housing stock and other economic conditions. Source: Energy Information Administration (EIA. 2016).</p>
Full article ">Figure 3
<p>An illustration of a typical 12 × 24 matrix provided by the utility representing the normalized value of energy during different times over the year with the highlighted green representing the target period window. Source: First Solar PVS Analytics.</p>
Full article ">Figure 4
<p>An illustration of Levelized Cost of Electricity (LCOE) vs. target period capacity factor (TPCF) showing the different system sizes (represented by dots) which can meet the TPCF requirement. Source: First Solar PVS Analytics.</p>
Full article ">Figure 5
<p>An example of LCOE vs TPCF portraying different technologies, system sizes and configurations which meet the TPCF requirement of 90% which is specified in the request for proposals (RFP). Source: First Solar PVS Analytics.</p>
Full article ">Figure 6
<p>(<b>a</b>) “Duck Curve” in CAISO and (<b>b</b>) Target Period window selected for analyzing the ability of PVS to provide high capacity factors during this duration. Source: CAISO, 2015.</p>
Full article ">Figure 6 Cont.
<p>(<b>a</b>) “Duck Curve” in CAISO and (<b>b</b>) Target Period window selected for analyzing the ability of PVS to provide high capacity factors during this duration. Source: CAISO, 2015.</p>
Full article ">Figure 7
<p>Single Line Diagram for AC-Coupled PVS showing the different nodes for power calculation purposes. Source: First Solar Plant Predict Resource Center.</p>
Full article ">Figure 8
<p>A typical profile for PV and Battery Energy Storage System (BESS) net energy for a typical summer day (27 June) in San Bernardino, California.</p>
Full article ">Figure 9
<p>PVS Net Energy Delivered in Year-1 (in MW<sub>DC</sub>) Heat Map for the 50 MW<sub>AC</sub> PV-60 MW/240 MWh BESS system for target period 7–10 p.m. in (<b>a</b>) San Bernardino, CA (<b>b</b>) San Luis Obispo, CA and (<b>c</b>) Sonoma, CA.</p>
Full article ">Figure 9 Cont.
<p>PVS Net Energy Delivered in Year-1 (in MW<sub>DC</sub>) Heat Map for the 50 MW<sub>AC</sub> PV-60 MW/240 MWh BESS system for target period 7–10 p.m. in (<b>a</b>) San Bernardino, CA (<b>b</b>) San Luis Obispo, CA and (<b>c</b>) Sonoma, CA.</p>
Full article ">Figure 10
<p>PVS Net Energy Delivered (MW<sub>DC</sub>) Heat Map for the 50 MWac PV-60 MW/240 MWh BESS system in Sonoma, CA for target period 6–10 p.m.</p>
Full article ">Figure 11
<p>Breakdown of Engineering Procurement &amp; Construction (EPC) Costs for DC (<b>top</b>) and AC Coupled (<b>bottom</b>) PVS Systems.</p>
Full article ">Figure 12
<p>Lifetime cost of operation (LCOO) comparison between PVS (a 50 MW<sub>AC</sub> PV with 60 MW/240 MWh BESS) and Conventional CT (70 MW) power plants in M<span>$</span> over twenty year system life without and with considering yearly environmental cost of operation due to emissions.</p>
Full article ">Figure 13
<p>Performance Ratio of the PV plant declines during the target period window leading to a proportional loss of PV production capability during the summer months (April–September).</p>
Full article ">Figure 14
<p>Future LCOO comparison between PVS and Combustion Turbine (CT) based on Li-Ion BESS cost reductions, natural gas price reduction and ITC ramp down.</p>
Full article ">
19 pages, 2721 KiB  
Article
Optimal Design of a Multibrid Permanent Magnet Generator for a Tidal Stream Turbine
by Khalil Touimi, Mohamed Benbouzid and Zhe Chen
Energies 2020, 13(2), 487; https://doi.org/10.3390/en13020487 - 19 Jan 2020
Cited by 7 | Viewed by 3780
Abstract
Tidal stream energy is acquiring more attention as a future potential renewable energy source. Considering the harsh submarine environment, the main challenges that face the tidal stream turbine (TST) industry are cost and reliability. Hence, simple and reliable technologies, especially considering the drivetrain, [...] Read more.
Tidal stream energy is acquiring more attention as a future potential renewable energy source. Considering the harsh submarine environment, the main challenges that face the tidal stream turbine (TST) industry are cost and reliability. Hence, simple and reliable technologies, especially considering the drivetrain, are preferred. The multibrid drivetrain configuration with only a single stage gearbox is one of the promising concepts for TST systems. In this context, this paper proposes the design optimization of a multibrid permanent magnet generator (PMG), the design of a planetary gearbox, and afterwards analyzes the multibrid concept cost-effectiveness for TST applications. Firstly, the system analytical model, which consists of a single-stage gearbox and a medium speed PMG, is presented. The optimization methodology is afterwards highlighted. Lastly, the multibrid system optimization results for different gear ratios including the direct-drive topology are discussed and compared where the suitable gear ratio (topology) is investigated. The achieved results show that the multibrid concept in TST applications seems more attractive than the direct-drive one especially for high power ratings. Full article
Show Figures

Figure 1

Figure 1
<p>OpenHydro/Naval Energies direct-drive tidal stream turbine [<a href="#B4-energies-13-00487" class="html-bibr">4</a>].</p>
Full article ">Figure 2
<p>The multibrid concept: (<b>a</b>) Schematic illustration [<a href="#B12-energies-13-00487" class="html-bibr">12</a>], (<b>b</b>) The AREVA Multibrid M5000 5 MW wind turbine nacelle [<a href="#B11-energies-13-00487" class="html-bibr">11</a>].</p>
Full article ">Figure 3
<p>Small scale multibrid tidal stream turbine [<a href="#B15-energies-13-00487" class="html-bibr">15</a>].</p>
Full article ">Figure 4
<p>Scheme of a grid-connected single stage permanent magnet generator-based tidal stream turbine.</p>
Full article ">Figure 5
<p>Tidal velocity in polar coordinates.</p>
Full article ">Figure 6
<p>Tidal velocity in the Ouessant Island.</p>
Full article ">Figure 7
<p>Tidal current energy distribution.</p>
Full article ">Figure 8
<p>Tidal current amplitude speed distribution.</p>
Full article ">Figure 9
<p>Harnessed energy rate versus power limitation rate.</p>
Full article ">Figure 10
<p>Illustration of a planetary gearbox with five planet gears.</p>
Full article ">Figure 11
<p>Basic dimensions of one pair of poles [<a href="#B8-energies-13-00487" class="html-bibr">8</a>].</p>
Full article ">Figure 12
<p>Flowchart describing the design optimization procedure.</p>
Full article ">Figure 13
<p>Tidal stream turbine (TST) estimated cost for different gear ratios and different power ratings.</p>
Full article ">Figure 14
<p>TST total estimated cost per kWh.</p>
Full article ">Figure 15
<p>Gearbox and generator costs for different gear ratios at the power rating of 1.5 MW.</p>
Full article ">Figure 16
<p>Generator active materials weight for different gear ratios at the power rating of 1.5 MW.</p>
Full article ">Figure 17
<p>Generator active materials cost for different gear ratios at the power rating of 1.5 MW.</p>
Full article ">Figure 18
<p>Active material geared generator cost compared to direct drive (DD) generator.</p>
Full article ">Figure 19
<p>View of the designed (3:1) geared generator at the power rating of 1.5 MW.</p>
Full article ">Figure 20
<p>Front and lateral view of the designed (3:1) geared generator at the power rating of 1.5 MW.</p>
Full article ">
19 pages, 13758 KiB  
Article
Prediction of Reservoir Quality from Log-Core and Seismic Inversion Analysis with an Artificial Neural Network: A Case Study from the Sawan Gas Field, Pakistan
by Zhang Qiang, Qamar Yasin, Naser Golsanami and Qizhen Du
Energies 2020, 13(2), 486; https://doi.org/10.3390/en13020486 - 19 Jan 2020
Cited by 47 | Viewed by 5655
Abstract
This paper presents a novel approach that aims to predict better reservoir quality regions from seismic inversion and spatial distribution of key reservoir properties from well logs. The reliable estimation of lithology and reservoir parameters at sparsely located wells in the Sawan gas [...] Read more.
This paper presents a novel approach that aims to predict better reservoir quality regions from seismic inversion and spatial distribution of key reservoir properties from well logs. The reliable estimation of lithology and reservoir parameters at sparsely located wells in the Sawan gas field is still a considerable challenge. This is due to three main reasons: (a) the extreme heterogeneity in the depositional environments, (b) sand-shale intercalations, and (c) repetition of textural changes from fine to coarse sandstone and very coarse sandstone in the reservoir units. In this particular study, machine learning (ML) inversion algorithm was selected to predict the spatial variations of acoustic impedance (AI), porosity, and saturation. While trained in a supervised mode, the support vector machine (SVM) inversion algorithm performed effectively in identifying and mapping individual reservoir properties to delineate and quantify fluid-rich zones. Meanwhile, the Sequential Gaussian Simulation (SGS) and Gaussian Indicator Simulation (GIS) algorithms were employed to determine the spatial variability of lithofacies and porosity from well logs and core analyses data. The calibration of the detailed spatial variations from post-stack seismic inversion using SVM and wireline logs data indicated an appropriate agreement, i.e., variations in AI is related to the variations in reservoir facies and parameters. From the current study, it was concluded that in a highly heterogeneous reservoir, the integration of SVM and GIS algorithms is a reliable approach to achieve the best estimation of the spatial distribution of detailed reservoir characteristics. The results obtained in this study would also be helpful to minimize the uncertainty in drilling, production, and injection in the Sawan gas field of Pakistan as well as other reservoirs worldwide with similar geological settings. Full article
(This article belongs to the Section H: Geo-Energy)
Show Figures

Figure 1

Figure 1
<p>(<b>a</b>) Structural map of Pakistan with the location of the Sawan Gas Field, (<b>b</b>) Generalized stratigraphic column of the Lower Indus Basin with highlighted C-sand interval.</p>
Full article ">Figure 2
<p>(<b>a</b>) Seismic base map with wells and Xline 932 (orange color), (<b>b</b>) Seismic to well tie at Xline 932.</p>
Full article ">Figure 3
<p>The diagrammatic representation of SVM for separation hyper-plane in 2D [<a href="#B22-energies-13-00486" class="html-bibr">22</a>].</p>
Full article ">Figure 4
<p>Estmated petrophysical properties of reservoir C-sand interval in the study wells; estimated values are shown in continuous lines while laboratory measurement results are indicated with red dots.</p>
Full article ">Figure 5
<p>Thin sections (<b>a</b>,<b>b</b>) and SEM images (<b>c</b>,<b>d</b>) of the C-sand interval demonstrate the depositional and diagenetic characteristics, where Chl = chlorite, Chlr = chlorite rim, Q = quartz, Fsp = dissolving of feldspar grains, VRF = alteration of volcanic rock fragments, Qo = Idiomorphic quartz outgrowths filling partially the pores after Chlr, Gc = Grain contact. P = porosity is 15.3 to 22% in wells C and D [<a href="#B24-energies-13-00486" class="html-bibr">24</a>].</p>
Full article ">Figure 6
<p>Typical T<sub>2</sub> relaxation curves observed for the core sample from the study area show various distributions of micro to macropores which result in the complex nature of the reservoir interval of the Sawan gas field, (<b>a</b>) sample S1, (<b>b</b>) sample S2, (<b>c</b>) sample S3, and (<b>d</b>) sample S4.</p>
Full article ">Figure 7
<p>Structural grid of the studied reservoir; (<b>a</b>) Low-resolution grid cells (71,760 cells), (<b>b</b>) High-resolution grid cells (72,649,401 cells).</p>
Full article ">Figure 8
<p>Comparison of facies analysis; Track B: Gamma-ray log, Track C: Facies analysis from well logs, Track D: Facies analysis from low-resolution up-scaling, Track F: Facies analysis from high-resolution up-scaling, and Track E: Potential reservoir zones.</p>
Full article ">Figure 9
<p>The petrophysical modeling and inversion of facies (<b>a</b>) 3D facies model, (<b>b</b>) cross-sectional view of the 3D facies model.</p>
Full article ">Figure 10
<p>Synthetic seismogram for Well-C, showing (left to right): TVD (true vertical depth) in meters, sonic transit time (µs/ft) and bulk density from logs (g/cm<sup>3</sup>), RC (reflection coefficient), synthetic seismogram, traces from a portion of seismic Xline 932, and TWT (two-way time in ms).</p>
Full article ">Figure 11
<p>The plots of the extracted wavelet (Time-Amplitude, Frequency-Power, Frequency-Phase) for correlation between seismic and well log in time and frequency at 2000–2300 ms.</p>
Full article ">Figure 12
<p>The inverted AI for Xline 932 using SVM seismic inversion. The impedance log of well-A, C, and D are not in good match with the inverted impedance surface.</p>
Full article ">Figure 13
<p>Comparison of original logs (blue) and SVM-inverted results (red) at Well-A, Well-C, and Well-D by SVM inversion method.</p>
Full article ">Figure 14
<p>A cross-plot of acoustic impedance (AI) and porosity with a color bar representing hydrocarbon saturation values (Sat. HC).</p>
Full article ">Figure 15
<p>The seismic inversion map of effective porosity for C-sand interval (2140 ms to 2280 ms).</p>
Full article ">Figure 16
<p>(<b>a</b>) 3D AI model, (<b>b</b>) seismic inversion of 3D effective porosity models within a specified time window (2160 to 2280 ms) covering target horizons.</p>
Full article ">Figure 17
<p>The cross-sectional analysis of 3D effective porosity model by well log simulation of C-sand interval.</p>
Full article ">Figure 18
<p>The time slice of porosity at various interval, (<b>a</b>) 2160 ms, (<b>b</b>) 2180 ms, (<b>c</b>) 2200 ms, and (<b>d</b>) 2220 ms.</p>
Full article ">Figure 19
<p>The time slices of (<b>a</b>) AI, (<b>b</b>) hydrocarbon saturation (Sat. HC) at Z = 2200, demonstrating the prospects within the target reservoir interval.</p>
Full article ">Figure 20
<p>3D Seismic attributes map, (<b>a</b>) low-frequency attenuation gradient, (<b>b</b>) amplitude above average at Z = 2200 ms, demonstrating the high amplitude (maroon) is being consistent with low-frequency (black).</p>
Full article ">
10 pages, 2650 KiB  
Article
Computer Simulation of Stochastic Energy Fluctuations in Tensile Test of Elasto-Plastic Porous Metallic Material
by Marcin Kamiński and Michał Strąkowski
Energies 2020, 13(2), 485; https://doi.org/10.3390/en13020485 - 19 Jan 2020
Viewed by 2222
Abstract
The main aim of this work is the computational implementation and numerical simulation of a metal porous plasticity model with an uncertain initial microdefects’ volume fraction using the Stochastic Finite Element Method (SFEM) based on the semi-analytical probabilistic technique. The metal porous plasticity [...] Read more.
The main aim of this work is the computational implementation and numerical simulation of a metal porous plasticity model with an uncertain initial microdefects’ volume fraction using the Stochastic Finite Element Method (SFEM) based on the semi-analytical probabilistic technique. The metal porous plasticity model applied here is based on Gurson–Tvergaard–Needleman theory and is included in the ABAQUS finite element system, while the external probabilistic procedures were programmed in the computer algebra system MAPLE 2017. Hybrid usage of these two computer systems enabled the determination of fluctuations in elastic and plastic energies due to initial variations in the ratio of the metal micro-voids, and the calculation of the first four probabilistic moments and coefficients of these energies due to Gaussian distribution of this ratio. A comparison with the Monte-Carlo simulation validated the numerical efficiency of the proposed approach for any level of input uncertainty and for the first four probabilistic characteristics traditionally seen in the experimental series. Full article
(This article belongs to the Special Issue Recent Advances in Stochastic Methods for Energy Analysis)
Show Figures

Figure 1

Figure 1
<p>Geometry, boundary conditions and discretization of the specimen.</p>
Full article ">Figure 2
<p>Deformed material specimen.</p>
Full article ">Figure 3
<p>Elastic energy variations during deformation.</p>
Full article ">Figure 4
<p>Plastic energy variations during deformation.</p>
Full article ">Figure 5
<p>Expected values of elastic energy fluctuations.</p>
Full article ">Figure 6
<p>Expected values of plastic energy fluctuations.</p>
Full article ">Figure 7
<p>Coefficients of variation of elastic energy fluctuations.</p>
Full article ">Figure 8
<p>Coefficients of variation of elastic energy fluctuations.</p>
Full article ">Figure 9
<p>Skewness of elastic energy fluctuations.</p>
Full article ">Figure 10
<p>Skewness of plastic energy fluctuations.</p>
Full article ">Figure 11
<p>Kurtosis of elastic energy fluctuations.</p>
Full article ">Figure 12
<p>Kurtosis of plastic energy fluctuations.</p>
Full article ">
19 pages, 2725 KiB  
Article
Electrical Insulator Fault Forecasting Based on a Wavelet Neuro-Fuzzy System
by Stéfano Frizzo Stefenon, Roberto Zanetti Freire, Leandro dos Santos Coelho, Luiz Henrique Meyer, Rafael Bartnik Grebogi, William Gouvêa Buratto and Ademir Nied
Energies 2020, 13(2), 484; https://doi.org/10.3390/en13020484 - 19 Jan 2020
Cited by 54 | Viewed by 3957
Abstract
The surface contamination of electrical insulators can increase the electrical conductivity of these components, which may lead to faults in the electrical power system. During inspections, ultrasound equipment is employed to detect defective insulators or those that may cause failures within a certain [...] Read more.
The surface contamination of electrical insulators can increase the electrical conductivity of these components, which may lead to faults in the electrical power system. During inspections, ultrasound equipment is employed to detect defective insulators or those that may cause failures within a certain period. Assuming that the signal collected by the ultrasound device can be processed and used for both the detection of defective insulators and prediction of failures, this study starts by presenting an experimental procedure considering a contaminated insulator removed from the distribution line for data acquisition. Based on the obtained data set, an offline time series forecasting approach with an Adaptive Neuro-Fuzzy Inference System (ANFIS) was conducted. To improve the time series forecasting performance and to reduce the noise, Wavelet Packets Transform (WPT) was associated to the ANFIS model. Once the ANFIS model associated with WPT has distinct parameters to be adjusted, a complete evaluation concerning different model configurations was conducted. In this case, three inference system structures were evaluated: grid partition, fuzzy c-means clustering, and subtractive clustering. A performance analysis focusing on computational effort and the coefficient of determination provided additional parameter configurations for the model. Taking into account both parametrical and statistical analysis, the Wavelet Neuro-Fuzzy System with fuzzy c-means showed that it is possible to achieve impressive accuracy, even when compared to classical approaches, in the prediction of electrical insulators conditions. Full article
(This article belongs to the Special Issue Machine Learning for Energy Systems)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Contaminated insulator removed from a 25 kV rural area distribution system.</p>
Full article ">Figure 2
<p>Comparison between original and rebuild signal using Wavelet Packets Transform (WPT).</p>
Full article ">Figure 3
<p>Adaptive Neuro-Fuzzy Inference System (ANFIS) structure for time series evaluation.</p>
Full article ">Figure 4
<p>Method flowchart for insulator withdrawal and model evaluation.</p>
Full article ">Figure 5
<p>Comparison between predicted and real data assuming an FCM structure for the testing phase.</p>
Full article ">
24 pages, 8301 KiB  
Article
Expansion of High Efficiency Region of Wind Energy Centrifugal Pump Based on Factorial Experiment Design and Computational Fluid Dynamics
by Wei Li, Leilei Ji, Weidong Shi, Ling Zhou, Hao Chang and Ramesh K. Agarwal
Energies 2020, 13(2), 483; https://doi.org/10.3390/en13020483 - 19 Jan 2020
Cited by 13 | Viewed by 3006
Abstract
The wind energy pump system is a new green energy technology. The wide high efficiency region of pump is of great significance for energy conservation of wind power pumping system. In this study, factorial experiment design (FED) and computational fluid dynamics (CFD) are [...] Read more.
The wind energy pump system is a new green energy technology. The wide high efficiency region of pump is of great significance for energy conservation of wind power pumping system. In this study, factorial experiment design (FED) and computational fluid dynamics (CFD) are employed to optimize the hydraulic design of wind energy centrifugal pump (WECP). The blade outlet width b2, blade outlet angle β2, and blade wrap angle ψ are chosen as factors of FED. The effect of the factors on the efficiency under the conditions of 0.6Qdes, 0.8Qdes, 1.0Qdes, and 1.4Qdes is systematically analyzed. The matching feature of various volute tongue angle with the optimized impeller is also investigated numerically and experimentally. After the optimization, the pump head changes smoothly during full range of flow condition and the high efficiency region is effectively improved. The weighted average efficiency of four conditions increases by 2.55%, which broadens the high efficiency region of WECP globally. Besides, the highest efficiency point moves towards the large flow conditions. The research results provide references for expanding the efficient operation region of WECP. Full article
(This article belongs to the Special Issue Design, Fabrication and Performance of Wind Turbines 2020)
Show Figures

Figure 1

Figure 1
<p>Definition of key parameters.</p>
Full article ">Figure 2
<p>Hydraulic model of Impeller 1.</p>
Full article ">Figure 3
<p>3D model of the calculation domain.</p>
Full article ">Figure 4
<p>Structured mesh for the numerical domains.</p>
Full article ">Figure 5
<p>Pump head with different mesh number.</p>
Full article ">Figure 6
<p>Efficiency curves of impeller under different schemes.</p>
Full article ">Figure 7
<p>Head curves of impeller under different schemes.</p>
Full article ">Figure 8
<p>Pareto diagram of efficiency under different operation conditions.</p>
Full article ">Figure 9
<p>Main effects diagram of different operation conditions.</p>
Full article ">Figure 10
<p>Interaction effect diagram of different flow rate point: (<b>a</b>) Interaction effect diagram of 0.6<span class="html-italic">Q</span><sub>des</sub>. (<b>b</b>) Interaction effect diagram of 0.8<span class="html-italic">Q</span><sub>des</sub>. (<b>c</b>) Interaction effect diagram of 1.0<span class="html-italic">Q</span><sub>des</sub>. (<b>d</b>) Interaction effect diagram of 1.4<span class="html-italic">Q</span><sub>des</sub>.</p>
Full article ">Figure 10 Cont.
<p>Interaction effect diagram of different flow rate point: (<b>a</b>) Interaction effect diagram of 0.6<span class="html-italic">Q</span><sub>des</sub>. (<b>b</b>) Interaction effect diagram of 0.8<span class="html-italic">Q</span><sub>des</sub>. (<b>c</b>) Interaction effect diagram of 1.0<span class="html-italic">Q</span><sub>des</sub>. (<b>d</b>) Interaction effect diagram of 1.4<span class="html-italic">Q</span><sub>des</sub>.</p>
Full article ">Figure 11
<p>Performance curves of optimized impeller and original impeller.</p>
Full article ">Figure 12
<p>3D model of different volute.</p>
Full article ">Figure 13
<p>Performance curves with different volutes.</p>
Full article ">Figure 14
<p>Velocity distribution on volute cross section under design flow rate condition with different volutes.</p>
Full article ">Figure 15
<p>Pressure distribution on volute cross section under design flow rate condition with different volutes.</p>
Full article ">Figure 16
<p>The positions of the cross sections.</p>
Full article ">Figure 17
<p>Vector distribution of four radial cross sections under design flow rate condition.</p>
Full article ">Figure 17 Cont.
<p>Vector distribution of four radial cross sections under design flow rate condition.</p>
Full article ">Figure 18
<p>Distribution of turbulent kinetic energy under design condition with different tongue angles.</p>
Full article ">Figure 19
<p>Model.</p>
Full article ">Figure 20
<p>Test bed. 1. Pipeline export. 2. Water tank. 3. Pipeline. 4. Electromagnetic flowmeter. 5. Discharge valve. 6. Test pump.</p>
Full article ">Figure 21
<p>Performance curves of repetitive testing.</p>
Full article ">Figure 22
<p>Comparison between experiment and simulation.</p>
Full article ">
22 pages, 1354 KiB  
Review
Factors Affecting the Performance of Membrane Osmotic Processes for Bioenergy Development
by Wen Yi Chia, Kuan Shiong Khoo, Shir Reen Chia, Kit Wayne Chew, Guo Yong Yew, Yeek-Chia Ho, Pau Loke Show and Wei-Hsin Chen
Energies 2020, 13(2), 481; https://doi.org/10.3390/en13020481 - 19 Jan 2020
Cited by 10 | Viewed by 8331
Abstract
Forward osmosis (FO) and pressure-retarded osmosis (PRO) have gained attention recently as potential processes to solve water and energy scarcity problems with advantages over pressure-driven membrane processes. These processes can be designed to produce bioenergy and clean water at the same time (i.e., [...] Read more.
Forward osmosis (FO) and pressure-retarded osmosis (PRO) have gained attention recently as potential processes to solve water and energy scarcity problems with advantages over pressure-driven membrane processes. These processes can be designed to produce bioenergy and clean water at the same time (i.e., wastewater treatment with power generation). Despite having significant technological advancement, these bioenergy processes are yet to be implemented in full scale and commercialized due to its relatively low performance. Hence, massive and extensive research has been carried out to evaluate the variables in FO and PRO processes such as osmotic membrane, feed solutions, draw solutions, and operating conditions in order to maximize the outcomes, which include water flux and power density. However, these research findings have not been summarized and properly reviewed. The key parts of this review are to discuss the factors influencing the performance of FO and PRO with respective resulting effects and to determine the research gaps in their optimization with the aim of further improving these bioenergy processes and commercializing them in various industrial applications. Full article
(This article belongs to the Special Issue Biomass Processing for Biofuels, Bioenergy and Chemicals)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Cross-section SEM micrographs of (<b>A</b>) commercial CTA-HTI and (<b>B</b>) TFC membrane. Reproduced with permission from Yip, Tiraferri, Phillip, Schiffman and Elimelech [<a href="#B16-energies-13-00481" class="html-bibr">16</a>], Copyright (1969) American Chemical Society.</p>
Full article ">Figure 2
<p>Hollow fiber membranes.</p>
Full article ">Figure 3
<p>Classification of draw solutes/solution.</p>
Full article ">Figure 4
<p>Factors affecting osmotic membrane processes.</p>
Full article ">
39 pages, 6922 KiB  
Article
Measuring Impact of Uncertainty in a Stylized Macroeconomic Climate Model within a Dynamic Game Perspective
by Valentijn Stienen and Jacob Engwerda
Energies 2020, 13(2), 482; https://doi.org/10.3390/en13020482 - 18 Jan 2020
Cited by 1 | Viewed by 2376
Abstract
In this paper, we present a stylized dynamic interdependent multi-country energy transition model. The goal of this paper is to provide a starting point for examining the impact of uncertainty in such models. To do this, we define a simple model based on [...] Read more.
In this paper, we present a stylized dynamic interdependent multi-country energy transition model. The goal of this paper is to provide a starting point for examining the impact of uncertainty in such models. To do this, we define a simple model based on the standard Solow macroeconomic growth model. We consider this model in a two-country setting using a non-cooperative dynamic game perspective. Total carbon dioxide (CO2) emission is added in this growth model as a factor that has a negative impact on economic growth, whereas production can be realized using either green or fossil energy. Additionally, a factor is incorporated that captures the difficulties of using green energy, such as accessibility per country. We calibrate this model for a two-player setting, in which one player represents all countries affiliated with the Organization for Economic Cooperation and Development (OECD) and the other player represents countries not affiliated with the OECD. It is shown that, in general, the model is capable to describe energy transitions towards quite different equilibrium constellations. It turns out that this is mainly caused by the choice of policy parameters chosen in the objective function. We also analyze the optimal response strategies of both countries if the model in equilibrium would be hit by a CO2 shock. Also, here we observe a quite natural response. As the model is quite stylized, a serious study is performed to the impact several model uncertainties have on the results. It turns out that, within the OECD/non-OECD framework, most of the considered uncertainties do not impact results much. However, the way we calibrate policy parameters does carry much uncertainty and, as such, influences equilibrium outcomes a lot. Full article
(This article belongs to the Special Issue Assessment of Energy–Environment–Economy Interrelations)
Show Figures

Figure 1

Figure 1
<p>Climate facts.</p>
Full article ">Figure 2
<p>Percentage green energy with <math display="inline"><semantics> <mrow> <mi>π</mi> <mo>=</mo> <mi>α</mi> <mo>·</mo> <msup> <mi>π</mi> <mi>initial</mi> </msup> </mrow> </semantics></math> with <math display="inline"><semantics> <mrow> <mi>α</mi> <mo>∈</mo> <mo>[</mo> <mn>1</mn> <mo>,</mo> <mn>2.5</mn> <mo>]</mo> </mrow> </semantics></math>.</p>
Full article ">Figure 3
<p>Percentage green energy with <math display="inline"><semantics> <mrow> <msub> <mi>e</mi> <mi>i</mi> </msub> <mo>+</mo> <mi>M</mi> <msub> <mi>e</mi> <mi>j</mi> </msub> </mrow> </semantics></math> with <math display="inline"><semantics> <mrow> <mi>M</mi> <mo>∈</mo> <mo>[</mo> <mn>0</mn> <mo>,</mo> <mn>1</mn> <mo>]</mo> </mrow> </semantics></math>.</p>
Full article ">Figure 4
<p>Control variables.</p>
Full article ">Figure 5
<p>State variables.</p>
Full article ">Figure 6
<p>Output variables.</p>
Full article ">Figure 7
<p>Control variables.</p>
Full article ">Figure 8
<p>State variables.</p>
Full article ">Figure 9
<p>Output variables.</p>
Full article ">Figure 10
<p>Estimated probability distribution <math display="inline"><semantics> <mi>ξ</mi> </semantics></math>.</p>
Full article ">Figure 11
<p>Equilibrium values, CO<sub>2</sub> lifetimes.</p>
Full article ">Figure 12
<p>Control variables.</p>
Full article ">Figure 13
<p>Equilibrium values, green energy use.</p>
Full article ">Figure 14
<p>Distribution of the social cost of carbon (<span class="html-italic">x</span>), according to the studies discussed (<span>$</span>).</p>
Full article ">Figure 15
<p>Control variables, simulation with <math display="inline"><semantics> <mi>π</mi> </semantics></math>.</p>
Full article ">Figure 16
<p>State variables, simulation with <math display="inline"><semantics> <mi>π</mi> </semantics></math>.</p>
Full article ">Figure 17
<p>Output variables, simulation with <math display="inline"><semantics> <mi>π</mi> </semantics></math>.</p>
Full article ">Figure 18
<p>Equilibrium values, added values.</p>
Full article ">Figure A1
<p>Population growth rate for OECD and (high-income) non-OECD members.</p>
Full article ">Figure A2
<p>Equilibrium values on grid, OECD countries: fitted plane.</p>
Full article ">Figure A3
<p>Equilibrium values on grid, non-OECD countries: fitted plane.</p>
Full article ">Figure A4
<p>Equilibrium values with stochastic <math display="inline"><semantics> <msub> <mi>π</mi> <mi>i</mi> </msub> </semantics></math>.</p>
Full article ">Figure A5
<p>Equilibrium values OECD countries with different <math display="inline"><semantics> <msub> <mi>π</mi> <mi>i</mi> </msub> </semantics></math>.</p>
Full article ">Figure A6
<p>Equilibrium values non-OECD countries with different <math display="inline"><semantics> <msub> <mi>π</mi> <mi>i</mi> </msub> </semantics></math>.</p>
Full article ">Figure A7
<p><span class="html-italic">e</span>, <span class="html-italic">f</span> and <span class="html-italic">g</span> for OECD countries.</p>
Full article ">Figure A8
<p><span class="html-italic">e</span>, <span class="html-italic">f</span> and <span class="html-italic">g</span> for non-OECD countries.</p>
Full article ">Figure A9
<p>Simulation with <math display="inline"><semantics> <mover accent="true"> <mi>k</mi> <mo>˙</mo> </mover> </semantics></math>, <math display="inline"><semantics> <mi mathvariant="sans-serif">Λ</mi> </semantics></math> stochastic: equilibrium values.</p>
Full article ">Figure A10
<p>Simulation with <math display="inline"><semantics> <mover accent="true"> <mi>k</mi> <mo>˙</mo> </mover> </semantics></math>, <math display="inline"><semantics> <mi mathvariant="sans-serif">Λ</mi> </semantics></math> stochastic: objective function values.</p>
Full article ">Figure A11
<p>Equilibrium values on grid, OECD countries: fitted planes.</p>
Full article ">Figure A12
<p>Equilibrium values on grid, non-OECD countries: fitted planes.</p>
Full article ">
10 pages, 2123 KiB  
Article
Preparation of Solid Fuel Hydrochar over Hydrothermal Carbonization of Red Jujube Branch
by Zhiyu Li, Weiming Yi, Zhihe Li, Chunyan Tian, Peng Fu, Yuchun Zhang, Ling Zhou and Jie Teng
Energies 2020, 13(2), 480; https://doi.org/10.3390/en13020480 - 18 Jan 2020
Cited by 18 | Viewed by 3296
Abstract
Biomass energy is becoming increasingly important, owing to the decreasing supply of fossil fuels and growing environmental problems. Hydrothermal carbonization (HTC) is a promising technology for producing solid biofuels from agricultural and forestry residues because of its lower fossil-fuel consumption. In this study, [...] Read more.
Biomass energy is becoming increasingly important, owing to the decreasing supply of fossil fuels and growing environmental problems. Hydrothermal carbonization (HTC) is a promising technology for producing solid biofuels from agricultural and forestry residues because of its lower fossil-fuel consumption. In this study, HTC was used to upgrade red jujube branch (RJB) to prepare hydrochar at six temperatures (220, 240, 260, 280, 300, and 320 °C) for 120 min, and at 300 °C for 30, 60, 90, and 120 min. The results showed that the energy recovery efficiency (ERE) reached maximum values of 80.42% and 79.86% at a residence time of 90 min and a reaction temperature of 220 °C, respectively. X-ray diffraction results and Fourier transform infrared spectroscopy measurements show that the microcrystal features of RJB were destroyed, whereas the hydrochar contained an amorphous structure and mainly lignin fractions at increased temperatures. Thermogravimetric analysis shows that the hydrochar had better fuel qualities than RJB, making hydrochar easier to burn. Full article
Show Figures

Figure 1

Figure 1
<p>Effect of residence time (300 °C) and reaction temperature (120 min) on the yield of solid, liquid, and gas. (<b>a</b>) Solid yield; (<b>b</b>) liquid yield; (<b>c</b>) gas yield; (<b>d</b>) solid yield; (<b>e</b>) liquid yield; (<b>f</b>) gas yield.</p>
Full article ">Figure 2
<p>Energy recovery efficiency and higher heating value of hydrochar after hydrothermal carbonization (HTC) treatment (residence time (300 °C) and reaction temperature (120 min). (<b>a</b>) Residence time; (<b>b</b>) reaction temperature; (<b>c</b>) residence time; (<b>d</b>) reaction temperature.</p>
Full article ">Figure 3
<p>XRD spectra for RJB and hydrochar produced by HTC. (<b>a</b>) Residence time; (<b>b</b>) reaction temperature.</p>
Full article ">Figure 4
<p>FT-IR spectra of RJB and its derived hydrochars. (<b>a</b>) Residence time; (<b>b</b>) reaction temperature.</p>
Full article ">Figure 5
<p>TGA of red jujube branch hydrochar. (<b>a</b>) TG of residence time; (<b>b</b>) derivative thermogravimetry (DTG) of residence time; (<b>c</b>) TG of reaction temperature; (<b>d</b>) DTG of reaction temperature. Note: The raw material is abbreviated as RM.</p>
Full article ">
12 pages, 2186 KiB  
Article
Production, Characterization, and Evaluation of Pellets from Rice Harvest Residues
by Cristina Moliner, Alberto Lagazzo, Barbara Bosio, Rodolfo Botter and Elisabetta Arato
Energies 2020, 13(2), 479; https://doi.org/10.3390/en13020479 - 18 Jan 2020
Cited by 16 | Viewed by 4129
Abstract
Pellets from residues from rice harvest (i.e., straw and husk) were produced and their main properties were evaluated. Firstly, rice straw pellets were produced at lab scale at varying operational conditions (i.e., load compression and wt % of feeding moisture content) to evaluate [...] Read more.
Pellets from residues from rice harvest (i.e., straw and husk) were produced and their main properties were evaluated. Firstly, rice straw pellets were produced at lab scale at varying operational conditions (i.e., load compression and wt % of feeding moisture content) to evaluate their suitability for palletization. Successively, rice straw and husk pellets were commercially produced. All the samples were characterized in terms of their main physical, chemical, and physico-chemical properties. In addition, axial/diametral compression and durability tests were performed to assess their mechanical performance. All the analyzed properties were compared with the established quality standards for non-woody pellets. In general, rice straw pellets presented suitable properties for their use as pelletized fuels. Rice husk pellets fell out of the standards in recommended size or durability and thus preliminary treatments might be required prior their use as fuels. Full article
Show Figures

Figure 1

Figure 1
<p>Particle size distribution of grounded rice straw.</p>
Full article ">Figure 2
<p>Pellets produced using the single press.</p>
Full article ">Figure 3
<p>N-Pico pelletizer (<b>a</b>) and pelletized rice straw (RS) (<b>b</b>).</p>
Full article ">Figure 4
<p>Axial compressive test.</p>
Full article ">Figure 5
<p>Fraction 0.3–0.85 mm (<b>a</b>) and 0.15–0.3 mm (<b>b</b>) of the straw powders at the magnification of 50× at optical microscope.</p>
Full article ">Figure 6
<p>RS (<b>a</b>) and RH (<b>b</b>) pellets.</p>
Full article ">Figure 7
<p>Axial compressive stress-strain curves for RS pellets (<b>a</b>) and RH pellets (<b>b</b>).</p>
Full article ">
16 pages, 5246 KiB  
Article
An Improved State of Charge and State of Power Estimation Method Based on Genetic Particle Filter for Lithium-ion Batteries
by Xingtao Liu, Chaoyi Zheng, Ji Wu, Jinhao Meng, Daniel-Ioan Stroe and Jiajia Chen
Energies 2020, 13(2), 478; https://doi.org/10.3390/en13020478 - 18 Jan 2020
Cited by 34 | Viewed by 3706
Abstract
In this paper, an improved method for estimating the state of charge (SOC) of lithium-ion batteries is proposed, which is developed from the particle filter (PF). An improved genetic particle filter (GPF), owing to the advantages of the PF and genetic algorithm, is [...] Read more.
In this paper, an improved method for estimating the state of charge (SOC) of lithium-ion batteries is proposed, which is developed from the particle filter (PF). An improved genetic particle filter (GPF), owing to the advantages of the PF and genetic algorithm, is proposed to overcome the disadvantage of the traditional particle filter: lacking the diversity of particles. Firstly, the relationship between SOC and open-circuit voltage (OCV) is identified on the low-current OCV test. Secondly, a first-order resistor and capacitance (RC) model is established, then, the least-squares algorithm is used to identify the model parameters via the incremental current test. Thirdly, GPF and the improved GPF (IGPF) are proposed to solve the problems of the PF. The method based on the IGPF is proposed to estimate the state of power (SOP). Finally, IGPF, GPF, and PF are employed to estimate the SOC on the federal urban driving schedule (FUDS). The results show that compared with traditional PF, the errors of the IGPF are 20% lower, and compared with GPF, the maximum error of the IGPF has declined 1.6% SOC. The SOC that is estimated by the IGPF is applied to estimate the SOP for battery, considering the restrictions from the peak SOC, the voltage, and the instruction manual. The result shows that the method based on the IGPF can successfully estimate SOP. Full article
(This article belongs to the Special Issue Testing and Management of Lithium-Ion Batteries)
Show Figures

Figure 1

Figure 1
<p>Resistor and capacitance equivalent circuit model.</p>
Full article ">Figure 2
<p>Current and voltage of low-current open-circuit voltage test.</p>
Full article ">Figure 3
<p>Curves at low-current OCV.</p>
Full article ">Figure 4
<p>Errors between the real OCV and the calculated OCV.</p>
Full article ">Figure 5
<p>Current and voltage of the incremental current test.</p>
Full article ">Figure 6
<p>Process of the improved genetic particle filter method.</p>
Full article ">Figure 7
<p>Current and voltage of federal urban driving schedule.</p>
Full article ">Figure 8
<p>SOC estimation for the FUDS driving cycle.</p>
Full article ">Figure 9
<p>Estimation errors for the FUDS driving cycle.</p>
Full article ">Figure 10
<p>SOP and charging SOP on the FUDS.</p>
Full article ">Figure 11
<p>(<b>a</b>) SOP and (<b>b</b>) charging SOP with different limitations on the FUDS.</p>
Full article ">
24 pages, 4013 KiB  
Article
An Energy Optimal Dispatching Model of an Integrated Energy System Based on Uncertain Bilevel Programming
by Xueying Song, Hongyu Lin, Gejirifu De, Hanfang Li, Xiaoxu Fu and Zhongfu Tan
Energies 2020, 13(2), 477; https://doi.org/10.3390/en13020477 - 18 Jan 2020
Cited by 10 | Viewed by 2366
Abstract
An integrated energy system (IES) involving a large number of decision-makers causes problems of bad coordination between energy sub-networks and the IES and it is not able to fully consider the multi-energy complementarity among multiple decision-makers. In this context, firstly, this paper constructs [...] Read more.
An integrated energy system (IES) involving a large number of decision-makers causes problems of bad coordination between energy sub-networks and the IES and it is not able to fully consider the multi-energy complementarity among multiple decision-makers. In this context, firstly, this paper constructs an energy optimal dispatching model of an IES based on uncertain bilevel programming. The upper model takes the transformation matrix of energy hubs as the upper decision-maker, taking the minimum operation cost of the IES in the form of confidence as the objective function; the lower model takes each optimal operation plan of the electric power sub-network, the thermal energy sub-network, and the gas energy sub-network as the lower decision-makers, aiming at the operation economy of each sub-network and considering their operation as necessary constraints. Secondly, a firefly algorithm with chaotic search and an improved light intensity coefficient is designed to improve the proposed model. An empirical analysis was conducted on a pilot area of an integrated energy system in Hebei Province. The results show the following: (1) The typical daily operating cost of the integrated energy system in winter is lower than that in summer; (2) under the same load level, the typical winter and summer running costs of the integrated energy system are lower than that of the traditional microgrid; (3) compared with the particle swarm optimization algorithm, the improved firefly algorithm proposed in the paper has obvious advantages both in terms of running cost and solution time; and (4) when the confidence of the objective function and the constraints increases, the operating cost of various schemes also increase. The above results validate the effectiveness of the energy optimal dispatching model of the IES and the economy of the system operation under the multiple decision-maker hierarchy. Full article
Show Figures

Figure 1

Figure 1
<p>Basic component module and energy flow relation between energy subgrids in micro-energy Internet.</p>
Full article ">Figure 2
<p>The information flow relation diagram of the bilevel programming model for an IES.</p>
Full article ">Figure 3
<p>A model solving flow chart based on an improved imperialist competitive algorithm.</p>
Full article ">Figure 4
<p>Wind power output, photovoltaic power output, power load, and heat load curves within a typical day in summer (Photovoltaic: PV).</p>
Full article ">Figure 5
<p>Wind power output, photovoltaic power output, power load, and cooling load curves within a typical day in winter (Photovoltaic: PV).</p>
Full article ">Figure 6
<p>An IES optimization operation scheme for a typical winter day (Micro-turbine: MT, FC: Fuel Cell, SB: Storage Battery).</p>
Full article ">Figure 7
<p>Operating cost curve of the IES for each period of a typical winter day.</p>
Full article ">Figure 8
<p>Optimal operation scheme of the IES for a typical summer day (Micro-turbine: MT, FC: Fuel Cell, SB: Storage Battery).</p>
Full article ">Figure 9
<p>Operating cost curve of the IES for each period of a typical summer day.</p>
Full article ">Figure 10
<p>The relationship between the planned operating cost and the confidence degree of the IES for a typical winter day.</p>
Full article ">Figure 11
<p>The relationship between the planned operating cost and the confidence degree of the IES for a typical summer day.</p>
Full article ">
Previous Issue
Back to TopTop