Flow characteristics and flow-induced forces of a stationary and rotating triangular cylinder with different incidence angles at low Reynolds numbers

In this paper, the problem of two-dimensional fluid flow past a stationary and rotationally oscillating equilateral triangular cylinder with a variable incident angle, Reynolds number, oscillating amplitude, and oscillating frequency is numerically investigated. The computations are carried out by using a two-step Taylor-characteristic-based Galerkin (TCBG) algorithm. For the stationary cases, simulations are conducted at various incident angles of α=0.0–60.0° and Reynolds numbers of Re=50–160. For the oscillation cases, the investigations are done at various oscillating amplitudes of θ_max=7.5–30.0° and oscillating frequencies of F_s/F_o=0.5–3.0 considering two different incidence angles (α=0.0°, 60.0°) and three different Reynolds numbers (Re=50, 100, 150). The results show that the influences of key parameters (incidence angle, Reynolds number, oscillating amplitude, and oscillating frequency) are significant on the flow pattern and hydrodynamic forces. For the stationary cases, at smaller angle of incidence (α≤30.0°), Reynolds number has a large impact on the position of the separation points. When α is between 30.0° and 60.0°, it was found that the separation points are located at the rear corners. From a topological point of view, the diagram of flow pattern is summarized, including two distinct patterns, namely, main separation and vortex merging. A deep analysis of the influence of Reynolds number and incidence angles on the mean pressure coefficient along the triangular cylinder surface is presented. Additionally, for the oscillating cases, the lock-on phenomenon is captured. The dominant flow patterns are 2S mode and P+S mode in lock-on region at α=0.0°. It is found at α=60.0°, however, that the flow pattern is predominantly 2S mode. Furthermore, except for the case of F_s/F_o=2.0, the mean drag decreases as the oscillating amplitude increases for each Reynolds number at α=0.0°. At α=60.0°, the minimum mean drag for F_s/F_o=1.5 is lower than that for stationary case, and occurs at θ_max=15.0° (Re=100) and θ_max=22.5° (Re=150), respectively. Finally, the effect of Reynolds number on a rotational oscillation cylinder is elucidated.     

Moving boundary in a non-Newtonian fluid

Moving boundary value problem in non-Newtonian fluid is considered. Exact analytical solution for the flow of second-grade fluid for a rigid moving plate oscillating in its own plane, is obtained. The Doppler effect has been observed due to the motion of the plate. The shearing stress on the plate is also calculated. It is concluded that the solutions for stationary porous boundaries can be obtained from the solutions of moving rigid boundaries.     

Hall effects on the magnetohydrodynamic shear flow past an infinite porous flat plate subjected to uniform suction or blowing

An analysis is made of Hall effects on the steady shear flow of a viscous incompressible electrically conducting fluid past an infinite porous plate in the presence of a uniform transverse magnetic field. It is shown that for suction at the plate, steady shear flow solution exists only when S²<Q, where S and Q are the suction and magnetic parameters, respectively. The primary flow velocity decreases with increase in Hall parameter m. But the cross-flow velocity first increases and then decreases with increase in m. Similar results are obtained for variation of the induced magnetic field with m. It is further found that for blowing at the plate, steady shear flow solution exists only when , where S₁ is the blowing parameter.     

Experimental study on the flow regimes past a confined prism undergoing self-sustained oscillations

An experimental study based on Particle Image Velocimetry (PIV) is presented with the objective of studying the flow regimes that appear in the flow past a confined prism undergoing self-sustained oscillations at low Reynolds numbers (Re). The square-section prism, placed inside a 3D square cross-section vertical channel with a confinement ratio of 1/2.5, was tethered to the channel walls and, therefore, it was allowed to move freely transverse to the incoming flow. Re (based on the prism cross-section height) was varied in the range from 100 to 700. Three different prism to fluid density ratios (m¹) were considered: 0.56, 0.70, and 0.91. These two parameters, Re and m¹, were used to map the results obtained. In particular, it was found that five different regimes appear: (1) steady prism with steady recirculation bubble, (2) steady prism with unsteady vortex shedding wake, (3) large amplitude low frequency oscillating prism with unsteady vortex shedding wake, (4) small amplitude high frequency oscillating prism with unsteady vortex shedding wake, and (5) irregular/chaotic motion of both the prism and the wake. The PIV results and associated numerical simulations were used to analyze the different prism and wake states.     

Identification of Fractional Systems Using an Output-Error Technique

An original method for modeling, simulation and identification of fractional systems in the time domain is presented in this article. The basic idea is to model the fractional system by a state-space representation, where conventional integration is replaced by a fractional one with the help of a non-integer integrator. This operator is itself approximated by a N-dimensional system composed of an integrator and of a phase-lead filter. An output-error technique is used in order to estimate the parameters of the model, including the fractional order N. Simulations exhibit the properties of the identification algorithm. Finally, this methodology is applied to the modeling of the dynamics of a real heat transfer system.     

