Next Article in Journal
The Application of Microfibrous Entrapped Activated Carbon Composite Material for the Sarin Simulant Dimethyl Methylphosphonate Adsorption
Previous Article in Journal
Zinc Oxide-Encapsulated Copper Nanowires for Stable Transparent Conductors
Previous Article in Special Issue
Exact Solutions for Non-Isothermal Flows of Second Grade Fluid between Parallel Plates
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Instability of Liquid Film with Odd Viscosity over a Non-Uniformly Heated and Corrugated Substrate

School of Mathematical Science, Inner Mongolia University, Hohhot 010021, China
*
Authors to whom correspondence should be addressed.
Nanomaterials 2023, 13(19), 2660; https://doi.org/10.3390/nano13192660
Submission received: 17 August 2023 / Revised: 19 September 2023 / Accepted: 26 September 2023 / Published: 28 September 2023
(This article belongs to the Special Issue Advances of Nanoscale Fluid Mechanics)

Abstract

:
The effect of odd viscosity on the instability of liquid film along a wavy inclined bottom with linear temperature variation is investigated. By utilizing the long-wave approximation, the non-linear evolution equation of the free surface is derived. By applying the normal mode method, the linear instability of thin film flow is investigated. With the help of multi-scale analysis methods, the weakly non-linear instability of thin film flow is also investigated. The results reveal that the Marangoni effect caused by non-uniform temperature distribution promotes the instability of the liquid film, while the odd viscosity has a stabilizing effect. In addition, for a positive local inclination angle θ , an increase in bottom steepness ζ inhibits the instability of the liquid film flow. In contrast, with a negative local inclination angle θ , increased bottom steepness ζ promotes the instability of the liquid film flow. The results of the temporal linear instability analysis and the weakly non-linear instability analysis have been substantiated through numerical simulations of the non-linear evolution equations.

1. Introduction

Liquid membranes find extensive applications in various industrial engineering fields, including thin film evaporators and liquid film dust collectors used in chemical equipment, as well as thin liquid film cooling of large-scale integrated circuits [1,2]. Additionally, the instability analysis of liquid films flowing along inclined or vertical planes is of significant relevance in the coating industry, where it plays a crucial role in microchip manufacturing, paper coating, magnetic film coating, and other related processes [3,4,5,6]. The complexities of liquid film flow have therefore garnered significant attention from researchers and practitioners alike, making it an important and constantly studied area of interest.
Benjamin [7] and Yih [8] were among the pioneers in studying the stability of liquid film flow over inclined planes. They focused on solving the linear stability problem for the fundamental flow of constant thickness and determined the critical Reynolds number R e c . When the Reynolds number exceeds R e c , the flow of liquid film becomes unstable. Samanta [9,10] conducted an analysis of the linear stability of thin liquid film flow on a non-uniformly heated inclined plane, neglecting the effect of inertial forces and employing the canonical mode method. Simultaneously, the long wave perturbation method was utilized to solve the control equation, leading to the derivation of a nonlinear surface wave equation. This enabled the analysis of linear, nonlinear, and sideband stability of liquid film flow on a vertically non-uniformly heated substrate. For additional studies on the stability of a thin film flow over the beveled and vertical planes, reference can be made to articles by Bauer and Kerczek [11], Hanratty [12] and Craster and Matar [13], among others. As for the stability of the film flow on a heated or inclined plane, research can be found in the work of Kalliadasis et al. [14], Sadiq et al. [15] and Mukhopadhyay et al. [16]. Gjevik [17] and Nakaya [18] extended the stability investigation of thin films to consider nonlinear stability.
Nevertheless, the aforementioned studies have primarily focused on liquid film flows over inclined or vertical substrates. In industrial equipment or other practical applications, certain substrates unavoidably have uneven surfaces. Solely investigating plane slopes may introduce bias when applied to non-level cases. In other words, the substrate’s structural characteristics significantly influence the stability of liquid film flow. The fluid dynamics of falling liquid films over undulating surfaces has attracted considerable attention, with early work by Pozrikidis [19] examining free surface Stokes flow along a sinusoidal base. Bielarz [20] discussed the stability of thin free surface liquid films flowing over structural surfaces, employing both one-dimensional and two-dimensional calculations. Wierschem and Aksel [21] performed a linear stability analysis of a Newtonian liquid film flowing down an inclined wavy plane. Their study investigated how wavy bottom variations, significantly longer than the film thickness, impact the stability of steady film flow compared to that over a flat inclined plane. Trifonov [22] studied the flow of a viscous liquid film along an inclined corrugated surface, using an integral model and employing Floquet theory to analyze the stability of nonlinear steady-state flows under various conditions. For research on the influence of substrate structure on film flow, references can be made to the works of Heining and Aksel [23], as well as Tougou [24]. These studies shed light on the crucial role played by the substrate’s unevenness or corrugations in determining the stability characteristics of a liquid film flow.
In recent years, odd viscosity has emerged as a prominent subject of investigation in the study of thin film flow stability. Odd viscosity [25,26] is the non-dissipative component of the viscosity tensor and is contained in its antisymmetric part. Avron [27,28] made a break-through discovery by demonstrating that in a classical fluid, when time-reversal symmetries are broken, either spontaneously or due to an external magnetic field or rotation, the viscosity tensor can have a non-zero odd part that gives rise to a dissipationless linear response coefficient known as odd or Hall viscosity. In natural situations where the time-reversal symmetry of a classical liquid is broken, odd viscosity effects are commonly observed in biological [29], granular [30] and colloidal [31] systems. Kirkinis and Andreev [32] explored the impact of odd viscosity on the thermocapillary instability of a viscous liquid film flowing along a uniformly heated solid substrate, while considering a fixed temperature gradient across the free surface. Their findings revealed that the initial wave of odd viscosity can suppress thermocapillary instability, leading to enhanced stability of the thin liquid film. A comprehensive exploration of odd viscosity in fluid dynamics was presented by Lapa et al. [33], providing a broader understanding of this intriguing phenomenon. Additionally, Zhao and Jian [34] investigated the effect of odd viscosity on the stability of a falling thin film in the presence of an electromagnetic field. Employing the lubrication approximation, they derived a new nonlinear evolution equation for the free surface that takes into account the influence of odd viscosity. Through linear and weakly nonlinear stability analyses, they found that odd viscosity has a stabilizing effect on the system. These studies collectively contribute to our understanding of the significant role that odd viscosity plays in influencing the stability characteristics of thin film flows, and highlight its relevance in various physical systems and applications.
This paper is centered around investigating the influence of odd viscosity on the instability of liquid film flowing down an undulated inclined plate with linear temperature variation. Specifically, the study examines the effects of odd viscosity, thermocapillary effects and bottom steepness on flow instability, leading to corresponding conclusions.

2. Mathematical Model

