Oscillatory Regime Research Articles

Transition from random self-propulsion to rotational motion in a non-Markovian microswimmer

Abstract We study the motion of an inertial microswimmer in a non-Newtonian environment with a finite memory and present the theoretical realization of an unexpected transition from its random self-propulsion to rotational (circular or elliptical) motion. Further, the rotational motion of the swimmer is followed by spontaneous local direction reversals yet with a steady state angular diffusion. Moreover, the advent of this behaviour is observed in the oscillatory regime of the inertia-memory parameter space of the dynamics. We quantify this unconventional rotational motion of microswimmer by measuring the time evolution of direction of its instantaneous velocity or orientation. By solving the generalized Langevin model of non-Markovian dynamics of an inertial active Ornstein-Uhlenbeck particle, we show that the emergence of the rotational (circular or elliptical) trajectory is due to the presence of both inertial motion and memory in the environment.

Analysis of flutter dynamics in thin flexible flags under streamlined and vortex-induced flows

In this study, numerical and experimental investigations are conducted to gain a thorough understanding of the aeroelastic behavior of flexible steel flags at low air speeds. Numerical simulations are conducted using a two-dimensional computational fluid dynamics solver to analyze the vortex shedding mechanism from a circular bluff body and to assess its effect on a downstream slender body placed at different gap distances. Experiments on a low-speed wind tunnel are performed to examine the flow-induced flutter dynamics of thin flexible flags. Four different flag shapes, namely, triangular, square, rectangular, and trapezoidal, are tested both with and without an upstream bluff body. Under streamlined flow conditions, flutter instability is characterized by subcritical bifurcation routes, while the presence of a cylindrical bluff body induces supercritical bifurcation scenarios. The aeroelastic response of the flag is inspected while varying the gap distance, separating it from the cylindrical bluff body from 0.5 to 15 times the cylinder's diameter. Experimental results show that the shape of the flag significantly influences its dynamic behavior, especially when it is placed in close proximity to the bluff body. The triangular flag exhibits the highest limit cycle oscillation (LCO) amplitudes, while the square flag has the lowest flutter onset speed. A regime of period-2 LCOs is observed in the triangular flag configuration, which is not seen elsewhere. For flags positioned farther away, the dynamics resembles those observed under streamlined flow, due to the vortices weakening with the increasing gap size. Finally, a comparative analysis demonstrates the flow speed regimes of high-amplitude oscillations for different flag shapes and bluff body positions, which can be explored for low speed energy harvesting applications.

Probing the contact time of droplet impacts: From the Hertz collision to oscillation regimes.

Droplet rebound on nonwetting surfaces is a common phenomenon. However, the underlying physics regulating the contact time remains unclear. In this work, we investigate droplet impacts on superamphiphobic surfaces through experiments and theoretical analyses. By analyzing the spreading and retraction of droplet impinging processes over a wide range of Weber numbers (We), it is revealed that droplet impacts experience three regimes as We is varied, which are denoted as the Hertz collision (We<1), transition (1<We<10), and oscillation (We>10) regimes. In the Hertz collision regime, the droplet impinging process is temporally symmetric, i.e., the spreading time, t_{S}, and the retraction time, t_{R}, are almost the same. Furthermore, t_{S} and t_{R} decrease with increasing We and follow a power-law dependence, which is different from previous theories. In the transition regime, t_{S} remains dominated by the Hertz collision, while t_{R} is governed by droplet oscillation. In the oscillation regime, t_{S}, t_{R}, and the total contact time, t_{C}, become independent of We. These three regimes are valid for both monophase and compound droplets. The findings in this work advance the understanding and offer a clear picture of droplet impact dynamics.

Capillary Orientation and Morphology Dictated Oscillatory Electro-magneto-imbibition of Viscoelastic Electrolytes.