Purely elastic flow asymmetries in flow-focusing devices

The flow of a viscoelastic fluid through a microfluidic flow-focusing device is investigated numerically with a finite-volume code using the upper-convected Maxwell (UCM) and Phan-Thien–Tanner (PTT) models. The conceived device is shaped much like a conventional planar “cross-slot” except for comprising three inlets and one exit arm. Strong viscoelastic effects are observed as a consequence of the high deformation rates. In fact, purely elastic instabilities that are entirely absent in the corresponding Newtonian fluid flow are seen to occur as the Deborah number (De) is increased above a critical threshold. From two-dimensional numerical simulations we are able to distinguish two types of instability, one in which the flow becomes asymmetric but remains steady, and a subsequent instability at higher De in which the flow becomes unsteady, oscillating in time. For the UCM model, the effects of the geometric parameters of the device (e.g. the relative width of the entrance branches, WR) and of the ratio of inlet average velocities (VR) on the onset of asymmetry are systematically examined. We observe that for high velocity ratios, the critical Deborah number is independent of VR (e.g. De_c ≈ 0.33 for WR = 1), but depends non-monotonically on the relative width of the entrance branches. Using the PTT model we are able to demonstrate that the extensional viscosity and the corresponding very large stresses are decisive for the onset of the steady-flow asymmetry.     

Visualization and structure of confined, milliscale, unsteady impinging slot jets and associated vortices

Flow characteristics of confined, laminar milliscale slot jets are investigated from visualizations, as they impinge upon a flat target plate, with a fully developed velocity profile at the nozzle exit. The effects of Reynolds number Re and normalized nozzle-to-plate distance H/B are considered for a nozzle width B of 1.0 mm. Transition from a stable symmetric jet to an unsteady oscillating jet is observed as the Reynolds number increases (with H/B constant), where the Reynolds number associated with this transition decreases as the normalized nozzle-to-plate distance H/B increases. Instantaneous visualizations show unsteady lateral distortions of jet columns at experimental conditions corresponding to the presence of continuous sinusoidal oscillations, intermittent oscillating motion of the jet column, and jet flow fluctuation/flapping motion. Also apparent in flow visualization sequences are smoke signatures associated with instantaneous vortex structures, which form as secondary flows develop in fluid which, initially, is just adjacent to and within the jet column. Associated jet and vortex structural changes are described as different modes of unsteadiness are present, including characterization of jet column unsteadiness using jet column oscillation frequency, and lateral and streamwise extents of jet distortion.     

Combined forced and free convective flow of water at 4°C past a semi-infinite vertical plate

Effects of buoyancy forces on forced and free convective flow of water at 4°C past a semi-infinite vertical plate at constant temperature are studied. Flow is assumed to be vertically upwards. Similarity solutions are derived and the resulting equations are solved numerically on a computer. Velocity and temperature profiles are shown graphically and numerical values of the skin friction and the rate of heat transfer are entered in tables. It is observed that the skin friction and the Nusselt number increase with increasing Gr/Re², where Gr is the Grashof number and Re is the Reynolds number     

Time periodic electroosmotic flow of the generalized Maxwell fluids between two micro-parallel plates

Analytical solutions are presented using method of separation of variables for the time periodic EOF flow of linear viscoelastic fluids between micro-parallel plates. The linear viscoelastic fluids used here are described by the general Maxwell model. The solution involves analytically solving the linearized Poisson–Boltzmann equation, together with the Cauchy momentum equation and the general Maxwell constitutive equation. By numerical computations, the influences of the electrokinetic width K denoting the characteristic scale of half channel width to Debye length, the periodic EOF electric oscillating Reynolds number Re and normalized relaxation time λ₁ω on velocity profiles and volumetric flow rates are presented. Results show that for prescribed electrokinetic width K, lower oscillating Reynolds number Re and shorter relaxation time λ₁ω reduces the plug-like EOF velocity profile of Newtonian fluids. For given Reynolds number Re and electrokinetic width K, longer relaxation time λ₁ω leads to rapid oscillating EOF velocity profiles with increased amplitude. With the increase of the K, the velocity variations are restricted to a very narrow region close to the EDL for small relaxation time. However, with the increase of the relaxation time, the elasticity of the fluid becomes conspicuous and the velocity variations can be expanded to the whole flow field. As far as volume flow rates are concerned, for given electrodynamic width K, larger oscillating Reynolds number Re results in a smaller volume flow rates. For prescribed oscillating Reynolds number Re, with the changes of relaxation time λ₁ω, volume flow rates will produce some peaks no matter how the electrodynamic width K varies. Moreover, the time periodic evolution of the velocity profiles provides a detail insight of the flow characteristic of this flow configuration.     

