Search

This study aims to investigate the wake of two tandem square cylinders based on the Proper Orthogonal Decomposition (POD) analyses of the PIV and hotwire data. The cylinder centre-to-centre spacing ratio L/w examined is from 1.2 to 4.2, covering the four flow regimes, i.e., extended body, reattachment, transition and co-shedding. The Reynolds number examined was 1.3 × 104. A novel Proper Orthogonal Decomposition (POD) technique (hereafter referred to as PODHW) is developed to analyse data from single point hotwire measurements, offering a new perspective compared to the conventional POD analysis (PODPIV) based on Particle Image Velocimetry (PIV) data. A key finding is the identification of two distinct states, reattachment and co-shedding, within the transition flow regime at L/w = 2.8, which PODPIV fails to capture due to the limited duration of the PIV data obtained. This study confirms, for the first time, the existence of these states as proposed by Zhou et al. (2024), highlighting the advantage of using PODHW for capturing intermittent flow phenomena. Furthermore, the analysis reveals how the predominant coherent structures contribute to the total fluctuating velocity energy in each individual regime. Other aspects of the flow are also discussed, including the Strouhal numbers, the contribution to the total fluctuating energy of the flow from the first four POD modes, and a comparison between different regimes.

The objective of this research is to study the fluid flow control allowing the reduction of aerodynamic drag around a square cylinder using two parallel partitions placed downstream of the cylinder using the lattice Boltzmann method with multiple relaxation times (MRT-LBM). In contrast to several existing investigations in the literature that study either the effect of position or the effect of length of a single horizontal or vertical plate, this work presents a numerical study on the effect of Reynolds number (Re), horizontal position (g), vertical position (a), and length (Lp) of the two control partitions. Therefore, this work will be considered as an assembly of several results presented in a single work. Indeed, the Reynolds numbers are selected from 20 to 300, the gap spacing (0 ≤ g ≤ 13), the vertical positions (0 ≤ a ≤ 0.8d), and the lengths of partitions (1d ≤ Lp ≤ 5d). To identify the different changes appearing in the flow and forces, we have conducted in this study a detailed analysis of velocity contours, lift and drag coefficients, and the root-mean-square value of the lift coefficient. The obtained results revealed three different flow regimes as the gap spacing was varied. Namely, the extended body regime for 0 ≤ g ≤ 3.9, the attachment flow regime for 4 ≤ g ≤ 5.5, and the completely developed flow regime for 6 ≤ g ≤ 13. A maximal percentage reduction in drag coefficient equal to 12.5%, is given at the critical gap spacing (gcr = 3.9). Also, at the length of the critical partition (Lpcr = 3d), a Cd reduction percentage of 12.95% was found in comparison with the case without control. Moreover, the position of the optimal partition was found to be equal to 0.8d i.e. one is placed on the top edge of the square cylinder and the second one is placed on the bottom edge. The maximum value of the lift coefficient is reached for a plate length Lp = 2d when the plates are placed at a distance g = 4. On the other hand, this coefficient has almost the same mean value for all spacings between the two plates. Similarly, the root means the square value of the lift coefficient (Clrms) admits zero values for low Reynolds numbers and then increases slightly until it reaches its maximum for Re = 300.

Airflow over a square cross-section at high Reynolds numbers and different angles of incidence is investigated with the aim of providing deeper insight into wind actions on elongated structures and, in particular, tall buildings. The flow around bluff bodies is characterized by separation at sharp corners, as well as possible flow reattachment at side surfaces. The alternate shedding of vortices is also generated in the wake of bluff bodies due to the unsteady nature of flow separation. Two-dimensional (2D) URANS numerical simulations were conducted in order to model transient flow and examine wind actions on a square used as a model of a typical cross-section of a tall building far from its roof and the ground. For validation purposes, the study’s numerical results on drag and lift coefficients, Strouhal numbers, as well as pressure coefficient distribution were found to be in good agreement with available experimental and numerical results in the literature for relatively low Reynolds numbers. The numerical study was then extended to higher Reynolds numbers, approaching values that are pertinent for wind flow around buildings, thus addressing the lack of such results in the literature. On the basis of these results, the impact of Reynolds numbers and angles of incidence on drag and lift coefficients, as well as the pressure coefficient distribution along the walls of the cross-section, is highlighted.

