Three dimensional simulations of penetrative convection in a porous medium with internal heat sources

The problem of penetrative convection in a fluid saturated porous medium heated internally is analysed. The linear instability theory and nonlinear energy theory are derived and then tested using three dimensions simulation.Critical Rayleigh numbers are obtained numerically for the case of a uniform heat source in a layer with two fixed surfaces. The validity of both the linear instability and global nonlinear energy stability thresholds are tested using a three dimensional simulation. Our results show that the linear threshold accurately predicts the onset of instability in the basic steady state. However, the required time to arrive at the basic steady state increases significantly as the Rayleigh number tends to the linear threshold.     

Horizontal thermal convection in a porous medium

Linearised instability and nonlinear stability bounds for thermal convection in a fluid-filled porous finite box are derived. A nonuniform temperature field in the steady state is generated by maintaining the vertical walls at different temperatures. The linearised instability threshold is shown to be well above the global stability boundary, which is strongly suggested by the lack of symmetry in the perturbed system. The numerical results are evaluated utilising a newly developed Legendre-polynomial-based spectral method.     

Weakly nonlinear EHD stability of slightly viscous jet

A weakly nonlinear approach is utilized here to study the electrohydrodynamic (EHD) instability of an incompressible viscous liquid jet stressed by an axial electric field. The linear motion equations is solved in the light of nonlinear boundary conditions. The viscosity is assumed to be small. The study takes into account both the shear and radial components of the stresses at the interface. In the linear theory, we discuss the breakup phenomena of liquid jets. Also, it is found that, the electrical shearing stresses have no effect at the linear marginal state, while the linear cutoff wavenumber depends on the electrical shearing stresses. A nonlinear perturbation method is introduced. This method can be described our problem precisely. The nonlinear stability is compared with the linear stability condition in the weak viscosity case. It is found that, the weak viscosity has effect on the nonlinear stability condition, in contrast with the linear analysis, whereas the nonlinear cutoff wavenumber doesn't depend on the weak viscosity in both the linear and nonlinear theory.     

Stability of triple diffusive convection in a viscoelastic fluid-saturated porous layer

The triple diffusive convection in an Oldroyd-B fluid-saturated porous layer is investigated by performing linear and weakly nonlinear stability analyses. The condition for the onset of stationary and oscillatory is derived analytically. Contrary to the observed phenomenon in Newtonian fluids, the presence of viscoelasticity of the fluid is to degenerate the quasiperiodic bifurcation from the steady quiescent state. Under certain conditions, it is found that disconnected closed convex oscillatory neutral curves occur, indicating the requirement of three critical values of the thermal Darcy-Rayleigh number to identify the linear instability criteria instead of the usual single value, which is a novel result enunciated from the present study for an Oldroyd-B fluid saturating a porous medium. The similarities and differences of linear instability characteristics of Oldroyd-B, Maxwell, and Newtonian fluids are also highlighted. The stability of oscillatory finite amplitude convection is discussed by deriving a cubic Landau equation, and the convective heat and mass transfer are analyzed for different values of physical parameters.     

An example of nonlinear toxic mass spreading

A one-dimensional, nonlinear problem of reproductive toxic mass spreading is studied in this paper. The nonlinearity is due to the difference of the reproduction rates in the toxic region and the nontoxic region. Multiple steady state solutions are found and their stability and instability are proved. Due to the instability, there may exist turning points (also called saddle-node bifurcation points), at which an infinitesimal perturbation of some parameters may cause a catastrophic change in the location of the steady state toxic front (the interface of the toxic region and the nontoxic region). For the time dependent case, the propagation of the toxic front is considered. An integral equation is derived to determine the propagation of the toxic front. Some numerical results are found for a specific example.     

Instability in gravity-driven flow over uneven permeable surfaces

