Nonlinear numerical simulation on extreme-wave kinematics

A fully nonlinear numerical model based on a time-domain higher-order boundary element method (HOBEM) is founded to simulate the kinematics of extreme waves. In the model, the fully nonlinear free surface boundary conditions are satisfied and a semi-mixed Euler-Lagrange method is used to track free surface; a fourth-order Runga-Kutta technique is adopted to refresh the wave elevation and velocity potential on the free surface at each time step; an image Green function is used in the numerical wave tank so that the integrations on the lateral surfaces and bottom are excluded. The extreme waves are generated by the method of wave focusing. The physical experiments are carried out in a wave flume. On the horizontal velocity of the measured point, numerical solutions agree well with experimental results. The characteristics of the nonlinear extreme-wave kinematics and the velocity distribution are studied here.     

Verification of A Numerical Harbour Wave Model

A numerical model for wave propagation in a harbour is verified by use of physical models.The extended time-dependent mild slope equation is employed as the governing equation,and the model is solved by use of ADI method containing the relaxation factor.Firstly,the reflection coefficient of waves in front of rubble-mound breakwaters under oblique incident waves is determined through physical model tests,and it is regarded as the basis for simulating partial reflection boundaries of the numerical model.Then model tests on refraction,diffraction and reflection of waves in a harbour are performed to measure wave height distribution.Comparative results between physical and numerical model tests show that the present numerical model can satisfactorily simulate the propagation of regular and irregular waves in a harbour with complex topography and boundary conditions.     

A study on dynamical features of air-sea coupling waves in the tropics

The dynamical features of air-sea coupling waves and their stabilities in a simple coupled air-sea model in the tropics have been studied with respect to interaction occurring among different types of the free waves in the o-cean and in the atmosphere. It is pointed out that there exist a stable and an unstable air-sea interaction modes in the tropical coupled system , respectively. The propagation of the unstable mode relies greatly on the zonal space scale, i. e. only for wave length ranging from 5 000 km to 10 000 km can the disturbance unstably move slowly eastward. The waves that slowly propagate unstably eastward agree well with the observational facts. Finally,it is also proposed that the interaction between Kelvin wave in one medium and Rossby wave in another medium is a necessary condition for the occurrence of destabilization of the coupled air-sea system in the tropics.     

Development of A Fully Nonlinear Numerical Wave Tank

A fully nonlinear numerical wave tank (NWT) based on the solution of the σ-transformed Navier-Stokes equation is developed in this study. The numerical wave is generated from the inflow boundary, where the surface elevation and/or velocity are specified by use of the analytical solution or the laboratory data. The Sommerfeld/Orlanski radiation condition in eonjunetion with an artificial damping zone is applied to reduce wave reflection from the outflow boundary. The whole numerical solution procedures are split into three steps, i.e., advection, diffusion and propagation, and a new method, the Lagrange-Euler Method, instead of the MAC or VOF method, is introduced to solve the free surface elevation at the new time step. Several typical wave cases, including solitary waves, regular waves and irregular waves, are simulated in the wave tank. The robustness and accuracy of the NWF are verified by the good agreement between the numerical results and the linear or nonlinear analytical solutions. This research will be further developed by study of wave-wave, wave-current, wave-structure or wave-jet interaction in the future.     

Wave-tide-surge coupled simulation for typhoon Maemi

The main task of this study focuses on studying the effect of wave-current interaction on currents, storm surge and wind wave as well as effects of current induced wave refraction and current on waves by using numerical models which consider the bottom boundary layer and sea surface roughness parameter for shallow and smooth bed area around Korean Peninsula. The coupled system (unstructured-mesh SWAN wave and ADCIRC) run on the same unstructured mesh. This identical and homogeneous mesh allows the physics of wave-circulation interactions to be correctly resolved in both models. The unstructured mesh can be applied to a large domain allowing all energy from deep to shallow waters to be seamlessly followed. There is no nesting or overlapping of structured wave meshes, and no interpolation is required. In response to typhoon Maemi (2003), all model components were validated independently, and shown to provide a faithful representation of the system’s response to this storm. The waves and storm surge were allowed to develop on the continental shelf and interact with the complex nearshore environment. The resulting modeling system can be used extensively for prediction of the typhoon surge. The result show that it is important to incorporate the wave-current interaction effect into coastal area in the wave-tide-surge coupled model. At the same time, it should consider effects of depth-induced wave breaking, wind field, currents and sea surface elevation in prediction of waves. Specially, we found that: (1) wave radiation stress enhanced the current and surge elevation otherwise wave enhanced nonlinear bottom boundary layer decreased that, (2) wind wave was significantly controlled by sea surface roughness thus we cautiously took the experimental expression. The resulting modeling system can be used for hindcasting (prediction) the wave-tide-surge coupled environments at complex coastline, shallow water and fine sediment area like areas around Korean Peninsula.     