In this work, we present a grid study oriented to capture 3D flow simulations around smooth and wrinkled cylinders that could have practical applications in various engineering areas. The study considers three Reynolds numbers, namely, a benchmark Re =2.14×104 and two orders of magnitude above and below it. The main contributions of the paper relate to the optimization of the computational mesh for the spanwise direction of the wind flow that results from the computational-mathematical framework employed, in addition to a novel visualization technique that unfolds features in the recording data that could otherwise be hidden when using traditional plots. We compare our benchmark results with those reported by other authors to conclude that the intermediate resolution grids employed with the widest spanwise provide acceptable results. Furthermore, the new visualization technique offers significant advantages compared to traditional pressure maps, regarding clarity for observing and interpreting local flow disturbances, making variations with Re clearer, and by enabling the detection of asymmetries.

The current manifestation is utilized to explicate the inspiration of combined convection flows in a cavity with a square cylinder placed at the center of the cavity having coordinates (0.5, 0.5). To be more specific, right- and left-sided vertical walls are kept cold, and the lower wall is to be heated in uniform and non-uniform manners, while the upper horizontal wall is moved with a constant velocity ULid and is thermally adiabatic. The obstacle is treated as cold as well as thermally adiabatic, and the no-slip velocity boundary conditions are specified at its surface. For the purpose of computing the velocity profile and temperature, a space including the quadratic polynomials (P2) is chosen; however, the pressure has been approximated using a linear ( P1 ) finite element space of functions. The Newton technique is used to perform the computations needed to solve the discrete systems of non-linear algebraic equations. The non-linear iterations are terminated at residual below 10–6, whereas the inner core linear solver is based on Gaussian elimination with special reordering of unknowns. To show the consistency of the implemented numerical technique, the parametric study is designed based on the most relevant non-dimensional parameters, namely, the Grashof number Gr ranging from 102 to 105, the Hartmann number Ha changing from 0 to 50, and the Reynolds number Re varying from 10 to 200. Computations in the forms of velocity streamlines and isotherm contour profiles are adorned. In addition, the production of temperature disparities is associated with an increase in Nuavg as Re increases. In contrast, a diminishing aptitude in kinetic energy is observed due to the creation of Lorentz forces.

The transonic flow around a square cylinder at Ma= 0.71 and Re= 4×105 has been studied by using the scale-adaptive simulation (SAS) method, and the characteristics of separated shear layer and wake have been analyzed in depth. To validate the SAS approach, the SAS results are compared with the existing numerical and experimental results. In the present transonic flow, the convective Mach number inside the shear layer is about 0.6. This indicates that the initial evolution of the separated shear layer is dominated by Kelvin-Helmholtz instability, and the roller spanwise eddies in the initial stage of the shear layer can be observed. In the regions near the shear layer and the wake, the doubling frequencies can be obtained indicative of the harmonic phenomenon inside the separated shear layer, which is closely related to the obvious merging of the vortices in the shear layer. Proper orthogonal decomposition of the pressure field shows that the transonic flow field of square cylinder is dominated by the antisymmetric mode, which is associated with the vortex shedding in the wake and the propagation of compression waves induced by the shear layer.

