Next Article in Journal
A Molecular View of the Surface Pressure/Area Per Lipid Isotherms Assessed by FTIR/ATR Spectroscopy
Next Article in Special Issue
Dynamic Light Scattering: A Powerful Tool for In Situ Nanoparticle Sizing
Previous Article in Journal
Structural Modification of Jackfruit Leaf Protein Concentrate by Enzymatic Hydrolysis and Their Effect on the Emulsifier Properties
Previous Article in Special Issue
The Contact Angle Hysteresis Puzzle
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Marangoni Patterns in a Non-Isothermal Liquid with Deformable Interface Covered by Insoluble Surfactant

by
Alexander B. Mikishev
1,* and
Alexander A. Nepomnyashchy
2
1
Department of Engineering Technology, Sam Houston State University, Huntsville, TX 77341, USA
2
Department of Mathematics, Technion-Israel Institute of Technology, Haifa 32000, Israel
*
Author to whom correspondence should be addressed.
Colloids Interfaces 2022, 6(4), 53; https://doi.org/10.3390/colloids6040053
Submission received: 14 August 2022 / Revised: 18 September 2022 / Accepted: 26 September 2022 / Published: 29 September 2022
(This article belongs to the Special Issue Feature Reviews in Colloids, Nanomaterials, and Interfaces)

Abstract

:
Marangoni patterns are created by instabilities caused by thermocapillary and solutocapillary stresses on the deformable free surface of a thin liquid layer. In the present paper, we consider the influence of the insoluble surfactant on the selection and modulational instability of stationary Marangoni patterns near their onset threshold. The basic governing parameters of the problem are the Biot number characterizing the heat-transfer resistances of and at the surface, the Galileo number indicating the role of gravity via viscous forces, and the elasticity number specifying the influence of insoluble surfactant on the interfacial dynamics of the liquid. The paper includes a review of the previous results obtained in that problem as well as new ones.

1. Introduction

Interfacial phenomena are significant in many natural and technological processes. A basic phenomenon is the interfacial pattern formation discovered by H. Bénard at the end of the 19th century [1]: under homogeneous external conditions, the system becomes spatially inhomogeneous. The correct theoretical explanation of these experimental results was proposed more than 50 years later in the 1950s of the 20th century by Block [2] and, independently, by Pearson [3]. They both explained the interfacial pattern formation as a result of surface tension variations of a liquid free surface due to temperature changes along the surface.
The interest in the pattern formation was awoken by the seminal paper of Turing [4], who predicted theoretically a spontaneous pattern generation in a chemical system with two reacting and diffusing chemicals (found experimentally in the 1990s [5]), and by the paper of Normand, Pomeau, and Velarde [6], who attracted the attention of physicists to the convective pattern formation. The most typical convective patterns are rolls, squares, and hexagons. Numerous examples of patterns can be found in the recent publications [7,8]. Now, we know that the process of pattern formation is typical for physical, chemical, and biological systems which are far from the thermodynamic equilibrium state.
The longwave Marangoni convection is the focus of the special interest of the experts on pattern formation because of a rich gallery of patterns that can be formed. There exist two kinds of longwave Marangoni instabilities: (i) the first is the already mentioned Pearson’s mode which is possible in the case of poorly conducting boundaries, and (ii) the mode described by Scriven and Sternling [9] existing in the case of a thin film. The evolution of longwave instability is described by a slowly evolving "active" variable that determines the nonlinear dynamics. The study of the pattern formation caused by the Pearson’s mode was started in the works of Sivashinsky [10] and Knobloch [11]. In this case, the temperature plays the role of the active variable. In the case of the mode of Scriven and Sternling, the active variable is the free-surface deformation. The nonlinear evolution of this mode was described by Davis [12].
Recently, a remarkable interest in the Marangoni patterns generated by the longwave Marangoni instability was motivated by the work of Shklyaev et al. [13]. The authors showed that in the interval of wavenumbers k O ( B i 1 / 2 ) ( B i 1 is the Biot number, defined as the ratio of heat-transfer resistances inside of and at the surface of the liquid), the temperature disturbance and the surface deformation are both active variables. Their interaction can generate large-scale monotonic and oscillatory instabilities, which produce stationary patterns and wave patterns [14].
In applications such as inject printing [15,16], spray cooling [17], the fabrication of DNA/RNA microarrays [18], photolithography [19], and various pulmonary therapies [20], one has to deal with the influence of surfactant contaminants. The nonuniform deposition of surfactant over the liquid surface can be a cause of special surfactant-induced solutal Marangoni stresses that significantly change the patterning profile of a liquid. Even a small concentration of impurities can cause the solutal stresses that interact with the thermal Marangoni flow [21,22]. Recently, the experimental and theoretical works in colloid science showed that the formation of unusual patterns (from rings and the “Marangoni ridge” to hexagonal arrays) are related to the Marangoni instabilities of deposited liquid films [23,24,25,26]. A highly ordered pattern is observed when a milimeter-size drop of dichloromethane spreads on an aqueous substrate under a change in Marangoni stresses [27]. The surfactant adsorption at a droplet interface is capable of modulating the flow patterns via the wetting shape, varying the liquid–vapor interfacial properties, and inducing superficial flows [28].
The conservation of the amount of the absorbed insoluble surfactant is the origin of an additional slow long-scale process: the redistribution of the adsorbed surfactant due to the surface flow and diffusion. The surfactant concentration is one more active variable. Though the presence of an insoluble surfactant on the surface provides no specific instability mechanism, it can significantly modify the development of instability generated by the thermocapillary instability.
In the present work, we review the latest results on the large-scale Marangoni convection in a liquid layer with an insoluble surfactant spread over a deformable liquid interface. The literature on the Marangoni convection with a soluble surfactant was reviewed in detail by Shklyaev and Nepomnyashchy in [8]. We present a weakly nonlinear analysis of the patterns on square and rhombic lattices in the Fourier space near the instability threshold. The modulation of these convective patterns is considered using the Newell, Whitehead, and Segel [29,30] approach, describing the interaction of the disturbances with various wavenumbers close to the critical one. As a result, this modulation is described by the Ginzburg–Landau-type equation for the amplitude of patterns coupled with equations for surface distortion and surfactant concentration. We discuss the role of 2D modulation on the roll stability. Several cases have been revealed. In each case, a specific stability map in terms of the Biot and Galileo numbers is constructed. In all these cases, the special role of insoluble surfactant is investigated. The results of the 1D modulation of rolls have been formerly published in [31]. All other results are presented here for the first time. At the end, the perspectives of the future research are discussed.

2. Description of Longwave Marangoni Convection with Insoluble Surfactant

2.1. Formulation of the Problem

