Flow structure beneath periodic waves with constant vorticity under normal electric fields M. V. Flamarion1 T. Gao2 R. Ribeiro-Jr3 A. Doak4

2025-05-06 0 0 4.28MB 14 页 10玖币
侵权投诉
Flow structure beneath periodic waves with constant vorticity
under normal electric fields
M. V. Flamarion1, T. Gao2, R. Ribeiro-Jr3& A. Doak4
1Unidade Acadêmica do Cabo de Santo Agostinho, UFRPE/Rural Federal University of Pernambuco, BR 101 Sul, Cabo de Santo
Agostinho-PE, Brazil, 54503-900
marcelo.flamarion@ufrpe.br
2School of Computing and Mathematical Sciences, University of Greenwich, London SE10 9LS, UK.
t.gao@gre.ac.uk
3UFPR/Federal University of Paraná, Departamento de Matemática, Centro Politécnico, Jardim das Américas, Caixa Postal
19081, Curitiba, PR, 81531-980, Brazil
robertoribeiro@ufpr.br
4Department of Mathematical Sciences, University of Bath, Bath BA2 7AY, UK.
add49@bath.ac.uk
Abstract
Waves with constant vorticity and electrohydrodynamics flows are two topics in fluid dynamics that have
attracted much attention from scientists for both the mathematical challenge and their industrial applica-
tions. The coupling of electric fields and vorticity is of significant research interest. In this paper, we study
the flow structure of steady periodic travelling waves with constant vorticity on a dielectric fluid under the
effect of normal electric fields. Through the conformal mapping technique combined with pseudo-spectral
numerical methods, we develop an approach that allows us to conclude that the flow can have zero, two or
three stagnation points according to variations in the voltage potential. We describe in detail the recircu-
lation zones that emerge together with the stagnation points. Besides, we show that the number of local
maxima of the pressure on the bottom boundary is intrinsically connected to the saddle points.
1 Introduction
The problem of surface water waves dates back to Stokes [40], who first considered surface gravity waves in
deep water with a very large lengthscale in kilometres. The research on this subject has been extended in many
directions. In particular, capillary waves, in which the surface tension is dominant with a typical wavelength
in millimetres, were studied first by Crapper [7] where an exact solution was derived, and then by Kinnersley
[28] who generalised the results to the fluid sheet problem. When gravity and surface tension are equally
important, it leads to a much more mathematical challenging problem named capillary-gravity wave, with a
small lengthscale about a centimetre, which plays a vital role in transferring energy momentum and material
fluxes across the ocean surface [4, 29]. A good monograph can be found in [44] for a review.
The aforementioned works assumed that the fluid is inviscid and the flow is incompressible and irrotational.
In practice, a theory with non-zero vorticity is often more physically realistic. The most common choice on the
vorticity is a constant value. A number of works have been achieved under such assumption for gravity waves.
For example, Constantin and Strauss[5] conducted a rigorous analytical study. Thomas et al. [43] investigated
the modulational instability of such waves in the framework of a nonlinear Schrödinger equation, and recently a
more comprehensive work of stability analysis was achieved by Francius and Kharif [13]. Numerically speaking,
a boundary integral equation method was employed to compute fully nonlinear travelling-wave solutions with
constant vorticity in this context by many authors [39, 42, 45, 46, 47]. In [2], time-dependent simulations
were conducted by a hodograph transformation method, first pioneered by Dyachenko et al. [11], to examine
the modulational instabilities numerically. In the presence of surface tension, i.e. capillary-gravity waves with
constant vorticity, there have been fruitful achievements by different authors. For example, steady-state solutions
were computed by Kang and Vanden-Broeck [26] via a boundary integral method. A Nonlinear Schrödinger
Equation was derived in [25] to study the modulational instabilities. Time-dependent wave interactions in this
context were investigated numerically in [15] in which various travelling-wave solutions were also computed.
All the studies in the last paragraph were interested solely in the shape of the free-surface, and integral
properties such as wave energy and momentum. More recent studies have considered the flow structure beneath
the waves [37, 36]. These structures are particularly interesting in the case of waves with constant vorticity. Of
1
arXiv:2210.10825v1 [physics.flu-dyn] 19 Oct 2022
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 xycoordinate 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
E
wall electrode
y=h
wall electrode
y= 0
V= 0
V=V0
V= 0
2φ= 0,2V= 0.
Dielectric fluid
Permittivity 1
Conducting gas
y=η(x, t)
g
y
x
Figure 1: Schematic of the problem.
X=xct,Y=y. It immediately follows that η(x, t) = η(X). The fluid motion can be characterised by the
potential function φand the voltage potential Vas follows
2φ= 0 in 1< Y < η(X),(2)
2V= 0 in 1< Y < η(X),(3)
X+ (φXγ(η+b)) ηX=φYfor Y=η(X),(4)
φY= 0 for Y=1,(5)
V= 0 for Y=η(X),(6)
V=1for Y=1,(7)
where we have chosen h,ph/g and V0as the reference length, time and voltage potential, respectively. Fur-
thermore, the continuity of pressure on the free-surface yields the dynamic boundary condition, which can be
written as
X+1
2(φ2
X+φ2
Y) + ηγ(η+b)φX+γψ τηXX
(1 + η2
X)3/2+Eb
2|∇V|2=B. (8)
Here, ψis the harmonic conjugate of φ. Parameter τand Ebare the non-dimensional Bond and Electric Bond
numbers respectively, given by
τ=T
ρgh2, Eb=1V2
0
ρgh3.(9)
The hydrodynamic pressure in fluid body is computed by the Bernoulli equation
p=X+1
2(φ2
X+φ2
Y) + Yγ(Y+b)φX+γψ B,(10)
and the stream function can be written as
ψs(X, Y ) := ψ(X, Y )γY 1
2Y+bcY. (11)
3 Linear theory
In this section, a linear theory is developed from the governing equations. A trivial solution takes the form of
φ0(X, Y )=0, V0(X, Y ) = Y, η0(X)=0,(12)
which is perturbed by a small disturbance, namely
η(X) = ˆη,
φ(X, Y ) = ˆ
φ,
V(X, Y ) = Y+ˆ
V ,
(13)
3
摘要:

FlowstructurebeneathperiodicwaveswithconstantvorticityundernormalelectriceldsM.V.Flamarion1,T.Gao2,R.Ribeiro-Jr3&A.Doak41UnidadeAcadêmicadoCabodeSantoAgostinho,UFRPE/RuralFederalUniversityofPernambuco,BR101Sul,CabodeSantoAgostinho-PE,Brazil,54503-900marcelo.amarion@ufrpe.br2SchoolofComputingandMat...

展开>> 收起<<
Flow structure beneath periodic waves with constant vorticity under normal electric fields M. V. Flamarion1 T. Gao2 R. Ribeiro-Jr3 A. Doak4.pdf

共14页,预览3页

还剩页未读, 继续阅读

声明:本站为文档C2C交易模式,即用户上传的文档直接被用户下载,本站只是中间服务平台,本站所有文档下载所得的收益归上传人(含作者)所有。玖贝云文库仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对上载内容本身不做任何修改或编辑。若文档所含内容侵犯了您的版权或隐私,请立即通知玖贝云文库,我们立即给予删除!

相关推荐

分类:图书资源 价格:10玖币 属性:14 页 大小:4.28MB 格式:PDF 时间:2025-05-06

开通VIP享超值会员特权

  • 多端同步记录
  • 高速下载文档
  • 免费文档工具
  • 分享文档赚钱
  • 每日登录抽奖
  • 优质衍生服务
/ 14
客服
关注