Dynamics of high-Deborah-number entry flows: a numerical study

2011 ◽  
Vol 677 ◽  
pp. 272-304 ◽  
Author(s):  
A. M. AFONSO ◽  
P. J. OLIVEIRA ◽  
F. T. PINHO ◽  
M. A. ALVES

High-elasticity simulations of flows through a two-dimensional (2D) 4 : 1 abrupt contraction and a 4 : 1 three-dimensional square–square abrupt contraction were performed with a finite-volume method implementing the log-conformation formulation, proposed by Fattal & Kupferman (J. Non-Newtonian Fluid Mech., vol. 123, 2004, p. 281) to alleviate the high-Weissenberg-number problem. For the 2D simulations of Boger fluids, modelled by the Oldroyd-B constitutive equation, local flow unsteadiness appears at a relatively low Deborah number (De) of 2.5. Predictions at higher De were possible only with the log-conformation technique and showed that the periodic unsteadiness grows with De leading to an asymmetric flow with alternate back-shedding of vorticity from pulsating upstream recirculating eddies. This is accompanied by a frequency doubling mechanism deteriorating to a chaotic regime at high De. The log-conformation technique provides solutions of accuracy similar to the thoroughly tested standard finite-volume method under steady flow conditions and the onset of a time-dependent solution occurred approximately at the same Deborah number for both formulations. Nevertheless, for Deborah numbers higher than the critical Deborah number, and for which the standard iterative technique diverges, the log-conformation technique continues to provide stable solutions up to quite (impressively) high Deborah numbers, demonstrating its advantages relative to the standard methodology. For the 3D contraction, calculations were restricted to steady flows of Oldroyd-B and Phan-Thien–Tanner (PTT) fluids and very high De were attained (De ≈ 20 for PTT with ϵ = 0.02 and De ≈ 10000 for PTT with ϵ = 0.25), with prediction of strong vortex enhancement. For the Boger fluid calculations, there was inversion of the secondary flow at high De, as observed experimentally by Sousa et al. (J. Non-Newtonian Fluid Mech., vol. 160, 2009, p. 122).

2021 ◽  
Vol 2119 (1) ◽  
pp. 012152
Author(s):  
D V Platonov ◽  
A V Minakov ◽  
A V Sentyabov

Abstract The paper presents a numerical study of the free discharge of water through the turbine with a braked runner. The simulation was carried out for a unit of a full-scale Francis turbine. The finite volume method was employed for unstructured meshes using the DES method. The simulation results show the flow structures, integral characteristics, and pressure pulsations in the flow path. The analysis of the applicability of this approach to real conditions is carried out.


2005 ◽  
Vol 15 (6) ◽  
pp. 398-405 ◽  
Author(s):  
Shewaferaw S. Shibeshi ◽  
William E. Collins

AbstractBlood flow rheology is a complex phenomenon. Presently there is no universally agreed upon model to represent the viscous property of blood. However, under the general classification of non-Newtonian models that simulate blood behavior to different degrees of accuracy, there are many variants. The power law, Casson and Carreau models are popular non-Newtonian models and affect hemodynamics quantities under many conditions. In this study, the finite volume method is used to investigate hemodynamics predictions of each of the models. To implement the finite volume method, the computational fluid dynamics software Fluent 6.1 is used. In this numerical study the different hemorheological models are found to predict different results of hemodynamics variables which are known to impact the genesis of atherosclerosis and formation of thrombosis. The axial velocity magnitude percentage difference of up to 2 % and radial velocity difference up to 90 % is found at different sections of the T-junction geometry. The size of flow recirculation zones and their associated separation and reattachment point’s locations differ for each model. The wall shear stress also experiences up to 12 % shift in the main tube. A velocity magnitude distribution of the grid cells shows that the Newtonian model is close dynamically to the Casson model while the power law model resembles the Carreau model.


Author(s):  
Zhenhua Wang ◽  
Bengt Sunden ◽  
Shikui Dong ◽  
Zhihong He ◽  
Weihua Yang ◽  
...  

In designing industrial cylindrical furnaces, it is important to predict the radiative heat flux on the wall with high accuracy. In this study, we consider CO2 and H2O which have strong absorption in the infrared range. The absorption coefficients of the gases are calculated by using the statistical narrow band (SNB) model. The spectrum is divided into 15 bands to cover all the absorption regions of the two non-gray gases. The radiative transfer equation is solved by the finite volume method (FVM) in cylindrical coordinates. To make the FVM more accurate, we discretize the solid angle into 80 directions with the S8 approximation which is found to be both efficient and less time consuming. Based on the existing species and temperature fields, which were modeled by the FLUENT commercial code, the radiative heat transfer in a cylinder combustor is simulated by an in-house code. The results show that the radiative heat flux plays a dominant part of the heat flux to the wall. Meanwhile, when the gas is considered as nongray, the computational time is very huge. Therefore, a parallel algorithm is also applied to speed up the computing process.


2014 ◽  
Vol 2014 ◽  
pp. 1-11 ◽  
Author(s):  
Yan Yu ◽  
Pingjian Ming ◽  
Song Zhou

Thermal ground testing is an accepted and frequently used method for simulating the aerodynamic heating of high speed flight vehicles. A numerical method based on a finite volume method for a quartz lamp heating system, used in thermal testing, is proposed. In this study, the unstructured finite-volume method (UFVM) for radiation has been formulated and implemented in a fluid flow solver GTEA on unstructured grids. For comparison and validation of the proposed method, a 2D furnace with cooling pipes was chosen. The results obtained from the proposed FVM agreed well with the exact solutions. Numerical results show that the quartz lamp can be simplified as a slat with the same temperature radiation source, and a simplified 2D thermal testing case was then simulated with the coupling effects of radiation, convection, and conduction heat transfer. Different temperature loading curves and ratios of intervals between the lamps and lamp length (l/s) were studied using the proposed method. The radiation heat flux on the metal surface was a wave-shaped curve. Comparing the different interval ratios, we found that the smaller the interval ratio, the larger the maximum value and the smaller the difference between the maximum and minimum heat flux.


Sign in / Sign up

Export Citation Format

Share Document