The presence of an insoluble surfactant changes the interfacial stresses on the deformable surface of the thin liquid layer. In the case of a non-isothermal liquid layer (the liquid is subjected to a transverse fixed temperature gradient a , ( a > 0 ) [32], the surfactant provides an additional factor, a solutocapillary mechanism, to generate Marangoni stresses on the liquid surface. The surface tension can be written as a linear function of both the temperature T and concentration of surfactant Γ ,
σ = σ 0 σ 1 ( T T 0 ) σ 2 ( Γ Γ 0 ) ,
and here, σ 0 is the reference value of the surface tension σ 1 = T σ , σ 2 = Γ σ , and T 0 and Γ 0 are the reference values of the corresponding variables at the surface in the absence of convection. The absorbed surfactant at the free surface is convected and diffuses over the free-liquid interface but not into the bulk. The parameter that describes the relation between surfactant diffusivity D 0 and thermal diffusivity χ of the liquid is called the Lewis number, L = D 0 / χ . Typically, the Lewis number is small, L 0.01 .
The parameter characterizing the influence of the insoluble surfactant on the interfacial dynamics of the liquid is the elasticity number, N = σ 2 d 0 Γ 0 / η χ . It determines the significance of the solutocapillary effect as compared to viscous forces. It is a changeable parameter. Here, d 0 is a mean thickness of the liquid layer and η is a dynamic viscosity ( η = ρ ν , ρ is a density of the liquid, and ν is a kinematic viscosity of the liquid).
Let us discuss the characteristic values of the surfactant parameters, N and L. While the elasticity number, N, which depends on the surfactant concentration, can vary in a relatively large interval, the value of the surfactant Lewis number has a significant limitation. As an example of an insoluble surfactant, let us consider C 12 E O 8 , produced by Nikko Chemical Co. (Tokyo, Japan); for its properties, see [33]. The typical surfactant values of the surface surfactant diffusion coefficient are (3.0–4.0) × 10 10 ( m 2 / s ) [34]. The surfactant is located at the air/water interface, and for other parameters, we take the physical values of the water at the room temperature. It means the value of the thermal diffusivity coefficient, χ , is 1.3 × 10 7 ( m 2 / s ) , and the kinematic viscosity, ν , is 0.9 × 10 6 ( m 2 / s ) . For the calculation, we fix the Lewis number at L = 0.003 .
The schematic of the problem is presented in Figure 1. Here, the liquid layer is bounded at the bottom by a rigid substrate and has an upper free surface. The coordinate axes ( x , y ) are in the plane of the substrate and the z axis points upward.
The typical additional dimensionless parameters of the instability problems are the Marangoni number, M = σ 1 a d 0 2 / η χ , the ratio of the thermocapillary surface force to the viscous force; the Galileo number, G = g d 0 3 / ν χ , the ratio between the gravity force and the viscous force; the inverse capillary number, Σ = σ 0 d 0 / η χ , the ratio between the surface tension force and the viscous force; and the Biot number, B i = q d 0 / Λ T (q is the heat-transfer coefficient and Λ T is the thermal conductivity), the ratio of the heat-transfer resistances inside of and at the surface of the liquid.
The mathematical description of the interfacial dynamics includes the surface transport equation for the surfactant. Traditionally, some authors use this equation in the form proposed by Levich [35]. However, it works only in the case of a nondeformable free surface. For a deformable surface, the equation modified by Wong et al. [36,37] needs to be used.
Let us describe the free-liquid interface as z = h ( x , y , t ) . Then, in the general case, it is necessary to distinguish between the liquid velocity v and the velocity along the surface, v s = v ( v · n ) n . Here, n = ( x h , y h , 1 ) ( 1 + ( x h ) 2 + ( y h ) 2 ) 1 / 2 is the unity vector normal to the interface. Then, the evolution equation for the surfactant concentration Γ can be written as
t Γ h ( e z · s ) Γ + s · ( v s Γ ) + ( s · n ) ( v · n ) Γ = D 0 s 2 Γ ,
where s = n ( n · ) denotes the operator of the gradient along the interface, = ( x , y , z ) , e z = ( 0 , 0 , 1 ) . The first two terms of Equation (1) describe the temporal change of Γ along the normal to the moving interface. The last term in the left-hand side describes the change in the surfactant concentration due to the temporal change in the surface area [38].

2.2. The Case of Perturbations with k O ( B i 1 / 4 )

It turns out that in the case B i 1 , there are two different characteristic spatial scales of disturbances, k = O ( B i 1 / 4 ) and k = O ( B i 1 / 2 ) . The former scaling has been discovered in the pioneering work of Pearson [3] for a nondeformable surface. It was shown that in the case B i = O ( ϵ 4 ) , ϵ 1 , the minimum of the neutral curve is shifted to the longwave region, k = O ( ϵ ) . We rescale the wavenumber as k = ϵ K . The analysis of the Marangoni convection with the insoluble surfactant adsorbed at the deformable free surface has been conducted in [39]. That analysis based on the longwave expansions shows the existence of two modes of instability—monotonic and oscillatory. At the leading order, the critical value of the monotonic Marangoni number is
M m 0 = 12 G ( 4 L + N ) L ( G + 72 ) .
The monotonic mode exists if
N < N * = 4 ( 216 + G ( G + 27 ) ) L 2 L ( 864 + G ( G + 36 ) ) G ( G + 72 ) .
The second-order correction has the form M m 2 = c 1 K 2 + c 2 K 2 , where c 1 and c 2 are constants. For N > N * , the instability becomes oscillatory with the critical Marangoni number
M o s c = 1 2 ( G + 27 + 48 L ) G 2 + 3 G ( 41 + 32 L + 5 N ) + 72 [ 5 N + 3 + 2 L ( L + N + 2 ) ] D ,
where
D = 144 { 4 ( 216 + G ( G + 27 ) ) L 2 G ( G + 72 ) N + ( 864 + G ( G + 36 ) ) L N } + + { G 2 9 G ( N 3 ) + 72 ( N + 3 + 2 L ( L + N ) ) } 2 .
The neutral oscillations have the frequency ω = K 2 Ω 0 with
Ω 0 2 = 1 288 { G 2 + 9 G ( N 3 ) 72 ( 3 + N + 2 L ( L + N ) ) + D } .
The deformability of the free surface is restricted by the conservation law of the volume that gives us the equation for the local layer thickness H. Two others relate to the conservation of energy described by the equation for the bottom temperature F and to the conservation of the mass of the surfactant described by the equation for the concentration of surfactant Γ . If all the parameters Σ , M, and G are of order of unity, at the leading order, one obtains the following system:
τ H = · G 3 H 3 H + M 2 H 2 θ + N 2 H 2 Γ · Q ˜ 1 ,
H τ F = · G H 4 8 H + M H 3 6 θ + N H 3 6 Γ + H F + Q ˜ 1 · θ 1 2 ( H ) 2 ,
τ Γ = · L Γ + Γ H G 2 H H + M θ + N Γ · Q ˜ 2 .
Here, θ = F H is a temperature perturbation on the free surface, = ( X , Y ) . Note that parameters Σ and B i do not appear in the leading-order equations. The obtained longwave system of equations is ill-posed [40], because the stabilizing action of the surface tension is not taken into account. To regularize the obtained problem, it is necessary to assume that the surface tension is strong: Σ = ϵ 2 S , S = O ( 1 ) . The strong surface tension suppresses shortwave disturbances, and the assumption of longwave instability is self-consistent.
The case of the finite value of the Biot number was considered in the problem on evaporation of the liquid layer with an insoluble surfactant [37], as well as in the problem of the formation of Faraday waves at the surfactant-covered free surface of a vertically vibrated liquid layer [41].
The influence of surfactants on convective instabilities in two-layer systems was discussed in [42,43,44].

2.3. The Case of Perturbations with k O ( B i 1 / 2 )

Another distinguished limit of the problem was discovered by Shklyaev et al. [13], where the authors showed that for the deformable interface, the disturbances in the interval k B i 1 / 2 without surfactant had two modes of instability, monotonic and oscillatory, that can be generated. This interval of wavenumbers determines the appropriate scaling for the development of possible nonlinear longwave structures in the liquid. One can rescale the spatial coordinates as
X = ϵ x , Y = ϵ y , 0 < ϵ 1 .
The characteristic time is proportional to k 2 and one can rescale the temporal coordinate as
τ = ϵ 2 t .
In the presence of an insoluble surfactant, in the framework of that scaling, we derived the following system of the longwave amplitude equations for the local thickness H ( X , Y , τ ) , perturbations of temperature F ( X , Y , τ ) , and surfactant concentration Γ ( X , Y , τ ) , see [45]:
τ H = · H 3 3 R + M H 2 2 θ + N H 2 2 Γ · Q 1 ,
H τ F = · H 4 8 R + M H 3 6 θ + N H 3 6 Γ + H F + Q 1 · θ 1 2 ( H ) 2 β θ ,
τ Γ = · Γ H H 2 R + M θ + N Γ + L Γ · Q 2 .
Here, R = G H S 2 H , = ( X , Y ) . The Biot number is B i = ϵ 2 β .
The system (7)–(9) describes the nonlinear dynamics of longwave perturbations and includes the following effects: Equation (7) includes the evolution of the surface deformation due to the bulk flow with the flow rate Q 1 generated by hydrostatic and Laplace pressures and both thermocapillary and solutocapillary effects. Equation (8) describes the advective heat transfer by the flow and heat conductivity in the longitudinal direction, and the last two terms of the equation describe the heat loss from the deformable interface. Finally, Equation (9) presents the evolution of the surfactant concentration due to the diffusion and surface flow caused by the gravity and surface tension, as well as by thermo- and solutocapillarity.
The linear analysis of (7)–(9) was performed in [45] and gave two instability modes: the monotonic one and oscillatory one. The wavelength of disturbances is rescaled as K = ϵ 1 k . The neutral stability curve is described for the monotonic mode as
M m ( K ) = ( 48 + 12 N / L ) ( G + K 2 S ) ( β + K 2 ) K 2 ( 72 + G + K 2 S ) .
Formula (2) is recovered in the limit K 2 , K 2 S 0 . The minimum of the marginal stability curve is reached at the critical wavenumber:
( K c ) 2 = G S β + 72 G S β ( 72 + G S β ) S ( 72 S β ) .
Note that without the loss of generality, one can choose S = 1 , which corresponds to the definition ϵ = Σ 1 / 2 .
In both limits, we arrive at a conclusion that in the considered problems, despite the appearance of a new oscillatory in the mode, increasing the concentration of the insoluble surfactant suppresses the instabilities creating the convection.

3. Pattern Selection in the Longwave Marangoni Convection

3.1. Square Lattice

We consider the nonlinear dynamics of small spatially periodic perturbations in the neighborhood of the threshold of the monotonic mode, M = M m + δ 2 M 2 (here, δ is a small parameter of supercriticality), corresponding to a square lattice in the Fourier space. At the present stage, only disturbances with the critical wavevector K c determined by (11) are taken into account. The plots showing basic wavevectors for all considered lattices are shown in Appendix A.
The problem is described by the set of Equations (7)–(9) for variables H, F, and Γ , which are expanded in powers of parameter δ :
H = 1 + δ h 1 + δ 2 h 2 + , F = 1 + δ f 1 + δ 2 f 2 + , Γ = 1 + δ γ 1 + δ 2 γ 2 + .
Substituting these expansions into (7)–(9), we obtain at the leading order the linear equations for h 1 , f 1 , and γ 1 . We can write their solution in the form
h 1 = A 1 ( τ 2 ) e i K c X + A 2 ( τ 2 ) e i K c Y + c . c . , f 1 = α 1 h 1 , γ 1 = α 2 h 1 ,
where α 1 = ( G K 2 + K 4 S 72 β ) 72 ( K 2 + β ) and α 2 = G + K 2 S 6 L .
Continuing the operations for higher powers of δ and applying the solvability conditions at the third order in δ , we obtain the set of Landau equations that describe the evolution of complex amplitudes A 1 and A 2 :
A 1 τ 2 = κ 0 A 1 + κ 1 | A 1 | 2 A 1 + κ 2 | A 2 | 2 A 1 ,
A 2 τ 2 = κ 0 A 2 + κ 1 | A 2 | 2 A 2 + κ 2 | A 2 | 2 A 2 .
Here,
κ 0 = K 4 L 2 ( 72 + G + K 2 S ) 2 M 2 12 ( K 2 D 1 β D 2 )
where
D 1 = N ( G + K 2 S ) ( 72 + G + K 2 S ) + L N ( 864 + 36 G + G 2 + 2 ( 18 + G ) K 2 S + K 4 S 2 ) + 4 L 2 ( 216 + 27 G + G 2 + ( 27 + 2 G ) K 2 S + K 4 S 2 ) , D 2 = 36 L N ( 24 + G + K 2 S ) + N ( G + K 2 S ) ( 72 + G + K 2 S ) + 36 L 2 ( 24 + 5 G + 5 K 2 S ) .
Other coefficients are cumbersome to be given here (see Supplemental Material to [45]).
We have two kinds of steady-state solutions: one of them describes rolls (if one of A 1 or A 2 is zero) and another one describes squares (both A 1 and A 2 are nonzero). Here, the pattern selection is determined by signs of κ 1 , κ 1 κ 2 , and κ 1 + κ 2 . The rolls and squares can be stable if κ 1 < 0 and κ 1 + κ 2 < 0 (supercritical bifurcation). The rolls are selected if κ 2 < κ 1 < 0 ; if κ 1 < κ 2 , κ 1 + κ 2 < 0 , the squares are stable, see Figure 2.
Note that the problem of selection between the roll and square patterns is mathematically similar to the classical problem of species competition [46]. The selection of rolls corresponds to the extinction of all species except one, while the selection of squares corresponds to the symbiosis.
The pattern selection domains are presented in Figure 2 for three different cases. The red color is for the case without surfactant, the blue color relates to N = 10 6 , and the black color is for N = 10 5 . It is seen that even at a low concentration of surfactant, the domains shrink. Figure 2b shows domains at small values of β .
It was found that the rolls are the most typical of the patterns. Note that this kind of structures is the only possible one in a rectangular container with one size much longer than another one or in a narrow annular container of a large radius.

3.2. Hexagonal Lattice

The most typical kind of non-equilibrium patterns near the threshold of stationary is the hexagonal pattern. Hexagons are formed by three sets of rolls, with wavenumbers satisfying the resonant condition k 1 + k 2 + k 3 = 0 , | k 1 | = | k 2 | = | k 3 | . Well-known examples are irregular hexagons appearing in the heated cooking oil and giraffe’s coat markings (see [14]); also, hexagonal patterns were observed in front solidification [47], in Faraday crispation [48], in a liquid-crystal-valve device [49], in ferrofluids [50], etc.
Here, we find how the surfactant concentration affects the type of hexagons in the longwave Marangoni convection. As in the previous section, we take the variable fields H ( X , Y , τ ) , F ( X , Y , τ ) , and Γ ( X , Y , τ ) near the threshold point M m ( K c ) (see (10) and (11)). The Marangoni number is presented as M = M m + δ M 1 + δ 2 M 2 + . We introduce the different time scales in the neighborhood of the threshold, i.e., τ 0 = τ , τ 1 = δ τ , τ 2 = δ 2 τ 2 , etc. Then, at the leading order of the expansions, at the hexagonal lattice in the Fourier space, we have the solution
h 1 = j = 1 3 A j ( τ 1 , τ 2 ) e i K c n j · X + c . c , f 1 = α 1 h 1 , γ 1 = α 2 h 1 .
Expression (17) corresponds to six peaks in the Fourier space around the points ± K c n j ( j = 1 , 2 , 3 ) , n 1 = ( 1 , 0 ) , n 2 = ( 1 / 2 , 3 / 2 ) , n 3 = ( 1 / 2 , 3 / 2 ) . The vectors n j are unit vectors parallel to pattern wavevectors. We also define mutually orthogonal to n j unit vectors t j , namely t 1 = ( 0 , 1 ) , t 2 = ( 3 / 2 , 1 / 2 ) , t 3 = ( 3 / 2 , 1 / 2 ) . Here, X = ( X , Y ) .
At the second order of the expansion, the solvability condition for h 1 , f 1 , and γ 1 gives
τ 1 A l = κ 0 A l + s 1 A m * A n * , { l , m , n } = { 1 , 2 , 3 } , { 2 , 3 , 1 } , { 3 , 1 , 2 } .
If s 1 = O ( 1 ) , then Equation (18) describes a two-sided bifurcation and unbounded growth of the selected type of hexagons, up-hexagons or down-hexagons, depending on the sign s 1 . However, there exists the line in the space of parameters, where s 1 = 0 ; hence, the term s 1 A m * A n * vanishes. The line s 1 ( G , β ) = 0 is shown in Figure 3. The solid line corresponds to the case without the surfactant, and the dashed line is plotted in the case with added surfactant on the surface ( N = 10 6 and L = 0.003 ). In all the calculations, here we choose S = 1 , which corresponds to the definition ϵ Σ 1 / 2 .
In the region where s 1 = O ( δ ) , we can denote s 1 = δ s ¯ 1 and relegate this term to the third order of expansion. In the vicinity of the line s 1 = 0 , we can set τ 1 = M 1 = 0 . Then, we obtain equations:
τ 2 A l = κ 0 A l + κ 1 | A l | 2 A l + κ 2 ( | A m | 2 + | A n | 2 ) A l + s ¯ 1 A m * A n * , { l , m , n } = { 1 , 2 , 3 } , { 2 , 3 , 1 } , { 3 , 1 , 2 } ,
where κ 1 is the same as in Equations (14) and (15), and κ 2 characterizes the interaction of rolls turned to the angle 2 π / 3 with each other. The growth of amplitudes is saturated if κ 1 < 0 and κ 1 + 2 κ 2 < 0 .

3.3. Rhombic Lattice

The patterns of a rhombic lattice are less widespread than rolls/squares or hexagons. However, they are known in different physical processes, such as [51,52,53,54]. These patterns are the result of the interaction of two external perturbations, oriented by angle θ one to another. In other words, we have an external one in rolls with respect to disturbances such that their wavevectors form angle θ with the basic wavevector of the roll. Thus, the solutions at the leading order of the variables can belong to a rhombic lattice in the Fourier space:
h 1 = A 1 ( τ 2 ) e i K c X + A 2 ( τ 2 ) e i ( K c X cos θ + K c Y sin θ ) + c . c . , f 1 = α 1 h 1 , γ 1 = α 2 h 1 .
The angle θ is an additional parameter of the problem. Without a loss of generality, we take for consideration θ < π / 2 , but θ does not equal π / 3 . The rhombic lattice at θ = π / 2 becomes a square one. Note that patterns with wavevectors on the rhombic lattice in the Fourier space look like rectangles in the real space.
Making the same procedures as in the case of the square lattice, we obtain for the complex amplitudes A 1 and A 2 the system of equations similar to (14) and (15):
A 1 τ 2 = κ 0 A 1 + κ 1 | A 1 | 2 A 1 + κ ˜ 2 ( θ ) | A 2 | 2 A 1 ,
A 2 τ 2 = κ 0 A 2 + κ 1 | A 2 | 2 A 2 + κ ˜ 2 ( θ ) | A 2 | 2 A 2 .
Here, κ 2 = κ ˜ 2 ( π / 2 ) . The growth of amplitudes on the rhombic lattice is saturated if κ 1 < 0 and κ 1 + κ ˜ 2 < 0 .
Figure 4 presents the pattern selection on the rhombic lattice for three values of angle θ . Figure 4a shows the case θ = 0.45 π . Here, the angle is closer to π / 2 and the map is similar to the patterns on the square lattice. In Figure 4b is the case θ = 0.4 π , and in Figure 4c, θ = 0.25 π .
Each panel includes two cases of different N; the black color lines are for N = 0 and the blue ones are for N = 10 6 . It is shown that surfactant diminishes the regions of stable rhombi. The calculations show that the most dangerous perturbations are directed by θ = π / 2 .

4. Modulational Instability of Stationary Rolls—The Case of 1D Disturbances

The stability analysis carried out in the previous section is incomplete because it does not include the stability of periodic patterns with respect to disturbances with wavenumbers different from K c . The most important class of disturbances are those with wavenumbers close to K c . Below, we present the analysis of the instability of roll patterns with respect to one-dimensional disturbances near the instability threshold ( M = M 0 + δ 2 M 2 , M 0 = M m ( K c ) ) carried out in [31]. In that case, the motionless state is unstable with respect to disturbances with the wavenumbers in the interval O ( δ ) around K c . Using the Newell–Whitehead–Segel approach, two spatial scales are introduced: X corresponding to the main wavelength of the pattern, and X 1 = δ X corresponding to the spatial modulation of the pattern. The time variable is rescaled as τ 2 = δ 2 τ . Thus, the solution at the leading order can be written as
h 1 = A ( X 1 , τ 2 ) e i K c X 0 + A ( X 1 , τ 2 ) * e i K c X 0 , f 1 = α 1 h 1 , γ 1 = α 2 h 1
(coefficients α 1 and α 2 are the same as in the previous section). At the next order, the solution can be written as
h 2 = h ¯ ( X 1 , τ 2 ) + B 1 ( X 1 , τ 2 ) e 2 i K c X 0 + B 1 ( X 1 , τ 2 ) * e 2 i K c X 0 , f 2 = f ¯ ( X 1 , τ 2 ) + B 2 ( X 1 , τ 2 ) e 2 i K c X 0 + B 2 ( X 1 , τ 2 ) * e 2 i K c X 0 , γ 2 = γ ¯ ( X 1 , τ 2 ) + B 3 ( X 1 , τ 2 ) e 2 i K c X 0 + B 3 ( X 1 , τ 2 ) * e 2 i K c X 0 .
The heat-transfer equation at the second order gives f ¯ h ¯ = K c 2 β | A | 2 . The calculation is performed taking into account that convection does not change the total volume of the liquid and the total amount of the surfactant one has < h ¯ > = 0 , < γ ¯ > = 0 , ( < > is averaging over X 1 ).
Finally, one obtains the following set of equations for the amplitude of the rolls A, for the large-scale disturbances of the surface deformation h 0 , and for surfactant concentration γ 0 :
τ 2 A = κ 0 A + κ 1 | A | 2 A + μ 1 X 1 X 1 A + μ 2 h ¯ A + μ 3 γ ¯ A ,
τ 2 h ¯ = a 1 X 1 X 1 | A | 2 + a 2 X 1 X 1 h ¯ + a 3 X 1 X 1 γ ¯ ,
τ 2 γ ¯ = b 1 X 1 X 1 | A | 2 + b 2 X 1 X 1 h ¯ + b 3 X 1 X 1 γ ¯ .
The nature of the terms and the expressions for the coefficients are described in [31]. The coefficients μ i ( i = 1 , 2 , 3 ) are cumbersome and not presented here. Coefficients a i and b i are described in Appendix B.
Stationary solution of Equations (24)–(26) without deformation and non-homogeneity of surfactant is a family of rolls A 0 = q 2 μ 1 κ 0 κ 1 e i q X 1 , h ¯ = 0 , γ ¯ = 0 ( κ 1 < 0 ) . Here, q is the deviation of the wavenumber from K c ( q 2 < q 0 2 = κ 0 / μ 1 ).
To consider the stability of the rolls, we disturb their amplitude and the phase in the form of normal modes. Obtaining the dispersion relation for the growth rate, λ , of the side-band instability, we find in the limit of longwave modulations (assuming λ = λ 0 + Λ k 2 + o ( k 2 ) ) that the growth rate of the mode related to the amplitude modulation, λ 0 , is negative and the other three eigenvalues tend to zero when k 0 . This equation for Λ describes the interaction of three Goldstone modes corresponding to the definite symmetries of the problem—phase disturbances, surface deformation, and surfactant concentration disturbances.
The analysis of the dispersion relation in the case without surfactant gives a quadratic equation for eigenvalue Λ 1 , 2 . Fixing parameters G and β (remember, S = 1 ), we can follow the evolution of the modulated rolls as a function of q 2 .
Figure 5 presents the stability map of the stationary rolls in the case without surfactant. The black solid line here shows the boundary between supercritical and subcritical domains for the rolls. We analyzed the modulational instability only in the region of supercritical rolls and found three regions:
  • Region 1. Here, all the rolls are unstable with respect to modulation.
  • Region 2. Here, the rolls are stable within the interval q 2 < q m 2 and monotonically unstable for q m 2 < q 2 < q 0 2 (monotonic Eckhaus instability).
  • Region 3. Here, the rolls are stable within interval q 2 < q o s c 2 , oscillatory unstable for q 2 slightly above q o s c 2 , and oscillatory or monotonically unstable for q o s c 2 < q 2 < q 0 2 .
The boundary of the monotonic Eckhaus instability is found as
q m 2 = κ 0 μ 1 a 2 κ 1 a 1 μ 2 3 a 2 κ 1 a 1 μ 2 .
For q o s c 2 , we found
q o s c 2 = κ 0 μ 1 a 2 κ 1 + κ 1 μ 1 a 1 μ 2 a 2 κ 1 + 3 κ 1 μ 1 a 1 μ 2 .
Even a small addition of surfactant changes the stability interval. The boundary of monotonic Eckhaus instability is
q m 2 = q 0 2 1 + C 3 + C ,
where q 0 2 = κ 0 / μ 1 is the width of the existence interval for periodic solutions and
C = ( a 3 b 1 a 1 b 3 ) μ 2 ( a 2 b 1 a 1 b 2 ) μ 3 ( a 2 b 3 a 3 b 2 ) κ 1 .
In the case without surface deformations and disturbances of surfactant concentration, we have C = 0 , and we obtain the classical result q m 2 = ( 1 / 3 ) q 0 2 .
Figure 6 presents the stability maps for N = 10 8 (a) and N = 10 7 (b). The shrinking of the region of stable supercritical rolls means that with an increasing surfactant concentration, the instability becomes subcritical in the larger region of the parameter.

5. Modulational Instability of Stationary Rolls—The Case of 2D Disturbances

Consider now the 2D modulation of the rolls. The roll pattern is typically selected by the analysis of competition of disturbances characterized by difference | K n K m | O ( 1 ) , i.e., the angle between these two wavevectors is O ( 1 ) . In the case of the problem with rotational symmetry, the set of wavevectors corresponding to unstable modes is a ring which contains a continuum of wavevectors with | K 2 1 | < δ , see Figure 7.
Let us take a roll pattern with K 0 = ( K c , 0 ) and consider the disturbances with K = K 0 + Δ around that wavevector. The most unstable modes have K = K c . When M a M a c O ( δ 2 ) , the instability region of the motionless state is a ring | | K | K c | = O ( δ ) .
Because
| K | = | K + Δ | = ( K c + Δ x ) 2 + Δ y 2 = K c + Δ x + Δ y 2 2 K c + ,
the instability of the motionless state takes place with respect to the disturbances with the wavevectors in the region
| Δ x | = O ( δ ) , | Δ y | = O ( δ 1 / 2 ) .
In addition, the coupling between the longwave roll modulation, the modulations of the surface distortion h, and the distribution of the surfactant concentration γ has some specific features in the region of | K y | = O ( δ ) . Therefore, one has to consider two regions of transversal modulation characterized by rescaled variables Y 1 / 2 = δ 1 / 2 Y and Y 1 = δ Y .

5.1. Transversal Modulation of Rolls, Y 1 = δ Y

Introduce the following rescaling of spatial and time coordinates
X 0 = X , X 1 = δ X , Y 1 = δ Y , τ 2 = δ 2 τ .
As in the previous section, applying the approach of Newell–Whitehead–Segel, one obtains the following set of evolution equations for amplitude A, for the surface deformation h 0 , and for the disturbance of surfactant concentration γ 0 :
τ 2 A = κ 0 A + κ 1 | A | 2 A + μ 1 X 1 X 1 A + μ 2 h 0 A + μ 3 γ 0 A ,
τ 2 h 0 = ( a 1 X 1 X 1 + a ˜ 1 Y 1 Y 1 ) | A | 2 + a 2 ( X 1 X 1 + Y 1 Y 1 ) h 0 + a 3 ( X 1 X 1 + Y 1 Y 1 ) γ 0 ,
τ 2 γ 0 = ( b 1 X 1 X 1 + b ˜ 1 Y 1 Y 1 ) | A | 2 + b 2 ( X 1 X 1 + Y 1 Y 1 ) h 0 + b 3 ( X 1 X 1 + Y 1 Y 1 ) γ 0 .
Coefficients κ 1 , μ 1 , μ 2 , and μ 3 are cumbersome and not given here; others are presented in Appendix A.
In the framework of the linear stability theory, the perturbations of the roll amplitude and its phase, surface deviation, and surfactant of concentration are taken proportional to e λ τ 2 + i K X X 1 + i K Y Y 1 . Consider a longwave modulation of the rolls and
| k | 2 = k 2 = K X 2 + K Y 2 = K X 2 ( 1 + χ 2 ) 1 .
Here, χ = K Y / K X . When | K Y | = 0 , we have the case described in the previous section. The boundary of the monotonic Eckhaus instability is described by (29). In the case of transversal modulation, | K Y | 0 , the parameter C in (29) is replaced by C ˜ = C + C ˜ χ 2 1 + χ 2 . Here, C ˜ is obtained from C by the replacement a 1 a ˜ 1 and b 1 b ˜ 1 . Because for any χ 2 , C ˜ is between C and C ˜ , it is sufficient to compare the criterion with the new one,
q 2 < q ˜ m 2 = q 0 2 1 + C ˜ 3 + C ˜ .
The rolls are stable with respect to disturbances with small K X 2 , K Y 2 , if q 2 < min ( q 0 2 , q ˜ m 2 ) . The stability maps have additional parameter, χ . We start plot for case without the surfactant, see Figure 8. When χ grows, the stability of stationary rolls increase, i.e., the longitudinal disturbances are more dangerous than disturbances with inclined wavevector. At χ 2 > 5.2 , the region “1” disappears.
Figure 9 presents the stability maps for two different values of N: panel (a) N = 10 8 , panel (b) N = 10 7 for different χ .

5.2. Transversal Modulation of Rolls, Y 1 / 2 = δ 1 / 2 Y

Here, we rescale coordinates as
X 0 = X , X 1 = δ X , Y 1 / 2 = δ 1 / 2 Y , τ 2 = δ 2 τ .
Now, Equations (34) and (35) in the third order become
Y 1 / 2 2 ( a ˜ 1 | A | 2 + a 2 h 0 + a 3 γ 0 ) ,
Y 1 / 2 2 ( b ˜ 1 | A | 2 + b 2 h 0 + b 3 γ 0 ) .
Thus, the system of algebraic equations for surface distortion h 0 and for distribution of the surfactant concentration γ 0 can be written as
a 2 h 0 + a 3 γ 0 = a ˜ 1 ( | A | 2 < | A | 2 > ) ,
b 2 h 0 + b 3 γ 0 = b ˜ 1 ( | A | 2 < | A | 2 > ) ,
where < > denotes averaging over both X 1 and Y 1 / 2 . Solving this system, we find the solutions for surface distortion h 0 and distribution of the surfactant concentration γ 0 can be written as
h 0 = c h ( | A | 2 < | A | 2 > ) , γ 0 = c γ ( | A | 2 < | A | 2 > ) .
Here,
c h = a 3 b ˜ 1 a ˜ 1 b 3 a 2 b 3 b 2 a 3 , c γ = a ˜ 1 b 2 a 2 b ˜ 1 a 2 b 3 b 2 a 3 .
In this case, we can rewrite the amplitude equation as the following modified Newell–Whitehead–Segel equation:
τ 2 A = κ 0 A + κ 1 | A | 2 A + κ ˜ 1 ( | A | 2 < | A | 2 > ) A + μ 1 X 1 i 2 2 K c Y 1 / 2 2 2 A .
Here, κ ˜ 1 = μ 2 c h + μ 3 c γ . Linearizing Equation (43) and substituting the disturbances of the rolls amplitude and its phase in the form exp ( λ τ 2 + i ( K X X 1 + K Y Y 1 / 2 ) , we obtain quadratic equation for the growth rate λ .
In the case when the perturbations vary only in Y direction (parallel to the rolls), K X = 0 the dispersion relation gives the growth rate for the phase perturbation
λ = K Y 2 q μ 1 K c K Y 4 μ 1 4 K c 2 .
Here, when q < 0 , we have the zigzag instability for small | K Y | ( K Y 2 < 4 K c q ). The zigzag instability caused by the phase disturbances is not influenced by κ ˜ 1 . The origin of the zigzag instability is fully “geometrical” for K < K c , there is a resonance of neutral disturbances with wavevectors ( ± K , K Y ) with zero linear growth rate, K 2 + K Y 2 = K c 2 :
( K , K Y ) + ( 2 K , 0 ) = ( K , K Y ) .

6. Distortions of Hexagons

Let us discuss now the longwave distortion of hexagons:
h 1 = j = 1 3 A j ( X 1 , τ 1 , τ 2 ) e i K c n j · X 0 + c . c , f 1 = α 1 h 1 , γ 1 = α 2 h 1 ,
where X 0 = X , X 1 = δ X . Distorted non-equilateral hexagons have been observed in chemical systems [55] and considered in theoretical models [56,57].
The general form of the amplitude equations for distorted patterns described by amplitudes A l ( l = 1 , 2 , 3 ) is as follows [58,59]:
τ 2 A l = κ 0 A l + κ 1 | A l | 2 A l + κ 2 ( | A m | 2 + | A n | 2 ) A l + μ 1 ( n l · ) 2 A l + μ 2 h ¯ A l + + μ 3 γ ¯ A l + s ¯ 1 A m * A n * + i s 2 [ A m * ( t n · ) A n * A n * ( t m · ) A m * ] + i s 3 [ A n * ( n m · ) A m * + + A m * ( n n · ) A n * ] ,
τ 2 h ¯ = j = 1 3 [ a 1 ( n j · ) 2 + a ˜ 1 ( t j · ) 2 ] | A j | 2 + a 2 2 h ¯ + a 3 2 γ ¯ ,
τ 2 γ ¯ = j = 1 3 [ b 1 ( n j · ) 2 + b ˜ 1 ( t j · ) 2 ] | A j | 2 + b 2 2 h ¯ + b 3 2 γ ¯ .
Here, X 1 . New coefficients are: a ˜ 1 = 1 3 ( G + K c 2 S ) K c 2 M m 2 β and b ˜ 1 = 1 36 L ( G + K c 2 S ) ( 12 L G K c 2 S ) K c 2 M m β . The terms with the coefficient s 2 describe the change in the nonlinear interaction by rotation of the wavevectors of rolls forming the hexagons, while the terms with s 3 correspond to that by dilatations of hexagons. The values of the coefficients s 2 and s 3 are computed directly from the original Equations (7)–(9), similarly to [60,61]. The terms μ 2 h ¯ A l and μ 3 γ ¯ A l describe the influence of surface deformation and surfactant concentration perturbation created by the pattern modulation on the growth rate of Marangoni instability.

7. Conclusions

In this paper, a review is presented of the recent works on the longwave Marangoni convection of the non-isothermal liquid layer covered by an insoluble surfactant. This analysis includes two different ranges of disturbances: one is with wavevectors k B i 1 / 4 and one with k B i 1 / 2 ( B i is the Biot number). While the first interval of wavevectors is well known and widely investigated, the second one, k B i 1 / 2 , is less known. Recent works elucidate the pattern selection at the interval of the disturbance wavenumbers: rolls/squares, hexagons, and rhombi were found. In the present work, we studied how the concentration of an insoluble surfactant, characterized by the elasticity parameter N, influences the roll stability, as well as that of rectangles and hexagons. The last pattern must be still investigated in more detail. The analysis of external instability is performed at a fixed value of the critical wavevector at which the monotonic instability mode of the Marangoni convection appears.
The special part of the paper is devoted to the modulational instability, when the 1D or 2D large-scale perturbations disturb the stable roll patterns, which are formed as a result of the stationary Marangoni convection in the liquid layer with/without insoluble surfactant. The result of this modulation is the appearance of the Eckhaus instability. The analysis is performed in the neighborhood of the critical wavenumber. The difference between the 1D and 2D modulation is revealed. In both cases, the stability maps are drawn in the plane of the Galileo number, G, and the Biot number, B i .
The existence of three regions of the supercritical rolls is found. Region “1” is where all the rolls are unstable with respect to modulation. Region “2” is where the rolls are stable within the interval q 2 < q m 2 and monotonically unstable for q m 2 < q 2 < q 0 2 . Region “3” is where the rolls are stable within the interval q 2 < q o s c 2 , oscillatory unstable for q 2 slightly above q o s c 2 , and oscillatory or monotonically unstable for any q o s c 2 < q 2 < q 0 2 . The zigzag instability boundary does not depend on the concentration of the insoluble surfactant.
The paper includes the known results published earlier in different papers, as well as new results never yet published. The authors anticipate that the discussed theory can motivate more experiments on the investigation of pattern formation in the system with surfactants. To estimate the parameters of potential experiments, we take the water layer of a thickness of 0.05 mm. As shown in [13], the longwave Marangoni convection is possible if B i Σ < 72 . Assuming the heat-transfer coefficient q 10 ( W / m 2 K ) and the thermal conductivity is 0.5 ( W / m K ) , one obtains B i 0.001 . The Galileo number corresponds to G 10 and the inverse capillary number Σ 3 · 10 4 .
The authors hope that this paper will be useful for experts in colloid physics and in surface active agents (surfactants).

Author Contributions

A.B.M.: software, formal analysis, investigation, writing—original draft, writing—review and editing, and visualization; A.A.N.: conceptualization, methodology, validation, writing—original draft and review and editing. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Acknowledgments

The authors acknowledge the support by the Israel Science Foundation (Grant No. 843/18).

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Basic Wavevectors for the Patterns

In this paper, we considered the patterns characterized by three different sets of basic wavevectors K i in the Fourier space.
For a square lattice, the basic wavevectors are perpendicular to another, Figure A1a. In the rhombic lattice, the basic wavevectors are oriented one to another by angle θ . In the real space, the rhombic structures look like rectangles, Figure A1b. The most typical structures are hexagonal patterns, Figure A1c, with three basic wavevectors oriented by 2 π / 3 . That pattern is supported by the resonance interaction of disturbances forming a triad, K 1 + K 2 + K 3 = 0 .
Figure A1. Basic wavevectors of disturbances. Panel (a) square/rolls lattice, (b) rhombic lattice, (c) hexagonal lattice.
Figure A1. Basic wavevectors of disturbances. Panel (a) square/rolls lattice, (b) rhombic lattice, (c) hexagonal lattice.
Colloids 06 00053 g0a1

Appendix B. Coefficients of (33)–(35)

a 1 = 1 3 ( G + 7 K 2 S ) K 2 M m 2 β , a ˜ 1 = 1 3 ( G + K 2 S ) K 2 M m 2 β , b 1 = 1 36 L [ G ( 12 L G ) + 4 K 2 ( G + 21 L ) S + 5 K 4 S 2 ] K 2 M m β , b ˜ 1 = ( G + K 2 S ) ( G 12 L + K 2 S ) 36 L K 2 M m β a 2 = G 3 , a 3 = N 2 , b 2 = G 2 , b 3 = L + N .

References

  1. Bénard, H. Les tourbillons cellulaires dans une nappe liquide. Rev. Gen. Sci Pures Appl. 1900, 11, 1261–1271. [Google Scholar]
  2. Block, M.J. Surface tension as a cause of Bénard cells and surface deformation in a liquid film. Nature 1956, 176, 650–651. [Google Scholar] [CrossRef]
  3. Pearson, J.R.A. On convection cells induced by surface tension. J. Fluid Mech. 1958, 4, 489–500. [Google Scholar] [CrossRef]
  4. Turing, A.M. The chemical basis of morphogenesis. Phil. Trans. R. Soc. Lond. B 1952, 237, 37–72. [Google Scholar]
  5. De Kepper, P.; Castets, V.; Dulos, E.; Boissonade, J. Turing-type chemical patterns in the chlorite-iodide malonic acid reaction. Physica D 1991, 49, 161–169. [Google Scholar] [CrossRef]
  6. Normand, C.; Pomeau, Y.; Velarde, M.G. Convective instability—Physicist’s approach. Rev. Mod. Phys. 1977, 49, 581–624. [Google Scholar] [CrossRef]
  7. Colinet, P.; Nepomnyashchy, A. (Eds.) Pattern Formation at Interfaces; Springer: Vienna, Austria, 2010; p. 304. [Google Scholar]
  8. Shklyaev, S.V.; Nepomnyashchy, A.A. Longwave Instabilities and Patterns in Fluids; Birkhäuser: New York, NY, USA, 2017; p. 453. [Google Scholar]
  9. Scriven, L.E.; Sternling, C.V. On cellular convection driven by surface-tension gradients: Effects of mean surface tension and surface viscosity. J. Fluid Mech. 1964, 19, 321–340. [Google Scholar] [CrossRef]
  10. Sivashinsky, G.I. Large cells in nonlinear Marangoni convection. Physica D 1982, 4, 227–235. [Google Scholar] [CrossRef]
  11. Knobloch, E. Pattern selection in long-wavelength convection. Physica D 1990, 41, 450–479. [Google Scholar] [CrossRef]
  12. Davis, S.H. Rapture of thin liquid films. In Waves on Fluid Interfaces; Mayer, R.E., Ed.; Academic Press: New York, NY, USA, 1983; p. 291. [Google Scholar]
  13. Shklyaev, S.; Alabuzhev, A.A.; Khenner, M. Long-wave Marangoni convection in a thin film heated from below. Phys. Rev. E 2012, 85, 016328. [Google Scholar] [CrossRef]
  14. Hoyle, R.B. Pattern Formation: An Introduction to Methods; Cambridge University Press: Cambridge, UK, 2006. [Google Scholar]
  15. Lohse, D. Fundamental fluid dynamics challenges in inkjet printing. Annu. Rev. Fluid Mech. 2022, 54, 349–382. [Google Scholar] [CrossRef]
  16. Singh, M.; Haverinnen, H.M.; Dhagat, P.; Jabbour, G.E. Inkjet printing-process and its applications. Adv. Mater. 2010, 22, 673–685. [Google Scholar] [CrossRef]
  17. Kim, J. Spray cooling heat transfer: The state of the art. Int. J. Heat Fluid Flow 2007, 28, 753–767. [Google Scholar] [CrossRef]
  18. Smalyukh, I.; Zribi, O.; Butler, J.; Lavrentovich, O.; Wong, G. Structure and dynamics of liquid crystalline pattern formation in drying droplets of DNA. Phys. Rev. Lett. 2006, 96, 177801. [Google Scholar] [CrossRef]
  19. Ziemelis, K. The future of microelectronics. Nature 2000, 406, 1021. [Google Scholar]
  20. Espinosa, F.F.; Shapiro, A.H.; Fredberg, J.J.; Kamm, R.D. Spreading of exogenous surfactant in an airway. J. Appl. Physiol. 1993, 75, 2028–2039. [Google Scholar] [CrossRef] [PubMed]
  21. Deegan, R.D. Pattern formation in drying drops. Phys. Rev. E 2000, 61, 475–485. [Google Scholar] [CrossRef]
  22. van Gaalen, R.T.; Wijshoff, H.M.A.; Kuerten, J.G.M.; Diddens, C. Competition between thermal and surfactant-induced Marangoni flow in evaporating sessile droplets. J. Coll. Int. Sci. 2022, 622, 892–903. [Google Scholar] [CrossRef]
  23. Karapetsas, G.; Sahu Kirti, C.; Matar, O.K. Evaporation of sessile droplets laden with particles and insoluble surfactants. Langmuir 2016, 32, 6871–6881. [Google Scholar] [CrossRef]
  24. Bhardwaj, R.; Fang, X.; Attinger, D. Pattern formation during the evaporation of a colloidal nanoliter drop: A numerical and experimental study. New J. Phys. 2009, 11, 075020. [Google Scholar] [CrossRef]
  25. Maillard, M.; Motte, L.; Ngo, A.T.; Pileni, M.P. Rings and hexagons of nanocrystals: A Marangoni effect. J. Phys. Chem. B 2000, 104, 11871–11877. [Google Scholar] [CrossRef]
  26. Sauleda, M.L.; Chu, H.C.W.; Tilton, R.D.; Garoff, S. Surfactant driven Marangoni spreading in the presence of predeposited insoluble surfactant monolayers. Langmuir 2021, 37, 3309–3320. [Google Scholar] [CrossRef] [PubMed]
  27. Wodlei, F.; Sebilleau, J.; Magnaudet, J.; Pimienta, V. Marangoni-driven flower-like patterning of a evaporating drop spreading on a liquid substrate. Nat. Commun. 2018, 9, 820. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Shao, X.; Duan, F.; Hou, Y.; Zhong, X. Role of surfactant in controlling the deposition pattern of a particle-laden droplet: Fundamentals and strategies. Adv. Coll. Int. Sci. 2020, 275, 102049. [Google Scholar] [CrossRef]
  29. Newell, A.C.; Whitehead, J.A. Finite bandwidth, finite amplitude convection. J. Fluid Mech. 1969, 38, 279–303. [Google Scholar] [CrossRef]
  30. Segel, L.A. Distant side-walls cause slow amplitude modulation of cellular convection. J. Fluid Mech. 1969, 38, 203–224. [Google Scholar] [CrossRef]
  31. Mikishev, A.B.; Nepomnyashchy, A.A. Patterns and their large-scale distortions in Marangoni convection with insoluble surfactant. Fluids 2021, 6, 282. [Google Scholar] [CrossRef]
  32. Gershuni, G.Z.; Zhukhovitsky, E.M. Convective Stability of Incompressible Fluid; Keter: Jerusalem, Israel, 1976. [Google Scholar]
  33. Rosen, M.J.; Zhen, H.Z. Enchancement of wetting properties of water-insoluble surfactant via solubilization. J. American Oil Chem. Soc. 1993, 70, 65–68. [Google Scholar] [CrossRef]
  34. Miller, R.; Aksenenko, E.V.; Fainerman, V.B. Dynamic interfacial tension of surfactant solutions. Adv. Colloid Interface Sci. 2017, 247, 115–129. [Google Scholar] [CrossRef]
  35. Levich, V.G. Physicochemical Hydrodynamics; Prentice Hall: Englewood Cliffs, NJ, USA, 1962. [Google Scholar]
  36. Wong, H.; Rumschitski, C.; Maldarelli, C. On the surfactant mass balance at a deforming fluid interface. Phys. Fluids 1996, 8, 3203–3204. [Google Scholar] [CrossRef]
  37. Mikishev, A.B.; Nepomnyashchy, A.A. Instabilities in evaporating liquid layer with insoluble surfactant. Phys. Fluids 2013, 25, 054109. [Google Scholar] [CrossRef]
  38. Stone, H. A simple derivation of the time-dependent convective diffusion equation for surfactant transport along a deformed interface. Phys. Fluids A 1990, 2, 111–112. [Google Scholar] [CrossRef]
  39. Mikishev, A.B.; Nepomnyashchy, A.A. Long-wavelength Marangoni convection in a liquid layer with insoluble surfactant: Linear theory. Microgravity Sci. Technol. 2010, 22, 415–423. [Google Scholar] [CrossRef]
  40. Mikishev, A.B.; Nepomnyashchy, A.A. Amplitude equations for large-scale Marangoni convection in a liquid layer with insoluble surfactant on deformable free surface. Microgravity Sci. Technol. 2011, 23, 59–63. [Google Scholar] [CrossRef]
  41. Mikishev, A.B.; Friedman, B.A.; Nepomnyashchy, A.A. Generation of transverse waves in a liquid layer with insoluble surfactant subjected to temperature gradient. Fluid Dyn. Res. 2016, 48, 061403. [Google Scholar] [CrossRef]
  42. Nepomnyashchy, A.A.; Simanovskii, I.B. Thermocapillary convection in two-layer systems in the presence of a surface active agent at the interface. Fluid Dyn. 1986, 21, 169–179. [Google Scholar] [CrossRef]
  43. Nepomnyashchy, A.A.; Simanovskii, I.B. Occurance of thermocapillary convection in a two-layer system in presence of a soluble surfactant. Fluid Dyn. 1988, 23, 302–306. [Google Scholar] [CrossRef]
  44. Nepomnyashchy, A.A.; Simanovskii, I.B. Onset of oscillatory convection in a two-layer system due to the presence of surfactant at the interface. Sov. Phys. Dokl. 1989, 34, 420–422. [Google Scholar]
  45. Mikishev, A.B.; Nepomnyashchy, A.A. Weakly nonlinear analysis of long-wave Marangoni convection in a liquid layer covered by insoluble surfactant. Phys. Rev. Fluids 2019, 4, 094002. [Google Scholar] [CrossRef]
  46. Murray, J.D. Mathematical Biology; Springer: Berlin, Germany, 1989; 767p. [Google Scholar]
  47. Wollkind, D.J.; Sriranganathan, R.; Oulton, D.B. Interfacial patterns during plane front alloy solidification. Physica D 1984, 12, 215–240. [Google Scholar] [CrossRef]
  48. Edwards, W.S.; Fauve, S. Patterns and quasi-patterns in the Faraday experiment. J. Fluid Mech. 1994, 278, 123–148. [Google Scholar] [CrossRef]
  49. Arecchi, F.T.; Boccaletti, S.; Ducci, S.; Pampoloni, E.; Ramazza, P.L.; Residori, S. The liquid crystal light valve with optical feed back: A case study in pattern formation. J. Nonlinear Opt. Phys. Mat. 2000, 9, 183–204. [Google Scholar] [CrossRef]
  50. Abou, B.; Westfreid, J.-E.; Roux, S. The normal field instability in ferrofluids: Hexagon-square transition mechanism and wavenumber selection. J. Fluid Mech. 2000, 416, 217–237. [Google Scholar] [CrossRef]
  51. Ouyang, Q.; Gunaratne, G.H.; Swinney, H.L. Rhombic patterns: Broken hexagonal symmetry. Chaos 1993, 3, 707–711. [Google Scholar] [CrossRef]
  52. Pansuwan, A.; Rattanakul, C.; Lenbury, Y.; Wollkind, D.J.; Harrison, L.; Rajapakse, I.; Cooper, K. Nonlinear stability analyses of pattern formation on solid surfaces during ion-sputtered erosion. Math. Comput. Modell. 2005, 41, 939–964. [Google Scholar] [CrossRef]
  53. Skeldon, A.C.; Silber, M. New stability results for patterns in a model of long-wavelength convection patterns. Physica D 1998, 122, 117–133. [Google Scholar] [CrossRef]
  54. Franz, J.; Grünebaum, J.; Schäfer, M.; Mulac, D.; Kehfeldt, F.; Langer, K.; Kramer, A.; Riethmüller, C. Rhombic organization of microvilli domains found in a cell model of the human intenstine. PLoS ONE 2018, 13, e0189970. [Google Scholar]
  55. Ouyang, Q.; Swinney, H.L. Transition from a uniform state to hexagonal and striped Turing patterns. Nature 1991, 352, 610–612. [Google Scholar] [CrossRef]
  56. Gunaratne, G.H.; Ouyang, Q.; Swinney, H.L. Pattern formation in the presence of symmetries. Phys. Rev. E 1994, 50, 2802–2820. [Google Scholar] [CrossRef]
  57. Nuz, A.; Nepomnyashchy, A.A.; Golovin, A.A.; Hari, A.A.; Pismen, L.M. Stability of rolls and hexagonal patterns in non-potential systems. Physica D 2000, 135, 233–262. [Google Scholar] [CrossRef]
  58. Brand, H.R. Envelope equations near the onset of a hexagonal pattern. Progr. Theor. Phys. Suppl. 1989, 99, 442–449. [Google Scholar] [CrossRef]
  59. Kuznetsov, E.A.; Nepomnyashchy, A.A.; Pismen, L.M. New amplitude equation for Boussinesq convection and non-equilateral hexagonal patterns. Phys. Let. A 1995, 2051, 261–265. [Google Scholar] [CrossRef]
  60. Golovin, A.A.; Nepomnyashchy, A.A.; Pismen, L.M. Nonlinear evolution and secondary instabilities of Marangoni convection in a liquid-gas system with deformable interface. J. Fluid Mech. 1997, 341, 317–341. [Google Scholar] [CrossRef]
  61. Bragard, J.; Velarde, M.G. Bénard-Marangoni convection: Planforms and related related theoretical predictions. J. Fluid Mech. 1998, 368, 165–194. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Schematic of the problem.
Figure 1. Schematic of the problem.
Colloids 06 00053 g001
Figure 2. Pattern selection domains for the square lattice. The solid lines are for κ 1 = 0 , the dashed lines are for κ 1 + κ 2 = 0 , and dot-dashed lines are for κ 2 κ 1 = 0 . Here, the red-colored set (# 1 ) is for N = 0 , the blue-colored set (# 2) is for N = 10 6 , and the black set (#3) is for N = 10 5 . The domains are marked: “R” for rolls, “S” for squares. Panel (b) is detailed panel (a) for small β .
Figure 2. Pattern selection domains for the square lattice. The solid lines are for κ 1 = 0 , the dashed lines are for κ 1 + κ 2 = 0 , and dot-dashed lines are for κ 2 κ 1 = 0 . Here, the red-colored set (# 1 ) is for N = 0 , the blue-colored set (# 2) is for N = 10 6 , and the black set (#3) is for N = 10 5 . The domains are marked: “R” for rolls, “S” for squares. Panel (b) is detailed panel (a) for small β .
Colloids 06 00053 g002
Figure 3. Pattern selection for regular hexagons ( q 0 = 0 ). Solid line: s 1 ( G , β ) = 0 in the case without surfactant ( N = 0 ), the dashed line: in the case with surfactant ( N = 10 6 , L = 0.003 ). Up-hexagons ( H + ) are selected in the region above the line s 1 = 0 , and down-hexagons ( H ) are selected below that line.
Figure 3. Pattern selection for regular hexagons ( q 0 = 0 ). Solid line: s 1 ( G , β ) = 0 in the case without surfactant ( N = 0 ), the dashed line: in the case with surfactant ( N = 10 6 , L = 0.003 ). Up-hexagons ( H + ) are selected in the region above the line s 1 = 0 , and down-hexagons ( H ) are selected below that line.
Colloids 06 00053 g003
Figure 4. Pattern selection domains for the rhombic lattice. The solid lines are for κ 1 = 0 , the dashed lines are for κ 1 + κ ˜ 2 ( θ ) = 0 . Here, the black color set (#1 ) is for N = 0 , the blue-colored set (#2 ) is for N = 10 6 . The stability domains are marked: “1Rh” is for N = 0 , “2Rh” for N = 10 6 . Panel (a): θ = 0.45 π , panel (b): θ = 0.4 π , and panel (c): θ = 0.25 π .
Figure 4. Pattern selection domains for the rhombic lattice. The solid lines are for κ 1 = 0 , the dashed lines are for κ 1 + κ ˜ 2 ( θ ) = 0 . Here, the black color set (#1 ) is for N = 0 , the blue-colored set (#2 ) is for N = 10 6 . The stability domains are marked: “1Rh” is for N = 0 , “2Rh” for N = 10 6 . Panel (a): θ = 0.45 π , panel (b): θ = 0.4 π , and panel (c): θ = 0.25 π .
Colloids 06 00053 g004
Figure 5. Stability map of the stationary rolls in the case without surfactant.
Figure 5. Stability map of the stationary rolls in the case without surfactant.
Colloids 06 00053 g005
Figure 6. Stability maps for N = 10 8 (a) and N = 10 7 (b). Other parameters L = 0.003 and S = 1 . Regions “1”, “2”, and “3” correspond to the same named in Figure 5.
Figure 6. Stability maps for N = 10 8 (a) and N = 10 7 (b). Other parameters L = 0.003 and S = 1 . Regions “1”, “2”, and “3” correspond to the same named in Figure 5.
Colloids 06 00053 g006
Figure 7. Ring of unstable wavevectors.
Figure 7. Ring of unstable wavevectors.
Colloids 06 00053 g007
Figure 8. Stability map of the stationary rolls in the case without surfactant. The solid black line is the boundary line between supercritical and subcritical rolls (“subcr. R”). Region “1” (below the dashed line) is the region of unstable rolls; in region “2”, rolls are stable within q 2 < q m 2 and monotonically unstable for q m 2 < q 2 < q 0 2 ; region “3” (above dot-dashed line) is the region of stable rolls for q 2 < q o s c 2 , oscillatory unstable for q 2 slightly above q o s c 2 , and oscillatory or monotonically unstable for any q o s c 2 < q 2 < q 0 2 . Dashed and solid lines: blue color for χ 2 = 2 , brown color for χ 2 = 5 , green color for χ 2 = 10 .
Figure 8. Stability map of the stationary rolls in the case without surfactant. The solid black line is the boundary line between supercritical and subcritical rolls (“subcr. R”). Region “1” (below the dashed line) is the region of unstable rolls; in region “2”, rolls are stable within q 2 < q m 2 and monotonically unstable for q m 2 < q 2 < q 0 2 ; region “3” (above dot-dashed line) is the region of stable rolls for q 2 < q o s c 2 , oscillatory unstable for q 2 slightly above q o s c 2 , and oscillatory or monotonically unstable for any q o s c 2 < q 2 < q 0 2 . Dashed and solid lines: blue color for χ 2 = 2 , brown color for χ 2 = 5 , green color for χ 2 = 10 .
Colloids 06 00053 g008
Figure 9. Stability maps for N = 10 8 (a) and N = 10 7 (b). Other parameters are L = 0.003 and S = 1 . Regions “1”, “2”, and “3” are the same as Figure 3. Dashed and solid lines: blue color for χ 2 = 2 , brown color for χ 2 = 5 , green color for χ 2 = 10 .
Figure 9. Stability maps for N = 10 8 (a) and N = 10 7 (b). Other parameters are L = 0.003 and S = 1 . Regions “1”, “2”, and “3” are the same as Figure 3. Dashed and solid lines: blue color for χ 2 = 2 , brown color for χ 2 = 5 , green color for χ 2 = 10 .
Colloids 06 00053 g009
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Mikishev, A.B.; Nepomnyashchy, A.A. Marangoni Patterns in a Non-Isothermal Liquid with Deformable Interface Covered by Insoluble Surfactant. Colloids Interfaces 2022, 6, 53. https://doi.org/10.3390/colloids6040053

AMA Style

Mikishev AB, Nepomnyashchy AA. Marangoni Patterns in a Non-Isothermal Liquid with Deformable Interface Covered by Insoluble Surfactant. Colloids and Interfaces. 2022; 6(4):53. https://doi.org/10.3390/colloids6040053

Chicago/Turabian Style

Mikishev, Alexander B., and Alexander A. Nepomnyashchy. 2022. "Marangoni Patterns in a Non-Isothermal Liquid with Deformable Interface Covered by Insoluble Surfactant" Colloids and Interfaces 6, no. 4: 53. https://doi.org/10.3390/colloids6040053

APA Style

Mikishev, A. B., & Nepomnyashchy, A. A. (2022). Marangoni Patterns in a Non-Isothermal Liquid with Deformable Interface Covered by Insoluble Surfactant. Colloids and Interfaces, 6(4), 53. https://doi.org/10.3390/colloids6040053

Article Metrics

Back to TopTop