No Arabic abstract
The effects of compressibility on Rayleigh-Taylor instability (RTI) are investigated by inspecting the interplay between thermodynamic and hydrodynamic non-equilibrium phenomena (TNE, HNE, respectively) via a discrete Boltzmann model (DBM). Two effective approaches are presented, one tracking the evolution of the emph{local} TNE effects and the other focussing on the evolution of the mean temperature of the fluid, to track the complex interfaces separating the bubble and the spike regions of the flow. It is found that, both the compressibility effects and the emph{global} TNE intensity show opposite trends in the initial and the later stages of the RTI. Compressibility delays the initial stage of RTI and accelerates the later stage. Meanwhile, the TNE characteristics are generally enhanced by the compressibility, especially in the later stage. The global or mean thermodynamic non-equilibrium indicators provide physical criteria to discriminate between the two stages of the RTI.
Rayleigh-Taylor (RT) instability widely exists in nature and engineering fields. How to better understand the physical mechanism of RT instability is of great theoretical significance and practical value. At present, abundant results of RT instability have been obtained by traditional macroscopic methods. However, research on the thermodynamic non-equilibrium (TNE) effects in the process of system evolution is relatively scarce. In this paper, the discrete Boltzmann method based on non-equilibrium statistical physics is utilized to study the effects of the specific heat ratio on compressible RT instability. The evolution process of the compressible RT system with different specific heat ratios can be analyzed by the temperature gradient and the proportion of the non-equilibrium region. Firstly, as a result of the competition between the macroscopic magnitude gradient and the non-equilibrium region, the average TNE intensity first increases and then reduces, and it increases with the specific heat ratio decreasing; the specific heat ratio has the same effect on the global strength of the viscous stress tensor. Secondly, the moment when the total temperature gradient in y direction deviates from the fixed value can be regarded as a physical criterion for judging the formation of the vortex structure. Thirdly, under the competition between the temperature gradients and the contact area of the two fluids, the average intensity of the non-equilibrium quantity related to the heat flux shows diversity, and the influence of the specific heat ratio is also quite remarkable.
In this paper, the coupled Rayleigh-Taylor-Kelvin-Helmholtz instability(RTI, KHI and RTKHI, respectively) system is investigated using a multiple-relaxation-time discrete Boltzmann model. Both the morphological boundary length and thermodynamic nonequilibrium (TNE) strength are introduced to probe the complex configurations and kinetic processes. In the simulations, RTI always plays a major role in the later stage, while the main mechanism in the early stage depends on the comparison of buoyancy and shear strength. It is found that, both the total boundary length $L$ of the condensed temperature field and the mean heat flux strength $D_{3,1}$ can be used to measure the ratio of buoyancy to shear strength, and to quantitatively judge the main mechanism in the early stage of the RTKHI system. Specifically, when KHI (RTI) dominates, $L^{KHI} > L^{RTI}$ ($L^{KHI} < L^{RTI}$), $D_{3,1}^{KHI} > D_{3,1}^{RTI}$ ($D_{3,1}^{KHI} < D_{3,1}^{RTI}$); when KHI and RTI are balanced, $L^{KHI} = L^{RTI}$, $D_{3,1}^{KHI} = D_{3,1}^{RTI}$. A second sets of findings are as below: For the case where the KHI dominates at earlier time and the RTI dominates at later time, the evolution process can be roughly divided into two stages. Before the transition point of the two stages, $L^{RTKHI}$ initially increases exponentially, and then increases linearly. Hence, the ending point of linear increasing $L^{RTKHI}$ can work as a geometric criterion for discriminating the two stages. The TNE quantity, heat flux strength $D_{3,1}^{RTKHI}$, shows similar behavior. Therefore, the ending point of linear increasing $D_{3,1}^{RTKHI}$ can work as a physical criterion for discriminating the two stages.
The dynamics of a thin liquid film on the underside of a curved cylindrical substrate is studied. The evolution of the liquid layer is investigated as the film thickness and the radius of curvature of the substrate are varied. A dimensionless parameter (a modified Bond number) that incorporates both geometric parameters, gravity, and surface tension is identified, and allows the observations to be classified according to three different flow regimes: stable films, films with transient growth of perturbations followed by decay, and unstable films. Experiments and theory confirm that, below a critical value of the Bond number, curvature of the substrate suppresses the Rayleigh-Taylor instability.
This work examines the effect of the embedded magnetic field strength on the non-linear development of the magnetic Rayleigh-Taylor Instability (RTI) (with a field-aligned interface) in an ideal gas close to the incompressible limit in three dimensions. Numerical experiments are conducted in a domain sufficiently large so as to allow the predicted critical modes to develop in a physically realistic manner. The ratio between gravity, which drives the instability in this case (as well as in several of the corresponding observations), and magnetic field strength is taken up to a ratio which accurately reflects that of observed astrophysical plasma, in order to allow comparison between the results of the simulations and the observational data which served as inspiration for this work. This study finds reduced non-linear growth of the rising bubbles of the RTI for stronger magnetic fields, and that this is directly due to the change in magnetic field strength, rather than the indirect effect of altering characteristic length scales with respect to domain size. By examining the growth of the falling spikes, the growth rate appears to be enhanced for the strongest magnetic field strengths, suggesting that rather than affecting the development of the system as a whole, increased magnetic field strengths in fact introduce an asymmetry to the system. Further investigation of this effect also revealed that the greater this asymmetry, the less efficiently the gravitational energy is released. By better understanding the under-studied regime of such a major phenomenon in astrophysics, deeper explanations for observations may be sought, and this work illustrates that the strength of magnetic fields in astrophysical plasmas influences observed RTI in subtle and complex ways.
We present a generic coarse-grained model to describe molecular motors acting on polymer substrates, mimicking, for example, RNA polymerase on DNA or kinesin on microtubules. The polymer is modeled as a connected chain of beads; motors are represented as freely diffusing beads which, upon encountering the substrate, bind to it through a short-ranged attractive potential. When bound, motors and polymer beads experience an equal and opposite active force, directed tangential to the polymer; this leads to motion of the motors along the polymer contour. The inclusion of explicit motors differentiates our model from other recent active polymer models. We study, by means of Langevin dynamics simulations, the effect of the motor activity on both the conformational and dynamical properties of the substrate. We find that activity leads, in addition to the expected enhancement of polymer diffusion, to an effective reduction of its persistence length. We discover that this effective softening is a consequence of the emergence of double-folded branches, or hairpins, and that it can be tuned by changing the number of motors or the force they generate. Finally, we investigate the effect of the motors on the probability of knot formation. Counter-intuitively our simulations reveal that, even though at equilibrium a more flexible substrate would show an increased knotting probability, motor activity leads to a marked decrease in the occurrence of knotted conformations with respect to equilibrium.