This work numerically investigates the effects of combined rotational and transverse oscillations of a square cylinder on the flow field and force coefficients. The primary non-dimensional parameters that were varied are frequency ratio fR (0.5, 0.8), Re (50-200), phase difference (ϕ) between the motions and rotational amplitude (θ0) with the influence of the last three parameters being discussed in detail. The amplitude of transverse oscillations is fixed at 0.2D, where D is the cylinder width. The study has been validated using the mean drag coefficient for stationary and transversely oscillating square cylinders from literature. Output data was obtained in the form of force coefficient (Cd), vorticity and pressure contours. The governing equations for the 2dimensional model were solved from an inertial frame of reference (overset meshing) using finite volume method. The interplay between the convective field and prescribed motion has been used to explain many of the results obtained. The relative dominance of cylinder motion over the flow stream was determined using Discrete Fast Fourier Transform. The influence of Re on Cd disappears when the motions are completely out of phase (ϕ = π). In general, the Cd for low Re flows exhibited low sensitivity to change in other parameters. Direct correlation has been observed between frontal area, vortex patterns and drag coefficient

Square cylinders are widely used in various fields. For example, they are common structures in fishways. The flow around square cylinders has been a common problem in various fields. However, reducing the flow drag of the square cylinder is a problem that remains unexplored. Many previous studies have reported the drag reduction of 2D square cylinders, which failed to reflect the drag of real structures. Also, some studies focus on the drag force of the inner wall of the square cylinder modified by the microstructure. Achieving drag reduction by microstructuring the surface of the 3D square cylinder is a challenging problem. This study applied a 3D numerical simulation and deep neural network to study the drag reduction performance of the square cylinder under different patch sizes. We studied the drag reduction performance of protrusion and pit-patched square cylinders and tried to find the rule between drag reduction performance and patch configuration. The results show that the square cylinder has better drag reduction performance in some cases. However, its drag reduction performance is greatly affected by the protrusion structure. Also, too large protrusions will increase the drag force of the structure. When the surface protrusion accounts for 10% of the total area of the square cylinder, the drag reduction performance is the best (22.1%). The pit patch structure demonstrated an insignificant drag reduction performance and even increased the drag in most cases. The DNN prediction model demonstrated the robustness of the numerical simulation data.

Under the action of incoming flow, the square cylinder can generate more intense vibration responses than the circular cylinder, which is beneficial for energy harvesting. Numerical simulations for FIV of the square-cylinder energy conversion system are carried out. URANS equations are used in conjunction with the shear stress transport k-ω turbulence model to predict the flow, and the equations for vibrations are solved by the Newmark-β algorithm. The present numerical method is validated against the published data with good consistency. The Reduced velocity Ur is varied from 1-20, with corresponding Reynolds numbers of 24 000-160 000. The numerical results indicate that the Reynolds number significantly affects the frequency response, amplitude response, vortex shedding mode, and energy conversion efficiency. The highest efficiency point locates at Re=88 000, with a value of 7.156%. When Re>120 000, the system transits from vortex-induced vibration into galloping, and its vibration responses as well as energy harvesting characteristics change sharply. Fully developed galloping motion occurs when Re>144 000.

Sınırlandırılmış bir kanal içerisinde yer alan tek kare silindir (KS) üzerinden iki boyutlu laminer sürekli akış için blokaj (=B/H) oranının ısı transferi ve akış karakteristiklerine olan etkisi incelenmiştir. Çalışmada Reynolds sayısı Re=40 değerinde sabit tutulurken blokaj oranı =0.125-0.8 değerleri arasında değiştirilmiştir. Hesaplamalarda ANSYS CFX 14.0 kullanılmıştır. Blokaj oranı etkisinin KS yüzeyleri üzerindeki sürükleme katsayısı (Cd), sürtünme katsayısı (Cf), boyutsuz yeniden birleşme uzunluğu (Lr/B) ve ortalama Nusselt sayısı (Nu) üzerine olan etkileri incelenmiştir. Blokaj oranı arttıkça sürükleme katsayısı (Cd), sürtünme katsayısı (Cf) ve ortalama Nusselt sayısı (Nu) değerlerinin arttığı ancak boyutsuz yeniden birleşme uzunluğu (Lr/B) değerinin azaldığı bulunmuştur. Sürükleme katsayısı (Cd), sürtünme faktörü (Cf), boyutsuz yeniden birleşme uzunluğu (Lr/B) ve ortalama Nusselt sayısı değerlerinin blokaj oranına göre değişimini veren bağıntılar elde edilmiştir.

