To achieve the carbon free electricity generation target for 2050, the penetration of renewable energy sources should further increase. To address the impacts of their unpredictable and intermittent characteristics on the future electricity grid, Pumped Hydro Energy Storage (PHES) plants should enhance their regulation capability by extending their continuous operating range far beyond the optimal normal working range. However, for the time being, the regulation capability of the new generation of PHES, equipped with reversible pump-turbines due to their cost-effectiveness, is limited at part load by instability problems. The aim of this paper is to analyse, during a pumping power reduction scenario, the onset and development of unsteady phenomena leading to unstable behaviour. A 3D transient numerical simulation was carried out on the first stage of a variable-speed two-stage pump-turbine from full load to the unstable operating zone by progressively reducing the speed from 100% to 88... [more]

Utilization of hydraulic-fracturing technology is dramatically increasing in exploitation of natural gas extraction. However the prediction of the configuration of propagated hydraulic fracture is extremely challenging. This paper presents a numerical method of obtaining the configuration of the propagated hydraulic fracture into discrete natural fracture network system. The method is developed on the basis of weighted fracture which is derived in combination of Dijkstra’s algorithm energy theory and vector method. Numerical results along with experimental data demonstrated that proposed method is capable of predicting the propagated hydraulic fracture configuration reasonably with high computation efficiency. Sensitivity analysis reveals a number of interesting observation results: the shortest path weight value decreases with increasing of fracture density and length, and increases with increasing of the angle between fractures to the maximum principal stress direction. Our method is... [more]

Two modeling approaches, the scaling-law and CFD (Computational Fluid Dynamics) approaches, are presented in this paper. To save on experimental cost of the pilot plant, the scaling-law approach as a low-computational-cost method was adopted and a small scale column operating under ambient temperature and pressure was built. A series of laboratory tests and computer simulations were carried out to evaluate the hydrodynamic characteristics of a pilot fluidized-bed biomass gasifier. In the small scale column solids were fluidized. The pressure and other hydrodynamic properties were monitored for the validation of the scaling-law application. In addition to the scaling-law modeling method, the CFD approach was presented to simulate the gas-particle system in the small column. 2D CFD models were developed to simulate the hydrodynamic regime. The simulation results were validated with the experimental data from the small column. It was proved that the CFD model was able to accurately predic... [more]

The authors wish to make the following corrections to the published paper [1].[...]

Parameterization of battery dynamics based on terminal operating data is a main concern in engineering applications of batteries. The key technology is designing an adequate test procedure and a data processing procedure to excite different inner dynamics and then estimate the parameters of a corresponding equivalent circuit model (ECM). This paper proposes a parameter identification method that utilizes the terminal voltage curves (TVC) under step-change current conditions and constant current conditions. With this method, I-V characteristics of battery’s Ohmic resistance, mass diffusion process, thermal process and SOC varying process are decoupled and parametric functions of an ECM are obtained. Experimental results show that the method is easy to be implemented and modeling accuracy is sufficient for applications.

The thermo-structural response of solid rocket motor nozzles is widely investigated in the design of modern rockets, and many factors related to the material properties have been considered. However, little work has been done to evaluate the effects of structure gaps on the generation of flame leaks. In this paper, a numerical simulation was performed by the finite element method to study the thermo-structural response of a typical nozzle with consideration of the structure gap. Initial boundary conditions for thermo-structural simulation were defined by a quasi-1D model, and then coupled simulations of different gap size matching modes were conducted. It was found that frictional interface treatment could efficiently reduce the stress level. Based on the defined flame leak criteria, gap size optimization was carried out, and the best gap matching mode was determined for designing the nozzle. Testing experiment indicated that the simulation results from the proposed method agreed well... [more]

In selecting a mathematical model for simulating physical behaviours, it is important to reach an acceptable compromise between analytical complexity and achievable precision. With the aim of helping researchers and designers working in the area of photovoltaic systems to make a choice among the numerous diode-based models, a criterion for rating both the usability and accuracy of one-diode models is proposed in this paper. A three-level rating scale, which considers the ease of finding the data used by the analytical procedure, the simplicity of the mathematical tools needed to perform calculations and the accuracy achieved in calculating the current and power, is used. The proposed criterion is tested on some one-diode equivalent circuits whose analytical procedures, hypotheses and equations are minutely reviewed along with the operative steps to calculate the model parameters. To assess the achievable accuracy, the current-voltage (I-V) curves at constant solar irradiance and/or cel... [more]

