• Journal of the Korean Society of Marine Environment & Safety
• /
• v.22 no.5
• /
• pp.564-575
• /
• 2016

The present study provides numerical simulations to predict the added resistance and ship motion of the KVLCC2 in regular waves using the unsteady Reynolds-Averaged Navier-Stokes (URANS) and 3-D potential methods. This numerical analysis is focused on added resistance and vertical ship motions (heave and pitch) under a wide range of wave conditions at three ship speeds (design, operating and zero speeds). Firstly, the characteristics of the CFD and 3-D potential flow methods are presented to predict added resistance and ship motions in regular waves taking into account various wave conditions at design speed to provide a validation study as well as at operating and zero speeds. Secondly, analyses of added resistance and ship motion with unsteady wave patterns and time history results as simulated by CFD were performed at each ship speed. Systematic validation and verifications of the numerical computations in this study were made against available Experimental Fluid Dynamics (EFD) data including grid convergence tests to demonstrate that reliable numerical results were obtained for the prediction of added resistance and ship motion in waves. Relationships between added resistance, vertical motion and changes in ship speeds were also found.

• Journal of the Society of Naval Architects of Korea
• /
• v.37 no.1
• /
• pp.67-81
• /
• 2000

An accurate and efficient numerical method for two-dimensional nonlinear radiation problem has been developed. The wave motion due to a moving body is described by the assumption of ideal fluid flow, and the governing Laplace equation can be effectively solved by the higher-order boundary element method with the help of the GMRES (Generalized Minimal RESidual) algorithm. The intersection or corner problem is resolved by utilizing the so-called discontinuous elements. The implicit trapezoidal rule is used in updating solutions at new time steps by considering stability and accuracy. Traveling waves caused by the oscillating body are absorbed downstream by the damping zone technique. It is demonstrated that the present method for time marching and radiation condition works efficiently for nonlinear radiation problem. To avoid the numerical instability enhanced by the local gathering of grid points, the regriding technique is employed so that all the grids on the free surface may be distributed with an equal distance. This makes it possible to reduce time interval and improve numerical stability. Special attention is paid to the local flow around the body during time integration. The nonlinear radiation force is calculated by the "acceleration potential technique". Present results show good agreement with other numerical computations and experiments.

• Journal of the Korean Society of Surveying, Geodesy, Photogrammetry and Cartography
• /
• v.16 no.1
• /
• pp.17-26
• /
• 1998

This paper deals with gravimetric geoid determination by Fast Hartely transform (FHT) technique in and around the Korean peninsula. A number of data files were compiled for this work, containing now more than 69, 001 point gravity data on land and ocean areas. Especially, regression was applied to estimate gravity anomalies in the northern area of peninsula. For evaluating accuracy of geoid obtained, GPS/Leveling data of 49 stations were prepared. EGM96 global geopotential model to degree 360 was used in order to determine the long wavelength effect of geoid undulations. By applying the remove-restore technique geoid undulations were determined by combining a geopotential model, free-air gravity anomalies. Fast Hartley Transform technique is a suitable solution that uses the advanced spectral technique on the sphere. It was applied to predict geoid undulations by Stokes's integral. Accuracy of geoid undulations was evaluated by comparing with results derived from GPS/Leveling. Standard deviation of differences is about 33 cm.

• Journal of the Society of Naval Architects of Korea
• /
• v.46 no.6
• /
• pp.587-594
• /
• 2009

In the paper the effect of a stern-plate attached to a ship was taken into account. The relationship between the trim angle of a ship and the wave-resistance coefficient induced by the a stern-plate was studied using the potential flow analysis method. Numerical algorithm was described using the panel method and the vortex lattice method(VLM) to simulate the flow phenomena around a ship. The non-linearity of the free surface boundary conditions were considered using the iterative method and the IGE-GMRES(Incomplete Gaussian Elimination-The Generalized Minimal RESidual) algorithm was adopted to solve the linear equation at each iterative step. Numerical calculations were carried out to investigate the validity of the adopted algorithm using KCS(KRISO 3600 TEU Container) hull. Possible cases for attachment of the plate were checked. The results showed that the numerical algorithm could be physically appropriate.

