Case Of Steady Flow Research Articles

Chaotic Behavior and Heat Transfer Analysis of Pulsating Flows with Different Frequencies

Due to the rapid development of electronic devices, the improvement of air–cooling methods for thermal management is urgent. In the present work, a promising way of coupling pulsating flow and vortex generators is proposed to eliminate the impacts of severe local high-temperature regions behind the vortex generators. Different from the traditional qualitative analysis, chaos is introduced to identify the quantitative impacts of the nonlinear dynamics of the coupling system. Particularly, the flow field was examined using the line stretching and the Poincaré map, and the chaotic characteristics of high-frequency pulsating flow were investigated. The results show that high-frequency pulsating flow enhances heat transfer mainly by promoting chaotic mixing. The chaotic fluid repeatedly stretches and folds, which increases the line stretching. Thus, the line stretching distribution can represent the heat transfer. Compared to the steady flow case, at Reynolds number of 1024 and pulsation frequency of 250 Hz, pulsating flow improves the Nusselt number by 5.82% and reduces friction factor by 2.62%, resulting in an improvement of thermal performance factor by 6.76%.

A one-dimensional computational model for blood flow in an elastic blood vessel with a rigid catheter.

Strokes are one of the leading causes of death in the United States. Stroke treatment involves removal or dissolution of the obstruction (usually a clot) in the blocked artery by catheter insertion. A computer simulation to systematically plan such patient-specific treatments needs a network of about 105 blood vessels including collaterals. The existing computational fluid dynamic (CFD) solvers are not employed for stroke treatment planning as they are incapable of providing solutions for such big arterial trees in a reasonable amount of time. This work presents a novel one-dimensional mathematical formulation for blood flow modeling in an elastic blood vessel with a centrally placed rigid catheter. The governing equations are first-order hyperbolic partial differential equations, and the hypergeometric function needs to be computed to obtain the characteristic system of these hyperbolic equations. We employed the Discontinuous Galerkin method to solve the hyperbolic system and validated the implementation by comparing it against a well-established 3D CFD solver using idealized vessels and a realistic truncated arterial network. The results showed clinically insignificant differences in steady flow cases, with overall variations between 1D and 3D models remaining below 10%. Additionally, the solver accurately captured wave reflection phenomena at domain discontinuities in unsteady cases. A primary advantage of this model over 3D solvers is its ease in obtaining a discretized geometry of complex vasculatures with multiple arterial branches. Thus, the 1D computational model offers good accuracy and applicability in simulating complex vasculatures, demonstrating promising potential for investigating patient-specific endovascular interventions in strokes.

Heat transfer characteristics of pulsating flow in a straight channel with rhombic-shaped expanding chamber: A numerical study

ABSTRACT The combination of passive and active heat transfer recovery techniques in channel flows has significant potential in terms of increasing thermal performance. Therefore, this study focused on the numerical examination of the flow and thermal behavior of pulsating flow in a straight duct containing a rhombus chamber. In the context of the study, iterations were solved using the finite volume method (FVM). The study was carried out for different pulsating amplitudes (A: 0.2, 0.4, 0.6, and 0.8), Strouhal numbers (St: 1, 2, 3, 4), and Reynolds numbers (200 ≤ Re ≤ 800). The surfaces other than the adiabatic lengths at the inlet and outlet of the duct were kept constant at Tw = 360 K. Results were compared with the results of the steady flow case. The effects of pulsating velocity and oscillatory parameters on Nusselt number, pressure drop, and performance factor were discussed. Velocity and temperature images were presented for different Reynolds numbers and pulsating components in the channel. The findings revealed that although the pulsating parameters significantly enhanced the heat transfer at increasing Reynolds numbers, they had a fairly low effect on heat transfer at the same Reynolds number. It was observed that at Re = 800, the Nusselt number formed a peak at St = 3 for all tested pulsating amplitudes. For the parameters of Re = 800, A = 0.6, and St = 2, heat transfer and performance factor in pulsating flow increased by 8.96 and 8.22 times, respectively, compared to steady flow conditions.

Nonlinear ice sheet/liquid interaction in a channel with an obstruction