We consider the flow of an incompressible viscous Newtonian fluid driven by gravity in two dimensions. The fluid flows down an inclined, corrugated substrate with uneven heating, as depicted in Figure 1. The Cartesian coordinate system e x ^ , e y ^ has an angle β with respect to the horizontal, and the base contour b ^ ( x ^ ) is periodic with amplitude a ^ and wavelength λ ^ , where x ^ is in the direction of the main flow. The substrate profile is undulated, so a local curvilinear coordinate system is introduced, and the when thickness of the thin films are thin enough compared to the curvature radius of the bottom, the flow ( u , v ) is still mainly parallel to the bottom [35]. Therefore, at every point of the bottom, a local coordinate system e x , e y is defined, where e x is the tangent and e y is normal to the base. For any point P in the fluid, the coordinates of the curve are the arc length x of the base and the distance y along the e y to the base. In e x ^ , e y ^ coordinates P = ( x ^ sin θ y , b ^ ( x ^ ) + cos θ y ) , where θ = θ ( x ^ ) = arctan ( b ^ ( x ^ ) / x ^ ) is the local inclination angle between e x ^ and e x .
This is the considered film flows along a moderate steepness undulating substrate and has a long length compared to the film thickness. The substrate curvature κ is denoted by
κ ( x ^ ) = 2 b ^ ( x ^ ) x ^ 2 1 + b ^ ( x ^ ) x ^ 2 3 / 2 .
For further details on the transformation to curvilinear coordinates, we refer to the study of Wierschem et al. [36].
We assume that such fluids are time reversal symmetry breaking, and the falling film problem should consider both even τ e and odd τ o viscosity coefficients. In this case, the Cauchy stress tensor τ is expressed as
τ = τ e + τ o ,
τ i j e = p δ i j + η e u i x j + u j x i ,
τ i j o = η o δ i 1 δ j 1 δ i 2 δ j 2 u 1 x 2 + u 2 x 1 + η o δ i 1 δ j 2 δ i 2 δ j 1 u 1 x 1 u 2 x 2 ,
where i, j = 1, 2, η e and η o denote the odd and even viscosity coefficients, respectively.
The dynamic properties of an incompressible fluid can be described by the continuity equation and the momentum equation. Referring to Wierschem et al. [36], the following equations are derived based on the curvilinear coordinate transformation considering the odd viscosity and the Marangoni effect. And, to simplify the notations, we set u = u 1 , v = u 2 ,   x = x 1 , and y = x 2 .
1 1 + κ y u x + κ v + v y = 0 ,
ρ u t + 1 1 + κ y u u x + κ v + v u y = 1 1 + κ y p x + ρ g sin ( β θ ) + η e 1 1 + κ y 3 κ x v y u x + 1 1 + κ y 2 2 u x 2 κ 2 u + 2 κ v x + 1 1 + κ y κ u y + 2 u y 2 η o 1 1 + κ y 3 κ x u + y v x + 1 1 + κ y 2 2 v x 2 κ 2 v 2 κ u x + 1 1 + κ y κ v y + 2 v y 2 ,
ρ v t + 1 1 + κ y u v x κ u + v v y = p y ρ g cos β θ + η e 1 1 + κ y 3 κ x u + y v x + 1 1 + κ y 2 2 v x 2 κ 2 v 2 κ u x + 1 1 + κ y κ v y + 2 v y 2 + η o 1 1 + κ y 3 κ x v y u x + 1 1 + κ y 2 2 u x 2 κ 2 u + 2 κ v x + 1 1 + κ y κ u y + 2 u y 2 ,
T t + 1 1 + κ y u T x + v T y = k c 1 1 + κ y 2 2 T x 2 1 1 + κ y 3 y κ x T x + 1 1 + κ y κ T y + 2 T y 2 ,
where ρ is the liquid density, p is the pressure, T is the temperature, g is the gravitational acceleration, and k c is the thermal diffusivity, which is assumed to be constant.
At the substrate y = 0 , the boundary conditions for fluid no-slip and no-penetration and the boundary condition for temperature are
u = v = 0 ,
T = T g + b x ,
where T g is the ambient temperature, b = Δ T / λ ^ is the linear rate of change of temperature, Δ T = T H T C , where T H and T C denote the temperatures at the hotter part and colder part of the substrate, respectively.
At the free surface y = h ( x , t ) , dynamic and kinematic boundary conditions are
1 1 + 1 1 + κ h h x 2 [ η e { ( 1 + κ h ) 2 h x 2 } 1 1 + κ h v x κ u + u y + 2 η e ( 1 + κ h ) v y u x + κ v h x + η o { ( 1 + κ h ) 2 h x 2 } 1 1 + κ h u x + κ v v y + 2 η o ( 1 + κ h ) u y + v x κ u h x ] = ( 1 + κ h ) σ x + h x σ y ,
p a p + 1 1 + 1 1 + κ h h x 2 [ 2 η e 1 1 + κ h 3 u x + κ v h x 2 + v y 1 1 + κ h 2 v x κ u × h x 1 1 + κ h u y h x + η o 1 1 1 + κ h h x 2 u y + 1 1 + κ h v x κ u 2 1 + κ h × 1 1 + κ h u x + κ v v y h x ] = σ 1 + κ h h x x h κ x h x κ 1 + κ h 2 + 2 h x 2 1 + κ h 2 + h x 2 3 / 2 ,
h t + 1 1 + κ h u h x v = 0 ,
where σ is the surface tension, and is assumed to vary linearly over a small temperature range
σ = σ 0 γ ( T T g ) ,
σ 0 is the surface tension at the reference temperature T g and γ = σ / T T = T g .
The balance between heat supply to and heat loss at the free surface y = h ( x , t ) is given by Newton’s law of cooling:
λ 1 + 1 1 + κ h h x 2 1 1 + κ h 2 h x T x + T y = k g T T g ,
where λ is the thermal conductivity, and k g is the heat transfer coefficient between the fluid and the air.
In order to investigate the effect of substrate undulation on the film flow, we use the thin film flow over a flat bottom as referenced. So the Nusselt velocity u 0 = ρ g h 0 2 sin θ / 3 η e , where h 0 is the constant film thickness, is also the length scale in the transverse direction, and λ ^ is the characteristic longitudinal length scale, which is very long compared to the film thickness.
The following dimensionless quantities are introduced to dimensionlessize the equation (indicated by the asterisk):
x = 2 π x λ ^ , y = y h 0 , h = h h 0 , u = u u 0 , v = λ ^ v 2 π h 0 u 0 , t = 2 π u 0 t λ ^ , p = p ρ u 0 2 , κ = λ ^ 2 κ 4 π 2 a ^ , T = 2 π T T g Δ T , x ^ x = 2 π x ^ x λ ^ , b ^ x ^ = b ^ a ^ λ ^ x ^ 2 π , θ = arctan ζ b ^ x ^ ,
where ζ = 2 π a ^ / λ ^ is the bottom steepness and α = 2 π h ^ / λ ^ is the aspect ratio.
Using the Equation (16) in the governing equations and boundary conditions, we arrive after dropping the asterisk as
u x + v y + α ζ κ v + y v y = 0 ,
α R e u t + u u x + v u y = α R e p x + 3 sin ( β θ ) sin β + 2 u y 2 + α κ ζ u y α μ 2 v y 2 + O α 2 ,
α R e ζ κ u 2 = R e p y 3 cos β θ sin β + α 2 v y 2 + μ α ζ κ u y + 2 u y 2 + O α 2 ,
α R e P r T t + u T x + v T y = 2 T y 2 + α ζ κ T y + O α 2 .
At the substrate y = 0 , we have
u = 0 , v = 0 , T = x .
At the free surface y = h ( x , t ) , we have
u y + α ζ κ 2 h u y u + α μ u x v y + 2 u y h x = M n T x + h x T y + O α 2 ,
p ¯ a p + 2 α R e v y u y h x + μ R e u y = α 2 W e ( 1 C a T ) h x x ξ κ + 2 ζ 2 κ 2 h + O α 2 ,
h t + u h x v α ζ κ h u h x + O α 2 = 0 ,
T y + B i T + O α 2 = 0 ,
where p ¯ a = p a ρ u 0 2 , μ = η o η e is the odd viscosity coefficient, R e = ρ u h η e is the Reynolds number, W e = σ 0 η e u 0 is the Weber number, M n = α γ Δ T η e u 0 is the Marangoni number, P r = η e ρ k c is the Prandtl number, C a = γ Δ T σ 0 is the Capillary number, and B i = k g h 0 λ is the Biot number, ξ = ζ α a ^ h ^ .

