Weakly Nonlinear Rayleigh–Taylor Instability in Cylindrically Convergent Geometry

The Rayleigh–Taylor instability(RTI) in cylindrical geometry is investigated analytically through a second-order weakly nonlinear(WN) theory considering the Bell–Plesset(BP) effect. The governing equations for the combined perturbation growth are derived. The WN solutions for an exponentially convergent cylinder are obtained. It is found that the BP and RTI growths are strongly coupled, which results in the bubble-spike asymmetric structure in the WN stage. The large Atwood number leads to the large deformation of the convergent interface. The amplitude of the spike grows faster than that of the bubble especially for large mode number m and large Atwood number A. The averaged interface radius is small for large mode number perturbation due to the mode-coupling effect.     

Rayleigh–Taylor instability at spherical interfaces of incompressible fluids

Rayleigh–Taylor instability(RTI) of three incompressible fluids with two interfaces in spherical geometry is derived analytically. The growth rate on the two interfaces and the perturbation feedthrough coefficients between two spherical interfaces are derived. For low-mode perturbation, the feedthrough effect from outer interface to inner interface is much more severe than the corresponding planar case, while the feedback from inner interface to the outer interface is smaller than that in planar geometry. The low-mode perturbations lead to the pronounced RTI growth on the inner interface of a spherical shell that are larger than the cylindrical and planar results. It is the low-mode perturbation that results in the difference between the RTI growth in spherical and cylindrical geometry. When the mode number of the perturbation is large enough, the results in cylindrical geometry are recovered.     

Cylindrical effects in weakly nonlinear Rayleigh Taylor instability

The classical Rayleigh–Taylor instability(RTI)at the interface between two variable density fluids in the cylindrical geometry is explicitly investigated by the formal perturbation method up to the second order.Two styles of RTI,convergent(i.e.,gravity pointing inward)and divergent(i.e.,gravity pointing outwards)configurations,compared with RTI in Cartesian geometry,are taken into account.Our explicit results show that the interface function in the cylindrical geometry consists of two parts:oscillatory part similar to the result of the Cartesian geometry,and non-oscillatory one contributing nothing to the result of the Cartesian geometry.The velocity resulting only from the non-oscillatory term is followed with interest in this paper.It is found that both the convergent and the divergent configurations have the same zeroth-order velocity,whose magnitude increases with the Atwood number,while decreases with the initial radius of the interface or mode number.The occurrence of non-oscillation terms is an essential character of the RTI in the cylindrical geometry different from Cartesian one.     

Nonlinear saturation amplitude of cylindrical Rayleigh Taylor instability

The nonlinear saturation amplitude （NSA） of the fundamental mode in the classical Rayleigh-Taylor instability with a cylindrical geometry for an arbitrary Atwood number is analytically investigated by considering the nonlinear corrections up to the third order. The analytic results indicate that the effects of the initial radius of the interface （r0） and the Atwood number （A） play an important role in the NSA of the fundamental mode. The NSA of the fundamental mode first increases gently and then decreases quickly with increasing A. For a given A, the smaller the ro/λ（λ is the perturbation wavelength）, the larger the NSA of the fundamental mode. When ro/λ is large enough （r0 〉〉 λ）, the NSA of the fundamental mode is reduced to the prediction in the previous literatures within the framework of the third-order perturbation theory.     

On the Second Harmonic Generation through Bell-Plesset Effects in Cylindrical Geometry

Generation of the second harmonic initiated by Bell-Plesset effects in a cylindrical geometry is studied analytically. For an initial single-mode velocity perturSation, the second-order mode-coupling formula is obtained by expanding the perturbation displacement and velocity potential up to the second-order accuracy. It is found that the initially symmetric interface evolves into a significant bubble-spike asymmetric pattern. The second-order solutions clearly show that the amplitude of the spike grows faster than that of the buSble. The temporal evolutions of the amplitudes of the 5ubSie and spike are dependent on the interface velocity Vo. The larger interface velocity leads to the smaller amplitude of the perturbation at an arbitrary interface position in a cylindrically convergent geometry.     

