2D Numerical Study on Wake Scenarios for a Flapping Foil

2021 ◽  
Author(s):  
Hui-li Xu ◽  
Marilena Greco ◽  
Claudio Lugni

Abstract Fishes are talented swimmers. Depending on the propulsion mechanisms many fishes can use flapping tails and/or fins to generate thrust, which seems to be connected to the formation of a reverse von Kármán wake. In the present work, the flow past a 2D flapping foil is simulated by solving the incompressible Navier-Stokes equations in the open-source OpenFOAM platform. A systematic study by varying the oscillating frequency, peak-to-peak amplitude and Reynolds number has been performed to analyze the transition of vorticity types in the wake as well as drag-thrust transition. The overset grid method is used herein to allow the pitching foil to move without restrictions. Spatial convergence tests have been carried out with respect to grid resolution and the size of overset mesh domain. Numerical results are compared with available experimental data and discussed. The results show that the adopted methodology can be well applied to simulate large amplitude motions of the flapping foil. The transitions in the types of wake are consistent with the benchmark experimental data, and the drag-thrust transition of the pitching foil does not coincide with von Kármán (vK)-reverse von Kármán (reverse-vK) wake transition and it is highly dependent on the Reynolds number.

1994 ◽  
Vol 116 (4) ◽  
pp. 202-208 ◽  
Author(s):  
K. Nakajima ◽  
Y. Kallinderis ◽  
I. Sibetheros ◽  
R. W. Miksad ◽  
K. Lambrakos

A numerical study of the nonlinear and random behavior of flow-induced forces on offshore structures and experimental verification of the results are presented. The numerical study is based on a finite-element method for the unsteady incompressible Navier-Stokes equations in two dimensions. The momentum equations combined with a pressure correction equation are solved employing fourth-order artificial dissipation with a nonstaggered grid, instead of the more commonly used staggered meshes. The solution is advanced in time with a combined explicit and implicit marching scheme. Emphasis is placed on study of reversing flows around a cylinder. Comparisons with experimental data evaluate accuracy and robustness of the method.


Author(s):  
S. Bhattacharyya ◽  
D. K. Maiti

Numerical study on the wake behind a square cylinder placed parallel to a wall has been made. Flow has been investigated in the laminar Reynolds number (based on the cylinder length) range. We have studied the flow field for different values of the non-dimensional gap length between cylinder and the wall. The case when the cylinder is placed on the wall has also been considered. The governing unsteady Navier-Stokes equations are discretised through the finite volume method on staggered grid system. A SIMPLER type of algorithm has been used to compute the discretised equations iteratively. Vortex shedding has been found to be influenced by the wall. Vortex shedding suppression occurs beyond a critical value of the gap length. Due to the shear, the drag experienced by the cylinder is found to increase with the reduction of gap length. The flow is found to be steady when the cylinder is placed on the wall at a range of Reynolds number.


1992 ◽  
Vol 114 (1) ◽  
pp. 14-25 ◽  
Author(s):  
R. X. Dai ◽  
Q. Dong ◽  
A. Z. Szeri

In this numerical study of the approximations that led Reynolds to the formulation of classical Lubrication Theory, we compare results from (1) the full Navier-Stokes equations, (2) a lubrication theory relative to the “natural,” i.e., bipolar, coordinate system of the geometry that neglects fluid inertia, and (3) the classical Reynolds Lubrication Theory that neglects both fluid inertia and film curvature. By applying parametric continuation techniques, we then estimate the Reynolds number range of validity of the laminar flow assumption of classical theory. The study demonstrates that both the Navier-Stokes and the “bipolar lubrication” solutions converge monotonically to results from classical Lubrication Theory, one from below and the other from above. Furthermore the oil-film force is shown to be invariant with Reynolds number in the range 0 < R < Rc for conventional journal bearing geometry, where Rc is the critical value of the Reynolds number at first bifurcation. A similar conclusion also holds for the off-diagonal components of the bearing stiffness matrix, while the diagonal components are linear in the Reynolds number, in accordance with the small perturbation theory of DiPrima and Stuart.


Author(s):  
Alexander I. Leontiev ◽  
Sergey A. Isaev ◽  
Nikolai V. Kornev ◽  
Yaroslav Chudnovsky ◽  
Egon Hassel

The paper presents a comprehensive analysis of conditions for numerical simulation and physical modeling of convective heat transfer in the vicinity of dimpled surface relief. Contradictory results, unreasonable assumptions, and non-justified conclusions are marked. Based on the analysis of physical experiments the correlation between the predictions and measured data is discussed. Detailed numerical study of turbulent air flow and heat transfer in the narrow channel with three types of dimples (spherical, conic and oval) was carried out. Various mathematical and discrete models, including, those based on solving Reynolds-averaged Navier-Stokes equations (RANS/URANS-SST), and also adaptive scale models (SAS-SST) are compared. The influence of flow parameters (Reynolds number) and geometric sizes (dimple diameter, depth, radius of rounding off of an edge, channel width and height) on local and integral characteristics of flow and heat transfer (total heat output and hydraulic losses) is determined. Special attention is given to reorganizing vortex structures and flow regime (with periodic fluctuations) with increasing relative dimple depth and Reynolds number. For the first time the influence of the scale factor of a constant cross-section channel is detailed. Thermal-hydraulic characteristics of various dimpled reliefs are compared, and the advantage of an oval dimple over a spherical one is shown.


1984 ◽  
Vol 106 (3) ◽  
pp. 591-596 ◽  
Author(s):  
R. S. Amano