Exact analytical solutions for axial flow of a fractional second grade fluid between two coaxial cylinders

The velocity field and the adequate shear stress corresponding to the longitudinal flow of a fractional second grade fluid, between two infinite coaxial circular cylinders, are determined by applying the Laplace and finite Hankel transforms. Initially the fluid is at rest, and at time t = 0⁺, the inner cylinder suddenly begins to translate along the common axis with constant acceleration. The solutions that have been obtained are presented in terms of generalized G functions. Moreover, these solutions satisfy both the governing differential equations and all imposed initial and boundary conditions. The corresponding solutions for ordinary second grade and Newtonian fluids are obtained as limiting cases of the general solutions. Finally, some characteristics of the motion, as well as the influences of the material and fractional parameters on the fluid motion and a comparison between models, are underlined by graphical illustrations.     

Mathematical modeling of fractional derivatives for magnetohydrodynamic fluid flow between two parallel plates by the radial basis function method

Investigations into the magnetohydrodynamics of viscous fluids have become more important in recent years, owing to their practical significance and numerous applications in astro-physical and geo-physical phenomena. In this paper, the radial base function was utilized to answer fractional equation associated with fluid flow passing through two parallel flat plates with a magnetic field. The magnetohydrodynamics coupled stress fluid flows between two parallel plates, with the bottom plate being stationary and the top plate moving at a persistent velocity. We compared the radial basis function approach to the numerical method (fourth-order Range-Kutta) in order to verify its validity. The findings demonstrated that the discrepancy between these two techniques is quite negligible, indicating that this method is very reliable. The impact of the magnetic field parameter and Reynolds number on the velocity distribution perpendicular to the fluid flow direction is illustrated. Eventually, the velocity parameter is compared for diverse conditions α, Reynolds and position (y), the maximum of which occurs at α = 0.4. Also, the maximum velocity values occur in α=0.4 and Re=1000 and the concavity of the graph is less for α=0.8.     

Laminar bubbly flow in an open capillary channel in microgravity

The study of a bubbly laminar two-phase flow in an open capillary channel under microgravity conditions was conducted aboard the sounding rocket, Texus-45. The channel consists of two parallel plates of width b = 25 mm and distance a = 10 mm. The flow along the length l = 80 mm is confined by a free surface on one side and a plate on the opposite side. The bubbles are injected at the nozzle of the capillary channel via six capillary tubes of 100 μm in inner diameter. Different liquid and gas flow rates were tested leading to different liquid free surface shape and bubble size.     

Similarity solution of mixed convection over a horizontal moving plate

Investigation to the mixed convective heat and mass transfer over a horizontal plate has been carried out. By applying transformation group theory to analysis of the governing equations of continuity, momentum, energy and diffusion, we show the existence of similarity solution for the problem provided that the temperature and concentration at the wall are proportional to x ^4/(7-5n) and that the moving speed of the plate is proportional to x ^(3-n)/(7-5n), and further obtain a similarity representation of the problem. The similarity equations have been solved numerically by a fourth-order Runge–Kutta scheme. The numerical results obtained for Pr=0.72 and various values of the parameters Sc, K ₁, K ₂ and K ₃ reveals the influence of the parameters on the flow, heat and mass transfer behavior.     

Visco-elastic flow due to a plate which starts oscillating in the presence of a parallel stationary plate

In this paper, the flow of a visco-elastic liquid between two parallel plates has been studied when one plate is stationary and the other plate suddenly starts oscillating. Both finite Fourier sine transform and Laplace transform technique have been employed to solve the basic differential equations. The flow phenomenon has been characterized by the parameters, and and the effects of these on the flow characteristics have been studied through several graphs.Late professor of the department, who died in an accident on 7th July 1978.     

Three-dimensional exact analysis of piezoelectric laminated plates via a sampling surfaces method

A paper focuses on the use of the efficient approach to three-dimensional (3D) exact solutions of electroelasticity for piezoelectric laminated plates. This approach is based on the new method of sampling surfaces (SaS) developed recently by the authors. We introduce inside the nth layer I_n not equally spaced SaS parallel to the middle surface of the plate and choose displacements of these surfaces as basic plate variables. Such an idea permits the representation of the proposed piezoelectric plate formulation in a very compact form. This fact gives the opportunity to derive the 3D exact solutions of electroelasticity for thick and thin piezoelectric laminated plates with a specified accuracy utilizing a sufficient number of SaS, which are located at interfaces and Chebyshev polynomial nodes.     

Numerical solution of the turbulent flow of an incompressible fluid in curved planar and axially symmetrical channels

