A New Upper Bound for the Largest Growth Rate of Linear Rayleigh--Taylor Instability

TL;DR

This paper establishes a new upper bound for the maximum growth rate of linear Rayleigh--Taylor instability in viscous fluids, demonstrating how surface tension influences the instability's development and confirming classical experimental observations.

Contribution

It introduces a novel upper bound for the growth rate and analyzes how surface tension affects the linear RT instability, extending understanding of the instability's limiting behavior.

Findings

Derived a new upper bound for the growth rate $\\Lambda$

Showed that $\\Lambda$ decreases to zero as surface tension approaches a critical threshold

Mathematically verified the limiting effect of surface tension on instability growth

Abstract

We investigate the effect of surface tension on the linear Rayleigh--Taylor (RT) instability in stratified incompressible viscous fluids with or without (interface) surface tension. The existence of linear RT instability solutions with largest growth rate $\Lambda$ is proved under the instability condition (i.e., the surface tension coefficient $\vartheta$ is less than a threshold $\vartheta_{\rm c}$) by modified variational method of PDEs. Moreover we find a new upper bound for $\Lambda$. In particular, we observe from the upper bound that $\Lambda$ decreasingly converges to zero, as $\vartheta$ goes from zero to the threshold $\vartheta_{\rm c}$. The convergence behavior of $\Lambda$ mathematically verifies the classical RT instability experiment that the instability growth is limited by surface tension during the linear stage.

Equations

