• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.35 no.2
• /
• pp.41-48
• /
• 2023

In this study, simplified linear numerical method that can simulate wave generation and transformation by a moving bottom is introduced. Numerical analysis is conducted in wave number domain after continuity equation, linear dynamic and kinematic free surface boundary conditions and linear kinematic bottom boundary condition are Fourier transformed, and the results are expressed in space domain by an inverse Fourier transform. In the wavenumber domain, the dynamic free water surface boundary condition and the kinematic free water surface boundary condition are numerically calculated, and the velocity potential in the mean water level (z = 0) satisfies the continuity equation and the kinematic bottom boundary condition. Wave generation and transformation are investigated when the triangular and rectangular shape of bottoms move periodically. The results of the simplified numerical method are compared with the results of previous analytical solutions and agree well with them. Stability of numerical results according to the calculation time interval (Δt) and the calculation wave number interval (Δk) was also investigated. It was found that the numerical results were appropriate when Δt ≤ T(period)/1000 and Δk ≤ π/100.

• The Journal of the Acoustical Society of Korea
• /
• v.25 no.3E
• /
• pp.131-136
• /
• 2006

The PML developed for the radio wave propagation is a powerful numerical domain truncation technique. We perform an analytic study on the reflection from the PML inserted in the ocean bottom. In the ocean bottom, we show the PML to have the improved performance but simultaneously the degeneration below the critical angle of the fast ocean bottom. The degeneration of the PML can be simply relaxed by stretching the thickness of the PML or putting the attenuation coefficient to the ocean bottom. As a better solution, we propose the improved truncation technique based on the PML and the non-local boundary condition. Finally, we apply the PML to the acoustic wave propagation using split-step Pade PE solver. For the problems of the ocean waveguide, the numerical efficiency of the PML is examined and the usefulness of the PML is confirmed.

• Journal of Korean Society of Water and Wastewater
• /
• v.25 no.5
• /
• pp.627-634
• /
• 2011

In the present study, the flow behaviors of square jets surface discharged and submerged discharged into shallow water were each simulated using computational fluid dynamics, and the results were compared. As for the verification of the models, the results of the hydraulic experiment conducted by Sankar, et al. (2009) were used. According to the results of the verification, the present application of computational fluid dynamics to the flow analysis of square jets discharged into shallow water was valid. As for the wall jet, which is one form of submerged discharges, at the bottom wall boundary, the peak velocity of the jet rapidly moved from the center of the jet to the bottom wall boundary due to the restriction of jet entrainment and the no-slip condition of the bottom wall boundary, and, as for the surface discharge, because jet entrainment is limited on the free water surface, the peak velocity of the jet moved from the center of the jet to the free water surface. This is because jet entrainment is restricted at the bottom wall boundary and the surface so that the momentum of the central core of the jet is preserved for considerable time at the bottom wall boundary and the surface. In addition, due to the effect of the bottom wall boundary and the free water surface, the jet discharged into shallow water had a smaller velocity diminution rate near the discharge outlet than did the free jet; at a location where it was so distant from the discharge outlet that the vertical profile of the velocity was nearly equal (b/x =20~30), moreover, it had a far smaller velocity diminution rate than did the free jet due to the effect of the finite depth.

• Proceedings of the KSME Conference
• /
• 2001.06e
• /
• pp.506-511
• /
• 2001

This paper presents how redevelopment of the boundary layer in a backward-facing step flow is affected by boundary conditions imposed on velocity at the inlet, top and exit of the flow. A two-dimensional, laminar, incompressible flow over a backward-facing step with an open top boundary has been computed by using numerical methods of second-order time and spatial accuracy and a fractional-step method that guarantees a divergence-free velocity field at all time. The inlet velocity profile above the step is of Blasius type. Along the top boundary, shear-tree and Dirichlet conditions on the streamwise velocity were considered and at the exit fully-developed and convective boundary conditions were examined. (The vertical velocity at all boundaries were assumed to be zero explicitly or implicitly.) From the computed flow fields, the reattachment on the bottom side of shear layer separated from the tip of the step and succeeding redevelopment of the boundary layer were investigated.

• Coupled systems mechanics
• /
• v.12 no.4
• /
• pp.365-389
• /
• 2023

Dynamic analysis of a typical concrete gravity dam-reservoir system is formulated by FE-(FE-TE) approach (i.e., Finite Element-(Finite Element-Truncation Element)). In this technique, dam and reservoir are discretized by plane solid and fluid finite elements. Moreover, the H-W (i.e., Hagstrom-Warburton) high-order condition imposed at the reservoir truncation boundary. This task is formulated by employing a truncation element at that boundary. It is emphasized that reservoir far-field is excluded from the discretized model. The formulation is initially reviewed which was originally proposed in a previous study. Thereafter, the response of Pine Flat dam-reservoir system is studied due to horizontal and vertical ground motions for two types of reservoir bottom conditions of full reflective and absorptive. It should be emphasized that study is carried out under high order of H-W condition applied on the truncation boundary. The initial part of study is focused on the time harmonic analysis. In this part, it is possible to compare the transfer functions against corresponding responses obtained by FE-(FE-HE) approach (referred to as exact method). Subsequently, the transient analysis is carried out. In that part, it is only possible to compare the results for low and high normalized reservoir length cases. Therefore, the sensitivity of results is controlled due to normalized reservoir length values.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.34 no.3
• /
• pp.58-71
• /
• 2022

