On Nonlinear Evolution Approach for Convective Flow during Alloy Solidification

We consider the problem of nonlinear buoyant flow in a horizontal mushy layer during alloy solidification. We study the nonlinear evolution of such flow based on a recently developed realistic model for the mushy layer. The evolution approach is based on a Landau type equation for the amplitude of the secondary nonlinear solution, which is derived in this article. Using both analytical and computational methods, we calculate the solution to the evolution equation for both subcritical and supercritical regimes. We find, in particular, that for a passive mush, where the permeability is constant, and supercritical regime, the primary solution is linearly unstable to the secondary solution which becomes a steady stable solution for sufficiently large time, while the secondary solution decays to zero for the subcritical regime. On the other hand, for a realistic reactive mush, where the permeability is variable, the secondary flow can break down in a finite time for either supercritical or subcritical regime, which indicates existence of some kind of bursting behavior. These results are then compared to the corresponding ones based on the weakly nonlinear theory.     

On Three-Dimensional Non-linear Buoyant Convection in Ternary Solidification

We consider the problem of three-dimensional non-linear buoyant convection in ternary solidification. Under the limit of large far-field temperature, the convective flow is modeled to be in a rectangular cube composing of a horizontal liquid layer above a primary mushy layer, which itself is over a secondary mushy layer. We first apply linear stability analysis to calculate the conditions at the onset of motion. Next, we carry out weakly non-linear analyses to determine solutions in the form of hexagons and their possible stability and to obtain information about tendency for chimney formation. We find that if the flow is driven either from both mushy layers with equal critical conditions at the onset of motion or only by the primary mushy layer, then the flow can be in the form of a double-cell structure vertically with down-hexagons below or above up-hexagons. There is tendency for vertically oriented chimney formation at different horizontal locations in each mushy layer. For the cases where only the critical conditions at the onset of motion are equal in both mushy layers and depending on the values of the mush Rayleigh numbers, the flow can be subcritical (or supercritical) in both mushy layers or mixed subcritical in one layer and supercritical in another layer.     

On non-linear magneto convection in a mushy layer with joule heating

This paper studies the flow pattern of non-linear magneto convection that can be realized in a horizontal mushy layer and in the presence of joule heating, which is the amount of heat produced by the induced magnetic field. We consider the appropriate system of equations and the associated boundary conditions for the flow in the mushy layer subjected to a vertical magnetic field of uniform strength. Under certain assumptions and conditions, we determine the stable finite-amplitude solutions of the resulting system using a perturbation approach and stability analysis. We find, in particular, that for a wide range of values of the joule heating parameter and sufficiently small amplitude of the flow, the only stable convective flow is in the form of subcritical down-hexagons with down-flow at the cells' centers and up-flow at the cells' boundaries. This result is in sharp contrast to the case in the absence of joule heating where instead the subcritical up-hexagons with up-flow at the cells' centers and down-flow at the cells' boundaries can be stable. In the presence of joule heating the stable subcritical down-hexagons were found to be enhanced with increasing the strength of the externally imposed magnetic field.     

On Mixed Oscillatory and Steady Modes of Nonlinear Convection in Mushy Layers

We consider the problem of mixed oscillatory and steady modes of nonlinear compositional convection in horizontal mushy layers during the solidification of binary alloys. Under a near-eutectic approximation and the limit of large far-field temperature, we determine a number of two- and three-dimensional weakly nonlinear mixed solutions, and the stability of these solutions with respect to arbitrary three-dimensional disturbances is then investigated. The present investigation is an extension of the problem of mixed oscillatory and steady modes of convection, which was investigated by Riahi (J Fluid Mech 517: 71–101, 2004), where some calculated results were inaccurate due to the presence of a singular point in the equation for the linear frequency. Here we resolve the problem and find some significant new results. In particular, over a wide range of the parameter values, we find that the properties of the preferred and stable solution in the form of particular subcritical mixed standing and steady hexagons appeared to be now in much better agreement with the available experimental results (Tai et al., Nature 359:406–408, 1992) than the one reported in Riahi (J Fluid Mech 517:71–101, 2004). We also determined a number of new types of preferred supercritical solutions, which can be preferred over particular values of the parameters and at relatively higher values of the amplitude of convection.     

Effect of permeability on the instability of a non-Newtonian film down a porous inclined plane