An experimental study was carried out to identify the effect of angle of attack on flow-induced vibration (FIV) of square cylinders. The experiment was conducted at the Aeronautical and Wind Engineering Laboratory (AEROLAB), UTM Kuala Lumpur using a wind tunnel that was free from external wind conditions. A supporting structure was designed and fabricated to conduct this experiment. The importance of this support structure was to enable the rigid cylinder to suspend and vibrate freely upon excitation of wind speed. The results were analysed through the response of amplitude and frequency of the rigid cylinder over a velocity range of 0.5m/s to 4.0m/s. The results showed that for a square cylinder of ?=0°, vortex-induced vibration (VIV) occurred at low reduced velocity (UR) in range of 5 ? UR ? 10 and galloping occurred at higher reduced velocity which started at UR=15. A tranquil zone was found between VIV and galloping in the reduced velocity range of 10 ? UR ? 15. As for ?=22.5° and 45°, only VIV response was found at low reduced velocity in range of 4? UR ? 9. ABSTRAK: Satu kajian eksperimentasi telah dilakukan bagi mengenal pasti pengaruh sudut serangan oleh getaran cetusan-aliran (FIV) dalam silinder persegi. Eksperimen ini dijalankan di Makmal Kejuruteraan Aeronautika dan Angin (AEROLAB), UTM Kuala Lumpur dengan menggunakan terowong angin yang bebas dari pengaruh angin luar. Struktur sokongan telah direka dan difabrikasi bagi tujuan eksperimen ini. Ini penting bagi membolehkan silinder pegun tergantung dan bergetar dengan bebas semasa ujian kelajuan angin. Dapatan kajian dianalisis melalui tindak balas amplitud dan frekuensi silinder pegun pada kadar halaju 0.5m/s sehingga 4.0m/s. Hasil kajian menunjukkan bahawa bagi silinder persegi ? = 0 °, getaran pengaruh-vorteks (VIV) berlaku pada halaju rendah (UR) dalam julat 5 ? UR ? 10 dan getaran lebih teruk telah ketara berlaku pada kadar halaju berkurang iaitu bermula pada UR = 15. Zon tenang dijumpai antara VIV dan getaran teruk pada kadar halaju berkurang 10 ? UR ? 15. Adapun pada ? = 22.5° dan 45°, hanya tindak balas VIV dijumpai pada halaju rendah dalam kadar 4? UR ? 9.

New data and review of the spanwise coherence length is provided for flows over cylinders of different cross-sections: circular of diameter d, and rectangular of sectional aspect ratios (breadth (b) to height (d) ratio AR = b/d) of 1, 2 and 3. In the present measurements, the body has both d and spanwise length of 70d fixed, and the Reynolds number (based on d) range 6000–27,000 is covered. Two-point data are obtained from two hot-wire probes, one fixed in the symmetry plane and the other moving on the corresponding spanwise axis. Their position in a cross plane are deduced from preliminary measurement of the mean flow with a single probe, allowing fair comparisons between the different geometries and the introduction of uncertainty bars on coherence length values. At all tested regimes, a very good agreement is noticed between velocity-based and pressure-based coherence experimental data. Coherence length definitions are revisited, and the aeroacoustically consistent, integral length definition is selected, allowing fair synthesis of literature data into a single chart and empirical functions. Definitions for coherence decay models (e.g. Gaussian or Laplacian) are also adapted so that coherence length and coherence integral shall be equivalent. This preliminary work on coherence data and its spanwise integration enables transparent regressions and model selection. Generally, the Gaussian model is relevant for the lift peak, while the coherence exhibits a Laplacian decay at harmonics. On average, at peak Strouhal number, the coherence length for the circular and square cylinders is of 5d while it is of the order of 15d for the rectangular sections. It is concluded that the flow over those latter geometries is still a two-dimensional dynamics at the tone frequency. These values are almost preserved over the tested Reynolds number range. Coherence length value at harmonics is extensively documented. Spanwise coherence length is also discussed as an ingredient of acoustic efficiency.

