New Weighted Total Acceleration on Momentum Euler Equation for Formulating Water Wave Dispersion Equation

In this present study, weighted total acceleration for Kinematic Free Surface Boundary Condition (KFSBC) and in momentum Euler equation was formulated. Furthermore, by using both aforementioned equations, the nonlinear water wave dispersion equation was then formulated. The wavelength obtained from dispersion equation is determined by weighting coefficient. The weighting coefficient value was determined by using the maximum wave height and critical wave steepness criteria which have been obtained from the previous studies.

Improvement of Non-Hydrostatic Hydrodynamic Solution Using a Novel Free-Surface Boundary Condition

Hydrodynamic models based on the RANS equation are well-established tools to simulate three-dimensional free surface flows in large aquatic ecosystems. However, when the ratio of vertical to horizontal motion scales is not small, a non-hydrostatic approximation is needed to represent these processes accurately. Increasing efforts have been made to improve the efficiency of non-hydrostatic hydrodynamic models, but these improvements require higher implementation and computational costs. In this paper, we proposed a novel free-surface boundary condition based on a fictional sublayer at the free-surface (FSFS). We applied the FSFS approach at a finite difference numerical discretization with a fractional step framework, which uses a Neumann type of boundary condition to apply a hydrostatic relation in the top layer. To evaluate the model performance, we compared the Classic Boundary Condition Approach (CBA) and the FSFS approach using two numerical experiments. The experiments tested the model’s phase error, capability in solving wave celerity and simulate non-linear wave propagation under different vertical resolution scenarios. Our results showed that the FSFS approach had a lower phase error (2 to 5 times smaller) than CBA with a little additional computational cost (ca. 7% higher). Moreover, it can better represent wave celerity and frequency dispersion with 2 times fewer layers and low mean computational cost (CBA δ t = 2.62 s and FSFS δ t = 1.22 s).

3D frequency-domain elastic wave modeling with the spectral element method using a massively parallel direct solver

ABSTRACT Complex topography, the free-surface boundary condition, and anelastic properties of media should be accounted for in the frame of onshore geophysical prospecting imaging, such as full-waveform inversion (FWI). In this context, an accurate and efficient forward-modeling engine is mandatory. We have performed 3D frequency-domain anisotropic elastic wave modeling by using the highly accurate spectral element method and a sparse multifrontal direct solver. An efficient approach similar to computing the matrix-vector product in the time domain is used to build the matrix. We validate the numerical results by comparing with analytical solutions. A parallel direct solver, the sparse direct multifrontal massively parallel solver (MUMPS), is used to solve the linear system. We find that a hybrid implementation of message passing interface and open multiprocessing is more efficient in flops and memory cost. The influence of the deformed mesh, free-surface boundary condition, and heterogeneity of media on MUMPS performance is negligible. Complexity analysis suggests that the memory complexity of MUMPS agrees with the theoretical order [Formula: see text] (or [Formula: see text] with an efficient matrix reordering method) for an [Formula: see text] grid when nontrivial topography is considered. With the available resources, we conduct a moderate scale modeling with a subset of the SEAM Phase II Foothills model, where 60 wavelengths in the [Formula: see text]-axis are propagated. Computing one gradient of FWI based on this model using the frequency-domain modeling is shown to require similar or fewer computational resources than what would be required for a time-domain solver, depending on the number of sources, while larger memory is necessary. An estimation of the increasing trend indicates that approximately 20 Tb of memory would be required for a [Formula: see text] wavelength modeling. The limit of MUMPS scalability hinders the application to larger scale applications.

Free-Surface Effects on Interaction of Multiple Ships Moving at Different Speeds

Ships often have to pass each other in proximity in harbor areas and waterways in dense shipping-traffic environment. Hydrodynamic interaction occurs when a ship is overtaking (or being overtaken) or encountering other ships. Such an interactive effect could be magnified in confined waterways, e.g., shallow and narrow rivers. Since Yeung published his initial work on ship interaction in shallow water, progress on unsteady interaction among multiple ships has been slow, though steady, over the following decades. With some exceptions, nearly all the published studies on ship-to-ship problem neglected free-surface effects, and a rigid-wall condition has often been applied on the water surface as the boundary condition. When the speed of the ships is low, this assumption is reasonably accurate as the hydrodynamic interaction is mainly induced by near-field disturbances. However, in many maneuvering operations, the encountering or overtaking speeds are actually moderately high (Froude number Fn > 0.2, where <inline-graphic xlink:href="josr10180089inf1.tif"/>, U is ship speed, g is the gravitational acceleration, and L is the ship length), especially when the lateral separation between ships is the order of ship length. Here, the far-field effects arising from ship waves can be important. The hydrodynamic interaction model must take into account the surface-wave effects. Classical potential-flow formulation is only able to deal with the boundary value problem when there is only one speed involved in the free-surface boundary condition. For multiple ships traveling with different speeds, it is not possible to express the free-surface boundary condition by a single velocity potential. Instead, a superposition method can be applied to account for the velocity field induced by each vessel with its own and unique speed. The main objective of the present article is to propose a rational superposition method to handle the unsteady free-surface boundary condition containing two or more speed terms, and validate its feasibility in predicting the hydrodynamic behavior in ship encountering. The methodology used in the present article is a three-dimensional boundary-element method based on a Rankine-type (infinite-space) source function, initially introduced by Bai and Yeung. The numerical simulations are conducted by using an in-house‐developed multibody hydrodynamic interaction program “MHydro.” Waves generated and forces (or moments) are calculated when ships are encountering or passing each other. Published model-test results are used to validate our calculations, and very good agreement has been observed. The numerical results show that free-surface effects need to be taken into account for Fn > 0.2.