A thin film of a power–law fluid flowing down a porous inclined plane is considered. It is assumed that the flow through the porous medium is governed by the modified Darcy’s law together with Beavers–Joseph boundary condition for a general power–law fluid. Under the assumption of small permeability relative to the thickness of the overlying fluid layer, the flow is decoupled from the filtration flow through the porous medium and a slip condition at the bottom is used to incorporate the effects of the permeability of the porous substrate. Applying the long-wave theory, a nonlinear evolution equation for the thickness of the film is obtained. A linear stability analysis of the base flow is performed and the critical condition for the onset of instability is obtained. The results show that the substrate porosity in general destabilizes the film flow system and the shear-thinning rheology enhances this destabilizing effect. A weakly nonlinear stability analysis reveals the existence of supercritical stable and subcritical unstable regions in the wave number versus Reynolds number parameter space. The numerical solution of the nonlinear evolution equation in a periodic domain shows that the fully developed nonlinear solutions are either time-dependent modes that oscillate slightly in the amplitude or time independent stable two-dimensional nonlinear waves with large amplitude referred to as ‘permanent waves’. The results show that the shape and the amplitude of the nonlinear waves are strongly influenced by the permeability of the porous medium and the shear-thinning rheology.     

Control of cross-flow-induced vibrations of square cylinders using linear and nonlinear delayed feedbacks

We investigate the effectiveness of linear and nonlinear time-delay feedback controls to suppress high amplitude oscillations of an elastically mounted square cylinder undergoing galloping oscillations. A representative model that couples the transverse displacement and the aerodynamic force is used. The quasi-steady approximation is used to model the galloping force. A linear analysis is performed to investigate the effect of linear time-delay controls on the onset speed of galloping and natural frequencies. It is demonstrated that a linear time-delay control can be used to delay the onset speed of galloping. The normal form of the Hopf bifurcation is then derived to characterize the type of the instability (supercritical or subcritical) and to determine the effects of the linear and nonlinear time-delay parameters on their outputs near the bifurcation. The results show that the nonlinear time-delay control can be efficiently implemented to significantly reduce the galloping amplitude and suppress any dangerous behavior by converting any subcritical Hopf bifurcation into a supercritical one.     

Nonlinear Analysis of a Cantilever Pipe Containing Pulsatile Flow

In this paper we investigate the nonlinear dynamics of a cantilever elastic pipe that contains pulsatile flow. The equation of motion was derived by using Hamiltonian action function. We use Galerkin's technique to include only finite number of spatial modes in the solution.The stability chart of the time-varying system was computed in the space of the relative perturbation amplitude of the flow velocity and dimensionless forcing frequency using an efficient numerical method based on Chebyshev polynomials. In the near of some critical regions bifurcation diagrams were also computed which show secondary Hopf bifurcations and phase locking followed by chaotic motion.     

三维方腔介电液体电对流的数值模拟研究

本文对封闭方腔内介电液体电对流进行了三维数值模拟研究.方腔的6个边界为固壁;4个侧边界为电绝缘边界;上下界面为两个电极.直流电场作用在从底部电极注入的自由电荷上,从而对液体施加库伦体积力并驱动流体流动形成电对流.为了求解这一物理问题,发展了一种二阶精度的有限体积法来求解完整的控制方程,包括Navier-Stokes方程和一组简化的Maxwell方程.考虑到电荷密度方程的强对流占优特性,采用了全逆差递减格式来求解该方程,获得了准确有界的解.通过研究发现,该流动在有限振幅区内的分叉类型为亚临界,即系统存在一个线性和非线性临界值,分别对应流动的开始和终止.由于非线性临界值比线性值小,因此两个临界值之间有一个迟滞回线.与无限大域中的自由对流相比,侧壁施加的额外约束改变了流场结构,使这两个临界值均有所增大.此外,还讨论了电荷密度和速度场的空间分布特征,发现电荷密度分布中存在电荷空白区.最后对更小空间尺寸情况计算结果表明,流动的线性分叉类型为超临界.本文的结果拓展了已有的二维有限空间内电对流的研究,并为三维电对流的线性和弱非线性理论分析提供参考.      

Nonlinear dynamics of extensible fluid-conveying pipes,supported at both ends

In this paper, the post-divergence behaviour of extensible fluid-conveying pipes supported at both ends is studied using the weakly nonlinear equations of motion of Semler, Li and Païdoussis. The two coupled nonlinear partial differential equations are discretized via Galerkin's method and the resulting set of ordinary differential equations is solved either by Houbolt's finite difference method or via AUTO. Typically, the pipe is stable at its original static equilibrium position up to the flow velocity where it loses stability by static divergence via a supercritical pitchfork bifurcation. The amplitude of the resultant buckling increases with increasing flow, but no secondary instability occurs beyond the pitchfork bifurcation. The effects of the system parameters on pipe behaviour as well as the possibility of a subcritical pitchfork bifurcation have also been studied.     