Power devices intended for high-voltage systems must be tested according to international standards, which includes the short-time withstand current test and peak withstand current test. However, these tests require very special facilities which consume huge amounts of electrical power. Therefore, mathematical tools to simulate such tests are highly appealing since they allow reproducing the electromagnetic and thermal behavior of the test object in a fast and economical manner. In this paper, a three-dimensional finite element method (3D-FEM) approach to simulate the transient thermal behavior of substation connectors is presented and validated against experimental data. To this end, a multiphysics 3D-FEM method is proposed, which considers both the connector and the reference power conductors. The transient and steady-state temperature profiles of both the conductors and connector provided by the 3D-FEM method prove its suitability and accuracy as compared to experimental data provid... [more]

This study presents the heat transfer and fluid flow characteristics in a rib-roughened SAH (solar air heater) channel. The artificial roughness of the rectangular channel was in the form of a thin circular wire in discrete multi V-pattern rib geometries. The effect of this geometry on heat transfer, fluid flow, and performance augmentation was investigated using the CFD (computational fluid dynamics). The roughness parameters were a relative discrete distance of 0.69, a relative rib height of 0.043, a relative rib pitch of 10, a relative rib width of 6.0, and a flow-attack-angle of 60°. The discrete width ratios and Reynolds numbers ranged from 0.5 to 2.0 and from 2000 to 20,000, respectively. The CFD results using the renormalization k-epsilon model were in good agreement with the empirical relationship. This model was used to investigate the heat transfer and fluid flow characteristics in the multi V-pattern rib roughened SAH channel. The thermo-hydraulic performance was found to be... [more]

In this study, variable injection-rate technology was numerically investigated in a pre-existing discrete fracture network (DFN) formation, the Tarim Basin in China. A flow-stress-damage (FSD) coupling model has been used in an initial attempt towards how reservoir response to variable injection-rates at different hydraulic fracturing stages. The established numerical model simultaneously considered the macroscopic and microscopic heterogeneity characteristics. Eight numerical cases were studied. Four cases were used to study the variable injection-rate technology, and the other four cases were applied for a constant injection-rate in order to compare with the variable injection-rate technology. The simulation results show that the variable injection-rate technology is a potentially good method to a form complex fracturing networks. The hydraulic fracturing effectiveness when increasing the injection-rate at each stage is the best, also, the total injected fluid is at a minimum. At the... [more]

A physics-based explicit mathematical model for the external voltage-dependent forward dark current in bulk heterojunction (BHJ) organic solar cells is developed by considering Shockley-Read-Hall (SRH) recombination and solving the continuity equations for both electrons and holes. An analytical model for the external voltage-dependent photocurrent in BHJ organic solar cells is also proposed by incorporating exponential photon absorption, dissociation efficiency of bound electron-hole pairs (EHPs), carrier trapping, and carrier drift and diffusion in the photon absorption layer. Modified Braun’s model is used to compute the electric field-dependent dissociation efficiency of the bound EHPs. The overall net current is calculated considering the actual solar spectrum. The mathematical models are verified by comparing the model calculations with various published experimental results. We analyze the effects of the contact properties, blend compositions, charge carrier transport properties... [more]

A novel modular arc-linear flux-switching permanent-magnet motor (MAL-FSPM) used for scanning system instead of reduction gearboxes and kinematic mechanisms is proposed and researched in this paper by the finite element method (FEM). The MAL-FSPM combines characteristics of flux-switching permanent-magnet motor and linear motor and can realize the direct driving and limited angular movement. Structure and operation principle of the MAL-FSPM are analyzed. Cogging torque model of the MAL-FSPM is established. The characteristics of cogging torque and torque ripple are investigated for: (1) distance (dend) between left end of rotor and left end of stator is more than two rotor tooth pitch (τp); and (2) dend is less than two rotor tooth pitch. Cogging torque is an important component of torque ripple and the period ratio of the cogging torque to the back electromotive force (EMF) equals one for the MAL-FSPM before optimization. In order to reduce the torque ripple as much as possible and af... [more]

Latent heat thermal energy storage (TES) plays an important role in the advocation of TES in contrast to sensible energy storage because of the large storage energy densities per unit mass/volume possible at a nearly constant thermal energy. In the current study, a heat exchange device with a zigzag configuration containing multiple phase-change-materials (m-PCMs) was considered, and an experimental system was built to validate the model for a single PCM. A two-dimensional numerical model was developed using the ANSYS Fluent 14.0 software program. The energy fractions method was put forward to calculate the average Ste number and the influence of Re and Ste numbers on the discharge process were studied. The influence of phase change temperature among m-PCMs on the solidification process has also been studied. A new boundary condition was defined to determine the combined effect of the Re and Ste numbers on the discharging process. The modelling results show that for a given input power... [more]

