• Journal of the Society of Naval Architects of Korea
• /
• v.47 no.3
• /
• pp.291-305
• /
• 2010

A three-dimensional time-domain calculation method is of crucial importance in prediction of the motions and wave loads of a ship advancing in a severe irregular sea. The exact solution of the free surface wave-ship interaction problem is very complicated because of the essentially nonlinear boundary conditions. In this paper, an approximate body nonlinear approach based on the three-dimensional time-domain forward-speed free-surface Green function has been presented. The Froude-Krylov force and the hydrostatic restoring force are calculated over the instantaneous wetted surface of the ship while the forces due to the radiation and scattering potentials over the mean wetted surface. The time-domain radiation and scattering potentials have been obtained from a time invariant kernel of integral equations for the potentials which are discretized according to the second-order boundary element method (Hong and Hong 2008). The diffraction impulse-response functions of the Wigley seakeeping model advancing in transient head waves at various Froude numbers have been presented. A simulation of coupled heave-pitch motion of a long rectangular barge advancing in regular head waves of large amplitude has been carried out. Comparisons between the linear and the approximate body nonlinear numerical results of motions and wave loads of the barge at a nonzero Froude number have been made.

• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.18 no.7
• /
• pp.819-827
• /
• 2007

A numerical algorithm for estimating precise backscattering coefficients of rice fields is proposed and its accuracy is verified in this paper. After a bunch of rice plants above water surface is modeled with a bunch of randomly oriented lossy dielectric bodies above an impedance surface and the equivalent volume currents of the lossy dielectrics are computed using the moment method. Then, the scattered fields of a rice field with many bunches are computed with a Monte Carlo method, and consequently the backscattering coefficient of the rice field is computed for various incidence angles and polarizations. Finally, the backscattering coefficient of a rice field is measured at 1.85 GHz using an R-band scatterometer system, and these experimental data are used to verify the numerical algorithm proposed in this paper. It is found that the numerical computation results agree well with the measurement data.

• The Journal of the Acoustical Society of Korea
• /
• v.18 no.6
• /
• pp.31-37
• /
• 1999

This paper describes the far-field estimation using the near-field measurement data. Measurement in far-field region gives us the acoustical characteristics of the source but in general measurement is made in near-field such as acoustic water tank or anechoic chamber, so far-field acoustical characteristics of the source should be predicted from near-field data. In this case, the number of measurement points in the near field which relates to the accuracy of the predicted field and the amount of data processing, should be optimized. Existing papers say that measurement points is proportional to kL and depends on geometry and directivity of the source. But they do not give us any definite criterion for the required number of measurement points. Boundary Collocation Method which is one of the far-field prediction methods, is analyzed based on Helmholtz integral equation and Green function and it has been found that the number of measurement points is optimized as 0.54kL which is about one half of the existing results.

• Journal of the Society of Naval Architects of Korea
• /
• v.28 no.1
• /
• pp.12-18
• /
• 1991

Three-dimensional nonlinear sloshing effects due to tank motions are simulated by solving boundary value problem using the panel method based on boundary integral technique. While Shinkai used boundary elements on which source strengths vary linearly between nodes, the source of constant strength is distributed on each triangular panel in the present study. The source strength at each time step is determined by solving the Fredholm integral equation of the second kind obtained from Green's theorem. To avoid cumulative numerical errors as time elapses, Adam-Bashforth-Moulton method is employed. Numerical examples for the case of partially filled spherical tank on board oscillating in harmonic sway mode or pitch mode are solved. The elevation of the free surface is compared with the result by Shinkai and confirmed in good agreement during early time. The input and the output energy are comparatively evaluated to check the overall accuracy of the present numerical scheme. Although some leakage of energy are found as time marches, it is plausible when we take into account nonlinearities of the problem and the number of panels of the model.

• Journal of the Society of Naval Architects of Korea
• /
• v.34 no.4
• /
• pp.31-41
• /
• 1997