Corner modifications can reduce wind loads acting on supertall buildings and modify the corresponding flow structures. The present study investigated the aerodynamic shape optimization of the corner recession square cylinders with multiple geometric parameters in a large design space via the GA-GRNN surrogate model updating-based multi-objective optimization framework. Six typical optimal aerodynamic shape sections M1~M6 were selected from the Pareto optimal front, and the effects of multiple geometric parameters of these sections on the aerodynamic performance and flow field were analyzed. The results showed that the present multi-objective optimization framework can significantly reduce the computational load and time cost, and significantly improve the optimization efficiency in solving complex engineering problems. The optimal corner recession sections can obviously reduce the mean drag coefficient CD and root mean square lift coefficient CσL while significantly increasing the Strouhal number St of the square cylinder, and it is concluded that the aerodynamic shape optimization can significantly improve the aerodynamic performance of square-sectional supertall buildings. When compared with the benchmark section, the CD and CσL of the optimal section M1 can be reduced up to 45.7% and 84.5%, respectively. Based on the analysis of the flow structures around the optimal sections, the flow mechanism can be attributed to the fact that the corner recession modifications postpone the flow separation, and deflect the separated shear layer towards the side surfaces and suppress the development of vortex shedding in the wake, which leads to significant elongation of the wake length and reduction of the width of the recirculation region. The proposed multi-objective optimization framework in this study can provide an important reference for the aerodynamic shape optimization of building structures and relevant studies.

Abstract Unsteady RANS (URANS), hybrid LES/RANS and IDDES simulations were conducted to numerically investigate the velocity field around, and pressures distribution and forces over a square cylinder immersed in a uniform, steady oncoming flow with Reynolds number Re = 21,400. The vortex shedding responses in terms of Strouhal number, the pressure distribution, the velocity profile and the velocity fluctuations obtained by numerical simulations are compared with experimental data. Compared with 2D URANS simulation, 3D simulations using hybrid LES/RANS and IDDES models provide more accurate prediction on the responses in the wake, including mean streamwise velocity profile and rms velocity fluctuations. This also results in more accurate prediction of time-averaged surface pressure coefficient on the rear surface obtained by 3D hybrid LES/RANS and IDDES simulations than by URANS simulation. When a hybrid LES/RANS model or IDDES model is used, a more accurate prediction for either pressure coefficient or velocity profile (especially in the far wake region) is not guaranteed by increasing the mesh resolution along the spanwise direction of the square cylinder.

The flow past a square cylinder under the influence of a one dimensional gust was investigated using computational fluid dynamics (CFD). The effect of upstream wind gusts of the same amplitude but different duration was investigated with respect to their effect on the flow, the vortex-shedding, and the pressure distribution around the square cylinder. For the computations, a very large eddy simulation (VLES) model was implemented in an in-house code and validated against numerical and experimental results from the literature. The gusts of different duration were found to have a distinctly different effect. The short-duration gust causes a lock-on behavior with cessation of the alternating vortex shedding, and a symmetric pair-vortex was created above and below the square cylinder. It was observed that the pressure distribution on the lateral sides of the cylinder has the same magnitude and phase, which resulted in a zero total lift coefficient. In terms of a free-standing structures, such as a building, this would lead to zero instantaneous forces and pressure difference in the lateral direction with obvious implications for dynamic response and cross ventilation.

