• Energy Research

In this paper different gas-steam combined-cycles fueled by syngas produced in a local downdraft gasifier, are analyzed. At first, the downdraft gasifier model is briefly described, where waste biomass is transformed into syngas, which can be used more efficiently than the original solid biomass to generate useful power, and can be transported much more easily. The gasifier model is able to estimate, with good approximation, the composition of the produced syngas, taking separately into account the biomass drying and the pyrolysis, oxidation and reduction processes. The gasifier operates at ambient pressure using air as gasification agent and biomass as input. Among others, pomace has been considered, since, in Italy (where the plant is supposed to be located) there are many regions, like Apulia, where this biomass is largely available. Three different plant configurations have been proposed and compared in terms of overall performance. The first two, named REXC (Regenerative cycle with EXternal Combustor) and CR (Conventional Regenerative cycle), burn the syngas in an internal combustor, whereas the third one, named SyEXC (Syngas External Combustion), considers an externally fired configuration for the syngas combustion. In the REXC cycle, a secondary external combustion system, fed by cellulosic biomass, is connected to a heat exchanger in order to increase the air temperature, as in a regenerative cycle. The combustion products pass through a primary heat exchanger placed in the external combustion system, heating the compressed air, which flows into the primary internal combustion chamber, where a defined quantity of syngas reacts with the compressed air. The turbine exhaust gas (TEG), before going back into the external combustion, partially transfers its enthalpy content to a Heat Recovery Steam Generator (HRSG) of a bottoming Rankine cycle. The second plant configuration (CR cycle) is characterized by a main internal combustor fueled by the produced syngas and implements a conventional regenerative process. The TEG is reheated by an externally fired combustion before flowing in the HRSG of a bottoming Rankine cycle. The last plant configuration (SyEXC) differs from the CR one only for the main externally fired gas turbine fueled by syngas, thus avoiding any costly cleanup operation and cooling.

A latent heat storage system for concentrated solar plants (CSP) is numerically examined by means of CFD simulations. This study aims at identifying the convective flows produced within the melted phase by temperature gradients and gravity. Simulations were carried out on experimental devices for applications to high temperature concentrated solar power plants. A shell-and-tube geometry composed by a vertical cylindrical tank, filled by a Phase Change Material (PCM) and an inner steel tube, in which the heat transfer fluid (HTF) flows, from the top to the bottom, is considered. The conjugate heat transfer process is examined by solving the unsteady Navier–Stokes equations for HTF and PCM and conduction for the tube. In order to take into account the buoyancy effects in the PCM tank the Boussinesq approximation is adopted. The results show that the enhanced heat flux, due to natural convective flow, reduce of about 30% the time needed to charge the heat storage. A detailed description of the convective motion in the melted phase and the heat flux distribution between the HTF and PCM are reported. The effect of the mushy zone constant is also investigated.

Sustainable development can no longer neglect the growth of those technologies that look at the recovery of any energy waste in industrial processes. For example, in almost every industrial plant it happens that pressure energy is wasted in throttling devices for pressure and flow control needs. Clearly, the recovery of this wasted energy can be considered as an opportunity to reach not only a higher plant energy efficiency, but also the reduction of the plant Operating Expenditures (OpEx). In recent years, it is getting common to replace throttling valves with turbine-based systems (tuboexpander) thus getting both the pressure control and the energy recovery, for instance, producing electricity. However, the wide range of possible operating conditions, technical requirements and design constrains determine highly customized constructions of these turboexpanders. Furthermore, manufacturers are interested in tools enabling them to rapidly get the design of their products. For these reasons, in this work we propose an optimization design procedure, which is able to rapidly come to the design of the turboexpander taking into account all the fluid dynamic and technical requirements, considering the already obtained achievements of the scientific community in terms of theory, experiments and numeric. In order to validate the proposed methodology, the case of a single stage axial impulse turbine is considered. However, the methodology extension to other turbomachines is straightforward. Specifically, the design requirements were expressed in terms of maximum allowable expansion ratio and flow coefficient, while achieving at least a minimum assigned value of the turbine loading factor. Actually, it is an iterative procedure, carried out up to convergence, made of the following steps: (i) the different loss coefficients in the turbine are set-up in order to estimate its main geometric parameters by means of a one dimensional (1D) study; (ii) the 2D blade profiles are designed by means of an optimization algorithm based on a “viscous/inviscid interaction” technique; (iii) 3D Computational Fluid Dynamic (CFD) simulations are then carried out and the loss coefficients are computed and updated. Regarding the CFD simulations, a preliminary model assessment has been performed against a reference case, chosen in the literature. The above-mentioned procedure is implemented in such a way to speed up the convergence, coupling analytical integral models of the 1D/2D approach with accurate local solutions of the finite-volume 3D approach. The method is shown to be able to achieve consistent results, allowing the determination of a turbine design respectful of the requirements more than doubling the minimum required loading factor.

Water Distribution Networks (WDNs) represent a noteworthy field for possible implementation of Small Hydropower (SHP), by replacing Pressure Reduction Valves (PRV) with turbomachines, in particular Pump as Turbines (PaTs), to control and regulate the pressure, while harvesting energy otherwise wasted. Different models were developed to predict the performance and select the positioning of the PaTs for the maximum energy recovery but most of them neglect practical aspect such as: power grid limitations and optimal harvesting strategy. In this framework, we intend to propose a new method to select a PaT, defining its optimal working point, by introducing an energy exploitation coefficient. The proposed methodology is based on the experimental results of a real PaT tested in the high capacity hydraulic laboratory at Polytechnic University of Bari. Firstly, the selected commercial centrifugal pump was tested in both pump and turbine modes. Then, three different approaches, for the Best Efficiency Point (BEP) selection, are described and compared in terms of energy exploitation and capacity factor for a WDN. The first consists of selecting the BEP at the average flow rate, the second one considers the probability distribution of the flow rate and the corresponding available hydraulic energy, whereas the latter is based on the highest energy harvesting. By applying energy production, economic and environmental analyses, the new proposed methodology, based on the third approach, shows a remarkable advantage in terms of exploited energy. Indeed a remarkable 60% energy recovery is achieved with 334 ton CO2/year avoided. Furthermore, the impact of the electrical motor on the maximum power generation (cut-off) is considered. Eventually, useful insights for the future PaT selection and installation are discussed.

