scholarly journals Stability of the Plane Bingham–Poiseuille Flow in an Inclined Channel

Fluids ◽  
2020 ◽  
Vol 5 (3) ◽  
pp. 141
Author(s):  
Paolo Falsaperla ◽  
Andrea Giacobbe ◽  
Giuseppe Mulone

We study the stability of laminar Bingham–Poiseuille flows in a sheet of fluid (open channel) down an incline with constant slope angle β∈(0,π/2). This problem has geophysical applications to the evolution of landslides. In this article, we apply to this problem recent results of Falsaperla et al. for laminar Couette and Poiseuille flows of Newtonian fluids in inclined channels. The stability of the basic motion of the generalised Navier–Stokes system for a Bingham fluid in a horizontal channel against linear perturbations has been studied. In this article, we study the flows of a Bingham fluid when the channel is oblique and we prove a stabilizing effect of the Bingham parameter B. We also study the stability of the linear system with an energy method (Lyapunov functions) and prove that the streamwise perturbations are always stable, while the spanwise perturbations are energy-stable if the Reynolds number Re is less than the critical Reynolds number Rc obtained solving a generalised Orr equation of a maximum variational problem.

Author(s):  
Jean-Yves Chemin ◽  
Benoit Desjardins ◽  
Isabelle Gallagher ◽  
Emmanuel Grenier

Let us now detail the stability properties of an Ekman layer introduced in Part I, page 11. First we will recall how to compute the critical Reynolds number. Then we will describe briefly what happens at larger Reynolds numbers. The first step in the study of the stability of the Ekman layer is to consider the linear stability of a pure Ekman spiral of the form where U∞ is the velocity away from the layer and ζ is the rescaled vertical component ζ = x3/√εν. The corresponding Reynolds number is Let us consider the Navier–Stokes–Coriolis equations, linearized around uE The problem is now to study the (linear) stability of the 0 solution of the system (LNSCε). If u=0 is stable we say that uE is linearly stable, if not we say that it is linearly unstable. Numerical results show that u=0 is stable if and only if Re<Rec where Rec can be evaluated numerically. Up to now there is no mathematical proof of this fact, and it is only possible to prove that 0 is linearly stable for Re<Re1 and unstable for Re>Re2 with Re1<Rec<Re2, Re1 being obtained by energy estimates and Re2 by a perturbative analysis of the case Re=∞. We would like to emphasize that the numerical results are very reliable and can be considered as definitive results, since as we will see below, the stability analysis can be reduced to the study of a system of ordinary differential equations posed on the half-space, with boundary conditions on both ends, a system which can be studied arbitrarily precisely, even on desktop computers (first computations were done in the 1960s by Lilly).


2001 ◽  
Vol 123 (4) ◽  
pp. 742-754 ◽  
Author(s):  
A. Pereira ◽  
G. McGrath ◽  
D. D. Joseph

The problem of predicting flow between rotating eccentric cylinders with axial throughput is studied. The system models a device used to test the stability of emulsions against changes in drop size distribution. The analysis looks for the major variation in flow properties which could put an emulsion at risk due to coalescence or breakage and finds the most likely candidate in the pressure gradient defined as the ratio of the difference between the maximum and minimum pressure to the arc length between the difference. The axial throughput is modeled by flow driven by a constant pressure gradient. The flow is calculated from the Navier-Stokes equation using the code SIMPLER (Patankar 1980). The effects of inertia at values typical for the device are studied. Several eccentricities and different rotational speeds are computed to sample the changes in flow and stress parameters in the idealized device for typical conditions. The numerical analysis is validated against the lubrication approximation in the low Reynolds number case. Conditions for stress induced cavitation are evaluated. The flow is completely determined by a Reynolds number, an eccentricity ratio and a dimensionless pressure gradient and all computed results are either presented or can be easily expressed in terms of these dimensionless parameters. The effect of inertia is to shift the eddy or re-circulation zone which develops in the more open region of the gap toward the region of low relative pressure; the zero of the relative pressure migrates away from the center and the distribution breaks the skew symmetry of the Stokes flow solution. The state of stress in the journal bearing is analyzed and a cavitation criterion based on the maximum tensile stress is compared with the traditional criterion based on pressure.


2012 ◽  
Vol 29 (1) ◽  
pp. 53-58 ◽  
Author(s):  
M. Bayareh ◽  
S. Mortazavi