A preliminary study is carried out to investigate the aerodynamic characteristics of a square cylinder with Savonius wind turbines and to explain the reason why this kind of structure can suppress wind-induced vibrations. A series of computational fluid dynamics simulations are performed for the square cylinders with stationary and rotating wind turbines at the cylinder corners. The turbine orientation and the turbine rotation speed are two key factors that affect aerodynamic characteristics of the cylinder for the stationary and rotating turbine cases, respectively. The numerical simulation results show that the presence of either the stationary or rotating wind turbines has a significant effect on wind forces acting on the square cylinder. For the stationary wind turbine cases, the mean drag and fluctuating lift coefficients decrease by 37.7% and 90.7%, respectively, when the turbine orientation angle is 45°. For the rotating wind turbine cases, the mean drag and fluctuating lift coefficients decrease by 34.2% and 86.0%, respectively, when the rotation speed is 0.2 times of vortex shedding frequency. Wind turbines installed at the corners of the square cylinder not only enhance structural safety but also exploit wind energy simultaneously.

In this paper, the entire process of the flow around a fixed square cylinder and the moving surface boundary-layer control (MSBC) at a low Reynolds number was numerically simulated. Two small rotating circular cylinders were located in each of the two rear corners of the square cylinder, respectively, to transfer momentum into the near wake behind the square cylinder. The rotations of the two circular cylinders were realized via dynamic mesh technology, when the two-dimensional incompressible Navier–Stokes equations for the flow around the square cylinder were solved. We analyzed the effects of different rotation directions, wind angles θ, and velocity ratios k (the ratio of the tangential velocity of the rotating cylinder to the incoming flow velocity) on the wake of flow around a square cylinder to evaluate the control effectiveness of the MSBC method. In the present work, the aerodynamic forces, the pressure distributions, and the wake patterns of the square cylinder are discussed in detail. The results show that the high suction areas near the surfaces of the rotating cylinders can delay or prevent the separation of the shear layer, reduce the wake width, achieve drag reduction, and eliminate the alternating vortex shedding. For a wind angle of 0°, the inward rotation of the small circular cylinders is the optimal arrangement to manipulate the wake vortex street behind the square cylinder, and k=2 is the optimal velocity ratio between the control effectiveness and external energy consumption.

Turbulent uniform flows past a two-dimensional square cylinder are investigated numerically. By varying the turbulence intensity and turbulence length scale of the approaching flow, the flow effect of the cylinder are compared to that in a laminar approaching-flow case. In addition, the variations of drag and lift coefficients with respect to the changes of turbulence intensity and turbulence length scale are analyzed on a systematic basis. In the large eddy simulations, the approaching-flow turbulence is generated by a spectral method according to Kármán spectrum. Two levels of turbulence intensities (5% and 10%) and three turbulence length scales (0.25, 0.50 and 1.0 times of the cylinder width) are selected in the study to examine the effect on the cylinder. Results show that the Strouhal number remains almost unchanged when the uniform approaching-flow changes from a laminar state to a turbulent one. The approaching-flow turbulence has noticeable effect in promoting the resulting drag and lift fluctuations. However, its effect on the mean drag appears negligible. In contrast, an increase of the approaching-flow turbulence length scale leads to mild increases of the mean and root-mean-square values of drag. On the other hand, the resulting lift fluctuation is insensitive to the change of the turbulence length scale.

The aerodynamic behavior of a square cylinder with rounded corner edges in steady flow regime in the range of Reynolds number (Re) 5–45; is predicted by Artificial Neural Network (ANN) using MATLAB. The ANN has trained by back propagation algorithm. The ANN requires input and output data to train the network, which is obtained from the commercial Computational Fluid Dynamics (CFD) software FLUENT in the present study. In FLUENT, all the governing equations are discretized by the finite volume method. Results from numerical simulation and back propagation based ANN have been compared. It has been discovered that the ANN predicts the aerodynamic behavior correctly within the given range of the training data. It is additionally observed that back propagation based ANN is an effective tool to forecast the aerodynamic behavior than simulation, that has very much longer computational time.