3. Approximate Solution of the Equations

The physical quantities u, v, p and T are expanded as power series of the small parameter α :
u = u 0 + α u 1 + , v = v 0 + α v 1 + , p = p 0 + α p 1 + , T = T 0 + α T 1 + .
Then we substitute the asymptotic Equation (26) into the dimensionless Equations (17)–(25) to obtain the zero-order governing equations and boundary conditions
u 0 x + v 0 y = 0 ,
3 sin β θ sin β + 2 u 0 y 2 = 0 ,
R e p y = 3 cos β θ sin β + μ 2 u 0 y 2 ,
2 T 0 y 2 = 0 .
At the substrate y = 0 ,
u 0 = v 0 = 0 , T 0 = x .
At the free surface y = h ( x , t ) ,
u 0 y = M n T 0 x + h x T 0 y ,
p ¯ a p 0 + μ R e u 0 y = α 2 W e h x x ξ κ + 2 ζ 2 κ 2 h ,
h t + u 0 h x v 0 = 0 ,
T 0 y = 0 .
The solutions at zeroth order can be found as
u 0 = 3 sin β θ sin β h y 1 2 y 2 M n y ,
v 0 = 3 sin β θ 2 sin β h x y 2 ,
p = 1 R e 3 cos β θ sin β + μ 3 sin β θ sin β h y + p ¯ a μ R e α 2 W e h x x ξ κ + 2 ζ 2 κ 2 h .
The first order governing equations and boundary conditions are obtained.
u 1 x + v 1 y + ζ κ v 0 + ζ κ y v 0 y = 0 ,
R e u 0 t + u 0 u 0 x + v 0 u 0 y = R e p 0 x + κ ζ u 0 y μ 2 v 0 y 2 + 2 u 1 y 2 ,
R e ζ κ u 0 2 = R e p 1 y + 2 v 0 y 2 + μ ζ κ u 0 y + 2 u 1 y 2 ,
R e P r T 0 t + u 0 T 0 x + v 0 T 0 y = 2 T 1 y 2 + ζ κ T 0 y .
At the substrate y = 0 ,
u 1 = v 1 = 0 , T 1 = 0 .
At the free surface y = h ( x , t ) ,
u 1 y + ζ κ 2 h u 0 y u 0 + μ u 0 x v 0 y + 2 u 0 y h x = M n T 1 x + h x T 1 y ,
p 1 + 2 R e v 0 y u 0 y h x + μ R e u 1 y = 0 ,
u 1 h t v 1 ζ κ h u 0 h x = 0 ,
T 1 y = 0 .
By solving the first order equation, we obtain the expression of T 1 and u 1 as
T 1 = R e P r { sin β θ sin β 1 8 y 4 + 1 2 h y 3 h 3 y M n 1 6 y 3 1 2 h 2 y } ,
u 1 = 1 2 R e sin β θ sin β y 3 3 h 2 y h t + M n R e P r 5 2 sin β θ sin β h 3 M n h 2 h x y + R e × sin β θ sin β 3 8 sin β θ sin β y 3 4 h 3 h y 1 8 M n y 3 4 h 3 y h x + R e 1 2 y 2 h y 3 cos β θ R e sin β h x α 2 W e h x x x ξ κ x + 2 ζ 2 κ 2 h x + 4 ζ 2 κ h κ x + ζ κ × sin β θ sin β 1 2 y 3 3 2 h y 2 + 3 h 2 y + 1 2 M n y 2 2 μ 3 sin β θ sin β h M n h x y .
From the continuity equation, we have
h t = 3 sin β θ sin β h 2 M n h h x .
Putting this in (49), we obtain
u 1 = M n R e P r 5 2 sin β θ sin β h 3 M n h 2 h x y + R e sin β θ sin β { sin β θ sin β 3 8 h y 4 3 2 h 2 y 3 + 3 h 4 y + M n 1 8 y 4 + 1 2 h y 3 h 3 y } h x + R e 1 2 y 2 h y { 3 cos ( β θ ) R e sin β h x α 2 W e h x x x ξ κ x + 2 ζ 2 κ 2 h x + 4 ζ 2 κ h κ x } + ζ κ { sin β θ sin β 1 2 y 3 3 2 h y 2 + 3 h 2 y + 1 2 M n y 2 } 2 μ 3 sin β θ sin β h M n h x y .
The local flow rate q ( x , t ) is defined in
q ( x , t ) = 0 h u ( x , y , t ) d y ,
where u ( x , y , t ) = u 0 ( x , y , t ) + α u 1 ( x , y , t ) + O ( α 2 ) . Solving for the integral on the right-hand side of Equation (50), we obtain
q ( x , t ) = sin β θ sin β h 3 1 2 M n h 2 + α [ M n R e 5 4 sin β θ sin β h 5 1 2 M n h 4 h x + R e sin β θ sin β × 6 5 sin β θ sin β h 6 2 5 M n h 5 h x R e 3 { 3 cos β θ R e sin β h x α 2 W e ( h x x x ξ κ x + 2 ζ 2 κ 2 h x + 4 ζ 2 κ h κ x ) } h 3 + ζ κ 9 8 sin ( β θ ) sin β h 4 + 1 6 M n h 3 μ 3 sin ( β θ ) sin β h 3 M n h 2 h x ] .
Using alternative form of the kinematic boundary conditions
h t + ( 1 α ζ κ h ) q x + O α 2 = 0 .
We obtain a non-linear evolution equation for the thickness of the film
h t + A ( h ) h x + α ( B ( h ) h x + C ( h ) h x x x ) x = 0 ,
where,
A ( h ) = 3 sin β θ sin β h M n h + α { α 2 W e R e ξ + 16 3 ζ 2 κ h h 2 κ x + 3 2 ζ κ h 2 sin ( β θ ) sin β h + M n } ,
B ( h ) = M n R e P r 5 4 sin β θ sin β h 5 1 2 M n h 4 + R e sin β θ sin β 6 5 sin β θ sin β h 6 2 5 M n h 5 cos β θ sin β 2 3 α 2 W e R e ζ 2 κ 2 h 3 μ 3 sin β θ sin β h 3 M n h 2 ,
C ( h ) = 1 3 α 2 W e R e h 3 .
We introduce a new parameter S, as defined by earlier researchers [37]:
S = α 2 W e .
From now on, for the sake of simplicity of the formula, we shall use the following abbreviation:
s = sin ( β θ ) sin β , c = cos ( β θ ) sin β .

