scholarly journals Wave-Structure Interaction for a Stationary Surface-Piercing Body Based on a Novel Meshless Scheme with the Generalized Finite Difference Method

Mathematics ◽  
2020 ◽  
Vol 8 (7) ◽  
pp. 1147
Author(s):  
Ji Huang ◽  
Hongguan Lyu ◽  
Chia-Ming Fan ◽  
Jiahn-Hong Chen ◽  
Chi-Nan Chu ◽  
...  

The wave-structure interaction for surface-piercing bodies is a challenging problem in both coastal and ocean engineering. In the present study, a two-dimensional numerical wave flume that is based on a newly-developed meshless scheme with the generalized finite difference method (GFDM) is constructed in order to investigate the characteristics of the hydrodynamic loads acting on a surface-piercing body caused by the second-order Stokes waves. Within the framework of the potential flow theory, the second-order Runge-Kutta method (RKM2) in conjunction with the semi-Lagrangian approach is carried out to discretize the temporal variable of governing equations. At each time step, the GFDM is employed to solve the spatial variable of the Laplace’s equation for the deformable computational domain. The results show that the developed numerical method has good performance in the simulation of wave-structure interaction, which suggests that the proposed “RKM2-GFDM” meshless scheme can be a feasible tool for such and more complicated hydrodynamic problems in practical engineering.

1989 ◽  
Vol 79 (4) ◽  
pp. 1210-1230
Author(s):  
C. R. Daudt ◽  
L. W. Braile ◽  
R. L. Nowack ◽  
C. S. Chiang

Abstract The Fourier method, the second-order finite-difference method, and a fourth-order implicit finite-difference method have been tested using analytical phase and group velocity calculations, homogeneous velocity model calculations for disperson analysis, two-dimensional layered-interface calculations, comparisons with the Cagniard-de Hoop method, and calculations for a laterally heterogeneous model. Group velocity rather than phase velocity dispersion calculations are shown to be a more useful aid in predicting the frequency-dependent travel-time errors resulting from grid dispersion, and in establishing criteria for estimating equivalent accuracy between discrete grid methods. Comparison of the Fourier method with the Cagniard-de Hoop method showed that the Fourier method produced accurate seismic traces for a planar interface model even when a relatively coarse grid calculation was used. Computations using an IBM 3083 showed that Fourier method calculations using fourth-order time derivatives can be performed using as little as one-fourth the CPU time of an equivalent second-order finite-difference calculation. The Fourier method required a factor of 20 less computer storage than the equivalent second-order finite-difference calculation. The fourth-order finite-difference method required two-thirds the CPU time and a factor of 4 less computer storage than the second-order calculation. For comparison purposes, equivalent runs were determined by allowing a group velocity error tolerance of 2.5 per cent numerical dispersion for the maximum seismic frequency in each calculation. The Fourier method was also applied to a laterally heterogeneous model consisting of random velocity variations in the lower half-space. Seismograms for the random velocity model resulted in anticipated variations in amplitude with distance, particularly for refracted phases.


Author(s):  
Harry B. Bingham ◽  
Allan P. Engsig-Karup

This contribution presents our recent progress on developing an efficient solution for fully nonlinear wave-structure interaction. The approach is to solve directly the three-dimensional (3D) potential flow problem. The time evolution of the wave field is captured by integrating the free-surface boundary conditions using a fourth-order Runge-Kutta scheme. A coordinate-transformation is employed to obtain a time-constant spatial computational domain which is discretized using arbitrary-order finite difference schemes on a grid with one stretching in each coordinate direction. The resultant linear system of equations is solved by the GMRES iterative method, preconditioned using a multigrid solution to the linearized, lowest-order version of the matrix. The computational effort and required memory use are shown to scale linearly with increasing problem size (total number of grid points). Preliminary examples of nonlinear wave interaction with variable bottom bathymetry and simple bottom mounted structures are given.


Sign in / Sign up

Export Citation Format

Share Document