$$\left\{{\begin{array}[]{ll}\rho_{\pm}(\partial_{t}v_{\pm}+v_{\pm}\cdot\nabla v_{\pm})+\mathrm{div}\mathcal{S}=-g\rho_{\pm}e_{3}&\mbox{ in }\Omega_{\pm}(t),\\[2.84526pt] \mathrm{div}v_{\pm}=0&\mbox{ on }\Omega_{\pm}(t),\\ d_{t}+v_{1}\partial_{1}d+v_{2}\partial_{2}d=v_{3}&\mbox{ on }\Sigma(t),\\ \llbracket v_{\pm}\rrbracket=0,\ \llbracket\mathcal{S}_{\pm}\nu\rrbracket=\vartheta\mathcal{C}\nu&\mbox{ on }\Sigma(t),\\ v_{\pm}=0&\mbox{ on }\Sigma_{\pm},\\ v_{\pm}|_{t=0}=v^{0}_{\pm}&\mbox{ in }\Omega_{\pm}(0),\\ d|_{t=0}=d^{0}&\mbox{ on }\Sigma(0).\end{array}}\right.$$

$${\Omega }_{+}(t):=\{(x_{\mathrm{h}},x_3)|x_{\mathrm{h}}:=(x_1,x_2)\in \mathbb{T},d(x_{\mathrm{h}},t)<x_3<h_{+}\},$$

$${\Omega }_{-}(t):=\{(x_{\mathrm{h}},x_3)|x_{\mathrm{h}}\in \mathbb{T},-h_{-}<x_3<d(x_{\mathrm{h}},t)\},$$

$$\Sigma (t):=\{(x_{\mathrm{h}},x_3)|x_{\mathrm{h}}\in \mathbb{T},x_3:=d(x_{\mathrm{h}},t)\}.$$

$${\mathcal{S}}_{\pm }:=p_{\pm }I-{\mu }_{\pm }\mathbb{D}v\text{mbox}with\mathbb{D}v=\nabla v_{\pm }+\nabla v_{\pm }^{\mathrm{T}}.$$

$$\mathcal{C}:=\frac{{\Delta }_{\mathrm{h}}d+({\partial }_{1}d{)}^2{\partial }_2^{2}d+({\partial }_{2}d{)}^2{\partial }_1^{2}d-2{\partial }_{1}d{\partial }_{2}d{\partial }_1{\partial }_{2}d}{(1+({\partial }_{1}d{)}^2+({\partial }_{2}d{)}^2{)}^{3/2}}.$$

$$\rho ={\rho }_{+}{\chi }_{{\Omega }_{+}(t)}+{\rho }_{-}{\chi }_{{\Omega }_{-}(t)},\mu ={\mu }_{+}{\chi }_{{\Omega }_{+}(t)}+{\mu }_{-}{\chi }_{{\Omega }_{-}(t)},$$

$$v=v_{+}{\chi }_{{\Omega }_{+}(t)}+v_{-}{\chi }_{{\Omega }_{-}(t)},p=p_{+}{\chi }_{{\Omega }_{+}(t)}+p_{-}{\chi }_{{\Omega }_{-}(t)},$$

$$v^0=v_{+}^0{\chi }_{{\Omega }_{+}(0)}+v_{-}^0{\chi }_{{\Omega }_{-}(0)},\mathcal{S}:=pI-\mu \mathbb{D}v.$$

$$\left\{{\begin{array}[]{ll}\rho(v_{t}+v\cdot\nabla v)+\mathrm{div}\mathcal{S}=-g\rho e_{3}&\mbox{ in }\Omega(t),\\[2.84526pt] \mathrm{div}v=0&\mbox{ in }\Omega(t),\\[2.84526pt] d_{t}+v_{1}\partial_{1}d+v_{2}\partial_{2}d=v_{3}&\mbox{ on }\Sigma(t),\\ \llbracket v\rrbracket=0,\ \llbracket\mathcal{S}\nu\rrbracket=\vartheta\mathcal{C}\nu&\mbox{ on }\Sigma(t),\\ v=0&\mbox{ on }\Sigma_{-}^{+},\\ v|_{t=0}=v^{0}&\mbox{ in }\Omega(0),\\ d|_{t=0}=d^{0}&\mbox{ on }\Sigma(0),\end{array}}\right.$$

$${\Omega }_{+}:=\{(y_{\mathrm{h}},y_3)\in {\mathbb{R}}^3|y_{\mathrm{h}}\in \mathbb{T},0<y_3<h_{+}\},$$

$${\Omega }_{-}:=\{(y_{\mathrm{h}},y_3)\in {\mathbb{R}}^3|y_{\mathrm{h}}\in \mathbb{T},-h_{-}<y_3<0\},$$

$${\zeta }_{\pm }^0:{\Omega }_{\pm }\to {\Omega }_{\pm }(0),$$

$$\Sigma (0)={\zeta }_{\pm }^0(\Sigma ),{\Sigma }_{\pm }={\zeta }_{\pm }^0({\Sigma }_{\pm })\text{mbox}and\det \nabla {\zeta }_{\pm }^0=1.$$

$$\left\{\begin{array}[]{ll}\partial_{t}\zeta(y,t)=v(\zeta(y,t),t)&\mbox{ in }\Omega_{-}^{+}\\ \zeta(y,0)=\zeta^{0}(y)&\mbox{ in }\Omega_{-}^{+},\end{array}\right.$$

$$y=\zeta (y,t)\text{mbox}on{\Sigma }_{-}^{+}.$$

$$(u,\sigma )(y,t)=(v,p+g{x}_3)(\zeta (y,t),t)\text{mbox}for(y,t)\in {\Omega }_{-}^{+}\times (0,\infty ),$$

$$\left\{\begin{array}[]{ll}\zeta_{t}=u&\mbox{ in }\Omega_{-}^{+},\\[2.84526pt] \rho u_{t}+\nabla_{\mathcal{A}}\sigma-\mu\Delta_{\mathcal{A}}u=0&\mbox{ in }\Omega_{-}^{+},\\[2.84526pt] \mathrm{div}_{\mathcal{A}}u=0&\mbox{ in }\Omega_{-}^{+},\\[2.84526pt] \llbracket\zeta\rrbracket=\llbracket u\rrbracket=0,\ \llbracket((\sigma-g\zeta_{3})I-\mathbb{D}_{\mathcal{A}}u)\vec{n}\rrbracket=\vartheta\mathcal{H}\vec{n}&\mbox{ on }\Sigma,\\ (\zeta,u)=(y,0)&\mbox{ on }\Sigma_{-}^{+},\\ (\zeta,u)|_{t=0}=(\eta^{0},u^{0}),&\mbox{ in }\Omega_{-}^{+},\end{array}\right.$$

$$\Sigma :=\mathbb{T}\times \{0\},\vec{n}:=\mathcal{A}{e}_3/|\mathcal{A}{e}_3|,{\mathbb{D}}_{\mathcal{A}}u={\nabla }_{\mathcal{A}}u+{\nabla }_{\mathcal{A}}{u}^{\mathrm{T}},$$

$$\mathcal{H}:=|{\partial }_1\zeta {|}^2{\partial }_2^2\zeta -2({\partial }_1\zeta \cdot {\partial }_2\zeta ){\partial }_1{\partial }_2\zeta +|{\partial }_2\zeta {|}^2{\partial }_1^2\zeta \cdot \vec{n}/(|{\partial }_1\zeta {|}^2|{\partial }_2\zeta {|}^2-|{\partial }_1\zeta \cdot {\partial }_2\zeta {|}^2).$$

$${\mathcal{A}}^{\mathrm{T}}=(\nabla \zeta {)}^{-1}:=({\partial }_j{\zeta }_i{)}_{3\times 3}^{-1},$$

$${\nabla }_{\mathcal{A}}w:=({\nabla }_{\mathcal{A}}{w}_1,{\nabla }_{\mathcal{A}}{w}_2,{\nabla }_{\mathcal{A}}{w}_3{)}^{\mathrm{T}}\text{mbox}and{\nabla }_{\mathcal{A}}{w}_i:=({\mathcal{A}}_{1k}{\partial }_k{w}_i,{\mathcal{A}}_{2k}{\partial }_k{w}_i,{\mathcal{A}}_{3k}{\partial }_k{w}_i{)}^{\mathrm{T}}$$

$${\mathrm{div}}_{\mathcal{A}}(f^1,f^2,f^3)=({\mathrm{div}}_{\mathcal{A}}{f}^1,{\mathrm{div}}_{\mathcal{A}}{f}^2,{\mathrm{div}}_{\mathcal{A}}{f}^3{)}^{\mathrm{T}}\text{mbox}and{\mathrm{div}}_{\mathcal{A}}{f}^i:={\mathcal{A}}_{lk}{\partial }_k{f}_l^i$$

$$\left\{\begin{array}[]{ll}\nabla\bar{p}_{\pm}=-\bar{\rho}_{\pm}ge_{3}&\mbox{ in }\Omega_{\pm},\\[2.84526pt] \llbracket\bar{p}_{\pm}\rrbracket e_{3}=0&\mbox{ on }\Sigma,\end{array}\right.$$

$$[[\rho ]]>0\text{mbox}on\Sigma .$$

$$\eta :=\zeta -y,u=u-0\text{mbox}andq=\sigma -(\bar{p}({\zeta }_3)+g\rho {\zeta }_3),$$

$$\left\{\begin{array}[]{ll}\eta_{t}=u&\mbox{ in }\Omega_{-}^{+},\\[2.84526pt] {\rho}u_{t}+\nabla q-\mu\Delta u=\mathcal{N}_{1}&\mbox{ in }\Omega_{-}^{+},\\[2.84526pt] \mathrm{div}u=\mathcal{N}_{2}&\mbox{ in }\Omega_{-}^{+},\\[2.84526pt] \llbracket\eta\rrbracket=\llbracket u\rrbracket=0,\ \llbracket((q-g\rho\eta_{3})I-\mu\mathbb{D}u)e_{3}\rrbracket=\vartheta\Delta_{\mathrm{h}}\eta_{3}e_{3}+\mathcal{N}_{3}&\mbox{ on }\Sigma,\\ (\eta,u)=0&\mbox{ on }\Sigma_{-}^{+},\\ (\eta,u)|_{t=0}=(\eta^{0},u^{0})&\mbox{ in }\Omega_{-}^{+},\end{array}\right.$$

$${\mathcal{N}}_1={\mathrm{div}}_{\tilde{\mathcal{A}}}{\nabla }_{\mathcal{A}}u+\mathrm{div}{\nabla }_{\tilde{\mathcal{A}}}u-{\nabla }_{\tilde{\mathcal{A}}}q,{\mathcal{N}}_2=-{\mathrm{div}}_{\tilde{\mathcal{A}}}u,$$

$${\mathcal{N}}_3={\mathbb{D}}_{\tilde{\mathcal{A}}}u+\vartheta \mathcal{H}\vec{n}-\vartheta {\Delta }_{\mathrm{h}}{\eta }_3{e}_3.$$

$$\left\{\begin{array}[]{ll}\eta_{t}=u&\mbox{ in }\Omega_{-}^{+},\\[2.84526pt] \rho u_{t}+\nabla q-\mu\Delta u=0&\mbox{ in }\Omega_{-}^{+},\\[2.84526pt] {\rm div}u=0&\mbox{ in }\Omega_{-}^{+},\\ \llbracket\eta\rrbracket=\llbracket u\rrbracket=0&\mbox{ on }\Sigma,\\ \llbracket((q-g\rho\eta_{3})I-\mu\mathbb{D}u)e_{3}\rrbracket=\vartheta\Delta_{\mathrm{h}}\eta_{3}e_{3}&\mbox{ on }\Sigma,\\ (\eta,u)=0&\mbox{ on }\Sigma_{-}^{+},\\ (\eta,u)|_{t=0}=(\eta^{0},u^{0})&\mbox{ in }\Omega_{-}^{+}.\end{array}\right.$$

$${\Lambda }_{\vartheta }⩽h_{-}g[[\rho ]]/4{\mu }_{-}\text{mbox}for{h}_{+}=1.$$

$${\Lambda }_{\vartheta }⩽m:=\min \left\{\frac{({\vartheta }_{\mathrm{c}}-\vartheta )}{4\max \{L_1^2,L_2^2\}}\min \left\{\frac{h_{+}}{{\mu }_{+}},\frac{h_{-}}{{\mu }_{-}}\right\},{\left(\frac{4(g[[\rho ]]({\vartheta }_{\mathrm{c}}-\vartheta ){)}^2}{{\vartheta }_{\mathrm{c}}^2\max \{{\rho }_{+}{\mu }_{+},{\rho }_{-}{\mu }_{-}\}}\right)}^{\frac{1}{3}}\right\}.$$

$$m⩽h_{-}g[[\rho ]]/4{\mu }_{-}.$$

$${\Lambda }_{\vartheta }\to 0\text{mbox}as\vartheta \to {\vartheta }_{\mathrm{c}}.$$

Full text

1 Introduction

Considering two completely plane-parallel layers of stratified (immiscible) fluids, the heavier one on top of the lighter one and both subject to the earth’s gravity, it is well-known that such equilibrium state is unstable to sustain small disturbances, and this unstable disturbance will grow and lead to a release of potential energy, as the heavier fluid moves down under the gravitational force, and the lighter one is displaced upwards. This phenomenon was first studied by Rayleigh [26] and then Taylor [27], and is called therefore the Rayleigh–Taylor (RT) instability. In the last decades, this phenomenon has been extensively investigated from mathematical, physical and numerical aspects, see [2, 29, 8] for instance. It has been also widely investigated how the RT instability evolves under the effects of other physical factors, such as elasticity [21, 23, 11, 30, 3], rotation [2, 4], (internal) surface tension [9, 34, 14], magnetic fields [20, 16, 17, 19, 31, 32, 22, 18] and so on. In this article, we are interested in the effect of surface tension on the linear RT instability in stratified incompressible viscous fluids. To conveniently introduce relevant mathematical progress and our main results, next we shall mathematically formulate our problem in details.

1.1 Motion equations in Eulerian coordinates

Let us first recall a mathematical model, which describes the horizontally periodic motion of stratified incompressible viscous fluids in an infinity layer domain [21]:

[TABLE]

The momentum equations in (1.1)1 describe the motion of the both upper heavier and lower lighter viscous fluids driven by the gravitational field along the negative $x_{3}$-direction, which occupy the two time-dependent disjoint open subsets $\Omega_{+}(t)$ and $\Omega_{-}(t)$ at time $t$, respectively. Moreover the fluids are incompressible due to (1.1)2. The two fluids interact with each other by the motion equation of a free interface (1.1)3 and the interfacial jump conditions in (1.1)4. The first jump condition in (1.1)4 represents that the velocity is continuous across the interface. The second jump in (1.1)4 represents that the jump in the normal stress is proportional to the mean curvature of the surface multiplied by the normal to the surface. The non-slip boundary condition of the velocities on the both upper and lower fixed flat boundaries are described by (1.1)5. (1.1)6 and (1.1)7 represents the initial status of the two fluids. Next we shall further explain the notations in (1.1) in details.

The subscripts $+$ resp. $-$ in the notations $f_{+}$ resp. $f_{-}$ mean that functions, parameters or domains $f_{+}$ resp. $f_{-}$ are relevant to the upper resp. lower fluids. For each given $t>0$, $d:=d(x_{\mathrm{h}},t):\mathbb{T}\mapsto(-h_{-},h_{+})$ is a height function of a point at the interface of stratified fluids, where $h_{-}$, $h_{+}>0$, $\mathbb{T}:=\mathbb{T}_{1}\times\mathbb{T}_{2}$, $\mathbb{T}_{i}=2\pi L_{i}(\mathbb{R}/\mathbb{Z})$, and $2\pi L_{i}$ ($i=1$, $2$) are the periodicity lengths. The domains $\Omega_{\pm}(t)$ and the interface $\Sigma(t)$ are defined as follows:

[TABLE]

In addition, $\Sigma_{+}=\mathbb{T}\times\{h_{+}\}$, $\Sigma_{-}=\mathbb{T}\times\{-h_{-}\}$ and we call $\Omega:=\mathbb{T}\times(-h_{-},h_{+})$ the domain of stratified fluids.

For given $t>0$, $v_{\pm}(x,t):\Omega_{\pm}(t)\mapsto\mathbb{R}^{3}$ are the velocities of the two fluids, and $\mathcal{S}_{\pm}$ the stress tensors enjoying the following expression:

[TABLE]

In the above expression the superscript $\mathrm{T}$ means matrix transposition and $I$ is the $3\times 3$ identity matrix. $\rho_{\pm}$ are the density constants, and the constants $\mu_{\pm}>0$ the shear viscosity coefficients. $g$ and $\vartheta$ represent the gravitational constant and the surface tension coefficient, reps. In addition, $e_{3}:=(0,0,1)^{\mathrm{T}}$.

For a function $f$ defined on $\Omega(t)$, we define $\llbracket f_{\pm}\rrbracket:=f_{+}|_{\Sigma(t)}-f_{-}|_{\Sigma(t)}$, where $f_{\pm}|_{\Sigma(t)}$ are the traces of the quantities $f_{\pm}$ on $\Sigma(t)$. $\nu$ is the unit outer normal vector at boundary $\Sigma(t)$ of $\Omega_{-}(t)$, and $\mathcal{C}$ the twice of the mean curvature of the internal surface $\Sigma(t)$, i.e.,

[TABLE]

Now we further introduce the indicator function $\chi_{\Omega_{\pm}(t)}$ and denote

[TABLE]

then the model (1.1) can be rewritten as follows:

[TABLE]

where we have defined that $\Omega(t):=\Omega_{+}(t)\cup\Omega_{-}(t)$, $\Sigma_{-}^{+}:=\Sigma_{-}\cup\Sigma_{+}$ and omitted the subscript $\pm$ in $\llbracket f_{\pm}\rrbracket$ for simplicity.

1.2 Reformulation in Lagrangian coordinates

Next we adopt the transformation of Lagrangian coordinates so that the interface and the domains stay fixed in time.

We define that

[TABLE]

and assume that there exist invertible mappings

[TABLE]

such that

[TABLE]

We further define $\zeta^{0}:=\zeta_{+}^{0}\chi_{\Omega_{+}}+\zeta_{-}^{0}\chi_{\Omega_{-}}$, and the flow map $\zeta$ as the solution to

[TABLE]

where $\Omega_{-}^{+}:=\Omega_{+}\cup\Omega_{-}$. We denote the Eulerian coordinates by $(x,t)$ with $x=\zeta(y,t)$, whereas the fixed $(y,t)\in\Omega_{-}^{+}\times\mathbb{R}^{+}$ stand for the Lagrangian coordinates.

In order to switch back and forth from Lagrangian to Eulerian coordinates, we shall assume that $\zeta_{\pm}(\cdot,t)$ are invertible and $\Omega_{\pm}(t)=\zeta_{\pm}(\Omega_{\pm},t)$, and since $v_{\pm}$ and $\zeta_{\pm}^{0}$ are all continuous across $\Sigma$, we have $\Sigma(t)=\zeta_{\pm}(\Sigma,t)$. In view of the non-slip boundary condition $v|_{\Sigma_{-}^{+}}=0$, we have

[TABLE]

Now we set the Lagrangian unknowns

[TABLE]

then the problem (1.3) can be rewritten as an initial-boundary value problem with an interface for $(\zeta,u)$ in Lagrangian coordinates:

[TABLE]

where we have defined that

[TABLE]

We shall introduce the notations involving $\mathcal{A}$. The matrix $\mathcal{A}:=(\mathcal{A}_{ij})_{3\times 3}$ is defined via

[TABLE]

where $\partial_{j}$ denote the partial derivative with respect to the $j$-th components of variables $y$. $\tilde{\mathcal{A}}:=\mathcal{A}-I$, and $I$ is the $3\times 3$ identity matrix. The differential operator $\nabla_{\mathcal{A}}$ is defined by

[TABLE]

for vector function $w:=(w_{1},w_{2},w_{3})$, and the differential operator $\mathrm{div}_{\mathcal{A}}$ is defined by

[TABLE]

for vector function $f^{i}:=(f_{1}^{i},f_{2}^{i},f_{3}^{i})^{\mathrm{T}}$. It should be noted that we have used the Einstein convention of summation over repeated indices. In addition, we define $\Delta_{\mathcal{A}}X:=\mathrm{div}_{\mathcal{A}}\nabla_{\mathcal{A}}X$.

1.3 Linearized motion

We choose a constant $\bar{d}\in(-l,\tau)$. Without loss of generality, we assume that $\bar{d}=0$. Then we consider an RT equilibrium state

[TABLE]

where $\rho$ satisfies the RT (jump) condition

[TABLE]

Let $\bar{p}:=\bar{p}_{+}\chi_{\Omega_{+}}+\bar{p}_{-}\chi_{\Omega_{-}}$. Then $(v,p)=(0,\bar{p})$ with $d=0$ is an RT equilibria solution of (1.3).

Denoting the perturbation in Lagrangian coordinates

[TABLE]

then subtracting (1.7) from (1.6) yields the perturbation RT problem in Lagrangian coordinates:

[TABLE]

where $\Delta_{\mathrm{h}}:=\partial_{1}^{2}+\partial_{2}^{2}$ and the nonlinear terms $\mathcal{N}_{1}$–$\mathcal{N}_{3}$ are defined as follows:

[TABLE]

Omitting the nonlinear terms in (1.9), we get a linearized RT problem:

[TABLE]

Of course, the motion equations of stratified viscous fluids in linear stage can be approximatively described by (1.10).

The inhibition of RT instability by surface tension was first analyzed by Bellman–Phennington [1] based on a linearized two-dimensional (2D) motion equations of stratified incompressible inviscid fluids defined on the domain $2\pi L_{1}\mathbb{T}_{1}\times(-h_{-},h_{+})$ (i.e., $\mu=0$ in the corresponding 2D case of (1.10)) in 1953. More precisely, they proved that the linear 2D stratified incompressible inviscid fluids is stable, resp. unstable for $\vartheta>g\llbracket\rho\rrbracket L_{1}^{2}$, resp. $\vartheta<g\llbracket\rho\rrbracket L_{1}^{2}$. The value $g\llbracket\rho\rrbracket L_{1}^{2}$ is a threshold of surface tension coefficient for linear stability and linear instability. Similar result was also found in the 3D viscous case, for example, Guo–Tice proved that $\vartheta_{\mathrm{c}}:={g\llbracket\rho\rrbracket\max\{L_{1}^{2},L_{2}^{2}\}}$ is a threshold of surface tension coefficient for stability and instability in the linearized 3D stratified compressible viscous fluids defined on $\Omega$ [9]. Next we further review the mathematical progress for the nonlinear case.

Pr$\mathrm{\ddot{u}}$ess–Simonett first proved that the RT equilibria solution of the stratified incompressible viscid fluids defined on the domain $\mathbb{R}^{3}$ is unstable based on a Henry instability method [25]. Later, Wang–Tice–Kim verified that the RT equilibria solution of stratified incompressible viscous fluids defined on $\Omega$ is stable, resp. unstable for $\vartheta>\vartheta_{\mathbb{T}}$, resp. $\vartheta\in[0,\vartheta_{\mathbb{T}})$ [34, 33]. Jang–Wang–Tice further obtained the same results of stability and instability in the corresponding compressible case [14, 13]. Recently, Wilke also proved there exists a threshold $\vartheta_{\mathrm{c}}$ for the stability and instability of stratified viscous fluids (with heavier fluid over lighter fluid) defined on a cylindrical domain with finite height [35]. Finally, we mention that the results of nonlinear RT instability in inhomogeneous fluid (without interface) were obtained based on the classical bootstrap instability method, see [12], resp. [15] for inviscid, resp. viscous cases.

2 Main result

In this paper, we investigate the effect of surface tension on the linear RT instability by the linearized motion (1.10). Wang–Tice used discrete Fourier transformation and modified variational method of ODEs to prove the existence of growth solutions with a largest growth rate $\Lambda_{\vartheta}$ for (1.10) with $h_{+}=1$ under the condition $\vartheta_{\vartheta}\in(0,\vartheta_{\mathbb{T}})$ [33]. Moreover, they provided an upper bound for $\Lambda_{\vartheta}$:

[TABLE]

In this paper, we exploit modified variational method of PDEs and existence theory of stratified (steady) Stokes problem to prove the existence of growth solutions with a largest growth rate $\Lambda_{\vartheta}$ for (1.10) under the instability condition $\vartheta\in[0,\vartheta_{\mathbb{T}})$. Moreover we find a new upper bound:

[TABLE]

It is easy to see that

[TABLE]

Therefore our upper bound is more precise than Wang–Tice’s one. Moreover, we see from (2.1) that

[TABLE]

In classical Rayleigh–Taylor (RT) experiments [6, 10], it has been shown the phenomenon of that the instability growth is limited by surface tension during the linear stage, where the growth is exponential in time. Obviously, the convergence behavior (2.2) mathematically verifies the phenomenon.

Before stating our main results in detail, we shall introduce some simplified notations throughout this article.

Basic notations: $I_{T}:=(0,T)$. $\mathbb{R}_{+}:=(0,\infty)$. The $j$-th difference quotient of size $h$ is $D_{j}^{h}w:=(w(y+he_{j})-w(y))/h$ for $j=1$ and $2$, and $D^{h}_{\mathrm{h}}w:=(D_{1}^{h}w_{1},D_{2}^{h}w_{2})$, where $|h|\in(0,1)$. $\Re f$, reps. $\Im f$ denote the real, resp. imaginary parts of the complex function $f$. $\nabla_{\mathrm{h}}^{k}f$ denotes a $(k+1)\times(k+1)$ matrix $(\partial_{1}^{i}\partial_{2}^{j}f)_{ij}$ for $k\geqslant 0$. $a\lesssim b$ means that $a\leqslant cb$ for some constant $c>0$, where the positive constant $c$ may depend on the domain $\Omega$, and known parameters such as $\rho_{\pm}$, $\mu_{\pm}$, $g$ and $\vartheta$, and may vary from line to line. 2. 2.

Simplified norms: $\|\cdot\|_{i}:=\|\cdot\|_{W^{i,2}}$, $|\cdot|_{s}:=\|\cdot|_{\Sigma}\|_{H^{s}(\mathbb{T})}$, where $s$ is a real number, and $i$ a non-negative integer. 3. 3.

Functionals: $\mathcal{E}(w):=\vartheta|\nabla_{\mathrm{h}}w_{3}|_{0}^{2}-g\llbracket\rho\rrbracket|w_{3}|_{0}^{2}$ and $\mathcal{F}(w,s):=-(\mathcal{E}(w)+s\|\sqrt{\mu}\mathbb{D}w\|^{2}_{0}/2)$.

In addition, we shall give the definition of the largest growth rate of RT instability in the linearized RT problem.

Definition 2.1**.**

We call $\Lambda>0$ the largest growth rate of RT instability in the linearized RT problem (1.10), if it satisfies the following two conditions:

1. (1)

For any strong solution $(\eta,u)\in C^{0}([0,T),H^{3}\cap H^{2}_{\sigma})\cap L^{2}(I_{T},H^{3}\cap H^{3}_{\sigma})$ of the linearized RT problem with $q$ enjoying the regularity $q\in C^{0}([0,T),H^{1})\cap L^{2}(I_{T},H^{2})$, then we have, for any $t\in[0,T)$,

[TABLE] 2. (2)

There exists a strong solution $(\eta,u)$ of the linearized RT problem in the form

[TABLE]

where $(\eta,{u})\in H^{2}$.

Now we state the first result on the existence of largest growth rate in the linearized RT problem.

Theorem 2.1**.**

Let $g>0$, $\rho>0$ and $\mu>0$ are given. Then, for any given

[TABLE]

there is an unstable solution

[TABLE]

to the linearized RT problem (1.10), where $(w,\beta)\in H^{\infty}$ solves the boundary value problem:

[TABLE]

with a largest growth rate $\Lambda>0$ satisfying

[TABLE]

Moreover,

[TABLE]

Next we briefly introduce how to prove Theorem 2.1 by modified variational method of PDEs and regularity theory of stratified (steady) Stokes problem. The detailed proof will be given in Section 4.

We assume a growing mode ansatz to the linearized problem:

[TABLE]

for some $\Lambda>0$. Substituting this ansatz into the linearized RT problem (1.10), we get a spectrum problem

[TABLE]

and then eliminating $\tilde{\eta}$ by using the first equation, we arrive at the boundary-value problem (2.5) for $w$ and $\beta$. Obviously, the linearized RT problem is unstable, if there exists a solution $(w,\beta)$ to the boundary-value problem (2.5) with $\Lambda>0$.

To look for the solution, we use a modified variational method of PDEs, and thus modify (2.5) as follows:

[TABLE]

where $s>0$ is a parameter. To emphasize the dependence of $s$ upon $\alpha$ and $\vartheta$, we will write $\alpha(s,\vartheta)=\alpha$.

Noting that the modified problem (2.8) enjoys the following variational identity

[TABLE]

Thus, by a standard variational approach, there exists a maximizer $w\in\mathcal{A}$ of the functional $\mathcal{F}$ defined on $\mathcal{A}$; moreover $w$ is just a weak solution to (2.8) with $\alpha$ defined by the relation

[TABLE]

see Proposition 4.1. Then we further use the method of difference quotients and the existence theory of the stratified (steady) Stokes problem to improve the regularity of the weak solution, and thus prove that $(w,\beta)\in H^{\infty}$ is a classical solution to the boundary-value problem (2.8), see Proposition 4.2.

In view of the definition of $\alpha(s,\vartheta)$ and the instability condition (2.4), we can infer that, for given $\vartheta$, the function $\alpha(s,\cdot)$ on the variable $s$ enjoys some good properties (see Proposition 4.3), which imply that there exists a $\Lambda$ satisfying the fixed-point relation

[TABLE]

Then we obtain a nontrivial solution $(w,\beta)\in H^{\infty}$ to (2.5) with $\Lambda$ defined by (2.10), and therefore the linear instability follows. Moreover, $\Lambda$ is the largest growth rate of RT instability in the linearized RT problem (see Proposition 4.4), and thus we get Theorem 2.1.

Next we turn to introduce the second result on the properties of largest growth rate constructed by (2.10).

Theorem 2.2**.**

The largest growth rate $\Lambda_{\vartheta}:=\Lambda$ in Theorem 2.1 enjoys the estimate (2.1). Moreover,

[TABLE]

In particular, we have $\Lambda_{\vartheta}\rightarrow 0$ as $\vartheta\rightarrow\vartheta_{\mathrm{c}}$.

The proof of Theorem 2.2 will be presented in Section 5. Here we briefly mention the idea of proof. We find that, for fixed $s$, $\alpha(\cdot,\vartheta)$ defined by (2.9) strictly decreases and is continuous with respect to $\vartheta$ (see Proposition 5.1). Thus, by the fixed-point relation (2.10) and some analysis based on the definition of continuity, we can show that $\Lambda_{\vartheta}:=\Lambda$ also inherits the monotonicity and continuity of $\alpha(\cdot,\vartheta)$. Finally, we derive (2.1) from (2.6) by some estimate techniques.

3 Preliminary

This section is devoted to the introduce of some preliminary lemmas, which will be used in the next two sections.

Lemma 3.1**.**

Difference quotients and weak derivatives: Let $D$ be $\Omega$, or $\mathbb{T}$.

1. (1)

Suppose $1\leqslant p<\infty$ and $w\in W^{1,p}(D)$. Then $\|D^{h}_{\mathrm{h}}w\|_{L^{p}(D)}\lesssim\|\nabla_{\mathrm{h}}w\|_{L^{p}(D)}$. 2. (2)

Assume $1<p<\infty$, $w\in L^{p}(D)$, and there exists a constant $c$ such that $\|D^{h}_{\mathrm{h}}w\|_{L^{p}(D)}\leqslant c$. Then $\nabla_{\mathrm{h}}w\in L^{p}(D)$ satisfies $\|\nabla_{\mathrm{h}}w\|_{L^{p}(D)}\leqslant c$ and $D^{-h_{k}}_{\mathrm{h}}w\rightharpoonup\nabla_{\mathrm{h}}w$ in $L^{p}(D)$ for some subsequence $-h_{k}\to 0$.

Proof 1**.**

Following the argument of [5, Theorem 3], and use the periodicity of $w$, we can easily get the desired conclusions. $\Box$

Lemma 3.2**.**

Existence theory of a stratified (steady) Stokes problem (see [34, Theorem 3.1]): let $k\geqslant 0$, $f^{\mathrm{S},1}\in H^{k}$ and $f^{\mathrm{S},2}\in H^{k+1/2}$, then there exists a unique solution $(u,q)\in H^{k+2}\times\underline{H}^{k+1}$ satisfying

[TABLE]

Moreover,

[TABLE]

Lemma 3.3**.**

Equivalent form of instability condition: the instability condition (2.4) is equivalent to the following integral version of instability condition:

[TABLE]

Proof 2**.**

The conclusion in Lemma 3.3 is obvious, if we have the assertion:

[TABLE]

Next we verify (3.7) by two steps. Without loss of generality, we assume that $L_{1}^{2}={\max\{L_{1}^{2},L_{2}^{2}\}}$.

(1) We first prove that $a\geqslant L_{1}^{2}$. We choose a non-zero function $\psi\in H_{0}^{2}(-h_{-},h_{+})$ such that $\psi(0)\neq 0$. We denote

[TABLE]

then ${w}\in H^{1}_{\sigma,\vartheta}$ and

[TABLE]

which yields $a\geqslant L_{1}^{2}$.

(2) We turn to the proof of $a\leqslant L_{1}^{2}$. It should be noted that

[TABLE]

In fact, let $w\in H_{\sigma,3}^{1}$. Since $\mathrm{div}w=0$, we have

[TABLE]

Thus, using Pocare’s inequality, we have

[TABLE]

which immediately implies the assertion (3.8).

Let $w\in H_{\sigma,3}^{1}$. Then $|\nabla_{\mathrm{h}}w_{3}|_{0}^{2}\neq 0$. Let $\hat{w}_{3}(\xi,y_{3})$ be the horizontal Fourier transform of $w_{3}(y)$, i.e.,

[TABLE]

where $\xi=(\xi_{1},\xi_{2})$, then $\widehat{\partial_{3}w_{3}}=\partial_{3}\widehat{w}_{3}$. We denote $\psi(\xi,y_{3}):=\psi_{1}(\xi,y_{3})+\mathrm{i}\psi_{2}(\xi,y_{3}):=\hat{w}_{3}(\xi,y_{3})$, where $\psi_{1}$ and $\psi_{2}$ are real functions. Noting that $\psi(0)=0$, by Parseval theorem (see [7, Proposition 3.1.16]), we have

[TABLE]

which immediately yields that $a\leqslant L^{2}_{1}$. The proof is complete. $\Box$

Lemma 3.4**.**

Friedrichs’s inequality (see [24, Lemma 1.42]): Let $1\leqslant p<\infty$, and $D$ be a bounded Lipschitz domain. Let a set $\Gamma\subset\partial D$ be measurable with respect to the $(N-1)$-dimensional measure ${\mu}:=\mathrm{meas}_{N-1}$ defined on $\partial D$ and let $\mathrm{meas}_{N-1}(\Gamma)>0$. Then

[TABLE]

for all $u\in W^{1,p}(D)$ satisfying that the trace of $u$ on $\Gamma$ is equal to [math] a.e. with respect to the $(N-1)$-dimensional measure ${\mu}$.

Remark 3.1**.**

By Friedrichs’s inequality and the fact

[TABLE]

we get the Korn’s inequality

[TABLE]

Lemma 3.5**.**

Trace estimates:

[TABLE]

Proof 3**.**

See [22, Lemma 9.7] for (3.11). Since $C_{\sigma}^{\infty}:=C^{\infty}_{0}(\mathbb{R}^{2}\times(-h_{-},h_{+}))\cap H_{\sigma}^{1}$ is dense in $H^{1}_{\sigma}$, it suffices to prove that (3.12) holds for any $w\in C_{\sigma}^{\infty}$ by (3.11).

Let $\hat{w}$ be the horizontal Fourier transformed function of $w\in C_{\sigma}^{\infty}$, and

[TABLE]

Then

[TABLE]

and $\psi(\cdot,y_{3})\in H^{2}_{0}(-h_{-},h_{+})$, because of $\mathrm{div}{w}=0$ and $w|_{\Sigma_{-}^{+}}=0$. Moreover,

[TABLE]

In addition, we can deduce from (3.13) that

[TABLE]

By (3.14) and the Fubini and Parseval theorems, one has

[TABLE]

and

[TABLE]

where

[TABLE]

Using (3.13), we find that

[TABLE]

which imply that

[TABLE]

for given $\xi\in(L^{-1}_{1}\mathbb{Z}\times L^{-1}_{2}\mathbb{Z})\backslash\{0\}$. Employing (3.15)–(3.17) and the relation

[TABLE]

we obtains

[TABLE]

Similarly, we also have

[TABLE]

which, together with (3.18), yields the desired conclusion. This completes the proof. $\Box$

Remark 3.2**.**

From the derivation of (3.12), we easily see that

[TABLE]

Lemma 3.6**.**

Negative trace estimate:

[TABLE]

Proof 4**.**

Estimate (3.20) can be derived by integration by parts and an inverse trace theorem [24, Lemma 1.47]. $\Box$

Lemma 3.7**.**

Let $X$ be a given Banach space with dual $X^{*}$ and let $u$ and $w$ be two functions belonging to $L^{1}((a,b),X)$. Then the following two conditions are equivalent

1. (1)

For each test function $\phi\in C_{0}^{\infty}(a,b)$,

[TABLE] 2. (2)

For each $\eta\in X^{*}$,

[TABLE]

in the scalar distribution sense, on $(a,b)$, where $<\cdot,\cdot>_{X\times X^{*}}$ denotes the dual pair between $X$ and $X^{*}$.

Proof 5**.**

See Lemma 1.1 in Chapter 3 in [28].

4 Linear instability

In this section, we will use modified variational method to construct unstable solutions for the linearized RT problem. The modified variational method was firstly used by Guo and Tice to construct unstable solutions to a class of ordinary differential equations arising from a linearized RT instability problem [9]. In this paper, we directly apply Guo and Tice’s modified variational method to the partial differential equations (2.5), and thus obtain a linear instability result of the RT problem by further using an existence theory of stratified Stokes problem. Next we prove Theorem 2.1 by four subsections.

4.1 Existence of weak solutions to the modified problem

In this subsection, we consider the existence of weak solutions to the modified problem

[TABLE]

where $s>0$ is any given. To prove the existence of weak solutions of the above problem, we consider the variational problem of the functional $F(\varpi,s)$:

[TABLE]

for given $s>0$, where we have defined that

[TABLE]

Sometimes, we denote $\alpha(s,\vartheta)$ and $F(\varpi,s)$ by $\alpha$ (or $\alpha(s)$) and $F(\varpi)$ for simplicity, resp.. Then we have the following conclusions.

Proposition 4.1**.**

Let $s>0$ be any given.

In the variational problem (4.2), $F(\varpi)$ achieves its supremum on $\mathcal{A}$. 2. 2.

Let $w$ be a maximizer and $\alpha:=\sup_{\varpi\in\mathcal{A}}F(\varpi)$, the $w$ is a weak solution the boundary problem (4.1) with given $\alpha$.

Proof 6**.**

Noting that

[TABLE]

thus, by Young’s inequality and Korn’s inequality (3.10), we see that $\{F(\varpi)\}_{\varpi\in\mathcal{A}}$ has an upper bound for any $\varpi\in\mathcal{A}$. Hence there exists a maximizing sequence $\{w^{n}\}_{n=1}^{\infty}\subset\mathcal{A}$, which satisfies $\alpha=\lim_{n\to\infty}F(w_{n})$. Moreover, making use of (4.3), the fact $\|\sqrt{{\rho}}w^{n}\|_{0}=1$, trace estimate (3.12) and Young’s and Korn’s inequalities, we have $\|w^{n}\|_{1}+\vartheta|\nabla_{\mathrm{h}}w^{n}_{3}|_{0}\leqslant c_{1}$ for some constant $c_{1}$, which is independent of $n$. Thus, by the well-known Rellich–Kondrachov compactness theorem and (4.3), there exist a subsequence, still labeled by $w^{n}$, and a function $w\in\mathcal{A}$, such that

[TABLE]

Exploiting the above convergence results, and the lower semicontinuity of weak convergence, we have

[TABLE]

Hence $w$ is a maximum point of the functional $F(\varpi)$ with respect to $\varpi\in\mathcal{A}$.

Obviously, $w$ constructed above is also a maximum point of the functional $F(\varpi)/\|\sqrt{\rho}\varpi\|^{2}_{0}$ with respect to $\varpi\in H_{\sigma,\vartheta}^{1}$. Moreover $\alpha=F(w)/\|\sqrt{\rho}w\|^{2}_{0}$. Thus, for any given $\varphi\in H^{1}_{\sigma,\vartheta}$, the point $t=0$ is the maximum point of the function

[TABLE]

Then, by computing out $I^{\prime}(0)=0$, we have the weak form:

[TABLE]

Noting that (4.4) is equivalent to

[TABLE]

The means that $w$ is a weak solution of the modified problem (4.1). $\Box$

4.2 Improving the regularity of weak solution

By Proposition 4.1, the boundary-value problem (4.1) admits a weak solution $w\in H^{1}_{\sigma,\vartheta}$. Next we further improve the regularity of $w$.

Proposition 4.2**.**

Let $w$ be a weak solution of the boundary-value problem (4.1). Then $w\in H^{\infty}$.

Proof 7**.**

To begin with, we shall establish the following preliminary conclusion:

For any $i\geqslant 0$, we have**

[TABLE]

and**

[TABLE]

Obviously, by induction, the above assertion reduces to verify the following recurrence relation:

For given $i\geqslant 0$, if $w\in H^{1,i}_{\sigma,\vartheta}$ satisfies (4.6) for any $\varphi\in H^{1}_{\sigma,\vartheta}$, then**

[TABLE]

and $w$ satisfies**

[TABLE]

Next we verify the above recurrence relation by method of difference quotients.

Now we assume that $w\in H^{1,i}_{\sigma,\vartheta}$ satisfies (4.6) for any $\varphi\in H^{1}_{\sigma,\vartheta}$. Noting that $\partial_{\mathrm{h}}^{i}w\in H^{1}_{\sigma,\vartheta}$, we can deduce from (4.6) that, for $j=1$ and $2$,

[TABLE]

and

[TABLE]

which yield that

[TABLE]

and

[TABLE]

resp..

By Korn’s inequality,

[TABLE]

thus, using (4.3), Young’s inequality, and the first conclusion in Lemma 3.1 , we further deduce from (4.10) that

[TABLE]

Thus, using (4.3), trace estimate (3.12) and the second conclusion in Lemma 3.1, there exists a subsequence of $\{-h\}_{h\in\mathbb{R}}$ (still denoted by $-h$) such that

[TABLE]

Using regularity of $w$ in (4.11) and the fact $w\in H_{\sigma,\vartheta}^{1,i}$, we have (4.7). In addition, exploiting the limit results in (4.11), we can deduce (4.8) from (4.9). This complete the proof of the recurrence relation, and thus (4.5) holds.

With (4.5) in hand, we can consider a stratified Stokes problem:

[TABLE]

where $k\geqslant 0$ is a given integer, and we have defined that

[TABLE]

Recalling the regularity (4.5) of $w$, we see that $\partial_{\mathrm{h}}^{k}w\in L^{2}$, and $\partial_{\mathrm{h}}^{k}\mathcal{L}^{1}\in H^{1}(\mathbb{T})$. Applying the existence theory of stratified Stokes problem (see Lemma 3.2), there exists a unique strong solution $(\omega^{k},\beta^{k})\in H^{2}\times\underline{H}^{1}$ of the above problem (4.12).

Multiplying (4.12)1 by $\varphi\in H^{1}_{\sigma,\vartheta}$ in $L^{2}$ (i.e., taking the inner product in $L^{2}$), and using the integration by parts and (4.12)2–(4.12)4, we have

[TABLE]

Subtracting the two identities (4.6) and (4.13) yields that

[TABLE]

Taking $\varphi:=\partial_{\mathrm{h}}^{k}w-\omega^{k}\in H^{1}_{\sigma,\vartheta}$ in the above identity, and using the Korn’s inequality, we find that $\omega^{k}=\partial_{\mathrm{h}}^{k}w$. Thus we immediately see that

[TABLE]

which implies $\partial_{\mathrm{h}}^{k}w\in H^{1}$, and $\partial_{\mathrm{h}}^{k}\mathcal{L}^{2}\in H^{2}(\mathbb{T})$ for any $k\geqslant 0$. Thus, applying the stratified Stokes estimate (3.5) to (4.12), we have

[TABLE]

Obviously, by induction, we can easily follow the improving regularity method from (4.14) to (4.15) to deduce that $w\in H^{\infty}$. In addition, we have $\beta:=\beta^{0}\in H^{\infty}$; moreover, $\beta^{k}$ in (4.12) is equal to $\partial^{k}_{\mathrm{h}}\beta$.

Finally, recalling the embedding $H^{k+2}\hookrightarrow C^{0}(\overline{\Omega})$ for any $k\geqslant 0$, we easily see that $(w,\beta)$ constructed above is indeed a classical solution to the modified problem (4.1).

4.3 Some properties of the function $\alpha(s)$

In this subsection, we shall derive some properties of the function $\alpha(s)$, which make sure the existence of fixed point of $\sqrt{\alpha(s)}$ in $\mathbb{R}^{+}$.

Proposition 4.3**.**

For given $\vartheta\in\mathbb{R}^{+}_{0}$, we have

[TABLE]

Proof 8**.**

To being with, we verify (4.16). For given $s_{2}>s_{1}$, then there exist $v^{s_{2}}\in\mathcal{A}$ such that $\alpha(s_{2})=F(v^{s_{2}},s_{2})$. Thus, by Korn’s inequality and the fact $\|\sqrt{\rho}v^{s_{2}}\|_{0}=1$,

[TABLE]

which yields (4.16).

Now we turn to prove (4.17). Choosing a bounded interval $[c_{4},c_{5}]\subset(0,\infty)$, then, for any $s\in[c_{4},c_{5}]$, there exists a function $v^{s}$ satisfying $\alpha(s)=F(v^{s},s)$. Thus, by the monotonicity (4.16), we have

[TABLE]

which yields

[TABLE]

Thus, for any $s_{1}$, $s_{2}\in[c_{4},c_{5}]$,

[TABLE]

and

[TABLE]

which immediately imply $|\alpha(s_{1})-\alpha(s_{2})|\leqslant\xi|{s_{2}}-s_{1}|$. Hence (4.17) holds.

Finally, (4.18) can be deduced from the definition of $\alpha$ by using Korn’s inequality and (4.3), while (4.19) is obvious by the definition of $\alpha$, Lemma 3.3 and (3.6). $\Box$

4.4 Construction of an interval for fixed point

Let $\mathfrak{I}:=\sup\{$all the real constant $s$, which satisfy that $\alpha(\tau)>0$ for any $\tau\in(0,s)\}$. By virtue of (4.18) and (4.19), $\mathfrak{I}\in\mathbb{R}_{+}$. Moreover, $\alpha(s)>0$ for any $s\in(0,\mathfrak{I})$, and, by the continuity of $\alpha(s)$,

[TABLE]

Using the monotonicity and the upper boundedness of $\alpha(s)$, we see that

[TABLE]

Now, exploiting (4.20), (4.21) and the continuity of ${\alpha}(s)$ on $(0,\mathfrak{I})$, we find by a fixed-point argument on $(0,\mathfrak{I})$ that there is a unique $\Lambda\in(0,\mathfrak{I})$ satisfying

[TABLE]

Thus we get a classical solution $(w,\beta)\in H^{\infty}$ to the boundary problem (2.5) with $\Lambda$ constructed by (4.22). Moreover,

[TABLE]

In addition, (2.7) directly follows (4.23) and the fact $w\in H_{\sigma}^{1}$.

4.5 Largest growth rate

Next we shall prove that $\Lambda$ constructed in previous section is the largest growth rate of RT instability in the linearized RT problem, and thus complete the proof of Theorem 2.1.

Proposition 4.4**.**

Under the assumptions of Theorem 2.1, $\Lambda>0$ constructed by (4.22) is the largest growth rate of RT instability in the linearized RT problem.

Proof 9**.**

Recalling the definition of largest growth rate, it suffices to prove that $\Lambda$ enjoys the first condition in Definition 2.1.

Let $u$ be strong solution to the linearized RT problem. Then we derive that, for a.e. $t\in I_{T}$ and all $w\in H_{\sigma}^{1}$,

[TABLE]

thus,

[TABLE]

Using regularity of $(\eta,u)$, we can show that the right hand side of (4.25) is bounded above by $A(t)(\|w\|_{1}+|w|_{1})$ for some positive function $A(t)\in L^{2}(I_{T})$. Then there exists a $f\in L^{2}(I_{T},H^{-1}_{\sigma})$ such that, for a.e. $t\in I_{T}$,

[TABLE]

Hence it follows from Lemma 3.7 that

[TABLE]

In addition, by a classical regularization method (referring to Theorem 3 in Chapter 5.9 in [5] and Lemma 6.5 in [24]), we have

[TABLE]

Therefore, we can derive from (4.26) and the above two identities that

[TABLE]

Then, integrating the above identity in time from [math] to $t$ yields that

[TABLE]

Using Newton–Leibniz’s formula and Young’s inequality, we find that

[TABLE]

In addition, by (2.6), we have

[TABLE]

Thus, we infer from (4.27)–(4.29) that

[TABLE]

Recalling that

[TABLE]

we further deduce from (4.30) the differential inequality:

[TABLE]

Applying Gronwall’s inequality [24, Lemma 1.2] to the above inequality, one concludes

[TABLE]

which, together with (4.30), yields

[TABLE]

Multiplying (1.10)2 by $u_{t}$ in $L^{2}$ and using the integral by parts, we get

[TABLE]

Exploiting (3.20), we can estimate that

[TABLE]

In addition, using (1.10)5 and trace estimate (3.12), we have

[TABLE]

Using the above two estimates, we can derive from (4.33) that

[TABLE]

which implies that

[TABLE]

By the above estimate and Korn’s inequality, we derive from (4.31) and (4.32) that

[TABLE]

Finally, from (1.10)1 we get

[TABLE]

By the two estimates above, we see that $\Lambda$ satisfies the first condition in Definition 2.1. The proof is complete. $\Box$

5 Effect of surface tension

5.1 Properties of $\alpha(s,\vartheta)$ with respect to $\vartheta$

To emphasize the dependence of $\Lambda$ and $\mathcal{G}$ upon $\vartheta$, we will denote them by $\Lambda_{\vartheta}$ and $\mathcal{G}_{\vartheta}$, respectively. To prove Theorem 2.2, we shall further derive the relations (2.1) and (2.11) of surface tension coefficient and the largest growth rate. To this end, we need the following auxiliary conclusions:

Proposition 5.1**.**

Let $g>0$, $\rho>0$ and $\mu>0$ are given.

Strict monotonicity: if $\vartheta_{1}$ and $\vartheta_{2}$ are constants satisfying $0\leqslant\vartheta_{1}<\vartheta_{2}$, then

[TABLE]

for any given $s>0$. Moreover, if $\vartheta_{2}$ further satisfies $\vartheta_{2}<\vartheta_{\mathrm{c}}$,

[TABLE]

where

[TABLE] 2. 2.

Continuity: for given $s>0$, $\alpha(s,\vartheta)\in C^{0,1}_{\mathrm{loc}}(\mathbb{R}_{+})$ with respect to the variable $\vartheta$.

Proof 10**.**

(1) Let $s>0$ be fixed, and $0\leqslant\vartheta_{1}<\vartheta_{2}$. Then there exist functions ${w}^{\vartheta_{i}}\in H^{\infty}\cap\mathcal{A}_{\vartheta_{i}}$, $i=1,2$, such that

[TABLE]

where ${E}({w}^{\vartheta_{i}}):=g\llbracket\rho\rrbracket|{w}^{\vartheta_{i}}_{3}|_{0}^{2}-{s}\|\sqrt{\mu}\mathbb{D}{w}^{\vartheta_{i}}\|_{0}^{2}/2$. Since ${w}^{\vartheta_{i}}\in\mathcal{A}_{\vartheta_{i}}$, by virtue of (2.7) and (3.8), we have

[TABLE]

and, thus

[TABLE]

This yields the desired conclusion (5.1).

Next we prove (5.2) by contradiction. If $\mathcal{G}_{\vartheta_{\mathrm{1}}}<\mathcal{G}_{\vartheta_{\mathrm{2}}}$, then we get from (5.1) and the strict monotonicity of $\alpha(s,\cdot)$ with respect to $s$ that

[TABLE]

which is a paradox. If $\mathcal{G}_{\vartheta_{\mathrm{1}}}=\mathcal{G}_{\vartheta_{\mathrm{2}}}$, exploiting (5.1), we have

[TABLE]

which is also a paradox. Thus we immediately get the desired conclusion.

(2) Let $s>0$ be fixed. We choose a bounded interval $[b_{1},b_{2}]\subset\mathbb{R}_{+}$. Then, for any given $\theta\in[b_{1}/2,b_{2}]$, there is a function ${w}^{\theta}\in\mathcal{A}_{\vartheta}$ satisfying $\alpha(s,\theta)={E}({w}^{\theta})-\theta|\nabla_{\mathrm{h}}{w}_{3}^{\theta}|_{0}^{2}$. Thus, in view of the monotonicity of $\alpha(\cdot,\theta)$, we know that

[TABLE]

which yields

[TABLE]

Thus, for any $\vartheta_{1}$, $\vartheta_{2}\in[b_{1},b_{2}]$,

[TABLE]

Reversing the role of the indices $1$ and $2$ in the derivation of the above inequality, we obtain the same boundedness with the indices switched. Therefore, we deduce that

[TABLE]

which yields $\alpha(s,\vartheta)\in C^{0,1}_{\mathrm{loc}}(\mathbb{R}_{+})$. This completes the proof. $\Box$

5.2 Proof of Theorem 2.2

First, we verify the monotonicity of $\Lambda_{\vartheta}$ with respect to the variable $\vartheta\in[0,\vartheta_{\mathrm{c}})$.

For given two constants $\vartheta_{1}$ and $\vartheta_{2}$ satisfying $0\leqslant\vartheta_{1}<\vartheta_{2}<\vartheta_{\mathrm{c}}$, then there exist two associated curve functions $\alpha(s,\vartheta_{1})$ and $\alpha(s,\vartheta_{2})$ defined in $(0,\vartheta_{\mathrm{c}})$. By the first assertion in Proposition 5.1.

[TABLE]

On the one hand, the fixed-point $\Lambda_{\vartheta_{i}}$ satisfying $\Lambda_{\vartheta_{i}}=\sqrt{\alpha(\Lambda_{\vartheta_{i}})}$ can be obtained from the intersection point of the two curves $y=s$ and $y=\sqrt{\alpha(s,\vartheta_{i})}$ on $(0,\mathcal{G}_{\vartheta_{i}})$ for $i=1$ and $2$. Thus we can immediately observe the monotonicity

[TABLE]

*Second, we prove the continuity for $\Lambda_{\vartheta}$. *

We choose a constant $\vartheta_{0}>0$ and an associated function $\alpha(s,\vartheta_{0})$. Noting that $\alpha(\Lambda_{\vartheta_{0}},\vartheta_{0})=\Lambda^{2}_{\vartheta_{0}}>0$ and $\alpha(\cdot,\vartheta)\in C^{0,1}_{\mathrm{loc}}[0,\vartheta_{\mathrm{c}})$ are strictly decreasing and continuous with respect to $\vartheta$, then, for any given $\varepsilon>0$, there exists a constant $\delta>0$, such that

[TABLE]

and

[TABLE]

In particular, we have

[TABLE]

By the monotonicity of $\Lambda_{\vartheta}$ with respect to $\vartheta$, we get

[TABLE]

Thus, using the monotonicity of $\alpha(s,\cdot)$ with respect to $s$, we obtain

[TABLE]

and

[TABLE]

Chaining the five inequalities above, we immediately get

[TABLE]

Then, for any $\vartheta\in(\vartheta_{0}-\delta,\vartheta_{0}+\delta)$, we arrive at $\Lambda_{\vartheta_{0}}-\varepsilon<\Lambda_{\vartheta}<\Lambda_{\vartheta_{0}}+\varepsilon$. Hence

[TABLE]

Now we study the limit of $\Lambda_{\vartheta}$ as $\vartheta\to 0$. For any $\varepsilon>0$, there exits a ${w}\in\mathcal{A}_{0}$ such that

[TABLE]

In addition,

[TABLE]

Thus, making use of (2.6), (5.7) and (5.8), there exists a sufficiently small constant $\vartheta_{1}\in(0,\vartheta_{\mathrm{c}})$ such that, for any $\vartheta\in(0,\vartheta_{1})$,

[TABLE]

Hence we get

[TABLE]

which, together with (5.6), yields that

[TABLE]

Finally, we derive the upper bound (2.1) for $\Lambda_{\vartheta}$.

Recalling the definition of $\vartheta_{\mathrm{c}}$, we see from (3.7) that

[TABLE]

Hence, by virtue of (2.6), for any given $\vartheta\in[0,\vartheta_{\mathrm{c}})$, there exists a ${w}^{\vartheta}\in\mathcal{A}_{\vartheta}$ such that

[TABLE]

which yields that

[TABLE]

By (3.9) and trace estimate (3.12), we can estimate that

[TABLE]

Similarly, we also have

[TABLE]

By the above two estimates, we derive from (5.11) that

[TABLE]

which yields that

[TABLE]

Noting that $\|\sqrt{\rho}{w}^{\vartheta}\|_{0}=1$, then, by (3.19),

[TABLE]

Putting the above estimate into (5.11), and then using Young’s inequality, we get

[TABLE]

which yields that

[TABLE]

Similarly, we also have

[TABLE]

Summing up the above two estimates yields that

[TABLE]

which, together with (5.12), implies that

[TABLE]

Consequently we complete the proof of Theorem 2.2 from (5.5), (5.10) and (5.16).

Acknowledgements. The research of Fei Jiang was supported by NSFC (Grant No. 11671086).