We develop a theoretical model for inclined free-surface flow over a porous surface exhibiting periodic undulations. The effect of bottom permeability is incorporated by imposing a slip condition that accounts for the nonplanar geometry of the fluid–porous medium interface. Under the assumption of shallow flow, equations of motion accounting for inertial effects are obtained by retaining in the Navier-Stokes equations terms that are up to second-order with respect to a small shallowness parameter. The explicit dependence on the cross-stream coordinate is eliminated from these equations by means of a weighted residual procedure. A linear stability analysis of the steady flow is performed in connection with Floquet–Bloch theory. The results predict that bottom permeability has a destabilizing influence on the flow. A physical explanation has been proposed which involves examining how permeability affects the steady-state flow. Conclusions are drawn regarding the combined effect of the surface tension of the fluid and the parameters describing the bottom surface including permeability, inclination and the amplitude and wavelength of the undulations that generate the bottom topography. A numerical scheme for solving the fully nonlinear governing equations is also outlined. The instability of particular steady flows is determined by conducting nonlinear simulations of the temporal evolution of the flow and comparisons are made with the predictions from the linear analysis. Comparisons with existing experimental data are also included.     

A Nonlinear Stability Analysis for Thermoconvective Magnetized Ferrofluid Saturating a Porous Medium

The generalized energy method is developed to study the nonlinear stability analysis for a magnetized ferrofluid layer heated from below saturating a porous medium, in the stress-free boundary case. The mathematical emphasis is on how to control the nonlinear terms caused by magnetic body force. By introducing a suitable generalized energy functional, we perform a nonlinear energy stability (conditional) analysis. It is found that the nonlinear critical stability magnetic thermal Rayleigh number does not coincide with that of linear instability analysis, and thus indicates that the subcritical instabilities are possible. However, it is noted that, in case of non-ferrofluid, global nonlinear stability Rayleigh number is exactly the same as that for linear instability. For lower values of magnetic parameters, this coincidence is immediately lost. The effect of magnetic parameter, M ₃, and medium permeability, Da, on subcritical instability region has also been analyzed. It is shown that with the increase of magnetic parameter (M ₃) and Darcy number (Da), the subcritical instability region between the two theories decreases quickly. We also demonstrate coupling between the buoyancy and magnetic forces in nonlinear energy stability analysis as well as in linear instability analysis.     

Effects of inertia and surface tension on a power-law fluid flowing down a wavy incline

We consider a thin film of a power-law liquid flowing down an inclined wall with sinusoidal topography. Based on the von Kármán–Pohlhausen method an integral boundary-layer model for the film thickness and the flow rate is derived. This allows us to study the influence of the non-Newtonian properties on the steady free surface deformation. For weakly undulated walls we solve the governing equation analytically by a perturbation approach and find a resonant interaction of the free surface with the wavy bottom. Furthermore, the analytical approximation is validated by numerical simulations. Increasing the steepness of the wall reveals that nonlinear effects like the resonance of higher harmonics grow in importance. We find that shear-thickening flows lead to a decrease while shear thinning flows lead to an amplification of the steady free surface. A linear stability analysis of the steady state shows that the bottom undulation has in most cases a stabilizing influence on the free surface. Shear thickening fluids enhance this effect. The open questions which occurred in the linear analysis are then clarified by a nonlinear stability analysis. Finally, we show the important role of capillarity and discuss its influence on the steady solution and on the stability.     

Numerical study of convection of a liquid heated from below