4. Linear Stability Analysis

To investigate the stability of the thin film flow, a small perturbation at the free interface is assumed. The film thickness h, which can be written as
h = 1 + η ( x , t ) ,
where η 1 denotes the dimensionless distance of free surface of the liquid film from free surface of the smooth laminar flow.
To eliminate the α in Equation (55), we set the conversion
x = α x ˜ , t = α t ˜ .
Substituting Equations (61) and (62) into Equation (55), retaining up to O ( η 3 ) , after dropping the cap sign, we obtain
η t + A 1 η x + B 1 η x x + C 1 η x x x x + A 1 η η x + B 1 η η x x + η x 2 + C 1 η η x x x x + η x η x x x + 1 2 A 1 η 2 η x + B 1 1 2 η 2 η x x + η η x 2 + C 1 1 2 η 2 η x x x x + η η x η x x x + O η 4 = 0 ,
where A 1 , B 1 , C 1 and their corresponding derivatives (denoted by primes) are the values corresponding to h = 1 . The linear response of the film is studied by assuming that the perturbation is in the form of a sinusoidal perturbation, that is
η ( x , t ) = Γ exp [ i ( k x ω t ) ] + c . c . ,
where Γ is the amplitude of the disturbance, k is the wave number, c . c . is the complex conjugate and ω = ω r + i ω i is the complex frequency, ω r and ω i are the linear growth rates of oscillation frequency and amplitude, respectively. Substituting Equation (64) into Equation (63) and considering the linear part, the dispersion relation can be obtained as
Disp ( ω , k ) = i ω + i A 1 k B 1 k 2 + C 1 k 4 = 0 .
In Equation (65), the real and the imaginary parts of ω are expressed as
ω r = A 1 k , ω i = B 1 k 2 C 1 k 4 .
The flow will be linearly unstable if the linear growth rate of amplitude is ω i > 0 , conversely, the flow is linearly stable, and the flow will be neutrally stable if ω i = 0 . At this point, the critical Reynolds number is
R e c = c 2 B o ζ 2 κ 2 + μ 3 s M n 1 4 M n P r 5 s 2 M n + 2 5 s 3 s M n ,
where B o = 4 π 2 σ 0 / ρ g λ ^ 2 sin β is the inverse Bond number and 3 B o = S R e .
As μ 0 , we obtained the same critical Reynolds number as that derived by Mukhopadhyay and Mukhopadhyay [35].

5. Weakly Non-Linear Analysis

We use the method of multiple scale and expend the interfacial perturbation η in the following form [38,39,40]
η ( x , x 1 , , t , t 1 , t 2 , ) = ε η 11 + ε 2 η 12 + ε 3 η 13 + ,
where
x 1 = ε x , t 1 = ε t , t 2 = ε 2 t , ,
here x and t represent the rapidly varying scales, while x 1 , t 1 , and so on represent the slowly varying scales. Assuming these variables are independent of each other, then the derivatives of time and space become as shown in
t t + ε t 1 + ε 2 t 2 ,
x x + ε x 1 .
Substituting the Equations (68)–(71) into the Equation (63), we obtain
( L 0 + ε L 1 + ε 2 L 2 + ) ( ε η 11 + ε 2 η 12 + ε 3 η 13 + ) = ε 2 N 2 ε 3 N 3 ,
where the operators L 0 , L 1 , L 2 , and the non-linear terms N 2 , N 3 in Equation (72) are shown below.
L 0 = t + A 1 x + B 1 2 x 2 + C 1 4 x 4 , L 1 = t 1 + A 1 x 1 + 2 B 1 2 x x 1 + 4 C 1 4 x 3 x 1 , L 2 = t 2 + B 1 2 x 1 2 + 6 C 1 4 x 2 x 1 2 , N 2 = A η 11 η 11 x + B 1 η 11 2 η 11 x 2 + η 11 x 2 + C 1 η 11 4 η 11 x 4 + η 11 x 3 η 11 x 3 , N 3 = A 1 η 11 η 12 x + η 11 x 1 + η 12 η 11 x + B 1 [ η 11 2 η 12 x 2 + 2 2 η 11 x x 1 + η 12 2 η 11 x 2 + 2 η 11 x η 12 x + η 11 x 1 ] + C 1 [ η 11 4 η 12 x 4 + 4 4 η 11 x 3 x 1 + η 12 4 η 11 x 4 + η 11 x × 3 η 12 x 3 + 3 3 η 11 x 2 x 1 + 3 η 11 x 3 η 11 x 1 + η 12 x ] + 1 2 A 1 η 11 2 η 11 x + B 1 1 2 η 11 2 2 η 11 x 2 + η 11 η 11 x 2 .
For the first-order equation of ε , we have
L 0 η 11 = 0 ,
the solution of this equation has the following form:
η 11 = Γ ( x 1 , t 1 , t 2 ) [ exp i Θ ] + c . c . ,
where Γ ( x 1 , t 1 , t 2 ) is the complex amplitude. Θ = k x ω r t , ω = ω r + i ω i is the complex frequency, because near the neutral curve ω i O ( ζ 2 ) . Thus, the function exp ( ω i t ) is slow and may be absorbed in Γ ( x 1 , t 1 , t 2 ) .
For the second-order equation of ε , we have
L 0 η 12 = L 1 η 11 N 2 .
Substituting the expression (75) into the Equation (76), we obtain
L 0 η 12 = i Disp ( ω r , k ) ω r Γ t 1 Disp ( ω r , k ) k Γ x 1 e i Θ Q 1 Γ 2 e 2 i Θ + c . c . ,
where Disp ( ω r , k ) is given by Equation (65) and
Q 1 = i A 1 k 2 B 1 k 2 + 2 C 1 k 4 .
The solution η 12 is obtained from Equation (77) in the following form:
η 12 = Q 1 Γ 2 e 2 i Θ Disp ( 2 ω r , 2 k ) + c . c . .
Introducing the coordinate transformation δ = x 1 c g t 1 , where c g = ( Disp k / Disp ω r ) , and using the solvability condition of the third-order equation, we obtain
Γ t 2 + J 1 2 Γ δ 2 ε 2 ω i Γ + J 2 + i J 3 | Γ | 2 Γ = 0 ,
where
J 1 = B 1 6 C 1 k 2 , J 2 = 1 2 C 1 k 4 B 1 k 2 + A 1 2 k 2 + 2 7 C 1 k 4 B 1 k 2 B 1 k 2 C 1 k 4 16 C 1 k 4 4 B 1 k 2 , J 3 = 1 2 A 1 k + A 1 k B 1 k 2 7 C 1 k 4 + 2 A 1 k B 1 k 2 C 1 k 4 16 C 1 k 4 4 B 1 k 2 .
The diffusion effect in Equation (80) is neglected to obtain
Γ t 2 ε 2 ω i Γ + J 2 + i J 3 | Γ | 2 Γ = 0 .
The solution of the equation can be written as
Γ = a exp [ i b ( t 2 ) t 2 ] .
Substituting Equation (83) into Equation (82), we obtain
a t 2 = ε 2 ω i J 2 a 2 a ,
[ b ( t 2 ) t 2 ] t 2 = J 3 a 2 .
The second term on the right-hand side of Equation (84) is a non-linear term that regulates the exponential change of the linear perturbation, depending on the sign of ω i and J 2 . When the right side of Equation (84) is 0, the equilibrium amplitude is solved as
ε a = ω i J 2 .
The term J 2 < 0 can lead to the instability of the system. According to the signs of ω i and J 2 , four non-linear regions are defined, which are: the supercritical stability region I ( ω i > 0 , J 2 > 0 ) , the subcritical instability region I I ( ω i < 0 , J 2 < 0 ) , the unconditionally stable region I I I ( ω i < 0 , J 2 > 0 ) and the explosive state region I V ( ω i > 0 , J 2 < 0 ) .

