• The Journal of the Acoustical Society of Korea
• /
• v.19 no.4
• /
• pp.84-92
• /
• 2000

In this paper, numerical scattering analysis for submerged scatterers is performed using finite and infinite elements. Unbounded domain is truncated into finite domain and finite elements are used in the domain. Infinite elements, So called Infinite Wave Envelope Elements (IWEE) which possess wave-like behavior, are used to take into account the infinite domain on the truncated boundary Scattering from rigid sphere is taken as an example and the effects of the order and mesh size of finite elements, size of finite element model and the order of IWEE are investigated. Quadratic finite element, refined mesh and higher order IWEE are recommended to improve the non-reflection boundary condition in the numerical scattering analysis.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.37 no.3
• /
• pp.187-195
• /
• 2024

This paper presents two-dimensional boundary value problems of the stress-based gradient elasticity within the smoothed finite element method (S-FEM) framework. Gradient elasticity is introduced to address the limitations of classical elasticity, particularly its struggle to capture size-dependent mechanical behavior at the micro/nano scale. The Ru-Aifantis theorem is employed to overcome the challenges of high-order differential equations in gradient elasticity. This theorem effectively splits the original equation into two solvable second-order differential equations, enabling its incorporation into the S-FEM framework. The present method utilizes a staggered scheme to solve the boundary value problems. This approach efficiently separates the calculation of the local displacement field (obtained over each smoothing domain) from the non-local stress field (computed element-wise). A series of numerical tests are conducted to investigate the influence of the internal length scale, a key parameter in gradient elasticity. The results demonstrate the effectiveness of the proposed approach in smoothing stress concentrations typically observed at crack tips and dislocation lines.

• Journal of Ocean Engineering and Technology
• /
• v.12 no.2 s.28
• /
• pp.97-110
• /
• 1998

The numerical damping and dispersion error characteristics associated with difference schemes and a panel shift method used for the calculation of steady free surface flows by a panel method are an analysed in this paper. First, 12 finite difference operators used for the double model flow by Letcher are applied to a two dimensional cylinder with the Kelvin free surface condition and the numerical errors with these schemes are compared with those by the panel shift method. Then, 3-D waves due to a submerged source are calculated by the difference schemes, the panel shift method and also by a higher order boundary element method(HOBEM). Finally, the waves and wave resistance for Wigley's hull are calculated with these three schemes. It is shown that the panel shift method is free of numerical damping and dispersion error and performs better than the difference schemes. However, it can be concluded that the HOBEM also free of the numerical damping and dispersion error is the most stable, accurate and efficient.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.8 no.1
• /
• pp.37-43
• /
• 1996

Breaking waves were generated in a 2-D flume. A piston-type wavemaker was operated in accordance with signals which consist of elementary harmonics with appropriate phase differences. These phase differences were estimated by using a linear wave theory so that wave crests were to be concentrated at the same position. The stroke of wavemaker was controlled to create plunging-type breaking waves. The signal with small amplitude could not generate breaking waves. In the case of moderate amplitudes, various breaking waves could be obtained. Most of breaking waves were spilling type. Only when the wavemaker was operated with appropriate amplitude, plunging-type breaking waves were generated. The parameters of breaking waves are the wave steepness and the frequency bandwidth. If the central frequency was low, breaking waves were not generated. Based on experimental data, we found that the wave height of breaking inception was H = 0.0113 gT$^2$. We also made computations by using a mixed Euler-Lagrangian scheme under the assumption of potential flow. The numerical results show good agreements with tank measurements.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.8 no.2
• /
• pp.174-184
• /
• 1996

In the present study, for fully-open rectangular harbors, sensitivity analyses are made for the major parameters which are relevant to the practical application of a hybrid element model widely used fur the analysis of harbor oscillation. The results show that it is desirable to extend the finite element region to the area in which depth change is not large and that it is appropriate to take the depth of the outer region for analytic solution as the average along the boundary between the two regions. It is expected that the number of Fourier components of the analytic solution may not be important for a constant-depth simple-shaped harbor but its significance may increase for harbors of varying depth and complex geometry. It is found that the effect of incident wave direction is not significant for the first resonance mode but its effect becomes important as the bottom slope increases, especially for the higher resonance modes.