The equilibrium of a liquid heated from below is stable only for small values of the vertical temperature gradient. With increase of the temperature gradient a critical equilibrium situation occurs, as a result of which convection develops. If the liquid fills a closed cavity, then there is a discrete sequence of critical temperature gradients (Rayleigh numbers) for which the equilibrium loses stability with respect to small characteristic disturbances. This sequence of critical gradients and motions may be found from the solution of the linear problem of equilibrium stability relative to small disturbances. If the temperature gradient exceeds the lower critical value, then (for steady-state heating conditions) there is established in the liquid a steady convective motion of a definite amplitude which depends on the magnitude of the temperature gradient. Naturally, the amplitude of the steady convective motion cannot be determined from linear stability theory; to find this amplitude we must solve the problem of convection with heating from below in the nonlinear formulation. A nonlinear study of the steady motion of a liquid in a closed cavity with heating from below was made in [1]. In that study it was shown that for Rayleigh numbers R which are less than the lower critical value R_c steady-state motions of the liquid are not possible. With R>R_c a steady convection arises, whose amplitude near the threshold is small and proportional to (R–R_c)^1/2 (the so-called soft instability)-this is in complete agreement with the results of the phenom-enological theory of Landau [2, 3].Primarily, various versions of the method of expansion in powers of the amplitude [4–8] have been used, and, consequently, the results obtained in those studies are valid only for values of R which are close to R_c, i. e., near the convection threshold.It is apparent that the study of developed convective motion far from the threshold can be carried out only numerically, with the use of digital computers. In [9, 10] the numerical methods have been successfully used for the study of developed convection in an infinite plane horizontal liquid layer.The present paper undertakes the numerical study of plane convective motions of a liquid in a closed cavity of square section. The complete nonlinear system of convection equations is solved by the method of finite differences on a digital computer for various values of the Rayleigh number, the maximal value exceeding by a factor of 40 the minimal critical value R_c. The numerical solution permits following the development of the steady motion which arises with R>R_c in the course of increase of the Rayleigh number and permits study of the oscillatory motions which occur at some value of the parameter R. The heat transfer through the cavity is studied. The corresponding linear problem on equilibrium stability is solved approximately by the Galerkin method.     

Instability of solution for a class of the third order nonlinear differential equation

ＩＮＳＴＡＢＩＬＩＴＹＯＦＳＯＬＵＴＩＯＮＦＯＲＡＣＬＡＳＳＯＦＴＨＥＴＨＩＲＤＯＲＤＥＲＮＯＮＬＩＮＥＡＲＤＩＦＦＥＲＥＮＴＩＡＬＥＱＵＡＴＩＯＮＬｕＤｅｙｕａｎ卢德渊（ＲｅｃｅｉｖｅｄＮｏｖ１０１９９４ＣｏｍｍｕｎｉｃａｔｅｄｂｙＺｈａｎｇＳｈｉ...     

Global properties of a class of virus infection models with multitarget cells

In this paper, we propose a class of virus infection models with multitarget cells and study their global properties. We first study three models with specific forms of incidence rate function, then study a model with a more general nonlinear incidence rate. The basic model is a (2n+1)-dimensional nonlinear ODEs that describes the population dynamics of the virus, n classes of uninfected target cells, and n classes of infected target cells. Model with exposed state and model with saturated infection rate are also studied. For these models, Lyapunov functions are constructed to establish the global asymptotic stability of the uninfected and infected steady states of these models. We have proven that if the basic reproduction number is less than unity then the uninfected steady state is globally asymptotically stable, and if the basic reproduction number is greater than unity then the infected steady state is globally asymptotically stable. For the model with general nonlinear incidence rate, we construct suitable Lyapunov functions and establish the sufficient conditions for the global stability of the uninfected and infected steady states of this model.     

进动圆筒内液体流动不稳定的非线性特征

在建立进动充液圆筒内液体偏差流动方程的基础上，结合液体惯性波和轴向二次流动线性解，通过对定常二次流动的线性稳定性分析，提出了函数空间表达的流动不稳定性非线性分岔分析方程. 对非惯性坐标系下液体流动的Navier-Stokes方程进行了数值求解，并对惯性波发生破裂(实验提供的3种主模态下得出的共振破裂现象)时的压力时间序列进行分析，得出了液体流动不稳定的基本非线性特征.     

Weakly nonlinear stability analysis of a falling film with countercurrent gas flow

