Oscillation Modes Research Articles

Coaxial air blast atomization of a particulate gel suspension jet

The objective of this study is to experimentally and theoretically explore the impact of particle volume fractions and nozzle thread structure on the coaxial air blast atomization of particulate gel suspension jets through the high-speed flow visualization technology. The four breakup modes of particulate gel suspension jets are confirmed, namely oscillation mode, membrane-type breakup, fiber-type breakup, and superpulsating submode. However, the fiber-type breakup of jets disappears as the particle concentration reaches 40%, primarily attributed to the numerous particles sharply promoting the instability on the free jet surface. The experimental statistical results demonstrate that an increase in particle concentration and the adoption of the threaded nozzle both result in an expansion of the jet spray angle. The breakup length exhibits a reduction trend with the increase of particle concentration. However, it remains nearly the same value as that in the case of volume fraction 20% when the particle concentration continues to increase, contributed by the surge of viscous dissipation counteracting the enhanced instability on the jet interface caused by a large number of particles. The breakup length of a pure gel jet is predicted through the linear instability analysis, which is further modified by the viscosity model of particle suspensions to derive the breakup lengths of particulate gel suspension jets with different concentrations. The theoretical predictions align well with the statistical results in the experiment. The research is poised to have potential implications for numerous industrial processes and engineering applications including gel propellants.

High-power in-phase and anti-phase mode emission from linear arrays of resonant-tunneling-diode oscillators in the 0.4-to-0.8-THz frequency range

Oscillators based on resonant tunneling diodes (RTDs) are able to reach the highest oscillation frequency among all electronic THz emitters. However, the emitted power from RTDs remains limited. Here, we propose linear RTD oscillator arrays capable of supporting coherent emission from both in-phase and anti-phase coupled modes. The oscillation modes can be selected by adjusting the mesa areas of the RTDs. Both the modes exhibit constructive interference at different angles in the far field, enabling high-power emission. Experimental demonstrations of coherent emission from linear arrays containing 11 RTDs are presented. The anti-phase mode oscillates at ∼450 GHz, emitting about 0.7 mW, while the in-phase mode oscillates at around 750 GHz, emitting about 1 mW. Moreover, certain RTD oscillator arrays exhibit dual-band operation: changing the bias voltage allows for controllable switching between the anti-phase and in-phase modes. Upon bias sweeping in both directions, a notable hysteresis feature is observed. Our linear RTD oscillator array represents a significant step forward in the realization of large arrays for applications requiring continuous-wave THz radiation with substantial power.

Bubble 36Ar and its new breathing modes

The bubble nuclei are important components of exotic nuclear structures characterized by special depletions of central densities. Focusing on bubble structures of 36Ar, the characterizations of bubble nuclei were explored with the framework of the extended quantum molecular dynamics model. Three density distribution modes were uncovered for the first time, i.e. micro-bubble, bubble, and cluster resonances, which show unique spectral signature compared to the monopole resonance spectrum as the excitation intensity was increased in bubbles. Of pivotal importance is the revelation that the bubble mode's oscillation frequency closely resembles macroscopic bubble dynamics, building a connection between classical macroscopic phenomena and the quantum complexity of the nuclear structure. The discovery marks a crucial step forward in deciphering the relationship between classical and quantum domains within the enigmatic world of atomic nuclei.

Oscillatory onset of rotating thermal convection subject to spatially varying magnetic fields and stable stratification

Convective fluid motions in deep planetary cores induce spatially and temporally varying magnetic fields observable as secular variation. Oscillatory instabilities occurring as onset modes in rapidly rotating thermal convection influenced by heterogeneous magnetic fields and stratification, investigated in this study, can be potentially linked to such observations. A plane layer approximation to near-polar and near-equatorial regions of the Earth's outer core is considered. In addition to penetrative convection and flow suppression effects, the simultaneous interaction of convective instabilities with asymmetrical magnetic fields and stable stratification is the focus of this study. The fundamental modes of magnetoconvection onset are preserved even under stable stratification, although the convection threshold is lowered irrespective of the parameter regime. The axial invariance of rapidly rotating columnar convection loses its midplane symmetry with intense localization of kinetic energy in the unstably stratified regions. The parameter regimes supporting the onset of the magnetically modified viscous oscillatory (mVO) modes are further extended toward Earthlike conditions such as high thermal conductivity (low Prandlt number, Pr) and magnetic dissipation (low Roberts number, q). Moreover, compared to fully unstable stratification, partial stable stratification enhances the range of imposed magnetic field length scale (δ) over which the onset of mVO modes is favored. The critical field length scale (δ⋆), above which the onset regime of the mVO modes is bounded by a minimum q, is determined. Irrespective of the rotation rate, below δ=δ⋆, the onset of the mVO mode is supported for asymptotically small q for Pr values close to the Earth's outer core.