The interaction between the flow in a channel with an obstruction on the bottom and an elastic sheet representing the ice covering the liquid is considered for the case of steady flow. The mathematical model based on the velocity potential theory and the theory of thin elastic shells fully accounts for the nonlinear boundary conditions at the elastic sheet/liquid interface and on the bottom of the channel. The integral hodograph method is employed to derive the complex velocity potential of the flow, which contains the velocity magnitude at the interface in explicit form. This allows one to formulate the coupled ice/liquid interaction problem and reduce it to a system of nonlinear equations in the unknown magnitude of the velocity at the interface. Case studies are carried out for a semi-circular obstruction on the bottom of the channel. Three flow regimes are studied: a subcritical regime, for which the interface deflection decays upstream and downstream; an ice supercritical and channel subcritical regime, for which two waves of different lengths may exist; and a channel supercritical regime, for which the elastic wave is found to extend downstream to infinity. All these regimes are in full agreement with the dispersion equation. The obtained results demonstrate a strongly nonlinear interaction between the elastic and the gravity wave near the first critical Froude number where their lengths approach each other. The interface shape, the bending moment and the pressure along the interface are presented for wide ranges of the Froude number and the obstruction height.

Effects of Flow Pulsation and Surface Geometry on Heat Transfer Performance in a Channel With Teardrop-Shaped Dimples Investigated by Large Eddy Simulation

Abstract Large eddy simulation was performed to investigate heat transfer performance of a pulsating flow over teardrop-shaped dimples. A total of six geometries of dimpled surfaces were examined for dimple arrangements of in-line/staggered/original and dimple inclination angle of 0–60 deg. Pulsating flows were generated by sinusoidally varying the volume-averaged velocity. The pulsation frequency and amplitude were changed for the Strouhal number of 0–0.60 and the root-mean-square velocity amplitude normalized by the bulk flow velocity of 0–0.14. The results showed that the surface-averaged Nusselt number and friction factor were larger for the pulsating flow case than for the steady flow case. The surface-averaged Nusselt number ratio and the friction factor increased with the Strouhal number up to the Strouhal number of 0.30. For the Strouhal number larger than 0.30, they decreased with the Strouhal number or stayed almost constant. Consequently, the heat transfer efficiency index increased with the Strouhal number. The increase in the local Nusselt number ratio due to the flow pulsation was observed at the leading-edge region of the dimples. The results of the streamlines near the dimple showed that the swirling separation bubble was located closer to the leading-edge region due to the pulsation, which resulted in the increase of the absolute values of the turbulent heat flux and the local Nusselt number ratio.

Acoustic streaming: insights across Reynolds numbers

When a fluid system is subjected to an acoustic wave (or another periodic actuation), the response of the fluid is not purely periodic, but is rather characterized by the combination of a periodic flow and a steady Stokes drift component, where the former is, in many cases, an acoustic wave and the latter is commonly referred to as acoustic streaming. Classical theories of acoustic streaming have focused on slow acoustic streaming, where the periodic flow is the leading-order flow, and is insensitive to the steady flow component which appears as a small correction and is characterized by a small hydrodynamic Reynolds number. In contrast, Dubrovski et al. (J. Fluid Mech. vol. 975, 2023, A4) tackle the fast acoustic streaming regime – conceived by Zarembo (Acoustic streaming. In High-Intensity Ultrasonic Fields, 1971, pp. 135–199. Springer) approximately fifty years ago – where both the periodic and steady flow components are of a similar order of magnitude such that the periodic flow both supports and is simultaneously impacted by the steady flow. They present a novel theoretical framework that accounts for the convection of momentum both within and between the periodic and steady flow to extend slow-streaming equations to the case of steady flow with arbitrary hydrodynamic Reynolds number. They leverage a scaling analysis of the resulting system of equations and a case study to demonstrate the compatibility of their equations with slow streaming theories and highlight the distinctive features of fast streaming.

Iterative dynamics-based mesh discretisation for multi-scale coastal ocean modelling