Numerical solution of a mathematical model for water wavesin large coastal areas

on the evolution equation for water waves,a mathematical model for wave propagation in large mild-slope areas is derived.The model is solved by the finite difference method with the staggered grid system.The computational results are in good agreement with experimental data and show that the model can obtain better results with relatively coarser grids.The model can be used to simulate water wave propagation in large coastal areas and can be efficiently solved without much programming effort.     

Investigation on the hydrodynamic performance of an ultra deep turret-moored FLNG system

Hydrodynamic performance of an ultra deep turret-moored Floating Liquefied Natural Gas(FLNG) system is investigated.Hydrodynamic modeling of a turret-moored FLNG system,in consideration of the coupling effects of the vessel and its mooring lines,has been addressed in details.Based on the boundary element method,a 3-D panel model of the FLNG vessel and the related free water surface model are established,and the first-order and second-order mean-drift wave loads and other hydrodynamic coefficients are calculated.A systematic model test program consisting of the white noise wave test,offset test and irregular wave test combined with current and wind,etc.is performed to verify the numerical model.Owing to the depth limit of the water basin,the model test is carried out for the hydrodynamics of the FLNG coupled with only the truncated mooring system.The numerical simulation model features well the hydrodynamic performance of the FLNG system obtained from the model tests.The hydrodynamic characteristics presented in both the numerical simulations and the physical model tests would serve as the guidance for the ongoing project of FLNG system.     

A Finite Element Solution of Wave Forces on Submerged Horizontal Circular Cylinders

In this paper, a numerical model is established for estimating the wave forces on a submerged horizontal circular cylinder. For predicting the wave motion, a set of two-dimensional Navier-Stokes equations is solved numerically with a finite element method. In order to track the moving non-linear wave surface boundary, the Navier-Stokes equations are discretized in a moving mesh system. After each computational time step, the mesh is modified according to the changed wave surface boundary. In order to stabilize the numerical procedure, a three-step finite element method is applied in the time integration. The water sloshing in a tank and wave propagation over a submerged bar are simulated for the first time to validate the present model. The computational results agree well with the analytical solution and the experimental data.Finally, the model is applied to the simulation of interaction between waves and a submerged horizontal circular cylinder.The effects of the KC number and the cylinder depth on the wave forces are studied.     

Power Maximization of A Point Absorber Wave Energy Converter Using Improved Model Predictive Control

This paper considers controlling and maximizing the absorbed power of wave energy converters for irregular waves. With respect to physical constraints of the system, a model predictive control is applied. Irregular waves’ behavior is predicted by Kalman filter method. Owing to the great influence of controller parameters on the absorbed power, these parameters are optimized by imperialist competitive algorithm. The results illustrate the method’s efficiency in maximizing the extracted power in the presence of unknown excitation force which should be predicted by Kalman filter.     

A 3-D time-domain coupled model for nonlinear waves acting on a box-shaped ship fixed in a harbor

A 3-D time-domain numerical coupled model is developed to obtain an efficient method for nonlinear waves acting on a box-shaped ship fixed in a harbor. The domain is divided into the inner domain and the outer domain. The inner domain is the area beneath the ship and the flow is described by the simplified Euler equations. The remaining area is the outer domain and the flow is defined by the higher-order Boussinesq equations in order to consider the nonlinearity of the wave motions. Along the interface boundaries between the inner domain and the outer domain, the volume flux is assumed to be continuous and the wave pressures are equal. Relevant physical experiment is conducted to validate the present model and it is shown that the numerical results agree with the experimental data. Compared the coupled model with the flow in the inner domain governed by the Laplace equation, the present coupled model is more efficient and its solution procedure is simpler, which is particularly useful for the study on the effect of the nonlinear waves acting on a fixed box-shaped ship in a large harbor.     

Wave prediction in a port using a fully nonlinear Boussinesq wave model

