Suppression of Irregular Frequency in Multi-Body Problem and Free-Surface Singularity Treatment

Author(s):  
Yujie Liu ◽  
Jeffrey M. Falzarano

Multibody operations are routinely performed in offshore activities. One classical example is the FLNG and LNGC side-by-side offloading case. To understand the phenomenon occurring inside the gap is of growing interest to the offshore industry. One important issue is the existence of the irregular frequency effect. The effect can be confused with the physical resonance. Thus it needs to be removed. An extensive survey of the previous approaches to the irregular frequency problem has been undertaken. The matrix formulated in the boundary integral equations will become nearly singular for some frequencies. The existence of numerical round-off errors will make the matrix still solvable by a direct solver, however will result in unreasonably large values in some aspects of the solution, namely the irregular frequency effect. The removal of the irregular effect is important especially for multi-body hydrodynamic analysis in identifying the physical resonances caused by the configuration of floaters. This paper will mainly discuss the lid method on the internal free surface. To reach a higher accuracy, the singularity resulting from the Green function needs special care. Each term in the wave Green function will be evaluated using the corresponding analysis methods. Specifically, an analytical integral method is proposed to treat the log singularity. Finally, results with and without irregular frequency removal will be shown to demonstrate the effectiveness of our proposed method. The validation cases include mini-boxbarge, boxbarge and cylindrical dock, which has apparent irregular frequency effect in their output results.

Author(s):  
Yujie Liu ◽  
Jeffrey M. Falzarano

Multibody operations are routinely performed in offshore activities, for example, the floating liquefied natural gas (FLNG) and liquefied natural gas carrier (LNGC) side-by-side offloading case. To understand the phenomenon occurring inside the gap is of growing interest to the offshore industry. One important issue is the existence of the irregular frequency effect. The effect can be confused with the physical resonance. Thus, it needs to be removed. An extensive survey of the previous approaches to the irregular frequency problem has been undertaken. The matrix formulated in the boundary integral equations will become nearly singular for some frequencies. The existence of numerical round-off errors will make the matrix still solvable by a direct solver, however, it will result in unreasonably large values in some aspects of the solution, namely, the irregular frequency effect. The removal of the irregular effect is important especially for multibody hydrodynamic analysis in identifying the physical resonances caused by the configuration of floaters. This paper will mainly discuss the lid method on the internal free surface. To reach a higher accuracy, the singularity resulting from the Green function needs special care. Each term in the wave Green function will be evaluated using the corresponding analysis methods. Specifically, an analytical integral method is proposed to treat the log singularity. Finally, results with and without irregular frequency removal will be shown to demonstrate the effectiveness of our proposed method.


Author(s):  
Philippe Helluy ◽  
Sylvain Maire ◽  
Patrice Ravel

Abstract A high order integration method is presented for regular or singular integrands over an integral. This method appears to be very useful to compute the integrals of the green function in the numerical resolution of boundary integral equations.


The question of non-uniqueness in boundary integral equation formu­lations of exterior problems for the Helmholtz equation has recently been resolved with the use of additional radiating multipoles in the definition of the Green function. The present note shows how this modification may be included in a rigorous formalism and presents an explicit choice of co­efficients of the added terms that is optimal in the sense of minimizing the least-squares difference between the modified and exact Green functions.


1999 ◽  
Vol 43 (01) ◽  
pp. 13-24
Author(s):  
M. Landrini ◽  
G. Grytøyr ◽  
O. M. Faltinsen

Fully nonlinear free-surface flows are numerically studied in the framework of the potential theory. The problem is formulated in terms of boundary integral equations which are solved by means of an arbitrary high-order boundary element method based on B-Spline representation of both the geometry and the fluid dynamic variables along the domain boundary. The solution is stepped forward in time either by following Lagrangian points attached to the free surface or by a less conventional scheme in which evolution equations for the B-Spline coefficients are integrated in time. Numerical examples for inner and outer free-surface flows are shown. The accuracy of the numerical solution is assessed either by checking mass and energy conservation or by comparing with reference solutions. Good results are generally obtained. Extended use of the developed algorithm to more applied problems in the context of naval hydrodynamics is now under development.


Author(s):  
Nikolai I. Makarenko

The fully nonlinear problem on unsteady two-dimensional water waves generated by elliptic cylinder, that is horizontally submerged beneath a free surface, is considered. An analytical boundary integral equations method using a version of Milne-Thomson transformation is developed. Boundary equations (the BEq system) determine immediately exact wave elevation and fluid velocity at free surface. Small-time solution expansion is obtained in the case of accelerated cylinder starting from rest.


Author(s):  
Chunmei Xie ◽  
Aurélien Babarit ◽  
François Rongère ◽  
Alain H. Clément

