Abstract

We are concerned with gravity-capillary waves propagating on the surface of a three-dimensional electrified liquid sheet under a uniform electric field parallel to the undisturbed free surface. For simplicity, we make an assumption that the permittivity of the fluid is much larger than that of the upper-layer gas; hence, this two-layer problem is reduced to be a one-layer problem. In this paper, we propose model equations in the shallow-water regime based on the analysis of the Dirichlet-Neumann operator. The modified Benney-Luke equation and Kadomtsev-Petviashvili equation will be derived, and the truly three-dimensional fully localized traveling waves, which are known as “lumps” in the literature, are numerically computed in the Benney-Luke equation.

1. Introduction

Interfacial electrohydrodynamic waves have many important applications in mechanical, chemical, and electrical industries, such as electrospray ionization, cooling systems, coating process, and electrowetting (see [110] and the references therein). The electric stresses have great impact on electrified free surfaces or interfaces, which can not only modify wave patterns but also change stability characteristics of the system. These influences are of great practical interest since they can lead to quantity (e.g., heat or mass) transfer enhancement (see, e.g., [5, 11, 12]). The study of the stability of interfacial electrohydrodynamic waves was initiated by Melcher [3] and Taylor and McEwan [13], and the role of interfacial stresses resulting from electrodes was reviewed by Melcher and Taylor [4]. Recent theoretical research in this field has focused on the nonlinear phenomena and corresponding mechanisms, such as nonlinear coherent structures (e.g., [2, 5, 8, 10, 11, 1420]) and touchdown singularities (e.g., [7, 21]). Considerable effects have been put into the modeling and numerical studies of nonlinear interfacial electrohydrodynamic waves.

For incompressible, inviscid, and irrotational flows, electrified interfacial waves were studied in the nonlinear regime under externally applied electric fields by different groups. The electric field is usually posed parallel or perpendicular to the undisturbed interface and is therefore called tangential or normal electric field, respectively. The orientation of the electric field with respect to the undisturbed interface plays an important role in the system. From the point of view of the linear theory, the dispersion relations are of the forms for two-dimensional gravity-capillary waves in deep water under normal and tangential electric field, respectively, where is the frequency, the wavenumber, the acceleration due to gravity, the surface tension coefficient, the strength of the electric field, the density of the liquid, and the ratio of permittivities. The forms of functions and depend on the specific problem; for example, for normal electric field acting on a perfect conducting fluid ([2, 8]) and for a dielectric fluid under a tangential electric field ([10]). Linearly speaking, a tangential electric field provides a dispersive contribution so as to stabilize the system, while a considerably strong normal electric field can provide energy to a certain range of wavenumbers to induce instability. On the nonlinear side, Tilley et al. [7] proposed a nonlinear long-wave system to study the dynamics of electrified liquid sheets and found that the tangential electric field can delay the formation of the film rupture. More than that, Barannyk et al. showed in [15] that the tangential electric field can even suppress the Rayleigh-Taylor instability in some situations. On the contrary, it can be deduced from the work by Gleeson et al. [16], Papageorgiou et al. [21], Lin et al. [2], and Wang [8] that the normal electric field has a destabilizing effect on the interface.

N. M. Zubarev and O. V. Zubareva [10, 20] and Tao and Guo [22] considered the electrified gas-fluid or vacuum-fluid interface so that they can make the assumption that the permittivity of the fluid was much larger compared to that of gas (the permittivities for pure water and air are 80 and 1, resp.). As a consequence, the actual two-layer problem could be reduced to one layer, and the theoretical and numerical techniques developed for free-surface water wave problems can be generalized to include the electric field.

On the numerical side, most researches focus on periodic waves (see, e.g., [5, 7, 10, 14, 21]); however the studies of electrified solitary waves have begun recently. Easwaran [23] first derived the Korteweg-de Vries (KdV) equation in the context of electrohydrodynamic waves which admits soliton solutions naturally. Hammerton [17] studied solitary waves in the KdV-Benjamin-Ono equation, a model arising from a conducting liquid sheet under a normal electric field generated by two electrodes with a sufficiently large separation distance. Barannyk et al. [15] computed interfacial solitary waves based on a system describing long waves between two immiscible dielectric fluids. It is noted that KdV-Benjamin-Ono equation was first derived in the context of interfacial waves between two immiscible fluids when surface tension is sufficiently strong (e.g., [24, 25]).