A wave forecasting system using FUNWAVE-TVD which is based on the fully nonlinear Boussinesq equations by Chen(2006) was developed to provide an accurate wave prediction in the Port of Busan, South Korea. This system is linked to the Korea Operational Oceanographic System(KOOS) developed by Park et al.(2015). The computational domain covers a region of 9.6 km×7.0 km with a grid size of 2 m in both directions, which is sufficient to resolve short waves and dominant sea states. The total number of grid points exceeds 16 millions,making the model computational expensive. To provide real-time forecasting, an interpolation method, which is based on pre-calculated results of FUNWAVE-TVD and SWAN forecasting results at the FUNWAVE-TVD offshore boundary, was used. A total of 45 cases were pre-calculated, which took 71 days on 924 computational cores of a Linux cluster system. Wind wave generation and propagation from the deep water were computed using the SWAN in KOOS. SWAN results provided a boundary condition for the FUNWAVE-TVD forecasting system. To verify the model, wave observations were conducted at three locations inside the port in a time period of more than 7 months. A model/model comparison between FUNWAVE-TVD and SWAN was also carried out. It is found that, FUNWAVE-TVD improves the forecasting results significantly compared to SWAN which underestimates wave heights in sheltered areas due to incorrect physical mechanism of wave diffraction, as well as large wave heights caused by wave reflections inside the port.     

An Improved Coupling of Numerical and Physical Models for Simulating Wave Propagation简

An improved coupling of numerical and physical models for simulating 2D wave propagation is developed in this paper. In the proposed model, an unstructured finite element model （FEM） based Boussinesq equations is applied for the numerical wave simulation, and a 2D piston-type wavemaker is used for the physical wave generation. An innovative scheme combining fourth-order Lagrange interpolation and Runge-Kutta scheme is described for solving the coupling equation. A Transfer function modulation method is presented to minimize the errors induced from the hydrodynamic invalidity of the coupling model and/or the mechanical capability of the wavemaker in area where nonlinearities or dispersion predominate. The overall performance and applicability of the coupling model has been experimentally validated by accounting for both regular and irregular waves and varying bathymetry. Experimental results show that the proposed numerical scheme and transfer function modulation method are efficient for the data transfer from the numerical model to the physical model up to a deterministic level.     

A Model for Nonlinear Irregular Waves in Coastal Area

Based on the refraction-diffraction theory of irregular waves in the waters of slowly-varying cur-rents and depths,and the generation dissipation theory of wind wave,a model for nonlinear irregular wavesin coastal area is developed.In light of the specific conditions of coastal wave character and engineering ap-plication,a practical mathematical model for the nonlinear irregular waves is presented.with directional spec-trum in coastal area.Coast effect,refraction,whitecapping.bottom friction.current,wind and nonlinearaction are considered in this model.The numerical methods and schemes for wave refraction ray,energy con-servation of propagation,energy balance of the generation and dissipation of wind waves have been studied.Finally,the model is used for the directional wave spectrum computation in the Daya Bay.Compared withthe measured data with 956 wave bouys in the Daya Bay,the model results are in good agreement with themeasured results.     

Three-dimensional numerical simulation of obfique wave action on vertical walls

A three-dimensional numerical model based on the potential theory was developed to study the oblique wave action on vertical walls. A source term inside the domain was used to generate incident waves and outgoing waves were dissipated by sponge layers and transmitted by radiation boundaries. The finite difference method was used to solve the governing equations and boundary conditions in the regular transformed domain in σ-coordinate. Satisfactory agreements between the numerical predictions and experimental results of wave force were obtained. It is concluded that the maximum wave force acting on the vertical walls is induced by the obliquely incident waves rather than the normally incident waves.     

Mooring tension and motion characteristics of a submerged fish reef with net in waves and currents using numerical analysis

A numerical model was used to analyze the motion response and mooring tension of a submerged fish reef system. The system included a net attached to a rigid structure suspended up from the bottom with a single, high tension mooring by fixed flotation. The analysis was performed by using a Morison equation type finite element model configured with truss elements. Input forcing parameters into the model consisted of both regular and irregular waves, with and without a steady current. Heave, surge and pitch dynamic calculations of the reef structure were made. Tension response results of the attached mooring line were also computed. Results were analyzed in both the time and frequency domain in which appropriate, linear transfer functions were calculated. The influence of the current was more evident in the tension and heave motion response data. This is most likely the result of the large buoyancy characteristics of the reef structure and the length of the mooring cable. Maximum mooring component tension was found to be 13.9 kN and occurred when the reef was subjected to irregular waves with a co-linear current of 1.0 m/s velocity. The results also showed that the system had little damping (in heave) with damped natural periods of 2.8 s. This combination of system characteristics promotes a possible resonating situation in typical open sea conditions with similar wave periods.     

Modeling of propagation and transformation of transient nonlinear waves on a current