AbstractThe equilibrium position of a deformable drop in Couette and Poiseuille flows is investigated numerically by solving the full Navier-Stokes equations using a finite difference/front tracking method. The objective of this work is to study the motion of a non-neutrally buoyant drop in Couette and Poiseuille flows with different density ratios at finite Reynolds numbers. Couette flow: The equilibrium position of the lighter drops is higher than the heavier drops at each particle Reynolds number. Also, the equilibrium position height increases with increasing the Reynolds number at a fixed density ratio. At this equilibrium distance from the wall, the migration velocity is zero, while the velocity of the drop in the flow direction and rotational velocity of the drop is finite. It is observed that the equilibrium position is independent of the initial position of the drops and depends on the density ratio and the shear Reynolds number. Poiseuille flow: When the drop is slightly buoyant, it moves to an equilibrium position between the wall and the centerline. The equilibrium position is close to the centerline if the drop lags the fluid but close to the wall if the drop leads the fluid. As the Reynolds number increases, the equilibrium position of lighter drops moves slightly closer to the wall and the equilibrium position of heavier drops moves towards the centerline.


2015 ◽  
Vol 112 (31) ◽  
pp. 9518-9523 ◽  
Author(s):  
Jianchun Wang ◽  
Qianxiao Li ◽  
Weinan E

The stability of the plane Poiseuille flow is analyzed using a thermodynamic formalism by considering the deterministic Navier–Stokes equation with Gaussian random initial data. A unique critical Reynolds number, Rec≈2,332, at which the probability of observing puffs in the solution changes from 0 to 1, is numerically demonstrated to exist in the thermodynamic limit and is found to be independent of the noise amplitude. Using the puff density as the macrostate variable, the free energy of such a system is computed and analyzed. The puff density approaches zero as the critical Reynolds number is approached from above, signaling a continuous transition despite the fact that the bifurcation is subcritical for a finite-sized system. An action function is found for the probability of observing puffs in a small subregion of the flow, and this action function depends only on the Reynolds number. The strategy used here should be applicable to a wide range of other problems exhibiting subcritical instabilities.


2015 ◽  
Vol 774 ◽  
pp. 1-4 ◽  
Author(s):  
Alexander J. Smits

Orlandi et al. (J. Fluid Mech., vol. 770, 2015, pp. 424–441) present direct numerical simulations over a very wide Reynolds number range for plane Couette and Poiseuille flows. The results reveal new information on the abrupt nature of transition in these flows, and the comparisons between Couette and Poiseuille flows help to provide a clearer picture of Reynolds number trends, especially with regard to inner/outer layer interactions. The stress distributions give strong support to Townsend’s attached eddy hypothesis, particularly for the wall-parallel component where there has been little experimental data available. The results pose some intriguing questions regarding the reconciliation of the present results with data at higher Reynolds numbers in different canonical flows.


2019 ◽  
Vol 14 (1) ◽  
pp. 52-58 ◽  
Author(s):  
A.D. Nizamova ◽  
V.N. Kireev ◽  
S.F. Urmancheev

The flow of a viscous model fluid in a flat channel with a non-uniform temperature field is considered. The problem of the stability of a thermoviscous fluid is solved on the basis of the derived generalized Orr-Sommerfeld equation by the spectral decomposition method in Chebyshev polynomials. The effect of taking into account the linear and exponential dependences of the fluid viscosity on temperature on the spectral characteristics of the hydrodynamic stability equation for an incompressible fluid in a flat channel with given different wall temperatures is investigated. Analytically obtained profiles of the flow rate of a thermovisible fluid. The spectral pictures of the eigenvalues of the generalized Orr-Sommerfeld equation are constructed. It is shown that the structure of the spectra largely depends on the properties of the liquid, which are determined by the viscosity functional dependence index. It has been established that for small values of the thermoviscosity parameter the spectrum compares the spectrum for isothermal fluid flow, however, as it increases, the number of eigenvalues and their density increase, that is, there are more points at which the problem has a nontrivial solution. The stability of the flow of a thermoviscous fluid depends on the presence of an eigenvalue with a positive imaginary part among the entire set of eigenvalues found with fixed Reynolds number and wavenumber parameters. It is shown that with a fixed Reynolds number and a wave number with an increase in the thermoviscosity parameter, the flow becomes unstable. The spectral characteristics determine the structure of the eigenfunctions and the critical parameters of the flow of a thermally viscous fluid. The eigenfunctions constructed in the subsequent works show the behavior of transverse-velocity perturbations, their possible growth or decay over time.


