Abstract
To date, the most industrially developed technology to produce green hydrogen is represented by alkaline water electrolysis (AWE). To improve on design and efficiency of these devices, ...however, multiphysics simulations based on Computational Fluid Dynamics (CFD) are needed, able to account for electrophysical phenomena and multiphase flows. Focusing on internal flow optimization, the requirements for CFD simulations are anyhow extremely challenging, since solving the gas bubbles’ motion implies the solution of a two-phase flow characterized by very low Reynolds numbers and a high fraction of dispersed gas. Despite some interesting studies have been presented in the literature so far, validation of CFD results with detailed experimental measurements is quite rare and, therefore, the reliability of the adopted modelling approaches is not assessed yet. This study presents the results of a multivariate CFD analysis of an electrochemical cell and its validation through a literature test case. Bubbles generation is introduced as a source term, thus overlooking for the moment the electrochemistry to focus on fluid-dynamics. In particular, attention is given to the Eulerian multiphase modelling, investigating the influence of both the inter-phase interaction sub-models’ settings (e.g., lift and drag forces, virtual-mass force) and the general settings of the simulation. The mean velocity field of the PIV-measured bubbles is considered to assess the accuracy of numerical predictions, while the available high-definition flow pictures allow a qualitative assessment of the bubbles size and location. CFD results are shown to be in decent agreement with experimental data and able to reproduce the key flow features such as the spreading of the bubble curtains and the gas shifting towards the inner part of the cell. The effect of the bubbles’ diameter and of source layer thickness is also discussed.
Abstract
Simulation of the complex, unsteady aerodynamics characterizing Darrieus rotors requires computational tools with a fidelity higher than the ubiquitous Blade Element Momentum (BEM) theory. ...Among them, the Actuator Line Method (ALM) stands out in terms of accuracy and computational cost. This approach, however, still fails to resolve the vortex-like structures shed at the blade ends, overestimating turbine performance at the higher rotational speeds. Moving from this background, in this study a comprehensive investigation on the ALM’s capability to simulate tip effects and their impact on rotor performance is carried out. To this end, the ALM tool developed by the authors in the ANSYS
®
FLUENT
®
environment (v. 20.2) and specifically tailored to the simulation of vertical-axis machines was employed. Both a steady finite wing and a fictitious one-blade Darrieus H-rotor, for which high-fidelity blade-resolved CFD data are available as benchmark, were considered as test cases. ALM simulations were first performed without any correction for different cell sizes and force projection radii, so that the limits of the original approach could be assessed. Then, two different sub-models were applied: the classical semi-empirical Glauert correction and a new methodology based on the Lifting Line Theory (LLT), which was recently proposed by Dağ and Sørensen (DS). The latter was here adapted to vertical-axis machines. Eventually, the blade spanwise load profiles coming from the three approaches were assessed and compared, proving the superior performance of the DS model.
Abstract
Aeroelasticity is recognized as the key enabler that allowed for the massive upscaling of wind turbines in the last decade, leading to long, slender, and flexible blades that equip rotors ...with lower specific power and unprecedented energy conversion capabilities. In this study, a selection of case studies with increasing size, specifically the NREL 5 MW, the DTU 10 MW, and the IEA 15 MW Reference Wind Turbines (RWTs) is considered to explore to what extent aeroelastic effects impact the functioning of these rotors. The rotors are not only different in size, with the 15 MW having blades nearly double the length of its 5 MW predecessor, but also in technological level as they belong to different phases of wind turbine development. To evaluate how these machines interact with complex inflow conditions, equivalent boundary conditions have been considered for all simulations, including a realistic atmospheric boundary layer, atmospheric turbulence, and wind shear. To carry out the comparative study, a new aero-servo-elastic module called CALMA is introduced herein, by coupling the engineering software OpenFAST with the CFD software CONVERGE, which solves the wind field and provides inflow conditions for the calculation of loads. Results show how not only aeroelasticity increasingly affects the power performance of the new-generation rotors, but is also a key driver of structural design, for example allowing for an alleviation of 1P load variations due to a sheared and turbulent inflow, thus proving that the inclusion of aeroelasticity in the design process is a key enabler to current and future upscaling trends.
Abstract
Weather routing (WR) systems are widely adopted in the maritime transport since safety of goods and saving of fuel are crucial for shipping companies. However, the need of protecting the ...local coasts and reducing CO
2
emission is making WR attractive even for the market of leisure boats, especially if comfort and safety are also accounted. In the present study, a prototype of a WR system is presented. The developed tool implements the Dijkstra’s algorithm to find the fuel-optimal path while ensuring comfort and safety. A computationally efficient digital twin of a planing boat, based on a 2D+t model and a powertrain mapped model, has been implemented for the estimation of the fuel consumption and the evaluation of vertical accelerations. The methodology for the acquisition of online weather data is presented, together with the strategy for the assessment of comfort and habitability. The application of the WR tool in a typical path of the Tuscan Archipelago shows improvements in fuel usage and comfort even in moderate weather conditions. The influence of each weather variable on fuel efficiency is assessed, highlighting the requirement of an accurate dynamic modelling. Criteria for the graph design are also outlined, showing a wrong estimation of comfort and fuel usage in case of low-definition graphs. The paper proves the potential and the effectiveness of the developed tool, moving toward a greener and more comfortable navigation in local seas.
Abstract
The realization of mathematical, multi-physics models for alkaline electrolyzers is crucial for advancing this technology. Lumped parameter models offer shorter simulation times, compared ...with other approaches, and practical industrial applications. If electrochemical models provide the polarization curve, the hydrogen production rate, and the device efficiency, thermal models solve the equations involving the electrolyte temperature. With a coupled approach, the two models can be linked together by considering the device voltage as temperature dependent. Despite the relevance of such models, few instances of their direct application on existing electrolyzers can be found in the literature, and combined electrochemical-thermal simulations are rare.
This study presents a multi-physics model applied on an alkaline electrolyzer and validated against measurements acquired on a dedicated experimental test bench.
The physics-based model accurately predicts the polarization curve, exhibiting a high precision match with experimental data. Additionally, it identifies material or geometrical imperfections in the electrolyzer, allowing for optimization in the design phase. The thermal model successfully converges to the desired electrolyte temperature of 72 °C under stationary conditions. Additional transient simulations demonstrate an average deviation of only 0.25% compared to the measured temperature trend. Finally, a sensitivity analysis explores the coupling of the electrolyzer with a wind turbine under different wind conditions.
The study showcases the effectiveness of the coupled electrochemical-thermal model in predicting electrolyzer performance, with a direct application to a wind turbine power output.
Simulation methods ensuring a level of fidelity higher than that of the ubiquitous Blade Element Momentum theory are increasingly applied to VAWTs, ranging from Lifting-Line methods, to Actuator Line ...or Computational Fluid Dynamics (CFD). The inherent complexity of these machines, characterised by a continuous variation of the angle of attack during the cycloidal motion of the airfoils and the onset of many related unsteady phenomena, makes nonetheless a correct estimation of the actual aerodynamics extremely difficult. In particular, a better understanding of the actual angle of attack during the motion of a VAWT is pivotal to select the correct airfoil and functioning design conditions. Moving from this background, a high-fidelity unsteady CFD model of a 2-blade H-Darrieus rotor was developed and validated against unique experimental data collected using Particle Image Velocimetry (PIV). In order to reconstruct the AoA variation during one rotor revolution, three different methods-detailed in the study-were then applied to the computed CFD flow fields. The resulting AoA trends were combined with available blade forces data to assess the corresponding lift and drag coefficients over one rotor revolution and correlate them with the most evident flow macro-structures and with the onset of dynamic stall.
The analysis of wind turbine wakes and their interaction with other machines installed in the same array or park has become a key element in the current wind energy research. If in case of siting of ...large rotors for energy production the use of high-fidelity CFD simulations is well established yet, there is still a lack of knowledge in the analysis of proper wind turbine siting for small wind turbines, which are typically installed in quite complex environments. The present study provides the description and validation of a hybrid BEM-CFD model for the analysis of wind turbine performance and wake structure. With respect to similar existing methods, the proposed one includes a specific correction of turbulence parameters able to make it compatible for use in combination with the standard settings of the turbulence models needed to properly describe the wind profile in the urban environment. The model was then used to carry out a demonstrative sensitivity analysis on the proper siting of a small wind turbine in the rooftop of a typical tall building in a densely built environment.
Small Darrieus Vertical-Axis Wind Turbines (VAWTs) are presently seen as a relevant research topic for the wind energy community, since they are thought to perform better than horizontal-axis rotors ...in the complex and highly-turbulent flow typical of the urban environment. Indeed, a preliminary wind tunnel test campaign on a H-Darrieus VAWT showed a significant increase of the performance for high turbulence levels. The present study analyses in detail the near wake of the turbine in the same turbulent conditions, enabling a better insight on the reasons of this power increase, and on the means to take advantage of it. Near-wake measurements are also benchmarked with a CFD simulation of the entire wake, in order to match the experimental wake measurements with the detachment of flow structures observed in CFD simulations. By doing so, a deeper and useful insight on the reasons why VAWTs perform better in turbulent environments is gained.
Light Detection and Ranging (LiDAR) technology can be a valuable tool for describing and quantifying vegetation structure. However, because of their size, extraction of leaf geometries remains ...complicated. In this study, the intensity data produced by the Terrestrial Laser System (TLS) FARO LS880 is corrected for the distance effect and its relationship with the angle of incidence between the laser beam and the surface of the leaf of a Conference Pear tree (Pyrus commmunis) is established. The results demonstrate that with only intensity, this relationship has a potential for determining the angle of incidence with the leaves surface with a precision of ±5° for an angle of incidence smaller than 60°, whereas it is more variable for an angle of incidence larger than 60°. It appears that TLS beam footprint, leaf curvatures and leaf wrinkles have an impact on the relationship between intensity and angle of incidence, though, this analysis shows that the intensity of scanned leaves has a potential to eliminate ghost points and to improve their meshing.
Computational fluid dynamics codes using the density-based compressible flow formulation of the Navier–Stokes equations have proven to be very successful for the analysis of high-speed flows. ...However, solution accuracy degradation and, for explicit solvers, reduction of the residual convergence rates occur as the local Mach number decreases below the threshold of 0.1. This performance impairment worsens remarkably in the presence of flow reversals at wall boundaries and unbounded high-vorticity flow regions. These issues can be resolved using low-speed preconditioning, but there exists an outstanding problem regarding the use of this technology in the strongly coupled integration of the Reynolds-averaged Navier–Stokes equations and two-equation turbulence models, such as the k−ω shear stress transport model. It is not possible to precondition only the RANS equations without altering parts of the governing equations, and there did not exist an approach for preconditioning both the RANS and the SST equations. This study solves this problem by introducing a turbulent low-speed preconditioner of the RANS and SST equations that does not require any alteration of the governing equations. The approach has recently been shown to significantly improve convergence rates in the case of a one-equation turbulence model. The study focuses on the explicit multigrid integration of the governing equations, but most algorithms are applicable also to implicit integration methods. The paper provides all algorithms required for implementing the presented turbulent preconditioner in other computational fluid dynamics codes. The new method is applicable to all low- and mixed-speed aeronautical and propulsion flow problems, and is demonstrated by analyzing the flow field of a Darrieus wind turbine rotor section at two operating conditions, one of which is characterized by significant blade/vortex interaction. Verification and further validation of the new method is also based on the comparison of the results obtained with the developed density-based code and those obtained with a commercial pressure-based code.