Flow in coastal waters contains multi-scale flow features that are generated by flow separation, shear-layer instabilities, bottom roughness and topographic form. Depending on the target application, the mesh design used for coastal ocean modelling needs to adequately resolve flow features pertinent to the study objectives. We investigate an iterative mesh design strategy, inspired by hydrokinetic resource assessment, that uses modelled dynamics to refine the mesh across key flow features, and a target number of elements to constrain mesh density. The method is solver-agnostic. Any quantity derived from the model output can be used to set the mesh density constraint. To illustrate and assess the method, we consider the cases of steady and transient flow past the same idealised headland, providing dynamic responses that are pertinent to multi-scale ocean modelling. This study demonstrates the capability of an iterative approach to define a mesh density that concentrates mesh resolution across areas of interest dependent on model forcing, leading to improved predictive skill. Multiple design quantities can be combined to construct the mesh density, refinement can be applied to multiple regions across the model domain, and convergence can be managed through the number of degrees of freedom set by the target number of mesh elements. To apply the method optimally, an understanding of the processes being model is required when selecting and combining the design quantities. We discuss opportunities and challenges for robustly establishing model resolution in multi-scale coastal ocean models.

Analytical solution for unsteady Walters-B fluid flow by a deforming surface with acceleration using OHAM based package BVPh2.0

Current study aims at simulating fluid flow due to a deformable heated surface in an otherwise static viscoelastic fluid obeying Walters-B model. Velocity of the surface is supposed to grow as time from its initiation of motion progress. Simulations in this work are based on the assumption of quadratic surface temperature distribution. Temperature rise attributed to the frictional heating effect is accounted for in the analysis. By choosing appropriate base functions, homotopy solutions are developed for reasonably large values of material fluid parameter. Reliability of the analytical results is established by computing averaged squared residual of the system. The contributions of the surface acceleration and elasticity on the boundary layer formation are enlightened through the plots of velocity components and temperature. Skin friction measuring the stress experienced by the surface is evaluated and examined under different controlling parameters. The paper also presents a numerical solution using NDSolve of MATHEMATICA in a special case of steady flow, and such solution agrees very well with the corresponding homotopy solution.

Effect of unsteady supersonic flow on detonation under different hot jet initiation conditions

To understand the actual operation of detonation engines, it is of great significance to simulate the detonation characteristics after being affected by unsteady supersonic flow. Adaptive mesh refinement method and single-step chemical reaction model were used for simulation calculations. The Navier–Stokes equations and hybrid sixth-order weighted essentially non-oscillatory centered difference scheme was adopted in this study. Three different amplitudes of disturbances (0.1, 0.2, and 0.3) and two hot jet injection times were analyzed. The results showed that a moderate disturbance can improve the stability of detonation wave under the ideal hot jet injection time. It allows detonation wave to stay in the combustion chamber for a longer time, which is beneficial for the release of thrust performance. Although the stability of detonation wave is significantly different between the low and high-amplitude disturbances, the residence time of detonation wave is close to the steady flow case under ideal initiation conditions. When the ignition time is insufficient, successful detonation can occur with low and medium-amplitude disturbances. However, the ignition delay was significantly prolonged with medium-amplitude disturbance. There was also greater instability, which is detrimental to the thrust release. When the inflow has a high-amplitude disturbance, the detonation wave fails to ignite due to insufficient energy.

Testing the complexity and chaotic nature of wave-dominated turbulent flows

This study employs chaos theory concepts to investigate the complexity and chaotic nature of surface-generated waves in comparison to steady-flow state. The flows are generated through laboratory-flume experiments. Two separate tests are carried out: steady-flow test; and test with the addition of waves to the steady-flow under identical flow conditions. For wave addition, two-different frequencies of wave (i.e., 1.1 Hz and 2.1 Hz) are considered. Two chaos theory-based approaches are employed to determine the complexity and chaotic nature of wave-current dynamics: False nearest neighbour (FNN) method; and Lyapunov exponent method. An effort is also made to confirm and understand the results from these chaos-theory methods with the results from the wavelet and Shannon entropy methods. The results from the chaos-theory methods suggest optimistic evidence of the presence of chaotic behaviour in the combined wave-current cases. The results show that greater complexity is found at the near-bed region with aperiodic nature of the eddy scales and the complexity becomes less when moving towards the near-surface region. These results are well supported by the tests with wavelet and Shannon entropy methods. The results further reveal that the complexity for the steady-flow case is high and the complexity decreases with the addition of waves.

Modelling Time-Dependent Flow through Railway Ballast

In this study, a numerical simulation of fluid flow through railway ballast in the time domain is presented, providing a model for unsteady-state flow. It is demonstrated that the position of the free surface with respect to time can also be used to solve the steady flow case. The effect of ballast fouling is included in the model to capture the realistic behavior of railway ballast, which is critical to understanding the impact of flooding. A thorough comparison with a range of previous studies, including theoretical and experimental approaches, is made, and very close agreement is obtained. The significant impact of ballast fouling on fluid flow and its potential consequences for railway infrastructure are highlighted by the simulation. Valuable insights into the behavior of water flow through porous media and its relevance to railway ballast management are offered by this study.