激励Stuart-Landau方程的研究--周期解、稳定性及流动控制

解析得出了有外部激励的Stuart-Landau（S-L）方程的频率锁定周期解,对这些解与外部激励振幅和频率的依赖关系做了详细研究,并用周期系统稳定性理论确定了解的稳定性边界.还对S-L方程所描述的流动控制效果进行了研究,发现由于外部激励的作用,稳定的锁频解可能比原来的饱和解能量减少了,外部的控制最多能使扰动能量减少为原来的一半.     

Effect of nonlinearity on the development of oscillatory magnetohydrodynamic Kelvin-Helmholtz instability in the weakly supercritical regime

The nonlinear stage of development of perturbations at a tangential magnetohydrodynamic discontinuity is investigated in the weakly subcritical and supercritical regimes. It is assumed that the fluid is incompressible and that the density and magnetic field, as well as the velocity, suffer a discontinuity. An equation describing the evolution of low-amplitude nonlinear perturbations is obtained. For periodic perturbations this equation reduces to an infinite system of ordinary differential equations for the amplitudes of the Fourier harmonics. The system is reduced to finite form by truncation and then integrated numerically. Calculations show that the evolution of an initially sinusoidal perturbation always ends with the appearance in the wave profile of an infinite derivative. This can take the form of either an infinitely sharp peak (knife-edge) or wave breaking.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 30–39, May–June, 1988.     

Horizontal mixing layer in shallow water flows

Horizontal-shear thin-layer homogeneous fluid flow in the open channel is considered. A one-dimensional mathematical model of the development and evolution of the horizontal mixing layer is derived within the framework of the three-layer scheme. The steady-state solutions of the equations of motion are constructed and investigated. In particular, supercritical (subcritical)-in-average flow concepts are introduced and the problem of the mixing layer structure is solved. The proposed model is verified on the basis of comparison with a numerical solution of two-dimensional equations of shallow water theory.     

Non-linear convection due to compositional and thermal buoyancy in Earth's outer core

The linear stability of convection due to compositional and thermal buoyancy in Earth's outer core has been investigated. We have obtained the values of Takens-Bogdanov bifurcation points by plotting graphs of neutral curves corresponding to stationary and oscillatory convection for different values of physical parameters. We have derived a non-linear two-dimensional Landau-Ginzburg equation with real coefficients near the onset of stationary convection at a supercritical pitchfork bifurcation and two non-linear one-dimensional coupled Landau-Ginzburg type equations with complex coefficients near the onset of oscillatory convection at a supercritical Hopf bifurcation. We have studied Nusselt number contribution from a Landau-Ginzburg equation at the onset of stationary convection. We have discussed the stability regions of standing and travelling waves. We have also discussed the occurrence of secondary instabilities such as Eckhaus, zigzag and Benjamin-Feir instabilities. We have also derived the non-linear amplitude equation near the Takens-Bogdanov bifurcation point.     

An analytical study of time-delayed control of friction-induced vibrations in a system with a dynamic friction model

We investigate the control of friction-induced vibrations in a system with a dynamic friction model which accounts for hysteresis in the friction characteristics. Linear time-delayed position feedback applied in a direction normal to the contacting surfaces has been employed for the purpose. Analysis shows that the uncontrolled system loses stability via. a subcritical Hopf bifurcation making it prone to large amplitude vibrations near the stability boundary. Our results show that the controller achieves the dual objective of quenching the vibrations as well as changing the nature of the bifurcation from subcritical to supercritical. Consequently, the controlled system is globally stable in the linearly stable region and yields small amplitude vibrations if the stability boundary is crossed due to changes in operating conditions or system parameters. Criticality curve separating regions on the stability surface corresponding to subcritical and supercritical bifurcations is obtained analytically using the method of multiple scales (MMS). We have also identified a set of control parameters for which the system is stable for lower and higher relative velocities but vibrates for the intermediate ones. However, the bifurcation is always supercritical for these parameters resulting in low amplitude vibrations only.     

Numerical Investigation of the Flow Past a Rotating Golf Ball and Its Comparison with a Rotating Smooth Sphere

Large-eddy simulations are conducted for a rotating golf ball and a rotating smooth sphere at a constant rotational speed at the subcritical, critical and supercritical Reynolds numbers. A negative lift force is generated in the critical regime for both models, whereas positive lift forces are generated in the subcritical and supercritical regimes. Detailed analysis on the flow separations on different sides of the models reveals the mechanism of the negative Magnus effect. Further investigation of the unsteady aerodynamics reveals the effect of rotating motion on the development of lateral forces and wake flow structures. It is found that the rotating motion helps to stabilize the resultant lateral forces for both models especially in the supercritical regime.     