Nonlinear Saturation Amplitude in Classical Planar Richtmyer–Meshkov Instability

The classical planar Richtmyer–Meshkov instability(RMI) at a fluid interface supported by a constant pressure is investigated by a formal perturbation expansion up to the third order,and then according to definition of nonlinear saturation amplitude(NSA) in Rayleigh–Taylor instability(RTI),the NSA in planar RMI is obtained explicitly.It is found that the NSA in planar RMI is affected by the initial perturbation wavelength and the initial amplitude of the interface,while the effect of the initial amplitude of the interface on the NSA is less than that of the initial perturbation wavelength.Without marginal influence of the initial amplitude,the NSA increases linearly with wavelength.The NSA normalized by the wavelength in planar RMI is about 0.11,larger than that corresponding to RTI.     

Interface Width Effect on the Weakly Nonlinear Rayleigh–Taylor Instability in Spherical Geometry

Interface width effect on the spherical Rayleigh–Taylor instability in the weakly nonlinear regime is studied by numerical simulations.For Legendre perturbation mode P_n with wave number K_n and interface half-width L,the commonly adopted empirical linear growth rate formula γ_n~(em)(L)=γn/√1 + K_nL is found to be sufficient in spherical geometry.At the weakly nonlinear stage,the interface width affects the mode coupling processes.The development of the mode P_(2n) is substantially influenced by the interface width.Moreover,the nonlinear saturation amplitude increases with the interface width.     

Molecular dynamics simulation of cylindrical Richtmyer-Meshkov instability

The microscopic-scale Richtmyer-Meshkov(RM) instability of a single-mode Cu-He interface subjected to a cylindrically converging shock is studied through the classical molecular dynamics simulation. An unperturbed interface is first considered to examine the flow features in the convergent geometry, and notable distortions at the circular inhomogeneity are observed due to the atomic fluctuation. Detailed processes of the shock propagation and interface deformation for the single-mode interface impacted by a converging shock are clearly captured. Different from the macroscopic-scale situation, the intense molecular thermal motions in the present microscale flow introduce massive small wavelength perturbations at the single-mode interface, which later significantly impede the formation of the roll-up structure. Influences of the initial conditions including the initial amplitude,wave number and density ratio on the instability growth are carefully analyzed. It is found that the late-stage instability development for interfaces with a large perturbation does not depend on its initial amplitude any more. Surprisingly, as the wave number increases from 8 to 12, the growth rate after the reshock drops gradually. The distinct behaviors induced by the amplitude and wave number increments indicate that the present microscopic RM instability cannot be simply characterized by the amplitude over wavelength ratio(η). The pressure history at the convergence center shows that the first pressure peak caused by the shock focusing is insensitive to η, while the second one depends heavily on it.     

A Weakly Nonlinear Model for Kelvin-Helmholtz Instability in Incompressible Fluids

A weakly nonlinear model is proposed for the Kelvin-Helmholtz instability in two-dimensional incompressible fluids by expanding the perturbation velocity potential to third order. The third-order harmonic generation effects of single-mode perturbation are analyzed, as well as the nonlinear correction to the exponential growth of the fundamental modulation. The weakly nonlinear results are supported by numerical simulations. Density and resonance effects exist in the development of mode coupling.     

Two-Dimensional Rayleigh-Taylor Instability in Incompressible Fluids at Arbitrary Atwood Numbers

The Rayleigh-Taylor instability in two-dimensional incompressible fluids at arbitrary Atwood numbers is studied by expanding the perturbation velocity potential to third order. The second and third harmonic generation effects of single-mode perturbation are analyzed, as well as the nonlinear correction to the exponential growth of the fundamental modulation. The mode coupling coefficients are dependent on the Atwood numbers. Our simulations support the weakly nonlinear results. We find that the ratio of the nonlinear saturation amplitude ηs and the perturbation wavelength λ is dependent on the Atwood number AT and the relation is ηs/λ=（1/π）[√2/5/√（1＋3AT2 ）].     

Influence of isothermal approximation on the phase-field simulation of directional growth in undercooled melt