A numerical study is reported on the flow and heat transfer in the channel with two right-angled bends. The modified hybrid scheme was employed to solve the steady full Navier-Stokes equations with the energy equation. The computations were performed for different step heights created in a long channel. The local heat transfer rate along the channel wall predicted by employing the present numerical model showed good agreement with the experimental data. The behavior of the flow and the heat transfer were investigated for the range of Reynolds number between 200 and 2000 and for step height ratios H/W = 1, 2, and 3. Finally, the correlation of the average Nusselt number in such channels as a function of Reynolds number is postulated.


Author(s):  
William Denner Pires Fonseca ◽  
Rafael Rosario Da Silva ◽  
Reinaldo Marcondes Orselli ◽  
Adson Agrico De Paula ◽  
Ricardo Galdino da Silva

In this work, a numerical study of flow around an airfoil with wavy leading edge is presented at a Reynolds number of 3X106. The flow is resolved by considering the RANS (Reynolds Average Navier-Stokes)equations. The baseline geometry is based on the NACA 0021 profile. The wavy leading edge has an amplitude of 3% and wavelength of 11%, both with respect to the airfoil chord. Cases without and with wavy leadingedges are simulated and compared. Initially, studies of the numerical sensitivity with respect to the obtained results, considering aspects such as turbulence modeling and mesh refinement, are carried out as well as bycomparison with corresponding results in the literature. Numerical data such as pressure distribution, shear stress lines on the wing surface, and aerodynamics coefficients are used to describe and investigate the flowfeatures around the wavy leading airfoil. Comparisons between the straight leading edge and the wavy leading edge cases shows an increase of the maximum lift coefficient as well as stall angle for the wavy leading edge configuration. In addition, at an angle of attack near the stall, the present numerical results shows an increase of the drag coefficient with the wavy leading edge airfoil when compared with the corresponding straight leading edge case.


Author(s):  
Bogdan Iwanowski ◽  
Marc Lefranc ◽  
Rik Wemmenhove

Numerical study of liquid dynamics in an LNG tank is presented. The available data from large scale (1:10) sloshing experiments of 2D section of an LNG carrier reveal large scatter in recorded values of peak pressures. The experimental data is analysed from statistical point of view in order to obtain distributions of the pressure peaks. Then the entire experimental data record is reproduced numerically by CFD simulations and it is shown that pressure peaks obtained numerically display scatter of values as well. A statistical description of the numerically obtained record is provided and compared with description derived from the experimental data. The applied CFD code ComFLOW solves Navier-Stokes equations and uses an improved Volume of Fluid (iVOF) method to track movement of fluid’s free surface. Two different fluid models, single-phase (liquid+void) and two-phase (liquid+compressible gas) can be applied, the latter model being capable of simulating bubbles and gas entrapped in liquid. For low tank filling rate discussed in the paper (10%) the single-phase approach is sufficient. Comparison of statistical properties of experimental and numerical records is offered.


1987 ◽  
Vol 177 ◽  
pp. 133-166 ◽  
Author(s):  
John Kim ◽  
Parviz Moin ◽  
Robert Moser

A direct numerical simulation of a turbulent channel flow is performed. The unsteady Navier-Stokes equations are solved numerically at a Reynolds number of 3300, based on the mean centreline velocity and channel half-width, with about 4 × 106 grid points (192 × 129 × 160 in x, y, z). All essential turbulence scales are resolved on the computational grid and no subgrid model is used. A large number of turbulence statistics are computed and compared with the existing experimental data at comparable Reynolds numbers. Agreements as well as discrepancies are discussed in detail. Particular attention is given to the behaviour of turbulence correlations near the wall. In addition, a number of statistical correlations which are complementary to the existing experimental data are reported for the first time.


1985 ◽  
Vol 157 ◽  
pp. 471-492 ◽  
Author(s):  
C.-Y. Lai ◽  
K. R. Rajagopal ◽  
A. Z. Szeri

In this paper we generalize the von Kármán solution for flow above a single rotating disk, to include non-axisymmetric solutions. These solutions contain an arbitrary parameter; for zero value of the parameter the asymmetric flow degenerates into the classical von Kármán solution. Thus the classical solution is never isolated when considered within the scope of the full Navier–Stokes equations; there are asymmetric solutions in every neighbourhood of the von Kármán solution. Calculations are reported here for s = 0, 0.02 and 0.06, where s represents the ratio of angular velocity of the fluid at infinity to the angular velocity of the disk. A subset of the solutions obtained here corresponds to flow induced by the rotation of a disk when the latter is placed in a fluid that is moving with a constant uniform velocity.


Processes ◽  
2021 ◽  
Vol 9 (5) ◽  
pp. 841
Author(s):  
Yuzhen Jin ◽  
Huang Zhou ◽  
Linhang Zhu ◽  
Zeqing Li

A three-dimensional numerical study of a single droplet splashing vertically on a liquid film is presented. The numerical method is based on the finite volume method (FVM) of Navier–Stokes equations coupled with the volume of fluid (VOF) method, and the adaptive local mesh refinement technology is adopted. It enables the liquid–gas interface to be tracked more accurately, and to be less computationally expensive. The relationship between the diameter of the free rim, the height of the crown with different numbers of collision Weber, and the thickness of the liquid film is explored. The results indicate that the crown height increases as the Weber number increases, and the diameter of the crown rim is inversely proportional to the collision Weber number. It can also be concluded that the dimensionless height of the crown decreases with the increase in the thickness of the dimensionless liquid film, which has little effect on the diameter of the crown rim during its growth.


Sign in / Sign up

Export Citation Format

Share Document