scholarly journals Singularity formation in three-dimensional vortex sheets

2003 ◽  
Vol 15 (1) ◽  
pp. 147-172 ◽  
Author(s):  
Thomas Y. Hou ◽  
Gang Hu ◽  
Pingwen Zhang
2002 ◽  
Vol 10 (02) ◽  
pp. 155-182 ◽  
Author(s):  
OLIVER V. ATASSI ◽  
AMR A. ALI

Inflow/Outflow conditions are formulated for time-harmonic waves in a duct governed by the Euler equations. These conditions are used to compute the propagation of acoustic and vortical disturbances and the scattering of vortical waves into acoustic waves by an annular cascade. The outflow condition is expressed in terms of the pressure, thus avoiding the velocity discontinuity across any vortex sheets. The numerical solutions are compared with the analytical solutions for acoustic and vortical wave propagation with and without the presence of vortex sheets. Grid resolution studies are also carried out to discern the truncation error of the numerical scheme from the error associated with numerical reflections at the boundary. It is observed that even with the use of exponentially accurate boundary conditions, the dispersive characteristics of the numerical scheme may result in small reflections from the boundary that slow convergence. Finally, the three-dimensional interaction of a wake with a flat plate cascade is computed and the aerodynamic and aeroacoustic results are compared with those of lifting surface methods.


1987 ◽  
Vol 178 ◽  
pp. 195-219 ◽  
Author(s):  
Douglas G. Dommermuth ◽  
Dick K. P. Yue

A numerical method is developed for nonlinear three-dimensional but axisymmetric free-surface problems using a mixed Eulerian-Lagrangian scheme under the assumption of potential flow. Taking advantage of axisymmetry, Rankine ring sources are used in a Green's theorem boundary-integral formulation to solve the field equation; and the free surface is then updated in time following Lagrangian points. A special treatment of the free surface and body intersection points is generalized to this case which avoids the difficulties associated with the singularity there. To allow for long-time simulations, the nonlinear computational domain is matched to a transient linear wavefield outside. When the matching boundary is placed at a suitable distance (depending on wave amplitude), numerical simulations can, in principle, be continued indefinitely in time. Based on a simple stability argument, a regriding algorithm similar to that of Fink & Soh (1974) for vortex sheets is generalized to free-surface flows, which removes the instabilities experienced by earlier investigators and eliminates the need for artificial smoothing. The resulting scheme is very robust and stable.For illustration, three computational examples are presented: (i) the growth and collapse of a vapour cavity near the free surface; (ii) the heaving of a floating vertical cylinder starting from rest; and (iii) the heaving of an inverted vertical cone. For the cavity problem, there is excellent agreement with available experiments. For the wave-body interaction calculations, we are able to obtain and analyse steady-state (limit-cycle) results for the force and flow field in the vicinity of the body.


1998 ◽  
Vol 9 (6) ◽  
pp. 607-634
Author(s):  
P. BASSANINI ◽  
C. M. CASCIOLA ◽  
M. R. LANCIA ◽  
R. PIVA

Steady incompressible inviscid flow past a three-dimensional multiconnected (toroidal) aerofoil with a sharp trailing edge TE is considered, adopting for simplicity a linearized analysis of the vortex sheets that collect the released vorticity and form the trailing wake. The main purpose of the paper is to discuss the uniqueness of the bounded flow solution and the role of the eigenfunction. A generic admissible flow velocity u has an unbounded singularity at TE; and the physical flow solution requires the removal of the divergent part of u (the Kutta condition). This process yields a linear functional equation along the trailing edge involving both the normal vorticity ω released into the wake, and the multiplicative factor of the eigenfunction, a1. Uniqueness is then shown to depend upon the topology of the trailing edge. If δTE=[empty ], as, for example, in an annular-aerofoil configuration, both ω and a1 are uniquely determined by the Kutta condition, and the bounded flow u is unique. If δTE≠[empty ], as, for example, in a connected-wing configuration, there is an infinity of bounded flows, parametrized by a1. Numerical results of relevance for these typical configurations are presented to show the different role of the eigenfunction in the two cases.


