Waves and turbulence on a beta-plane

1975 ◽  
Vol 69 (3) ◽  
pp. 417-443 ◽  
Author(s):  
Peter B. Rhines

Two-dimensional eddies in a homogeneous fluid at large Reynolds number, if closely packed, are known to evolve towards larger scales. In the presence of a restoring force, the geophysical beta-effect, this cascade produces a field of waves without loss of energy, and the turbulent migration of the dominant scale nearly ceases at a wavenumber kβ = (β/2U)½ independent of the initial conditions other than U, the r.m.s. particle speed, and β, the northward gradient of the Coriolis frequency.The conversion of turbulence into waves yields, in addition, more narrowly peaked wavenumber spectra and less fine-structure in the spatial maps, while smoothly distributing the energy about physical space.The theory is discussed, using known integral constraints and similarity solutions, model equations, weak-interaction wave theory (which provides the terminus for the cascade) and other linearized instability theory. Computer experiments with both finite-difference and spectral codes are reported. The central quantity is the cascade rate, defined as \[ T = 2\int_0^{\infty} kF(k)dk/U^3\langle k\rangle , \] where F is the nonlinear transfer spectrum and 〈k〉 the mean wavenumber of the energy spectrum. (In unforced inviscid flow T is simply U−1d〈k〉−1/dt, or the rate at which the dominant scale expands in time t.) T is shown to have a mean value of 3·0 × 10−2 for pure two-dimensional turbulence, but this decreases by a factor of five at the transition to wave motion. We infer from weak-interaction theory even smaller values for k [Lt ] kβ.After passing through a state of propagating waves, the homogeneous cascade tends towards a flow of alternating zonal jets which, we suggest, are almost perfectly steady. When the energy is intermittent in space, however, model equations show that the cascade is halted simply by the spreading of energy about space, and then the end state of a zonal flow is probably not achieved.The geophysical application is that the cascade of pure turbulence to large scales is defeated by wave propagation, helping to explain why the energy-containing eddies in the ocean and atmosphere, though significantly nonlinear, fail to reach the size of their respective domains, and are much smaller. For typical ocean flows, $k_{\beta}^{-1} = 70\,{\rm km} $, while for the atmosphere, $k_{\beta}^{-1} = 1000\,{\rm km}$. In addition the cascade generates, by itself, zonal flow (or more generally, flow along geostrophic contours).

1984 ◽  
Vol 146 ◽  
pp. 21-43 ◽  
Author(s):  
James C. Mcwilliams

A study is made of some numerical calculations of two-dimensional and geostrophic turbulent flows. The primary result is that, under a broad range of circumstances, the flow structure has its vorticity concentrated in a small fraction of the spatial domain, and these concentrations typically have lifetimes long compared with the characteristic time for nonlinear interactions in turbulent flow (i.e. an eddy turnaround time). When such vorticity concentrations occur, they tend to assume an axisymmetric shape and persist under passive advection by the large-scale flow, except for relatively rare encounters with other centres of concentration. These structures can arise from random initial conditions without vorticity concentration, evolving in the midst of what has been traditionally characterized as the ‘cascade’ of isotropic, homogeneous, large-Reynolds-number turbulence: the systematic elongation of isolines of vorticity associated with the transfer of vorticity to smaller scales, eventually to dissipation scales, and the transfer of energy to larger scales. When the vorticity concentrations are a sufficiently dominant component of the total vorticity field, the cascade processes are suppressed. The demonstration of persistent vorticity concentrations on intermediate scales - smaller than the scale of the peak of the energy spectrum and larger than the dissipation scales - does not invalidate many of the traditional characterizations of two-dimensional and geostrophic turbulence, but I believe it shows them to be substantially incomplete with respect to a fundamental phenomenon in such flows.


2008 ◽  
Vol 602 ◽  
pp. 303-326 ◽  
Author(s):  
E. PLAUT ◽  
Y. LEBRANCHU ◽  
R. SIMITEV ◽  
F. H. BUSSE

