Pressure-driven flow in a two-dimensional channel with porous walls

2009 ◽  
Vol 631 ◽  
pp. 1-21 ◽  
Author(s):  
QUAN ZHANG ◽  
ANDREA PROSPERETTI

The finite-Reynolds-number two-dimensional flow in a channel bounded by a porous medium is studied numerically. The medium is modelled by aligned cylinders in a square or staggered arrangement. Detailed results on the flow structure and slip coefficient are reported. The hydrodynamic force and couple acting on the cylinder layer bounding the porous medium are also evaluated as a function of the Reynolds number. In particular, it is shown that, at finite Reynolds numbers, a lift force acts on the bodies, which may be significant for the mobilization of bottom sediments.

2011 ◽  
Vol 679 ◽  
pp. 77-100 ◽  
Author(s):  
QIANLONG LIU ◽  
ANDREA PROSPERETTI

The finite-Reynolds-number three-dimensional flow in a channel bounded by one and two parallel porous walls is studied numerically. The porous medium is modelled by spheres in a simple cubic arrangement. Detailed results on the flow structure and the hydrodynamic forces and couple acting on the sphere layer bounding the porous medium are reported and their dependence on the Reynolds number illustrated. It is shown that, at finite Reynolds numbers, a lift force acts on the spheres, which may be expected to contribute to the mobilization of bottom sediments. The results for the slip velocity at the surface of the porous layers are compared with the phenomenological Beavers–Joseph model. It is found that the values of the slip coefficient for pressure-driven and shear-driven flow are somewhat different, and also depend on the Reynolds number. A modification of the relation is suggested to deal with these features. The Appendix provides an alternative derivation of this modified relation.


1972 ◽  
Vol 39 (4) ◽  
pp. 873-878 ◽  
Author(s):  
J. W. McDonald ◽  
V. E. Denny ◽  
A. F. Mills

Numerical solutions of the Navier-Stokes equations are obtained for steady two-dimensional flow in the inlet region of both a tube and a channel. The entering flow is considered to be either uniform (u = constant, v = 0) or irrotational (u = constant, ω = 0). Values of Reynolds number Re = u0a/ν range from 0.75 to 2 × 106. An improved method for solving the stream function-vorticity equations of hydrodynamics has been developed. The method is stable at all Reynolds numbers and appears to be computationally superior to previous methods.


2002 ◽  
Vol 473 ◽  
pp. 167-190 ◽  
Author(s):  
DWIGHT BARKLEY ◽  
M. GABRIELA M. GOMES ◽  
RONALD D. HENDERSON

Results are reported from a three-dimensional computational stability analysis of flow over a backward-facing step with an expansion ratio (outlet to inlet height) of 2 at Reynolds numbers between 450 and 1050. The analysis shows that the first absolute linear instability of the steady two-dimensional flow is a steady three-dimensional bifurcation at a critical Reynolds number of 748. The critical eigenmode is localized to the primary separation bubble and has a flat roll structure with a spanwise wavelength of 6.9 step heights. The system is further shown to be absolutely stable to two-dimensional perturbations up to a Reynolds number of 1500. Stability spectra and visualizations of the global modes of the system are presented for representative Reynolds numbers.


2010 ◽  
Vol 1 (1-2) ◽  
pp. 15-20 ◽  
Author(s):  
B. Bolló

Abstract The two-dimensional flow around a stationary heated circular cylinder at low Reynolds numbers of 50 < Re < 210 is investigated numerically using the FLUENT commercial software package. The dimensionless vortex shedding frequency (St) reduces with increasing temperature at a given Reynolds number. The effective temperature concept was used and St-Re data were successfully transformed to the St-Reeff curve. Comparisons include root-mean-square values of the lift coefficient and Nusselt number. The results agree well with available data in the literature.


Author(s):  
Francine Battaglia ◽  
George Papadopoulos