The Study of a Parametrically Excited Nonlinear Mechanical System with the Continuation Method

The dynamic behavior of a one-degree-of-freedom, parametrically excited nonlinear system is investigated. The Galerkin method is applied to the principal and fundamental parameteric resonance of the system. The continuation method is used to study the change of harmonic oscillation with respect to the variation of excitation frequency. The numerical stability analysis of the trivial solution is carried out and the stable and unstable regions of the trivial solution are given. They are found to agree with the results obtained by the analytical method of Galerkin. Periodic solutions are traced and the coexistence of multi-periodic solutions is observed With the change of excitation frequency the large amplitude periodic-2 oscillation is found to be in the same closed branch with the small amplitude periodic-2 solution. In addition, the bifurcation pattern of the trivial solution is found to change from subcritical Hopf bifurcation into supercritical Hopf bifurcation with the increase of excitation amplitude. Combined with the conventional numerical integration method, new complex dynamic behavior is detected.     

Effect of Low Rotation Rate on Steady Convection During the Solidification of a Ternary Alloy

We consider the problem of steady convective flow during the directional solidification of a horizontal ternary alloy system rotating at a constant and low rate about a vertical axis. Under the limit of large far-field temperature, the flow region is modeled to be composed of two horizontal mushy layers, which are referred to here as a primary layer over a secondary layer. We first determine the basic state solution and then carry out linear stability analysis to calculate the neutral stability boundary and the critical conditions at the onset of motion. We find, in particular, that there are two flow solutions and each solution exhibits two neutral stability boundaries, and the flow can be multi-modal in the low rotating rate case with local minima on each neutral boundary. The critical Rayleigh number and the wave number as well as the vertical volume flux increase with the rotation rate, but the flow is found to be less stabilizing as compared to the binary alloy counterpart flow. The effects of low rotation rate increase the solid fraction and the liquid fraction at certain vertically oriented fluid lines, and the highest value of such increase is at a horizontal level close to the interface between the two mushy layers.     

Finite-amplitude stability analysis of liquid films down a vertical wall with and without interfacial phase change

The generalized kinematic equation for film thickness, taking into account the effect of phase change at the interface, is used to investigate the nonlinear stability of film flow down a vertical wall. The analysis shows that supercritical stability and subcritical instability are both possible for the film flow system. Applications of the result to isothermal, condensate and evaporate film flow show that mass transfer into (away from) the liquid phase will stabilize (destabilize) the film flow. Finally, we find that supercritical filtered waves are always linearly stable with regard to side-band disturbance.     

Numerical Investigation on Marginal Stability and Convection with and without Magnetic Field in a Mushy Layer

In this article, an investigation is conducted to analyze the marginal stability with and without magnetic field in a mushy layer. During alloy solidification, such mushy layer, which is adjacent to the solidification front and composed of solid dendrites and liquid, is known to produce vertical chimneys. Here, we carry out numerical investigation for particular range of parameter values, which cover those of available experimental studies, to determine the convective flow at the onset of motion. The governing coupled non-linear partial differential equations are non-dimensionalised and solved to get the steady basic-state solution. The thickness of the mushy layer is determined as a part of the solution. Using multiple shooting technique, we determine the steady-state solutions in a range of critical Rayleigh number. We analyse the effect of several parameters, Chandrasekhar number Q, and Robert’s number τ on the problem. It was found that an increase in Q has a stabilizing effect on solidification and the critical Rayleigh number increases on increasing Q. It was also found that for moderate or small values of Robert’s number τ the critical Rayleigh number is mostly insensitive.     

Nonlinear flutter behavior of a plate with motion constraints in subsonic flow

This paper is aimed at presenting the nonlinear flutter peculiarities of a cantilevered plate with motion-limiting constraints in subsonic flow. A non-smooth free-play structural nonlinearity is considered to model the motion constraints. The governing nonlinear partial differential equation is discretized in space and time domains by using the Galerkin method. The equilibrium points and their stabilities are presented based on qualitative analysis and numerical studies. The system loses its stability by flutter and undergoes the limit cycle oscillations (LCOs) due to the nonlinearity. A heuristic analysis scheme based on the equivalent linearization method is applied to theoretical analysis of the LCOs. The Hopf and two-multiple semi-stable limit cycle bifurcation bifurcations are supercritical or subcritical, which is dependent on the location of the motion constraints. For some special cases the bifurcations are, interestingly, both supercritical and subcritical. The influence of varying parameters on the dynamics is discussed in detail. The results predicted by the analysis scheme are in good agreement with the numerical ones.