Formation of columnar baroclinic vortices in thermally stratified nonlinear spin-up

2012 ◽  
Vol 702 ◽  
pp. 265-285 ◽  
Author(s):  
J. R. Pacheco ◽  
R. Verzicco

AbstractWe investigate the mechanisms that affect the formation of columnar vortices for spin-up in a cylinder where the temperatures at the horizontal walls are prescribed. Numerical results from the three-dimensional Navier–Stokes equations show that a short-lived instability, suppressed by the combined effect of rotation and stratification, generates small temperature variations in the azimuthal direction. Temperature-gradient anomalies produce vorticity, and these vortices stir the fluid at the interface of the central vortex core thus reinforcing the temperature gradients. For sufficiently strong temperature gradients, the central vortex core breaks up into several columnar vortices. It is found, in particular, that small aspect ratios (height over radius of the cylindrical fluid layer) $\Gamma = 1, 2$ tend to inhibit the instability, while larger ones, $\Gamma = 3. 3$, have the opposite effect. The main source of instability is the baroclinic vorticity production and not the presence of a solid sidewall since, counter-intuitively, the flow is more unstable for a free-slip boundary than for a no-slip one. Finally the effect of the temperature boundary conditions (isothermal versus adiabatic) on the horizontal boundaries has been investigated. The adiabatic boundaries help to preserve for longer times the sloping density interfaces that feed, with their potential energy, the baroclinic vorticity production; this results in more unstable flows.

2014 ◽  
Vol 740 ◽  
Author(s):  
Mohamed Ali ◽  
Malek Abid

AbstractWe report a self-similar behaviour of solutions (obtained numerically) of the Navier–Stokes equations behind a single-blade rotor. That is, the helical vortex core in the wake of a rotating blade is self-similar as a function of its age. Profiles of vorticity and azimuthal velocity in the vortex core are characterized, their similarity variables are identified and scaling laws of these variables are given. Solutions of incompressible three-dimensional Navier–Stokes equations for Reynolds numbers up to $Re= 2000$ are considered.


1998 ◽  
Vol 185 ◽  
pp. 217-218
Author(s):  
Kwing L. Chan ◽  
Y.C. Kim

Currently, the most successful direct simulation of the solar granules (and the convection/radiation transition layer) is the three-dimensional (3D) model computed by Stein and Nordlund (1989). So far, there is no other similar 3D models available for comparison [however, see Ludwig et al. (1997) for a recent 2D calculation]. We are developing an alternative numerical approach to simulate the 3D radiation hydrodynamics of this layer. In this approach, the Eddington approximation is used to handle the radiation rather than solving the radiative transfer equations along rays, and the ADISM method (Chan and Wolff 1982) which solves the Navier Stokes equations in conservative forms is used to speed up the thermal relaxation of the fluid layer. We are in the process of testing the numerical accuracy of the codes. This paper summarizes the results of a test that illustrate the effects of vertical space resolution on the mean profiles of some important quantities.


1986 ◽  
Vol 108 (1) ◽  
pp. 93-102 ◽  
Author(s):  
B. C. Weinberg ◽  
R.-J. Yang ◽  
H. McDonald ◽  
S. J. Shamroth

The multidimensional, ensemble-averaged, compressible, time-dependent Navier-Stokes equations have been used to study the turbulent flow field in two and three-dimensional turbine cascades. The viscous regions of the flow were resolved and no-slip boundary conditions were utilized on solid surfaces. The calculations were performed in a constructive ‘O’-type grid which allows representation of the blade rounded trailing edge. Converged solutions were obtained in relatively few time steps (∼ 80–150) and comparisons for both surface pressure and heat transfer showed good agreement with data. The three-dimensional turbine cascade calculation showed many of the expected flow-field features.


1977 ◽  
Vol 83 (1) ◽  
pp. 1-31 ◽  
Author(s):  
G. D. Mallinson ◽  
G. De Vahl Davis

The solution of the steady-state Navier–Stokes equations in three dimensions has been obtained by a numerical method for the problem of natural convection in a rectangular cavity as a result of differential side heating. In the past, this problem has generally been treated as though it were two-dimensional. The solutions explore the three-dimensional motion generated by the presence of no-slip adiabatic end walls. For Ra = 104, the three-dimensional motion is shown to be the result of the inertial interaction of the rotating flow with the stationary walls together with a contribution arising from buoyancy forces generated by longitudinal temperature gradients. The inertial effect is inversely dependent on the Prandtl number, whereas the thermal effect is nearly constant. For higher values of Ra, multiple longitudinal flows develop which are a delicate function of Ra, Pr and the cavity aspect ratios.