We explore the intricate dynamics of imbibition by a viscoelastic electrolyte within an arbitrarily oriented, nonuniform microcapillary, while under the stimulus of external electromagnetic fields and internal electroviscous forces stemming from streaming potential. The microcapillary walls are envisaged to be tapered relative to each other, with the entire system inclined with respect to the horizontal plane. The rheological behavior of the electrolyte is characterized using the Phan-Thien-Tanner (PTT) model. To manipulate the imbibition dynamics, external transverse magnetic and electric fields are imposed. Incorporating all contributing forces, we obtain semianalytical formulations for the velocity and shear stress distributions. We identify distinct stages during the imbibition process: (i) initial, (ii) filling, (iii) oscillatory, and (iv) stagnation stages. Moreover, we scrutinize the impact of four pivotal parameters, namely, the Weissenberg number (Wi), the Hartmann number (Ha), the transverse electric to viscous force (EVF) number (S), and the electric Reynolds number (Ree), on the imbibition dynamics across different inclinations and taper angles. We also delineate the parameter space for these four parameters, identifying where the onset of oscillations occurs. Finally, through scaling analysis, we establish the existence of four distinct regimes corresponding to the aforementioned stages: (i) the linear regime, (ii) the Washburn regime, (iii) the oscillatory regime, and (iv) the equilibrium regime. Our findings are anticipated to enhance the understanding of capillary imbibition under such complex flow conditions and contribute significantly to the advancement of capillarity-driven microfluidic devices.

Tailoring femtosecond LSP resonance and lifetime in a nanoresonator via phase retardation.

The manipulation of femtosecond plasmon resonance and lifetime in a nanoantenna is crucial for the realization of integrated and miniaturized plasmonic circuits. Here, we have used FDTD simulations to study the plasmonic resonance and lifetime variation of the far-end (PA) and near-end (PB) hotspots of size-gradient nanoresonators. We found that the near-field spectrum of PA is red-shifted compared to PB due to the phase retardation effect. By capturing the ultrafast dynamics of both hotspots, we confirm that these phenomena are governed by the transient evolution of the plasmonic field in the forced oscillation regime. Furthermore, the lifetimes of plasmonic hotspots scale directly with their near-field intensities. Meanwhile, the lifetime τA is always larger than that of τB in the same nanoresonator because of the smaller non-radiative damping of hotspot PA. Our results provide a basis for the miniaturization of plasmonic nanoresonators.

Operation regimes of spinal circuits controlling locomotion and the role of supraspinal drives and sensory feedback

Locomotion in mammals is directly controlled by the spinal neuronal network, operating under the control of supraspinal signals and somatosensory feedback that interact with each other. However, the functional architecture of the spinal locomotor network, its operation regimes, and the role of supraspinal and sensory feedback in different locomotor behaviors, including at different speeds, remain unclear. We developed a computational model of spinal locomotor circuits receiving supraspinal drives and limb sensory feedback that could reproduce multiple experimental data obtained in intact and spinal-transected cats during tied-belt and split-belt treadmill locomotion. We provide evidence that the spinal locomotor network operates in different regimes depending on locomotor speed. In an intact system, at slow speeds (&lt;0.4 m/s), the spinal network operates in a non-oscillating state-machine regime and requires sensory feedback or external inputs for phase transitions. Removing sensory feedback related to limb extension prevents locomotor oscillations at slow speeds. With increasing speed and supraspinal drives, the spinal network switches to a flexor-driven oscillatory regime and then to a classical half-center regime. Following spinal transection, the model predicts that the spinal network can only operate in the state-machine regime. Our results suggest that the spinal network operates in different regimes for slow exploratory and fast escape locomotor behaviors, making use of different control mechanisms.

Operation regimes of spinal circuits controlling locomotion and the role of supraspinal drives and sensory feedback.

Locomotion in mammals is directly controlled by the spinal neuronal network, operating under the control of supraspinal signals and somatosensory feedback that interact with each other. However, the functional architecture of the spinal locomotor network, its operation regimes, and the role of supraspinal and sensory feedback in different locomotor behaviors, including at different speeds, remain unclear. We developed a computational model of spinal locomotor circuits receiving supraspinal drives and limb sensory feedback that could reproduce multiple experimental data obtained in intact and spinal-transected cats during tied-belt and split-belt treadmill locomotion. We provide evidence that the spinal locomotor network operates in different regimes depending on locomotor speed. In an intact system, at slow speeds (<0.4 m/s), the spinal network operates in a non-oscillating state-machine regime and requires sensory feedback or external inputs for phase transitions. Removing sensory feedback related to limb extension prevents locomotor oscillations at slow speeds. With increasing speed and supraspinal drives, the spinal network switches to a flexor-driven oscillatory regime and then to a classical half-center regime. Following spinal transection, the model predicts that the spinal network can only operate in the state-machine regime. Our results suggest that the spinal network operates in different regimes for slow exploratory and fast escape locomotor behaviors, making use of different control mechanisms.

Learning the latent dynamics of fluid flows from high-fidelity numerical simulations using parsimonious diffusion maps

