note, stagnation points beneath the surface – points within the bulk of the fluid that travels with the same speed
as the wave. Such configurations can be also observed in waves with a point vortex [22, 38, 8, 48]. Another
peculiarity of flows with constant vorticity is the emergence of pressure anomalies, such as the maxima and
minima of pressure on the bottom boundary being attained at locations other than beneath the crests and
troughs of the free surface respectively.
In some practical situations is of interest to couple the fluid motion with electric fields. This subject of study
is known as Electrohydrodynamics (EHD) and it can be applied widely in chemical engineering such as cooling
systems in conducting pumps [16], coating process [21] and etc. A good review paper was recently published
by Papageorgiou [31]. Beginning with the early work and experiments of Taylor and McEwan [41], researchers
have been interested in the effect of electric fields on fluids with varying electrical properties. In that particular
paper, normal electric fields were shown to be capable to destabilise the interface between two fluids. It was
quickly followed by work on tangential electric fields by Melcher and Schwarz [30], where, on the contrary, the
interfacial waves can be stabilised due to such effects. Kelvin-Helmholtz and Rayleigh-Taylor instabilities were
shown to be controlled and suppressed by tangential electric fields in [51] and [1, 3] respectively. Besides many
works were achieved theoretically by multi-scale technique, e.g. [12, 23, 20, 24, 35, 32, 27] in which weakly
nonlinear model equations were derived for the two-dimensional problem. The readers are referred to the work
by Wang [50] for a comprehensive review. Travelling waves were computed by a boundary integral method, e.g.
[33, 34, 10, 9], for different configurations. Of note, in [27, 9], the electric fields were shown to be capable of
changing the type of the associated nonlinear Schrödinger equation from focusing to defocusing such that dark
solitary waves bifurcate at the phase speed minimum. Nonlinear wave interactions in this context were studied
theoretically and numerically by Gao et al. [19]. Unsteady simulations were performed based on the time-
dependent conformal mapping technique in [14, 18] for a special case of capillary-gravity waves on a dielectric
fluid of finite and infinite depth respectively covered above by a conducting gas layer. More recently, Gao et al.
[17] followed to investigate the singularity formation of the capillary-gravity wave due to the effect of normal
electric fields.
In this work, we consider a dielectric fluid in a two-dimensional space bounded below by an electrode, and
above by a conducting passive gas, which in turn is bounded above by another electrode. A potential difference
is applied between the electrode such that the electric field is imposed in the direction perpendicular to the
undisturbed fluid surface. This configuration is well known for destabilising the interface between the fluids. The
fluid is assumed to be under the effect of constant vorticity. Our aim is to study numerically the flow structure
beneath waves with constant vorticity under the influence of normal electric fields. More specifically, we make
use of the conformal mapping technique and pseudo-spectral numerical methods to compute the locations of
stagnation points and the pressure. The rest of this paper is structured as follows. The mathematical formulation
is given in section 2. The linear theory is studied in section 3. The numerical scheme is presented in section 4.
Section 5 is devoted to the results on the full nonlinear governing equations. Finally, concluding remarks are
given in section 6.
2 Mathematical Formulation
An inviscid and incompressible dielectric fluid with density ρof mean depth hand electric permittivity 1is
considered in a two-dimensional space. It is bounded below by a wall electrode and surrounded by a region
occupied by perfectly conducting gas and enclosed by another wall electrode atop. In the electrostatic limit
of Maxwell’s equations, the induced magnetic fields are negligible so the electric field is irrotational due to
Faraday’s law. A potential function Vof the electric fields ~
Ecan be then introduced such that ~
E=∇V. It
follows that Vsatisfies the Laplace Equation in the dielectric fluid layer. The potential difference between the
two electrodes is denoted by V0. The voltage potential is invariant in the conducting gas layer. Without losing
generality, we set V=−V0on the bottom boundary and V= 0 in the gas layer. A schematic of the problem is
presented in Figure 1 . The gravitational acceleration and the surface tension coefficient are denoted by gand
T, respectively. A Cartesian x−ycoordinate system is introduced with the gravity pointing in the negative
y-direction and y= 0 being the undisturbed free-surface. It follows that the bottom boundary of the fluid is
at y=−h, and the upper boundary is free to move and denoted by η(x, t). We consider a periodic wavetrain
propagating at a constant speed in the positive x-direction. Also, we let x= 0 be at a wave crest such that the
wave profile is symmetric by the y-axis.
The velocity field is assumed to be an irrotational perturbation of a shear flow, namely
~
U= ˆu+~u, (1)
where the vorticity is constant that equals γ, and ˆu= (−γ(y+b),0) with bbeing another constant. The velocity
field ~u is irrotational so that there exists a potential function φsatisfying ~u =∇φ. We seek steady waves of
wavelength λin a frame of reference moving with the wave speed c. This is achieved by a change of variables:
2