• 한국해안·해양공학회논문집
• /
• 제31권3호
• /
• pp.180-198
• /
• 2019

최근, 유럽이나 일본 등은 해빈유실방지대책공법의 하나이며, 잠제의 대안으로 여겨지는 저천단구조물(LCS)에 관한 많은 연구를 수행하였고, 그의 결과들을 집약하여 설계매뉴얼까지 편찬하였다. 지금까지 LCS에 관한 연구는 2차원적인 파랑전달율과 피복재의 안정중량산정에 치우쳐 있으며, 이들은 주로 실험에 기초하여 검토 논의되었다. 본 연구에서는 투과성의 LCS를 대상으로 3차원수치해석을 수행한다. 수치해석에 Navier-Stokes 운동방정식에 기반한 오픈소스 CFD Code인 olaFlow를 적용하였으며, 이는 쇄파와 난류해석까지도 가능한 강비선형해석법이다. 이로부터 수위, 흐름 및 난류운동에너지 등의 분포특성을 검토 논의하였으며, 동시에 잠제의 경우와도 비교 검토하였다. 이로부터 해안선 근방에서 연안류의 흐름패턴과 평균난류운동에너지의 연안방향 및 종단방향의 공간분포에 관해 잠제와 LCS의 경우에 각각 차이가 발생하는 것을 알 수 있었다. 이러한 결과의 차이는 모래이동에서의 차이로 이어질 것으로 판단된다.

• 전자공학회논문지B
• /
• 제31B권11호
• /
• pp.53-62
• /
• 1994

본 논문에서는 湖光立體視法(Photometric Stereo Method : PSM)을 적용하여 3차원 물체의 표면 방향정보를 추정할 때 조명원이 이상적인 평행광이 아닌 범용 조명환경인데서 기인한 표면 방향추정 오차의 해석 방법을 제시하였다. 제시한 방법에서는 원반 형태의 조사 대상 면을 포함한 조명환경에 대한 구체적인 모델링과 오차분포에 대한 직접적인 시뮬레이션을 통하여 오차해석을 하도록 하였다. 범용조명원은 beam angle을 갖는 점광원으로 가정하고, 조사 대상면에서의 밝기분포는 분산값이 변하는 가우스함수 형태로 모델링하였다. 설정된 조명 모델에서 분산을 여러값으로 변화시켜서 이에따라 PSM에서 나타나는 오차분포를 얻기 위한 목적으로, 임의의 세 방향으로 위치한 조사 대상면의 평균 밝기값을 계산하여 고정된 조명원의 위치벡터를 얻는 시뮬레이션 알고리듬을 제시하였다. PSM과 제안된 조명원의 위치추정 알고리듬 사이에는 상사관계(analogy)가 성립하기 때문에, PSM 적용시 비이상적인 조명모델에 의해 나타나는 물체표면방향 추정오차는 제시한 조명원 위치추정 시뮬레이션으로부터 직접적으로 구할 수 있다. 또한 설정된 조사대상면 모델과 같은 원방형태의 계측기구를 제작하여 임의의 방향벡터를 제공하도록 하고, 이를 측정과 시뮬레이션에 사용함으로써 오차해석과 보정을 체계적으로 수행할 수 있도록 하였다. 제시한 방법에 따라 방향 벡터를 제공하도록 하고, 이를 측정과 시뮬레이션에 사용함으로써 오차해석과 보정을 체계적으로 수행할 수 있도록 하였다. 제시한 방법에 따라 방향 벡터값들의 조합을 공간상에서 다양하게 변화시키면서, PSM에서 나타나는 방향추정 오차분포를 구하는 컴퓨터 시뮬레이션을 수행하였다. 시뮬레이션과 같은 방법으로 실측에 의한 방향추정 오차를 구하고 시뮬레이션에서 얻은 결과와 비교 검토하였다. 그 결과, 실제 측정실험에서 얻은 복합적인 오차중 상당부분은 비이상적인 범용조명원 모델로부터 기인한 것임을 확인할 수 있었다.

• 한국해양공학회지
• /
• 제15권4호
• /
• pp.14-19
• /
• 2001

In this paper, the wave and current forces acting on cylinders are investigated by theoretical and experimental methods. The models used are one-cylinder, four-cylinder and semi-submersible types. The theoretical investigations are carried out by the Morison equation and three dimensional source distribution method to calculate exciting forces in waves with and without currents. The experimental investigations are carried out in the wave tank which can generate currents in both directions. In these tests, the models have been exposed to the regular waves with and without currents. It is shown that the exciting forces acting on the one-cylinder or four-cylinders can be approximately estimated by the Morison equation and also by the diffraction theory. However, the Morison equation seems to be not appropriate to estimate the exciting forces on the present type of semi-submersible.