We use parsimonious diffusion maps (PDMs) to discover the latent dynamics of high-fidelity Navier–Stokes simulations with a focus on the two-dimensional (2D) fluidic pinball problem. By varying the Reynolds number Re, different flow regimes emerge, ranging from steady symmetric flows to quasi-periodic asymmetric and chaos. The proposed non-linear manifold learning scheme identifies in a crisp manner the expected intrinsic dimension of the underlying emerging dynamics over the parameter space. In particular, PDMs estimate that the emergent dynamics in the oscillatory regime can be captured by just two variables, while in the chaotic regime, the dominant modes are three as anticipated by the normal form theory. On the other hand, proper orthogonal decomposition/principal component analysis (POD/PCA), most commonly used for dimensionality reduction in fluid mechanics, does not provide such a crisp separation between the dominant modes. To validate the performance of PDMs, we also compute the reconstruction error, by constructing a decoder using geometric harmonics (GHs). We show that the proposed scheme outperforms the POD/PCA over the whole Re number range. Thus, we believe that the proposed scheme will allow for the development of more accurate reduced order models for high-fidelity fluid dynamics simulators, relaxing the curse of dimensionality in numerical analysis tasks such as bifurcation analysis, optimization, and control.

Investigating ocean circulation dynamics through data assimilation: A mathematical study using the Stommel box model with rapid oscillatory forcings.

We investigate ocean circulation changes through the lens of data assimilation using a reduced-order model. Our primary interest lies in the Stommel box model, which reveals itself to be one of the most practicable models that has the ability of reproducing the meridional overturning circulation. The Stommel box model has at most two regimes: TH (temperature driven circulation with sinking near the north pole) and SA (salinity driven with sinking near the equator). Currently, the meridional overturning is in the TH regime. Using box-averaged Met Office EN4 ocean temperature and salinity data, our goal is to provide a probability that a future regime change occurs and establish how this probability depends on the uncertainties in initial conditions, parameters, and forcings. We will achieve this using data assimilation tools and DAPPER within the Stommel box model with fast oscillatory regimes.

Marginal stability analyses for thermochemical convection and its implications for the dynamics of continental lithosphere and core–mantle boundary regions

SUMMARY Significant compositional differences may exist in the lithospheric mantle and above the core–mantle boundary (CMB) relative to the ambient mantle. The intrinsic density differences may affect the development of thermal boundary layer (TBL) instabilities associated with lithospheric delamination and formation of thermochemical plumes. In this study, we explored the instability of two-layer thermochemical fluid using two different techniques: marginal stability analysis with a propagator-matrix method and finite element modelling. We investigated both the instabilities in lithospheric mantle (i.e. lithospheric instability) and the mantle above the CMB (i.e. plume-forming instability) using a background temperature Tbg(z) with the TBL. For lithospheric instability, we found that two-layer fluid with free-slip boundary conditions mainly undergoes the same three different convective modes (i.e. two oscillatory convection modes and one layered convection regime) as that with no-slip boundary conditions reported in Jaupart et al. However, with free-slip boundary conditions, the transitions between these convection modes occur at larger values of buoyancy number B. Free-slip boundary conditions lead to smaller critical Rayleigh number Rac, but larger convective wavelength and oscillation frequency ωc, compared with those with no-slip boundary conditions. Our numerical modelling results demonstrate that Rac and ωc predicted from the classical marginal stability analyses using Tbg(z) with TBL temperature may have significant errors when the oscillatory period is comparable with or larger than the timescale of lithospheric thermal diffusion that causes Tbg(z) to vary with time significantly. In this case, using a more gently sloped background temperature profile ignoring the TBL temperature, the stability analysis predicts more accurate stability conditions, thus presenting an effective remedy to the stability analysis. For plume-forming instability, because of the reduced viscosity in the hot and compositionally dense bottom layer, the transition to the layered convection occurs at significantly smaller B values, and in the oscillatory convection regime, Rac is larger but ωc is smaller, compared with those for lithospheric instability. Finally, our study provides a successful benchmark of numerical models of thermochemical convection by comparing Rac and ωc from numerical models with those from the marginal stability analysis.

Optimally truncated WKB approximation for the 1D stationary Schrödinger equation in the highly oscillatory regime

