URANS Approach Research Articles

Numerical study of a fish swimming in hydrokinetic turbine wake

The environmental effects of hydrokinetic turbines are still under investigation, reflecting the emerging status of this technology. This study investigates the interaction between hydrokinetic rotor wakes and fish swimming, revealing insights into fish biomechanics in complex flows and assessing the environmental implications of marine energy solutions. We conducted numerical simulations with the URANS approach and k−ω−SST turbulence closure model to predict three-dimensional turbulent flow in the OpenFOAM software. The hydrokinetic rotor wake was simulated employing the actuator line method, providing a computationally efficient alternative to full geometry simulations. For accurate replication of the motion of a fish-like tuna (Thunnus atlanticus), dynamic adaptive mesh discretization was employed. The results offer a comparative analysis of fish swimming performance within the wake rotor, particularly when immersed in the tip blade vortex, contrasted with scenarios where fish swim in undisturbed flow conditions. The analysis encompasses three-dimensional wake structures, force generation, efficiency, and equilibrium states (balancing drag and thrust) across varying Swimming numbers (Sw). Key findings include the enhanced attachment of the leading-edge vortex due to the caudal fin’s interaction with the tip blade vortex, resulting in improved auto-propulsive force production; a reduced tail stride frequency observed in fish swimming downstream of the rotor to achieve longitudinal force balance compared to unperturbed flow; and transverse hydrodynamic forces pushing fish radially away from the wake’s influence zone, potentially mitigating the risk of collision with turbine blades.

Numerical Investigation of Dimensionless Parameters in Carangiform Fish Swimming Hydrodynamics.

Research into how fish and other aquatic organisms propel themselves offers valuable natural references for enhancing technology related to underwater devices like vehicles, propellers, and biomimetic robotics. Additionally, such research provides insights into fish evolution and ecological dynamics. This work carried out a numerical investigation of the most relevant dimensionless parameters in a fish swimming environment (Reynolds Re, Strouhal St, and Slip numbers) to provide valuable knowledge in terms of biomechanics behavior. Thus, a three-dimensional numerical study of the fish-like lambari, a BCF swimmer with carangiform kinematics, was conducted using the URANS approach with the k-ω-SST transition turbulence closure model in the OpenFOAM software. In this study, we initially reported the equilibrium Strouhal number, which is represented by St∗, and its dependence on the Reynolds number, denoted as Re. This was performed following a power-law relationship of St∝Re(-α). We also conducted a comprehensive analysis of the hydrodynamic forces and the effect of body undulation in fish on the production of swimming drag and thrust. Additionally, we computed propulsive and quasi-propulsive efficiencies, as well as examined the influence of the Reynolds number and Slip number on fish performance. Finally, we performed a vortex dynamics analysis, in which different wake configurations were revealed under variations of the dimensionless parameters St, Re, and Slip. Furthermore, we explored the relationship between the generation of a leading-edge vortex via the caudal fin and the peak thrust production within the motion cycle.

Highly resolved LES and URANS computations of a differentially heated square cavity

This investigation focusses on the use of two wall-resolved LES approaches for the simulation of the buoyancy driven flow in a differentially heated square cavity, aiming to advance understanding of the flow phenomena involved. The study also assesses the predictive accuracy and cost-effectiveness of LES for natural convection flows. This assessment is pursued through comparisons of the resulting predictions, and also of their computational requirements, with those for DNS, produced in an earlier study, and those of a number of URANS models produced during the current investigation. URANS models tested here include both high- and low-Reynolds-number schemes. In the case of the high-Reynolds-number models a recently developed numerical variant of the Analytical Wall Function (AWF) has been employed to account for the effects of near-wall turbulence. Both LES approaches result in predictions that are close to those of the DNS but have considerably more modest resource requirements. The URANS approaches, which have the lowest resource requirements, for key parameters such as the local Nusselt number are not as accurate as the LES approaches and their performance is highly dependent on the turbulence model employed. Surprisingly, the simplest of the low-Re models employed, the Launder–Sharma two-equation model, is reasonably close to the DNS data. The approach that requires the least computational resource, the high-Re model with the recent version of the AWF, also offers a promising URANS route. Simulations were also conducted to investigate the flow behaviour at a Rayleigh number that is an order of magnitude higher than that of the available DNS. At such high Rayleigh number, at which the phenomena are of close relevance to cooling loops of nuclear reactors, comparisons are performed between the LES and the high-Re AWF combination. The findings show that as the Rayleigh number increases, the turbulence levels become stronger, but these are confined to more refined structures within the thinner near-wall boundary layers.