Investigation of Free Surface Damping Models With Applications to Gap Resonance Problems

In this paper two methods for modelling the damping in a narrow gap are investigated. The first method is called the Pressure Damping Model. This method has been used in studies of wave energy devices. An attractive feature of this model is that the modified input is directly related to the energy dissipation in the gap, which means that if this dissipation is estimated the input to the model can be obtained directly. The idea of the method is to add a pressure input in the gap to suppress the resonant motion. A challenge with the method is that it contains a non-linear term. The second method is the Newtonian Cooling damping model. The method is based on introducing a dissipation term in the free surface boundary condition. This dissipation term contains a coefficient which is not directly related to the energy dissipation. Hence this method is not so easy to relate directly to the estimated energy dissipation. An advantage with this method is that it is linear and hence can be expected to be more robust. In the first part of the paper a 2-dimensional problem is investigated using both methods. In addition to the numerical performance and robustness, much focus is put on investigation of the energy balance in the solution, and we attempt to relate both models to the energy dissipation in the gap. In the second part the Newtonian cooling method is implemented in a 3-dimensional potential flow solver and it is shown that the method provides a robust way to handle the resonance problem. The method will give rise to a modified set of equations which are described. Two different problems are investigated with the 3D solver. First we look at a side-by-side problem, where the 3D results are also compared with 2D results. Finally, the moonpool problem is investigated by two different 3D solvers, a classical Green’s function based method and a Rankine solver. It is also shown how this damping model can be combined with a similar model on the internal waterplane to remove irregular frequencies.

Inherently unstable internal gravity waves due to resonant harmonic generation

Here we show that there exist internal gravity waves that are inherently unstable, that is, they cannot exist in nature for a long time. The instability mechanism is a one-way (irreversible) harmonic-generation resonance that permanently transfers the energy of an internal wave to its higher harmonics. We show that, in fact, there are a countably infinite number of such unstable waves. For the harmonic-generation resonance to take place, the nonlinear terms in the free surface boundary condition play a pivotal role, and the instability does not occur in a linearly stratified fluid if a simplified boundary condition, such as a rigid lid or a linearized boundary condition, is employed. Harmonic-generation resonance presented here provides a mechanism for the transfer of internal wave energy to the higher-frequency part of the spectrum hence affecting, potentially significantly, the evolution of the internal waves spectrum.

The Vertical Mode Decomposition of Surface and Internal Tides in the Presence of a Free Surface and Arbitrary Topography

The method of decomposing surface and internal tides determines the expression for internal tide energy, energy flux, and energy conversion. The de facto standard is to define surface tides as depth-averaged pressure and horizontal velocity and internal tides as the residuals. This decomposition, which is equivalent to projecting motion onto vertical modes that obey a rigid lid, is known to produce spurious energy conversion CS through movement of the free surface. Here, motion is instead projected onto modes that obey a linear, free-surface boundary condition. The free-surface modes are shown to obey a more complicated orthogonality condition than rigid-lid modes but are still straightforward to calculate numerically. The resulting decomposition (i) completely eliminates spurious energy conversion CS and (ii) leads to a more precise expression for topographic internal tide generation C, which now depends on horizontal gradients in the vertical structure of the surface tide. Numerical simulations and rough global estimates indicate that corrections to C are a maximum of a few percent. However, CS produces spurious energy flux divergences/convergences in the open ocean, which are the same order of magnitude [O(1–10) mW m−2] as open-ocean internal tide energy dissipation.

Improvement of Rankine Panel Method for Seakeeping Prediction of a Ship in Low Frequency Region

Rankine panel methods have been studied for solving 3D seakeeping problems of a ship with forward speed and oscillatory motions. Nevertheless, there is a drawback in the numerical method for satisfying the radiation condition of outgoing waves at low frequencies, because the waves generated ahead of a ship reflect from the outward computational boundary and smear the flow around the ship. The so-called panel shift technique has been adopted in the frequency-domain Rankine panel method, which is effective when the generated waves propagate downstream of a ship. In this paper, in addition to this conventional panel shift method, Rayleigh’s artificial friction is introduced in the free-surface boundary condition to suppress longer wave components in a computational region apart from the ship. With this practical new method, it is shown that there is no prominent wave reflection from the side and/or upstream computational boundaries even in the range of low frequencies. As a consequence, the unsteady pressure, hydrodynamic forces, wave-induced ship motions, added resistance are computed with reasonable accuracy even in following waves and in good agreement with measured results in the experiment using a bulk carrier model which is also conducted for the present study.