Longwave Marangoni convection in a surfactant solution between poorly conducting boundaries

2013 ◽  
Vol 718 ◽  
pp. 428-456 ◽  
Author(s):  
Sergey Shklyaev ◽  
Alexander A. Nepomnyashchy

AbstractWe consider Marangoni convection in a heated layer of a binary liquid. The solute is a surfactant, which is present in both surface and bulk phases; the bulk gradient of the concentration is formed due to the Soret effect. Linear stability analysis demonstrates a well-pronounced stabilization of the layer due to the adsorption kinetics and advection of the surface phase. We derive nonlinear amplitude equations for longwave perturbations in the case of fast sorption kinetics (small Langmuir number) and demonstrate that with increase in the effect of the adsorption, subcritical excitation occurs. In the case of a finite Langmuir number, the weakly nonlinear problem is ill-posed. A physical mechanism of subcritical bifurcation is discussed.

Fluids ◽  
2021 ◽  
Vol 6 (8) ◽  
pp. 282
Author(s):  
Alexander B. Mikishev ◽  
Alexander A. Nepomnyashchy

Nonlinear dynamics of patterns near the threshold of long-wave monotonic Marangoni instability of conductive state in a heated thin layer of liquid covered by insoluble surfactant is considered. Pattern selection between roll and square planforms is analyzed. The dependence of pattern stability on the heat transfer from the free surface of the liquid characterized by Biot number and the gravity described by Galileo number at different surfactant concentrations is studied. Using weakly nonlinear analysis, we derive a set of amplitude equations governing the large-scale roll distortions in the presence of the surface deformation and the surfactant redistribution. These equations are used for the linear analysis of modulational instability of stationary rolls.


1996 ◽  
Vol 06 (09) ◽  
pp. 1665-1671 ◽  
Author(s):  
J. BRAGARD ◽  
J. PONTES ◽  
M.G. VELARDE

We consider a thin fluid layer of infinite horizontal extent, confined below by a rigid plane and open above to the ambient air, with surface tension linearly depending on the temperature. The fluid is heated from below. First we obtain the weakly nonlinear amplitude equations in specific spatial directions. The procedure yields a set of generalized Ginzburg–Landau equations. Then we proceed to the numerical exploration of the solutions of these equations in finite geometry, hence to the selection of cells as a result of competition between the possible different modes of convection.


1994 ◽  
Vol 04 (05) ◽  
pp. 1085-1094 ◽  
Author(s):  
MICHAEL BESTEHORN

Pattern formation in fluids with a free flat upper surface is examined. On that surface, the Marangoni effect provides an additional instability mechanism. Based on amplitude equations it is shown that phase instabilities confine the region of stable hexagons to a narrow band of wavelengths. On the other hand we developed a numerical scheme that allows for a direct integration of the fully three-dimensional hydrodynamic equations. There we show the evolution of random patterns and the creation and stabilization of defects as well as the instability of hexagonal patterns lying outside the stable band of wave vectors.


The amplitude equations that govern the motion of wavetrains near the critical point of unstable dispersive, weakly nonlinear physical systems are considered on slow time and space scales T m ═ ε m t ; X m ═ ε m x ( m ═ 1, 2,...). Such systems arise when the dispersion relation for the harmonic wavetrain is purely real and complex conjugate roots appear when a control parameter ( μ ) is varied. At the critical point, when the critical wavevector k c is non-zero, a general result for this general class of unstable systems is that the typical amplitude equations are either of the form ( ∂/∂ T 1 + c 1 ∂/∂ X 1 ) (∂/∂ T 1 + c 2 ∂/∂ X 1 ) A ═ ±α A ─ β AB , ( ∂/∂ T 1 + c 2 ∂/∂ X 1 ) B ═ (∂/∂ T 1 + c 1 ∂/∂ X 1 ) | A | 2 , or of the form ( ∂/∂ T 1 + c 1 ∂/∂ X 1 ) (∂/∂ T 1 + c 2 ∂/∂ X 1 ) A ═ ±α A - β A | A | 2 . The equations with the AB -nonlinearity govern for example the two-layer model for baroclinic instability and self-induced transparency (s. i. t.) in ultra-short optical pulse propagation in laser physics. The second equation occurs for the two-layer Kelvin-Helmholtz instability and a problem in the buckling of elastic shells. This second type of equation has been considered in detail by Weissman. The AB -equations are particularly important in that they are integrable by the inverse scattering transform and have a variety of multi-soliton solutions. They are also reducible to the sine-Gordon equation ϕ ξƬ ═ ± sin ϕ when A is real. We prove some general results for this type of instability and discuss briefly their applications to various other examples such as the two-stream instability. Examples in which dissipation is the dominant mechanism of the instability are also briefly considered. In contrast to the dispersive type which operates on the T 1 -time scale, this type operates on the T 2 -scale.