Numerical investigation of active flow control using zero-net-mass-flux jets around a high-lift morphing cambered wing-flap system

PurposeThis study aims to investigate the physical mechanisms of the use of active flow control (AFC) for a high-lift wing-flap configuration.Design/methodology/approachBy means of high-fidelity numerical simulations, the flow dynamics around a high-lift wing-flap system at high Reynolds number (Re/c = 4.6 million) is studied. Adapted turbulence models based on the URANS approach are used to capture the flow separation and the subsequent development of coherent structures. The present study focuses on the use of AFC using a synthetic jet known as zero-net-mass-flux (ZNMF) using the blowing–suction approach. Different parameters (geometry, frequency and velocity) of a ZNMF placed at the cambered flap’s chord are optimized to obtain the most efficient parameter settings to suppress the flow separation.FindingsA synthetic jet with the optimal shape and orientation enforces the flow reattachment on the wing-flap surface. This leads to an improvement of the aerodynamic performance of the system. The wake thickness was reduced by 30%, and an increase of 17.6% in lift-to-drag ratio was obtained. Concerning the ZNMF location, they should be installed upstream of the separation point to achieve the best performance.Originality/valueThe effectiveness of ZNMF devices integrated on a high-lift wing-flap configuration was studied in real flight conditions at high Reynolds number. A detailed analysis of the wake dynamics explains how AFC forces the reattachment of the boundary layer and attenuates the predominant wake instabilities up to −20 dB.

Turbulence-induced vibrations prediction through use of an anisotropic pressure fluctuation model

In nuclear fuel rod bundles, turbulence-induced pressure fluctuations caused by an axial flow can create small but significant vibrations in the fuel rods, which in turn can cause structural effects such as material fatigue and fretting wear. Fluid-structure interaction simulations can be used to model these vibrations, but for affordable simulations based on the URANS approach, a model for the pressure fluctuations must be utilised. Driven by the goal to improve the current state-of-the-art pressure fluctuation model, AniPFM (Anisotropic Pressure Fluctuation Model) was developed. AniPFM can model velocity fluctuations based on anisotropic Reynolds stress tensors, with temporal correlation through the convection and decorrelation of turbulence. From these velocity fluctuations and the mean flow properties, the pressure fluctuations are calculated. The model was applied to several test cases and shows promising results in terms of reproducing qualitatively similar flow structures, as well as predicting the root-mean-squared pressure fluctuations. While further validation is being performed, the AniPFM has already demonstrated its potential for affordable simulations of turbulence-induced vibrations in industrial nuclear applications.

Atmosphere homogenization induced by vertical jets in large enclosures

To theoretically investigate atmosphere mixing in an enclosure, three experiments of interaction of a vertical air jet with a helium-rich layer performed in SPARC, PANDA and MiniPanda experimental facilities, were simulated using the URANS approach. For the purposes of this work, firstly, a local Froude number is introduced, which can be directly used in computational fluid dynamics calculations to specify the regions with possible occurrence of Kelvin-Helmholtz instability. Namely, it is shown that in these regions, turbulent diffusion is underpredicted and turbulence model fails to replicate some phenomena. Secondly, a model for dynamically prescribe turbulent Schmidt and Prandtl numbers is proposed. The proposed model improves the results in isothermal case with possible occurrence of Kelvin-Helmholtz instability, while its effects are negligible in cases with higher temperature jets with lower local Froude number values, where the conditions for the Kelvin-Helmholtz instability are not established.

Flow instability in a volute-free centrifugal fan subjected to non-axisymmetric pre-swirl flow from upstream bended inflow tube

