Discretized Simulation of Vortex Sheet Evolution with Buoyancy and Surface Tension Effects

AIAA Journal ◽  
1976 ◽  
Vol 14 (11) ◽  
pp. 1517-1523 ◽  
Author(s):  
Robert G. Zalosh
Keyword(s):  
1998 ◽  
Vol 19 (5) ◽  
pp. 1737-1766 ◽  
Author(s):  
Gregory Baker ◽  
André Nachbin
Keyword(s):  

1981 ◽  
Vol 108 ◽  
pp. 159-170 ◽  
Author(s):  
P. N. Shankar

The initial-value problem for the evolution of the interface η(x, t) separating two unbounded, inviscid streams is considered in the framework of linearized analysis. Given the initial shape y = εη0(x) of the interface at t = 0 the objective is to calculate the interface shape η(x,t) for later times. First, it is shown that, if the vortex sheet is of infinite extent, if surface tension is absent and if the two streams are of the same density, the evolution is given by \[ \eta(x,t) = \epsilon(1-\alpha)^{-1}{\rm Re}[\{(1-\alpha)+(1+\alpha)i\}\eta_0\{x-{\textstyle\frac{1}{2}}((1+\alpha)+(1-\alpha)i)t\}], \] where α (≠ 1) is the ratio of the speeds of the streams, provided the initial interface shape εη0(x) is analytic and its Fourier transform decays sufficiently rapidly. An interesting consequence is that it is possible, under certain circumstances, for the interface to develop singularities after a finite time. Next it is shown that when the two streams move at the same speed (α = 1) the growth of η is given by \[ \eta(x,t) = \epsilon\eta_0(x-t)+\epsilon t\,d\eta_0(x-t)/dx \] with mild restrictions on η0(x). The major effect of surface tension, it is found, is to prevent the occurrence of singularities after a finite time, a distinct possibility in its absence. Finally the vortex sheet shed by a semi-infinite flat plate is considered. The unsteady mixed boundary-value problem is formally solved by using parabolic coordinates and Fourier-Laplace transforms.


1982 ◽  
Vol 119 ◽  
pp. 507-532 ◽  
Author(s):  
D. I. Pullin

We consider the behaviour of an interface between two immiscible inviscid incompressible fluids of different density moving under the action of gravity, inertial and interfacial tension forces. A vortex-sheet model of the exact nonlinear two-dimensional motion of this interface is formulated which includes expressions for an appropriate set of integral invariants. A numerical method for solving the vortex-sheet initial-value equations is developed, and is used to study the nonlinear growth of finite-amplitude normal modes for both Kelvin-Helmholtz and Rayleigh-Taylor instability. In the absence of an interfacial or surface-tension term in the integral-differential equation that describes the evolution of the circulation distribution on the vortex sheet, it is found that chaotic motion of, or the appearance of curvature singularities in, the discretized interface profiles prevent the simulations from proceeding to the late-time highly nonlinear phase of the motion. This unphysical behaviour is interpreted as a numerical manifestation of possible ill-posedness in the initial-value equations equivalent to the infinite growth rate of infinitesimal-wavelength disturbances in the linearized stability theory. The inclusion of an interfacial tension term in the circulation equation (which stabilizes linearized short-wavelength perturbations) was found to smooth profile irregularities but only for finite times. While coherent interfacial motion could then be followed well into the nonlinear regime for both the Kelvin-Helmholtz and Rayleigh-Taylor modes, locally irregular behaviour eventually reappeared and resisted subsequent attempts at numerical smoothing or suppression. Although several numerical and/or physical mechanisms are discussed that might produce irregular behaviour of the discretized interface in the presence of an interfacial-tension term, the basic cause of this instability remains unknown. The final description of the nonlinear interface motion thus awaits further research.


2008 ◽  
Vol 18 (4) ◽  
pp. 463-484 ◽  
Author(s):  
F. de la Hoz ◽  
M. A. Fontelos ◽  
L. Vega

Sign in / Sign up

Export Citation Format

Share Document