Thermal Performance of a Counter-Flow Double-Pass Solar Air Heater With The Steady and Pulsating Flow

The present work studied the influence of pulsating flow as an active method on the thermal performance of a double-pass solar air heater with a tubular solar absorber. A ball valve has been used as a pulse generator mounted at the downstream flow of the solar air heater. The experiments were under indoor conditions with a constant heat flux of 1000 W/m2, and different air mass flow rates ranged from 0.01 to 0.03 kg/s. Moreover, the study covered three pulsation frequencies varied from 1 to 3 Hz. Based on the experimental outcomes, it can be observed that the heat transfer rate is enhanced by applying the pulsating flow, where it was found that the outlet temperature in the case of applying the pulsating flow rises by about 25.6 - 27% as compared with the steady flow case. Moreover, pulsating flow offers a higher effective thermal performance by about 15.2 % at the maximum air mass flow rate compared with the steady flow. In addition, the findings pointed out that varying the pulsation frequency from 1 to 3 Hz produces an enhancement in heat transfer rate and in solar heater effective efficiency, where it was found that when changing the frequency from 1 to 3, the increment of effective efficiency ranged from 3.8 to 6.9 % depending on the air mass flow value.

Linear instabilities of pulsatile plane channel flow between compliant walls

The linear dynamics of perturbations developing in an infinite channel with compliant walls is investigated for pulsatile flow conditions. Two-dimensional modal perturbations are considered for Womersley-type pulsating base flows and the wall motion is only allowed in the normal direction. It is found that the flow dynamics is mainly governed by four control parameters: the Reynolds number $Re$ , the reduced velocity $V_R$ , the Womersley number $Wo$ and the amplitude of the base-flow modulation $\tilde Q$ . Linear stability analyses are carried out within the framework of Floquet theory, implementing an efficient approach for removing spurious eigenmodes. The characteristics of flow-based (Tollmien–Schlichting) and wall-based (both travelling-wave flutter and divergence) modes are investigated over a large control-parameter space. It is shown that travelling-wave flutter (TWF) modes are predominantly influenced by the reduced velocity and that the Reynolds number has only a marginal effect. The critical reduced velocity (corresponding to onset of linear instability) is demonstrated to depend both on the Womersley number and modulation amplitude for a given set of wall parameters. Similarly to the steady flow case, the Tollmien–Schlichting (TS) mode is also found to be only weakly affected by the flexibility of the wall in pulsatile flow conditions. Finally, the classification given by Benjamin (J. Fluid Mech., vol. 16, 1963, pp. 436–450) is found to be too restrictive in the case of pulsatile base flows. In particular, a new type of mode is identified that shares characteristics of two distinct Floquet eigenmodes: TS and TWF modes. Due to coupling of the different Floquet harmonics, a phenomenon specific to time-periodic base flows, this two-wave mode exhibits a beating over the intracyclic dynamics.

Assessment of polymer feedback coupling approaches in simulation of viscoelastic fluids using the lattice Boltzmann method

The lattice Boltzmann method (LBM) has emerged as a popular approach to simulate complex non-Newtonian fluid mechanics problems, such as those involving viscoelastic fluids. In these cases, the LBM solver often requires modelling additional non-linear and non-ideal hydrodynamic contributions to the flow resulting from the addition of polymer molecules. These additional contributions, commonly referred to as the polymer feedback, are either incorporated directly in the Boltzmann stress (pressure method) or alternatively in the momentum field as a forcing contribution (forcing method). The choice of coupling approach has in previous works been strongly speculated to impact the numerical accuracy and stability of LBM solvers for viscoelastic fluids, however no formal comparison has been conducted. In this work, we compare the two popular polymer feedback coupling approaches, as well as the various schemes commonly used in the forcing method, namely, the explicit forcing (EF) scheme (He et al., 1998), the Guo et al. scheme (Guo et al., 2002), and the exact-difference-method (EDM) (Kupershtokh et al., 2009). First, through a theoretical comparison in recovering the macroscopic equations, we show that both the EF and Guo et al. forcing schemes retain the exact hydrodynamic representation up to second-order, whereas the EDM scheme includes additional lattice error terms. Similarly, the pressure method retains additional unphysical terms, which are known to violate Galilean invariance. Through numerical tests, involving the steady and time-dependent Poiseuille flow cases, as well as the four-roll mill problem under both steady and transient regimes, it is shown that the forcing method has enhanced numerical accuracy and stability at varying elastic effects. Alternatively, owing to the violation of Galilean invariance, the pressure method proved to be unstable even for moderate elastic effects. The obtained results are useful in determining the optimal polymer feedback coupling approach when simulating viscoelastic fluids using LBM, as well as further testing the generality of previous LBM findings on non-ideal coupling approaches for viscoelastic fluids.