The effect of three-dimensionality on low Reynolds number flows past a symmetric sudden expansion in a channel was investigated. The geometric expansion ratio of in the current study was 2:1 and the aspect ratio was 6:1. Both experimental velocity measurements and two- and three-dimensional simulations for the flow along the centerplane of the rectangular duct are presented for Reynolds numbers in the range of 150 to 600. Comparison of the two-dimensional simulations with the experiments revealed that the simulations fail to capture completely the total expansion effect on the flow, which couples both geometric and hydrodynamic effects. To properly do so requires the definition of an effective expansion ratio, which is the ratio of the downstream and upstream hydraulic diameters and is therefore a function of both the expansion and aspect ratios. When the two-dimensional geometry was consistent with the effective expansion ratio, the new results agreed well with the three-dimensional simulations and the experiments. Furthermore, in the range of Reynolds numbers investigated, the laminar flow through the expansion underwent a symmetry-breaking bifurcation. The critical Reynolds number evaluated from the experiments and the simulations was compared to other values reported in the literature. Overall, side-wall proximity was found to enhance flow stability, helping to sustain laminar flow symmetry to higher Reynolds numbers in comparison to nominally two-dimensional double-expansion geometries. Lastly, and most importantly, when the logarithm of the critical Reynolds number from all these studies was plotted against the reciprocal of the effective expansion ratio, a linear trend emerged that uniquely captured the bifurcation dynamics of all symmetric double-sided planar expansions.


2017 ◽  
Vol 825 ◽  
pp. 631-650 ◽  
Author(s):  
Francesco Romanò ◽  
Arash Hajisharifi ◽  
Hendrik C. Kuhlmann

The topology of the incompressible steady three-dimensional flow in a partially filled cylindrical rotating drum, infinitely extended along its axis, is investigated numerically for a ratio of pool depth to radius of 0.2. In the limit of vanishing Froude and capillary numbers, the liquid–gas interface remains flat and the two-dimensional flow becomes unstable to steady three-dimensional convection cells. The Lagrangian transport in the cellular flow is organised by periodic spiralling-in and spiralling-out saddle foci, and by saddle limit cycles. Chaotic advection is caused by a breakup of a degenerate heteroclinic connection between the two saddle foci when the flow becomes three-dimensional. On increasing the Reynolds number, chaotic streamlines invade the cells from the cell boundary and from the interior along the broken heteroclinic connection. This trend is made evident by computing the Kolmogorov–Arnold–Moser tori for five supercritical Reynolds numbers.


Micromachines ◽  
2021 ◽  
Vol 12 (9) ◽  
pp. 1055
Author(s):  
Alireza Farahinia ◽  
Jafar Jamaati ◽  
Hamid Niazmand ◽  
Wenjun Zhang

One approach to achieve a homogeneous mixture in microfluidic systems in the quickest time and shortest possible length is to employ electroosmotic flow characteristics with heterogeneous surface properties. Mixing using electroosmotic flow inside microchannels with homogeneous walls is done primarily under the influence of molecular diffusion, which is not strong enough to mix the fluids thoroughly. However, surface chemistry technology can help create desired patterns on microchannel walls to generate significant rotational currents and improve mixing efficiency remarkably. This study analyzes the function of a heterogeneous zeta-potential patch located on a microchannel wall in creating mixing inside a microchannel affected by electroosmotic flow and determines the optimal length to achieve the desired mixing rate. The approximate Helmholtz–Smoluchowski model is suggested to reduce computational costs and simplify the solving process. The results show that the heterogeneity length and location of the zeta-potential patch affect the final mixing proficiency. It was also observed that the slip coefficient on the wall has a more significant effect than the Reynolds number change on improving the mixing efficiency of electroosmotic micromixers, benefiting the heterogeneous distribution of zeta-potential. In addition, using a channel with a heterogeneous zeta-potential patch covered by a slip surface did not lead to an adequate mixing in low Reynolds numbers. Therefore, a homogeneous channel without any heterogeneity would be a priority in such a range of Reynolds numbers. However, increasing the Reynolds number and the presence of a slip coefficient on the heterogeneous channel wall enhances the mixing efficiency relative to the homogeneous one. It should be noted, though, that increasing the slip coefficient will make the mixing efficiency decrease sharply in any situation, especially in high Reynolds numbers.