6. Numerical Simulations

We numerically solve Equation (55) within a periodic domain to comprehend the evolution of finite amplitude perturbations. To achieve this, we transform the problem into a set of ordinary differential equations (ODE) by discretizing the spatial variables into a series of wave numbers through the Fourier transform. Subsequently, we employ the fast Fourier transform algorithm to compute these wave numbers as the numerical solution to the nonlinear evolution Equation (55) [41].
Fourier transform of the Equation (55) in the x domain
h ^ t = i k · F { s F 1 [ h ^ ] 3 1 2 M n F 1 [ h ^ ] 2 + ζ κ 3 8 s F 1 [ h ^ ] 3 S R e ( 1 3 ξ F 1 [ h ^ ] 3 4 3 ζ 2 κ × F 1 [ h ^ ] 4 ) κ x + M n R e P r 5 4 s F 1 [ h ^ ] 5 1 2 M n F 1 [ h ^ ] 4 F 1 [ i k h ^ ] + R e s 6 5 s × F 1 [ h ^ ] 6 2 5 M n F 1 [ h ^ ] 5 F 1 [ i k h ^ ] c F 1 [ h ^ ] 3 F 1 [ i k h ^ ] + 2 3 S R e ζ 2 κ 2 F 1 [ h ^ ] 3 × F 1 [ i k h ^ ] μ 3 s F 1 [ h ^ ] 3 M n F 1 [ h ^ ] 2 F 1 [ i k h ^ ] + 1 3 S R e F 1 [ h ^ ] 3 F 1 [ ( i k ) 3 h ^ ] } .
Initially, a finite-amplitude monochromatic disturbance is given as
h ( x , 0 ) = 1 + 0.03 cos k x ,
where k is the wave number. The computation is performed on a uniform grid with the number of spatial grid points varying N = 200 and Δ t = 0.1∼0.2.

7. Specific Case Study

We choose the substrate profile as
b ^ ( x ^ ) = a ^ sin 2 x ^ / λ ^ ,
where λ ^ = 3 is the wavelength and a ^ = 0.15 is the amplitude of the wavy bottom profile [35]. x ^ = 10 π / 3 is a point on the “uphill” and x ^ = 5 π / 3 is a point on the “downhill”. We will understand more easily according to Figure 1. Moreover, the basal steepness ζ is a fixed quantity rather than a perturbation parameter. Since both the basal curvature κ ( x ^ ) and the local inclination θ ( x ^ ) are functions of x ^ , the critical Reynolds number is also a function of x ^ .

8. Results and Discussion

8.1. Linear Stability Analysis

Figure 2 depicts the temporal growth rate curve for various values of the Marangoni number. It is evident that the Marangoni number exerts a destabilizing effect on the system. The destabilizing role of the Marangoni number M n can be explained from a physical perspective as follows: when the liquid film flows downward along the substrate, the thermocapillary force appearing at the interface acts in the opposite direction of the gravitational acceleration as M n increases and then enhances the growth of the surface instabilities [42].
Figure 3 displays the temporal growth rate curve for different values of odd viscosity. Notably, the odd viscosity μ has a stabilizing effect on the system. This is because the odd viscosity μ comes into play effectively through shear stress on the free surface as the liquid film flows downward along the substrate, causing additional stress on the free surface, which increases the critical Reynolds number giving the stabilizing effect [32]. Equation (67) further supports this observation, showing that the critical Reynolds number R e c rises as the odd viscosity μ increases. As a result, higher values of odd viscosity lead to a more stable fluid flow, reducing the likelihood of surface perturbations and instability. In addition, we compare Figure 2a,b, or Figure 3a,b: interestingly, the fluid exhibits greater stability in the “uphill” part under the same conditions. This distinction is attributed to the difference in the sign of the local inclination angle θ in these regions. The variation in θ plays a significant role in influencing the stability characteristics of the fluid flow along the substrate.
In Figure 4, it is evident that under same conditions, the critical Reynolds number of the “uphill” point increases with an increase in the bottom steepness ζ . In contrast, the critical Reynolds number of the “downhill” point decreases as the bottom steepness ζ increases. From Equation (67), for a positive local inclination angle θ , an increase in bottom steepness ζ inhibits the instability of the liquid film flow. Conversely, with a negative local inclination angle θ , increased bottom steepness ζ promotes the instability of the liquid film flow.
Moreover, as M n 0 , μ 0 and ζ 0 , we find from Figure 5 that the value of the critical Reynolds number is equal to R e c = ( 5 / 6 ) cot β 0.48 , which was originally obtained by Benjamin [7] and Yih [8].

8.2. Weakly Non-Linear Stability Analysis

From Figure 6, we observe that, under same conditions, the subcritical instability region II and unconditional stability region III gradually enlarge as the odd viscosity μ increases. In contrast, the supercritical stability region I and the explosive region IV gradually shrink. Moreover, the critical Reynolds number R e c increases with the odd viscosity μ . This implies that the odd viscosity has a stabilizing effect on the flow, as it leads to larger stability regions and higher critical Reynolds numbers.
Similarly, Figure 7 reveals that, under same conditions, the subcritical instability region II and unconditional stability region III gradually shrink as the Marangoni number M n increases. Conversely, the supercritical stability region I and the explosive region IV gradually enlarge. Additionally, the critical Reynolds number R e c decreases with the increase in the Marangoni number. Hence, the Marangoni effect destabilizes the flow, resulting in reduced stability regions and lower critical Reynolds numbers.

8.3. Numerical Simulations

