3,577 research outputs found

    Convergence of a Boundary Integral Method for Water Waves

    Get PDF
    We prove nonlinear stability and convergence of certain boundary integral methods for time-dependent water waves in a two-dimensional, inviscid, irrotational, incompressible fluid, with or without surface tension. The methods are convergent as long as the underlying solution remains fairly regular (and a sign condition holds in the case without surface tension). Thus, numerical instabilities are ruled out even in a fully nonlinear regime. The analysis is based on delicate energy estimates, following a framework previously developed in the continuous case [Beale, Hou, and Lowengrub, Comm. Pure Appl. Math., 46 (1993), pp. 1269–1301]. No analyticity assumption is made for the physical solution. Our study indicates that the numerical methods must satisfy certain compatibility conditions in order to be stable. Violation of these conditions will lead to numerical instabilities. A breaking wave is calculated as an illustration

    Asymptotic expansions and fast computation of oscillatory Hilbert transforms

    Full text link
    In this paper, we study the asymptotics and fast computation of the one-sided oscillatory Hilbert transforms of the form H+(f(t)eiΟ‰t)(x)=βˆ’int0∞eiΟ‰tf(t)tβˆ’xdt,Ο‰>0,xβ‰₯0,H^{+}(f(t)e^{i\omega t})(x)=-int_{0}^{\infty}e^{i\omega t}\frac{f(t)}{t-x}dt,\qquad \omega>0,\qquad x\geq 0, where the bar indicates the Cauchy principal value and ff is a real-valued function with analytic continuation in the first quadrant, except possibly a branch point of algebraic type at the origin. When x=0x=0, the integral is interpreted as a Hadamard finite-part integral, provided it is divergent. Asymptotic expansions in inverse powers of Ο‰\omega are derived for each fixed xβ‰₯0x\geq 0, which clarify the large Ο‰\omega behavior of this transform. We then present efficient and affordable approaches for numerical evaluation of such oscillatory transforms. Depending on the position of xx, we classify our discussion into three regimes, namely, x=O(1)x=\mathcal{O}(1) or x≫1x\gg1, 0<xβ‰ͺ10<x\ll 1 and x=0x=0. Numerical experiments show that the convergence of the proposed methods greatly improve when the frequency Ο‰\omega increases. Some extensions to oscillatory Hilbert transforms with Bessel oscillators are briefly discussed as well.Comment: 32 pages, 6 figures, 4 table
    • …
    corecore