A high-fidelity real-time simulation code based on a lumped, nonlinear representation of gas turbine components is presented. The code is a general-purpose simulation software environment useful for setting up and testing control equipments. The mathematical model and the numerical procedure are specially developed in order to efficiently solve the set of algebraic and ordinary differential equations that describe the dynamic behavior of gas turbine engines. For high-fidelity purposes, the mathematical model takes into account the actual composition of the working gases and the variation of the specific heats with the temperature, including a stage-by-stage model of the air-cooled expansion. The paper presents the model and the adopted solver procedure. The code, developed in Matlab-Simulink using an object-oriented approach, is flexible and can be easily adapted to any kind of plant configuration. Simulation tests of the transients after load rejection have been carried out for a single-shaft heavy-duty gas turbine and a double-shaft aero-derivative industrial engine. Time plots of the main variables that describe the gas turbine dynamic behavior are shown and the results regarding the computational time per time step are discussed.

An object-oriented program for the dynamic simulation of gas turbines has been developed using Matlab-Simulink® utilities. The advantages given by the object-oriented program are the flexibility and the user-friendly interface. The components of the gas turbine are modeled as blocks that can be assembled like an engineering drawing in which the connections between two elements represent either a mechanical power transfer or a fluid transfer. In a fluid-type connection the data regarding the fluid composition, the thermodynamic properties and the fluid velocity are simultaneously transferred. A library of gas turbine components has been developed: each component appears as a block that can be placed into the drawing by means of click-and-drag operations. In order to obtain an accurate description of the physical processes characterizing the gas turbine, the components are described by a set of non-linear algebraic equations and ordinary differential equations. Each component block is created using the graphical programming tools of Simulink, following a block diagram approach. Specific blocks have been created for the evaluation of the thermodinamic properties of the working fluids. The paper provides the description of the mathematical model adopted for the simulation and the some examples of the Matlab-Simulink formulation adopted for implementing the model. The accuracy of the program has been verified considering test cases for which experimental data are available in the public literature. The computational times, tested for the transients of a single-shaft and a double-shaft gas turbine, resulted to be very short, even making use of a personal computer.

Abstract A theoretical model to predict the complete melting time of a Latent Heat Thermal Energy Storage (LHTES) considering a shell-and-tube configuration is presented and applied to several geometries. The shape of the device has been modified according to the internal radius ( r i ), external radius ( r e ) and total height (L) retaining constant the volume of the storage and the heat exchange area. The model has been validated by means of multiphase numerical simulations of the charging phase. The numerical simulations have been performed considering four configurations, with a radii ratio, r e / r i , equal to 1.5, 2, 4.375 and 6. The comparison between the model predictions and the numerical simulations confirms the reliability of the theoretical model in terms of melting time within the range investigated. The study reveals that, even considering the same storage volume, heat exchange area and wall temperature, for low values of radii ratio ( r e / r i ), the shape of the device is able to reduce the charging time of the LHTES up to 50% for a radii ratio r e / r i = 1.5 with respect to r e / r i = 4.375 . Increasing further the radii ratio from r e / r i = 4.375 to 6, the melting time decreases. The unsteady numerical simulations support the prediction of the theoretical model. Thus, in the here studied geometrical configurations the proposed approach represents a simplified and accurate design tool to predict the charging time of a LHTES shell-and-tube device.

In the present paper the new wind tunnel located in the Fluid-dynamic Laboratory of the Dipartimento di Ingegneria Meccanica e Gestionale (DIMeG) of the Bari Polytechnic will be shortly described and the first experimental measurements on a vertical axis wind turbine (VAWT) will be shown. The DIMeG wind tunnel has been designed by the research group on wind energy of the Department. It is a subsonic, closed loop, wind tunnel with a transparent test part where small scale models can be analyzed. A four bladed axial fan is driven by an asynchronous three phase electric motor, which is connected to an inverter in order to change the wind speed. Angular blades have been inserted at the two curves between the fan and the test section in order to increase the uniformity of the velocity profile after the two curves. An optimization fluid-dynamic study has been carried out in order to find the best blade profile. A honeycomb has been also inserted upstream the test section in order to destroy the still existing small vorticity generated by the fan and the curves. A three-axis traversing, called Cartesian robot, has been designed and built above the test section, in order to move the hot wire probe, for wind speed measurements, by means of four step by step electric motors controlled by a personal computer. A data acquisition system has been set up. All the principal commands and controls can be performed by a dedicated personal computer, which has been programmed using LabVIEW® G-programming language. The first experimental activity has been performed on a VAWT model, of the Giromill type with parallel blades. The turbine has been connected to an AC brushless servo, able to control the braking torque. Experimental results of the flow field in two horizontal planes have been set up using a two component hot wire probe (Dantec 55R51) calibrated with the manual system Dantec 54H10. The measurement grid adopted is formed by 20 nodes in the Y direction (main flow direction) and 10 nodes in the X direction.