By using the phase-field approach,we have simulated the directional growth of alloys in undercooled moten states under the isothermal and nonisothermal conditions.The influences of the isothermal approximation on simulation results are discussed.We found that for undercooling greater than 25K,the isothermal approximation overestimates the interface growth velocity and reduces a critical velocity for an absolute stable planar interface,thus in this simulation,the uinterface morphology shows the plane-cell-plane transition with increasing initial undercooling of the mele,and the planar interface obtained under a large undercooling is absolutely stable.Whereas in the nonisothermal simulation,only plane-cell transition occures in the same range of the initial undercoolings of the melt,and the planar interface tends to be destabilized and evolve into cells.     

Interface coupling effects of weakly nonlinear Rayleigh–Taylor instability with double interfaces

Taking the Rayleigh–Taylor instability with double interfaces as the research object,the interface coupling effects in the weakly nonlinear regime are studied numerically.The variation of Atwood numbers on the two interfaces and the variation of the thickness between them are taken into consideration.It is shown that,when the Atwood number on the lower interface is small,the amplitude of perturbation growth on the lower interface is positively related with the Atwood number on the upper interface.However,it is negatively related when the Atwood number on the lower interface is large.The above phenomenon is quantitatively studied using an analytical formula and the underlying physical mechanism is presented.     

The numerical study of shock-induced hydrodynamic instability and mixing

Based on multi-fluid volume fraction and piecewise parabolic method (PPM), a multi-viscosity-fluid hydrodynamic code MVPPM (Multi-Viscosity-Fluid Piecewise Parabolic Method) is developed and applied to the problems of shock-induced hydrodynamic interfacial instability and mixing. Simulations of gas/liquid interface instability show that the influences of initial perturbations on the fluid mixing zone (FMZ) growth are significant, especially at the late stages, while grids have only a slight effect on the FMZ width, when the interface is impulsively accelerated by a shock wave passing through it. A numerical study of the hydrodynamic interfacial instability and mixing of gaseous flows impacted by re-shocks is presented. It reveals that the numerical results are in good agreement with the experimental results and the mixing growth rate strongly depends on initial conditions. Ultimately, the jelly layer experiment relevant to the instability impacted by exploding is simulated. The shape of jelly interface, position of front face of jelly layer, crest and trough of perturbation versus time are given; their simulated results are in good agreement with experimental results.     

Competition of the multiple Grtler modes in hypersonic boundary layer flows

Competition of multiple Grtler modes in hypersonic boundary layer flows are investigated with the local and marching methods.The wall-layer mode(mode W)and the trapped-layer mode(mode T)both occur in the compressible boundary layer where there exists a temperature adjustment layer near the upper edge.The mode T has the largest growth rate at a lower Grtler number while the mode W dominates at larger Grtler numbers.These two modes are both responsible for the flow transition in the hypersonic flows especially when Grtler number is in the high value range in which the crossover of these two modes takes place.Such high Grtler numbers are virtually far beyond the neutral regime.The nonparallel base flows,therefore,cease to influence the stability behavior of the Grtler modes.The effects of the Mach number on the multiple Grtler modes are studied within a chosen Mach number of 0.95,2,4 and 6.When the flow Mach number is sufficiently large,e.g.,Ma 4,the growth rate crossover of the mode T and mode W occurs both in the conventional G-βmap as well as on the route downstream for a fixed wavelength disturbance.Four particular regions(Region T,T-W,W-T and W)around the crossover point are highlighted with the marching analysis and the result matches that of the local analysis.The initial disturbance of a normal mode maintains the shape in its corresponding dominating region while a shape-transformation occurs outside this region.     

Physiological Flow of Jeffrey Six Constant Fluid Model due to Ciliary Motion