All the aforementioned papers confine themselves to two-dimensional (2D) problems, but there are relatively fewer studies on the three-dimensional (3D) case. Recent study on 3D electrohydrodynamic modeling has focused on conducting fluids under normal electric fields: Hunt et al. [18] proposed a one-way model with weak transverse variations which is called the Benjamin-Ono-Kadomtsev-Petviashvili equation, Aliev and Yurchenko [14] established a reduced system for the same setup, and Wang [8] developed fully nonlinear numerical models and weakly nonlinear theories for electrohydrodynamic surface waves in the Hamiltonian framework based on analyses of the Dirichlet-Neumann operators. In this paper, we are concerned with a dielectric liquid under a tangential electric field, and restoring forces due to gravity, surface tension, and Maxwell stresses are all taken into account. We will propose weakly nonlinear models (including both unidirectional and bidirectional equations) under Zubarev’s approximation and numerically search for truly 3D solitary waves (which is known as “lumps” in the literature) within these models.

The rest of the paper is structured as follows. In Section 2, the governing equations are described and simplified by making the assumption that the permittivity of the fluid is large. In Section 3, the Dirichlet-Neumann operator is introduced, which plays an important role in deriving weakly nonlinear models. In Section 4, a bidirectional model (the modified Benney-Luke equation) and a unidirectional model (the modified Kadomtsev-Petviashvili equation) are derived under the long-wave approximation. Then the main numerical results are followed, including the typical profiles of solitary waves, the bifurcation diagrams, and the formation of lumps resulting from the transverse instability of plane solitary waves.

2. Formulation

2.1. Governing Equations

We consider an incompressible, inviscid, and irrotational fluid, moving on a solid wall, under a horizontal electric field with uniform strength in the far field. We introduce the Cartesian coordinates, such that and are horizontal variables, and points upwards with being the undisturbed air-water interface (see Figure 1). We denote the displacement of the air-water interface by ; then the velocity potential of the fluid, which is denoted by , is governed by And the voltage potentials within and above the fluid domain, denoted by and , respectively, are also governed by Laplace’s equation: The kinematic and dynamic boundary conditions at the interface are given by (see, e.g., [10, 22, 26]) where , , and is the ratio of the permittivity of water, , to that of air, . The detailed derivation of the electrostatic pressure at the interface between an ideal dielectric liquid and air (or vacuum) in the presence of free electric charges can be found in [26] or [5]. At the interface , the continuity condition of voltage potentials and the normal stress induced by electricity are, respectively, Following Barannyk et al. [15], the kinematic boundary conditions for and at are of the same form Finally, the far-field conditions for the electric field, complete the whole system.

2.2. Simplification

In order to simplify the problem, we make assumption that the quantity is small. This is a reasonable assumption in many physical situations when the upper layer is occupied by air () and the lower layer is a fluid with much lager electric permittivity such as water (). This assumption could decouple the system. We then nondimensionalize the system by choosing as the time, length, voltage potential, and velocity potential scales, respectively. Therefore, following [10, 19, 20] the dynamic boundary condition can be rewritten as where , and the rescaled surface tension coefficient becomes the Bond number . The rigid bottom is now represented by , and as . Meanwhile, the normal stress due to the electric field becomes

Finally, we introduce the variable by , and, as a consequence, as and at . It also follows at the interface that

3. Dirichlet-Neumann Operators

The Dirichlet-Neumann operator (DNO), which maps the Dirichlet boundary condition to normal derivatives on the boundary via solving Laplace’s equation, plays an essential role in surface water wave problems. Coifman and Meyer [27] initially proved that if the Lipschitz-norm of wave profiles is smaller than a certain constant then DNO is analytic in two-dimensional problems. Subsequently, the analyticity of this operator was extended to three dimensions by Craig et al. [28], and Craig and Nicholls [29] and Nicholls and Reitich [30] delivered further results. If we define as the surface velocity potential, the DNO, denoted by , can be defined by In the same vein, we define the surface voltage . Since also satisfies the Laplace equation as well as the boundary condition at , as a consequence, we can use the same Dirichlet-Neumann operator to define Using the boundary condition of the voltage potential (12), we have If we define a new pseudo-differential operator , it follows that . Here we remark that strictly speaking is not a well-defined operator due to the singularity at the origin in the Fourier space; however, in the present problem always appears with (see (17) and (19)) which removes the singularity. We then reformulate the problem by use of the DNO. It is not difficult to conduct the following calculation:Similarly, we can haveUsing the above formulae, the kinematic and dynamic boundary conditions can be rewritten as follows:Linearizing the initial system (18)-(19) and seeking solutions in the form lead to the dispersion relation where is the wave frequency and is the wavenumber with its component corresponding to - and -directions, respectively.

4. Long-Wave Models

