Vortex merger in rotating stratified flows

2002 ◽  
Vol 455 ◽  
pp. 83-101 ◽  
Author(s):  
DAVID G. DRITSCHEL

This paper describes the interaction of symmetric vortices in a three-dimensional quasi-geostrophic fluid. The initial vortices are taken to be uniform-potential-vorticity ellipsoids, of height 2h and width 2R, and with centres at (±d/2; 0, 0), embedded within a background flow having constant background rotational and buoyancy frequencies, f/2 and N respectively. This problem was previously studied by von Hardenburg et al. (2000), who determined the dimensionless critical merger distance d/R as a function of the height-to-width aspect ratio h/R (scaled by f/N). Their study, however, was limited to small to moderate values of h/R, as it was anticipated that merger at large h/R would reduce to that for two columnar two-dimensional vortices, i.e. d/R ≈ 3.31. Here, it is shown that no such two-dimensional limit exists; merger is found to occur at any aspect ratio, with d ∼ h for h/R [Gt ] 1.New results are also found for small to moderate values of h/R. In particular, our numerical simulations reveal that asymmetric merger is predominant, despite the initial conditions, if one includes a small amount of random noise. For small to moderate h/R, decreasing the initial separation distance d first results in a weak exchange of material, with one vortex growing at the expense of the other. As d decreases further, this exchange increases and leads to two dominant but strongly asymmetric vortices. Finally, for yet smaller d, rapid merger into a single dominant vortex occurs – in effect the initial vortices exchange nearly all of their material with one another in a nearly symmetrical fashion.

1990 ◽  
Vol 216 ◽  
pp. 71-91 ◽  
Author(s):  
Bryan Travis ◽  
Peter Olson ◽  
Gerald Schubert

The stability of two-dimensional thermal convection in an infinite-Prandtl-number fluid layer with zero-stress boundaries is investigated using numerical calculations in three-dimensional rectangles. At low Rayleigh numbers (Ra < 20000) calculations of the stability of two-dimensional rolls to cross-roll disturbances are in agreement with the predictions of Bolton & Busse for a fluid with a large but finite Prandtl number. Within the range 2 × 104 < Ra [les ] 5 × 105, steady rolls with basic wavenumber α > 2.22 (aspect ratio < 1.41) are stable solutions. Two-dimensional rolls with basic wavenumber α < 1.96 (aspect ratio > 1.6) are time dependent for Ra > 4 × 104. For every case in which the initial condition was a time-dependent large-aspect-ratio roll, two-dimensional convection was found to be unstable to three-dimensional convection. Time-dependent rolls are replaced by either bimodal or knot convection in cases where the horizontal dimensions of the rectangular box are less than twice the depth. The bimodal planforms are steady states for Ra [les ] 105, but one case at Ra = 5 × 105 exhibits time dependence in the form of pulsating knots. Calculations at Ra = 105 in larger domains resulted in fully three-dimensional cellular planforms. A steady-state square planform was obtained in a 2.4 × 2.4 × 1 rectangular box. started from random initial conditions. Calculations in a 3 × 3 × 1 box produced steady hexagonal cells when started from random initial conditions, and a rectangular planform when started from a two-dimensional roll. An hexagonal planform started in a 3.5 × 3.5 × 1 box at Ra = 105 exhibited oscillatory time dependence, including boundary-layer instabilities and pulsating plumes. Thus, the stable planform in three-dimensional convection is sensitive to the size of the rectangular domain and the initial conditions. The sensitivity of heat transfer to planform variations is less than 10%.


2020 ◽  
Vol 25 (5) ◽  
pp. 1101-1123 ◽  
Author(s):  
Lidong Fang ◽  
Apala Majumdar ◽  
Lei Zhang

We study nematic equilibria on rectangular domains, in a reduced two-dimensional Landau–de Gennes framework. These reduced equilibria carry over to the three-dimensional framework at a special temperature. There is one essential model variable, [Formula: see text], which is a geometry-dependent and material-dependent variable. We compute the limiting profiles exactly in two distinguished limits: the [Formula: see text] 0 limit relevant for macroscopic domains and the [Formula: see text] limit relevant for nanoscale domains. The limiting profile has line defects near the shorter edges in the [Formula: see text] limit, whereas we observe fractional point defects in the [Formula: see text] 0 limit. The analytical studies are complemented by some bifurcation diagrams for these reduced equilibria as a function of [Formula: see text] and the rectangular aspect ratio. We also introduce the concept of ‘non-trivial’ topologies and study the relaxation of non-trivial topologies to trivial topologies mediated via point and line defects, with potential consequences for non-equilibrium phenomena and switching dynamics.


2015 ◽  
Vol 61 (228) ◽  
pp. 702-712 ◽  
Author(s):  
Tong Zhang ◽  
Lili Ju ◽  
Wei Leng ◽  
Stephen Price ◽  
Max Gunzburger

AbstractFor many regions, glacier inaccessibility results in sparse geometric datasets for use as model initial conditions (e.g. along the central flowline only). In these cases, two-dimensional (2-D) flowline models are often used to study glacier dynamics. Here we systematically investigate the applicability of a 2-D, first-order Stokes approximation flowline model (FLM), modified by shape factors, for the simulation of land-terminating glaciers by comparing it with a 3-D, ‘full’-Stokes ice-flow model (FSM). Based on steady-state and transient, thermomechanically uncoupled and coupled computational experiments, we explore the sensitivities of the FLM and FSM to ice geometry, temperature and forward model integration time. We find that, compared to the FSM, the FLM generally produces slower horizontal velocities, due to simplifications inherent to the FLM and to the underestimation of the shape factor. For polythermal glaciers, those with temperate ice zones, or when basal sliding is important, we find significant differences between simulation results when using the FLM versus the FSM. Over time, initially small differences between the FLM and FSM become much larger, particularly near cold/temperate ice transition surfaces. Long time integrations further increase small initial differences between the two models. We conclude that the FLM should be applied with caution when modelling glacier changes under a warming climate or over long periods of time.