An unsteady two-dimensional forced convection over a square cylinder with sharp and rounded corner edge is numerically analyzed for the low Reynolds number laminar flow regime. In this study, the analysis is carried out for Reynolds number (Re) in the range of 80 to 180 with Prandtl number (Pr) variation from 0.01 to 1000 for various corner radius (r=0.50, 0.51, 0.54, 0.59, 0.64 and 0.71). The lateral sides of the computational domain are kept constant to maintain the blockage as 5%. Heat transfer due to unsteady forced convection has been predicted by Artificial Neural network (ANN). The present ANN is trained by the input and output data which has been acquired from the numerical simulation, performed in finite volume based Computational Fluid Dynamics (CFD) commercial software FLUENT. The heat transfer characteristics over the sharp and rounded corner square cylinder are evaluated by analyzing the local Nusselt number (Nulocal), average Nusselt number (Nuavg) at various Reynolds number, Prandtl numbers and for various corner radii. It is found that the heat transfer rate of a circular cylinder can be enhanced by 12% when Re is varying and 14% when Prandtl number is varying by introducing a new cylinder geometry of corner radius r=0.51. It is found that the unsteady forced convection heat transfer over a cylinder can be predicted appropriately by ANN. It is also observed that the back propagation ANN can predict the heat transfer characteristics of forced convection very quickly compared to a standard CFD method.

A quasi two-dimensional numerical study is performed to analyze the thermo-magneto-convective transport of liquid metal around a square cylinder in a square duct subjected to a strong externally imposed axial magnetic field. The channel bottom wall is considered heated while the top wall is maintained at the free stream temperature keeping the cylinder adiabatic. The Reynolds and Hartmann numbers are kept in the range and . The flow dynamics in the aforementioned range of parameters reveals the existence of four different regimes out of which the first three ones are similar to the classical non-MHD 2-D cylinder wakes while the fourth one is characterized by the vortices evolved from the duct side walls due to the boundary layer separation which strongly disturbs the Kármán vortex street. The flow dynamics and heat transfer rate from the heated channel wall are observed to depend on the imposed magnetic field strength. With increasing magnetic field, the flow becomes stabilized resulting in a degradation in the forced convection heat transfer. A special case at a very high Reynolds number with Ha = 2160 is also considered to show the development of a Kelvin–Helmholtz-type instability that substantially affects the heat transfer rate.

Flow around an extended triangular solid (thorn) attached to a square cylinder is investigated numerically. The numerical analysis is carried out at low Reynolds number, Re = 100 & 180 for different non-dimensional thorn lengths (l΄ = 0. 2, 0.4 & 0.6), different inclination angles (θ = 5°, 10°, 15° and 20°) and two different thorn positions. It is found that drag and lift reduction can be achieved by attaching the thorn on a square cylinder. It is observed that the fluctuation of the drag force as well as the lift force is reduced and there is a comparatively large variation of drag and lift when the thorn is placed at the front stagnation point instead of placing at rear stagnation point. The reduction of drag and lift coefficient are directly proportional to thorn length and thorn inclination angle. It is found that the drag and lift are minimized by 16% & 46% for Re = 100 respectively, and 22% & 60% for Re = 180 compared to a square model (without thorn).

Numerical analysis of two dimensional steady flows over a triangular extended solid (thorn) attached to square cylinder positioned at front stagnation point and at rear stagnation point separately is reported at low Reynolds number (Re = 40). The variation of thorn length (l′ = 0.2, 0.4 & 0.6) & inclination angle (Θ = 5°, 10°, 15° & 20°) and its effect on the drag, pressure, shear stress, boundary layer as well as on the inertia force and viscous force are of interest. There is comparatively large variation on drag when the thorn is placed at the front side instead of placing at rear. The recirculation length is remained constant by varying the length and inclination of the thorn irrespective of its position. The variation of drag is comparatively less by changing thorn inclination. It is found that the drag is minimized by 2–3% compared to square model.