Analytical solutions for two-dimensional velocities of regular waves generated by bottom wave makers in a flume were derived in this study. Triangular and rectangular bottom wave makers were adopted. The velocity potential was derived based on the linear wave theory with the bottom moving boundary condition, kinematic and dynamic free surface boundary conditions. Then, analytical solutions of two-dimensional particle velocities were derived from the velocity potential. The velocity potential and two-dimensional particle velocities which were derived as complex integral equations were numerically calculated. The solutions showed physically valid results as velocities of regular waves generated by bottom wave makers in a flume.

• Journal of Advanced Research in Ocean Engineering
• /
• v.1 no.3
• /
• pp.157-167
• /
• 2015

In order to provide simple and accurate wave theory in design of offshore structure, an analytical approximation is introduced in this paper. The solution is limited to flat bottom having a constant water depth. Water is considered as inviscid, incompressible and irrotational. The solution satisfies the continuity equation, bottom boundary condition and non-linear kinematic free surface boundary condition exactly. Error for dynamic condition is quite small. The solution is suitable in description of breaking waves. The solution is presented with closed form and dispersion relation is also presented with closed form. In the last century, there have been two main approaches to the nonlinear problems. One of these is perturbation method. Stokes wave and Cnoidal wave are based on the method. The other is numerical method. Dean's stream function theory is based on the method. In this paper, power series method was considered. The power series method can be applied to certain nonlinear differential equations (initial value problems). The series coefficients are specified by a nonlinear recurrence inherited from the differential equation. Because the non-linear wave problem is a boundary value problem, the power series method cannot be applied to the problem in general. But finite number of coefficients is necessary to describe the wave profile, truncated power series is enough. Therefore the power series method can be applied to the problem. In this case, the series coefficients are specified by a set of equations instead of recurrence. By using the set of equations, the nonlinear wave problem has been solved in this paper.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1991.04a
• /
• pp.8-14
• /
• 1991

The finite element technique incorporating infinite elements is applied to analyzing the general three dimensional wave-structure interaction problems within the limits of linear wave theory. The hydrodynamic farces are assumed to be inertially dominated, and viscous effects are neglected. In order to analyze the corresponding boundary value problems efficiently, two types of elements are developed. One is the infinite element for modeling the radiation condition at infinity, and the other is the fictitious bottom boundary element for the case of deep water. To validate those elements, numerical analyses are performed for several floating structures. Comparisons with the results from culler available solution methods show that the present method incorporating tile infinite and the fictitious bottom boundary elements gives good results.

• Earthquakes and Structures
• /
• v.3 no.2
• /
• pp.97-115
• /
• 2012

The energy-transmitting boundary, which is used in the well-known finite element method (FEM) program FLUSH, is quite efficient for the earthquake response analysis of buildings considering soil-structure interaction. However, it is applicable only in the frequency domain. The author proposed methods for transforming frequency dependent impedance into the time domain, and studied the time domain transform of the boundary. In this paper, first, the estimation methods for both the halfspace condition under the bottom of the soil model and the pseudo three-dimensional effect were studied with the time domain transmitting boundary. Next, response behavior when using the boundary was studied in detail using a practical soil and building model. The response accuracy was compared with those using viscous boundary, and the boundary that considers the excavation force. Through these studies, the accuracy and efficiency of the proposed time domain transmitting boundary were confirmed.

• Structural Engineering and Mechanics
• /
• v.40 no.5
• /
• pp.665-688
• /
• 2011

This paper presents a theoretical analysis for the free vibration of rectangular tanks partially filled with an ideal liquid. Wet dynamic displacements of the tanks are approximated by combining the orthogonal polynomials satisfying the boundary conditions, since the rectangular tanks are composed of four rectangular plates. The classical boundary conditions of the tanks at the top and bottom ends are considered, such as clamped, simply supported, and clamped-free boundary conditions. As the facing rectangular plates are assumed to be geometrically and structurally identical, the vibration modes of the facing plates of the tanks can be divided into two categories: symmetric and antisymmetric modes with respect to the planes passing through the center of the tanks and perpendicular to the free liquid surface. The liquid displacement potentials satisfying the Laplace equation and liquid boundary conditions are derived, and the wet dynamic modal functions of a quarter of the tanks can be expanded by the finite Fourier transform for compatibility requirements along the contacting surfaces between the tanks and liquid. An eigenvalue problem is derived using the Rayleigh-Ritz method. Consequently, the wet natural frequencies of the rectangular tanks can be extracted. The proposed analytical method is verified by observing an excellent agreement with three-dimensional finite element analysis results. The effects of the liquid level and boundary condition at the top and bottom edges are investigated.