The characteristics of internal flow and performance of a centrifugal fan is greatly dependent on the inflow pattern. As the fan is subjected to incoming flow from an upstream tube, the size and geometry of the tube affect the three-dimensional motion of local flow and possibly degrades the aerodynamic performance of the fan. In this work, we performed a numerical investigation on the internal flow in a centrifugal fan subjected to incoming flow from an upstream bended inflow tube of various radii using the steady and unsteady Reynolds-averaged Navier-Stokes (RANS and URANS) simulation approaches. The effects of the non-axisymmetric pre-swirl flow generated due to the curvature of the bended inflow tube are demonstrated by analyzing the internal flow characteristics of the fan, including the spatial distributions and temporal variations of pressure field and streamlines, pressure fluctuations in the upstream tube, the inflow and outflow sections of the impeller, and the circumferential distributions of velocity and pressure in the impeller. The numerical results reveal that as the inflow tube is curved with larger curvature (smaller radius of the bended section), the pre-swirl inflow is strong and deteriorates the static pressure rise and static pressure efficiency of the centrifugal fan more remarkably, and the circumferential non-uniformity of pressure and velocity distributions appears inside of the channels of the fan. As the radius of the bended section increases, the instability of the internal flow gets more pronounced, as represented by the stronger pressure fluctuations at the inflow and outflow sections. The prediction capabilities of RANS and URANS approaches are also analyzed based on the numerical data and we found that the latter is more reliable in predicting the performance of the fan.

Analysis and Numerical Simulation of no Reacting Swirling Flow by Using URANS Approach

In this paper, we investigate the no-reacting swirling flow by using the numerical simulation based to the unsteady Reynolds-averaged Navier-Stokes approach. The numerical simulation was realized by using a computational fluid dynamics CFD code. The governing equations are solved by using the finite volume method with two classical models of turbulence K-epsilon and Shear Stress K-ω. The objective of this paper is therefore to evaluate the performance of the two models in predicting the recirculation zones in a swirled turbulent flow. The current models are validated by comparing the numerical results of the axial, radial and tangential velocities to the experimental data from literature.

Wind Flow Characterisation over a PV Module through URANS Simulations and Wind Tunnel Optical Flow Methods

Despite their simplicity, photovoltaic (PV) modules are often arranged in structures that can be affected by severe and complex wind loads: in this context, the wind flow and the dynamic excitation induced by vortex shedding can introduce unexpected aeroelastic responses. This work introduces a novel wind tunnel application of experimental techniques to address this issue by the use of flow visualisation and video postprocessing, through the optical flow algorithm. Numerical simulations based on unsteady Reynolds-averaged Navier–Stokes (RANS) models are performed and compared against the experimental wind tunnel tests on a PV panel that was also instrumented with pressure taps. A setup with a 65∘ tilt angle was examined because, based on preliminary analyses, it was considered interesting for the free flow–wake transition associated with the dynamic response of the PV panel. The comparison of the experimental and numerical average wind fields supported that the proposed optical flow method was appropriate for characterising the wake of the panel, because there was enough seeding to perform the video postprocessing. Experiments and numerical predictions were compared as regards the average pressure distribution on the panel surfaces, and the average percentage was in the error of 7%; this supports that the URANS method was capable of reproducing the average behaviour of the panel, as well as for the selected configuration, which is particularly challenging. Furthermore, the simulated and measured power spectral densities of the wind speed were compared, and this resulted in the numerical model quite faithfully reproducing the frequency of the peak at 5 m/s, while the error was in the order of 20% for the 10 m/s case; this supports that, despite the URANS approach being affected by well-known critical points regarding the simulation of instantaneous quantities, it can be employed to elaborate information that can be particularly useful for the structural design of the panel. This kind of result can be considered as a first step, obtained with simplified and affordable methods, towards a characterisation of the dynamic behaviour of a PV panel in a real-world setup.

On the Accuracy of uRANS and LES-Based CFD Modeling Approaches for Rotor and Wake Aerodynamics of the (New) MEXICO Wind Turbine Rotor Phase-III