1997 ◽  
Vol 335 ◽  
pp. 1-28 ◽  
Author(s):  
MELVIN E. STERN ◽  
ERIC P. CHASSIGNET ◽  
J. A. WHITEHEAD

The previously observed spatial evolution of the two-dimensional turbulent flow from a source on the vertical wall of a shallow layer of rapidly rotating fluid is strikingly different from the non-rotating three-dimensional counterpart, insofar as the instability eddies generated in the former case cause the flow to separate completely from the wall at a finite downstream distance. In seeking an explanation of this, we first compute the temporal evolution of two-dimensional finite-amplitude waves on an unstable laminar jet using a finite difference calculation at large Reynolds number. This yields a dipolar vorticity pattern which propagates normal to the wall, while leaving some of the near-wall vorticity (negative) of the basic flow behind. The residual far-field eddy therefore contains a net positive circulation and this property is incorporated in a heuristic point-vortex model of the spatial evolution of the instability eddies observed in a laboratory experiment of a flow emerging from a source on a vertical wall in a rotating tank. The model parameterizes the effect of Ekman bottom friction in decreasing the circulation of eddies which are periodically emitted from the source flow on the wall. Further downstream, the point vortices of the model merge and separate abruptly from the wall; the statistics suggest that the downstream separation distance scales with the Ekman spin-up time (inversely proportional to the square root of the Coriolis parameter f) and with the mean source velocity. When the latter is small and f is large, qualitative support is obtained from laboratory experiments.


Author(s):  
Abhishek Agrawal ◽  
Amit Agrawal

Three-dimensional lattice Boltzmann method based simulations of a microduct have been undertaken in this paper. The objective is to understand the different physical phenomena occurring at these small scales and to investigate when the flow can be treated as two-dimensional. Towards this end, the Knudsen number and aspect ratio (depth to width ratio) are varied for a fixed pressure ratio. The pressure in the microduct is non-linear with the non-linearity in pressure reducing with an increase in Knudsen number. The pressure and velocity behaves somewhat similar to two-dimensional microchannels even when the aspect ratio is unity. The slip velocity at the impenetrable wall has two components: along and perpendicular to the flow. Our results show that the streamwise velocity near the centerline is relatively invariant along the depth for aspect ratio more than three, suggesting that the microduct can be modeled as a two-dimensional microchannel. However, the velocity component along the depth is never identically zero, implying that the flow is not truly two-dimensional. A curious change in vector direction in a plane normal to the flow direction is observed around aspect ratio of four. These first set of three-dimensional results are significant because they will help in theoretical development and flow modeling at micro scales.


1990 ◽  
Vol 45 (9-10) ◽  
pp. 1059-1073 ◽  
Author(s):  
G. Knorr ◽  
J. P. Lynov ◽  
H. L. Pécseli

Abstract The three-dimensional incompressible Euler equations are expanded in eigenflows of the curl operator, which represent positive and negative helicity flows in a particularly simple and convenient way. Four different basic types of interactions between eigenflows are found. Two represent an "inverse cascade", the interaction familiar from the two-dimensional Euler equations, in which only modes of the same sign of the helicity interact. The other two interactions mix positive and negative helicity modes. Only these interactions can transport all of the available energy to higher wave numbers. Initial conditions, which lead to the appearance of structures and self-organization, are discussed.


2014 ◽  
Vol 58 (01) ◽  
pp. 1-19
Author(s):  
Michael J. Hughes ◽  
Young T. Shen

The behavior of the force on a rudder changes significantly after the inception of stall, requiring different mathematical formulae to compute rudder forces prior-and poststall. Determining the inception angle at which stall occurs is important for predicting the rudder force on a maneuvering ship. A method to compute the inception angle of stall on a rudder is presented in this article. The theoretical formulation is based on a flow similarity approach, which relates three-dimensional rudder stall inception with two-dimensional airfoil data. Rudders are low-aspect ratio wings, and the three-dimensional lift is based on the low-aspect ratio wing theory. The two-dimensional airfoil stall data are obtained from National Advisory Committee for Aeronautics (NACA) reports. The derived theory is first validated with wind tunnel data from foils with a NACA 0015 profile of aspect ratios 1, 2, and 3. The theory is also validated with data from foils with a NACA 0012 profile and an aspect ratio of 2, 3, and 6.


2010 ◽  
Vol 658 ◽  
pp. 492-508 ◽  
Author(s):  
ZIV KIZNER ◽  
RUVIM KHVOLES ◽  
DAVID A. KESSLER

A theory of viscous evolution and selection of symmetric two-dimensional dipoles is suggested, based on a combination of numerical simulations and an asymptotic analysis, where the slow time scale associated with the vorticity diffusion due to viscosity is incorporated. It is shown that viscosity first brings a dipole to an intermediate asymptotic state, which is independent of the initial conditions, and then slowly takes the dipole away from this state. We demonstrate that, among the variety of possible ideal-fluid dipole solutions, viscosity going to zero selects a unique solution, which is described to high accuracy by the elliptical dipole solution with a separatrix aspect ratio of 1.037.


Sign in / Sign up

Export Citation Format

Share Document