Stochastically driven instability in rotating shear flows

The origin of hydrodynamic turbulence in rotating shear flows is investigated, with particular emphasis on the flows whose angular velocity decreases but whose specific angular momentum increases with the increasing radial coordinate. Such flows are Rayleigh stable, but must be turbulent in order to explain the observed data. Such a mismatch between the linear theory and the observations/experiments is more severe when any hydromagnetic/magnetohydrodynamic instability and then the corresponding turbulence therein is ruled out. This work explores the effect of stochastic noise on such hydrodynamic flows. We essentially concentrate on a small section of such a flow, which is nothing but a plane shear flow supplemented by the Coriolis effect. This also mimics a small section of an astrophysical accretion disc. It is found that such stochastically driven flows exhibit large temporal and spatial correlations of perturbation velocities and hence large energy dissipations of perturbation, which presumably generate the instability. A range of angular velocity (Omega) profiles of the background flow, starting from that of a constant specific angular momentum (lambda = Omega r(2); r being the radial coordinate) to a constant circular velocity (v(phi) = Omega r), is explored. However, all the background angular velocities exhibit identical growth and roughness exponents of their perturbations, revealing a unique universality class for the stochastically forced hydrodynamics of rotating shear flows. This work, to the best of our knowledge, is the first attempt to understand the origin of instability and turbulence in three-dimensional Rayleigh stable rotating shear flows by introducing additive noise to the underlying linearized governing equations. This has important implications to resolve the turbulence problem in astrophysical hydrodynamic flows such as accretion discs.

Magnetohydrodynamic stability of stochastically driven accretion flows

We investigate the evolution of magnetohydrodynamic (or hydromagnetic as coined by Chandrasekhar) perturbations in the presence of stochastic noise in rotating shear flows. The particular emphasis is the flows whose angular velocity decreases but specific angular momentum increases with increasing radial coordinate. Such flows, however, are Rayleigh stable but must be turbulent in order to explain astrophysical observed data and, hence, reveal a mismatch between the linear theory and observations and experiments. The mismatch seems to have been resolved, at least in certain regimes, in the presence of a weak magnetic field, revealing magnetorotational instability. The present work explores the effects of stochastic noise on such magnetohydrodynamic flows, in order to resolve the above mismatch generically for the hot flows. We essentially concentrate on a small section of such a flow which is nothing but a plane shear flow supplemented by the Coriolis effect, mimicking a small section of an astrophysical accretion disk around a compact object. It is found that such stochastically driven flows exhibit large temporal and spatial autocorrelations and cross-correlations of perturbation and, hence, large energy dissipations of perturbation, which generate instability. Interestingly, autocorrelations and cross-correlations appear independent of background angular velocity profiles, which are Rayleigh stable, indicating their universality. This work initiates our attempt to understand the evolution of three-dimensional hydromagnetic perturbations in rotating shear flows in the presence of stochastic noise.

A DYNAMO MODEL OF MAGNETIC ACTIVITY IN SOLAR-LIKE STARS WITH DIFFERENT ROTATIONAL VELOCITIES

We attempt to provide a quantitative theoretical explanation for the observations that Ca II H/K emission and X-ray emission from solar-like stars increase with decreasing Rossby number (i.e., with faster rotation). Assuming that these emissions are caused by magnetic cycles similar to the sunspot cycle, we construct flux transport dynamo models of 1M(circle dot) stars rotating with different rotation periods. We first compute the differential rotation and the meridional circulation inside these stars from a mean-field hydrodynamics model. Then these are substituted in our dynamo code to produce periodic solutions. We find that the dimensionless amplitude f(m) of the toroidal flux through the star increases with decreasing rotation period. The observational data can be matched if we assume the emissions to go as the power 3-4 of f(m). Assuming that the Babcock-Leighton mechanism saturates with increasing rotation, we can provide an explanation for the observed saturation of emission at low Rossby numbers. The main failure of our model is that it predicts an increase of the magnetic cycle period with increasing rotation rate, which is the opposite of what is found observationally. Much of our calculations are based on the assumption that the magnetic buoyancy makes the magnetic flux tubes rise radially from the bottom of the convection zone. Taking into account the fact that the Coriolis force diverts the magnetic flux tubes to rise parallel to the rotation axis in rapidly rotating stars, the results do not change qualitatively.