This work presents a comparison study of the CFD modeling with two different turbulence modeling approaches viz. unsteady RANS and LES, on a full-scale model of the (New) MEXICO rotor wind turbine. The main emphasis of the paper is on the rotor and wake aerodynamics. Simulations are carried out for the three wind speeds considered in the MEXICO experiment (10, 15, and 24 ms−1). The results of uRANS and LES are compared against the (New) MEXICO experimental measurements of pressure distributions, axial, radial, and azimuth traverse of three velocity components. The near wake characteristics and vorticity are also analyzed. The pressure distribution results show that the LES can predict the onset of flow separation more accurately than uRANS when the turbine operates in the stall condition. The LES can compute the flow structures in wake significantly better than the uRANS for the stall condition of the blade. For the design condition, the mean absolute error in axial and radial velocity components along radial traverse is less than 10% for both the modeling approaches, whereas tangential component error is less than 2% from the LES approach. The results also reveal that wake recovers faster in the uRANS approach, requiring further research of the far wake region using both CFD modeling approaches.

An URANS Simulation of the Kelvin-Helmholtz Aerodynamic Effect over the Ahmed Body

The reduction of energy consumption of cars is always a significant issue in automotive design. Turbulent flow around a car body is very difficult to simulate accurately due to the complexity of the flow conditions around the body, such as complex flow separation and laminar to turbulent flow transition. In particular, flow over the Ahmed body with a rear angle of 25o is considered a challenging problem for the RANS approach with two-equation turbulence models. In this study, we aim to analyze the Kelvin-Helmholtz instability associated with this flow with a URANS approach. Methodology for utilizing the URANS method is fully discussed. The predicted velocity profiles and drag coefficient are com-pared with experimental results. Three turbulence models, such as the k-ε, k-ω and SST models, are assessed and validated with experimental data. The aim of the study is to evaluate the performance of these models for the study of the Kel-vin-Helmholtz instability over the Ahmed body and for car bodies generally us-ing experimental data for their validations. It is found that the URANS approach with the turbulence models with proper numerical treatment can perform as well as or even better than the LES. And the SST model shows the best performance compared with other turbulence models.

Development of the Wind Tunnel CFD Model for Clarifying Rotor Blade Airfoil Unsteady Aerodynamic Researches

A computational model of a wind tunnel (WT) with a special experimental rig has been developed directly for studying unsteady aerodynamic characteristics of a wing section with a helicopter blade profile. The model was constructed using CFD (Computational Fluid Dynamics) methods based on the URANS approach (Unsteady Reynolds Averaged Navier Strokes Equations). Therefore, the aerodynamic characteristics of the airfoil can be determined, taking into account influence of the WT walls such as walls perforation and various configurations of the experimental rig. Simulation of the flow around the wing section with the SC1095 airfoil in steady and unsteady settings is performed. The flow features in the working section of the WT and the experimental rig are analyzed. For a particular case, the calculation method was validated in a 2D formulation on the basis of available experimental data. The developed model can be used to refine the methods of processing experimental data, taking into account the individual characteristics of the WT and the experimental rig configuration.

UNSTEADY TWO-PHASE GAS-PARTICLE FLOWS IN BLADE CASCADES

Unsteady flows of a gas with solid particles in a system of two plane “rotor–stator” cascades are investigated. The carrier-gas flow is modeled using the Navier–Stokes equations (pseudo DNS approach) and, for comparison, the Reynolds equations (URANS approach) with the Menter SST k–ω model of turbulence. In both cases, the equations are solved using a second-order finite- volume method. In the absence of particle collisions, the admixture motion is modeled using a Lagrangian approach and the colliding particles are modeled by the Monte Carlo method. The feedback effect of the particles on the carrier gas flow is taken into account. The influence of different factors of random nature (particle distribution over sizes, particle scattering after collisions with the blades, particle–particle collisions) on the admixture flow pattern and particle concentration profiles at the outlet of the stator cascade is analyzed.

Simulation of transonic buffet with an automated zonal DES approach