• 동력기계공학회지
• /
• 제18권3호
• /
• pp.36-41
• /
• 2014

The structure of a variable liquid column oscillator(a VLCO) is analogous to that of the tuned liquid column damper used to suppress oscillatory motion in large structures like tall buildings and cargo ships. The VLCO is a system absorbing high kinetic energy of accelerated motions of the multiple floating bodies in the effect of air springs occurred by installation of inner air chambers. Thus, VLCO can improve the efficiency of energy than wave energy converters of the activating object type made in Pelamis Company. In this research, the experiment was performed in two models of same draft. The one is that weights were filled, and the other is that water was filled. The numerical results were estimated by assuming that do not exist internal flow, and the results were compared with the results of experiments.

• 수산해양기술연구
• /
• 제37권4호
• /
• pp.313-320
• /
• 2001

In the presence of incident waves with different frequencies, there are second order sum and difference frequency wane exciting forces due to the nonlinearty of the incident waves. Although the magnitudes of these nonlinear forces are small, they act at sum and difference frequencies away from those of the incident waves. So, the second order sum and difference frequency wane loads occurring close to the natural frequencies of TLPs often give greater contributions to high and low frequency resonant responses. The components of the second order forces which depend on first order quantities have been evaluated using the three dimensional source distribution method. The numerical results of time domain motion analysis for the nonlinear wave exciting forces in regular waves are compared with the numerical ones of frequency domain analysis. The results of comparison confirmed the validity of the proposed approach.

• 한국해양공학회지
• /
• 제20권4호
• /
• pp.70-75
• /
• 2006

A numerical procedure is described for estimating the effects of the multi-directional irregular waves on the motion responses and tension variations of the ISSC-TLP. The numerical approach is based on a three-dimensional source distribution method and a spectral analysis technique of directional waves. The spectral description for the linear system of ISSC-TLP in the frequency domain is sufficient to completely define the motion responses and tension variations. This is because both the wave inputs and responses are stationary Gaussian random processes, of which the statistical properties in the amplitude domain are well known. The numerical results for the linear motion responses and tension variations in regular waves are compared with the experimental and numerical ones, which are obtained in the literature. The results of comparison confirmed the validity of the proposed approach.

• 한국에너지공학회:학술대회논문집
• /
• 한국에너지공학회 2000년도 추계 학술발표회 논문집
• /
• pp.115-120
• /
• 2000

The analysis of flux distribution under stead-state in large power reactors with assymetry reactivity insertions requires the use of three-dimensional diffusion calculations. For the purpose, consistently formulated modern nodal methods based on higher order interface techniques have become popular tools for flux distributions in large commercial nuclear reactors. Among the earlier developments, the nodal Green's function method obtains its nodal interface equation from the transverse-integrated integral diffusion equation using a finite-medium Green's function. In this method, the outgoing current from a node surface is formulated as a response of the incoming currents and the spatially integrated neutron source within the same node. The well-known nodal expansion method is also based on an interface partial current formulation. Nodal methods high-level interface variables, i.e., interface net current and flux, may be more computationally efficient than the nodal Green's function method because they have one fewer unknown per interface. The Analytic Nodal Method(ANM), which can be classified as an interface net current technique and, was faster in solving some standard benchmark problems than the other two methods.(omitted)

• 한국의학물리학회지:의학물리
• /
• 제22권4호
• /
• pp.190-197
• /
• 2011