This paper is dedicated to the efficient numerical computation of solutions to the 1D stationary Schrödinger equation in the highly oscillatory regime. We compute an approximate solution based on the well-known WKB-ansatz, which relies on an asymptotic expansion w.r.t. the small parameter ɛ. Assuming that the coefficient in the equation is analytic, we derive an explicit error estimate for the truncated WKB series, in terms of ɛ and the truncation order N. For any fixed ɛ, this allows to determine the optimal truncation order Nopt which turns out to be proportional to ɛ−1. When chosen this way, the resulting error of the optimally truncated WKB series behaves like O(exp(−r/ɛ)), with some parameter r>0. The theoretical results established in this paper are confirmed by several numerical examples.

Second register production on the clarinet: Nonlinear losses in the register hole as a decisive physical phenomenon.

This study investigates the role of localized nonlinear losses in the register hole of a clarinet in producing second-register notes. First, an experiment is conducted to study the ability of the opening of a register hole to trigger a jump in oscillatory regime from the first to the second register. A cylindrical tube is drilled with holes of increasing diameter: five at the register hole level and five at the thumb hole level of a B-flat clarinet. Clarinetists are asked to play with constant parameters, blindfolded, beginning with all holes closed. The operator randomly opens one of the ten holes, noting the resulting register. The experiment is replicated numerically by time integration of two different models. The first is the model from Taillard, Silva, Guillemain, and Kergomard [(2018). Appl. Acoust. 141, 271-280] based on the modal decomposition of the input impedance. The second accounts for localized nonlinear losses in the register hole, through the model from Dalmont, Nederveen, Dubos, Ollivier, Méserette, and Sligte [(2002). Acta Acust. united Ac. 88, 567-575]. These losses are handled through variable modal coefficients. For the first model, simulations never produce the second register for any of the open holes. For the second, the proportion of second-register production is close to the experiment for upstream holes, but remains at zero for downstream holes.

Operation regimes of spinal circuits controlling locomotion and role of supraspinal drives and sensory feedback.

Locomotion in mammals is directly controlled by the spinal neuronal network, operating under the control of supraspinal signals and somatosensory feedback that interact with each other. However, the functional architecture of the spinal locomotor network, its operation regimes, and the role of supraspinal and sensory feedback in different locomotor behaviors, including at different speeds, remain unclear. We developed a computational model of spinal locomotor circuits receiving supraspinal drives and limb sensory feedback that could reproduce multiple experimental data obtained in intact and spinal-transected cats during tied-belt and split-belt treadmill locomotion. We provide evidence that the spinal locomotor network operates in different regimes depending on locomotor speed. In an intact system, at slow speeds (< 0.4 m/s), the spinal network operates in a non-oscillating state-machine regime and requires sensory feedback or external inputs for phase transitions. Removing sensory feedback related to limb extension prevents locomotor oscillations at slow speeds. With increasing speed and supraspinal drives, the spinal network switches to a flexor-driven oscillatory regime and then to a classical half-center regime. Following spinal transection, the model predicts that the spinal network can only operate in the state-machine regime. Our results suggest that the spinal network operates in different regimes for slow exploratory and fast escape locomotor behaviors, making use of different control mechanisms.

Effect of the dynamic contact angle on electromagnetically driven flows in free liquid films

We study the effect of a dynamic contact angle on an electromagnetically driven flow in a horizontal free electrolyte film stretched between two coaxial cylindrical electrodes and placed in a uniform magnetic field. The flow dynamics and film deformation are described using the reduced hydrodynamic model derived in the lubrication approximation in [A. Pototsky and S. A. Suslov, J. Fluid Mech., 984, A75 (2024)]. The linearized molecular kinetic model is used to relate the dynamic and static contact angles to the wetting-line friction coefficient. Steady azimuthal flow is found for arbitrary static contact angles. Linear stability of the base azimuthal flow with respect to azimuthally invariant perturbations is studied using the numerical continuation method. We find that the flow stability is highly sensitive to variations of the wetting-line friction coefficient and the static contact angle. The least stable azimuthal flow is found for the frictionless contact line corresponding to a free film that remains perpendicular to the surface of the electrodes. The azimuthal flows are found to experience a supercritical Hopf bifurcation upon which they transition to a stable oscillatory dynamic regime characterized by alternating squeezing and swelling of the film near the inner and outer electrodes.

Complex viscosity and cohesive energy density of mixtures of sodium hyaluronate and hypromellose for medical devices