A study of transonic buffet on the NASA Common Research Model at flight Reynolds numbers is presented. The ability of two different hybrid RANS/LES models as well as the URANS approach for resolving three-dimensional buffet motion was evaluated by means of spectral analysis. Automated Zonal DES and URANS simulations show similar results in terms of buffet frequency and spanwise propagation of buffet cells, whereas the Delayed Detached Eddy Simulation results indicate a strong interaction between flow separation and shock motion. The extracted characteristic frequencies which are associated with transonic buffet are located in a range of Sr = 0.2–0.65 for URANS and AZDES and are therefore in accordance with findings from related recent research. Furthermore, the simulation time series were investigated and a structure of spanwise moving buffet cells with varying convection speed and wavelength could be observed.

On the accuracy of predicting cavitation impact loads on marine propellers

Predicting the cavitation impact loads on a propeller surface using numerical tools is becoming essential, as the demand for more efficient designs, stretched to the limit, is increasing. One of the possible design limits is governed by cavitation erosion. The accuracy of estimating such loads, using a URANS approach, has been investigated. We follow the energy balance approach by (Schenke and van Terwisga, 2019), (Schenke et al., 2019), where we take account of the focusing of the potential energy into the collapse center before it is radiated as shock wave energy in the domain. In complex flows, satisfying the total energy balance, when reconstructing the radiated energy, has always been an issue in the past. Therefore, in this study, we investigate different considerations for the vapor reduction rate, in order to minimize the numerical errors, when estimating the local surface impact power. We show that when the vapor volume reduction rate is estimated using the mass transfer source term, then all the energy is conserved and the total energy balance is satisfied. The model is verified on a single cavitating bubble collapse, and it is further validated on a model propeller test case. The obtained surface impact distribution agrees well with the experimental paint test results, illustrating the potential for practical use of our fully conservative method to predict cavitation implosion loads on propeller blades.

Biomechanics of fish swimming in a hydrokinetic turbine wake

In the present paper, the assessment of the hydrodynamics of the fish swimming in undisturbed and altered wake-flow is carried out identifying hydrodynamical forces and flow patterns of the fish swimming wake. URANS approach with k-ω/SST turbulence model are employed combining fish and turbine in the same simulation. Fish motion is realized using dynamically adaptive mesh. The actuator line method is employed to induce the wake of a hydrokinetic turbine, which is a simplified method that requires lower computational cost than full geometry simulations. This work brings a new numerical approach involving fish and turbine wake highlighting that fish swimming in the wake presents higher thrust forces than in the undisturbed flow, due to the x-component to velocity in the wake is lower than free flow velocity.

Numerical Analysis of a Surface Ship in Waves

Abstract Ship's resistance and engine power to sustain ship's speed in seaways are augmented due to complex non-linear interactions between the ship and the ambient sea (waves). Ship designers, in early design stage, use an ad hoc "sea margin" to account for the effects of seaways in selecting propeller and engine. A numerical tool capable of accurately predicting added resistance and power of a ship cruising in waves would greatly help select its powering (margin) requirement and determine the optimal operating point that can maximize the energy efficiency. For seakeeping analysis, strip theory-based methods have long been used. More recently, nonlinear time-domain three-dimensional (3D) panel methods have started being used widely. A more physics-based avenue to seakeeping analysis is offered by coupled solutions of two-phase unsteady Reynolds-Averaged Navier-Stokes equations and six degrees-of-freedom rigid-body motion (RBM) equations. The URANS approach also avails itself of including the effects of propulsors, either explicitly or approximately. By accounting for all the nonlinear effects in hydrodynamic forces and moments and the resulting ship motions, and the effects of fluid viscosity and turbulence, the coupled URANS-RBM method is believed not only able to predict added resistance and speed loss more accurately, but also to provide valuable insights into the physical mechanisms underlying added resistance and power. The objectives of this study are: (1) to validate a coupled URANS-RBM solver developed for high-fidelity prediction of added resistance, speed loss and added power on ships cruising in regular head sea and irregular waves, and (2) to conduct a detailed analysis of the interactions among ship hull, propeller and waves for a 1/49 scaled model of the ONR Tumblehome (ONRT) (Model 5613) in order to shed light on the physical mechanisms leading to added resistance, speed loss and added power. Figure 1 depicts the ONRT self-propellers with two 4-bladed propellers in regular waves. The main flow features such as the free surface, the hub vortices and blade-tip vortices from the propeller, as well as vortices generated by the sonar dome, shafts, shaft brackets and bilge keels are captured.