Author(s):  
Johannes Ruhland ◽  
Christian Breitsamter

AbstractThis study presents two-dimensional aerodynamic investigations of various high-lift configuration settings concerning the deflection angles of droop nose, spoiler and flap in the context of enhancing the high-lift performance by dynamic flap movement. The investigations highlight the impact of a periodically oscillating trailing edge flap on lift, drag and flow separation of the high-lift configuration by numerical simulations. The computations are conducted with regard to the variation of the parameters reduced frequency and the position of the rotational axis. The numerical flow simulations are conducted on a block-structured grid using Reynolds Averaged Navier Stokes simulations employing the shear stress transport $$k-\omega $$ k - ω turbulence model. The feature Dynamic Mesh Motion implements the motion of the oscillating flap. Regarding low-speed wind tunnel testing for a Reynolds number of $$0.5 \times 10^{6}$$ 0.5 × 10 6 the flap movement around a dropped hinge point, which is located outside the flap, offers benefits with regard to additional lift and delayed flow separation at the flap compared to a flap movement around a hinge point, which is located at 15 % of the flap chord length. Flow separation can be suppressed beyond the maximum static flap deflection angle. By means of an oscillating flap around the dropped hinge point, it is possible to reattach a separated flow at the flap and to keep it attached further on. For a Reynolds number of $$20 \times 10^6$$ 20 × 10 6 , reflecting full scale flight conditions, additional lift is generated for both rotational axis positions.


Author(s):  
Michael Leschziner ◽  
Ning Li ◽  
Fabrizio Tessicini

This paper provides a discussion of several aspects of the construction of approaches that combine statistical (Reynolds-averaged Navier–Stokes, RANS) models with large eddy simulation (LES), with the objective of making LES an economically viable method for predicting complex, high Reynolds number turbulent flows. The first part provides a review of alternative approaches, highlighting their rationale and major elements. Next, two particular methods are introduced in greater detail: one based on coupling near-wall RANS models to the outer LES domain on a single contiguous mesh, and the other involving the application of the RANS and LES procedures on separate zones, the former confined to a thin near-wall layer. Examples for their performance are included for channel flow and, in the case of the zonal strategy, for three separated flows. Finally, a discussion of prospects is given, as viewed from the writer's perspective.


Aerospace ◽  
2021 ◽  
Vol 8 (8) ◽  
pp. 216
Author(s):  
Emanuel A. R. Camacho ◽  
Fernando M. S. P. Neves ◽  
André R. R. Silva ◽  
Jorge M. M. Barata

Natural flight has consistently been the wellspring of many creative minds, yet recreating the propulsive systems of natural flyers is quite hard and challenging. Regarding propulsive systems design, biomimetics offers a wide variety of solutions that can be applied at low Reynolds numbers, achieving high performance and maneuverability systems. The main goal of the current work is to computationally investigate the thrust-power intricacies while operating at different Reynolds numbers, reduced frequencies, nondimensional amplitudes, and mean angles of attack of the oscillatory motion of a NACA0012 airfoil. Simulations are performed utilizing a RANS (Reynolds Averaged Navier-Stokes) approach for a Reynolds number between 8.5×103 and 3.4×104, reduced frequencies within 1 and 5, and Strouhal numbers from 0.1 to 0.4. The influence of the mean angle-of-attack is also studied in the range of 0∘ to 10∘. The outcomes show ideal operational conditions for the diverse Reynolds numbers, and results regarding thrust-power correlations and the influence of the mean angle-of-attack on the aerodynamic coefficients and the propulsive efficiency are widely explored.


1973 ◽  
Vol 59 (2) ◽  
pp. 391-396 ◽  
Author(s):  
N. C. Freeman ◽  
S. Kumar

It is shown that, for a spherically symmetric expansion of a gas into a low pressure, the shock wave with area change region discussed earlier (Freeman & Kumar 1972) can be further divided into two parts. For the Navier–Stokes equation, these are a region in which the asymptotic zero-pressure behaviour predicted by Ladyzhenskii is achieved followed further downstream by a transition to subsonic-type flow. The distance of this final region downstream is of order (pressure)−2/3 × (Reynolds number)−1/3.


Sign in / Sign up

Export Citation Format

Share Document