In the present study binary formulations of sodium hyaluronate (▪) and hypromellose (▪) were proposed as medical devices in the field of eye surgeries. Viscoelasticity is tested directly on binary solutions while the rheological behavior in the human eye is mimicked by monitoring the viscoelasticity of ternary mixtures ▪/▪/▪. Both binary and ternary formulations were rheologically characterized at 25∘C and 37∘C by means of measurements in the oscillatory regime. Amplitude sweep tests revealed that the range of linear viscoelasticity is strongly dependent on the concentration of ▪. For the proposed formulations, it was ascertained that the cohesive energy density at 37∘C decreases as ▪ concentration surges; in contrast, at 25∘C we detected an increase of cohesive energy were detected. At ▪ mass fraction of α1, the two obtained values were comparable. Frequency sweep tests revealed the pseudoplastic character of the binary formulations. In order to fit the complex viscosity modules, the Carreau-Yasuda model was modified. Numerical analysis of the results indicates that the best-fit to the model is obtained when the characteristic parameters take constant values, i.e. n=0 and m=1, regardless of the composition of the formulation. This implies that the modulus of complex viscosity decays hyperbolically. Frequency sweep tests on the ternary formulations were analyzed in terms of both the modulus of complex viscosity to confirm the previous scientific evidence together with an investigation of the phase shift angle between the viscous and elastic components. The results indicate that the addition of water has a different impact on the structure and the cohesion of the mixture at different temperatures. These results suggest that the ternary mixture structure can be sensitively modulated by temperature and the content of water.

Heat Flow Overshooting in Suspended Graphene

ABSTRACT In this work we found that when the hydrodynamic heat conduction in suspended graphene is in the under-damped free oscillation regime the disturbance of heat flux at the boundaries of graphene and the joint action of the initial heat flux and its time change rate can cause the heat flow overshooting which refers to the excess heat flux established in the heat conduction medium when two heat flow wave-fronts meet. The overshooting can make the maximum amplitude of the heat flux more than twice that of the initial heat flux, which poses a great danger to the next high-frequency and high-efficiency electronic devices based on graphene, such as the graphene field-effect transistors. Interestingly, results show that when the scale ratio of graphene with rectangular shape in the heat flow direction to its vertical direction is greater than a certain value, the overshooting caused by the disturbance of the heat flux at the boundaries no longer occurs. For the overshooting induced by the uneven distribution of the initial heat flux and its time change rate, it is shown that if the length in the region with the higher heat flux and time change rate in the direction of heat flow is less than a critical value the overshooting phenomenon can be avoided. These findings provide a new way for thermal management of electronic devices based on graphene.

Unsupervised detection of large-scale weather patterns in the northern hemisphere via Markov State Modelling: from blockings to teleconnections

Detecting recurrent weather patterns and understanding the transitions between such regimes are key to advancing our knowledge of the low-frequency variability of the atmosphere and have important implications in terms of weather and climate-related risks. We adopt an analysis pipeline inspired by Markov State Modelling and detect in an unsupervised manner the dominant winter mid-latitude Northern Hemisphere weather patterns in the Atlantic and Pacific sectors. The daily 500 hPa geopotential height fields are first classified in about 200 microstates. The weather dynamics are then represented on the basis of these microstates and the slowest decaying modes are identified from the spectral properties of the transition probability matrix. These modes are defined on the basis of the nonlinear dynamical processes of the system and not as tentative metastable states, as often done in Markov state analysis. When focusing on a shifting longitudinal window of 60∘, we find that the longitude-dependent estimate of the longest relaxation time is smaller where stronger baroclinic activity is found. In the Atlantic and Pacific sectors slow relaxation processes are mainly related to transitions between blocked regimes and zonal flow. We also find strong evidence of a dynamical regime associated with the simultaneous Atlantic-Pacific blocking. When the analysis is performed on a broader geographical region of the Atlantic sector, we discover that the slowest relaxation modes of the system are associated with transitions between dynamical regimes that resemble teleconnection patterns like the North Atlantic Oscillation and weather regimes like the Scandinavian and Greenland blocking, yet have a much stronger dynamical foundation than classical methods based e.g. on EOF analysis. Our method clarifies that, as a result of the lack of a time-scale separation in the atmospheric variability of the mid-latitudes, there is no clear-cut way to represent the atmospheric dynamics in terms of few, well-defined modes of variability. The approach proposed here can be seamlessly applied across different regions of the globe for detecting regional modes of variability, and has a great potential for intercomparing climate models and for assessing the impact of climate change on the low-frequency variability of the atmosphere.