Weakly nonlinear stability analysis of a falling film with countercurrent gas–liquid flow has been investigated. A normal mode approach and the method of multiple scales are employed to carry out the linear and nonlinear stability solutions for the film flow system. The results show that both supercritical stability and subcritical instability are possible for a film flow system when the gas flows in the countercurrent direction. The stability characteristics of the film flow system are strongly influenced by the effects of interfacial shear stress when the gas flows in the countercurrent direction. The effect of countercurrent gas flow in a falling film is to stabilize the film flow system.     

Convective Stability of a Horizontal Binary-Mixture Layer in a Modulated External Force Field

The thermovibrational instability of a plane horizontal layer of incompressible binary mixture is analyzed with account for the Soret thermodiffusion effect. To find the stability threshold in the linear approximation, Floquet's theory is used. It is shown that, depending on the amplitude and frequency, vibration can stabilize an unstable ground state or destabilize fluid equilibrium. Apart from the synchronous or subharmonic responses to external influences, instability may be related to quasi-periodic perturbations. The behavior of the threshold values at the low-frequency limit is considered.     

Görtler Vortices with System Rotation

The steady primary instability of Görtler vortices developing along a curved Blasius boundary layer subject to spanwise system rotation is analysed through linear and nonlinear approaches, to clarify issues of vortex growth and wavelength selection, and to pave the way to further secondary instability studies.A linear marching stability analysis is carried out for a range of rotation numbers, to yield the (predictable) result that positive rotation, that is rotation in the sense of the basic flow, enhances the vortex development, while negative rotation dampens the vortices. Comparisons are also made with local, nonparallel linear stability results (Zebib and Bottaro, 1993) to demonstrate how the local theory overestimates vortex growth. The linear marching code is then used as a tool to predict wavelength selection of vortices, based on a criterion of maximum linear amplification.Nonlinear finite volume numerical simulations are performed for a series of spanwise wave numbers and rotation numbers. It is shown that energy growths of linear marching solutions coincide with those of nonlinear spatially developing flows up to fairly large disturbance amplitudes. The perturbation energy saturates at some downstream position at a level which seems to be independent of rotation, but that increases with the spanwise wavelength. Nonlinear simulations performed in a long (along the span) cross section, under conditions of random inflow disturbances, demonstrate that: (i) vortices are randomly spaced and at different stages of growth in each cross section; (ii) upright vortices are the exception in a universe of irregular structures; (iii) the average nonlinear wavelengths for different inlet random noises are close to those of maximum growth from the linear theory; (iv) perturbation energies decrease initially in a linear filtering phase (which does not depend on rotation, but is a function of the inlet noise distribution) until coherent patches of vorticity near the wall emerge and can be amplified by the instability mechanism; (v) the linear filter represents the receptivity of the flow: any random noise, no matter how strong, organizes itself linearly before subsequent growth can take place; (vi) the Görtler number, by itself, is not sufficient to define the state of development of a vortical flow, but should be coupled to a receptivity parameter; (vii) randomly excited Görtler vortices resemble and scale like coherent structures of turbulent boundary layers.A.Z. has been supported, during his stay at EPFL, by an ERCOFTAC Visitor Grant. A.B. acknowledges the Swiss National Fund, Grant No. 21-36035.92, for travel support associated with this research. This work was also supported by the Swedish Board of Technical Development (NUTEK), the Swedish Technical-Scientific Council (TFR), and an ERCOFTAC Visitor Grant, through which the stay of B.G.B.K. at the EPFL was made possible. Cray-2 computing time for this research was generously provided by the Service Informatique Centrale of EPFL.     

Nonlinear Stability of Convection in a Porous Layer with Solid Partitions

We show that for many classes of convection problem involving a porous layer, or layers, interleaved with finite but non-deformable solid layers, the global nonlinear stability threshold is exactly the same as the linear instability one. The layer(s) of porous material may be of Darcy type, Brinkman type, possess an anisotropic permeability, or even be such that they are of local thermal non-equilibrium type where the fluid and solid matrix constituting the porous material may have different temperatures. The key to the global stability result lies in proving the linear operator attached to the convection problem is a symmetric operator while the nonlinear terms must satisfy appropriate conditions.     