• Korean Journal of Soil Science and Fertilizer
• /
• v.36 no.4
• /
• pp.181-192
• /
• 2003

We developed a mathematical simulation model to portray the vertical distribution of soil water from the measured weather data and the known soil hydraulic properties, and then compared simulation results with the periodically measured soil water profiles obtained on Jungdong sandy loam to verify the model, In this model, we solved potential-based Richards' equation by the implicit finite difference method superimposed on the predictor-corrector scheme. We presumed that: soil hydraulic properties are homogeneous; soil water flows isothermally; hysteresis is not considered; no vapor flows; no heat transfers into the soil profiles; and water added to soil surface is distributed along the soil profile following partial displacement principle. The input data were broadly classified into two groups: (1) daily weather data such as rainfall, maximum and minimum air temperatures, relative humidity and solar radiation and (2) soil hydraulic data to approximate unsaturated hydraulic conductivity and water retention. Each hydraulic polynomial function approximated using the Chebyshev polynomial and least square difference technique in tandem showed a fairly good fit of the given set of data. Vertical distribution of soil water as approximations to the Richards' equation subject to changing surface and phreatic boundaries was solved numerically during 53 days with a comparatively large time increment, and this pattern agreed well with field neutron scattering data, except for the surface 0.1 m slab.

• Journal of Korean Society for Geospatial Information Science
• /
• v.16 no.3
• /
• pp.117-126
• /
• 2008

Earth geopotential models were used to determine the reference surface in geoid modelling and many global geopotential models were developed after 1980's. Nowadays, EGM96 and EIGEN-GL04C geopotential models have been most widely used in the world, but what so called EGM2008 earth geopotential model were developed in 2008 by NGA. In this paper, we intended to compare the results of spherical harmonic analyses using the three geopotential model, EGM96, EIGEN-GL04C and EGM2008. So, the spherical harmonic analyses were performed up to degree and order 360(in case of EGM2008, up to degree and order 720, 1440, 2190 in addition), on each $1'{\times}1'$grid point in and around Korean peninsula. Geometric geoid were calculated at 464 GPS/leveling points for accuracy evaluation and then the results of three geopotential models were compared to geometric geoid. The results show that the accuracy of EGM2008 is improved considerablely compared to EGM96 and EIGEN-GL04C and it is possible to calculate geoidal heights within 14cm standard deviation and 5.5cm standard deviation after LSC fitting in and around Korean peninsula using EGM2008 geopotential model.

• Progress in Medical Physics
• /
• v.10 no.1
• /
• pp.1-7
• /
• 1999

In this study, the dose distributions of a $^{32}$ p uniform cylindrical volume source and a surface source, a pure $\beta$emitter, were calculated in order to obtain information relevant to the utilization of a balloon catheter and a radioactive stent. The dose distributions of $^{32}$ p were calculated by means of the EGS4 code system. The sources are considered to be distributed uniformly in the volume and on the surface in the form of a cylinder with a radius of 1.5 mm and length of 20 mm. The energy of $\beta$particles emitted is chosen at random in the $\beta$ energy spectrum evaluated by the solution of the Dirac equation for the Coulomb potential. Liquid water is used to simulate the particle transport in the human body. The dose rates in a target at a 0.5mm radial distance from the surface of cylindrical volume and surface source are 12.133 cGy/s per GBq (0.449 cGy/s per mCi, uncertainty: 1.51%) and 24.732 cGy/s per GBq (0.915 cGy/s per mCi, uncertainty: 1.01%), respectively. The dose rates in the two sources decrease with distance in both radial and axial direction. On the basis of the above results, the determined initial activities were 29.69 mCi and 1.2278 $\mu$Ci for the balloon catheter and the radioactive stent using $^{32}$ P isotope, respectively. The total absorbed dose for optimal therapeutic regimen is considered to be 20 Gy and the treatment time in the case of the balloon catheter is less than 3 min. Absorbed doses in targets placed in a radial direction for the two sources were also calculated when it expressed initial activity in a 1 mCi/ml volume activity density for the cylindrical volume source and a 0.1 mCi/cm$^2$ area activity density for the surface source. The absorbed dose distribution around the $^{32}$ P cylindrical source with different size can be easily calculated using our results when the volume activity density and area activity density for the source are known.