The main purpose of this article is to present a mathematical model of ciliary motion in an annulus. In this analysis, the peristaltic motion of non-Newtonian Jeffrey six constant fluid is observed in an annulus with ciliated tips in the presence of heat and mass transfer. The effects of viscous dissipation are also considered. The flow equations of non-Newtonian fluid for the two-dimensional tube in cylindrical coordinates are simplified using the low Reynolds number and long wave-length approximations. The main equations for Jeffrey six constant fluid are considered in cylindrical coordinates system. The resulting nonlinear problem is solved using the regular perturbation technique in terms of a variant of small dimensionless parameter α. The results of the solutions for velocity, temperature and concentration field are presented graphically. B_k is Brinkman number, ST is soret number, and SH is the Schmidth number. Outcome for the longitudinal velocity, pressure rise, pressure gradient and stream lines are represented through graphs. In the history, the viscous-dissipation effect is usually represented by the Brinkman number.     

Role of Viscosity Stratification and Insoluble Surfactant in Instability of Two-Layer Channel Flow

We study, for the case of the two layer plane Poiseuille flow, the effect of viscosity stratification and interracial surfactant on the flow instability. Considering a normal mode of the streamwise wave number α, both the linear and energy analyses are presented. The expressions of perturbation energy supplied at the interface are derived. The result demonstrates that the jumps of horizontal velocity and tangential stress of the perturbed flow across the interface could be induced by the presence of viscosity stratification and surfactant. This is expected to be responsible for the Yih and Marangoni instability.     

Toroidicity and shape dependence of peeling mode growth rates in axisymmetric toroidal plasmas

The growth rate of the peeling mode instability with large toroidal mode number is calculated for general axisymmetric toroidal plasmas, including tokamaks and the spherical torus(ST) equilibia by using formalism presented by Connor et al.Analytic equilibia with non-zero edge current density and quasi-uniform current profiles are assumed. It is found that in sharp D-shape tokamak plasma, the derivative of the safety factor with respect to the poloidal flux becomes very large,making the perturbed poloidal motion very large, in turn making a significant reduction of the growth rate of the peeling mode, similar to the X-point effect in diverted plasma. The large aspect ratio effect is also studied, which reduces the growth rate further.     

Effect of far-field flow on a columnar crystal in the convective undercooled melt

The growth behavior of a columnar crystal in the convective undercooled melt affected by the far-field uniform flow is studied and the asymptotic solution for the interface evolution of the columnar crystal is derived by means of the asymptotic expansion method.The results obtained reveal that the far-field flow induces a significant change of the temperature around the columnar crystal and the convective flow caused by the far-field flow accelerates the growth velocity of the interface of the growing columnar crystal in the upstream direction and inhibits its growth velocity in the downstream direction.Our results are similar to the experimental data and numerical simulations.     

Compressibility Effects on the Rayleigh--Taylor Instability Growth Rates

Effects of two compressibility parameters, i.e. the ratio of specific heats and the equilibrium pressure at the interface, on the Rayleigh-Taylor instability （RTI） growth rates are studied under the same initial conditions, which include the mass, pressure profile, and density profile of the two superposed fluids. The results obtained reconcile the stabilizing and destabilizing effects of compressibility reported in the literature. The influences of the ratio of specific heats on the RTI growth rates are not only stabilized but also destabilized. The effects of the equilibrium pressure at the interface on the growth rates are destabilized.     

Polar interface and surface optical vibration spectra in multi-layer wurtzite quantum wires： transfer matrix method

The polar interface optical （IO） and surface optical （SO） phonon modes and the corresponding Froehlich electron phonon-interaction Hamiltonian in a freestanding multi-layer wurtzite cylindrical quantum wire （QWR） are derived and studied by employing the transfer matrix method in the dielectric continuum approximation and Loudon＇s uniaxial crystal model. A numerical calculation of a freestanding wurtzite GaN/AlN QWR is performed. The results reveal that for a relatively large azimuthal quantum number m or wave-number kz in the free z-direction, there exist two branches of IO phonon modes localized at the interface, and only one branch of SO mode localized at the surface in the system. The degenerating behaviours of the IO and SO phonon modes in the wurtzite QWR have also been clearly observed for a small kz or m. The limiting frequency properties of the IO and SO modes for large kz and m have been explained reasonably from the mathematical and physical viewpoints. The calculations of electron-phonon coupling functions show that the high-frequency IO phonon branch and SO mode play a more important role in the electron phonon interaction.