A general reformulation of the Reynolds stresses created by two-dimensional waves breaking a translational or a rotational invariance is described. This reformulation emphasizes the importance of a geometrical factor: the slope of the separatrices of the wave flow. Its physical relevance is illustrated by two model systems: waves destabilizing open shear flows; and thermal Rossby waves in spherical shell convection with rotation. In the case of shear-flow waves, a new expression of the Reynolds–Orr amplification mechanism is obtained, and a good understanding of the form of the mean pressure and velocity fields created by weakly nonlinear waves is gained. In the case of thermal Rossby waves, results of a three-dimensional code using no-slip boundary conditions are presented in the nonlinear regime, and compared with those of a two-dimensional quasi-geostrophic model. A semi-quantitative agreement is obtained on the flow amplitudes, but discrepancies are observed concerning the nonlinear frequency shifts. With the quasi-geostrophic model we also revisit a geometrical formula proposed by Zhang to interpret the form of the zonal flow created by the waves, and explore the very low Ekman-number regime. A change in the nature of the wave bifurcation, from supercritical to subcritical, is found.


1969 ◽  
Vol 91 (3) ◽  
pp. 353-358 ◽  
Author(s):  
W. A. Gustafson ◽  
I. Pelech

The two-dimensional, incompressible laminar boundary layer on a strongly curved wall in a converging channel is investigated for the special case of potential velocity inversely proportional to the distance along the wall. Similarity solutions of the momentum equation are obtained by two different methods and the differences between the methods are discussed. The numerical results show that displacement and momentum thickness increase linearly with curvature while skin friction decreases linearly.


2021 ◽  
Vol 932 ◽  
Author(s):  
Pierre Ricco ◽  
Claudia Alvarenga

The development and growth of unsteady three-dimensional vortical disturbances entrained in the entry region of a circular pipe is investigated by asymptotic and numerical methods for Reynolds numbers between $1000$ and $10\,000$ , based on the pipe radius and the bulk velocity. Near the pipe mouth, composite asymptotic solutions describe the dynamics of the oncoming disturbances, revealing how these disturbances are altered by the viscous layer attached to the pipe wall. The perturbation velocity profiles near the pipe mouth are employed as rigorous initial conditions for the boundary-region equations, which describe the flow in the limit of low frequency and large Reynolds number. The disturbance flow is initially primarily present within the base-flow boundary layer in the form of streamwise-elongated vortical structures, i.e. the streamwise velocity component displays an intense algebraic growth, while the cross-flow velocity components decay. Farther downstream the disturbance flow occupies the whole pipe, although the base flow is mostly inviscid in the core. The transient growth and subsequent viscous decay are confined in the entrance region, i.e. where the base flow has not reached the fully developed Poiseuille profile. Increasing the Reynolds number and decreasing the frequency causes more intense perturbations, whereas small azimuthal wavelengths and radial characteristic length scales intensify the viscous dissipation of the disturbance. The azimuthal wavelength that causes the maximum growth is found. The velocity profiles are compared successfully with available experimental data and the theoretical results are helpful to interpret the only direct numerical dataset of a disturbed pipe-entry flow.


2013 ◽  
Vol 2013 ◽  
pp. 1-13
Author(s):  
Motoaki Saruwatari ◽  
Masa-aki Hashimoto ◽  
Ryohei Fukuda ◽  
Shin-ichiro Fujimoto

We investigate the r-process nucleosynthesis during the magnetohydrodynamical (MHD) explosion of a supernova in a helium star of 3.3 M⊙, where effects of neutrinos are taken into account using the leakage scheme in the two-dimensional (2D) hydrodynamic code. Jet-like explosion due to the combined effects of differential rotation and magnetic field is able to erode the lower electron fraction matter from the inner layers. We find that the ejected material of low electron fraction responsible for the r-process comes out from just outside the neutrino sphere deep inside the Fe-core. It is found that heavy element nucleosynthesis depends on the initial conditions of rotational and magnetic fields. In particular, the third peak of the distribution is significantly overproduced relative to the solar system abundances, which would indicate a possible r-process site owing to MHD jets in supernovae.


2006 ◽  
Vol Volume 5, Special Issue TAM... ◽  
Author(s):  
Mario Lefebvre

International audience A two-dimensional controlled stochastic process defined by a set of stochastic differential equations is considered. Contrary to the most frequent formulation, the control variables appear only in the infinitesimal variances of the process, rather than in the infinitesimal means. The differential game ends the first time the two controlled processes are equal or their difference is equal to a given constant. Explicit solutions to particular problems are obtained by making use of the method of similarity solutions to solve the appropriate partial differential equation. On considère un processus stochastique commandé bidimensionnel défini par un ensemble d'équations différentielles stochastiques. Contrairement à la formulation la plus fréquente, les variables de commande apparaissent dans les variances infinitésimales du processus, plutôt que dans les moyennes infinitésimales. Le jeu différentiel prend fin lorsque les deux processus sont égaux ou que leur différence est égale à une constante donnée. Des solutions explicites à des problèmes particuliers sont obtenues en utilisant la méthode des similitudes pour résoudre l'équation aux dérivées partielles appropriée.