To improve vehicle fuel economy whilst enhancing road handling and ride comfort, power generating suspension systems have recently attracted increased attention in automotive engineering. This paper presents our study of a regenerative hydraulic shock absorber system which converts the oscillatory motion of a vehicle suspension into unidirectional rotary motion of a generator. Firstly a model which takes into account the influences of the dynamics of hydraulic flow, rotational motion and power regeneration is developed. Thereafter the model parameters of fluid bulk modulus, motor efficiencies, viscous friction torque, and voltage and torque constant coefficients are determined based on modelling and experimental studies of a prototype system. The model is then validated under different input excitations and load resistances, obtaining results which show good agreement between prediction and measurement. In particular, the system using piston-rod dimensions of 50⁻30 mm achieves recovera... [more]

In a future sustainable energy vision, in which diversified conversion of renewable energies is essential, vertical axis wind turbines (VAWTs) exhibit some potential as a reliable means of wind energy extraction alongside conventional horizontal axis wind turbines (HAWTs). Nevertheless, there is currently a relative shortage of scientific, academic and technical investigations of VAWTs as compared to HAWTs. Having this in mind, in this work, we aim to, for the first time, study the wake of a single VAWT placed in the atmospheric boundary layer using large-eddy simulation (LES). To do this, we use a previously-validated LES framework in which an actuator line model (ALM) is incorporated. First, for a typical three- and straight-bladed 1-MW VAWT design, the variation of the power coefficient with both the chord length of the blades and the tip-speed ratio is analyzed by performing 117 simulations using LES-ALM. The optimum combination of solidity (defined as N c / R , where N is t... [more]

At present, the common practice in the power system of China is to represent the small and medium hydropower generator group as a negative load. This paper presents a method to build a dynamic equivalent model of the hydropower generator group using a 3rd order generator model and a static characteristic load model. Based on phasor measurements in the tie line which connects to the modeled hydropower generator group, the dynamic multi-swarm particle swarm optimizer (DMS-PSO) algorithm is used to obtain parameters of the equivalent model. The proposed method is verified in the small and medium hydropower generator group of Sichuan power grid with both simulation and actual data. The results show that the dynamic responses and the transient stability are consistent before and after the equivalence. The proposed method can be used for modeling a group of small and medium hydropower generators whose structures and parameters are unknown.

The paper proposes two mathematical models of a photo-voltaic (PV) cell—the complete model and the simplified model—which can be used also for modeling a PV module or a PV string under any environmental condition. Both of them are based on the well-known five-parameters model, while the approach allows to write a new descriptive equation, whose terms are functions of the information always available in the modern datasheet of a PV module’s manufacturer. This implies that no pre-processing of the datasheet parameters is needed to use the proposed model, whichever the solar irradiance and the cell/module temperature are. Moreover, these models are interpreted from a circuital point of view, providing the electrical circuits constituted only by basic electrical components. Particularly, in order to take into account the variability of the environment parameters, several variable resistors and voltage-controlled sources are used. The proposed models are tested with the datasheet parameters... [more]

Multiple-factor analysis and optimization play a critical role in the the ability to maximizethe stimulated reservoir volume (SRV) and the success of economic shale gas production. In this paper, taking the typical continental naturally fractured silty laminae shale in China as anexample, response surface methodology (RSM) was employed to optimize multiple hydraulic fracturing parameters to maximize the stimulated area in combination with numerical modeling based on the coupled flow-stress-damage (FSD) approach. This paper demonstrates hydraulic fracturing effectiveness by defining two indicesnamelythe stimulated reservoir area (SRA) and stimulated silty laminae area (SLA). Seven uncertain parameters, such as laminae thickness, spacing, dip angle, cohesion, internal friction angle (IFA), in situ stress difference (SD), and an operational parameter-injection rate (IR) with a reasonable range based on silty Laminae Shale, Southeastern Ordos Basin, are used to fit a response of SRA and SL... [more]

Groundwater flow is one of the most important factors for the design of a ground heat exchanger (GHEX) since the thermal environment of the ground around the buried GHEX is significantly affected by heat convection due to the groundwater flow. Several preceding studies have been conducted to develop analytical solutions to the heat transfer model of GHEX with consideration of groundwater flow. One of these solutions is the combined heat transfer model of conduction and convection. However, the developed combined analytical models are inapplicable to all of the configurations of ordinary GHEXs because these solutions assume that the inner part of the borehole is thermally inert or consists of the same material as that of the surrounding ground. In this paper, the applicability of the combined solid cylindrical heat source model, which is the most suitable to energy piles until now, was evaluated by performing a series of numerical analyses. In the numerical analysis, the inner part of t... [more]

If a turbine loses its electrical load, it will rotate freely and increase speed, eventually achieving that rotational speed which produces zero net torque. This is known as a runaway situation. Unlike many other types of turbine, a marine current turbine will typically overshoot the final runaway speed before slowing down and settling at the runaway speed. Since the hydrodynamic forces acting on the turbine are dependent on rotational speed and acceleration, turbine behaviour during runaway becomes important for load analyses during turbine design. In this article, we consider analytical and numerical models of marine current turbine runaway behaviour in one dimension. The analytical model is found not to capture the overshoot phenomenon, while still providing useful estimates of acceleration at the onset of runaway. The numerical model incorporates turbine wake build-up and predicts a rotational speed overshoot. The predictions of the models are compared against measurements of runaw... [more]

Cavitation is a common flow phenomena in most hydraulic turbines and has the potential to cause vibration, blade surface damage and performance loss. Despite the fact that crossflow turbines have been used in small-scale hydropower systems for a long time, cavitation has not been studied in these turbines. In this paper, we present the findings of a computational study on cavitation inception in crossflow turbines. Cavitation inception was assessed using three-dimensional (3D) Reynolds-Averaged Navier⁻Stokes (RANS) computations. A homogeneous, free-surface two-phase flow model was used. Pressure distributions on the blades were examined for different flow rates, heads and impeller speeds to assess cavitation inception. The results showed that cavitation occurs in the second stage of the turbine and was observed on the suction side near the inner edge of the blades. For the particular turbine studied, cavitation always occurred at shaft speeds greater than that, giving the maximum effic... [more]

The evaluation of the thermodynamic performance of the mutual transformation of different kinds of exergy linked to the intensive thermodynamic parameters of the flow inside the ejector of a refrigeration system is undertaken. Two thermodynamic metrics, exergy produced and exergy consumed, are introduced to assess these transformations. Their calculation is based on the evaluation of the transiting exergy within different ejector sections taking into account the temperature, pressure and velocity variations. The analysis based on these metrics has allowed pinpointing the most important factors affecting the ejector’s performance. A new result, namely the temperature rise in the sub-environmental region of the mixing section is detected as an important factor responsible for the ejector’s thermodynamic irreversibility. The overall exergy efficiency of the ejector as well as the efficiencies of its sections are evaluated based on the proposed thermodynamic metrics.

In order to properly manage lithium-ion batteries of electric vehicles (EVs), it is essential to build the battery model and estimate the state of charge (SOC). In this paper, the fractional order forms of Thevenin and partnership for a new generation of vehicles (PNGV) models are built, of which the model parameters including the fractional orders and the corresponding resistance and capacitance values are simultaneously identified based on genetic algorithm (GA). The relationships between different model parameters and SOC are established and analyzed. The calculation precisions of the fractional order model (FOM) and integral order model (IOM) are validated and compared under hybrid test cycles. Finally, extended Kalman filter (EKF) is employed to estimate the SOC based on different models. The results prove that the FOMs can simulate the output voltage more accurately and the fractional order EKF (FOEKF) can estimate the SOC more precisely under dynamic conditions.

This paper investigates the performance analysis and mathematical modeling of a staggered-tooth transverse-flux permanent magnet linear synchronous machine (STTF-PMLSM), which is characterized by simple structure and low flux leakage. Firstly, the structure advantages and operation principle of the STTF-PMLSM are introduced, and a simplified one phase model is established to investigate the performance of the machine in order to save the computation time. Then, the electromagnetic characteristics, including no-load flux linkage, electromotive force (EMF), inductance, detent force and thrust force, are simulated and analyzed in detail. After that, the theoretical analysis of the detent force, thrust force, and power factor are carried out. And the theoretical analysis results are validated with 3-D finite-element method (FEM). Finally, an improved mathematical model of the machine based on d-q rotating coordinate system is proposed, in which inductance harmonics and coupling between d-... [more]

To capture the hydraulic fractures in heterogeneous and layered rocks, a numerical code that can consider the coupled effects of fluid flow, damage, and stress field in rocks is presented. Based on the characteristics of a typical thin and inter-bedded sedimentary reservoir, China, a series of simulations on the hydraulic fracturing are performed. In the simulations, three points, i.e., (1) confining stresses, representing the effect of in situ stresses, (2) strength of the interfaces, and (3) material properties of the layers on either side of the interface, are crucial in fracturing across interfaces between two adjacent rock layers. Numerical results show that the hydrofracture propagation within a layered sequence of sedimentary rocks is controlled by changing in situ stresses, interface properties, and lithologies. The path of the hydraulic fracture is characterized by numerous deflections, branchings, and terminations. Four types of potential interaction, i.e., penetration, arres... [more]