REEF3D::FNPF: A Flexible Fully Nonlinear Potential Flow Solver

2019 ◽  
Author(s):  
Hans Bihs ◽  
Weizhi Wang ◽  
Tobias Martin ◽  
Arun Kamath

Abstract In situations where the calculation of ocean wave propagation and impact on offshore structures is required, fast numerical solvers are desired in order to find relevant wave events in a first step. After the identification of the relevant events, Computational Fluid Dynamics (CFD) based Numerical Wave Tanks (NWT) with an interface capturing two-phase flow approach can be used to resolve the complex wave structure interaction, including breaking wave kinematics. CFD models emphasize detail of the hydrodynamic physics, which makes them not the ideal candidate for the event identification due to the large computational resources involved. In the current paper a new numerical wave model is represented that solves the Laplace equation for the flow potential and the nonlinear kinematic and dynamics free surface boundary conditions. This approach requires reduced computational resources compared to CFD based NWTs. In contrast to existing approaches, the resulting fully nonlinear potential flow solver REEF3D::FNPF uses a σ-coordinate grid for the computations. Solid boundaries are incorporated through a ghost cell immersed boundary method. The free surface boundary conditions are discretized using fifth-order WENO finite difference methods and the third-order TVD Runge-Kutta scheme for time stepping. The Laplace equation for the potential is solved with Hypres stabilized bi-conjugated gradient solver preconditioned with geometric multi-grid. REEF3D::FNPF is fully parallelized following the domain decomposition strategy and the MPI communication protocol. The model is successfully tested for wave propagation benchmark cases for shallow water conditions with variable bottom as well as deep water.

Author(s):  
Hans Bihs ◽  
Weizhi Wang ◽  
Csaba Pakozdi ◽  
Arun Kamath

Abstract In situations where the calculation of ocean wave propagation and impact on structures are required, fast numerical solvers are desired in order to find relevant wave events. Computational fluid dynamics (CFD)-based numerical wave tanks (NWTs) emphasize on the hydrodynamic details such as fluid–structure interaction, which make them less ideal for the event identification due to the large computational resources involved. Therefore, a computationally efficient numerical wave model is needed to identify the events both for offshore deep-water wave fields and coastal wave fields where the bathymetry and coastline variations have strong impact on wave propagation. In the current paper, a new numerical wave model is represented that solves the Laplace equation for the flow potential and the nonlinear kinematic and dynamics free surface boundary conditions. This approach requires reduced computational resources compared to CFD-based NWTs. The resulting fully nonlinear potential flow solver REEF3D::FNPF uses a σ-coordinate grid for the computations. This allows the grid to follow the irregular bottom variation with great flexibility. The free surface boundary conditions are discretized using fifth-order weighted essentially non-oscillatory (WENO) finite difference methods and the third-order total variation diminishing (TVD) Runge–Kutta scheme for time stepping. The Laplace equation for the potential is solved with Hypre’s stabilized bi-conjugated gradient solver preconditioned with geometric multi-grid. REEF3D::FNPF is fully parallelized following the domain decomposition strategy and the message passing interface (MPI) communication protocol. The numerical results agree well with the experimental measurements in all tested cases and the model proves to be efficient and accurate for both offshore and coastal conditions.


2020 ◽  
Vol 159 ◽  
pp. 103717
Author(s):  
Nikta Iravani ◽  
Peyman Badiei ◽  
Maurizio Brocchini

Author(s):  
Zaibin Lin ◽  
Ling Qian ◽  
Wei Bai ◽  
Zhihua Ma ◽  
Hao Chen ◽  
...  

Abstract A 3-Dimensional numerical wave tank based on the fully nonlinear potential flow theory has been developed in OpenFOAM, where the Laplace equation of velocity potential is discretized by Finite Volume Method. The water surface is tracked by the semi-Eulerian-Lagrangian method, where water particles on the free surface are allowed to move vertically only. The incident wave is generated by specifying velocity profiles at inlet boundary with a ramp function at the beginning of simulation to prevent initial transient disturbance. Additionally, an artificial damping zone is located at the end of wave tank to sufficiently absorb the outgoing waves before reaching downstream boundary. A five-point smoothing technique is applied at the free surface to eliminate the saw-tooth instability. The proposed wave model is validated against theoretical results and experimental data. The developed solver could be coupled with multiphase Navier-Stokes solvers in OpenFOAM in the future to establish an integrated versatile numerical wave tank for studying efficiently wave structure interaction problems.


2004 ◽  
Vol 126 (4) ◽  
pp. 287-296
Author(s):  
R. E. Baddour ◽  
W. Parsons

We are studying numerically the problem of generation and propagation of long-crested gravity waves in a tank containing an incompressible inviscid homogeneous fluid initially at rest with a horizontal free surface of finite extent and of infinite depth. A nonorthogonal curvilinear coordinate system, which follows the free surface, is constructed and the full nonlinear kinematic and dynamic free surface boundary conditions are utilized in the algorithm. “Wavemakers” are modeled using both the Dirichlet and Neumann lateral boundary conditions and a full comparison is given. Overall, the Dirichlet model was more stable than the Neumann model, with an upper limit steepness S=2A/λ of 0.08 using good resolution compared with the Neumann’s maximum of 0.05.


Sign in / Sign up

Export Citation Format

Share Document