2009 ◽  
Vol 2009 ◽  
pp. 1-22 ◽  
Author(s):  
Edson D. Leonel

A phase transition from integrability to nonintegrability in two-dimensional Hamiltonian mappings is described and characterized in terms of scaling arguments. The mappings considered produce a mixed structure in the phase space in the sense that, depending on the combination of the control parameters and initial conditions, KAM islands which are surrounded by chaotic seas that are limited by invariant tori are observed. Some dynamical properties for the largest component of the chaotic sea are obtained and described in terms of the control parameters. The average value and the deviation of the average value for chaotic components of a dynamical variable are described in terms of scaling laws, therefore critical exponents characterizing a scaling function that describes a phase transition are obtained and then classes of universality are characterized. The three models considered are: The Fermi-Ulam accelerator model, a periodically corrugate waveguide, and variant of the standard nontwist map.


Author(s):  
D.Y. Ivanov ◽  

Here we consider the initial-boundary value problems in a homogeneous cylindrical domain YI Ω ×+ ( Ω+ is an open two-dimensional bounded simply connected domain with a boundary 5 ∂Ω ∈C , 2 \ Ω≡ Ω − + R is the open exterior of the domain Ω+ , [0, ] YI ≡ Y is the height of the cylinder) on a time interval [0, ] TI ≡ T . The initial conditions and the boundary conditions on the bases of the cylinder are zero, and the boundary conditions on the lateral surface of the cylinder are given by the function 1 2 wx x yt ( , , ,) ( 1 2 (, ) x x ∈∂Ω , Y y ∈ I , T t I ∈ ). An approximate solution of such problems is obtained through the combined use of the Fourier method and the collocation boundary element method based on piecewise quadratic interpolation (PQI). The solution to the problem in the cylinder is expanded in a Fourier series in terms of eigenfunctions of the operator 2 By yy ≡ ∂ with the corresponding zero boundary conditions. The coefficients of such a Fourier series are solutions of problems for two-dimensional heat equations 2 2 t ∇ =∂ + u u ku . With a low smoothness of the functions w in the variable y, the weight of solutions at large values of k increases and the accuracy of solving the problem in the cylinder decreases. To maintain accuracy on a uniform grid, the step of discretization of the boundary function w with respect to the variable y is decreased by a factor of j. Here j is an averaged value of the quantity Y k π depending on the function w. In addition, the steps of discretization of functions ( ) 2 exp − τ k with respect to the variable τ in domains τ≤πT k are reduced by a factor of 2 2 k π . The steps in the remaining ranges of values τ and the steps by the other variables remain unchanged. The approximate solutions obtained on the basis of this procedure converge stably to exact solutions in the 2 ( ) LI I Y T × -norm with a cubic velocity uniformly with respect to sets of functions w, bounded by norm of functions with low smoothness in the variable y, uniformly along the length of the generatrix of the cylinder Y , and uniformly in the domain Ω . The latter is also associated with the use of PQI along the curve ∂Ω over the variable 2 2 ρ≡ − r d , which is carried out at small values of r ( d and r are the distances from the observed point of the domain Ω to the boundary ∂Ω and to the current point of integration along ∂Ω , respectively). The theoretical conclusions are confirmed by the results of the numerical solution of the problem in a circular cylinder, where the dependence of the boundary functions w on y is given by the normalized eigenfunctions of the differential operator By which vary in a sufficiently large range of values of k .


2021 ◽  
Vol 31 (15) ◽  
Author(s):  
Penghe Ge ◽  
Hongjun Cao

The existence of chaos in the Rulkov neuron model is proved based on Marotto’s theorem. Firstly, the stability conditions of the model are briefly renewed through analyzing the eigenvalues of the model, which are very important preconditions for the existence of a snap-back repeller. Secondly, the Rulkov neuron model is decomposed to a one-dimensional fast subsystem and a one-dimensional slow subsystem by the fast–slow dynamics technique, in which the fast subsystem has sensitive dependence on the initial conditions and its snap-back repeller and chaos can be verified by numerical methods, such as waveforms, Lyapunov exponents, and bifurcation diagrams. Thirdly, for the two-dimensional Rulkov neuron model, it is proved that there exists a snap-back repeller under two iterations by illustrating the existence of an intersection of three surfaces, which pave a new way to identify the existence of a snap-back repeller.


Sign in / Sign up

Export Citation Format

Share Document