Modelling of soldier fly halteres for gyroscopic oscillations

Nature has evolved a beautiful design for small-scale vibratory rategyro in the form of dipteran halteres that detect body rotations via Coriolis acceleration. In most Diptera, including soldier fly, Hermetia illucens, halteres are a pair of special organs, located in the space between the thorax and the abdomen. The halteres along with their connecting joint with the fly's body constitute a mechanism that is used for muscle-actuated oscillations of the halteres along the actuation direction. These oscillations lead to bending vibrations in the sensing direction (out of the haltere's actuation plane) upon any impressed rotation due to the resulting Coriolis force. This induced vibration is sensed by the sensory organs at the base of the haltere in order to determine the rate of rotation. In this study, we evaluate the boundary conditions and the stiffness of the anesthetized halteres along the actuation and the sensing direction. We take several cross-sectional SEM (scanning electron microscope) images of the soldier fly haltere and construct its three dimensional model to get the mass properties. Based on these measurements, we estimate the natural frequency along both actuation and sensing directions, propose a finite element model of the haltere's joint mechanism, and discuss the significance of the haltere's asymmetric cross-section. The estimated natural frequency along the actuation direction is within the range of the haltere's flapping frequency. However, the natural frequency along the sensing direction is roughly double the haltere's flapping frequency that provides a large bandwidth for sensing the rate of rotation to the soldier flies.

On the control of rapidly rotating convection by an axially varying magnetic field

The magnetic field in rapidly rotating dynamos is spatially inhomogeneous. The axial variation of the magnetic field is of particular importance because tall columnar vortices aligned with the rotation axis form at the onset of convection. The classical picture of magnetoconvection with constant or axially varying magnetic fields is that the Rayleigh number and wavenumber at onset decrease appreciably from their non-magnetic values. Nonlinear dynamo simulations show that the axial lengthscale of the self-generated azimuthal magnetic field becomes progressively smaller as we move towards a rapidly rotating regime. With a small-scale field, however, the magnetic control of convection is different from that in previous studies with a uniform or large-scale field. This study looks at the competing viscous and magnetic mode instabilities when the Ekman number E (ratio of viscous to Coriolis forces) is small. As the applied magnetic field strength (measured by the Elsasser number Lambda) increases, the critical Rayleigh number for onset of convection initially increases in a viscous branch, reaches an apex where both viscous and magnetic instabilities co-exist, and then falls in the magnetic branch. The magnetic mode of onset is notable for its dramatic suppression of convection in the bulk of the fluid layer where the field is weak. The viscous-magnetic mode transition occurs at Lambda similar to 1, which implies that small-scale convection can exist at field strengths higher than previously thought. In spherical shell dynamos with basal heating, convection near the tangent cylinder is likely to be in the magnetic mode. The wavenumber of convection is only slightly reduced by the self-generated magnetic field at Lambda similar to 1, in agreement with previous planetary dynamo models. The back reaction of the magnetic field on the flow is, however, visible in the difference in kinetic helicity between cyclonic and anticyclonic vortices.

Onset of plane layer magnetoconvection at low Ekman number

We study the onset of magnetoconvection between two infinite horizontal planes subject to a vertical magnetic field aligned with background rotation. In order to gain insight into the convection taking place in the Earth's tangent cylinder, we target regimes of asymptotically strong rotation. The critical Rayleigh number Ra-c and critical wavenumber a(c) are computed numerically by solving the linear stability problem in a systematic way, with either stress-free or no-slip kinematic boundary conditions. A parametric study is conducted, varying the Ekman number E (ratio of viscous to Coriolis forces) and the Elsasser number. (ratio of the Lorentz force to the Coriolis force). E is varied from 10(-9) to 10(-2) and. from 10(-3) to 1. For a wide range of thermal and magnetic Prandtl numbers, our results verify and confirm previous experimental and theoretical results showing the existence of two distinct unstable modes at low values of E-one being controlled by the magnetic field, the other being controlled by viscosity (often called the viscous mode). It is shown that oscillatory onset does not occur in the range of parameters we are interested in. Asymptotic scalings for the onset of these modes are numerically confirmed and their domain of validity is precisely quantified. We show that with no-slip boundary conditions, the asymptotic behavior is reached for E < 10(-6) and establish a map in the (E, Lambda) plane. We distinguish regions where convection sets in either through the magnetic mode or through the viscous mode. Our analysis gives the regime in which the transition between magnetic and viscous modes may be observed. We also show that within the asymptotic regime, the role played by the kinematic boundary conditions is minimal. (C) 2015 AIP Publishing LLC.

