scholarly journals DynamFluid: Development and Validation of a New GUI-Based CFD Tool for the Analysis of Incompressible Non-Isothermal Flows

Processes ◽  
2019 ◽  
Vol 7 (11) ◽  
pp. 777
Author(s):  
Héctor Redal ◽  
Jaime Carpio ◽  
Pablo A. García-Salaberri ◽  
Marcos Vera

A computational fluid dynamics software (DynamFluid) based on the application of the finite element method with the characteristic-based-split algorithm is presented and validated. The software is used to numerically integrate the steady and unsteady Navier–Stokes equations for both constant-density and Boussinesq non-isothermal flows. Benchmark two-dimensional computations carried out with DynamFluid show good agreement with previous results reported in the literature. Test cases used for validation include (i) the lid-driven cavity flow, (ii) mixed convection flow in a vertical channel with asymmetric wall temperatures, (iii) unsteady incompressible flow past a circular cylinder, and (iv) steady non-isothermal flow past a circular cylinder with negligible buoyancy effects. The new software is equipped with a graphical user interface that facilitates the definition of the fluid properties, the discretization of the physical domain, the definition of the boundary conditions, and the post-processing of the computed velocity, pressure and temperature fields.

1990 ◽  
Vol 220 ◽  
pp. 459-484 ◽  
Author(s):  
H. M. Badr ◽  
M. Coutanceau ◽  
S. C. R. Dennis ◽  
C. Ménard

The unsteady flow past a circular cylinder which starts translating and rotating impulsively from rest in a viscous fluid is investigated both theoretically and experimentally in the Reynolds number range 103 [les ] R [les ] 104 and for rotational to translational surface speed ratios between 0.5 and 3. The theoretical study is based on numerical solutions of the two-dimensional unsteady Navier–Stokes equations while the experimental investigation is based on visualization of the flow using very fine suspended particles. The object of the study is to examine the effect of increase of rotation on the flow structure. There is excellent agreement between the numerical and experimental results for all speed ratios considered, except in the case of the highest rotation rate. Here three-dimensional effects become more pronounced in the experiments and the laminar flow breaks down, while the calculated flow starts to approach a steady state. For lower rotation rates a periodic structure of vortex evolution and shedding develops in the calculations which is repeated exactly as time advances. Another feature of the calculations is the discrepancy in the lift and drag forces at high Reynolds numbers resulting from solving the boundary-layer limit of the equations of motion rather than the full Navier–Stokes equations. Typical results are given for selected values of the Reynolds number and rotation rate.


2015 ◽  
Vol 2015 ◽  
pp. 1-9 ◽  
Author(s):  
Yingnan Fu ◽  
Xizeng Zhao ◽  
Xinggang Wang ◽  
Feifeng Cao

Viscous flow past an upstream in-line forced oscillating circular cylinder with a stationary cylinder downstream at Reynolds number of 100 is investigated using a CIP model. The model is established in a Cartesian coordinate system using a high-order difference method to discretise the Navier-Stokes equations. The fluid-structure interaction is treated as a multiphase flow with fluid and solid phases solved simultaneously. An immersed boundary method is used to deal with the fluid-body coupling. The CFD model is firstly applied to the computation of flow past a fixed circular cylinder for its validation; then flow over two stationary tandem cylinders is investigated and good agreements are obtained comparing with existing ones. Computations are then performed with flow past two tandem cylinders with an upstream in-line oscillating cylinder with a small spacingL=2D. Considerable attention is paid to the spectrum characteristics and vortex modes.


1973 ◽  
Vol 60 (1) ◽  
pp. 105-127 ◽  
Author(s):  
W. M. Collins ◽  
S. C. R. Dennis

An accurate method is described for integrating the Navier-Stokes equations numerically for the time-dependent flow past an impulsively started circular cylinder. Results of integrations over the range of Reynolds numbers, based on the diameter of the cylinder, from 5 to ∞ are presented and compared with previous numerical, theoretical and experimental results. In particular, the growth of the length of the separated wake behind the cylinder has been calculated forR= 40, 100 and 200 and is found to be in very good agreement with the results of recent experimental measurements. The calculated pressure distribution over the surface of the cylinder forR= 500 is also found to be in reasonable agreement with experimental measurements for the caseR= 560.For Reynolds numbers up to 100 the equations were integrated until most of the features of the flow showed a close approximation to steady-state conditions. The results obtained are in good agreement with previous calculations of the steady flow past a circular cylinder. ForR> 100 the integrations were continued until the implicit method of integration broke down by reason of its failure to converge. A secondary vortex appeared on the surface of the cylinder in the caseR= 500, but for higher Reynolds numbers, including the caseR= ∞, the procedure broke down before the appearance of a secondary vortex. In all cases the flow was assumed to remain symmetrical.


2013 ◽  
Vol 730 ◽  
pp. 19-51 ◽  
Author(s):  
Jiten C. Kalita ◽  
Shuvam Sen