A new acceleration technique for the computation of first order hydrodynamic coefficients for floating bodies in frequency domain and in deep water is proposed. It is based on the classical boundary element method (BEM) which requires solving a boundary integral equation for distributions of sources and/or dipoles and evaluating integrals of Kelvin’s Green function and its derivatives over panels. The Kelvin’s Green function includes two Rankine sources and a wave term. In present study, for the two Rankine sources, analytical integrations of strongly singular kernels are adopted for the linear density distributions. It is shown that these analytical integrations are more accurate and faster than numerical integrations. The wave term is obtained by solving Clément’s ordinary differential equations (ODEs) [1] and an adaptive numerical quadrature is performed for integrations over the panels. It is shown here that the computational time of the wave term by solving the ODEs is greatly reduced compared to the classical integration method [7].


Author(s):  
Yingyi Liu ◽  
Ying Gou ◽  
Bin Teng

The present study aims to develop an efficient numerical method for computing the diffraction and radiation of water waves with horizontal long cylindrical structures, such as floating breakwaters. A higher-order scheme is used to discretize geometry of the structure as well as the relevant physical quantities. As the kernel of this method, Wehausen’s free-surface Green function is calculated by a newly-developed Gauss-Kronrod adaptive quadrature algorithm after elimination of its Cauchy-type singularities. To improve computational efficiency, a Chebyshev approximation approach is applied to a fast calculation of the Green function that needs evaluation thousands of times. In addition, OpenMP parallel technique is used to the formation of influence coefficient matrix, which significantly reduces CPU time. Finally, computations are performed on wave exciting forces and hydrodynamic coefficients for the long cylindrical structures, either floating or submerged. Comparison with other numerical and analytical methods demonstrates good performance of the present method.


Author(s):  
Petr Denisov ◽  
◽  
Anna Balaban ◽  

The article proposes the modification of a technique for assessing the magnetization of permanent magnets from the known field pattern. The identification method is based on solving an ill-conditioned system of linear algebraic equations by the Tikhonov regularization method. The method of boundary integral equations based on scalar potentials is used to compile the matrix of coefficients. The article presents the algorithm that uses parallel computations when performing the most time-consuming operations to reduce the time for solving the inverse problem. In order to check the proposed method, a program was developed that allows to simulate the measurement process: to calculate the direct problem and find the magnetic induction at the points of the air gap, then introduce the error into the "measurement results" and solve the inverse problem. The results of nu-merical experiments that allow us to evaluate the advantages of parallel implementation using the capabilities of modern multi-core processors are presented.


Author(s):  
Hui Li ◽  
Hao Lizhu ◽  
Huilong Ren ◽  
Xiaobo Chen

The solution of hydrodynamic problem with forward speed still has some well-known problems such as high oscillation and slow convergence of the wave term when using a moving and oscillating source as the Green function. Recently, Ten and Chen (2010) has come up with a new method to benefit the merits of both the Rankine source and moving and oscillating source by taking a hemisphere as the control surface which separates the fluid region into two domains, but some troubles have been induced in the process of solution. Therefore, in this paper, a cylindrical surface instead of a hemisphere is selected to be the control surface to make the solution easy, and in this method, the control surface isn’t divided into panels. In the interior domain near the ship, the Rankin Green function is used to simplify the calculation. In the exterior domain some distance from the ship, there is no panels representing the free surface by using the Green function which satisfy the free surface boundary condition. The whole fluid region matches by the condition that the velocity potentials and their normal derivatives in the interior domain and exterior domain are equal on the control surface separately. In this paper, we have validated the Rankine-Kelvin hybrid method is applicable by adopting it to solve the zero speed problem in this work.


Author(s):  
Oscar P. Bruno ◽  
Stephen P. Shipman ◽  
Catalin Turc ◽  
Stephanos Venakides

This work, part I in a two-part series, presents: (i) a simple and highly efficient algorithm for evaluation of quasi-periodic Green functions, as well as (ii) an associated boundary-integral equation method for the numerical solution of problems of scattering of waves by doubly periodic arrays of scatterers in three-dimensional space. Except for certain ‘Wood frequencies’ at which the quasi-periodic Green function ceases to exist, the proposed approach, which is based on smooth windowing functions, gives rise to tapered lattice sums which converge superalgebraically fast to the Green function—that is, faster than any power of the number of terms used. This is in sharp contrast to the extremely slow convergence exhibited by the lattice sums in the absence of smooth windowing. (The Wood-frequency problem is treated in part II.) This paper establishes rigorously the superalgebraic convergence of the windowed lattice sums. A variety of numerical results demonstrate the practical efficiency of the proposed approach.


Sign in / Sign up

Export Citation Format

Share Document