Концентрационная конвекция в замкнутой области пористой среды при заданном вертикальном перепаде концентрации и учете иммобилизации примеси

This work is devoted to the study of stability of horizontal filtration flow of a mixture through a closed porous domain taking into account impurity immobilization. The instability arises due to the vertical concentration difference of heavy impurities, which creates unstable density stratification. A general mathematical model describing the transport of an impurity through a porous medium is presented. The equations are simplified for the case of low impurity concentration. Simplification made it possible to analytically obtain the solution corresponding to homogeneous horizontal filtration and to study its stability. It is known that, in the narrow regions of porous medium and at weak intensities of the external flow, convection is excited in a monotonous manner. On the contrary, in the case of an infinite horizontal layer, oscillatory instability is observed. A study of the transition between the instability modes is presented. It is shown that the oscillatory regime is observed in long regions or at significant intensity of the external horizontal flow. At low flow intensities, convective cells do not move relative to the region and, hence, there is no reason for oscillations. It has been established that the range of flow intensity values, in which oscillations are observed, grows with increasing length of the domain. Impurity immobilization leads to the stabilization of horizontal filtration with respect to convective perturbations. Critical curves and stability maps are obtained in a wide range of problem parameters and then analyzed. For limiting cases, a comparison is made with the known results obtained for an infinite layer.

A Holocene history of climate, fire, landscape evolution, and human activity in northeastern Iceland

Abstract. Paleoclimate reconstructions across Iceland provide a template for past changes in climate across the northern North Atlantic, a crucial region due to its position relative to the global northward heat transport system and its vulnerability to climate change. The roles of orbitally driven summer cooling, volcanism, and human impact as triggers of local environmental changes in the Holocene of Iceland remain debated. While there are indications that human impact may have reduced environmental resilience during late Holocene summer cooling, it is still difficult to resolve to what extent human and natural factors affected Iceland's late Holocene landscape instability. Here, we present a continuous Holocene fire record of northeastern Iceland from proxies archived in Stóra Viðarvatn sediment. We use pyrogenic polycyclic aromatic hydrocarbons (pyroPAHs) to trace shifts in fire regimes, paired with continuous biomarker and bulk geochemical records of soil erosion, lake productivity, and human presence. The molecular composition of pyroPAHs and a wind pattern reconstruction indicate a naturally driven fire signal that is mostly regional. Generally low fire frequency during most of the Holocene significantly increased at 3 ka and again after 1.5 ka BP before known human settlement in Iceland. We propose that shifts in vegetation type caused by cooling summers over the past 3 kyr, in addition to changes in atmospheric circulation, such as shifts in North Atlantic Oscillation (NAO) regime, led to increased aridity and biomass flammability. Our results show no evidence of faecal biomarkers associated with human activity during or after human colonisation in the 9th century CE. Instead, faecal biomarkers follow the pattern described by erosional proxies, pointing toward a negligible human presence and/or a diluted signal in the lake's catchment. However, low post-colonisation levels of pyroPAHs, in contrast to an increasing flux of erosional bulk proxies, suggest that farming and animal husbandry may have suppressed fire frequency by reducing the spread and flammability of fire-prone vegetation (e.g. heathlands). Overall, our results describe a fire frequency heavily influenced by long-term changes in climate through the Holocene. They also suggest that human colonisation had contrasting effects on the local environment by lowering its resilience to soil erosion while increasing its resilience to fire.

Vortex-induced rotation of a square cylinder under the influence of Reynolds number and density ratio

Numerical simulations are carried out on the vortex-induced rotations of a freely rotatable rigid square cylinder in a two-dimensional uniform cross-flow. A range of Reynolds numbers between 40 and 150 and density ratios between 0.1 and 10 are considered. Results show eight different characteristic regimes, expanding the classification of Ryu &amp; Iaccarino (J. Fluid Mech., vol. 813, 2017, pp. 482–507). New regimes include the transition and wavy rotation regimes; in the ${\rm \pi}$ -limited oscillation regime we observe multipeak subregimes. Moment-generating mechanisms of these regimes and subregimes are further elucidated. A phenomenon related to the influence of density ratio is the tooth-like shape of the ${\rm \pi} /2$ -limit oscillation regime observed in the regime map, which is explained as a result of the imbalance relation between the main frequencies of rotation response and the vortex shedding frequency. In addition, existence of multiple regimes and multistable states are discussed, indicating multiple stable attractive structures in phase space.