Figure 8 and Figure 9 illustrate the short-time evolution of the free surface for different values of odd viscosity μ and the Marangoni number M n in the subcritical instability region, respectively. As shown in both figures, the disturbance amplitude gradually decreases over time, indicating that the liquid film stabilizes as the disturbances dampen. This behavior is characteristic of the subcritical instability region, where disturbances do not grow unbounded but rather stabilize over time. Furthermore, when comparing Figure 8 and Figure 9, we can observe that the odd viscosity μ has a stabilizing effect. Conversely, the Marangoni number M n has a destabilizing effect.
In Figure 10a, we can see that under the same conditions, the amplitude of the free surface disturbance first becomes smaller and then larger as the bottom steepness ζ increases in the subcritical instability region. This observation supports the previous conclusion that the critical Reynolds number of the “uphill” point increases with the bottom steepness ζ , while the “downhill” point decreases with increasing ζ . The varying amplitudes of the free surface disturbance further highlight the contrasting roles of the bottom steepness ζ in influencing the stability of the fluid flow at different regions of the substrate. We can observe wave changes within the explosive region, as depicted in Figure 10b and Figure 11. The obtained results align with those within the subcritical unstable region.

9. Conclusions

We have mainly studied the effect of odd viscosity on the instability of falling liquid film over a non-uniformly heated inclined corrugated substrate. To simplify the analysis, we neglect the evaporation effect by assuming a nonvolatile fluid. Additionally, for the sake of convenience, we consider the free surface to be adiabatic. The mechanism of thermocapillary helps move the fluid from a warmer region to a colder neighborhood. This study is conducted on general periodic bottom contours and is analyzed and discussed in the specific case of sinusoidal bottoms.
In our investigation, we carefully considered the impact of various factors, including odd viscosity, thermocapillary effects, and bottom steepness. By taking these factors into account, we derived the non-linear evolution equation. Temporal linear stability analysis are performed based on Equation (55), and we find that the odd viscosity has a stabilizing effect, while the Marangoni number always has a destabilizing effect. We observed that regardless of the presence of odd viscosity or the value of the Marangoni number, the fluid flow is consistently more stable at the “uphill” point compared to the “downhill” point. This distinction is a consequence of different signs of the local inclination angle.
In the investigation of weakly non-linear stability, we employed the multiple scales method to derive the Ginzburg–Landau equation. By analyzing the signs of ω i and J 2 , we have identified four distinct nonlinear regions. We find that the subcritical instability region II and the unconditional stability region III enlarge with odd viscosity μ increases, and conversely, the supercritical stability region I and the explosive region IV shrink. While the results for the Marangoni effect are opposite to the odd viscosity effect.
To gain a deeper understanding of the stability of liquid film flow and investigate the influence of different parameters, we employ the fast Fourier transform method to solve the non-linear Equation (55). The results obtained from the numerical simulations are in agreement with the findings of the previous linear stability analysis. Numerical simulations indicate that, increasing odd viscosity diminishes perturbations and decreases the wave’s maximum height, which is precisely the opposite to the Marangoni effect. While with the increase in the bottom steepness, the height of the wave amplitude lowers first and then rises.

Author Contributions

Conceptualization, D.X. and Z.D.; methodology, Q.L.; software, R.Z.; validation, D.X. and Z.D.; formal analysis, Q.L. and R.Z.; writing—original draft preparation, D.X.; writing—review and editing, D.X.; supervision, Z.D. All authors have read and agreed to the published version of the manuscript.

Funding