AbstractStudies on the$\alpha $- and$\beta $-phenomena, terms coined by Bouard & Coutanceau (J. Fluid Mech., vol. 101, 1980, pp. 583–607) for the flow past an impulsively started circular cylinder, have been confined only to the very early stages of the flow. In this paper, besides making a comprehensive in-depth analysis of these phenomena for a much longer period of time, we report the existence of some tertiary vortex phenomena for the first time, which we term the sub-$\alpha $- and sub-$\beta $-phenomena. The mechanism of unsteady flow separation at high Reynolds numbers for the flow past a circular cylinder developed in the last two decades has been used to understand these flow phenomena. The flow is computed using a recently developed compact finite difference method for the biharmonic form of the two-dimensional Navier–Stokes equations for the range of Reynolds number$500\leq \mathit{Re}\leq 10\hspace{0.167em} 000$. We specifically choose$\mathit{Re}= 5000$to describe the interplay among the primary, secondary and tertiary vortices leading to these interesting vortex dynamics. We also report a$\beta $-like phenomenon which is very similar to the$\beta $-phenomenon, but slightly differs in details. We offer a new perception of the$\alpha $-phenomenon by defining its existence in a strong and weak sense along with a clearer characterization of the$\beta $-phenomenon. Apart from numerical computation, a detailed theoretical characterization using topological aspects of the boundary layer separation leading to the secondary and tertiary vortex phenomena has also been carried out. We compare our numerical results with established experimental and numerical results wherever available and an excellent match with the experimental results is obtained in all cases.


2013 ◽  
Vol 736 ◽  
pp. 414-443 ◽  
Author(s):  
Y. Ueda ◽  
T. Kida ◽  
M. Iguchi

AbstractThe long-time viscous flow about two identical rotating circular cylinders in a side-by-side arrangement is investigated using an adaptive numerical scheme based on the vortex method. The Stokes solution of the steady flow about the two-cylinder cluster produces a uniform stream in the far field, which is the so-called Jeffery’s paradox. The present work first addresses the validation of the vortex method for a low-Reynolds-number computation. The unsteady flow past an abruptly started purely rotating circular cylinder is therefore computed and compared with an exact solution to the Navier–Stokes equations. The steady state is then found to be obtained for $t\gg 1$ with ${\mathit{Re}}_{\omega } {r}^{2} \ll t$, where the characteristic length and velocity are respectively normalized with the radius ${a}_{1} $ of the circular cylinder and the circumferential velocity ${\Omega }_{1} {a}_{1} $. Then, the influence of the Reynolds number ${\mathit{Re}}_{\omega } = { a}_{1}^{2} {\Omega }_{1} / \nu $ about the two-cylinder cluster is investigated in the range $0. 125\leqslant {\mathit{Re}}_{\omega } \leqslant 40$. The convection influence forms a pair of circulations (called self-induced closed streamlines) ahead of the cylinders to alter the symmetry of the streamline whereas the low-Reynolds-number computation (${\mathit{Re}}_{\omega } = 0. 125$) reaches the steady regime in a proper inner domain. The self-induced closed streamline is formed at far field due to the boundary condition being zero at infinity. When the two-cylinder cluster is immersed in a uniform flow, which is equivalent to Jeffery’s solution, the streamline behaves like excellent Jeffery’s flow at ${\mathit{Re}}_{\omega } = 1. 25$ (although the drag force is almost zero). On the other hand, the influence of the gap spacing between the cylinders is also investigated and it is shown that there are two kinds of flow regimes including Jeffery’s flow. At a proper distance from the cylinders, the self-induced far-field velocity, which is almost equivalent to Jeffery’s solution, is successfully observed in a two-cylinder arrangement.


2016 ◽  
Vol 21 (3) ◽  
pp. 667-681 ◽  
Author(s):  
K.D. Singh

Abstract An unsteady mixed convection flow of a visco-elastic, incompressible and electrically conducting fluid in a hot vertical channel is analyzed. The vertical channel is filled with a porous medium. The temperature of one of the channel plates is considered to be fluctuating span-wise cosinusoidally, i.e., $T^* \left( {y^* ,z^* ,t^* } \right) = T_1 + \left( {T_2} - {T_ 1} \right)\cos \left( {{{\pi z^* } \over d} - \omega ^* t^* } \right)$ . A magnetic field of uniform strength is applied perpendicular to the planes of the plates. The magnetic Reynolds number is assumed very small so that the induced magnetic field is neglected. It is also assumed that the conducting fluid is gray, absorbing/emitting radiation and non-scattering. Governing equations are solved exactly for the velocity and the temperature fields. The effects of various flow parameters on the velocity, temperature and the skin friction and the Nusselt number in terms of their amplitudes and phase angles are discussed with the help of figures.


1994 ◽  
Vol 271 ◽  
pp. 1-16 ◽  
Author(s):  
Peter Y. Huang ◽  
Jimmy Feng ◽  
Daniel D. Joseph

We do a direct two-dimensional finite-elment simulation of the Navier–Stokes equations and compute the forces which turn an ellipse settling in a vertical channel of viscous fluid in a regime in which the ellipse oscillates under the action of vortex shedding. Turning this way and that is induced by large and unequal values of negative pressure at the rear separation points which are here identified with the two points on the back face where the shear stress vanishes. The main restoring mechanism which turns the broadside of the ellipse perpendicular to the fall is the high pressure at the ‘stagnation point’ on the front face, as in potential flow, which is here identified with the one point on the front face where the shear stress vanishes.


Sign in / Sign up

Export Citation Format

Share Document