A novel theoretical approach is applied to predict the propagation and transformation of transient nonlinear waves on a current. The problem was solved by applying an eigenfunction expansion method and the derived semi-analytical solution was employed to study the transformation of wave profile and the evolution of wave spectrum arising from the nonlinear interactions of wave components in a wave train which may lead to the formation of very large waves. The results show that the propagation of wave trains is significantly affected by a current. A relatively small current may substantially affect wave train components and the wave train shape. This is observed for both opposing and following current. The results demonstrate that the application of the nonlinear model has a substantial effect on the shape of a wave spectrum. A train of originally linear and very narrow-banded waves changes its one-peak spectrum to a multi-peak one in a fairly short distance from an initial position. The discrepancies between the wave trains predicted by applying the linear and nonlinear models increase with the increasing wavelength and become significant in shallow water even for waves with low steepness. Laboratory experiments were conducted in a wave flume to verify theoretical results. The free-surface elevations recorded by a system of wave gauges are compared with the results provided by the nonlinear model. Additional verification was achieved by applying a Fourier analysis and comparing wave amplitude spectra obtained from theoretical results with experimental data. A reasonable agreement between theoretical results and experimental data is observed for both amplitudes and phases. The model predicts fairly well multi-peak spectra, including wave spectra with significant nonlinear wave components.     

Numerical Modelling of Standing Waves with Three-Dimensional Non-Linear Wave Propagation Model

Based on the theoretical high-order model with a dissipative term for non-linear and dispersive wave in water of varying depth, a 3-D mathematical model of non-linear wave propagation is presented. The model, which can be used to calculate the wave particle velocity and wave pressure, is suitable to the complicated topography whose relative depth ratio of the characteristic water depth to the characteristic wavelength in deep-water) is equal to or smaller than one. The governing equations are discretized with the improved 2-D Crank-Nicolson method in which the first-order derivatives are corrected by Taylor series expansion, .and the general boundary conditions with an arbitrary reflection coefficient and phase shift are adopted in the model. The surface elevation, horizontal and vertical velocity components and wave pressure of standing waves are numerically calculated. The results show that the numerical model can effectively simulate the complicated standing waves, and the general boundary conditions     

A Numerical Study on Water Waves Generated by a A Submerged Moving Body in a A Two-Llayer Fluid System

This is a numerical study on the time development of surface waves generated by a submerged body moving steadily in a two-layer fluid system, in which a layer of water is underlain by a layer of viscous mud. The fully nonlinear Navier–Stokes equations are solved on FLUENT with the Volume-of-Fluid (VOF) multiphase scheme in order to simulate the free surface waves as well as the water–mud interface waves as functions of time. The numerical model is validated by mimicking a reported experiment in a one-layer system before it is applied to a two-layer system. It is found that the presence of bottom mud in a water layer can lead to large viscous damping of the surface waves. To investigate the problem systematically, the effects of the Froude number and the mud layer thickness, density and viscosity relative to those of water are evaluated and discussed in detail.     

Numerical Simulation of Sloshing Using the MPS-FSI Method with Large Eddy Simulation

A numerical model has been developed to study sloshing of turbulent flow in a tank with elastic baffles. The Moving-Particle Semi-implicit method(MPS) is a kind of meshless Lagrangian calculation method. The large eddy simulation(LES) approach is employed to model the turbulence by using the Smagorinsky Sub-Particle Scale(SPS)closure model. This paper uses MPS-FSI method with LES to simulate the interaction between free surface flow and a thin elastic baffle in sloshing. Then, the numerical model is validated, and the numerical solution has good agreement with experimental data for sloshing in a tank with elastic baffles. Furthermore, under external excitations,the MPS is applied to viscous laminar flow and turbulent flow, with both the deformation of elastic baffles and the wave height of the free surface are compared with each other. Besides, the impact pressure with/without baffles and wave height of free surface are investigated and discussed in detail. Finally, preliminary simulations are carried out in the damage problem of elastic baffles, taking the advantage of the MPS-FSI method in computations of the fluid–structure interaction with large deformation.     

Dynamic Responses of Pontoons Connecting Process of Mobile Offshore Base

The dynamic responses of two pontoons while connected with each other in irregular waves are calculated by means of three-dimensional (3-D) potential flow theory.The computation is to find the optimal status for connection at a certain sea state.On the basis of the relative motion of two pontoons in irregular waves,Visual FORTRAN programming language,as well as OpenGL (Open Graphics Library),is used to develop a set of virtual reality system of the relative motion of two pontoons,which is fully interactive with realistic effect.The transfinite interpolation scheme is applied for the mesh generation of wave surface,and the wave motion is simulated by surface elevation and calculated by 3-D potential flow theory.