• Journal of the Society of Naval Architects of Korea
• /
• v.35 no.3
• /
• pp.1-13
• /
• 1998

In the present paper, the hydrodynamic characteristics of three dimensional hydrofoils moving with a constant speed below the free surface using a higher-order boundary element method based on 9-node Lagrangian curvilinear elements are investigated. A bi-quadratic spline scheme is employed to improve the numerical results on the free surface. To validate the present scheme, the calculated results are compared with the analytic solutions for a submerged sphere and a spheroid showing a good agreement. For the validation of the hydrofoil study, the computed lift and drag of a hydrofoil having $NACA64_{1}A412$ section with aspect ratio(A.R.) of 4 are compared with the experimental data by Wadlin et al.[28]. The comparison covers a number of variations of angle of attack and submergence depth. Then, using an A.R. hydrofoil with NACA0012 section, the free surface on the lift and drag are investigated and these are compared with the previous results. The wave elevations and patterns created by the aforementioned submerged bodies are also investigated with Froude numbers and submergences.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.13 no.3
• /
• pp.1352-1358
• /
• 2012

The dynamic analysis of delaminated composite structures is carried out based on the higher order plate theory. In the finite element (FE) formulation, the seven degrees of freedom per each node are used with transformations in order to fit the displacement continuity conditions at the delamination region. The boundaries of the instability regions are determined using the method proposed by Bolotin. The numerical results obtained for skew plates are in good agreement with those reported by other investigators. The new results for delaminated skew plate structures in this study mainly show the effect of the interactions between the geometries and other various parameters.

• The Journal of the Acoustical Society of Korea
• /
• v.31 no.1
• /
• pp.19-28
• /
• 2012

The time-domain solution from the Kirchhoff integral equation for an exterior problem is not unique at certain eigen-frequencies associated with the fictitious internal modes as happening in frequency-domain analysis. One of the solution methods is the CHIEF (Combined Helmholtz Integral Equation Formulation) approach, which is based on employing additional zero-pressure constraints at some interior points inside the body. Although this method has been widely used in frequency-domain boundary element method due to its simplicity, it was not used in time-domain analysis. In this work, the CHIEF approach is formulated appropriately for time-domain acoustic boundary element method by constraining the unknown surface pressure distribution at the current time, which was obtained by setting the pressure at the interior point to be zero considering the shortest retarded time between boundary nodes and interior point. Sound radiation of a pulsating sphere was used as a test example. By applying the CHIEF method, the low-order fictitious modes could be damped down satisfactorily, thus solving the non-uniqueness problem. However, it was observed that the instability due to high-order fictitious modes, which were beyond the effective frequency, was increased.

• Journal of the Society of Naval Architects of Korea
• /
• v.40 no.5
• /
• pp.53-59
• /
• 2003

Ocean space utilization using VLFS(Very Large Floating Structures) can provide environmental impact free space by allowing sea water flow freely through the floating structure. Use of Pontoon type VLFS for that purpose needs employment of breakwaters for reduction of wave effects. Therefore, in order to maximize advantage of environmental impact free structure, the breakwater should be the one that can allow water flow freely through it, too. In this paper hydroelastic response of a pontoon type structure is analyzed considering breakwaters which allow water flow through its opening at bottom of the breakwaters. Mode superposition technique is used for solving equation of flexible body while interactions between the pontoon and breakwaters is considered based on generalized mode concept. Bi-quadratic nine node higher-order boundary element method is adopted for more accurate numerical treatment near sharp edged body shape. Performance of various combinations of breakwaters is investigated.

• Journal of Ocean Engineering and Technology
• /
• v.25 no.1
• /
• pp.8-13
• /
• 2011

A compact array system of small buoys is used for wave energy extraction. To evaluate the performance of this system, hydrodynamic analysis is carried out in regular waves using the higher order boundary element method. The motion response of each buoy is calculated considering hydrodynamic interactions caused by other buoys. The effect of energy extraction device is modeled as a linear damping load. The efficiencies of energy conversion are compared using the various sizes and arrangements of the array system and the damping coefficients for energy extraction. The increase in size or the packing ratio of the system gives better efficiency. However, the wave condition and the cost for the system should be considered to optimize performance from the perspective of engineering and economics. The proposed nondimensionalized damping coefficient for energy extraction is 0.1~0.5.