scholarly journals The effects of surface tension on the initial development of a free surface adjacent to an accelerated plate

2015 ◽  
Vol 776 ◽  
pp. 37-73 ◽  
Author(s):  
J. Uddin ◽  
D. J. Needham

When a vertical rigid plate is uniformly accelerated horizontally from rest into an initially stationary layer of inviscid incompressible fluid, the free surface will undergo a deformation in the locality of the contact point. This deformation of the free surface will, in the early stages, cause a jet to rise up the plate. An understanding of the local structure of the free surface in the early stages of motion is vital in many situations, and has been developed in detail by King & Needham (J. Fluid Mech., vol. 268, 1994, pp. 89–101). In this work we consider the effects of introducing weak surface tension, characterized by the inverse Weber number $\mathscr{W}$, into the problem considered by King & Needham. Our approach is based upon matched asymptotic expansions as $\mathscr{W}\rightarrow 0$. It is found that four asymptotic regions are needed to describe the problem. The three largest regions have analytical solutions, whilst a numerical method based on finite differences is used to solve the time-dependent harmonic boundary value problem in the last region. Our results identify the local structure of the jet near the vicinity of the contact point, and we highlight a number of key features, including the height of this jet as well as its thickness and strength. We also present some preliminary experimental results that capture the spatial structure near the contact point, and we then show promising comparisons with the theoretical results obtained within this paper.

2007 ◽  
Vol 578 ◽  
pp. 67-84 ◽  
Author(s):  
D. J. NEEDHAM ◽  
J. BILLINGHAM ◽  
A. C. KING

The free-surface deformation and flow field caused by the impulsive horizontal motion of a rigid vertical plate into a horizontal strip of inviscid incompressible fluid, initially at rest, is studied in the small time limit using the method of matched asymptotic expansions. It is found that three different asymptotic regions are necessary to describe the flow. There is a main, O(1) sized, outer region in which the flow is singular at the point where the free surface meets the plate. This leads to an inner region, centred on the point where the free surface initially meets the plate, with size of O(-t log t). To resolve the singularities that arise in this inner region, it is necessary to analyse further the flow in an inner-inner region, with size of O(t), again centred upon the wetting point of the nascent rising jet. The solutions of the boundary value problems in the two largest regions are obtained analytically. The solution of the parameter-free nonlinear boundary value problem that arises in the inner-inner region is obtained numerically.


1969 ◽  
Vol 35 (1) ◽  
pp. 77-96 ◽  
Author(s):  
L. M. Perko

This paper considers the large-amplitude symmetric and asymmetric irrota-tional motion of an inviscid incompressible fluid with a liquid—vapour interface in an accelerating container of revolution. A combined analytical—numerical method which involves no linearizations in the hydrodynamical equations and applies to all but surface-tension dominated motions is used to compute a variety of such motions. One important aspect of this non-linear method is that it accurately determines the initial development of surface instabilities such as breakers near the wall of the container.


2001 ◽  
Vol 12 (3) ◽  
pp. 233-252 ◽  
Author(s):  
S. K. WILSON ◽  
B. R. DUFFY ◽  
S. H. DAVIS

In this paper two similarity solutions describing a steady, slender, symmetric dry patch in an infinitely wide liquid film draining under gravity down an inclined plane are obtained. The first solution, which predicts that the dry patch has a parabolic shape and that the transverse profile of the free surface always has a monotonically increasing shape, is appropriate for weak surface-tension effects and far from the apex of the dry patch. The second solution, which predicts that the dry patch has a quartic shape and that the transverse profile of the free surface has a capillary ridge near the contact line and decays in an oscillatory manner far from it, is appropriate for strong surface-tension effects (in particular, when the plane is nearly vertical) and near (but not too close) to the apex of the dry patch. With the average volume flux per unit width (or equivalently with the uniform height of the layer far from the dry patch) prescribed, both solutions contain a free parameter. For each value of this parameter there is a unique solution in the first case and either no solution or a one-parameter family of solutions in the second case. The solutions capture some of the qualitative features observed in experiments.