Limiting States Of Internal Wave Rays

Existing models of baroclinic tides are based upon the "traditional approximation'', i. e., neglect of the horizontal component of the Earth's rotation, leading to a well- known conclusion that no freely propagating internal waves can exist beyond the critical latitude and the wave rays are symmetric to the vertical. However, recent studies have contended that the situation may change if both the vertical and horizontal components of the Earth's rotation are taken into account. With the full account of the Coriolis force, characteristics of the internal wavefield generated by tidal flow over uneven topography are investigated. It is found that "nontraditional effects'' profoundly change not only the dynamics of internal waves but also the rate at which the barotropic tidal energy is fed into the internal wavefield. Discarding the traditional approximation, internal waves are proved to be able to generate poleward of the critical latitude, rays of which are no longer symmetric and the limiting values of ray angles become greater or less than 90 degrees, depending on the local latitude and the direction of ray. More importantly, in contrast to the predictions of models based upon the traditional approximation, a substantial conversion occurs in the situations when stratification is so weak that the buoyancy frequency is below the tidal one.

旋转流体侵入型异重流实验研究

文中介绍了对旋转系统中侵入型异重流特征的若干实验室观察结果,在Coriolis力作用下,异重流取沿岸(壁)流形式,其宽度与其Rossby变形半径成比例,在实验的参数范围内,这种沿岸流均为不稳定的,在其演化过程中产生出气旋——反气旋涡对,乃至与岸壁分离与再附的大弯曲,本文还初步考察了波长、相对厚度与Frou de数之间的关系。

论旋转液体星的稳定

本文应用文献[1]中的方法,证明了旋转液体星(作为自吸引、不可压粘性液体处理)由所有振动本征模决定的稳定临界和长期稳定临界是一致的.同时,以Columbus问题作为这个结论的实验依据.这是Thomson和Tait在1883年预见而没有证明的重要结论,由此可以重新肯定旋转液体星准稳态演化过程必须按长期稳定临界判断;而对于近三十年来被否定的Jeans-Darwin双星分裂理论,必须重新进行评价.

熔液法晶体生长流场的数值模拟(英文)

本文给出了对二种熔液法生长晶体的液相流场进行的计算机模拟的结果,对CZ法晶体生长的定常流场数值计算,采用了二阶迎风QUICK差分格式的有限控制体积法,通过对12种不同参数的流场情况所作的计算表明,二阶QUICK格式具有更高的计算精度和数值稳定性。对离心机环境下的水平Bridgman晶体生长的流场进行了三维非定常的数值模拟。结果表明,在一定的旋转方向和旋转速率条件下,离心机环境中的Coriolis力对熔体振荡流动有致稳作用。计算结果符合以往实验的发现,并优于二维非定常模拟结果。

Nuclear hyperfine interactions in W^(182), W^(183), and Pt^(195) as studies by the Mossbauer effect

