Coupled FSI Simulations of the Interaction of a Flexible Hydrofoil With Large Scale Unsteady Flows

Author(s):  
Abe H. Lee ◽  
Robert L. Campbell ◽  
Brent A. Craven ◽  
Stephen A. Hambric

Fluid-structure interaction (FSI) effects must be considered when flexible structures are subjected to unsteady flows. Large-scale unsteady flows can excite structural vibrations significantly and cause the fluid flow to be altered by the large amplitude vibrations. In this work, an in-house finite-element structural code FEANL is tightly coupled with the open-source computational-fluid dynamics (CFD) library package OpenFOAM to simulate the interaction of a backward-skewed, flexible hydrofoil with vortical flow structures shed from a large upstream rigid cylinder in the Penn State-ARL 12” water tunnel. To simulate the turbulent flow at a moderate computational cost, hybrid LES-RANS approaches, i.e. Delayed-Detached-Eddy-Simulation (DDES) and k–ω SST-SAS, are used. The hybrid approaches have been widely employed to simulate massively-separated flows at moderately high Reynolds numbers. Both of the turbulence models are used for a coarse mesh CFD-only case (no FSI effects by assuming a rigid structure) to test their capabilities, and the results of the two models are compared. DDES is chosen to simulate a fine mesh CFD-only case to conduct a mesh convergence study, and it is then used for final FSI simulations. The purpose of this work is focused on obtaining computational results; detailed comparisons against experimental data will be made in future work.

Author(s):  
David Forbes ◽  
Gary Page ◽  
Martin Passmore ◽  
Adrian Gaylard

This study is an evaluation of the computational methods in reproducing experimental data for a generic sports utility vehicle (SUV) geometry and an assessment on the influence of fixed and rotating wheels for this geometry. Initially, comparisons are made in the wake structure and base pressures between several CFD codes and experimental data. It was shown that steady-state RANS methods are unsuitable for this geometry due to a large scale unsteadiness in the wake caused by separation at the sharp trailing edge and rear wheel wake interactions. unsteady RANS (URANS) offered no improvements in wake prediction despite a significant increase in computational cost. The detached-eddy simulation (DES) and Lattice–Boltzmann methods showed the best agreement with the experimental results in both the wake structure and base pressure, with LBM running in approximately a fifth of the time for DES. The study then continues by analysing the influence of rotating wheels and a moving ground plane over a fixed wheel and ground plane arrangement. The introduction of wheel rotation and a moving ground was shown to increase the base pressure and reduce the drag acting on the vehicle when compared to the fixed case. However, when compared to the experimental standoff case, variations in drag and lift coefficients were minimal but misleading, as significant variations to the surface pressures were present.


2021 ◽  
Vol 157 (A1) ◽  
Author(s):  
Z Q Leong ◽  
D Ranmuthugala ◽  
I Penesis ◽  
H D Nguyen

Computational Fluid Dynamics (CFD) simulations using Reynolds Averaged Navier-Stokes (RANS) equations are increasingly adopted as an analysis tool to predict the hydrodynamic coefficients of underwater vehicles. These simulations have shown to offer both a high degree of accuracy comparable to experimental methods and a greatly reduced computational cost compared to Large Eddy Simulation (LES) and Direct Numerical Simulation (DNS). However, one of the major challenges faced with CFD simulations is that the results can vary greatly depending on the numerical model settings. This paper uses the DARPA SUBOFF hull form undergoing straight-line and rotating arm manoeuvres at different drift angles to analyse the hydrodynamic forces and moments on the vehicle against experimental data, showing that the selection of the boundary conditions and turbulence models, and the quality of the mesh model can have a considerable and independent effect on the computational results. Comparison between the Baseline Reynolds Stress Model (BSLRSM) and Shear Stress Transport with Curvature Correction (SSTCC) were carried out for both manoeuvres, showing that with a sufficiently fine mesh, appropriate mesh treatment, and simulation conditions matching the experiments; the BSLRSM predictions offer good agreement with experimental measurements, while the SSTCC predictions are agreeable with the longitudinal force but fall outside the experimental uncertainty for both the lateral force and yawing moment.