Exact Solutions to Navier–Stokes Equations Describing a Gradient Nonuniform Unidirectional Vertical Vortex Fluid Flow

The paper announces a family of exact solutions to Navier–Stokes equations describing gradient inhomogeneous unidirectional fluid motions (nonuniform Poiseuille flows). The structure of the fluid motion equations is such that the incompressibility equation enables us to establish the velocity defect law for nonuniform Poiseuille flow. In this case, the velocity field is dependent on two coordinates and time, and it is an arbitrary-degree polynomial relative to the horizontal (longitudinal) coordinate. The polynomial coefficients depend on the vertical (transverse) coordinate and time. The exact solution under consideration was built using the method of indefinite coefficients and the use of such algebraic operations was for addition and multiplication. As a result, to determine the polynomial coefficients, we derived a system of simplest homogeneous and inhomogeneous parabolic partial equations. The order of integration of the resulting system of equations was recurrent. For a special case of steady flows of a viscous fluid, these equations are ordinary differential equations. The article presents an algorithm for their integration. In this case, all components of the velocity field, vorticity vector, and shear stress field are polynomial functions. In addition, it has been noted that even without taking into account the thermohaline convection (creeping current) all these fields have a rather complex structure.

Experimental and Theoretical Study of Surge Behavior in a Boil-Off Gas Centrifugal Compressor on an LNG Carrier

In this study, we conducted experiments and numerical analysis on a 270 kW industrial-scale centrifugal compressor used in the fuel supply system of an LNG carrier in order to improve the lumped-parameter surge model by considering the viscosity in the pipeline and to confirm whether the improved model would be applicable. The steady and unsteady compressor performance curves were constructed using measurements and predictions, respectively. The flow through the pipeline was assumed to be both steady and unsteady, and each governing equation under the assumptions was derived in accordance with the lumped-parameter model. In the steady flow case, the surge behavior of the modified model was in a good agreement with the lumped-parameter model at surge parameter B = 4.8716. In the unsteady flow case, however, the modified model showed a deviation from the lumped-parameter model, and the simulation from the modified model described the surge behavior 5% more accurately than the lumped-parameter model. Through experiments and numerical analysis, this study showed that the present models are useful and applicable for describing the surge behavior of an industrial-scale single-stage centrifugal compressor.

Numerical investigation on buoyancy-driven flow over a circular cylinder in a channel with nonparallel walls

The flow over a stationary cylinder is destabilized mainly in the separated shear layer (SSL) and wake regions and transverse shedding forms downstream of the cylinder. It has been revealed that aiding/opposing buoyancy or diverging/converging channel could stabilize/destabilize the separated and wake flows, while their joint effect has not been investigated. In this numerical investigation, the low-Re buoyancy-driven flow over a circular cylinder placed in a channel consists of diverging or converging walls is studied. We aim to reveal the effects of two parameters, i.e. Richardson number Ri representing buoyancy and channel angle α, on the flow dynamics and heat transfer. The cases of aiding (Ri > 0) and opposing (Ri < 0) buoyancy, and diverging (α > 0) and converging (α < 0) channels are considered in this work. The effects of Ri and α are presented in terms of the instantaneous flow field, global characteristic quantities, heat transfer between the fluid and cylinder surface, separation and shearing of wake flow, and temperature variations of the channel walls. We found that the flow transits to steady subjected to aiding buoyancy at a lower value of Ri where the fluctuations of force and heat transfer rate reduces to zero; however, the mean drag and Nusselt number do not exhibit monotonic variations with Ri. The convection heat transfer is prone to be notably affected by Ri rather than α on the leeward side of the cylinder, and its fluctuation is significant in the separation zone downstream of the separation point. Depending on the value of Ri, the steady and unsteady wake flow present different velocity profile as well as the shearing; the deficit wake flow recovers more rapidly for the steady flow cases.