The Mössbauer technique has been used to study the nuclear hyperfine interactions and lifetimes in W¹⁸² (2⁺ state) and W¹⁸³ (3/2^- and 5/2^- states) with the following results: g(5/2^-)/g(2⁺) = 1.40 ± 0.04; g(3/2^- = -0.07 ± 0.07; Q(5/2^-)/Q(2⁺) = 0.94 ± 0.04; T_1/2(3/2^-) = 0.184 ± 0.005 nsec; T_1/2(5/2^-) >̰ 0.7 nsec. These quantities are discussed in terms of a rotation-particle interaction in W¹⁸³ due to Coriolis coupling. From the measured quantities and additional information on γ-ray transition intensities magnetic single-particle matrix elements are derived. It is inferred from these that the two effective g-factors, resulting from the Nilsson-model calculation of the single-particle matrix elements for the spin operators ŝ_z and ŝ₊, are not equal, consistent with a proposal of Bochnacki and Ogaza.

The internal magnetic fields at the tungsten nucleus were determined for substitutional solid solutions of tungsten in iron, cobalt, and nickel. With g(2⁺) = 0.24 the results are: |H_eff(W-Fe)| = 715 ± 10 kG; |H_eff(W-Co)| = 360 ± 10 kG; |H_eff(W-Ni)| = 90 ± 25 kG. The electric field gradients at the tungsten nucleus were determined for WS₂ and WO₃. With Q(2⁺) = -1.81b the results are: for WS₂, eq = -(1.86 ± 0.05) 10¹⁸ V/cm²; for WO₃, eq = (1.54 ± 0.04) 10¹⁸ V/cm² and ƞ = 0.63 ± 0.02.

The 5/2^- state of Pt¹⁹⁵ has also been studied with the Mössbauer technique, and the g-factor of this state has been determined to be -0.41 ± 0.03. The following magnetic fields at the Pt nucleus were found: in an Fe lattice, 1.19 ± 0.04 MG; in a Co lattice, 0.86 ± 0.03 MG; and in a Ni lattice, 0.36 ± 0.04 MG. Isomeric shifts have been detected in a number of compounds and alloys and have been interpreted to imply that the mean square radius of the Pt¹⁹⁵ nucleus in the first-excited state is smaller than in the ground state.

Movement of deposited water on turbomachinery rotor blade surfaces

A theoretical approach for calculating the movement of liquid water following deposition onto a turbomachine rotor blade is described. Such a situation can occur during operation of an aero-engine in rain. The equation of motion of the deposited water is developed on an arbitrarily oriented plane triangular surface facet. By dividing the blade surface into a large number of facets and calculating the water trajectory over each one crossed in turn, the overall trajectory can be constructed. Apart from the centrifugal and Coriolis inertia effects, the forces acting on the water arise from the blade surface friction, and the aerodynamic shear and pressure gradient. Non- dimensionalisation of the equations of motion provides considerable insight and a detailed study of water flow on a flat rotating plate set at different stagger angles demonstrates the paramount importance of blade surface friction. The extreme cases of low and high blade friction are examined and it is concluded that the latter (which allows considerable mathematical generalisation) is the most likely in practice. It is also shown that the aerodynamic shear force, but not the pressure force, may influence the water motion. Calculations of water movement on a low-speed compressor blade and the fan blade of a high bypass ratio aero-engine suggest that in low rotational speed situations most of the deposited water is centrifuged rapidly to the blade tip region. Copyright © 2006 by ASME.

Spatial structure of internal tidal waves generated on weak topographies

The generation of internal gravity waves by barotropic tidal flow passing over a two-dimensional topography is investigated. Rather than calculating the conversion of tidal energy, this study focuses on delineating the geometric characteristics of the spatial structure of the resulting internal wave fields (i.e., the configurations of the internal beams and their horizontal projections) which have usually been ignored. it is found that the various possible wave types can be demarcated by three characteristic frequencies: the tidal frequency, wo; the buoyancy frequency, N; and the vertical component of the Coriolis vector or earth's rotation.f. When different possibilities arising from the sequence of these frequencies are considered, there occur 12 kinds of wave structures in the full 3D space in contrast to the 5 kinds identified by the 2D theory. The constant wave phase lines may form as ellipses or hyperbolic lines on the horizontal plane, provided the buoyancy frequency is greater or less than the tidal frequency. The effect that stems from the consideration of the basic flow is also found, which not only serves as the reason for the occurrence of higtter harmonics but also increases the wave strength in the direction of basic flow. (C) 2009 Elsevier B.V. All rights reserved.

Two modes of dipole events in tropical Indian Ocean

By analyzing the distributions of subsurface temperature and the surface wind stress anomalies in the tropical Pacific and Indian Oceans during the Indian Ocean Dipole (IOD) events, two major modes of the IOD and their formation mechanisms are revealed. (1) The subsurface temperature anomaly (STA) in the tropical Indian Ocean during the IOD events can be described as a "<" -shaped and west-east-oriented dipole pattern; in the east side of the "<" pattern, a notable tongue-like STA extends westward along the equator in the tropical eastern Indian Ocean; while in the west side of the "<" pattern, the STA has opposite sign with two centers (the southern one is stronger than the northern one in intensity) being of rough symmetry about the equator in the tropical mid-western Indian Ocean. (2) The IOD events are composed of two modes, which have similar spatial pattern but different temporal variabilities due to the large scale air-sea interactions within two independent systems. The first mode of the IOD event originates from the air-sea interaction on a scale of the tropical Pacific-Indian Ocean and coexists with ENSO. The second mode originates from the air-sea interaction on a scale of the tropical Indian Ocean and is closely associated with changes in the position and intensity of the Mascarene high pressure. The strong IOD event occurs when the two modes are in phase, and the IOD event weakens or disappears when the two modes are out of phase. Besides, the IOD events are normally strong when either of the two modes is strong. (3) The IOD event is caused by the abnormal wind stress forcing over the tropical Indian Ocean, which results in vertical transports, leading to the upwelling and pileup of seawater. This is the main dynamic processes resulting in the STA. When the anomalous easterly exists over the equatorial Indian Ocean, the cold waters upwell in the tropical eastern Indian Ocean while the warm waters pileup in the tropical western Indian Ocean, hence the thermocline in the tropical Indian Ocean is shallowed in the east and deepened in the west. The off-equator component due to the Coriolis force in the equatorial area causes the upwelling of cold waters and the shallowing of the equatorial India Ocean thermocline. On the other hand, the anomalous anticyclonic circulations and their curl fields located on both sides of the equator, cause the pileup of warm waters in the central area of their curl fields and the deepening of the equatorial Indian Ocean thermocline off the equator. The above three factors lead to the occurrence of positive phase IOD events. When anomalous westerly dominates over the tropical Indian Ocean, the dynamic processes are reversed, and the negative-phase IOD event occurs.

The effects of random surface waves on the steady Ekman current solutions

The response of near-surface current profiles to wind and random surface waves are studied based on the approach of Jenkins [1989. The use of a wave prediction model for driving a near surface current model. Dtsch. Hydrogr. Z. 42,134-149] and Tang et al. [2007. Observation and modeling of surface currents on the Grand Banks: a study of the wave effects on surface currents. J. Geophys. Res. 112, C10025, doi:10.1029/2006JC004028]. Analytic steady solutions are presented for wave-modified Ekman equations resulting from Stokes drift, wind input and wave dissipation for a depth-independent constant eddy viscosity coefficient and one that varies linearly with depth. The parameters involved in the solutions can be determined by the two-dimensional wavenumber spectrum of ocean waves, wind speed, the Coriolis parameter and the densities of air and water, and the solutions reduce to those of Lewis and Belcher [2004. Time-dependent, coupled, Ekman boundary layer solutions incorporating Stokes drift. Dyn. Atmos. Oceans. 37, 313-351] when only the effects of Stokes drift are included. As illustrative examples, for a fully developed wind-generated sea with different wind speeds, wave-modified current profiles are calculated and compared with the classical Ekman theory and Lewis and Belcher's [2004. Time-dependent, coupled, Ekman boundary layer solutions incorporating Stokes drift. Dyn. Atmos. Oceans 37, 313-351] modification by using the Donelan and Pierson [1987. Radar scattering and equilibrium ranges in wind-generated waves with application to scatterometry. J. Geophys. Res. 92, 4971-5029] wavenumber spectrum, the WAM wave model formulation for wind input energy to waves, and wave energy dissipation converted to currents. Illustrative examples for a fully developed sea and the comparisons between observations and the theoretical predictions demonstrate that the effects of the random surface waves on the classical Ekman current are important, as they change qualitatively the nature of the Ekman layer. But the effects of the wind input and wave dissipation on surface current are small, relative to the impact of the Stokes drift. (C) 2008 Elsevier Ltd. All rights reserved.