scholarly journals Direct numerical simulation of the oscillatory flow around a sphere resting on a rough bottom

2017 ◽  
Vol 822 ◽  
pp. 235-266 ◽  
Author(s):  
Marco Mazzuoli ◽  
Paolo Blondeaux ◽  
Julian Simeonov ◽  
Joseph Calantoni

The oscillatory flow around a spherical object lying on a rough bottom is investigated by means of direct numerical simulations of the continuity and Navier–Stokes equations. The rough bottom is simulated by a layer/multiple layers of spherical particles, the size of which is much smaller that the size of the object. The period and amplitude of the velocity oscillations of the free stream are chosen to mimic the flow at the bottom of sea waves and the size of the small spherical particles falls in the range of coarse sand/very fine gravel. Even though the computational costs allow only the simulation of moderate values of the Reynolds number characterizing the bottom boundary layer, the results show that the coherent vortex structures, shed by the spherical object, can break up and generate turbulence, if the Reynolds number of the object is sufficiently large. The knowledge of the velocity field allows the dynamics of the large-scale coherent vortices shed by the object to be determined and turbulence characteristics to be evaluated. Moreover, the forces and torques acting on both the large spherical object and the small particles, simulating sediment grains, can be determined and analysed, thus laying the groundwork for the investigation of sediment dynamics and scour developments.

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.


2014 ◽  
Vol 752 ◽  
pp. 602-625 ◽  
Author(s):  
Kengo Deguchi ◽  
Philip Hall

AbstractOur concern in this paper is with high-Reynolds-number nonlinear equilibrium solutions of the Navier–Stokes equations for boundary-layer flows. Here we consider the asymptotic suction boundary layer (ASBL) which we take as a prototype parallel boundary layer. Solutions of the equations of motion are obtained using a homotopy continuation from two known types of solutions for plane Couette flow. At high Reynolds numbers, it is shown that the first type of solution takes the form of a vortex–wave interaction (VWI) state, see Hall & Smith (J. Fluid Mech., vol. 227, 1991, pp. 641–666), and is located in the main part of the boundary layer. On the other hand, here the second type is found to support an equilibrium solution of the unit-Reynolds-number Navier–Stokes equations in a layer located a distance of $\def \xmlpi #1{}\def \mathsfbi #1{\boldsymbol {\mathsf {#1}}}\let \le =\leqslant \let \leq =\leqslant \let \ge =\geqslant \let \geq =\geqslant \def \Pr {\mathit {Pr}}\def \Fr {\mathit {Fr}}\def \Rey {\mathit {Re}}O(\ln \mathit{Re})$ from the wall. Here $\mathit{Re}$ is the Reynolds number based on the free-stream speed and the unperturbed boundary-layer thickness. The streaky field produced by the interaction grows exponentially below the layer and takes its maximum size within the unperturbed boundary layer. The results suggest the possibility of two distinct types of streaky coherent structures existing, possibly simultaneously, in disturbed boundary layers.


2021 ◽  
Vol 930 ◽  
Author(s):  
Kartik P. Iyer ◽  
Katepalli R. Sreenivasan ◽  
P.K. Yeung

Using direct numerical simulations performed on periodic cubes of various sizes, the largest being $8192^3$ , we examine the nonlinear advection term in the Navier–Stokes equations generating fully developed turbulence. We find significant dissipation even in flow regions where nonlinearity is locally absent. With increasing Reynolds number, the Navier–Stokes dynamics amplifies the nonlinearity in a global sense. This nonlinear amplification with increasing Reynolds number renders the vortex stretching mechanism more intermittent, with the global suppression of nonlinearity, reported previously, restricted to low Reynolds numbers. In regions where vortex stretching is absent, the angle and the ratio between the convective vorticity and solenoidal advection in three-dimensional isotropic turbulence are statistically similar to those in the two-dimensional case, despite the fundamental differences between them.


2021 ◽  
Vol 8 (3) ◽  
pp. 418-424
Author(s):  
Syed Fazuruddin ◽  
Seelam Sreekanth ◽  
G. Sankara Sekhar Raju