In Vitro Pressure Measurements Across an Interatrial Shunt for HFpEF Treatment.

Preserved ejection fraction heart failure (HFpEF) can be treated by installing a shunt in the interatrial septum, which relieves excess pressure in the left atrium by allowing blood to flow from left to right. This technique has proven effective in clinical trials, but the details of the flow through the shunted heart are not well understood. The current study aims to collect quantitative data on the relationship between pressure and flow rate in such shunts. An in vitro, shunted double atrium flow phantom was fabricated and used to investigate the relationship between pressure drop and flow across an interatrial shunt. Flow rate was controlled and the resulting pressure drop across the shunt was measured for a variety of flow cases, including steady and pulsatile flow, flow rate waveforms typical of healthy and failing hearts, and low and high heart rates. The results show a positive relationship between shunt flow rate and pressure drop which is more pronounced in steady flow than in pulsatile flow. Increasing heart rate increases the time-averaged pressure drop across the shunt but not the maximum pressure drop. For steady-flow cases, large changes in pressure drop resulting from moderate changes in flow rate suggest a flow regime transition during parts of the cardiac cycle. Comparison of time-averaged pulsatile flow pressure measurements with steady-flow measurements and two analytical plate-orifice models suggests that none approximate pulsatile flow accurately. The flow rate/pressure drop relationship across an in vitro model of an interatrial shunt has been measured for a variety of physiologically relevant cases. Among other things, the results suggest that steady flow approximations to the heart's pulsatile flow should be used with caution and simplified theoretical models do not approximate the flow rate/pressure drop relationship accurately.

On the steady laminar boundary-layer flow over the family of rotating tori

We present the laminar boundary layer flow over the standard tori rotating in an otherwise quiescent incompressible fluid. The mean flow equations are derived by using a specific form of the toroidal–poloidal coordinate system. A pseudo-spectral collocation method mixed with the first-order Euler scheme is employed to solve the related steady mean flow equations. A fixed aspect ratio R is defined as the ratio of radii of the torus and its tube that distinguishes each fixed torus within the general family of standard tori. Tori of spindle type turn into a sphere when the aspect ratio R=0 and thus the mean flow results for the rotating sphere are consistently reproduced. The influence of aspect ratio R on the steady mean velocity components and the wall shear stresses are presented in detail. In a precise manner, the mean flow results for each rotating torus are explainable in a consistent way to the already established case of steady laminar boundary layer flow over the rotating sphere. The significance of the current work is discussed in terms of the hydrodynamic instabilities of boundary layer flows over the family of rotating tori.

Experimental investigation of gradually-varied unsteady flow passed a circular pile

The flow field around a vertically mounted circular pile exposed to gradually-varied unsteady flow is experimentally investigated. Experiments were conducted in a 30 m-long and 1 m-wide recirculating flume equipped with a variable discharge pump. A circular cylinder with 9 cm diameter was used as the model pile. To understand the influence of accelerating and decelerating flow conditions, three unsteady cases with different unsteadiness degrees were tested as well as a reference steady flow case. The spatial and temporal variations of Reynolds-averaged velocity and turbulence characteristics around the pile, as well as undisturbed flow, were analyzed. Findings show that there are distinct differences between the tested gradually-varied unsteady flow cases and the reference steady flow case. The Reynolds-averaged velocity vs. TKE plots of undisturbed flow indicated a hysteresis effect, such that larger turbulence is generated during the falling stage of the flow compared to the rising stage. This hysteresis was considerably reduced in the pile wake, and even reversed hysteresis was seen at certain cases. The spatial variation of Reynolds-averaged velocity and turbulence in the peak instant of unsteady flow was qualitatively similar to that of steady flow, but quantitatively, turbulence, flow contraction, and velocity deficit in the near-wake region were smaller in the case of unsteady flow. Contrarily, the unsteady flow generated remarkably higher turbulence levels at further downstream in the pile wake. It is concluded that in the case of unsteady flow the pile behaves as if it has a more streamlined shape. The results were also interpreted from structure-bed interactions perspective, explaining the differences between the pile scour induced by steady and unsteady flow conditions.