The authors acknowledge financial support provided by the National Natural Science Foundation of China (Nos. 11902165, 12272188, 12102205 and 12262025), the National Science Foundation for Distinguished Young Scholars of the Inner Mongolia Autonomous Region of China (No. 2023JQ16), the Program for Young Talents of Science and Technology in Universities of Inner Mongolia Autonomous Region (NJYT23098), the Scientific Startin and the Innovative Research Team in Universities of Inner Mongolia Autonomous Region of China (No. NMGIRT2208).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Oron, A.; Davis, S.H.; Bankoff, S.G. Long-scale evolution of thin liquid films. Rev. Mod. Phys. 1997, 69, 931. [Google Scholar] [CrossRef]
  2. Tseluiko, D.; Papageorgiou, D.T. Wave evolution on electrified falling films. J. Fluid Mech. 2006, 556, 361–386. [Google Scholar] [CrossRef]
  3. Kalliadasis, S.; Kiyashko, A.; Demekhin, E.A. Marangoni instability of a thin liquid film heated from below by a local heat source. J. Fluid Mech. 2003, 475, 377–408. [Google Scholar] [CrossRef]
  4. Chen, C.I.; Chen, C.K.; Yang, Y.T. Weakly nonlinear stability analysis of thin viscoelastic film flowing down on the outer surface of a rotating vertical cylinder. Int. J. Eng. Sci. 2003, 41, 1313–1336. [Google Scholar] [CrossRef]
  5. Thiele, U.; Goyeau, B.; Velarde, M.G. Stability analysis of thin film flow along a heated porous wall. Phys. Fluids 2009, 21, 014103. [Google Scholar] [CrossRef]
  6. Closa, F.; Ziebert, F.; Raphael, E. Effects of In-plane Elastic Stress and Normal External Stress on Viscoelastic Thin Film Stability. Math. Model. Nat. Phenom. 2012, 7, 6–19. [Google Scholar] [CrossRef]
  7. Benjamin, T.B. Wave formation in laminar flow down an inclined plane. J. Fluid Mech. 1957, 2, 554–573. [Google Scholar] [CrossRef]
  8. Yih, C.S. Stability of liquid flow down an inclined plane. Phys. Fluids 1963, 6, 321–334. [Google Scholar] [CrossRef]
  9. Samanta, A. Stability of liquid film falling down a vertical non-uniformly heated wall. Physica D 2008, 237, 2587–2598. [Google Scholar] [CrossRef]
  10. Samanta, A. Stability of inertialess liquid film flowing down a heated inclined plane. Phys. Lett. A 2008, 372, 6653–6657. [Google Scholar] [CrossRef]
  11. Bauer, R.J.; Kerczek, C.V. Stability of liquid film flow down an oscillating wall. J. Appl. Mech. 1991, 58, 278–282. [Google Scholar] [CrossRef]
  12. Hanratty, T.J. Interfacial instabilities caused by air flow over a thin liquid layer. In Waves on Fluid Interfaces; Meyer, R.E., Ed.; Academic Press: Cambridge, MA, USA, 1983; pp. 221–259. [Google Scholar]
  13. Craster, R.V.; Matar, O.K. Dynamics and stability of thin liquid films. Rev. Mod. Phys. 2009, 81, 1131. [Google Scholar] [CrossRef]
  14. Trevelyan, P.M.J.; Kalliadasis, S. Wave dynamics on a thin-liquid film falling down a heated wall. J. Eng. Math. 2004, 50, 177–208. [Google Scholar] [CrossRef]
  15. Sadiq, I.M.R.; Usha, R.; Joo, S.W. Instabilities in a liquid film flow over an inclined heated porous substrate. Chem. Eng. Sci. 2010, 65, 4443–4459. [Google Scholar] [CrossRef]
  16. Mukhopadhyay, A. Stability of a thin viscous fluid film flowing down a rotating non-uniformly heated inclined plane. Acta Mech. 2011, 216, 225–242. [Google Scholar] [CrossRef]
  17. Gjevik, B. Occurrence of Finite-Amplitude Surface Waves on Falling Liquid Films. Phys. Fluids 1970, 13, 1918–1925. [Google Scholar] [CrossRef]
  18. Nakaya, C. Long waves on a thin fluid layer flowing down an inclined planes. Phys. Fluids 1975, 18, 1407–1412. [Google Scholar] [CrossRef]
  19. Pozrikidis, C. The flow of a liquid film along a periodic wall. J. Fluid Mech. 1988, 188, 275–300. [Google Scholar] [CrossRef]
  20. Bielarz, C.; Kalliadasis, S. Time-dependent free-surface thin film flows over Topography. Phys. Fluids 2003, 15, 2512–2524. [Google Scholar] [CrossRef]
  21. Wierschem, A.; Aksel, N. Instability of a liquid film flowing down an inclined wavy plane. Physica D 2003, 186, 221–237. [Google Scholar] [CrossRef]
  22. Trifonov, Y.Y. Viscous liquid film flow down an inclined corrugated surface. Calculation of the flow stability to arbitrary perturbations using an integral method. J. Appl. Mech. Tech. Phys. 2016, 57, 195–201. [Google Scholar] [CrossRef]
  23. Heining, C.; Aksel, N. Bottom reconstruction in thin-film flow over topog-raphy: Steady solution and linear stability. Phys. Fluids 2009, 21, 083605. [Google Scholar] [CrossRef]
  24. Tougou, H. Long waves on a film flow of a viscous fluid down an inclined uneven wall. J. Phys. Soc. Jpn. 1978, 44, 1014–1019. [Google Scholar] [CrossRef]
  25. Fruchart, M.; Scheibner, C.; Vitelli, V. Odd Viscosity and Odd Elasticity. Annu. Rev. Condens. Matter Phys. 2023, 14, 471–510. [Google Scholar] [CrossRef]
  26. Granero-Belinchón, R.; Ortega, A. On the Motion of Gravity–Capillary Waves with Odd Viscosity. J. Nonlinear Sci. 2022, 32, 28. [Google Scholar] [CrossRef]
  27. Avron, J.E. Odd viscosity. J. Stat. Phys. 1998, 92, 543–557. [Google Scholar] [CrossRef]
  28. Avron, J.E.; Elgart, A. Adiabatic theorem without a gap condition: Two-level system coupled to quantized radiation field. Phys. Rev. A 1998, 58, 4300. [Google Scholar] [CrossRef]
  29. Sumino, Y.; Nagai, K.H.; Shitaka, Y.; Tanaka, D.; Yoshikawa, K.; Chate, H.; Oiwa, K. Large-scale vortex lattice emerging from collectively moving micro-tubules. Nature 2012, 483, 448–452. [Google Scholar] [CrossRef]
  30. Tsai, J.C.; Ye, F.; Rodriguez, J.; Gollub, J.P.; Lubensky, T.C. A chiral granular gas. Phys. Rev. Lett. 2005, 94, 214301. [Google Scholar] [CrossRef]
  31. Maggi, C.; Saglimbeni, F.; Dipalo, M.; Angelis, F.D.; Leonardo, R.D. Micromotors with asymmetric shape that efficiently convert light into work by thermocapillary effects. Nat. Commun. 2015, 6, 7855. [Google Scholar] [CrossRef]
  32. Kirkinis, E.; Andreev, A.V. Odd-viscosity-induced stabilization of viscous thin liquid films. J. Fluid Mech. 2019, 878, 169–189. [Google Scholar] [CrossRef]
  33. Lapa, M.F.; Hughes, T.L. Swimming at low reynolds number in fluids with odd, or hall, viscosity. Phys. Rev. E 2014, 89, 043019. [Google Scholar] [CrossRef] [PubMed]
  34. Zhao, J.; Jian, Y. Effect of odd viscosity on the stability of a falling thin film in presence of electromagnetic field. Fluid Dyn. Res. 2021, 53, 015510. [Google Scholar] [CrossRef]
  35. Mukhopadhyay, S.; Mukhopadhyay, A. Hydrodynamics and instabilities of falling liquid film over a non-uniformly heated inclined wavy bottom. Phys. Fluids 2020, 32, 074103. [Google Scholar] [CrossRef]
  36. Wierschem, A.; Lepski, C.; Aksel, N. Effect of long undulated bottoms on thin gravity-driven films. Acta Mech. 2005, 179, 41–66. [Google Scholar] [CrossRef]
  37. Miladinova, S.; Slavtchev, S.; Lebon, G.; Legros, J.C. Long-wave instabilities of non-uniformly heated falling films. J. Fluid Mech. 2002, 453, 153–175. [Google Scholar] [CrossRef]
  38. Dandapat, B.S.; Mukhopadhyay, A. Finite amplitude long wave instability of a film of conducting fluid flowing down an inclined plane in presence of electromagnetic field. Int. J. Appl. Mech. Eng. 2003, 8, 379–393. [Google Scholar]
  39. Mukhopadhyay, A.; Dandapat, B.S. Nonlinear stability of conducting viscous film flowing down an inclined plane at moderate Reynolds number in the presence of a uniform normal electric field. J. Phys. D Appl. Phys. 2004, 38, 138–143. [Google Scholar] [CrossRef]
  40. Mukhopadhyay, A.; Dandapat, B.S.; Mukhopadhyay, A. Stability of conducting liquid flowing down an inclined plane at moderate Reynolds number in the presence of constant electromagnetic field. Int. J. Non-Linear Mech. 2008, 43, 632–642. [Google Scholar] [CrossRef]
  41. Chattopadhyay, S.; Subedar, G.Y.; Gaonkar, A.K.; Barua, A.K.; Mukhopadhyay, A. Effect of odd-viscosity on the dynamics and stability of a thin liquid film flowing down on a vertical moving plate. Int. J. Non-Linear Mech. 2022, 140, 103905. [Google Scholar] [CrossRef]
  42. Mukhopadhyay, A.; Haldar, S. Long-wave instabilities of viscoelastic fluid film flowing down an inclined plane with linear temperature variation. Z. Naturforsch. A 2010, 65, 618–632. [Google Scholar] [CrossRef]