Cross-Scale Phase Relationship of the Ca II K Index with Solar Wind Parameters: A Space Climate Focus

The solar wind, representing one of the most impacting phenomena in the circum-terrestrial space, constitutes one of the several manifestations of the magnetic activity of the Sun. With the aim of shedding light on the scales beyond the rotational period of the Sun (i.e., Space Climate scales), this study investigates the phase relationship of a solar activity physical proxy, the Ca II K index, with solar wind properties measured near the Earth, over the whole space era (last five solar cycles). Using a powerful tool such as the Hilbert–Huang transform, we investigate the dependence of their phase coherence on the obtained time scale components. Phase coherence at the same time scales is found between all the components and is also preserved between adjacent components with time scales ≳ 2 yrs. Finally, given the availability of the intrinsic modes of oscillation, we explore how the relationship of Ca II K index with solar wind parameters depends on the time scale considered. According to our results, we hypothesize the presence of a bifurcation in the phase-space Ca II K index vs. solar wind speed (dynamic pressure), where the time scale seems to act as a bifurcation parameter. This concept may be pivotal for unraveling the complex interplay between solar activity and solar wind, bearing implications from the prediction and the interpretation point of view in Space Climate studies.

Experimental evaluation of the influence of the diffuser inlet to nozzle exit cross sectional area ratio on pressure oscillation in a high-altitude test facility

The ratio of diffuser inlet to nozzle exit area (Ad/Ae) is one of the most important parameters in designing and evaluating the performance of a high-altitude test facility (HATF). Typically, at motor pressures near start or operating pressure of the diffuser, the vacuum chamber pressure oscillates at a high Ad/Ae, which may be dangerous for the HATF. In the present study, the effect of the Ad/Ae on the starting and breakdown performance of a second throat exhaust diffuser (STED) has been investigated experimentally using a HATF supplied by high-pressure cold air apparatus. Upon examining the overall performance of the diffuser, it is observed that with an increase in Ad/Ae, the oblique sealing shock wave at the diffuser inlet undergoes significant changes in location and strength. This leads to a further reduction of the flow total pressure and a stronger separation of the flow along the diffuser. Consequently, flow instability arises during starting or breakdown phases near the STED starting or operating pressure. By frequency analysis, it is observed that as Ad/Ae increases, the number of oscillatory modes of the vacuum chamber pressure increases and the dominant frequency of the oscillations in these transient situations becomes larger. Additionally, there is a significant decrease in nozzle pressure related to the onset of pressure oscillations. Furthermore, with an increase in Ad/Ae from 1.27 to 7.81, the minimum starting pressure and breakdown pressure of STED increase by 20.33% and 28.47%, respectively. Also, the hysteresis range in performance of STED vanishes at high Ad/Ae values.

Numerical simulation of the influence of wave parameters on the horizontal cylinder in the nonlinear coupling of the wave–current

Considering that marine structures are frequently subjected to the combined effects of waves and currents, the wave–current coupling environment largely determines the structural load and the surrounding water–air mixed flow, which is a typical feature of offshore structures, such as ships and offshore platforms. This study focuses on the interaction between a horizontal cylinder and a free surface in a wave–current nonlinear coupled environment using numerical simulations. The numerical method is based on the incompressible Navier–Stokes equation, with the volume-of-fluid method used to capture the free surface. Based on the second-order Stokes wave theory, we studied the impact of wave height and steepness on the cylinder force, vorticity field, free-surface deformation, and spatial–temporal distribution characteristics of entrainment bubbles. The results showed that the wave height and wave steepness have opposite effects on the root mean square (RMS) value of the force and influence the amplitude and period of the force curve. The stretching of the negative vortex led to varying degrees of double-frequency oscillation modes in the force curves. The main sources of bubbles in the wake are the breaking of the free surface and the entrainment caused by the cylinder vortex, and the bubbles caused by the former account for the majority. Compared with the wave height, an increase in wave steepness can cause a more severe interface breaking, resulting in more entrainment bubbles.