The uniqueness of the solution of boundary-value problems of a thermal model of discharge

The paper is devoted to a nonlinear analysis of superheating [1, 2] instability of an electric discharge stabilized by electrodes [3] in the framework of a thermal model [4] where the stability of the discharge relative to the long-wave and short-wave perturbations is proved in a linear approximation. Similar boundary-value problems arise in the theories of chemically and biologically reacting mixtures [5–7], thermal breakdown of dielectrics [8], thermal explosion [9], in the investigation of nonlinear waves in semiconductors and superconductors [10, 11], and in the investigation of Couette flow with variable viscosity [12]. The uniqueness of the one-dimensional steady solutions of the thermal model of discharge and the stability relative to the small spatial perturbations, respectively, for the exponential and step dependence of the electrical conductivity on the temperature are proved in [3, 13]. The uniqueness of the solutions in the one-dimensional case for the same electrode temperature and arbitrary dependences of the electrical and thermal conductivity on the temperature is established in paper [14]. In the present paper, the existence and uniqueness of steady solutions of the thermal model of discharge in a three-dimensional formulation for arbitrary fairly smooth electrical and thermal conductivity functions of the temperature in the case of isothermal isopotential electrodes are proved analytically.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 140–145, January–February, 1986.The author expresses his gratitude to A. G. Kulikovskii and A. A. Barmin for the formulation of the problem and their discussions.     

Marangoni instability in a liquid layer bounded by two coaxial cylinder surfaces

A liquid layer, confined between two coaxial cylinder surfaces, has either a gasliquid interface on the inside and is heated from the outer (solid) boundary, or it has a gas-liquid interface on the outside and is heated from the inner (solid) boundary. Neglecting gravity and using a standard normal-mode approach, we analyse surface-tension driven instability (Marangoni instability) of the motionless steady state in which the temperature depends on the radial coordinate only. Numerical results for the critical Marangoni number and corresponding wave-number pair are presented for various values of the curvature of the interface. This curvature turns out to exert a significant influence on the onset of Marangoni convection flows. Further, the stability behaviour of the system is found to be quite variable, depending on whether the interface is on the inside or on the outside of the layer and whether it is well-conducting or nearly-isolated.     

磁悬浮列车系统的随机最优控制

本文根据磁悬浮列车和车行道的结构特点，将总系统模型按分块原则分成列车、磁执行环节和车行道系统，并在平衡点附近对非线性方程线性化处理，形成末加控制的总系统的状态方程，它是一组考虑外干扰情况下线性时变系统模型。而基于电磁关系原则建立的磁悬浮列车系统模型在末加控制状态下是不稳定的，为了保证列车的行驶舒适性、稳定性及可靠性，承重磁铁与导向磁铁必须加以控制。附加控制方程后，就形成了被控制的总系统的状态方程，从而实现车、磁及车行道模型的有机组合。对于实际工程问题，被控制的总系统的动力学性质由于维数较高，直接计算比较困难，本文采用计算机进行数值仿真，利用随机最优控制理论，对系统悬浮气隙和垂向加速度的变化规律进行了研究，并通过实例给出时变系统的仿真结果。     

The nonlinear dynamics of elastic tubes conveying a fluid

The Kirchhoff equations for elastic tubes are modified to include the effect of fluid flow. Using the techniques of linear and nonlinear analysis specially developed for the Kirchhoff equations, the effect of the fluid flow on the basic twist-to-writhe instability is investigated. The results suggest an intriguing modification of the bifurcation threshold due to the flow. Beyond threshold the buckled tube acquires a slight curvature which modifies the flow rate and results in a correction to nonlinearity of the amplitude equation governing the deformation dynamics.