1996 ◽  
Vol 63 (2) ◽  
pp. 467-473
Author(s):  
F. Y. Huang ◽  
C. D. Mote

A new model of the thin viscous fluid film, constrained between two translating, flexible surfaces, is presented in this paper: The unsteady inertia of the film is included in the model. The derivation starts with the reduced three-dimensional Navier-Stokes equations for an incompressible viscous fluid with a small Reynolds number. By introduction of an approximate velocity field, which satisfies the continuity equation and the no-slip boundary conditions exactly, into weighted integrals of the three-dimensional equations over the film thickness, a two-dimensional thin film equation is obtained explicitly in a closed form. The 1th thin film equation is obtained when the velocity field is approximated by 21th order polynominals, and the three-dimensional viscous film is described with increasing accuracy by thin film equations of increasing order. Two cases are used to illustrate the coupling of the film to the vibration of the structure and to show that the second thin film equation can be applied successfully to the prediction of a coupled film-structure response in the range of most applications. A reduced thin film equation is derived through approximation of the second thin film equation that relates the film pressure to transverse accelerations and velocities, and to slopes and slope rates of the two translating surfaces.


2017 ◽  
Vol 139 (8) ◽  
Author(s):  
K. Srinivasan ◽  
P. M. V. Subbarao ◽  
S. R. Kale

The present work investigates the extension of Navier–Stokes equations from slip-to-transition regimes with higher-order slip boundary condition. To achieve this, a slip model based on the second-order slip boundary condition was derived and a special procedure was developed to simulate slip models using FLUENT®. The boundary profile for both top and bottom walls was solved for each pressure ratio by the customized user-defined function and then passed to the FLUENT® solver. The flow characteristics in microchannels of various aspect ratios (a = H/W = 0.002, 0.01, and 0.1) by generating accurate and high-resolution experimental data along with the computational validation was studied. For that, microchannel system was fabricated in silicon wafers with controlled surface structure and each system has several identical microchannels of same dimensions in parallel and the processed wafer was bonded with a plane wafer. The increased flow rate reduced uncertainty substantially. The experiments were performed up to maximum outlet Knudsen number of 1.01 with nitrogen and the second-order slip coefficients were found to be C1 = 1.119–1.288 (TMAC = 0.944–0.874) and C2 = 0.34.


2017 ◽  
Vol 828 ◽  
pp. 837-866 ◽  
Author(s):  
Jeffrey Tithof ◽  
Balachandra Suri ◽  
Ravi Kumar Pallantla ◽  
Roman O. Grigoriev ◽  
Michael F. Schatz

We present a combined experimental and theoretical study of the primary and secondary instabilities in a Kolmogorov-like flow. The experiment uses electromagnetic forcing with an approximately sinusoidal spatial profile to drive a quasi-two-dimensional (Q2D) shear flow in a thin layer of electrolyte suspended on a thin lubricating layer of a dielectric fluid. Theoretical analysis is based on a two-dimensional (2D) model (Suri et al., Phys. Fluids, vol. 26 (5), 2014, 053601), derived from first principles by depth-averaging the full three-dimensional Navier–Stokes equations. As the strength of the forcing is increased, the Q2D flow in the experiment undergoes a series of bifurcations, which is compared with results from direct numerical simulations of the 2D model. The effects of confinement and the forcing profile are studied by performing simulations that assume spatial periodicity and strictly sinusoidal forcing, as well as simulations with realistic no-slip boundary conditions and an experimentally validated forcing profile. We find that only the simulation subject to physical no-slip boundary conditions and a realistic forcing profile provides close, quantitative agreement with the experiment. Our analysis offers additional validation of the 2D model as well as a demonstration of the importance of properly modelling the forcing and boundary conditions.


Author(s):  
S. Itoh ◽  
A. Tani

The initial-boundary value problem for the non-homogeneous Navier-Stokes equations including the slipping on the solid boundary is considered. The unique solvability is established locally in time for the three-dimensional problem and globally in time for the two-dimensional problem without so-called smallness restrictions.


Sign in / Sign up

Export Citation Format

Share Document