Compact inertial sensors for measuring external disturbances of physics experiments

Compact, high-precision inertial sensors are needed in the control schemes of many modern physics experiments to isolate them from disturbances caused by seismic motion. We present an inertial sensor whose mechanical oscillator fits on a one-inch diameter optic. The oscillators achieve a mechanical Quality factor of a fundamental oscillation mode of 600,000 and a resonance frequency of 50 Hz, giving them a suspension thermal noise floor lower than all commercially available inertial sensors. The motion of this fundamental mode is suitable to encode inertial motion into the sensor readout. The oscillator is combined with an optical resonator readout scheme that achieves a displacement noise of 100 fm/Hz\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\sqrt{\ ext{Hz}}$$\\end{document} above 0.2 Hz. We validate the sensors’ noise floor using a huddle test. Below 20 Hz, the sensor offers comparable performance to the best inertial sensors available today while being a fraction of the size. Above 20 Hz, the sensor is, to the author’s knowledge, the best demonstrated in the literature to date for such a sensor, with a self-noise floor of 0.1 ng/Hz\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\sqrt{\ ext{Hz}}$$\\end{document}. The excellent performance of the sensors across seismically relevant frequencies, vacuum compatibility, and compact size make it a prime candidate for integration into sophisticated seismic isolation schemes, such as those used by gravitational wave detectors.

Small signal stability analysis of a wind power system integrated with a multi-band supplementary damping controller

Abstract With the integration of wind turbines represented by doubly-fed induction generators (DFIG) into the power system, the inertia level of the system decreases sharply, and small disturbances can easily cause low-frequency oscillation (LFO) in the system. In this paper, we propose a multi-band supplementary damping controller (MBSDC) directly applied to the active power control link of the rotor-side converter in the wind power generation system. We establish a linearized model for DFIG integrated with the power system, calculate eigenvalues, and screen out multiple electromechanical oscillation modes to be suppressed. Based on the calculation results, an MBSDC is designed to suppress LFO, and we propose a parameter-setting algorithm to enhance the suppression effect. Taking the New England 10-generator 39-bus system containing DFIG as an example for simulation verification, the results show that the MBSDC and parameter setting algorithm designed in this paper can effectively suppress LFO.

Experimental study on liquid piston Stirling engine combined with self-rectifying turbine.

A liquid piston Stirling engine is an external combustion engine that uses air and water under atmospheric pressure as its working fluids. Resulting from its uncomplicated design and the capacity to operate under relatively low temperature differentials of less than 100 °C, it has attracted considerable attention in recent years. This paper presents the fundamental characteristics of the liquid piston engine combined with a self-rectifying turbine for the advancement of thermal generators. When the turbine is installed in the water region rather than in the air region, it exhibits unidirectional rotation with a rotational speed directly proportional to the velocity amplitude of the reciprocating axial flow. Additionally, the acoustic impedance within the duct section containing the turbine is determined, demonstrating that the real part of impedance rises with increasing axial velocity, indicating a loss mechanism similar to the minor loss. Furthermore, the installation of the turbine results in a breakdown of symmetry in the engine oscillation mode. To maintain symmetry and improve system design, future developments must consider the installation of a turbine in each unit. These findings can pave the way to the design of liquid piston Stirling engines and their applications in thermal energy conversion.

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