Author(s):  
Chih-Hua Wu ◽  
Shengwei Ma ◽  
Chang-Wei Kang ◽  
Teck-Bin Arthur Lim ◽  
Rajeev Kumar Jaiman ◽  
...  

Bluff body structures exposed to ocean current can undergo vortex-induced motion (VIM) for certain geometric and physical conditions. Recently, the study of VIM has been gaining attention for many engineering applications, in particular offshore structures such as buoys, FPSOs, semi-submersibles, Spars and TLPs. The present work is a part of a systematic effort to investigate the VIM response of multi-columns floating platform. In real sea condition, floating platforms are in high Reynolds numbers region and flow patterns around structures are turbulent in nature. For the purpose of investigating and simulating accurately the nonlinear dynamic processes of vortex shedding, transport and wake interactions with the bluff body, the fundamental study of VIM around a square column at moderate Reynolds numbers (1500 ≤ Re ≤ 14000) is firstly investigated. In the present work, the transient flow pattern around a free vibrating square cylinder at moderate Reynolds numbers is numerically investigated by an open source CFD toolbox, OpenFOAM. Good consistency and agreement are found between the present numerical findings and that of experiments. The cylinder, with a blockage area of 4.2%, is mounted on an elastic support for free vibration in the transverse direction. Hybrid RANS-LES turbulence models are considered here for accurate prediction of massively separated turbulent wake flow while maintaining the reasonable computational cost. Three hybrid turbulence models, the DDES (Delayed Detached Eddy Simulation, the k-ω SST-DES (Detached Eddy Simulation), and the k–ω SST-SAS (Scale Adaptive Simulation), are studied and their results are compared with the reported experimental measurements. It is shown that the result of simulation with the k–ω SST-SAS model is closer to the reported literature than the other two and therefore, subsequently adopted for all the simulations of our study in this paper. The scaling effect of cylinder length in the spanwise direction is also studied with the objective to reduce the computational cost. From the comparison with the recent experimental measurements, the discrepancy between the present simulations of reducing cylinder length and the experiment increases only when Re ≥ 4000. This might stem from the increase in wavelength of some vortex shedding modes and turbulence intensity variation in the spanwise direction near the cylinder as Re ≥ 4000. The detailed flow patterns, 3D vortex structures (based on Q-criterion) and vortex-shedding modes are presented in this work as well.


2017 ◽  
Vol 139 (4) ◽  
Author(s):  
Abe H. Lee ◽  
Robert L. Campbell ◽  
Brent A. Craven ◽  
Stephen A. Hambric

Fluid–structure interaction (FSI) is investigated in this study for vortex-induced vibration (VIV) of a flexible, backward skewed hydrofoil. An in-house finite element structural solver finite element analysis nonlinear (FEANL) is tightly coupled with the open-source computational fluid dynamics (CFD) library openfoam to simulate the interaction of a flexible hydrofoil with vortical flow structures shed from a large upstream rigid cylinder. To simulate the turbulent flow at a moderate computational cost, hybrid Reynolds-averaged Navier–Stokes–large eddy simulation (RANS–LES) is used. Simulations are first performed to investigate key modeling aspects that include the influence of CFD mesh resolution and topology (structured versus unstructured mesh), time-step size, and turbulence model (delayed-detached-eddy-simulation and k−ω shear stress transport-scale adaptive simulation). Final FSI simulations are then performed and compared against experimental data acquired from the Penn State-ARL 12 in water tunnel at two flow conditions, 2.5 m/s and 3.0 m/s, corresponding to Reynolds numbers of 153,000 and 184,000 (based on the cylinder diameter), respectively. Comparisons of the hydrofoil tip-deflections, reaction forces, and velocity fields (contours and profiles) show reasonable agreement between the tightly coupled FSI simulations and experiments. The primary motivation of this study is to assess the capability of a tightly coupled FSI approach to model such a problem and to provide modeling guidance for future FSI simulations of rotating propellers in crashback (reverse propeller operation).


Author(s):  
Juraj Králik

Abstract Aim of this contribution is to provide insight view into analysis focused on obtaining external pressure coefficients on isolated two storey low-rise building with 15° elevation gable roof using Computer Fluid Dynamics simulation and these are compared to values that offering Eurocodes. Final Volume Model consisting of polyhedral mesh will be used for analysis with two different turbulence models. Mesh was created with respect to y+ parameter, where desired value was below one which leads us to fine mesh type. Secondary aim of this contribution is to compare performance of selected turbulence models. For this purpose were chosen Detached Eddy Simulation and Large Eddy Simulation which are part of the Scale Resolving Simulation turbulence models.


2021 ◽  
Author(s):  
Anup Zope ◽  
Avery Schemmel ◽  
Xiao Wang ◽  
Shanti Bhushan ◽  
Prashant Singh ◽  
...  

Abstract In this study, we have assessed performance of URANS model, various hybrid RANS/LES turbulence models such as detached eddy simulation, Nichols-Nelson HRLES model, dynamic HRLES (DHRL) model, as well as LES for two classes of problems: (a) heat transfer due to subsonic swirling flow subjected to a sudden expansion leading to cylindrical chamber, and (b) flow separation due to oblique shock wave-turbulent boundary layer interaction (STBLI). The results are assessed using the heat transfer characteristics, separation and reattachment characteristics, and capability to predict flow unsteadiness. The study indicates that URANS can predict large scale flow features reasonably well. However, it fails to resolve turbulence. PANS improves TKE prediction, hence, improves heat transfer prediction. Among the hybrid RANS/LES models, DHRL coupled with ILES is capable of providing accurate prediction of flow separation/reattachment characteristics for boundary layer flows. For free-shear dominated flows, implicit LES performs better compared to the explicit LES models.


Aerospace ◽  
2019 ◽  
Vol 6 (8) ◽  
pp. 86 ◽  
Author(s):  
Amne ElCheikh ◽  
Michel ElKhoury

Numerical simulations are crucial for fast and accurate estimations of the flow characteristics in many engineering applications such as atmospheric boundary layers around buildings, external aerodynamics around vehicles, and pollutant dispersion. In the simulation of flow over urban-like obstacles, it is crucial to accurately resolve the flow characteristics with reasonable computational cost. Therefore, Large Eddy Simulations on non-uniform grids are usually employed. However, an undesirable accumulation of energy at grid-refinement interfaces was observed in previous studies using non-uniform grids. This phenomenon induced oscillations in the spanwise velocity component, mainly on fine-to-coarse grid interfaces. In this study, the two challenging test cases of flow over urban-like cubes and flow over a 3-D circular cylinder were simulated using three different scale-resolving turbulence models. Simulations were performed on uniform coarse and fine grids on one hand, and a non-uniform grid on the other, to assess the effect of mesh density and mesh interfaces on the models’ performance. Overall, the proposed One-Equation Scale-Adaptive Simulation (One-Equation SAS) showed the least deviation from the experimental results in both tested cases and on all grid sizes and types when compared to the Shear Stress Transport-Improved Delayed Detached Eddy Simulation (IDDES) and the Algebraic Wall-Modeled Large Eddy Simulation (WMLES).


Fluids ◽  
2019 ◽  
Vol 4 (3) ◽  
pp. 148 ◽  
Author(s):  
Chunhui Zhang ◽  
Charles Patrick Bounds ◽  
Lee Foster ◽  
Mesbah Uddin

In today’s road vehicle design processes, Computational Fluid Dynamics (CFD) has emerged as one of the major investigative tools for aerodynamics analyses. The age-old CFD methodology based on the Reynolds Averaged Navier–Stokes (RANS) approach is still considered as the most popular turbulence modeling approach in automotive industries due to its acceptable accuracy and affordable computational cost for predicting flows involving complex geometries. This popular use of RANS still persists in spite of the well-known fact that, for automotive flows, RANS turbulence models often fail to characterize the associated flow-field properly. It is even true that more often, the RANS approach fails to predict correct integral aerodynamic quantities like lift, drag, or moment coefficients, and as such, they are used to assess the relative magnitude and direction of a trend. Moreover, even for such purposes, notable disagreements generally exist between results predicted by different RANS models. Thanks to fast advances in computer technology, increasing popularity has been seen in the use of the hybrid Detached Eddy Simulation (DES), which blends the RANS approach with Large Eddy Simulation (LES). The DES methodology demonstrated a high potential of being more accurate and informative than the RANS approaches. Whilst evaluations of RANS and DES models on various applications are abundant in the literature, such evaluations on full-car models are relatively fewer. In this study, four RANS models that are widely used in engineering applications, i.e., the realizable k - ε two-layer, Abe–Kondoh–Nagano (AKN) k - ε low-Reynolds, SST k - ω , and V2F are evaluated on a full-scale passenger vehicle with two different front-end configurations. In addition, both cases are run with two DES models to assess the differences between the flow predictions obtained using RANS and DES.


2020 ◽  
Vol 142 (12) ◽  
Author(s):  
David Romera ◽  
Roque Corral

Abstract This paper presents an efficient method of approximating unsteady flows using a block-wise discrete spatial Fourier series for the modeling of three-dimensional non-axisymmetric flows without making any hypothesis about its temporal periodicity. The method aims at capturing the long-wavelength flow patterns that are present in many unsteady problems of industrial interest, such as compressor stability, with a drastic reduction in computational resources. The method is intended to be used to compute flows exhibiting large-scale instabilities and where the fundamental frequency of the problem is not known beforehand. The approach discretizes the domain using a finite number of blocks or passages, where the flow variables at the supposedly periodic boundaries are continuously updated using the spatial Fourier coefficients of a uniformly spaced set of reduced-passage domains. The NASA rotor 67 under the effect of distorted inflow conditions has been used as verification case to demonstrate the effectiveness and viability of the method. The comparison between the passage-spectral method and the full-annulus solution shows that sound solutions can be obtained with a low number of harmonics. The new method has also been applied to investigate the rotating stall inception of the NASA rotor 67 for distorted inlet flows near stall operating conditions. The method is shown to accurately reproduce the full-annulus solution with a few spatial harmonics, capturing the characteristic features of the complex flow induced by the tip leakage vortex breakdown. The computational cost in this application has been reduced by a factor of between three and seven. This number heavily depends on the ratio between the number of retained harmonics to the number of blades.


Author(s):  
David Romera ◽  
Roque Corral

Abstract This paper presents an efficient method of approximating unsteady flows using a blockwise discrete spatial Fourier series for the modeling of three-dimensional non-axisymmetric flows without making any hypothesis about its temporal periodicity. The method aims at capturing the long wavelength flow patterns which are present in many unsteady problems of industrial interest, such as compressor stability, with a drastic reduction in computational resources. The method is intended to be used to compute flows exhibiting large-scale instabilities and where the fundamental frequency of the problem is not known beforehand. The approach discretizes the domain using a finite number of blocks or passages, where the flow variables at the supposedly periodic boundaries are continuously updated using the spatial Fourier coefficients of a uniformly spaced set of reduced-passage domains. The NASA rotor 67 under stall conditions has been used as verification validation case to demonstrate the effectiveness and viability of the proposed modeling strategy. The comparison between the solutions obtained with the discrete Fourier series and the full-annulus solution shows that accurate solutions can be obtained with a low number of harmonics. The new method has been applied to investigate the rotating stall inception of the NASA rotor 67 for clean and distorted inlet flow near stall operating conditions. The method is shown to accurately reproduce the full-annulus solution with a few spatial harmonics, capturing the characteristic features of the complex flow induced by the tip leakage vortex breakdown. The computational cost in this application has been reduced by a factor of between three and seven, although this number heavily depends on the ratio between the number of retained harmonics and the number of blades.


Sign in / Sign up

Export Citation Format

Share Document