This paper describes a potential-based panel method formulated for the prediction of the steady performance of a waterjet propulsor. The method employs normal dipoles and sources distributed on the solid surfaces such as the impeller/stator blades, hub and duct, and normal dipoles in the shed wakes trailing the impeller and stator to represent the potential flow around the waterjet propulsor. To define a closed boundary surface, the inlet and outlet open boundary surfaces are introduced where the sources and dipoles are distributed. The kinematic boundary condition on the solid boundary surface is satisfied by requiring that the normal component of the total velocity should vanish. On the inlet surface, the total inflow flux into the duct is specified, and on the outlet surface the conservation of mass principle is applied to evaluate the source strength. The solid surfaces are discretized into a set of quadrilateral panel elements and the strengths of sources and dipoles are assumed constant at each panel. Applying this approximation to the boundary conditions leads to a set of simultaneous equations. Systematic numerical tests show that the present numerical method is fast and stable. In order to validate the present method, sample computations are carried out first for the case of a conventional axial flow fan which has a similar geometry as the waterjet propulsor, and then for the case of a waterjet propulsor on which experiments are carried out at KRISO(Korea Research Institute of Ships and Ocean Engineering).

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.12 no.1
• /
• pp.1-8
• /
• 1992

A method analizing contact pressure between plate and elastic half space is presented by using F.E.M. With the method, the pressure intensities at surface nodes of half space cae be directly calculated by using flexibility matrix of half space. The method is originally presented by Y.K. Cheung et al.(3) Insted of Y.K. Cheung's method, which use a conception of equi-contact pressure area around each surface nodes of half space in the noded rectanqular element area. We use the equi-contact pressure area around the Gaussian integration points of half space surface in the noded isoparametric element area. Numarical examples are presented and compared with other's studies.

• Journal of Korea Water Resources Association
• /
• v.30 no.6
• /
• pp.597-610
• /
• 1997

One-dimensional numerical models using finite difference methods for unsteady sediment transport on alluvial river channel are developed. The Preissmann implicit scheme and the Lax-Wendroff two-step explicit scheme with the Method of Characteristics for water motion and a forward time centered space explicit scheme for sediment motion are developed to simulate the sediment transport rate and the variation of channel bed level. The program correctness of each model is successfully verified using volume conservation tests. The sensitivity studies show that higher peak stage level, steeper channel slope and longer flooding duration produce more channel bed erosion. and median grain size, $D_{50}=0.4mm$ give maximum volume loss in this study. Finally, the numerical models are found to produce reasonable results from the various sensitivity tests which reveal that the numerical models have properly responded to the changes of each model parameter.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.10 no.4
• /
• pp.141-150
• /
• 1990

Numerical experiments of sballow water equations are performed under various boundary conditions by finite element method to simulate the circulation in estuaries and coastal areas. Galerkin method is employed to discretize spatial domain, and for time integration, finite difference method (Crank-Nicolson scheme) is used. This method is tested in five problems, in which first four cases have analytic solutions. The computed values are well in agreement with the analytic solutions in four experiments and the result of the last 2-dimensional ease is resonable. Implicit and two step Lax-Wendroff schemes in time domain are compared, and the results when using four node bilinear and triangular elements are presented. Consequently it takes very long time for complex problems requiring many elements to integrate all the time steps using the implicit schemes. And the explicit scheme requires careful consideration in selecting the time step and the grid size to obtain the desired accuracy.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.20 no.4
• /
• pp.493-499
• /
• 2007

The Taylor expansion expresses a differentiable function and its coefficients provide good approximations for the given function and its derivatives. In this study, m-th order Taylor Polynomial is constructed and the coefficients are computed by the Moving Least Squares method. The coefficients are applied to the governing partial differential equation for solid problems including crack problems. The discrete system of difference equations are set up based on the concept of point collocation. The developed method effectively overcomes the shortcomings of the finite difference method which is dependent of the grid structure and has no approximation function, and the Galerkin-based meshfree method which involves time-consuming integration of weak form and differentiation of the shape function and cumbersome treatment of essential boundary.

• Journal of the Society of Naval Architects of Korea
• /
• v.38 no.4
• /
• pp.39-47
• /
• 2001

In this paper, the hydroelastic behavior of a mat-type large floating structure is analyzed in time domain by using mode superposition method. The time-memory function is estimated by Fourier transforming the wave damping coefficients, which are computed by a higher-order boundary element method based on potential theory. Meanwhile, the structural response is obtained by time integrating the eigenmodes of the structure. Numerical examples are made for three test cases on the scaled model of a mat-type large floating structure ; weight pull-up case, weight drop case and weight moving case. In all three cases, the numerical results coincide well with experimental data.