고선량률 근접치료에 사용되는 상업용 선원과 치료계획 시스템들은 AAPM TG 43에서 권고하는 점 및 선 선원에 의해 선량분포를 계산한다. 하지만, 근접치료용 선원에 대한 인체 내의 정확한 선량계산을 위해서 3차원 부피의 선원을 고려하는 MC 기반의 선량계산 방법이 필요하다. 본 연구에서는 microSelectron HDR Ir-192 선원을 작은 부분으로 분할하여 계산하는 미소선원 적분법을 이용하여 기하학적 인수를 계산하였다. 또한, 범용 방사선 수송코드인 MCNPX를 사용하여 30 cm 직경의 구형 물 팬텀 내에서 선원의 선량률을 계산하여 비등방성함수와 반경선량함수를 구하였다. 그 결과를 MC 기반 광자 수송코드인 MCPT를 사용하여 계산한 Williamson의 결과와 비교 및 분석하였다. 미소선원 적분법과 선 선원 근사법에 따른 기하학적 인수는 $r{\geq}0.5cm$에서는 0.2% 이내에서 일치하였고 r=0.1 cm일 때 1.33%의 차이를 보였다. 본 연구에서 계산된 비등방성함수와 반경선량함수가 Williamson의 계산된 결과의 차이는 비등방성함수의 경우 r=0.25 cm에 서 2.33%의 가장 큰 R-RMSE를 보였고 $r{\geq}0.5cm$에서는 1% 미만의 R-RMSE를 보였다. 반경선량함수의 경우는 r=0.1~14.0 cm에서 0.46%의 R-RMSE를 보였다. 미소선원 적분법과 선 선원 근사법으로 계산한 기하학적 인수는 $r{\geq}0.1cm$에서 잘 일치하지만 3차원의 Ir-192 선원을 적용하여 계산한 미소선원 적분법이 실제 기하학적 인수를 잘 반영할 것으로 생각된다. r=0.25 cm에서 비등방성함수를 제외하고는 MCPT와 MCNPX의 몬테칼로 코드를 이용하여 얻어진 비등방성함수와 반경선량함수는 각각의 몬테칼로 코드에 대한 불확실성 이내에서 잘 일치함을 확인하였다. 따라서 MCNPX 전산모사 결과를 통해 TG-43의 선량 계산식에 사용된 인자를 Williamson 등의 결과와 비교 및 검증함으로써, 추후 다른 종류의 선원에 대해서도 Monte Carlo 기반의 연구가 가능할 것으로 기대된다.

• 한국진공학회:학술대회논문집
• /
• 한국진공학회 2014년도 제46회 동계 정기학술대회 초록집
• /
• pp.221-221
• /
• 2014

The present work proposes an improved numerical simulator for design and modification of large area capacitively coupled plasma (CCP) processing chamber. CCP, as notoriously well-known, demands the tremendously huge computational cost for carrying out transient analyses in realistic multi-dimensional models, because electron dissociations take place in a much smaller time scale (${\Delta}t{\approx}10-8{\sim}10-10$) than time scale of those happened between neutrals (${\Delta}t{\approx}10-1{\sim}10-3$), due to the rf drive frequencies of external electric field. And also, for spatial discretization of electron flux (Je), exponential scheme such as Scharfetter-Gummel method needs to be used in order to alleviate the numerical stiffness and resolve exponential change of spatial distribution of electron temperature (Te) and electron number density (Ne) in the vicinity of electrodes. Due to such computational intractability, it is prohibited to simulate CCP deposition in a three-dimension within acceptable calculation runtimes (<24 h). Under the situation where process conditions require thickness non-uniformity below 5%, however, detailed flow features of reactive gases induced from three-dimensional geometric effects such as gas distribution through the perforated plates (showerhead) should be considered. Without considering plasma chemistry, we therefore simulated flow, temperature and species fields in three-dimensional geometry first, and then, based on that data, boundary conditions of two-dimensional plasma discharge model are set. In the particular case of SiH4-NH3-N2-He CCP discharge to produce deposition of SiNxHy thin film, a cylindrical showerhead electrode reactor was studied by numerical modeling of mass, momentum and energy transports for charged particles in an axi-symmetric geometry. By solving transport equations of electron and radicals simultaneously, we observed that the way how source gases are consumed in the non-isothermal flow field and such consequences on active species production were outlined as playing the leading parts in the processes. As an example of application of the model for the prediction of the deposited thickness uniformity in a 300 mm wafer plasma processing chamber, the results were compared with the experimentally measured deposition profiles along the radius of the wafer varying inter-electrode gap. The simulation results were in good agreement with experimental data.

• 대한조선학회논문집
• /
• 제37권2호
• /
• pp.57-69
• /
• 2000

기하학적 형상과 유동의 해를 B-스플라인으로 표현하는 3차원 고차 패널법은 프로펠러 주위의 유동을 해석하기 위해 개발되었다. 제어점이 패널내에 놓이는 경우, 고차의 다이폴과 쏘오스에 의해 유기되는 자기 유기 포텐셜의 특이 거동은 2차 변환(quadratic transformation)을 통하여 제거하였으며, 특이 부분은 해석적인 적분으로 비특이 부분은 정도 높은 Gauss 구적법으로 계산함으로써 유기 포텐셜을 정도 높게 구할 수 있음을 보였다. 또한, 날개 뒷날에서의 압력 점프의 값이 명시적으로 영이 되도록하는 동역학적 Kutta 조건식을 도입하고, 이의 적용이 안정된 해를 보장함을 확인하였다. 수치 실험을 통하여, 제안된 수치해석 기법이 안정적이고 정확한 해를 줌을 확인하였으며, 특히 저차 패널법과 비교하여 적은 수의 패널로 동일한 정도의 해를 유지할 수 있음을 보였다.