Numerical-experimental analysis of the slipstream produced by a high speed train

The recent development of high-speed trains over the last decade led to a growing interest in their aerodynamics. A train at full operational speed generates a strong induced airflow that may damage infrastructures near the trackside and endanger people near the rails. Nowadays the tests required for the train homologation are based on full scale measurements of the airspeed taken in specified positions close to the railway line. An important part of the work was initially dedicated to the data analysis of an extensive experimental campaign performed on the Italian high-speed line, in order to obtain reliable results to be used for comparison with CFD results. These ones have been performed on the full-scale train geometry using an URANS approach with a wall treatment based on wall functions, checking mesh and turbulence models dependencies. The full-scale simulations highlighted the great complexity of the problem being the comparison with the experimental measurements performed close to the shear layer of the train boundary layer. Encouraging results are obtained from the numerical analysis on the full-scale train geometry indicating the ability of the URANS coupled with the SST turbulence model in the prediction of both the train induced flow and the turbulent structures around it.

Simulation of turbulent wakes in model wind farm with arbitrary location for wind turbines

Wind energy is an important part of renewable energy sources in the Russian Federation. A new modern wind farm with 28 wind turbines was built on the territory of the Russian Federation, in the Ulyanovsk region in 2017-2018. When designing and operating a wind farm, it is necessary to assess correctly the power of each individual wind turbine. In this paper we used the open SOWFA library based on OpenFOAM 2.4.0 and URANS approach with k-e model of turbulence. The Actuator Line Model for the calculation of the forces from the wind turbines was used in pisoFoamTurbine.ALM solver. To set the boundary conditions at the input, the Precursor method was used. The unstructured grids with different number of cells were built. As a result of the calculation for 14 simulators of wind turbines, the values of the power coefficients Cp and thrust coefficients Ct were also determined. The maximum Cp value was in agreement with the Bets-Zhukovsky law for an ideal wind turbine. The new numerical case with 28 wind turbines was prepared for future simulations. The calculations were performed on a cluster of the UNICFD web-laboratory of ISP RAS. Each calculation was performed on 48-72 cpu cores.

Numerical and experimental study of the flow around a 4:1 rectangular cylinder at moderate Reynolds number

This paper presents the results of investigations into the flow around a rectangular cylinder with a chord-to-depth ratio equal to 4. The studies are performed through wind tunnel dynamic pressure measurements along a cross section combined with two-dimensional Unsteady Reynolds-Averaged Navier-Stokes (urans) and three-dimensional Delayed-Detached Eddy Simulation (ddes). These experimental and numerical studies are complementary and combining them allows a better understanding of the unsteady dynamics of the flow. These studies aim mainly at determining the effects of the rectangle incidence and freestream velocity on the variation of the flow topology and the aerodynamic loads, and at assessing the capability of the industrially affordable urans and ddes approaches to provide a sufficiently accurate estimation of the flow for different incidences. The comparison of experimental and numerical data is performed using statistics and Dynamic Mode Decomposition. It is shown that the rectangular cylinder involves complex separation-reattachment phenomena that are highly sensitive to the Reynolds number. In particular, the time-averaged lift slope c¯lα increases rapidly with the Reynolds number in the range 7.8×103≤Re≤1.9×104 due to the modification of the time-averaged vortex strength, thickness and distance from the surface. Additionally, it is shown that both urans and ddes simulations fail to accurately predict the flow at all the different incidence angles considered. The urans approach is able to qualitatively estimate the spatio-temporal variations of vortices for incidences below the stall angle α=4∘. Nonetheless, urans does not predict stall, while ddes correctly identifies the stall angle observed experimentally.