Author(s):  
Xuemin Ye ◽  
Weiping Yan ◽  
Chunxi Li

When liquid film is under evaporating or condensing conditions, the flow stability is clearly different to that under isothermal condition due to thermal non-equilibrium effect at interface, especially under lower Reynolds number. The universal linear temporal and spatial stability formulations of the two-dimensional surface waves on evaporating or isothermal or condensing liquid films are established in present paper with the collocation method based on the boundary layer theory and complete boundary conditions. The models include the effects of Reynolds number, thermocapillarity, inclination angle, liquid property, evaporation, isothermal or condensation. The effects of above factors are investigated with the neutral stability curves at different Reynolds numbers, and stabilities characteristics are fully indicated in theory for evaporating or condensing films.


1975 ◽  
Vol 70 (3) ◽  
pp. 529-542 ◽  
Author(s):  
Edmond Szechenyi

In wind-tunnel tests on bluff bodies the Reynolds number is often limited to values that are very much smaller than those of the flows being simulated. In such cases the experiments may have no practical significance whatsoever since both the fluctuating and the steady aerodynamic phenomena can vary considerably with Reynolds number.This difficulty was encountered in an investigation of supercritical incompressible flow over cylinders, and an attempt at artificially increasing the Reynolds number by means of surface roughness was made. In order to evaluate this simulation technique, the influence of various grades of surface roughness on the aerodynamic forces acting on cylinders of different diameters was studied over a wide range of Reynolds numbers in two very different wind tunnels. The results allow very positive conclusions to be drawn.


2011 ◽  
Vol 681 ◽  
pp. 411-433 ◽  
Author(s):  
HEMANT K. CHAURASIA ◽  
MARK C. THOMPSON

A detailed numerical study of the separating and reattaching flow over a square leading-edge plate is presented, examining the instability modes governing transition from two- to three-dimensional flow. Under the influence of background noise, experiments show that the transition scenario typically is incompletely described by either global stability analysis or the transient growth of dominant optimal perturbation modes. Instead two-dimensional transition effectively can be triggered by the convective Kelvin–Helmholtz (KH) shear-layer instability; although it may be possible that this could be described alternatively in terms of higher-order optimal perturbation modes. At least in some experiments, observed transition occurs by either: (i) KH vortices shedding downstream directly and then almost immediately undergoing three-dimensional transition or (ii) at higher Reynolds numbers, larger vortical structures are shed that are also three-dimensionally unstable. These two paths lead to distinctly different three-dimensional arrangements of vortical flow structures. This paper focuses on the mechanisms underlying these three-dimensional transitions. Floquet analysis of weakly periodically forced flow, mimicking the observed two-dimensional quasi-periodic base flow, indicates that the two-dimensional vortex rollers shed from the recirculation region become globally three-dimensionally unstable at a Reynolds number of approximately 380. This transition Reynolds number and the predicted wavelength and flow symmetries match well with those of the experiments. The instability appears to be elliptical in nature with the perturbation field mainly restricted to the cores of the shed rollers and showing the spatial vorticity distribution expected for that instability type. Indeed an estimate of the theoretical predicted wavelength is also a good match to the prediction from Floquet analysis and theoretical estimates indicate the growth rate is positive. Fully three-dimensional simulations are also undertaken to explore the nonlinear development of the three-dimensional instability. These show the development of the characteristic upright hairpins observed in the experimental dye visualisations. The three-dimensional instability that manifests at lower Reynolds numbers is shown to be consistent with an elliptic instability of the KH shear-layer vortices in both symmetry and spanwise wavelength.


Sign in / Sign up

Export Citation Format

Share Document