Deriving long-wave models using the DNO was initiated by Craig and Groves [31]. It is well-known that the DNO is an analytical operator if the -norm of is smaller than a constant. Therefore, the Taylor series of the DNO converge; namely, can be expressed as , and the first three terms are given by It is worth mentioning that, in contrast to [32], we expand the DNO in this paper based on the results of Nicholls and Reitich [33], who provided the joint analyticity and analytic continuation results for the DNO. We assume is the typical length scale in -direction. Since the bottom of the fluid is at , we define . If we consider the isotropic problem; namely, the typical wavelength in -direction, , is similar to ; therefore, and we can expandFollowing [34], we can also expend as In the meantime, we obtain

4.1. Benney-Luke Scaling

The isotropic Boussinesq regime includes the additional expectation of small amplitude motions. First of all we assume and the Bond number . The scales for the variables are as follows: , ; the time scale is of order , and . It follows that We proceed by expanding (18)-(19) in terms of the small parameter and truncating at Taking derivative with respect to time for (27) yieldsSubstituting (26) into (28) and retaining terms valid to yieldEquation (27) implies that ; therefore, substituting (29) yieldsIf the original problem is independent of the transverse direction, (27) reduces to And its traveling-wave solutions satisfy which admits KdV soliton

4.2. Kadomtsev-Petviashvili Scaling

If the variation in the -direction is much slower than that in the -direction, we choose and the scaling of the other variables is the same as the case of Benney-Luke. In the same vein, it is not difficult to obtain the following two equations: Following the same argument stated in Section 4.1, we eliminate and arrive atWe change variables as Therefore, we can recast (35) as Truncating at and returning back to the original scale, we arrive at Taking derivative with respect to and letting , the KP equation takes the form

4.3. Lumps

Fully confined three-dimensional traveling waves, known as lumps, do not exist in pure gravity waves ([35]); however, if the surface tension is also included and strong enough, lumps do exist and were first predicted by the celebrated KP equation. Equation (39) is the KP equation whose coefficient of the term is modified by the electric field, and it admits analytical fully confined traveling-wave solutions, which is called the “lumps” in the literature. Pego and Quintero [36] mathematically proved the existence of lumps in the Benney-Luke equation (30). In the present paper, we numerically search for these solutions in (30). We assume translates with speed in the -direction, namely, ; therefore, (30) becomesThe numerical scheme adopted to find these lump-type solutions was first proposed by Petviashvili [37], and here we follow the procedure given by Wang [8]. Applying the Fourier transform to (40) leads to where is the vector of wavenumber. Following [8], we propose the numerical iteration scheme as where is defined asIn Figure 2 we present the bifurcation diagram and typical profile of lumps in (30) for . Figure 2(a) shows the speed-amplitude curve for fixed , and Figure 2(b) shows the -amplitude curve for fixed , both of which are monotonic. For positive , only depression lumps which feature a negative value at the center (see Figure 2(c)) were found.

For time-dependent simulation of the Benney-Luke equation (30), we introduce a new variable ; then it follows that which is a first-order dispersive system and can be integrated using pseudo-spectral method (see, e.g., [8, 15, 38, 39]). We study the transverse instability of plane solitary-wave solutions (33), which gives mechanics of generation of lumps. We denote a plane solution (33) by and perturb it using a long wave in the transverse direction; namely,is taken as the initial profile. The subsequent evolution of this initial data shows a focusing behaviour which is arrested by the formation of a lump propagating almost steadily (the snapshots are shown in Figure 3).

5. Conclusion

In the present paper, we have studied the electro-capillary-gravity waves with a horizontal electric field. For simplicity, we make assumption that the permittivity of the fluid is much larger than that of the upper-layer gas. Therefore, as the first-order approximation, the whole system could be reduced to a one-layer problem. The Dirichlet-Neumann operator, which transfers the Dirichlet boundary condition to the Neumann derivative on the boundary via solving the Laplace equation, plays an important role in this problem. Based on the long-wave approximation and symbolical calculations, we have derived the unidirectional KP-like model and bidirectional Benney-Luke-like model. Then the bifurcation diagram and typical profiles for three-dimensional solitary waves (i.e., lumps) were found in the modified Benney-Luke equation based on the Petviashvili method.

Very recently considerable attention has been paid to asymmetric gravity-capillary waves. Wang et al. [40] first found asymmetric plane solitary waves in the full Euler equations in [40] and asymmetric lumps in a model equation in [41]. Shimizu and Shoji [42] and Gao et al. [43] thoroughly investigated asymmetric periodic waves. Gao et al. [44] found asymmetric generalized solitary waves in the full Euler equations. These findings raise the question as to the existence and stability of asymmetric waves in electrified waves, which will be left for future studies.

Finally we remark that if the electric field is strong enough so that gravity and surface tension are negligible, linearizing the original system yields Combining the above equations, we end up with a wave equation for arbitrary depth of fluid, which is commonly true in shallow water only in water-wave problems.

Conflicts of Interest

The author declares that there are no conflicts of interest regarding the publication of this paper.