In this research, an experimental study was conducted to high light the impact of the exterior shape of a cylindrical body on the forced and free convection heat transfer coefficients when the body is hold in the entrance of an air duct. The impact of changing the body location within the air duct and the air speed are also demonstrated. The cylinders were manufactured with circular, triangular and square sections of copper for its high thermal conductivity with appropriate dimensions, while maintaining the surface area of all shapes to be the same. Each cylinder was heated to a certain temperature and put inside the duct at certain locations. The temperature of the cylinder was then monitored. The heat transfer coefficient were then calculated for forced convection for several Reynolds number (4555-18222).The study covered free convection impact for values of Rayleigh number ranging between (1069-3321). Imperical relationships were obtained for all cases of forced and free convection and compared with equations of circular cylindrical shapes found in literature. These imperical equations were found to be in good comparison with that of other sources.

A detailed numerical investigation of the flow behind a square cylinder at a Reynolds number of 21,400 is conducted to assess the ability of the delayed detached-eddy simulation (DDES) modeling approach to accurately predict the velocity recovery in the wake of a bluff body. Three-dimensional unsteady Reynolds-averaged Navier–Stokes (URANS) and DDES simulations making use of the Spalart–Allmaras turbulence model are carried out using the open-source computational fluid dynamics (CFD) toolbox OpenFOAM-2.1.x, and are compared with available experimental velocity measurements. It is found that the DDES simulation tends to overestimate the averaged streamwise velocity component, especially in the near wake, but a better agreement with the experimental data is observed further downstream of the body. The velocity fluctuations also match reasonably well with the experimental data. Moreover, it is found that the spanwise domain length has a significant impact on the flow, especially regarding the fluctuations of the drag coefficient. Nonetheless, for both the averaged and fluctuating velocity components, the DDES approach is shown to be superior to the URANS approach. Therefore, for engineering purposes, it is found that the DDES approach is a suitable choice to simulate and characterize the velocity recovery in a wake.

The flow of a non-Newtonian, power-law fluid, directed normally to a horizontal cylinder with square cross-section (two-dimensional flow) is considered in the present paper. The problem is investigated numerically with a very large calculation domain in order that the flow could be considered unconfined. The investigation covers the power-law index from 0.1 up to 2 and the Reynolds number ranges from 60 to 160. Over this range of Reynolds numbers the flow is unsteady. It is found that the drag coefficient and the Strouhal number are higher in a confined flow compared to those of an unconfined flow. In addition some flow characteristics are lost in a confined flow. Complete results for the drag coefficient and Strouhal number in the entire shear-thinning and shear-thickening region have been produced. In shear-thinning fluids chaotic structures exist which diminish at higher values of power-law index. This study represents the first investigation of unsteady, non-Newtonian power-law flow past a square cylinder in an unconfined field.

In this research, an experimental study was conducted to high light the impact of the exterior shape of a cylindrical body on the forced and free convection heat transfer coefficients when the body is hold in the entrance of an air duct. The impact of changing the body location within the air duct and the air speed are also demonstrated. The cylinders were manufactured with circular, triangular and square sections of copper for its high thermal conductivity with appropriate dimensions, while maintaining the surface area of all shapes to be the same. Each cylinder was heated to a certain temperature and put inside the duct at certain locations. The temperature of the cylinder was then monitored. The heat transfer coefficient were then calculated for forced convection for several Reynolds number (4555-18222).The study covered free convection impact for values ​​of Rayleigh number ranging between (1069-3321). Imperical relationships were obtained for all cases of forced and free convection and compared with equations of circular cylindrical shapes found in literature. These imperical equations were found to be in good comparison with that of other sources.

The paper is devoted to the problem of boundary conditions influence on the quality of the solution obtained with use of k-ε turbulence models. There are calculation results for different boundary conditions and two methods: standard k-ε and RNG k-ε in the paper. The flow parameters obtained from the calculation are compared with our own measurement results. Moreover, the influence of input data on the inflow edge on sensitivity coefficients is shown and analysed in the paper. The research is performed for components of velocity and turbulence kinetic energy.