The governing equations for axially symmetric flow, where the Reynolds stresses are expressed by scalar turbulent viscosity, are the Reynolds equations. The turbulence model k, ? is used in the well-known form for fully developed turbulent flow.The numerical method, a continuation of the MAC system¹, is adapted so that even for high Reynolds cell numbers precision (δx²) can be achieved for the steady flow. Irregular cells join the rectangular network on the curved surface. Von Neumann's stability condition of the linearised numerical system is investigated. Special problems concerning the numerical solution of the turbulence model equations are stated and a special procedure is worked out to ensure that the fields k, ? do not converge to physically meaningless values. The program for the computer is universal in that the boundary problems can be assigned by input data.As an example, an axially symmetrical diffuser with an area ratio of widening 1.40 is computed. Fields of velocity and pressure at the wall as well as fields v_T and k are assessed. The results are compared with an experiment. The conclusion is that this method is suitable for the problems mentioned in this study as well as for unsteady flow.     

动失速型非定常分离涡结构的控制

本文对动失速型非定常分离涡结构的控制方法，在低速风洞中应用相平均测压技术进行了实验研究。在二元平板模型中部安装一作俯仰振荡的扰流板产生动失速型分离涡，在其上游安装另一用作控制的小扰流板。实验结果表明，应用前置的振荡小扰流板可影响并改变动失速分离涡的强度和对流特性。在最有利的控制相位下，涡吸力峰可降低４８％，涡对流时间可以推迟０．１９周期。对于间歇式振荡扰流板，采用相位提前控制方式比相位滞后控制方式更有效。     

3-D measurements of separated flow over rigid plates with spanwise periodic cambering at low Reynolds number

Measurements of the flow field around a flat plate and rigid plates with spanwise periodic cambering were performed using volumetric three-component velocimetry (V3V) at a Reynolds numbers of 28,000 at α=12° where the flow is fully separated. The Reynolds normal and shear stresses, and the streamwise, spanwise and normal components of the vorticity vector are investigated for three-dimensionality. Flow features are discussed in context of the periodic cambering and corresponding aerodynamic force measurements. The periodic cambering results in spanwise variation in the reversed-flow region, Reynolds stresses and spanwise vorticity. These spanwise variations are induced by streamwise and normal vortices of opposite directions of rotation. Moreover, measurements were carried out for the cambered plates at α=8°, where a long separation bubble exists, to further understand the behavior of the streamwise and normal vortices. These vortices become more organized and increase in strength and size at the lower angle of attack. It is also speculated that these vortices contribute to the increase in lift at and beyond the onset of stall angle of attack.     

A low Reynolds number variant of partially-averaged Navier-Stokes model for turbulence

A low Reynolds number (LRN) formulation based on the Partially Averaged Navier-Stokes (PANS) modelling method is presented, which incorporates improved asymptotic representation in near-wall turbulence modelling. The effect of near-wall viscous damping can thus be better accounted for in simulations of wall-bounded turbulent flows. The proposed LRN PANS model uses an LRN k-ε model as the base model and introduces directly its model functions into the PANS formulation. As a result, the inappropriate wall-limiting behavior inherent in the original PANS model is corrected. An interesting feature of the PANS model is that the turbulent Prandtl numbers in the k and ε equations are modified compared to the base model. It is found that this modification has a significant effect on the modelled turbulence. The proposed LRN PANS model is scrutinized in computations of decaying grid turbulence, turbulent channel flow and periodic hill flow, of which the latter has been computed at two different Reynolds numbers of Re = 10,600 and 37,000. In comparison with available DNS, LES or experimental data, the LRN PANS model produces improved predictions over the standard PANS model, particularly in the near-wall region and for resolved turbulence statistics. Furthermore, the LRN PANS model gives similar or better results - at a reduced CPU time - as compared to the Dynamic Smagorinsky model.     

Unsteady flow of a fourth-grade fluid due to an oscillating plate

The governing non-linear high-order, sixth-order in space and third-order in time, differential equation is constructed for the unsteady flow of an incompressible conducting fourth-grade fluid in a semi-infinite domain. The unsteady flow is induced by a periodically oscillating two-dimensional infinite porous plate with suction/blowing, located in a uniform magnetic field. It is shown that by augmenting additional boundary conditions at infinity based on asymptotic structures and transforming the semi-infinite physical space to a bounded computational domain by means of a coordinate transformation, it is possible to obtain numerical solutions of the non-linear magnetohydrodynamic equation. In particular, due to the unsymmetry of the boundary conditions, in numerical simulations non-central difference schemes are constructed and employed to approximate the emerging higher-order spatial derivatives. Effects of material parameters, uniform suction or blowing past the porous plate, exerted magnetic field and oscillation frequency of the plate on the time-dependent flow, especially on the boundary layer structure near the plate, are numerically analysed and discussed. The flow behaviour of the fourth-grade non-Newtonian fluid is also compared with those of the Newtonian fluid.