In this paper, we investigate convective modes of pseudoplastic fluid in a square cavity with solid perfectly heat-conducting boundaries under reduced gravity. The cavity performs vertical linearly polarized high-frequency vibrations. Between the vertical walls of the cavity a temperature drop is prescribed in the direction perpendicular to vibrations and the field of gravity. The fluid rheology is described using the Williamson model. The problem is solved based on the averaged thermovibration convection equations for nonlinear viscous fluids. The vibration intensity is determined by the vibration parameter V, which is proportional to the ratio of the amplitude of the vibration acceleration to the free-fall acceleration and does not depend on the temperature difference. The problem is found to have two types of solutions, which we call Newtonian and non-Newtonian modes. These solutions correspond to different convective structures, for which the dependences of the maximum of the stream function and the Nusselt number on the Grashof number are derived. We use these dependences for determining at different values of rheological parameters the threshold values of the Grasgof numbers, corresponding to the change in the modes of stationary averaged convection and the critical Grasgof numbers, corresponding to the loss of stability of stationary averaged flow and the occurrence of oscillatory modes of averaged convection. The structures of different modes of averaged stationary and oscillatory convection are studied. The results obtained for the Newtonian mode have shown that at small values of the Grasgoff number, a slow one-vortex stationary convective flow is realized in the cavity. With increasing the Grasgoff number it transforms into a three-vortex flow. A further increase in the Grasgoff number results in the appearance of a four-vortex convective flow, which eventually transforms again into a three-vortex motion. At even greater Gr, the stationary averaged motion becomes unstable and the oscillatory modes appear in the cavity. For the non-Newtonian mode, a stationary convective flow is the basic flow pattern detected in the cavity; the oscillatory modes are not observed in the whole investigated range of the Grasgoff numbers.

Anisotropic dark energy stars within vanishing complexity factor formalism: Hydrostatic equilibrium, radial oscillations, and observational implications

We investigate the structure and radial oscillations of anisotropic compact stars composed of dark energy, using the vanishing complexity factor formalism within general relativity. This novel approach establishes a direct link between the energy density and anisotropic factor, providing a robust framework for studying these exotic stellar objects. Employing an Extended Chaplygin Gas equation of state, we numerically compute interior solutions for both isotropic and anisotropic stars, revealing distinct differences in their properties. Additionally, we examine the oscillation modes and frequencies of these stars, highlighting the impact of anisotropy on their pulsational behaviour. Our results reveal distinct differences in the stellar properties, such as the metric potentials, pressure, speed of sound, and relativistic adiabatic index, between the two cases. Furthermore, we calculate the large frequency separation for the fundamental and first excited modes, offering insights relevant for future asteroseismology studies. Our findings shed light on the complex interplay of gravity, matter, and anisotropy in compact stars, providing a new perspective on dark energy's role in their structure and dynamics.

An experimental study of thermoacoustic instabilities in a RP-3 fueled multi-swirl stabilized lean direct injection combustor

An investigation of self-excited combustion instabilities in a RP-3 fueled, multi-swirl stabilized, lean direct injection (LDI) combustor was performed and investigated with different operating conditions (inlet Reynolds number and fuel–air ratio (FAR)). Pressure fluctuation and heat release data were simultaneously collected using a 3-channel synchronous acquisition system. The results show that as the FAR increases, combustion transitions from a low oscillation mode to a high oscillation mode. During the modal transition, the oscillating flame structure, formed by vortex shedding in the shear layer and the recirculation zone, is then transformed into contraction and expansion movements from the front of the central recirculation zone towards the downstream and radial directions of the combustor. When combustion is in a low oscillation mode, the driving capability of the flame is relatively low. Upon entering a high oscillation mode, the overall driving energy of the flame can maintain combustion in a limit cycle oscillation state.

An enhanced proportional resonance controller design for the PMSM based electric vehicle drive system

Permanent magnet synchronous machine (PMSM) has proven to be a more economical traction drive system for electric vehicle (EV) applications owing to increased efficiency and high-power density. However, the drive system requires more efficient control schemes to deliver better dynamic performance irrespective of dynamic changes in the motor speed, machine parameters and disturbances. Hence, to tackle the dynamic changes, to enhance the wider operating speed, to achieve precise speed tracking capability, and improved efficiency, a novel control algorithm for the PMSM based EV is proposed in this paper. The control algorithm is implemented by adopting the merits of conventional proportional resonance (PR) and proportional integral (PI) controller. The proposed control strategy is designed with an outer PI speed regulator and the inner enhanced PR (EPR) current regulator. The uniqueness of the proposed EPR controller is that the controller is designed to damp the torsional mode oscillation owing to dynamic changes such as speed and torque regulation evading the additional control loop. The effectiveness of the control scheme is tested in MATLAB Simulink and hardware-in-loop (HIL) real time simulator RT5700. To validate the effectiveness of the proposed control scheme the results are compared with the conventional control schemes. The results presented show that the proposed control technique successfully enhances the static and dynamic performance, and resilience of the EV system. Also, the proposed scheme significantly reduces the flux ripples, torque ripples, current jitter, peak overshoot, undershoot compared to the conventional current controllers.