2014 ◽  
Vol 759 ◽  
pp. 295-320 ◽  
Author(s):  
J. Mizushima ◽  
G. Hatsuda

AbstractNonlinear interactions between the two wakes behind a pair of square cylinders, which are placed side by side in a uniform flow, are investigated by the linear and weakly nonlinear stability analyses and numerical simulations. It is known from the linear stability analysis that the flow past a pair of cylinders becomes unstable to a symmetric or an antisymmetric mode of disturbance, depending on the gap ratio, the ratio of the gap distance between the two cylinders to the cylinder diameter. The antisymmetric mode gives the critical condition for smaller gap ratios than a threshold value, and for larger gap ratios the symmetric mode becomes the most unstable. We focus on the flow pattern arising through the nonlinear interactions of the two modes of disturbance for gap ratios around the threshold value when both modes are growing. We derive a couple of amplitude equations for the two modes to properly describe the nonlinear interaction between them by applying the weakly nonlinear stability theory. The amplitude equations are shown to have three equilibrium solutions except the null solution such as a mixed-mode solution, symmetric and antisymmetric single-mode solutions. Examination of the stability of each equilibrium solution leads to a conclusion that the mixed-mode solution exchange its stability with both the symmetric and the antisymmetric single-mode solutions simultaneously. In the case where the mixed-mode solution is stable, both the symmetric and antisymmetric modes have finite amplitudes, and the resultant flow has an asymmetric flow pattern comprising of finite amplitudes of the two modes of disturbance superposed on the steady symmetric flow. While in the case where both the single-mode solutions are stable, either of the symmetric- and antisymmetric-mode solutions survives, overwhelming the other. Then, if the symmetric mode attains at an equilibrium finite amplitude and the antisymmetric mode vanishes, the resultant flow is symmetric, and if the antisymmetric mode survives and the symmetric mode decays out, the flow becomes asymmetric with the antisymmetric mode of disturbance superposed on the steady symmetric flow. Thus, the flow appearing due to instability differs depending on the initial condition, not uniquely determined, when both single-mode solutions are stable. We numerically delineated the region in the parameter space of the gap ratio and the Reynolds number where the mixed-mode solution is stable. The theoretical results obtained from the weakly nonlinear stability analyses are confirmed by numerical simulations. The conclusion derived from the stability analysis of the equilibrium solutions of the amplitude equations is widely applicable also to other double Hopf bifurcation problems.


Author(s):  
K Fujimura ◽  
S Yamada

On a weakly nonlinear basis, we revisit the pattern formation problem in the Boussinesq convection, for which nonlinear terms of the quadratic order are known to vanish from amplitude equations. It is thus necessary to proceed to the quintic-order approximation in order for the amplitude equations to be generic. By deriving the quintic amplitude equations from the governing PDEs, we examined the bifurcation of steady solutions under rigid–free, rigid–rigid and free–free boundary conditions. Right above the criticality, all the axial solutions are obtained including up- and down-hexagons under the asymmetric boundary conditions and hexagons and regular triangles under the symmetric conditions. Hexagons and regular triangles are unstable whereas rolls are stable as has already been predicted by the cubic-order amplitude equations. Irrespective of the boundary conditions, quintic-order terms stabilize hexagons except near the criticality; rolls and hexagons thus coexist stably in an open region. This suggests that amplitude equations of higher order are possible to predict re-entrant hexagons .


2017 ◽  
Vol 27 (08) ◽  
pp. 1730026 ◽  
Author(s):  
Biao Liu ◽  
Ranchao Wu ◽  
Naveed Iqbal ◽  
Liping Chen

Turing instability and pattern formation in the Lengyel–Epstein (L–E) model with superdiffusion are investigated in this paper. The effects of superdiffusion on the stability of the homogeneous steady state are studied in detail. In the presence of superdiffusion, instability will occur in the stable homogeneous steady state and more complex dynamics will exist. As a result of Turing instability, some patterns are formed. Through stability analysis of the system at the equilibrium point, conditions ensuring Turing and Hopf bifurcations are obtained. To further explore pattern selection, the weakly nonlinear analysis and multiple scale analysis are employed to derive amplitude equations of the stationary patterns. Then complex dynamics of amplitude equations, such as the existence of homogeneous solutions, stripe and hexagon patterns, mixed structure patterns, their stability, interaction and transition between them, are analyzed. Then different patterns occur immediately. Finally, the numerical simulations are presented to show the effectiveness of theoretical analysis and patterns are identified numerically. Whereas in the existing results of such model with normal diffusion, no amplitude equations are derived and patterns are only identified through numerical simulations.


Sign in / Sign up

Export Citation Format

Share Document