2018 ◽  
Vol 841 ◽  
pp. 146-166 ◽  
Author(s):  
M. T. Gallagher ◽  
D. J. Needham ◽  
J. Billingham

We consider the problem of a rigid plate, inclined at an angle $\unicode[STIX]{x1D6FC}\in (0,\unicode[STIX]{x03C0}/2)$ to the horizontal, accelerating uniformly from rest into, or away from, a semi-infinite strip of inviscid, incompressible fluid under gravity. Following on from Gallagher et al. (J. Fluid Mech., vol. 841, 2018, pp. 109–145) (henceforth referred to as GNB), it is of interest to analyse the well-posedness and stability of the principal flow with respect to perturbations in the initially horizontal free surface close to the plate contact point. In particular we find that the solution to the principal unperturbed problem, denoted by [IBVP] in GNB, is well-posed and stable with respect to perturbations in initial data in the region of interest, localised close to the contact point of the free surface and the plate, when the plate is accelerated with dimensionless acceleration $\unicode[STIX]{x1D70E}\geqslant -\cot \,\unicode[STIX]{x1D6FC}$, while the solution to [IBVP] is ill-posed with respect to such perturbations in the initial data, when the plate is accelerated with dimensionless acceleration $\unicode[STIX]{x1D70E}<-\cot \,\unicode[STIX]{x1D6FC}$. The physical source of the ill-posedness of the principal problem [IBVP], when $\unicode[STIX]{x1D70E}<-\cot \,\unicode[STIX]{x1D6FC}$, is revealed to be due to the leading-order problem in the innermost region localised close to the initial contact point being in the form of a local Rayleigh–Taylor problem. As a consequence of this mechanistic interpretation we anticipate that, when the plate is accelerated with $\unicode[STIX]{x1D70E}<-\cot \,\unicode[STIX]{x1D6FC}$, the inclusion of weak surface tension effects will restore well-posedness of the problem [IBVP] which will, however, remain temporally unstable.


2018 ◽  
Vol 841 ◽  
pp. 109-145 ◽  
Author(s):  
M. T. Gallagher ◽  
D. J. Needham ◽  
J. Billingham

The free surface and flow field structure generated by the uniform acceleration (with dimensionless acceleration $\unicode[STIX]{x1D70E}$) of a rigid plate, inclined at an angle $\unicode[STIX]{x1D6FC}\in (0,\unicode[STIX]{x03C0}/2)$ to the exterior horizontal, as it advances ($\unicode[STIX]{x1D70E}>0$) or retreats ($\unicode[STIX]{x1D70E}<0$) from an initially stationary and horizontal strip of inviscid incompressible fluid under gravity, are studied in the small-time limit via the method of matched asymptotic expansions. This work generalises the case of a uniformly accelerating plate advancing into a fluid as studied by Needham et al. (Q. J. Mech. Appl. Maths, vol. 61 (4), 2008, pp. 581–614). Particular attention is paid to the innermost asymptotic regions encompassing the initial interaction between the plate and the free surface. We find that the structure of the solution to the governing initial boundary value problem is characterised in terms of the parameters $\unicode[STIX]{x1D6FC}$ and $\unicode[STIX]{x1D707}$ (where $\unicode[STIX]{x1D707}=1+\unicode[STIX]{x1D70E}\tan \unicode[STIX]{x1D6FC}$), with a bifurcation in structure as $\unicode[STIX]{x1D707}$ changes sign. This bifurcation in structure leads us to question the well-posedness and stability of the governing initial boundary value problem with respect to small perturbations in initial data in the innermost asymptotic regions, the discussion of which will be presented in the companion paper Gallagher et al. (J. Fluid Mech. vol. 841, 2018, pp. 146–166). In particular, when $(\unicode[STIX]{x1D6FC},\unicode[STIX]{x1D707})\in (0,\unicode[STIX]{x03C0}/2)\times \mathbb{R}^{+}$, the free surface close to the initial contact point remains monotone, and encompasses a swelling jet when $(\unicode[STIX]{x1D6FC},\unicode[STIX]{x1D707})\in (0,\unicode[STIX]{x03C0}/2)\times [1,\infty )$ or a collapsing jet when $(\unicode[STIX]{x1D6FC},\unicode[STIX]{x1D707})\in (0,\unicode[STIX]{x03C0}/2)\times (0,1)$. However, when $(\unicode[STIX]{x1D6FC},\unicode[STIX]{x1D707})\in (0,\unicode[STIX]{x03C0}/2)\times \mathbb{R}^{-}$, the collapsing jet develops a more complex structure, with the free surface close to the initial contact point now developing a finite number of local oscillations, with near resonance type behaviour occurring close to a countable set of critical plate angles $\unicode[STIX]{x1D6FC}=\unicode[STIX]{x1D6FC}_{n}^{\ast }\in (0,\unicode[STIX]{x03C0}/2)$ ($n=1,2,\ldots$).