RHEOLOGICAL PROPERTIES OF LOW-TEMPERATURE GEL-FORMING COMPOSITIONS

In this work, the rheological properties of low-temperature versions of gel-forming compositions developed at the IPC SB RAS are studied. A composition on an inorganic basis and a composition comprising a polymer and inorganic components are intended to increase oil recovery from low temperature reservoirs. The measurements were carried using a Haake Viscotester iQ rheometer with a CC25 coaxial cylinder measuring system in oscillatory mode with controlled strain at a frequency of 1 Hz. During the oscillation test, the dependences of the elastic modulus G' and the viscosity modulus G'' on the strain amplitude were recorded, the ranges of linear viscoelasticity of the samples under study were established, and the strain value for kinetic experiments was chosen. At a given strain value (0.01), the kinetic dependences of G' were obtained and the dynamics of the rheological properties and the strength of the formed structure were evaluated. The time of gelation onset has been found. It was 2300 s at 38 °C, 3900 s at 30 °C, and 21000 s at 20 °C for the inorganic composition and 200 s at 70 °C, 570 s at 60 °C, and 1920 s at 50 °C for the polymer-inorganic composition. The time required to form a gel of maximum strength also was found to be increased with increasing temperature. It has been established that for all systems under study the time of gel formation was reduced with increasing temperature, while the strength of the formed structure also increased. A slight increase in the maximum value of the modulus G' (33490 and 32805 Pa) was observed for the composition on an inorganic basis at 38 and 30 °C respectively as compared with the maximum of the dependence obtained at 20 °C (24660 Pa). The major increase in the value of the modulus G' has been recorded for the composition on a polymer-inorganic basis with increase in temperature. In this case the value of the modulus G' increased from 8970 Pa at 50 °C to 11680 Pa at 60 °C and 14590 Pa at 70 °C, while the general form of the rheokinetic dependence for different compositions and temperatures remained the same. For citation: Kozhevnikov I.S., Bogoslovskii А.V. Rheological properties of low-temperature gel-forming compositions. ChemChemTech [Izv. Vyssh. Uchebn. Zaved. Khim. Khim. Tekhnol.]. 2024. V. 67. N 8. P. 22-28. DOI: 10.6060/ivkkt.20246708.1t.

Linear operator theory of phase mixing

ABSTRACT We study solutions of the collisionless Boltzmann equation (CBE) in a functional Koopman representation. This facilitates the use of linear spectral techniques characteristic of the analysis of Schrödinger-type equations. For illustrative purposes, we consider the classical phase mixing of a non-interacting distribution function in a quartic potential. Solutions are determined perturbatively relative to a harmonic oscillator. We impose a form of coarse-graining by choosing a finite-dimensional basis to represent the distribution function and time evolution operators, which sets a minimum length-scale on phase space structure. We observe a relationship between the dimension of the representation and the multiplicity of the harmonic oscillator eigenvalues. System dynamics are understood in terms of degenerate subspaces of the linear operator spectra. Each subspace is associated with a mode of the harmonic oscillator, the first two being bending and breathing structures. The quartic potential splits the degenerate eigenvalues within each subspace. This facilitates the formation of spiral structure as deformations from the harmonic oscillator modes. We ultimately argue that this construction provides a promising avenue for study of self-interacting systems experiencing phase mixing, which is an outstanding problem in the context of the Gaia DR2 vertical phase space spirals.

A Comparative Analysis between FFT, EMD, and EEMD for Epilepsy Detection

Context: Epilepsy is a neurological disease that affects more than 50 million people worldwide, causing recurrent seizures, with a significant impact on patients' quality of life due to abnormally synchronized neuronal activity. Method: This article discusses three methods used for signal analysis in patients diagnosed with epilepsy. Conventional signal decomposition methods, such as the fast Fourier transform, widely used in signal analysis based on time series techniques, have some issues when analyzing nonlinear and non-stationary signals, in addition to difficulties in detecting low-order frequencies. Results: To overcome these limitations, alternatives such as empirical mode decomposition and one of its variants, called ensemble empirical mode decomposition, have been developed. These techniques allow observing different oscillation modes through intrinsic mode functions and instantaneous frequencies. Conclusions: In this study, the results obtained through the aforementioned techniques were compared, revealing the impact of nonlinear methods on the reconstruction of brain activity.