2013 ◽  
Vol 723 ◽  
pp. 289-317 ◽  
Author(s):  
Andrew L. Stewart ◽  
Paul J. Dellar

AbstractWe analyse the hyperbolicity of our multilayer shallow water equations that include the complete Coriolis force due to the Earth’s rotation. Shallow water theory represents flows in which the vertical shear is concentrated into vortex sheets between layers of uniform velocity. Such configurations are subject to Kelvin–Helmholtz instabilities, with arbitrarily large growth rates for sufficiently short-wavelength disturbances. These instabilities manifest themselves through a loss of hyperbolicity in the shallow water equations, rendering them ill-posed for the solution of initial value problems. We show that, in the limit of vanishingly small density difference between the two layers, our two-layer shallow water equations remain hyperbolic when the velocity difference remains below the same threshold that also ensures the hyperbolicity of the standard shallow water equations. Direct calculation of the domain of hyperbolicity becomes much less tractable for three or more layers, so we demonstrate numerically that the threshold for the velocity differences, below which the three-layer equations remain hyperbolic, is also unchanged by the inclusion of the complete Coriolis force. In all cases, the shape of the domain of hyperbolicity, which extends outside the threshold, changes considerably. The standard shallow water equations only lose hyperbolicity due to shear parallel to the direction of wave propagation, but the complete Coriolis force introduces another mechanism for loss of hyperbolicity due to shear in the perpendicular direction. We demonstrate that this additional mechanism corresponds to the onset of a transverse shear instability driven by the non-traditional components of the Coriolis force in a three-dimensional continuously stratified fluid.


1977 ◽  
Vol 99 (1) ◽  
pp. 187-196
Author(s):  
M. Ribaut

This study deals with the calculation of the three-dimensional steady flow of a compressible and viscous fluid through a turbomachine. Using the potential theory, a quasi three-dimensional method is developed, which reduces the simplifying assumptions and yields high numerical accuracy. The method includes all main effects induced by the primary vortex system, in particular the influence of the casing and blade boundary layers on the cascade circulation. The use of two vortex sheets, representing the boundary surface, makes the calculation of flows having strong developed or separated boundary layers possible and allows a direct calculation of the secondary vorticity. Computed divergent and rotational flows are presented and compared with exact solutions or experimental data.


1995 ◽  
Vol 300 ◽  
pp. 339-366 ◽  
Author(s):  
Takashi Ishihara ◽  
Yukio Kaneda

The evolution of a small but finite three-dimensional disturbance on a flat uniform vortex sheet is analysed on the basis of a Lagrangian representation of the motion. The sheet at time t is expanded in a double periodic Fourier series: R(λ1, λ2, t) = (λ1, λ2, 0) + Σn,mAn,m exp[i(nλ1 + δmλ2)], where λ1 and λ2 are Lagrangian parameters in the streamwise and spanwise directions, respectively, and δ is the aspect ratio of the periodic domain of the disturbance. By generalizing Moore's analysis for two-dimensional motion to three dimensions, we derive evolution equations for the Fourier coefficients An,m. The behaviour of An,m is investigated by both numerical integration of a set of truncated equations and a leading-order asymptotic analysis valid at large t. Both the numerical integration and the asymptotic analysis show that a singularity appears at a finite time tc = O(lnε−1) where ε is the amplitude of the initial disturbance. The singularity is such that An,0 = O(tc−1) behaves like n−5/2, while An,±1 = O(εtc) behaves like n−3/2 for large n. The evolution of A0,m(spanwise mode) is also studied by an asymptotic analysis valid at large t. The analysis shows that a singularity appears at a finite time t = O(ε−1) and the singularity is characterized by A0,2k ∝ k−5/2 for large k.


Sign in / Sign up

Export Citation Format

Share Document