Water Waves ◽  
2021 ◽  
Author(s):  
M. D. Groves

AbstractIn the applied mathematics literature solitary gravity–capillary water waves are modelled by approximating the standard governing equations for water waves by a Korteweg-de Vries equation (for strong surface tension) or a nonlinear Schrödinger equation (for weak surface tension). These formal arguments have been justified by sophisticated techniques such as spatial dynamics and centre-manifold reduction methods on the one hand and variational methods on the other. This article presents a complete, self-contained account of an alternative, simpler approach in which one works directly with the Zakharov–Craig–Sulem formulation of the water-wave problem and uses only rudimentary fixed-point arguments and Fourier analysis.


2021 ◽  
Vol 126 (1) ◽  
Author(s):  
Alex Doak ◽  
Jean-Marc Vanden-Broeck

AbstractThis paper concerns the flow of fluid exiting a two-dimensional pipe and impacting an infinite wedge. Where the flow leaves the pipe there is a free surface between the fluid and a passive gas. The model is a generalisation of both plane bubbles and flow impacting a flat plate. In the absence of gravity and surface tension, an exact free streamline solution is derived. We also construct two numerical schemes to compute solutions with the inclusion of surface tension and gravity. The first method involves mapping the flow to the lower half-plane, where an integral equation concerning only boundary values is derived. This integral equation is solved numerically. The second method involves conformally mapping the flow domain onto a unit disc in the s-plane. The unknowns are then expressed as a power series in s. The series is truncated, and the coefficients are solved numerically. The boundary integral method has the additional advantage that it allows for solutions with waves in the far-field, as discussed later. Good agreement between the two numerical methods and the exact free streamline solution provides a check on the numerical schemes.


2000 ◽  
Vol 406 ◽  
pp. 337-346 ◽  
Author(s):  
L. ENGEVIK

The instabilities of a free surface shear flow are considered, with special emphasis on the shear flow with the velocity profile U* = U*0sech2 (by*). This velocity profile, which is found to model very well the shear flow in the wake of a hydrofoil, has been focused on in previous studies, for instance by Dimas & Triantyfallou who made a purely numerical investigation of this problem, and by Longuet-Higgins who simplified the problem by approximating the velocity profile with a piecewise-linear profile to make it amenable to an analytical treatment. However, none has so far recognized that this problem in fact has a very simple solution which can be found analytically; that is, the stability boundaries, i.e. the boundaries between the stable and the unstable regions in the wavenumber (k)–Froude number (F)-plane, are given by simple algebraic equations in k and F. This applies also when surface tension is included. With no surface tension present there exist two distinct regimes of unstable waves for all values of the Froude number F > 0. If 0 < F [Lt ] 1, then one of the regimes is given by 0 < k < (1 − F2/6), the other by F−2 < k < 9F−2, which is a very extended region on the k-axis. When F [Gt ] 1 there is one small unstable region close to k = 0, i.e. 0 < k < 9/(4F2), the other unstable region being (3/2)1/2F−1 < k < 2 + 27/(8F2). When surface tension is included there may be one, two or even three distinct regimes of unstable modes depending on the value of the Froude number. For small F there is only one instability region, for intermediate values of F there are two regimes of unstable modes, and when F is large enough there are three distinct instability regions.


Sign in / Sign up

Export Citation Format

Share Document