Weak time-scale separation at the onset of oscillatory magnetoconvection in rapidly rotating fluids

Convective instabilities are one of the integral parts of the dynamics of flows driven by thermal buoyancy. Naturally, physical phenomena exhibit a wide disparity in the length and timescales of the field variables in numerical simulations and experimental observations. Such variations are not represented in the traditional normal mode stability analysis attempting to understand the onset of convection. This study attempts to incorporate different time constants for different perturbation variables in the linear stability analysis with the help of a Taylor series expansion. The infinite horizontal layer model is chosen for simplicity. Apart from the classical Rayleigh-Bénard system, additional physical effects such as background rotation and magnetic field have been considered with plausible implications for geophysical flow applications. The time scale separation is implemented by considering a slight change in the frequency of temperature perturbation compared to that for other physical quantities. Both analytical and numerical methods have been utilised for the investigation. The threshold buoyancy force is reduced when the temperature perturbation has a smaller frequency than the frequencies of other variables. Besides that, the onset wavenumber and frequency of the oscillatory modes are modified due to weak scale separation from the onset characteristics of the reference case. In particular, enhanced frequency of temperature perturbations leads to smaller-scaled magnetically controlled convective rolls and larger-scaled viscously controlled instabilities at the onset. A robust dependence of the onset characteristics with the parameter quantifying the timescale separation is obtained. Additionally, two transitions in the convective onset modes with scale separation have been identified.

Mechanism of vortex oscillation around a hemisphere–cylinder body

Previous studies have shown that low-frequency vortex oscillations occur around a hemisphere–cylinder body at different angles of attack, but the underlying mechanism is still unclear. In this study, we examine the origin of the vortex oscillation using numerical simulations and global linear stability analysis. The vortex oscillation is reproduced using numerical simulations, and the oscillatory modes are computed through dynamic mode decomposition (DMD). We obtain the base flow through a selective frequency damping method, which exhibits a pair of steady leeward vortices over the body. The four unstable modes are computed using a modified Arnoldi iteration. The antisymmetric mode with a Strouhal number of 0.105 is discovered to be responsible for the alternate oscillation of the vortex pair, and the mode with a Strouhal number of 0.220 corresponds to the in-phase vortex oscillation. Their frequencies have good agreement with the modes of DMD. The other two unstable modes with higher frequencies, one antisymmetric and one symmetric, are harmonic frequencies of the above two modes. The study conclusively verifies that the vortex oscillation over a hemisphere–cylinder body originates from a global flow instability.

Analysis of hovering flight stability of an insect-like flapping-wing robot in Martian condition

In efforts to search for life-form on a planet other than Earth, many projects have been launched to measure and analyze the surface and atmospheric conditions on Mars. To extend the exploration of Mars, scientists have attempted to perform flight in the extremely thin air environment of Mars. Inspired by the successful controlled flight of NASA's Helicopter on Mars, in this work, we investigate the six-degrees-of-freedom hovering flight dynamics of an insect-like flapping-wing robot, called KUBeetle. To hover the KUBeetle on Mars, using the scaling ratios of air density and gravity on Earth and Mars, we have found that the flapping frequency should be about 113.48 Hz, which is 4.93 times faster than that on Earth, and the required lift on Mars is about 38% of that on Earth. The computational fluid dynamics (CFD) analysis suggests that the lift produced by the flapping wings at the frequency on Mars is close to that predicted by the scaling ratios. A series of simulations are used to compute the stability derivatives to complete the equations of motion. The eigenvalue and eigenvector analyses have identified one fast subsidence mode, one slow subsidence mode, and one divergence oscillation mode in the longitudinal and lateral motions on Mars, which are identical to those on Earth. The dynamic responses of the robot under three initial disturbances in the translational speeds indicate that due to the 94% smaller real part of the complex roots in the divergence oscillation modes of the longitudinal and lateral motions on Mars, the growth of flight disturbances is much slower than on Earth. Therefore, the KUBeetle is capable of hovering flight on Mars, since it has enough time to use the feedback controls to stabilize the system. This work can be used to support the studies of the flight control of flapping wing robots in the Martian atmospheric condition, and to generalize such flight control for other non-terrestrial applications in the atmospheres of other planets and/or their moons.