Incompressible 2-D Navier-stokes equations for various values of Reynolds number with and without partial slip conditions are studied numerically. The Lid-Driven cavity (LDC) with uniform driven lid problem is employed with vorticity - Stream function (VSF) approach. The uniform mesh grid is used in finite difference approximation for solving the governing Navier-stokes equations and developed MATLAB code. The numerical method is validated with benchmark results. The present work is focused on the analysis of lid driven cavity flow of incompressible fluid with partial slip conditions (imposed on side walls of the cavity). The fluid flow patterns are studied with wide range of Reynolds number and slip parameters.


1998 ◽  
Vol 120 (2) ◽  
pp. 257-262 ◽  
Author(s):  
Peter Gerlinger ◽  
Dieter Bru¨ggemann

A multigrid method for convergence acceleration is used for solving coupled fluid and turbulence transport equations. For turbulence closure a low-Reynolds-number q-ω turbulence model is employed, which requires very fine grids in the near wall regions. Due to the use of fine grids, convergence of most iterative solvers slows down, making the use of multigrid techniques especially attractive. However, special care has to be taken on the strong nonlinear turbulent source terms during restriction from fine to coarse grids. Due to the hyperbolic character of the governing equations in supersonic flows and the occurrence of shock waves, modifications to standard multigrid techniques are necessary. A simple and effective method is presented that enables the multigrid scheme to converge. A strong reduction in the required number of multigrid cycles and work units is achieved for different test cases, including a Mack 2 flow over a backward facing step.


2018 ◽  
Vol 856 ◽  
Author(s):  
M. Borgnino ◽  
G. Boffetta ◽  
F. De Lillo ◽  
M. Cencini

We study the dynamics and the statistics of dilute suspensions of gyrotactic swimmers, a model for many aquatic motile microorganisms. By means of extensive numerical simulations of the Navier–Stokes equations at different Reynolds numbers, we investigate preferential sampling and small-scale clustering as a function of the swimming (stability and speed) and shape parameters, considering in particular the limits of spherical and rod-like particles. While spherical swimmers preferentially sample local downwelling flow, for elongated swimmers we observe a transition from downwelling to upwelling regions at sufficiently high swimming speed. The spatial distribution of both spherical and elongated swimmers is found to be fractal at small scales in a wide range of swimming parameters. The direct comparison between the different shapes shows that spherical swimmers are more clusterized at small stability and speed numbers, while for large values of the parameters elongated cells concentrate more. The relevance of our results for phytoplankton swimming in the ocean is briefly discussed.


Author(s):  
Bakhtier Farouk ◽  
Murat K. Aktas

Formation of vortical flow structures in a rectangular enclosure due to acoustic streaming is investigated numerically. The oscillatory flow field in the enclosure is created by the vibration of a vertical side wall of the enclosure. The frequency of the wall vibration is chosen such that a standing wave forms in the enclosure. The interaction of this standing wave with the horizontal solid walls leads to the production of Rayleigh type acoustic streaming flow patterns in the enclosure. All four walls of the enclosure considered are thermally insulated. The fully compressible form of the Navier-Stokes equations is considered and an explicit time-marching algorithm is used to explicitly track the acoustic waves. Numerical solutions are obtained by employing a highly accurate flux corrected transport (FCT) algorithm for the convection terms. A time-splitting technique is used to couple the viscous and diffusion terms of the full Navier-Stokes equations. Non-uniform grid structure is employed in the computations. The simulation of the primary oscillatory flow and the secondary (steady) streaming flows in the enclosure is performed. Streaming flow patterns are obtained by time averaging the primary oscillatory flow velocity distributions. The effect of the amount of wall displacement on the formation of the oscillatory flow field and the streaming structures are studied. Computations indicate that the nonlinearity of the acoustic field increases with increasing amount of the vibration amplitude. The form and the strength of the secondary flow associated with the oscillatory flow field and viscous effects are found to be strongly correlated to the maximum displacement of the vibrating wall. Total number of acoustic streaming cells per wavelength is also determined by the strength and the level of the nonlinearity of the sound field in the resonator.


Sign in / Sign up

Export Citation Format

Share Document