Figure 1. The schematic diagram of the physical model.
Figure 1. The schematic diagram of the physical model.
Nanomaterials 13 02660 g001
Figure 2. Temporal growth rate curve derived from ω i = B 1 k 2 C 1 k 4 , (a) x = 10 π / 3 ; (b) x = 5 π / 3 , for different values of Marangoni number when R e = 2 , ζ = 0.1 π , S = 4.5 , P r = 7 , μ = 0.4 , β = π / 3 .
Figure 2. Temporal growth rate curve derived from ω i = B 1 k 2 C 1 k 4 , (a) x = 10 π / 3 ; (b) x = 5 π / 3 , for different values of Marangoni number when R e = 2 , ζ = 0.1 π , S = 4.5 , P r = 7 , μ = 0.4 , β = π / 3 .
Nanomaterials 13 02660 g002
Figure 3. Temporal growth rate curve derived from ω i = B 1 k 2 C 1 k 4 , (a) x = 10 π / 3 ; (b) x = 5 π / 3 , for different values of odd viscosity when R e = 2 , ζ = 0.1 π , S = 4.5 , P r = 7 , M n = 0.4 , β = π / 3 .
Figure 3. Temporal growth rate curve derived from ω i = B 1 k 2 C 1 k 4 , (a) x = 10 π / 3 ; (b) x = 5 π / 3 , for different values of odd viscosity when R e = 2 , ζ = 0.1 π , S = 4.5 , P r = 7 , M n = 0.4 , β = π / 3 .
Nanomaterials 13 02660 g003
Figure 4. Critical Reynolds number as a function of bottom steepness, (a) x = 10 π / 3 ; (b) x = 5 π / 3 , for different values of odd viscosity when R e = 2 , ζ = 0.1 π , S = 4.5 , P r = 7 , M n = 0.4 , β = π / 3 .
Figure 4. Critical Reynolds number as a function of bottom steepness, (a) x = 10 π / 3 ; (b) x = 5 π / 3 , for different values of odd viscosity when R e = 2 , ζ = 0.1 π , S = 4.5 , P r = 7 , M n = 0.4 , β = π / 3 .
Nanomaterials 13 02660 g004
Figure 5. Comparison of R e c as a function of bottom steepness for isothermal bottom at a point on the “downhill” ( x = 5 π / 3 ) and at a point on the “uphill” ( x = 10 π / 3 ) portion when B o = 1 , P r = 7 , β = π / 3 .
Figure 5. Comparison of R e c as a function of bottom steepness for isothermal bottom at a point on the “downhill” ( x = 5 π / 3 ) and at a point on the “uphill” ( x = 10 π / 3 ) portion when B o = 1 , P r = 7 , β = π / 3 .
Nanomaterials 13 02660 g005
Figure 6. Stability curves, (a) μ = 0 ; (b) μ = 0.4 ; (c) μ = 0.8 ; (d) μ = 1.2 , for different odd viscosity μ when ζ = 0.1 π , S = 4.5 , P r = 7 , M n = 0.4 , β = π / 3 . I–IV represent the supercritical stability region, the subcritical instability region, unconditional stability region and the explosive region, respectively.
Figure 6. Stability curves, (a) μ = 0 ; (b) μ = 0.4 ; (c) μ = 0.8 ; (d) μ = 1.2 , for different odd viscosity μ when ζ = 0.1 π , S = 4.5 , P r = 7 , M n = 0.4 , β = π / 3 . I–IV represent the supercritical stability region, the subcritical instability region, unconditional stability region and the explosive region, respectively.
Nanomaterials 13 02660 g006
Figure 7. Stability curves, (a) M n = 0 ; (b) M n = 0.2 ; (c) M n = 0.4 ; (d) M n = 0.6 , for different Marangoni number M n when ζ = 0.1 π , S = 4.5 , P r = 7 , μ = 0.2 , β = π / 3 . I–IV represent the supercritical stability region, the subcritical instability region, unconditional stability region and the explosive region, respectively.
Figure 7. Stability curves, (a) M n = 0 ; (b) M n = 0.2 ; (c) M n = 0.4 ; (d) M n = 0.6 , for different Marangoni number M n when ζ = 0.1 π , S = 4.5 , P r = 7 , μ = 0.2 , β = π / 3 . I–IV represent the supercritical stability region, the subcritical instability region, unconditional stability region and the explosive region, respectively.
Nanomaterials 13 02660 g007aNanomaterials 13 02660 g007b
Figure 8. Film thickness at different times, (a) μ = 0 ; (b) μ = 0.4 , when R e = 2 , k = 1.75 , ζ = 0.1 π , S = 4.5 , P r = 7 , M n = 0.5 , β = π / 3 .
Figure 8. Film thickness at different times, (a) μ = 0 ; (b) μ = 0.4 , when R e = 2 , k = 1.75 , ζ = 0.1 π , S = 4.5 , P r = 7 , M n = 0.5 , β = π / 3 .
Nanomaterials 13 02660 g008
Figure 9. Film thickness at different times, (a) M n = 0 ; (b) M n = 0.4 , when R e = 1 , k = 1.7 , ζ = 0.1 π , S = 4.5 , P r = 7 , μ = 0.2 , β = π / 3 .
Figure 9. Film thickness at different times, (a) M n = 0 ; (b) M n = 0.4 , when R e = 1 , k = 1.7 , ζ = 0.1 π , S = 4.5 , P r = 7 , μ = 0.2 , β = π / 3 .
Nanomaterials 13 02660 g009
Figure 10. Variations of waves with bottom steepness through numerical simulation, (a) when R e = 1 , k = 1.5 , S = 4.5 , P r = 7 , μ = 0.2 , M n = 0.4 , β = π / 3 ; (b) R e = 3 , k = 1.48 , S = 4.5 , P r = 7 , μ = 0.2 , M n = 0.4 , β = π / 3 .
Figure 10. Variations of waves with bottom steepness through numerical simulation, (a) when R e = 1 , k = 1.5 , S = 4.5 , P r = 7 , μ = 0.2 , M n = 0.4 , β = π / 3 ; (b) R e = 3 , k = 1.48 , S = 4.5 , P r = 7 , μ = 0.2 , M n = 0.4 , β = π / 3 .
Nanomaterials 13 02660 g010
Figure 11. Variations of waves with different odd viscosity and Marangoni number through numerical simulations, (a) odd viscosity μ when R e = 3 , k = 1.4 , S = 4.5 , P r = 7 , M n = 0.5 , β = π / 3 ; (b) Marangoni number M n when R e = 3 , k = 0.67 , S = 4.5 , P r = 7 , μ = 0.2 , β = π / 3 .
Figure 11. Variations of waves with different odd viscosity and Marangoni number through numerical simulations, (a) odd viscosity μ when R e = 3 , k = 1.4 , S = 4.5 , P r = 7 , M n = 0.5 , β = π / 3 ; (b) Marangoni number M n when R e = 3 , k = 0.67 , S = 4.5 , P r = 7 , μ = 0.2 , β = π / 3 .
Nanomaterials 13 02660 g011
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Xue, D.; Zhang, R.; Liu, Q.; Ding, Z. Instability of Liquid Film with Odd Viscosity over a Non-Uniformly Heated and Corrugated Substrate. Nanomaterials 2023, 13, 2660. https://doi.org/10.3390/nano13192660

AMA Style

Xue D, Zhang R, Liu Q, Ding Z. Instability of Liquid Film with Odd Viscosity over a Non-Uniformly Heated and Corrugated Substrate. Nanomaterials. 2023; 13(19):2660. https://doi.org/10.3390/nano13192660

Chicago/Turabian Style

Xue, Danting, Ruigang Zhang, Quansheng Liu, and Zhaodong Ding. 2023. "Instability of Liquid Film with Odd Viscosity over a Non-Uniformly Heated and Corrugated Substrate" Nanomaterials 13, no. 19: 2660. https://doi.org/10.3390/nano13192660

APA Style

Xue, D., Zhang, R., Liu, Q., & Ding, Z. (2023). Instability of Liquid Film with Odd Viscosity over a Non-Uniformly Heated and Corrugated Substrate. Nanomaterials, 13(19), 2660. https://doi.org/10.3390/nano13192660

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop